SFHMarEvent.h

Go to the documentation of this file.

// Class SFHMarEvent
//
// Author     : J. List
// Created    : 06/2005
// Last update: $Date: 2005/07/28 09:14:27 $
//          by: $Author: boehmej $
// Comment: Contains definition of event class

#ifndef __SFHMAREVENT
#define __SFHMAREVENT


class NTFloatFun;
class NTBaseCut;
class FloatFun;
class BaseCut;

#include "Marana/TMarEvent.h"

class SFHMarEvent :  public TMarEvent
{
public :
//... which particles to define 
  NTFloatFun& fffind_ElePho;
  NTFloatFun& fffind_Muon;
  NTFloatFun& fffind_Jet;
  NTFloatFun& fffind_NP;

// //--- pointer to Mods Particles
//   H1ExclHfsIterator ModsPartHFS;
// 
// //...mods particles after calibrations
// 
//   vector<TMarBody> PartCalEm;
//   vector<TMarBody> PartCalMuon;
//   vector<TMarBody> PartCalJet;
//   TLorentzVector HadCalJetVec;  
//   TLorentzVector HadCalNoJetVec;
// 
// //...mods particles after calibrations and systematics
//   vector<TMarBody>  PartScaled[3];
//   vector<TMarBody>* PartEm;
//   vector<TMarBody>* PartMuon;
//   vector<TMarBody>* PartJet;
//   TLorentzVector HadNoJetVec;
//   TLorentzVector HadJetVec;
//   TLorentzVector TotalVec;
// 
// //...List of selected particles
// 
//   TClonesArray *electron;
//   TClonesArray *photon;
//   TClonesArray *muon;
//   TClonesArray *jet;
//   TClonesArray *np;
// 
// 
// //... List of generated particles
//   TClonesArray *GenElectron;
//   TClonesArray *GenPhoton;
//   TClonesArray *GenMuon;
//   TClonesArray *GenTau;
//   TClonesArray *GenNp;
//   TClonesArray *GenHfs;
//   TLorentzVector GenScatEle;
// 
  NTFloatFun& ffGenVtxZ;
  NTFloatFun& ffyhgen;
  NTFloatFun& ffGenRad; // flag pour les photons rad
//  
// //... to cut MC on generated quantities
//   list<TMarCutMCList> genCut;
// 
//... type of systematic to apply (see TMarSyst)
  NTFloatFun& ffSyst_Type;
  NTFloatFun& ffSyst_Num;
  NTFloatFun& ffSyst_Sign;
//  
// //... class for the calibration of jets 
//   H1JetCalibration* JetCalib;
// //... class for the calibration of electrons 
//   TMarH1ElecCalibration* ElecCalib;
//   H1ElecCalibration* H1ElecCalib;
//   
// //... Type and year of data
//   
//   enum Type {Data=0,  NC_Djo, NC_Rap, CC_Djo, EE_Grape,  MM_Grape,
//               TT_Grape, EG_Wabgen, GP_Pythia, W_Epvec, QCD_Ins, Htm, Anotop, PsCC, UNKNOWN } RunType;
//   enum Year {Y9497 = 3, Y9899 = 4 , Y9900 = 5, Y0203 = 6, Y0304 = 7, Y05 = 8} RunYear;
//   
// //... lumi object 
//   TMarLumi* lumi;    
  NTFloatFun& ffpolar;
  
//... general quantities 
  NTFloatFun& ffVertex_x; 
  NTFloatFun& ffVertex_y; 
  NTFloatFun& ffVertex_z; 
  NTFloatFun& ffEpz;         //total E-Pz    
  NTFloatFun& ffPtmiss;      //total Ptmiss
  NTFloatFun& ffEt;

//... scattered electron
  NTFloatFun& ffIndexElScat;//scattered electron =higest Pt electron if filled from MODS 
  NTFloatFun& ffPte;        // == caracteristics of scattered electron, independant of TMARBody electronR
  NTFloatFun& ffEe;
  NTFloatFun& ffPze;
  NTFloatFun& ffThe;
  NTFloatFun& ffPhe;
  NTFloatFun& ffQ2e;
  NTFloatFun& ffye;
  NTFloatFun& ffxe;
  

//... hfs quantities  
  NTFloatFun& ffEtotCalo;  // energy sums in LAr + Spacal

  NTFloatFun& ffPth;    // Pt of all hadronic final state
  NTFloatFun& ffThh;
  NTFloatFun& ffPhh; 
  NTFloatFun& ffEh;
  NTFloatFun& ffPzh;
  NTFloatFun& ffQ2h;
  NTFloatFun& ffyh;
  NTFloatFun& ffxh;
  NTFloatFun& ffGammah;

//... Calorimeter quantities ....
  NTFloatFun& ffptIronfrac;
  NTFloatFun& ffEIronfrac;
  NTFloatFun& ffptSpacalfrac;
  NTFloatFun& ffESpacalfrac;


//... Kinematic variables (recalculated from all Particles candidate on Mods)
  NTFloatFun& ffPtda;  //da method
  NTFloatFun& ffQ2da;
  NTFloatFun& ffyda;
  NTFloatFun& ffxda;
  NTFloatFun& ffEda;

  NTFloatFun& ffPts;  //sigma method
  NTFloatFun& ffQ2s;
  NTFloatFun& ffys;
  NTFloatFun& ffxs;

  NTFloatFun& ffQ2es;  //electron-sigma method
  NTFloatFun& ffxes;
  NTFloatFun& ffyes;


//... Timing variables
  NTFloatFun& fftCJC;  // timing CJC  
  NTFloatFun& fftLAr;  // timing tLAr

//... Weigths to apply to events
  NTFloatFun& ffMCWeight; // total MC Weight
  NTFloatFun& ffLumiWeight; // Weight to sum MC
  NTFloatFun& fffTriggerWeight; // Weight from triggers
  NTFloatFun& fffQCDWeight; // Weight from QCD (pdf's) or LO->NLO MC reweigthing

//... Vap/Vp calculation
  NTFloatFun& ffVparl;
  NTFloatFun& ffVanti;

//... Triggers arrays

//... cuts for the definition of selected particles (TMarBody)
//...Cuts for electron definition  
  NTFloatFun& ffPtMinEle;
  NTFloatFun& ffThetaMinEle;  
  NTFloatFun& ffThetaMaxEle;  
  NTFloatFun& ffConeEle;

//...Cuts for muon definition 
  NTFloatFun& ffPtMinMu;
  NTFloatFun& ffThetaMinMu;  
  NTFloatFun& ffThetaMaxMu;  
  NTFloatFun& ffDTrTrMinMu;

//...Cuts for jets definition 
  NTFloatFun& ffPtMinJet;
  NTFloatFun& ffThetaMinJet;  
  NTFloatFun& ffThetaMaxJet;
  NTFloatFun& ffRMinJet;
  NTFloatFun& ffEmFracJet;

//...Cuts for Neutral particle definition 
  NTFloatFun& ffPtMinNp;
  NTFloatFun& ffThetaMinNp;  
  NTFloatFun& ffThetaMaxNp;

// //... Rejection handler: 
//    list<string> RejectReason;
// 
// //... Trigger efficiencies class
//   TTriggerReweight* fTrigReweight;
//   
// //... MC reweight class
//   TMCXSectionReweight* fMCReweight;  


// //...from Marana/TTools.h 
//   // Is Phi Crack
//   NTBaseCut& bcIsPhiCrack(Float_t Phi);
//   NTBaseCut& bcIsLArCrack(Float_t z,Float_t phi);
//   NTBaseCut& bcIsLArCrack(Float_t z,Float_t zBBE,Float_t phi);
// 
//   // Is z Crack
//   NTBaseCut& bcIsZCrack(Float_t z);
// 
// //--- Is timing ok
//   NTBaseCut& bcIsTimingOk();
//   NTBaseCut& bcIs2004CalibratedTimingOK(void);
// 

// ... MY stuff
  FloatFun& ffNSelJets;
  FloatFun& ffPtSelJet1;
  FloatFun& ffThSelJet1;
  
  BaseCut& bcTimingOK;

//... Standard constructor and destructor
  SFHMarEvent(const TRun *run);
  virtual ~SFHMarEvent();

//  ClassDef(SFHMarEvent,1) // Make class known to root 

};

#endif

---