Top.hh

Go to the documentation of this file.

#pragma once

#include "boca/fastjet/Particles.hh"

#include "standardmodel/analysis/StandardModel.hh"
#include "standardmodel/tagger/WHadronic.hh"
#include "standardmodel/tagger/WLeptonic.hh"

namespace standardmodel
{

enum class TopTagger
{
    boca,
    hep,
    hep2
};

std::string Name(TopTagger top_tagger);

namespace analysis
{

template<typename Tagger_>
class Top : public StandardModel<Tagger_>
{

public:

    Top() {
        this->PreCuts().PtLowerCut().Set(Id::top, this->LowerPtCut());
        this->PreCuts().PtUpperCut().Set(Id::top, this->UpperPtCut());
        this->PreCuts().TrackerMaxEta().Set(Id::top, Settings::TrackerEtaMax());
    }

    static Decay TopDecay() {
        return Decay::hadronic;
        return Decay::leptonic;
    }

    static standardmodel::TopTagger TopTagger() {
        return standardmodel::TopTagger::boca;
        return standardmodel::TopTagger::hep;
        return standardmodel::TopTagger::hep2;
    }

private:

    std::string Name() const override {
      return boca::Name(this->Collider()) + "-" + boca::units::Name(this->LowerPtCut()) + "-" + boca::Name(TopDecay()) + (TopDecay() == Decay::hadronic ? "-" + standardmodel::Name(TopTagger()) : "") + "new-est";
    }

    void SetFiles(Phase const& phase) override {
        switch (phase.Tag()) {
        case Tag::signal :
            if (TopDecay() == Decay::hadronic || this->template TaggerIs<tagger::Bottom>()) this->AddSignal(Process::tt_had);
            if (this->template TaggerIs<tagger::Bottom>()) this->AddSignal(Process::hh);
            if (this->template TaggerIs<tagger::Bottom>()) this->AddSignal(Process::bb);
            if (TopDecay() == Decay::leptonic || this->template TaggerIs<tagger::Bottom>()) this->AddSignal(Process::tt_lep);
            if (this->template TaggerIs<tagger::WHadronic>() || this->template TaggerIs<tagger::WLeptonic>()) this->AddSignal(Process::ww);
            break;
        case Tag::background :
            if (TopDecay() == Decay::leptonic && !this->template TaggerIs<tagger::Bottom>()) this->AddBackground(Process::tt_had);
            if (!this->template TaggerIs<tagger::Bottom>()) this->AddBackground(Process::hh);
            if (!this->template TaggerIs<tagger::Bottom>()) this->AddBackground(Process::bb);
            this->AddBackground(Process::cc);
            this->AddBackground(Process::gg);
            if (!this->template TaggerIs<tagger::WLeptonic>()) this->AddBackground(Process::qq);
            if ((TopDecay() == Decay::hadronic) && !this->template TaggerIs<tagger::Bottom>()) this->AddBackground(Process::tt_lep);
            this->AddBackground(Process::zz);
            if (!this->template TaggerIs<tagger::WHadronic>() && !this->template TaggerIs<tagger::WLeptonic>()) this->AddBackground(Process::ww);
            break;
        }
    }

    bool PassPreCut(boca::Event const& event) const override {
        auto particles = SortedByPt(event.GenParticles());
        particles = CopyIfDrellYan(particles);
        particles = RemoveIfOutsidePtWindow(particles, this->LowerPtCut(), this->UpperPtCut());
//         std::cout << "pre cut " << particles.size() << '\n';
        return particles.size() == 1 ? true : false;
    }

};

}

}