SusyNtuple/Event.h

// Dear emacs, this is -*- c++ -*-
#ifndef SUSYNTUPLE_EVENT_H
#define SUSYNTUPLE_EVENT_H

#include "SusyNtuple/SusyNtSys.h"
#include "SusyNtuple/SusyDefs.h" // DataStream

#include "TBits.h"
#include "TObject.h"

namespace Susy
{

class Event: public TObject
{
public:
    Event() :
        trigBits(m_nTriggerBits)
        {
            clear();
        }
    virtual ~Event(){};

    unsigned int run;   
    unsigned int eventNumber; 
    unsigned int lb;          
    DataStream stream;        

    bool isMC;                
    unsigned int mcChannel;   
    float w;                  
    uint64_t initialNumberOfEvents;  
    double sumOfEventWeights;        
    double sumOfEventWeightsSquared; 

    unsigned int larError;    

    unsigned int nVtx;        
    float avgMu;              

    int hfor;                 

    int susyFinalState;       
    int hDecay;               
    bool eventWithSusyProp;   

    int susySpartId1;         
    int susySpartId2;         

    float mllMcTruth;         
    bool passMllForAlpgen;    

    //unsigned int trigFlags; ///< Event level trigger bits
    long long trigFlags;      

    TBits               trigBits;
    static const size_t m_nTriggerBits=64;


    bool passTrig(long long mask, bool requireAll=true) const {
      if(requireAll) return (trigFlags & mask) == mask;
      else return mask == 0 || (trigFlags & mask) != 0;
    }

    // Event Flag to check for LAr, bad jet, etc. List found in SusyDefs.h under EventCheck
    // Use cutFlags instead
    //int evtFlag[NtSys_N];

    unsigned int cutFlags[NtSys::SYS_UNKNOWN];

    // Reweighting and scaling
    float wPileup;            
    float wPileup_up;         
    float wPileup_dn;         
    float wPileupAB3;         
    float wPileupAB;          
    float wPileupIL;          
    float wPileupAE;          
    float xsec;               
    float errXsec;            
    float sumw;               
    float pdfSF;              

    int pdf_id1;
    int pdf_id2;
    float pdf_x1;
    float pdf_x2;
    float pdf_scale;

    float eventScale;
    float alphaQCD;
    float alphaQED;

    void print() const;

    void clear(){
      run = eventNumber = lb = 0;
      stream = Stream_Unknown;
      isMC = false;
      mcChannel = w = 0;
      initialNumberOfEvents = sumOfEventWeights = sumOfEventWeightsSquared = 0;
      larError = 0;
      nVtx = avgMu = trigFlags = 0;
      hfor = -1;
      susyFinalState = 0;
      mllMcTruth = -1.0;
      hDecay = -1;
      eventWithSusyProp = false;
      susySpartId1 = susySpartId2 = 0;
      passMllForAlpgen = true;
      //memset(evtFlag,0,sizeof(evtFlag));
      memset(cutFlags,0,sizeof(cutFlags));
      wPileup = wPileup_up = wPileup_dn = wPileupAB3 = wPileupAB = wPileupIL = wPileupAE = 0;
      xsec = errXsec = sumw = pdfSF = 0;
      pdf_id1 = pdf_id2 = pdf_x1 = pdf_x2 = pdf_scale = 0;
      eventScale = alphaQCD = alphaQED = 0;
    }

    ClassDef(Event, 31);
  };
} // Susy
#endif