Jet.h

Go to the documentation of this file.

#pragma once
 
#include "Core/Objects/interface/PhysicsObject.h"
 
namespace DAS {
 
struct GenJet : public PhysicsObject {
    static inline const char * const ScaleVar = "GenJetScales", 
                             * const WeightVar = "GenJetWgts";  
 
    int nBHadrons = -1, 
        nCHadrons = -1; 
    int partonFlavour = 0; 
 
    inline int hadronFlavour () const
    {
             if (nBHadrons > 0) return 5;
        else if (nCHadrons > 0) return 4;
        else                    return 0;
    }
 
private:
    std::string_view scale_group () const override { return ScaleVar; }
    std::string_view weight_group () const override { return WeightVar; }
};
 
struct RecJet : public GenJet {
    static inline const char * const ScaleVar = "RecJetScales", 
                             * const WeightVar = "RecJetWgts"; 
 
    // for jet energy scale corrections
    float area = -1, 
          puID = -1; 
 
    struct Tagger {
        float probc = -1,    
              probb = -1,    
              probbb = -1,   
              problepb = -1, 
              probuds = -1,  
              probg = -1;    
 
        inline float B    () const { return probb+probbb+problepb; }
        inline float CvsL () const { return probc/(probc+probuds+probg); }
        inline float CvsB () const { return probc/(probc+probb+probbb+problepb); }
    } DeepJet; 
 
    RecJet () = default; 
 
private:
    std::string_view scale_group () const final { return ScaleVar; }
    std::string_view weight_group () const final { return WeightVar; }
};
 
} // end of DAS namespace
 
inline std::ostream& operator<< (std::ostream& s, const DAS::GenJet& jet)
{
    return s << jet.p4 << ' ' << jet.nCHadrons << ' ' << jet.nBHadrons;
}
 
#if defined(__ROOTCLING__)
#pragma link C++ class DAS::GenJet +;
#pragma link C++ class std::vector<DAS::GenJet> +;
 
#pragma link C++ class DAS::RecJet +;
#pragma link C++ class DAS::RecJet::Tagger +;
#pragma link C++ class std::vector<DAS::RecJet> +;
#endif