TopHadronic.hh

Go to the documentation of this file.

#pragma once

#include <functional>
#include "boca/multiplets/Triplet.hh"
#include "standardmodel/tagger/WHadronic.hh"
#include "standardmodel/branch/TopHadronic.hh"

namespace boca
{
class MomentumRange;
}

namespace standardmodel
{

namespace tagger
{

class TopHadronic : public Tagger<Triplet, branch::TopHadronic>
{

public:

    TopHadronic(Id id = Id::top);

    int Train(boca::Event const& event, PreCuts const& pre_cuts, Tag tag) override;

    using Tagger::Multiplets;

    std::vector<Triplet> Multiplets(boca::Event const& event, PreCuts const& pre_cuts, TMVA::Reader const& reader) override;

    std::string Name() const override;

    latex::String LatexName() const override;

    virtual std::vector<Particle> Particles(boca::Event const& event) const;

    using Function = std::function<boost::optional<Triplet>(Triplet&, std::vector<boca::Jet> const&)>;

    std::vector<Triplet> Triplets(boca::Event const& event, Function const& function);

private:

    std::vector<Triplet> ThreeJets(const std::vector< boca::Jet >& jets, const std::vector< boca::Lepton >& leptons, const Function& function, const boca::MomentumRange& range);

    std::vector<Triplet> TwoJets(const std::vector< boca::Jet >& jets, boca::Jet const& jet, const std::vector< boca::Lepton >& leptons, const Function& function, const boca::MomentumRange& range);

    std::vector<Triplet> ThreeSubJets(boca::Jet const& jet, const std::vector< boca::Lepton >& leptons, const Function& function, const boca::MomentumRange& range);

    std::vector<Triplet> TwoSubJets(boca::Jet const& jet, const std::vector< boca::Lepton >& leptons, const Function& function, const boca::MomentumRange& range);

    boost::optional<Triplet> HighlyBoosted(boca::Jet const& jet, std::vector<Lepton> const& leptons, Function const& function) const;

    boost::optional<Triplet> Tripple(Doublet const& doublet, boca::Jet const& jet, std::vector<Lepton> const& leptons, Function const& function, MomentumRange const& range, bool check_overlap = false) const;

    boost::optional<Triplet> Tripple(Triplet& triplet, std::vector<Lepton> const& leptons, PreCuts const& pre_cuts, Tag tag) const;

    std::vector<Triplet> Triplets(std::vector<Doublet> const& doublets, std::vector<boca::Jet> const& jets, std::vector<Lepton> const& leptons, Function const& function, MomentumRange const& range) const;

    std::vector<Triplet> Triplets(Doublet const& doublet, std::vector<boca::Jet> const& jets, std::vector<Lepton> const& leptons, Function const& function, MomentumRange const& range) const;

    boost::optional<Triplet> Multiplet(Triplet& triplet, std::vector<Lepton> const& leptons, PreCuts const& pre_cuts, TMVA::Reader const& reader);

    bool Problematic(Triplet const& triplet, PreCuts const& pre_cuts, Tag tag) const;

    bool Problematic(Triplet const& triplet, PreCuts const& pre_cuts) const;

    Reader<tagger::Bottom> bottom_reader_;

    Reader<tagger::WHadronic> w_hadronic_reader_;

    Mass top_mass_window_;

    Id id_;

};

}

}