Event.h

Go to the documentation of this file.

#pragma once
 
#include <vector>
 
#include "Core/Objects/interface/Variation.h"
#include "Core/Objects/interface/Weight.h"
 
#include <TRandom3.h>
 
namespace DAS {
 
class AbstractEvent {
protected:
    virtual std::string_view weight_group () const = 0;
 
    AbstractEvent () = default;
    virtual ~AbstractEvent () = default;
 
public:
    Weights weights; 
 
    inline double Weight (const Uncertainties::Variation& v) const 
    {
        if (v.group == weight_group()) {
            const DAS::Weight& w = weights.at(v.index);
            if (w.i == v.bit) return w;
        }
        return weights.front();
    }
};
 
struct GenEvent : public AbstractEvent {
    static inline const char * const WeightVar = "GenEvtWgts"; //< Name of gen weights in variations
 
    float hard_scale; 
 
    void clear (); 
 
private:
    std::string_view weight_group () const override { return WeightVar; }
};
 
struct RecEvent : public AbstractEvent {
    static inline const char * const WeightVar = "RecEvtWgts"; //< Name of rec weights in variations
 
    int runNo, 
        lumi;  
    unsigned long long evtNo; 
 
    void clear (); 
 
private:
    std::string_view weight_group () const override { return WeightVar; }
};
 
struct Trigger {
    std::vector<bool> Bit; 
    std::vector<int> PreHLT,   
                     PreL1min, 
                     PreL1max; 
 
    void clear (); 
};
 
struct MET {
    float Et,    
          SumEt, 
          Pt,    
          Phi;   
    std::vector<bool> Bit; 
 
    void clear (); 
};
 
struct PileUp {
 
    static float MBxsec, MBxsecRelUnc;
 
    float rho;      
    int nVtx;       
    float trpu;     
    int intpu;      
    std::vector<float> pthats; 
    std::vector<unsigned long long> MBevents; 
 
    float GetTrPU (const char v = '0') const; 
    float GetInTimePU (const char v = '0', const bool useDefault = true) const; 
    static TRandom3 r3;
    void clear (); 
};
 
struct PrimaryVertex {
    float Rho,  
          z,    
          chi2, 
          ndof; 
    bool fake; 
 
    float distGen; 
 
    void clear (); 
};
 
}
 
#if defined(__ROOTCLING__)
#pragma link C++ class DAS::AbstractEvent +;
#pragma link C++ class DAS::GenEvent +;
#pragma link C++ class DAS::RecEvent +;
 
#pragma link C++ class DAS::Trigger +;
 
#pragma link C++ class DAS::PrimaryVertex +;
 
#pragma link C++ class DAS::PileUp +;
 
#pragma link C++ class DAS::MET +;
 
#endif