Event.h

Go to the documentation of this file.

#pragma once
 
// STD
#include <vector>
 
// Darwin
#include <GenericObject.h>
#include <Variation.h>
#include <Weight.h>
 
// ROOT
#include <TRandom3.h>
#include <Math/Point3D.h>
 
namespace Darwin::Physics {
 
class GenEvent : public AnonymousWeightableObject {
 
    std::string_view WeightGroup () const override { return WeightVar; }
 
public:
    static inline const char * const WeightVar = "GenEvtWgts"; //< Name of gen weights in variations
 
    GenEvent () = default;
    virtual ~GenEvent () = default;
};
 
class RecEvent : public AnonymousWeightableObject {
 
    std::string_view WeightGroup () const override { return WeightVar; }
 
public:
    static inline const char * const WeightVar = "RecEvtWgts"; //< Name of rec weights in variations
 
    int fill = 0,  
        runNo = 0, 
        lumi = 0;  
    unsigned long long evtNo = 0; 
 
    RecEvent () = default;
    virtual ~RecEvent () = default;
 
    void Clear () override 
    {
        fill = 0;
        runNo = 0;
        lumi = 0;
        evtNo = 0;
        AnonymousWeightableObject::Clear();
    }
};
 
class GenPV : public AnonymousWeightableObject {
 
    std::string_view WeightGroup () const override { return WeightVar; }
 
public:
    static inline const char * const WeightVar = "GenVtxWgts"; //< Name of gen weights in variations
 
    ROOT::Math::RhoZPhiPointF location = ROOT::Math::RhoZPhiPointF(std::nan("rho"), std::nan("z"), std::nan("phi"));
 
    GenPV () = default;
    virtual ~GenPV () = default;
 
    void Clear () override
    {
        AnonymousWeightableObject::Clear();
        location.SetCoordinates(std::nan("rho"), std::nan("z"), std::nan("phi"));
    }
};
 
class RecPV : public GenPV { 
 
    std::string_view WeightGroup () const override { return WeightVar; }
 
public:
    static inline const char * const WeightVar = "RecVtxWgts"; //< Name of gen weights in variations
 
    float chi2 = std::nan("chi2"); 
    int ndf = -1; 
 
    RecPV () = default;
    virtual ~RecPV () = default;
 
    void Clear () override
    {
        GenPV::Clear();
        chi2 = std::nan("chi2");
        ndf = -1;
    }
};
 
} // end of Darwin::Physics
 
inline std::ostream& operator<< (std::ostream& s, const Darwin::Physics::GenPV& pv)
{
    return s << '(' << pv.location.rho() << ','
                    << pv.location.z  () << ','
                    << pv.location.phi() << ')';
}
 
inline std::ostream& operator<< (std::ostream& s, const Darwin::Physics::RecEvent& event)
{
    return s << event.fill << ' ' << event.runNo << ' ' << event.lumi << ' ' << event.evtNo;
}