Jet.h

Go to the documentation of this file.

#pragma once
 
// STD
#include <cstdint>
#include <ios>
 
// Darwin
#include <GenericObject.h>
 
namespace Darwin::Physics {
 
class GenJet : public GenericObject {
 
    std::string_view scale_group () const override { return ScaleVar; }
    std::string_view weight_group () const override { return WeightVar; }
 
public:
    static inline const char * const ScaleVar = "GenJetScales",
                             * const WeightVar = "GenJetWgts";
 
    int nBmesons = -1, 
        nDmesons = -1; 
    int partonFlavour = 0; 
 
    GenJet () = default;
    virtual ~GenJet () = default;
};
 
class RecJet : public GenJet {
 
    std::string_view scale_group () const override { return ScaleVar; }
    std::string_view weight_group () const override { return WeightVar; }
 
public:
    static inline const char * const ScaleVar = "RecJetScales",
                             * const WeightVar = "RecJetWgts";
 
    float area = M_PI*0.4*0.4; 
    std::uint32_t selectors = 0; 
 
    RecJet () = default;
    virtual ~RecJet () = default;
};
 
} // end of Darwin::Physics namespace
 
inline std::ostream& operator<< (std::ostream& s, const Darwin::Physics::GenJet& jet)
{
    return s << jet.CorrP4() << ' ' << jet.nDmesons << ' ' << jet.nBmesons;
}
 
inline std::ostream& operator<< (std::ostream& s, const Darwin::Physics::RecJet& jet)
{
    return s << jet.CorrP4() << " 0x" << std::hex << jet.selectors << std::dec;
}