Photon.h

Go to the documentation of this file.

#pragma once
 
#include "Core/Objects/interface/PhysicsObject.h"
 
#include <TString.h>
 
#include <cstdint>
 
namespace DAS {
 
struct GenPhoton : public PhysicsObject {
    static inline const char * const ScaleVar = "GenPhotonScales", //< Name of photon rec scales in variations
                             * const WeightVar = "GenPhotonWgts"; //< Name of photon rec weights in variations
 
    bool zAncestor = false, 
         prompt = false; 
 
private:
    std::string_view scale_group () const override { return ScaleVar; }
    std::string_view weight_group () const override { return WeightVar; }
};
 
struct RecPhoton : public GenPhoton {
    static inline const char * const ScaleVar = "RecPhotonScales", //< Name of photon rec scales in variations
                             * const WeightVar = "RecPhotonWgts"; //< Name of photon rec weights in variations
 
    enum EnergyVariation {
        Nominal = 0,    
        ScaleUp,        
        ScaleDown,      
        SigmaPhiUp,     
        SigmaPhiDown,   
        SigmaRhoUp,     
        SigmaRhoDown    
    };
    static inline const std::vector<TString> uncs {"Scale", "SigmaPhi", "SigmaRho"};
 
    enum Identification {
        CutBasedLoose              = 0b0000001, 
        CutBasedMedium             = 0b0000010, 
        CutBasedTight              = 0b0000100, 
        MVAWorkingPoint80          = 0b0001000, 
        MVAWorkingPoint90          = 0b0010000, 
        ConversionSafeElectronVeto = 0b0100000, 
        PixelSeedVeto              = 0b1000000, 
    };
 
    float scEta, 
          sigmaIEtaIEta, 
          hOverE, 
          chargedIsolation, 
          neutralHadronIsolation, 
          photonIsolation, 
          worstChargedIsolation; 
 
    float ecalEnergyErrPostCorr = -1; 
 
    std::uint32_t selectors = 0; 
 
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::GenPhoton& photon)
{
    return s << photon.p4 << ' ' << photon.zAncestor << ' ' << photon.prompt;
}
 
inline std::ostream& operator<< (std::ostream& s, const DAS::RecPhoton& photon)
{
    return s << photon.p4 << ' ' << photon.scales.size() << ' ' << photon.weights.size();
}
 
#if defined(__ROOTCLING__)
#pragma link C++ class DAS::GenPhoton +;
#pragma link C++ class std::vector<DAS::GenPhoton> +;
 
#pragma link C++ class DAS::RecPhoton +;
#pragma link C++ class std::vector<DAS::RecPhoton> +;
#endif