Multiplet.hh

Go to the documentation of this file.

#pragma once

#include "boca/multiplets/Singlet.hh"
#include "boca/fastjet/Particle.hh"
#include "boca/multiplets/ClosestLepton.hh"
#include "boca/fastjet/SubJettiness.hh"
#include "boca/event/EventShapes.hh"

namespace boca
{

class Multiplet : public Identification
{

    template<typename Multiplet_>
    using NotJet = typename std::enable_if < !std::is_same<Multiplet_, boca::Jet>::value && !std::is_same<Multiplet_, boca::PseudoJet>::value && !std::is_same<Multiplet_, boca::Particle>::value >::type;

public:

    void SetClosestLepton(std::vector<boca::Lepton> const &leptons);

    ClosestLepton Lepton() const;

    boca::Jet Jet() const;

    boca::Mass Mass() const;

    Momentum Pt() const;

    Angle Rap() const;

    Angle Phi() const;

    template <typename Multiplet_, typename = NotJet<Multiplet_>>
    Angle DeltaPhiTo(Multiplet_ const &multiplet) const
    {
        return Jet().DeltaPhiTo(multiplet.Jet());
    }
    Angle DeltaPhiTo(PseudoJet const &jet) const;

    template <typename Multiplet_, typename = NotJet<Multiplet_>>
    Angle DeltaRapTo(Multiplet_ const &multiplet) const
    {
        return Jet().DeltaRapTo(multiplet.Jet());
    }
    Angle DeltaRapTo(PseudoJet const &jet) const;

    template <typename Multiplet_, typename = NotJet<Multiplet_>>
    Angle DeltaRTo(Multiplet_ const &multiplet) const
    {
        return Jet().DeltaRTo(multiplet.Jet());
    }
    Angle DeltaRTo(PseudoJet const &jet) const;

    template <typename Multiplet_, typename = NotJet<Multiplet_>>
    Vector2<Angle> DeltaTo(Multiplet_ const &multiplet) const
    {
        return Jet().DeltaTo(multiplet.Jet());
    }
    Vector2<Angle> DeltaTo(PseudoJet const &jet) const;

    Vector2<Angle> Angles(bool wrap_phi = false) const;

    template <typename Multiplet_/*, typename = NotJet<Multiplet_>*/>
    Vector2<Angle> AnglesMinTo(Multiplet_ const &multiplet) const
    {
        return Jet().AnglesMinTo(multiplet.Jet());
    }
    Vector2<Angle> AnglesMinTo(PseudoJet const& jet) const;

    boca::Singlet ConstituentJet() const;

    bool HasConstituents() const;

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

    boca::SubJettiness SubJettiness() const;

    void SetExtraInfo(double extra_info);

    double ExtraInfo() const;

protected:

    virtual std::string Name() const;

    virtual Singlet GetConstituentJet() const = 0;

    virtual boca::Jet GetJet() const = 0;

    virtual std::vector<boca::Jet> Jets() const = 0;

    virtual std::vector<LorentzVector<Momentum>> LorentzVectors() const = 0;

    Mutable<boca::EventShapes> event_shapes_;

private:

    ClosestLepton closest_lepton_;

    Mutable<boca::Singlet> constituent_jet_;

    Mutable<boca::Jet> jet_;

    double extra_info_;

};

}