SingleLeptonic.hh

Go to the documentation of this file.

#pragma once

#include "toppartner/analysis/TopPartner.hh"
#include "toppartner/tagger/EventSingleLeptonic.hh"
#include "toppartner/tagger/TopPartnerHadronicNeutral.hh"
#include "boca/generic/DEBUG_MACROS.hh"

namespace toppartner
{

namespace analysis
{

template<typename Tagger_>
class SingleLeptonic : public TopPartner<Tagger_>
{

protected:

    std::string Name() const override {
        INFO0;
        return "Single-Leptonic-" + boca::Name(Settings::Collider()) + "-" + boca::units::Name(this->Mass()) + "-new-higgs";
    }

    void SetFiles(Phase const& phase) override {
        INFO0;
        switch (phase.Tag()) {
        case Tag::signal :
            if (this->template TaggerIs<tagger::VetoTopPartnerHadronic>() || this->template TaggerIs<tagger::TopPartnerHadronicNeutral>()) this->AddSignal(Process::TT);
            else this->AddSignal(Process::TthLep);
            break;
        case Tag::background :
            if (this->template TaggerIs<tagger::VetoTopPartnerHadronic>() || this->template TaggerIs<tagger::TopPartnerHadronicNeutral>()) this->AddBackground(Process::TthLep);
            else if (!this->template TaggerIs<tagger::TopPartnerLeptonicNeutral>()) this->AddBackground(Process::TT);
            if (!this->template TaggerIs<tagger::VetoTopPartnerHadronic>() || this->template TaggerIs<tagger::TopPartnerHadronicNeutral>()) {
                this->AddBackground(Process::ttBB);
                this->AddBackground(Process::ttBjj);
            }
            break;
        }
    }

private:
//   TopPartnerLeptonicNeutralTagger partner_tagger_;

    bool PassPreCut(boca::Event const&) const override {
        INFO0;
//     if(tag == Tag::signal){
//      std::vector<Jet> partner = partner_tagger.Particles(event);
//       if(partner.empty()) return false;
//     }
//        std::vector<Jet> jets = SortedByPt(event.Jets());
//         if (jets.size() < 3) return false;
//         if (jets.at(2).Pt() < this->JetPreCut()) return false;
//
//
//        std::vector<Jet> particles = event.GenParticles();
//         particles = RemoveIfSoft(particles, this->PreCut());
//        std::vector<Jet> tops = CopyIfParticle(particles, Id::top);
//        std::vector<Jet> higgs = CopyIfParticle(particles, Id::higgs);
//        std::vector<Jet>vectors = CopyIfParticles(particles, Id::Z, Id::W);
//        std::vector<Jet> partner = CopyIfParticle(particles, Id::top_partner);
//         if (tag == Tag::signal && partner.size() != 1) {
//             return false;
//         }
//         if (tops.size() < 2 || (higgs.empty() && vectors.empty())) {
//             return false;
//         }
        return true;
    }

};

}

}