o First Version of TRDnSigma implementation (Xianguo) o still requires some catching...
[u/mrichter/AliRoot.git] / PWGLF / totEt / AliAnalysisEtCuts.cxx
1 //_________________________________________________________________________
2 //  Utility Class for transverse energy studies
3 //  Base class for ESD & MC analysis
4 //  - cuts for reconstruction and MonteCarlo 
5 // implementation file
6 //  
7 //*-- Authors: Oystein Djuvsland (Bergen), David Silvermyr (ORNL)
8 //_________________________________________________________________________
9
10 #include "AliAnalysisEtCuts.h"
11
12 #include <iostream>
13 #include "TNamed.h"
14
15 using namespace std;
16
17 ClassImp(AliAnalysisEtCuts);
18
19
20 AliAnalysisEtCuts::AliAnalysisEtCuts() : 
21   TNamed()
22                                     //
23   ,fCommonEtaCut(0.5)
24   ,fCommonClusterEnergyCut(0.15)
25   ,fCommonTrackPtCut(0.0)
26   ,fCommonSingleCell(1)
27   ,fEmcalTrackDistanceCut(15.0)
28   ,fEmcalTrackDxCut(5.0)
29   ,fEmcalTrackDzCut(5.0)
30   ,fPhosTrackDistanceCut(10.0)  
31   ,fPhosTrackDxCut(8.0)
32   ,fPhosTrackDzCut(3.0)
33   ,fPhosTrackRCut(5.0)
34   ,fPhosBadDistanceCut(3.0)
35   
36   ,fGeometryPhosEtaAccCut(0.12)
37   ,fGeometryPhosPhiAccMinCut(-100)//260.0)
38   ,fGeometryPhosPhiAccMaxCut(-40)//320.0)
39   ,fGeometryPhosDetectorRadius(460.0)
40                                     //
41   ,fGeometryEmcalEtaAccCut(0.6)
42   ,fGeometryEmcalPhiAccMinCut(80.0) // rad 1.4
43   ,fGeometryEmcalPhiAccMaxCut(120.0) // rad 2.1
44                                     //,fGeometryEmcalPhiAccMaxCut(180.0) // rad 3.14
45   ,fGeometryEmcalDetectorRadius(440.0)
46                                     //
47   ,fReconstructedVertexXCut(0.5)
48   ,fReconstructedVertexYCut(0.5)
49   ,fReconstructedVertexZCut(12.0)
50   ,fReconstructedIPxyCut(1.5)
51   ,fReconstructedIPzCut(1.5)
52   ,fReconstructedNTpcClustersCut(30)
53   ,fReconstructedNItsClustersCut(3)
54   ,fReconstructedPidCut(0.0)
55                                     //
56   ,fReconstructedPhosClusterType(-1)
57   ,fReconstructedPhosClusterEnergyCut(0.25)// GeV
58   ,fReconstructedPhosSingleCellEnergyCut(0.5)
59   ,fReconstructedPhosTrackDistanceTightCut(3.0)
60   ,fReconstructedPhosTrackDistanceMediumCut(5.0)
61   ,fReconstructedPhosTrackDistanceLooseCut(15.0)
62                                     //
63   ,fReconstructedEmcalClusterType(1)
64   ,fReconstructedEmcalClusterEnergyCut(0.30) // GeV
65   ,fReconstructedEmcalSingleCellEnergyCut(0.5)
66   ,fReconstructedEmcalTrackDistanceTightCut(5.0)
67   ,fReconstructedEmcalTrackDistanceMediumCut(10.0)
68   ,fReconstructedEmcalTrackDistanceLooseCut(15.0)
69   
70   ,fMonteCarloSingleChargedParticle(3)
71   ,fMonteCarloNeutralParticle(0)
72
73   ,fHistMakeTree(kFALSE)
74   ,fHistMakeTreeDeposit(kFALSE)
75   ,fHistNbinsMult(2000)
76   ,fHistMinMult(-0.5)
77   ,fHistMaxMult(1999.5)
78   ,fHistNbinsTotEt(10000)
79   ,fHistMinTotEt(0.000)
80   ,fHistMaxTotEt(1000)
81   ,fHistNbinsParticleEt(5000)
82   ,fHistMinParticleEt(0)
83   ,fHistMaxParticleEt(500)
84   ,fHistNbinsParticlePt(200) 
85   ,fHistMinParticlePt(0)
86   ,fHistMaxParticlePt(20)
87   
88   ,fPrimaryVertexCutXY(4.0)
89   ,fPrimaryVertexCutZ(20.0)
90 { // ctor
91 }
92
93 AliAnalysisEtCuts::~AliAnalysisEtCuts()
94 { // dtor
95 }
96
97 void AliAnalysisEtCuts::SetPbPbDefaults()
98 { // just increase seom histogram max values for now
99   // enough to multiply conservative p+p defaults by a factor 100?
100   fHistMaxMult = 20000;
101   fHistMaxTotEt = 10000;
102 }
103
104