Photon.h

Go to the documentation of this file.

#pragma once
 
#include "Core/Objects/interface/PhysicsObject.h"
 
#include <TString.h>
 
#include <array>
#include <cstdint>
 
namespace DAS {
 
struct GenPhoton : public PhysicsObject {
    static inline const char * ScaleVar = "GenPhotonScales",   //< Name of photon rec scales in variations
                             * WeightVar = "GenPhotonWgts"; //< Name of photon rec weights in variations
 
    bool zAncestor = false, 
         prompt = false; 
 
    GenPhoton() = default; 
 
private:
    std::string_view scale_group () const final { return ScaleVar; }
    std::string_view weight_group () const final { return WeightVar; }
};
 
struct RecPhoton {
    static inline const char * ScaleVar = "RecPhotonScales",   //< Name of photon rec scales in variations
                             * WeightVar = "RecPhotonWgts"; //< Name of photon rec weights in variations
 
    enum EnergyVariation {
        Nominal = 0,    
        ScaleUp,        
        ScaleDown,      
        SigmaPhiUp,     
        SigmaPhiDown,   
        SigmaRhoUp,     
        SigmaRhoDown,   
        VariationsCount 
    };
    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, 
    };
 
    FourVector p4; 
    Weights weights = {{1.,0}}; 
 
    float scEta, 
          sigmaIEtaIEta, 
          hOverE, 
          chargedIsolation, 
          neutralHadronIsolation, 
          photonIsolation, 
          worstChargedIsolation; 
 
    std::array<float, VariationsCount> scales;
    float ecalEnergyErrPostCorr = -1; 
 
    std::uint32_t selectors = 0; 
 
    RecPhoton () = default; 
    ~RecPhoton () = default; 
 
    inline float Rapidity () const { return       p4.Rapidity() ; };
    inline float AbsRap   () const { return std::abs(Rapidity()); };
 
    inline FourVector CorrP4 (EnergyVariation var = Nominal) const { return p4 * scales[var]; }
};
 
} // 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