TMarEvent.h

Go to the documentation of this file.

// Class TMarEvent
//
// Author     : Emmanuel Sauvan/Matti Peez/F. Cassol Brunner
// Created    : 02/2002
// Last update: $Date: 2005/06/30 17:08:00 $
//          by: $Author: sauvan $
// Comment: Contains definition of event class

#ifndef __TMAREVENT
#define __TMAREVENT

#include <string>
#include <vector>

// Include ROOT libs
#include <TTree.h>
#include <TBenchmark.h>
#include <TLorentzVector.h>


//------ H1 OO includes ------------
// Access general oo objects
#include "H1Skeleton/H1Tree.h"
#include "H1Skeleton/H1EventList.h"
#include "H1Skeleton/H1PointerManager.h"
#include "H1Pointers/H1ObjectPtr.h"
#include "H1Pointers/H1FloatPtr.h"
#include "H1Pointers/H1IntPtr.h"
#include "H1Pointers/H1BytePtr.h"
#include "H1JetCalibration/H1JetCalibration.h"
#include "H1HfsCalibration/H1HfsCalibEmHad.h"
#include "Marana/TMarH1ElecCalibration.h" 
#include "Marana/TTriggerReweight.h"
#include "Marana/TMCXSectionReweight.h"
#include "H1ElecCalibration/H1ElecCalibration.h"

// Access ODS
  

#include "H1OOBanks/H1BankPtr.h"


// Access mods
#include "H1Mods/H1PartCand.h"
#include "H1Mods/H1PartCandArrayPtr.h"
#include "H1Mods/H1PartEmArrayPtr.h"
#include "H1Mods/H1PartEm.h"
#include "H1Mods/H1PartMuonArrayPtr.h"
#include "H1Mods/H1PartMuon.h"
#include "H1Mods/H1PartJetArrayPtr.h"
#include "H1Mods/H1PartJet.h"
#include "H1Mods/H1PartSelTrackArrayPtr.h"
#include "H1Mods/H1PartSelTrack.h"
#include "H1Mods/H1PartCandArrayPtr.h"
#include "H1Mods/H1PartCand.h"


#include "H1Mods/H1ExclHfsIterator.h"
#include "H1Mods/H1PartMCArrayPtr.h"
#include "H1Mods/H1PartMC.h"


// Access tracks
#include "H1Tracks/H1VertexFittedTrack.h"
#include "H1Tracks/H1Vertex.h"
#include "H1Tracks/H1Track.h"
#include "H1Tracks/H1ReconstructedVertex.h"
#include "H1Steering/h1debug.h"

// Include Body definition
#include "Marana/TMarBody.h"


#include "Marana/TTools.h"
#include "Marana/TMarLumi.h"
#include "Marana/TRunHisto.h"
#include "Marana/TMarCutMC.h"
#include "Marana/TMarSyst.h" 

// Rad => Deg 
#define MYR2D 180.0/3.14159265359

class TMarEvent :  public TTools
{
//     class for useful event variables 
public :

//... which particles to define 
  Bool_t find_ElePho;
  Bool_t find_Muon;
  Bool_t find_Jet;
  Bool_t find_NP;



//--- pointer to Mods Particles
  H1PartCandArrayPtr ModsPartCand; 
  H1PartEmArrayPtr   ModsPartEm;
  H1PartMuonArrayPtr ModsPartMuon;
  H1ExclHfsIterator ModsPartHFS;
  H1PartJetArrayPtr  ModsPartJet;
  H1PartSelTrackArrayPtr ModsPartSelTrack;
  H1PartMCArrayPtr ModsPartMC;

//... Access user tree
  H1PartJetArrayPtr ModsGenKtJets;
  H1PartSelTrackArrayPtr ModsDTRATracks;
  H1PartSelTrackArrayPtr ModsDTNVTracks;

//...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;

  Float_t GenVtxZ;
  Double_t yhgen;
  Short_t  GenRad; // flag pour les photons rad
 
//... to cut MC on generated quantities
  list<TMarCutMCList> genCut;

//... type of systematic to apply (see TMarSyst)
  Int_t Syst_Type;
  Short_t Syst_Num;
  Short_t Syst_Sign;
 
//... class for the calibration of jets 
  H1JetCalibration* JetCalib;
//... class for the calibration of hfs (low Q2, h1oo test module) 
  H1HfsCalibEmHad* HfsCalib;
//... 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, POUDS, RAPD_UDS, UNKNOWN } RunType;
  enum Year {Y9497 = 3, Y9899 = 4 , Y9900 = 5, Y0203 = 6, Y0304 = 7, Y05 = 8} RunYear;
  
//... lumi object 
  TMarLumi* lumi;    
  Double_t polar;
  
//... general quantities 
  Double_t Vertex[3]; 
  Double_t Epz;         //total E-Pz    
  Double_t Ptmiss;      //total Ptmiss
  Double_t Et;

//... scattered electron

  Int_t IndexElScat;  //scattered electron =higest Pt electron if filled from MODS 
  Double_t Pte;   // == caracteristics of scattered electron, independant of TMARBody electronR
  Double_t Ee;
  Double_t Pze;
  Double_t The;
  Double_t Phe;
  Double_t Q2e;
  Double_t ye;
  Double_t xe;
  

//... hfs quantities  
  Double_t EtotCalo;  // energy sums in LAr + Spacal

  Double_t Pth;    // Pt of all hadronic final state
  Double_t Thh;
  Double_t Phh; 
  Double_t Eh;
  Double_t Pzh;
  Double_t Q2h;
  Double_t yh;
  Double_t xh;
  Double_t Gammah;

//... Calorimeter quantities ....
  Double_t ptIronfrac;
  Double_t EIronfrac;
  Double_t ptSpacalfrac;
  Double_t ESpacalfrac;


//... Kinematic variables (recalculated from all Particles candidate on Mods)
  Double_t Ptda;  //da method
  Double_t Q2da;
  Double_t yda;
  Double_t xda;
  Double_t Eda;

  Double_t Pts;  //sigma method
  Double_t Q2s;
  Double_t ys;
  Double_t xs;

  Double_t Q2es;  //electron-sigma method
  Double_t xes;
  Double_t yes;


//... Timing variables
  Double_t tCJC;  // timing CJC  
  Double_t tLAr;  // timing tLAr

//... Weigths to apply to events
  Double_t MCWeight; // total MC Weight
  Double_t LumiWeight; // Weight to sum MC
  Double_t fTriggerWeight; // Weight from triggers
  Double_t fQCDWeight; // Weight from QCD (pdf's) or LO->NLO MC reweigthing

//... Vap/Vp calculation
  Float_t Vparl;
  Float_t Vanti;

//... Triggers arrays
  Byte_t Il1ac[128]; // array of 128 l1 actual subtrigger numbers
  Byte_t Il1rw[128]; // array of 128 l1 raw subtrigger numbers
  Byte_t Il1te[192]; // array of 192 l1 trigger elements


//... cuts for the definition of selected particles (TMarBody)
//...Cuts for electron definition  
  Float_t PtMinEle;
  Float_t ThetaMinEle;  
  Float_t ThetaMaxEle;  
  Float_t ConeEle;

//...Cuts for muon definition 
  Float_t PtMinMu;
  Float_t ThetaMinMu;  
  Float_t ThetaMaxMu;  
  Float_t DTrTrMinMu;

//...Cuts for jets definition 
  Float_t PtMinJet;
  Float_t ThetaMinJet;  
  Float_t ThetaMaxJet;
  Float_t RMinJet;
  Float_t EmFracJet;

//...Cuts for Neutral particle definition 
  Float_t PtMinNp;
  Float_t ThetaMinNp;  
  Float_t ThetaMaxNp;


//... Rejection handler: 
   list<string> RejectReason;

//... Trigger efficiencies class
  TTriggerReweight* fTrigReweight;
  
//... MC reweight class
  TMCXSectionReweight* fMCReweight;  

//... Clear up lists
  virtual void Clear(void);
  virtual void ClearSelectedParticles(void);


//... Calibrate/copy mods particles
  virtual void DefineInputParticles();

//... Apply the systematic shifts 
  virtual void ApplySystematicShifts(TMarSyst &);

//... Calculate all Derived Variables after caliabration ans scaling
  virtual void CalculateDerivedVariables();

//... Fill event information
  virtual Bool_t FillFromMods(){return kTRUE;};
  virtual Bool_t FillFromHat();
  virtual Bool_t FillGeneratedParticlesAndCut();
  virtual Bool_t MCGenSelection(TMarCutMCList& List);

//... initialize the event
  virtual Bool_t Init();  

//... next event 
  virtual Bool_t Next(TRunHisto& run);  

//... Find ele/pho 
  virtual Int_t FindElePho();

//... Find muon
  virtual Int_t FindMuon();
 
//... Find np
  virtual Int_t FindNP();

//... Find jet
  virtual Int_t FindJet();

//... retrieve generated particles
  virtual void FindGeneratedParticles();


//... apply hfs calibrations
  virtual void ApplyHfsCalibration();
 virtual void ApplyLowQ2HfsCalibration();
 
//...select particles
  virtual void SelectParticles();

//... return Cluster/Tracks Pth and Eh
  void GetPthEhTracksClusters(double&,double&,double&,double&);


//... print event informations
  virtual void PrintEventInfo();


//... Calulate min distance to jets
  Float_t MinDisttoJets(Float_t Theta,Float_t Phi);
  Float_t MinDisttoTracks(Float_t Theta,Float_t Phi);

//... matching tracks-clusters
  Float_t CalcDistTracks(H1PartSelTrackArrayPtr Tracks,const H1PartEm* EMParticle);
  Float_t CalcDistDTNV(const H1PartEm* EMParticle);
  void CalcDCA(TVector3 ClusStart, TVector3 TrackEnd,
        Float_t Theta, Float_t Phi,Float_t &Dca,Float_t &Alpha,Float_t &Dlon) const;


//... Some anti-background cut for CC
  virtual Float_t ApplyCCBackgroundCuts();


//... Functions for trigger efficiency
  virtual Float_t CCTriggerEfficiency();
  virtual Float_t JJTriggerEfficiency();
  virtual Float_t MuTriggerEfficiency();
  virtual Float_t CCAndLeptonTriggerEfficiency();
  virtual void ApplyTriggerEfficiencies();
  
//... for MC xsection reweighting
  virtual void ApplyMCXSectionReweight();

//... Add particle to list of preselected bodies
  Int_t AddParticleToList(TClonesArray *List,Int_t IndexCand,Int_t ID,Float_t E,Float_t Pt,Float_t Eta,Float_t Phi);
  Int_t AddParticleToList(TClonesArray *List,TMarBody *body);

//...sort particle list in Pt
  void SortPt(TClonesArray *list);

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

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

};

#endif

---