Jet.hh

Go to the documentation of this file.

#pragma once

#include "boca/fastjet/PseudoJet.hh"
#include "boca/fastjet/JetInfo.hh"

namespace boca
{

class Jet : public PseudoJet
            , public boost::additive<Jet>
            , public boost::additive<Jet,  Particle>
            , public boost::multiplicative<Jet, double>
{

public:

    Jet();

    Jet(const boca::Momentum& x, const boca::Momentum& y, const boca::Momentum& z, const boca::Energy& e);

    Jet(PseudoJet const& jet);

    Jet(fastjet::PseudoJet const& jet);

    Jet(fastjet::PseudoJet const& jet, JetInfo const& info);

    Jet(boca::Vector3<boca::Momentum> const& spatial, const boca::Energy &e);

    Jet(Vector2<boca::Momentum> const& transverse, Momentum const& z, boca::Energy const& e);

    Jet(boca::LorentzVector<Momentum> const& lorentz_vector);

    Jet(TLorentzVector const& lorentz_vector);

    Jet(TLorentzVector const& lorentz_vector, ::delphes::Jet const& jet);

    Jet(TLorentzVector const& lorentz_vector, Constituent const& constituent);

    Jet(TLorentzVector const& lorentz_vector, std::vector<Constituent> const& constituents);

    Jet(TLorentzVector const& lorentz_vector, int charge);

    Jet(exroot::Electron const& electron);

    Jet(exroot::GenJet const& gen_jet);

    Jet(exroot::Jet const& jet);

    Jet(exroot::Muon const& muon);

    Jet(exroot::Photon const& photon);

    Jet(exroot::Tau const& tau);

    JetInfo const& Info() const;

    JetInfo& Info();

    void SetDelphesTags(::delphes::Jet const& delphes_jet);

    bool HasConsitutents() const;

    std::vector<Jet> Constituents() const;

    double Bdt() const;

    void SetInfo(JetInfo const& info = JetInfo());

    Jet& operator+=(Jet const& jet);

    Jet& operator-=(Jet const& jet);

    Jet& operator*=(double scalar);

    Jet& operator/=(double scalar);

    Jet& operator+=(Particle const& particle);

    Jet& operator-=(Particle const& particle);


    bool HasJetInfo() const;

protected:

    void ResetInfo(JetInfo const& info);

private:

};

using Lepton = Jet;

using Photon = Jet;

using MissingEt = Jet;

std::vector<Jet> JetVector(std::vector<fastjet::PseudoJet> const& pseudo_jets);

std::vector<fastjet::PseudoJet> PseudoJetVector(std::vector<Jet> const& jets);

Jet Join(std::vector<Jet> const& jets);

Jet Join(Jet const& jet_1, Jet const& jet_2);

Jet Join(Jet const& jet_1, Jet const& jet_2, Jet const& jet_3);

}