1 //_________________________________________________________________________
2 // Utility Class for transverse energy studies
3 // Base class for ESD & MC analysis
4 // - cuts for reconstruction and MonteCarlo
7 //*-- Authors: Oystein Djuvsland (Bergen), David Silvermyr (ORNL)
8 //_________________________________________________________________________
10 #include "AliAnalysisEtCuts.h"
17 ClassImp(AliAnalysisEtCuts);
20 AliAnalysisEtCuts::AliAnalysisEtCuts() :
24 ,fCommonClusterEnergyCut(0.0)
25 ,fCommonTrackPtCut(0.0)
28 ,fGeometryPhosEtaAccCut(0.12)
29 ,fGeometryPhosPhiAccMinCut(260.0)
30 ,fGeometryPhosPhiAccMaxCut(320.0)
31 ,fGeometryPhosDetectorRadius(460.0)
33 ,fGeometryEmcalEtaAccCut(0.7)
34 ,fGeometryEmcalPhiAccMinCut(80.0) // rad 1.4
35 ,fGeometryEmcalPhiAccMaxCut(120.0) // rad 2.1
36 ,fGeometryEmcalDetectorRadius(440.0)
38 ,fReconstructedVertexXCut(0.5)
39 ,fReconstructedVertexYCut(0.5)
40 ,fReconstructedVertexZCut(12.0)
41 ,fReconstructedIPxyCut(1.5)
42 ,fReconstructedIPzCut(1.5)
43 ,fReconstructedNTpcClustersCut(30)
44 ,fReconstructedNItsClustersCut(3)
45 ,fReconstructedPidCut(0.0)
47 ,fReconstructedPhosClusterType(-1)
48 ,fReconstructedPhosClusterEnergyCut(0.0)
49 ,fReconstructedPhosSingleCellEnergyCut(0.5)
50 ,fReconstructedPhosTrackDistanceCut(15.0)
52 ,fReconstructedEmcalClusterType(1)
53 ,fReconstructedEmcalClusterEnergyCut(0.1) // GeV
54 ,fReconstructedEmcalSingleCellEnergyCut(0.5)
55 ,fReconstructedEmcalTrackDistanceCut(15.0)
57 ,fMonteCarloSingleChargedParticle(3)
58 ,fMonteCarloNeutralParticle(0)
60 ,fHistMakeTree(kFALSE)
61 ,fHistMakeTreeDeposit(kTRUE)
65 ,fHistNbinsTotEt(1000)
66 ,fHistMinTotEt(0.0001)
68 ,fHistNbinsParticleEt(1000)
69 ,fHistMinParticleEt(0)
70 ,fHistMaxParticleEt(10)
71 ,fHistNbinsParticlePt(200)
72 ,fHistMinParticlePt(0)
73 ,fHistMaxParticlePt(20)
77 AliAnalysisEtCuts::~AliAnalysisEtCuts()
81 void AliAnalysisEtCuts::SetPbPbDefaults()
82 { // just increase seom histogram max values for now
83 // enough to multiply conservative p+p defaults by a factor 100?
85 fHistMaxTotEt = 10000;