]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGGA/CaloTrackCorrelations/AliAnaParticleIsolation.h
Add centrality control histogram plus some cosmetics
[u/mrichter/AliRoot.git] / PWGGA / CaloTrackCorrelations / AliAnaParticleIsolation.h
CommitLineData
1a31a9ab 1#ifndef ALIANAPARTICLEISOLATION_H
2#define ALIANAPARTICLEISOLATION_H
3/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
1a31a9ab 5
6//_________________________________________________________________________
7
8// Class for the analysis of particle isolation
9// Input is selected particles put in AOD branch (AliAODPWG4ParticleCorrelation)
10//
11// Class created from old AliPHOSGammaJet
12// (see AliRoot versions previous Release 4-09)
13
14//-- Author: Gustavo Conesa (INFN-LNF)
15
16// --- ROOT system ---
17class TH2F;
18class TList ;
19class TObjString;
20
21// --- ANALYSIS system ---
745913ae 22#include "AliAnaCaloTrackCorrBaseClass.h"
1a31a9ab 23class AliAODPWG4Particle;
24class AliAODPWG4ParticleCorrelation ;
25
26
745913ae 27class AliAnaParticleIsolation : public AliAnaCaloTrackCorrBaseClass {
1a31a9ab 28
29 public:
30 AliAnaParticleIsolation() ; // default ctor
803d06a8 31 virtual ~AliAnaParticleIsolation() { ; } //virtual dtor
1a31a9ab 32
1a31a9ab 33 // Main general methods
803d06a8 34
1a31a9ab 35 TObjString * GetAnalysisCuts();
36
37 TList * GetCreateOutputObjects();
38
803d06a8 39 void InitParameters();
40
1a31a9ab 41 void MakeAnalysisFillAOD() ;
42
43 void MakeAnalysisFillHistograms() ;
44
45 void Print(const Option_t * opt)const;
46
1a31a9ab 47 //Analysis specific methods
803d06a8 48
1a31a9ab 49 void MakeSeveralICAnalysis(AliAODPWG4ParticleCorrelation * ph);
50
1a31a9ab 51 // Analysis Setters and Getters
52
53 TString GetCalorimeter() const { return fCalorimeter ; }
54 Int_t GetNCones() const { return fNCones ; }
55 Int_t GetNPtThresFrac() const { return fNPtThresFrac ; }
56 Float_t GetConeSizes(Int_t i) const { return fConeSizes[i] ; }
57 Float_t GetPtThresholds(Int_t i) const { return fPtThresholds[i] ; }
58 Float_t GetPtFractions(Int_t i) const { return fPtFractions[i] ; }
59
803d06a8 60 void SetCalorimeter(TString & det) { fCalorimeter = det ; }
61 void SetNCones(Int_t ncs) { fNCones = ncs ; }
62 void SetNPtThresFrac(Int_t npt) { fNPtThresFrac = npt ; }
63 void SetConeSizes(Int_t i, Float_t r) { fConeSizes[i] = r ; }
64 void SetPtThresholds(Int_t i, Float_t pt) { fPtThresholds[i] = pt ; }
65 void SetPtFractions(Int_t i, Float_t pt) { fPtFractions[i] = pt ; }
1a31a9ab 66
67 Bool_t IsReIsolationOn() const { return fReMakeIC ; }
09273901 68 void SwitchOnReIsolation() { fReMakeIC = kTRUE ; }
69 void SwitchOffReIsolation() { fReMakeIC = kFALSE ; }
1a31a9ab 70
71 Bool_t IsSeveralIsolationOn() const { return fMakeSeveralIC ; }
72 void SwitchOnSeveralIsolation() { fMakeSeveralIC = kTRUE ; }
73 void SwitchOffSeveralIsolation() { fMakeSeveralIC = kFALSE ; }
09273901 74
75 void SwitchOnTMHistoFill() { fFillTMHisto = kTRUE ; }
76 void SwitchOffTMHistoFill() { fFillTMHisto = kFALSE ; }
77
78 void SwitchOnSSHistoFill() { fFillSSHisto = kTRUE ; }
79 void SwitchOffSSHistoFill() { fFillSSHisto = kFALSE ; }
80
1a31a9ab 81 //Histogrammes setters and getters
82
83 virtual void SetHistoPtSumRangeAndNBins(Float_t min, Float_t max, Int_t n){
84 fHistoNPtSumBins = n ; fHistoPtSumMax = max ; fHistoPtSumMin = min ; }
85
86 Int_t GetHistoNPtSumBins() const { return fHistoNPtSumBins ; }
87 Float_t GetHistoPtSumMin() const { return fHistoPtSumMin ; }
88 Float_t GetHistoPtSumMax() const { return fHistoPtSumMax ; }
89
90 virtual void SetHistoPtInConeRangeAndNBins(Float_t min, Float_t max, Int_t n) {
91 fHistoNPtInConeBins = n ; fHistoPtInConeMax = max ; fHistoPtInConeMin = min ; }
92
803d06a8 93 Int_t GetHistoNPtInConeBins() const { return fHistoNPtInConeBins; }
94 Float_t GetHistoPtInConeMin() const { return fHistoPtInConeMin ; }
95 Float_t GetHistoPtInConeMax() const { return fHistoPtInConeMax ; }
1a31a9ab 96
97 private:
98
99 TString fCalorimeter ; // Calorimeter where neutral particles in cone for isolation are;
100 Bool_t fReMakeIC ; // Do isolation analysis
101 Bool_t fMakeSeveralIC ; // Do analysis for different IC
09273901 102 Bool_t fFillTMHisto; // Fill track matching plots
103 Bool_t fFillSSHisto; // Fill Shower shape plots
104
803d06a8 105 // Analysis data members for multiple cones and pt thresholds
106 Int_t fNCones ; //! Number of cone sizes to test
107 Int_t fNPtThresFrac ; //! Number of ptThres and ptFrac to test
108
109 Float_t fConeSizes[5] ; //! Array with cones to test
110 Float_t fPtThresholds[5] ; //! Array with pt thresholds to test
111 Float_t fPtFractions[5] ; //! Array with pt thresholds to test
112
113 TH1F* fhPtThresIsolated[5][5] ; //! Isolated particle with pt threshold
114 TH1F* fhPtFracIsolated[5][5] ; //! Isolated particle with pt threshold
115 TH2F* fhPtSumIsolated[5] ; //! Isolated particle with threshold on cone pt sum
1a31a9ab 116
117 //Histograms
118
0fb69ade 119 TH1F * fhEIso ; //! Number of isolated particles
1a31a9ab 120 TH1F * fhPtIso ; //! Number of isolated particles
121 TH2F * fhPhiIso ; //! Phi of isolated particles
122 TH2F * fhEtaIso ; //! eta of isolated particles
0fb69ade 123 TH2F * fhEtaPhiIso ; //! eta vs phi of isolated particles
1a31a9ab 124 TH1F * fhPtNoIso ; //! Number of not isolated leading particles
803d06a8 125 TH1F * fhPtDecayIso ; //! Number of isolated Pi0 decay particles (invariant mass tag)
126 TH1F * fhPtDecayNoIso ; //! Number of not isolated Pi0 decay leading particles (invariant mass tag)
1a31a9ab 127 TH2F * fhConeSumPt ; //! Sum Pt in the cone
128 TH2F * fhPtInCone ; //! Particle Pt in the cone
129 TH2F * fhFRConeSumPt ; //! Sum Pt in the forward region cone (phi +90)
130 TH2F * fhPtInFRCone ; //! Particle Pt in the forward region cone (phi +90 )
131
1a31a9ab 132
133 //MC
134 TH1F * fhPtIsoPrompt; //! Number of isolated prompt gamma
135 TH2F * fhPhiIsoPrompt; //! Phi of isolated prompt gamma
136 TH2F * fhEtaIsoPrompt; //! eta of isolated prompt gamma
137 TH1F * fhPtThresIsolatedPrompt[5][5]; //! Isolated prompt gamma with pt threshold
138 TH1F * fhPtFracIsolatedPrompt[5][5]; //! Isolated prompt gamma with pt frac
139 TH2F * fhPtSumIsolatedPrompt[5]; //! Isolated prompt gamma with threshold on cone pt sume
140 TH1F * fhPtIsoFragmentation; //! Number of isolated fragmentation gamma
141 TH2F * fhPhiIsoFragmentation; //! Phi of isolated fragmentation gamma
142 TH2F * fhEtaIsoFragmentation; //! eta of isolated fragmentation gamma
143 TH1F * fhPtThresIsolatedFragmentation[5][5]; //! Isolated fragmentation gamma with pt threshold
144 TH1F * fhPtFracIsolatedFragmentation[5][5]; //! Isolated fragmentation gamma with pt frac
145 TH2F * fhPtSumIsolatedFragmentation[5]; //! Isolated fragmentation gamma with threshold on cone pt sume
803d06a8 146 TH1F * fhPtIsoPi0Decay; //! Number of isolated pi0 decay gamma
147 TH2F * fhPhiIsoPi0Decay; //! Phi of isolated pi0 decay gamma
148 TH2F * fhEtaIsoPi0Decay; //! eta of isolated pi0 decay gamma
149 TH1F * fhPtThresIsolatedPi0Decay[5][5]; //! Isolated pi0 decay gamma with pt threshold
150 TH1F * fhPtFracIsolatedPi0Decay[5][5]; //! Isolated pi0 decay gamma with pt frac
151 TH2F * fhPtSumIsolatedPi0Decay[5]; //! Isolated pi0 decay gamma with threshold on cone pt sume
152 TH1F * fhPtIsoEtaDecay; //! Number of isolated eta decay gamma
153 TH2F * fhPhiIsoEtaDecay; //! Phi of isolated eta decay gamma
154 TH2F * fhEtaIsoEtaDecay; //! eta of isolated eta decay gamma
155 TH1F * fhPtThresIsolatedEtaDecay[5][5]; //! Isolated eta decay gamma with pt threshold
156 TH1F * fhPtFracIsolatedEtaDecay[5][5]; //! Isolated eta decay gamma with pt frac
157 TH2F * fhPtSumIsolatedEtaDecay[5]; //! Isolated eta fecay gamma with threshold on cone pt sume
158 TH1F * fhPtIsoOtherDecay; //! Number of isolated other decay gamma
159 TH2F * fhPhiIsoOtherDecay; //! Phi of isolated other decay gamma
160 TH2F * fhEtaIsoOtherDecay; //! eta of isolated other decay gamma
1a31a9ab 161 TH1F * fhPtThresIsolatedOtherDecay[5][5]; //! Isolated OtherDecay gamma with pt threshold
162 TH1F * fhPtFracIsolatedOtherDecay[5][5]; //! Isolated OtherDecay gamma with pt frac
163 TH2F * fhPtSumIsolatedOtherDecay[5]; //! Isolated OtherDecay gamma with threshold on cone pt sume
164 TH1F * fhPtIsoConversion; //! Number of isolated Conversion gamma
165 TH2F * fhPhiIsoConversion; //! Phi of isolated Conversion gamma
166 TH2F * fhEtaIsoConversion; //! eta of isolated Conversion gamma
167 TH1F * fhPtThresIsolatedConversion[5][5]; //! Isolated Conversion gamma with pt threshold
168 TH1F * fhPtFracIsolatedConversion[5][5]; //! Isolated Conversion gamma with pt frac
169 TH2F * fhPtSumIsolatedConversion[5]; //! Isolated Conversion gamma with threshold on cone pt sume
170 TH1F * fhPtIsoUnknown; //! Number of isolated Unknown
171 TH2F * fhPhiIsoUnknown; //! Phi of isolated Unknown
172 TH2F * fhEtaIsoUnknown; //! eta of isolated Unknown
173 TH1F * fhPtThresIsolatedUnknown[5][5]; //! Isolated Unknown gamma with pt threshold
174 TH1F * fhPtFracIsolatedUnknown[5][5]; //! Isolated Unknown gamma with pt frac
175 TH2F * fhPtSumIsolatedUnknown[5]; //! Isolated Unknown gamma with threshold on cone pt sume
176
803d06a8 177 TH1F * fhPtNoIsoPi0Decay; //! Number of not isolated leading pi0 decay gamma
178 TH1F * fhPtNoIsoEtaDecay; //! Number of not isolated leading eta decay gamma
179 TH1F * fhPtNoIsoOtherDecay; //! Number of not isolated leading other decay gamma
1a31a9ab 180 TH1F * fhPtNoIsoPrompt; //! Number of not isolated leading prompt gamma
181 TH1F * fhPtIsoMCPhoton; //! Number of isolated leading gamma
182 TH1F * fhPtNoIsoMCPhoton; //! Number of not isolated leading gamma
0fb69ade 183 TH1F * fhPtNoIsoConversion; //! Number of not isolated leading conversion gamma
184 TH1F * fhPtNoIsoFragmentation; //! Number of not isolated leading fragmentation gamma
185 TH1F * fhPtNoIsoUnknown; //! Number of not isolated leading hadrons
1a31a9ab 186
09273901 187 // Track matching studies
188 TH2F * fhTrackMatchedDEta ; //! Eta distance between track and cluster vs cluster E
189 TH2F * fhTrackMatchedDPhi ; //! Phi distance between track and cluster vs cluster E
190 TH2F * fhTrackMatchedDEtaDPhi ; //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV
31ae6d59 191 TH2F * fhdEdx ; //! matched track dEdx vs cluster E
192 TH2F * fhEOverP; //! matched track E cluster over P track vs cluster E, after dEdx cut
193 TH2F * fhTrackMatchedMCParticle; //! Trace origin of matched particle
194
09273901 195 // Shower Shape histograms
196 TH2F * fhELambda0; //! Shower shape of isolated photons
197 TH2F * fhELambda1; //! Shower shape of isolated photons
198
1a31a9ab 199 //Histograms settings
200 Int_t fHistoNPtSumBins; // Number of bins in PtSum histograms
201 Float_t fHistoPtSumMax; // PtSum maximum in histogram
202 Float_t fHistoPtSumMin; // PtSum minimum in histogram
203 Int_t fHistoNPtInConeBins; // Number of bins in PtInCone histogram
204 Float_t fHistoPtInConeMax; // PtInCone maximum in histogram
205 Float_t fHistoPtInConeMin; // PtInCone maximum in histogram
31ae6d59 206
09273901 207 AliAnaParticleIsolation( const AliAnaParticleIsolation & iso) ; // cpy ctor
208 AliAnaParticleIsolation & operator = (const AliAnaParticleIsolation & iso) ; // cpy assignment
c5693f62 209
31ae6d59 210 ClassDef(AliAnaParticleIsolation,8)
1a31a9ab 211} ;
212
213
214#endif //ALIANAPARTICLEISOLATION_H
215
216
217