]>
Commit | Line | Data |
---|---|---|
1c5acb87 | 1 | #ifndef ALIANAPHOTON_H |
2 | #define ALIANAPHOTON_H | |
3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
4 | * See cxx source for full Copyright notice */ | |
1c5acb87 | 5 | |
6 | //_________________________________________________________________________ | |
7 | // | |
8 | // Class for the photon identification. | |
9 | // Clusters from calorimeters are identified as photons | |
10 | // and kept in the AOD. Few histograms produced. | |
6175da48 | 11 | // Produces input for other analysis classes like AliAnaPi0, |
12 | // AliAnaParticleHadronCorrelation ... | |
1c5acb87 | 13 | // |
14 | ||
15 | //-- Author: Gustavo Conesa (INFN-LNF) | |
16 | ||
17 | // --- ROOT system --- | |
18 | class TH2F ; | |
123fc3bd | 19 | class TH1F; |
1c5acb87 | 20 | class TString ; |
0c1383b5 | 21 | class TObjString; |
5812a064 | 22 | class TList ; |
1c5acb87 | 23 | |
24 | // --- ANALYSIS system --- | |
745913ae | 25 | #include "AliAnaCaloTrackCorrBaseClass.h" |
1c5acb87 | 26 | |
745913ae | 27 | class AliAnaPhoton : public AliAnaCaloTrackCorrBaseClass { |
1c5acb87 | 28 | |
78219bac | 29 | public: |
5812a064 | 30 | AliAnaPhoton() ; // default ctor |
31 | virtual ~AliAnaPhoton() { ; } // virtual dtor | |
0c1383b5 | 32 | |
6175da48 | 33 | //--------------------------------------- |
34 | // General analysis frame methods | |
35 | //--------------------------------------- | |
c4a7d28a | 36 | |
0c1383b5 | 37 | TObjString * GetAnalysisCuts(); |
6175da48 | 38 | |
0c1383b5 | 39 | TList * GetCreateOutputObjects(); |
c4a7d28a | 40 | |
6175da48 | 41 | void Init(); |
6639984f | 42 | |
6175da48 | 43 | void InitParameters(); |
44 | ||
45 | void MakeAnalysisFillAOD() ; | |
46 | ||
47 | void MakeAnalysisFillHistograms() ; | |
1c5acb87 | 48 | |
6175da48 | 49 | void Print(const Option_t * opt)const; |
521636d2 | 50 | |
3d5d5078 | 51 | |
52 | // Analysis methods | |
53 | ||
521636d2 | 54 | Bool_t ClusterSelected(AliVCluster* cl, TLorentzVector mom) ; |
1c5acb87 | 55 | |
3d5d5078 | 56 | void FillAcceptanceHistograms(); |
57 | ||
3d5d5078 | 58 | void FillShowerShapeHistograms( AliVCluster* cluster, const Int_t mcTag) ; |
59 | ||
60 | void SwitchOnFillShowerShapeHistograms() { fFillSSHistograms = kTRUE ; } | |
61 | void SwitchOffFillShowerShapeHistograms() { fFillSSHistograms = kFALSE ; } | |
62 | ||
4bfeae64 | 63 | void FillTrackMatchingResidualHistograms(AliVCluster* calo, const Int_t cut); |
64 | ||
65 | void SwitchOnTMHistoFill() { fFillTMHisto = kTRUE ; } | |
66 | void SwitchOffTMHistoFill() { fFillTMHisto = kFALSE ; } | |
67 | ||
3d5d5078 | 68 | |
6175da48 | 69 | // Analysis parameters setters getters |
c4a7d28a | 70 | |
521636d2 | 71 | TString GetCalorimeter() const { return fCalorimeter ; } |
72 | void SetCalorimeter(TString & det) { fCalorimeter = det ; } | |
73 | ||
6175da48 | 74 | // ** Cluster selection methods ** |
75 | ||
c4a7d28a | 76 | void SetMinDistanceToBadChannel(Float_t m1, Float_t m2, Float_t m3) { |
521636d2 | 77 | fMinDist = m1; fMinDist2 = m2; fMinDist3 = m3; } |
6175da48 | 78 | |
c4a7d28a | 79 | void SetTimeCut(Double_t min, Double_t max) { fTimeCutMin = min; |
521636d2 | 80 | fTimeCutMax = max ; } |
81 | Double_t GetTimeCutMin() const { return fTimeCutMin ; } | |
82 | Double_t GetTimeCutMax() const { return fTimeCutMax ; } | |
1e86c71e | 83 | |
521636d2 | 84 | void SetNCellCut(Int_t n) { fNCellsCut = n ; } |
85 | Double_t GetNCellCut() const { return fNCellsCut ; } | |
c4a7d28a | 86 | |
87 | Bool_t IsTrackMatchRejectionOn() const { return fRejectTrackMatch ; } | |
88 | void SwitchOnTrackMatchRejection() { fRejectTrackMatch = kTRUE ; } | |
89 | void SwitchOffTrackMatchRejection() { fRejectTrackMatch = kFALSE ; } | |
09273901 | 90 | |
f66d95af | 91 | void FillNOriginHistograms(Int_t n) { fNOriginHistograms = n ; |
92 | if(n > 14) fNOriginHistograms = 14; } | |
93 | void FillNPrimaryHistograms(Int_t n) { fNPrimaryHistograms= n ; | |
94 | if(n > 7) fNPrimaryHistograms = 7; } | |
95 | ||
3d5d5078 | 96 | // For histograms in arrays, index in the array, corresponding to a particle |
c5693f62 | 97 | enum mcTypes { kmcPhoton = 0, kmcPi0Decay = 1, kmcOtherDecay = 2, |
98 | kmcPi0 = 3, kmcEta = 4, kmcElectron = 5, | |
99 | kmcConversion = 6, kmcOther = 7, kmcAntiNeutron = 8, | |
100 | kmcAntiProton = 9, kmcPrompt = 10, kmcFragmentation = 11, | |
101 | kmcISR = 12, kmcString = 13 }; | |
41121cfe | 102 | |
c5693f62 | 103 | enum mcPTypes { kmcPPhoton = 0, kmcPPi0Decay = 1, kmcPOtherDecay = 2, kmcPOther = 3, |
104 | kmcPPrompt = 4, kmcPFragmentation = 5, kmcPISR = 6 }; | |
f66d95af | 105 | |
c5693f62 | 106 | enum mcssTypes { kmcssPhoton = 0, kmcssOther = 1, kmcssPi0 = 2, |
107 | kmcssEta = 3, kmcssConversion = 4, kmcssElectron = 5 }; | |
3d5d5078 | 108 | |
1c5acb87 | 109 | private: |
110 | ||
6175da48 | 111 | TString fCalorimeter ; // Calorimeter where the gamma is searched; |
112 | Float_t fMinDist ; // Minimal distance to bad channel to accept cluster | |
113 | Float_t fMinDist2; // Cuts on Minimal distance to study acceptance evaluation | |
114 | Float_t fMinDist3; // One more cut on distance used for acceptance-efficiency study | |
115 | Bool_t fRejectTrackMatch ; // If PID on, reject clusters which have an associated TPC track | |
09273901 | 116 | Bool_t fFillTMHisto; // Fill track matching plots |
6175da48 | 117 | Double_t fTimeCutMin ; // Remove clusters/cells with time smaller than this value, in ns |
118 | Double_t fTimeCutMax ; // Remove clusters/cells with time larger than this value, in ns | |
119 | Int_t fNCellsCut ; // Accept for the analysis clusters with more than fNCellsCut cells | |
c4a7d28a | 120 | Bool_t fFillSSHistograms ; // Fill shower shape histograms |
f66d95af | 121 | Int_t fNOriginHistograms; // Fill only NOriginHistograms of the 14 defined types |
122 | Int_t fNPrimaryHistograms; // Fill only NPrimaryHistograms of the 7 defined types | |
521636d2 | 123 | |
2244659d | 124 | //Histograms |
fc195fd0 | 125 | TH1F * fhClusterCuts[9]; //! control histogram on the different photon selection cuts |
c4a7d28a | 126 | TH2F * fhNCellsE; //! number of cells in cluster vs E |
f66d95af | 127 | TH2F * fhMaxCellDiffClusterE; //! Fraction of energy carried by cell with maximum energy |
f15c25da | 128 | TH2F * fhTimeE; //! time of cluster vs E |
129 | ||
20218aea | 130 | TH1F * fhEPhoton ; //! Number of identified photon vs energy |
6175da48 | 131 | TH1F * fhPtPhoton ; //! Number of identified photon vs transerse momentum |
132 | TH2F * fhPhiPhoton ; //! Azimuthal angle of identified photon vs transerse momentum | |
133 | TH2F * fhEtaPhoton ; //! Pseudorapidity of identified photon vs transerse momentum | |
134 | TH2F * fhEtaPhiPhoton ; //! Pseudorapidity vs Phi of identified photon for transerse momentum > 0.5 | |
135 | TH2F * fhEtaPhi05Photon ; //! Pseudorapidity vs Phi of identified photon for transerse momentum < 0.5 | |
123fc3bd | 136 | |
521636d2 | 137 | //Shower shape |
f66d95af | 138 | |
521636d2 | 139 | TH2F * fhDispE; //! cluster dispersion vs E |
140 | TH2F * fhLam0E; //! cluster lambda0 vs E | |
141 | TH2F * fhLam1E; //! cluster lambda1 vs E | |
7c65ad18 | 142 | |
521636d2 | 143 | TH2F * fhDispETRD; //! cluster dispersion vs E, SM covered by TRD |
144 | TH2F * fhLam0ETRD; //! cluster lambda0 vs E, SM covered by TRD | |
145 | TH2F * fhLam1ETRD; //! cluster lambda1 vs E, SM covered by TRD | |
7c65ad18 | 146 | |
b5dbb99b | 147 | TH2F * fhDispETM; //! cluster dispersion vs E, cut on Track Matching residual |
148 | TH2F * fhLam0ETM; //! cluster lambda0 vs E, cut on Track Matching residual | |
149 | TH2F * fhLam1ETM; //! cluster lambda1 vs E, cut on Track Matching residual | |
150 | ||
151 | TH2F * fhDispETMTRD; //! cluster dispersion vs E, SM covered by TRD, cut on Track Matching residual | |
152 | TH2F * fhLam0ETMTRD; //! cluster lambda0 vs E, SM covered by TRD, cut on Track Matching residual | |
153 | TH2F * fhLam1ETMTRD; //! cluster lambda1 vs E, SM covered by TRD, cut on Track Matching residual | |
154 | ||
155 | ||
521636d2 | 156 | TH2F * fhNCellsLam0LowE; //! number of cells in cluster vs lambda0 |
157 | TH2F * fhNCellsLam1LowE; //! number of cells in cluster vs lambda1 | |
158 | TH2F * fhNCellsDispLowE; //! number of cells in cluster vs dispersion | |
159 | TH2F * fhNCellsLam0HighE; //! number of cells in cluster vs lambda0, E>2 | |
160 | TH2F * fhNCellsLam1HighE; //! number of cells in cluster vs lambda1, E>2 | |
161 | TH2F * fhNCellsDispHighE; //! number of cells in cluster vs dispersion, E>2 | |
162 | ||
521636d2 | 163 | TH2F * fhEtaLam0LowE; //! cluster eta vs lambda0, E<2 |
164 | TH2F * fhPhiLam0LowE; //! cluster phi vs lambda0, E<2 | |
165 | TH2F * fhEtaLam0HighE; //! cluster eta vs lambda0, E>2 | |
166 | TH2F * fhPhiLam0HighE; //! cluster phi vs lambda0, E>2 | |
167 | TH2F * fhLam0DispLowE; //! cluster lambda0 vs dispersion, E<2 | |
168 | TH2F * fhLam0DispHighE; //! cluster lambda0 vs dispersion, E>2 | |
169 | TH2F * fhLam1Lam0LowE; //! cluster lambda1 vs lambda0, E<2 | |
170 | TH2F * fhLam1Lam0HighE; //! cluster lambda1 vs lambda0, E>2 | |
171 | TH2F * fhDispLam1LowE; //! cluster disp vs lambda1, E<2 | |
172 | TH2F * fhDispLam1HighE; //! cluster disp vs lambda1, E>2 | |
7c65ad18 | 173 | |
4c8f7c2e | 174 | //Fill MC dependent histograms, Origin of this cluster is ... |
175 | ||
5812a064 | 176 | TH2F * fhMCDeltaE[14] ; //! MC-Reco E distribution coming from MC particle |
177 | TH2F * fhMCDeltaPt[14] ; //! MC-Reco pT distribution coming from MC particle | |
178 | TH2F * fhMC2E[14] ; //! E distribution, Reco vs MC coming from MC particle | |
179 | TH2F * fhMC2Pt[14] ; //! pT distribution, Reco vs MC coming from MC particle | |
4c8f7c2e | 180 | |
5812a064 | 181 | TH1F * fhMCE[14]; //! Number of identified photon vs cluster energy coming from MC particle |
182 | TH1F * fhMCPt[14]; //! Number of identified photon vs cluster pT coming from MC particle | |
183 | TH2F * fhMCPhi[14]; //! Phi of identified photon coming from MC particle | |
184 | TH2F * fhMCEta[14]; //! eta of identified photon coming from MC particle | |
3d5d5078 | 185 | |
5812a064 | 186 | TH1F * fhEPrimMC[7]; //! Number of generated photon vs energy |
187 | TH1F * fhPtPrimMC[7]; //! Number of generated photon vs pT | |
188 | TH2F * fhPhiPrimMC[7]; //! Phi of generted photon | |
189 | TH2F * fhYPrimMC[7]; //! Rapidity of generated photon | |
3d5d5078 | 190 | |
5812a064 | 191 | TH1F * fhEPrimMCAcc[7]; //! Number of generated photon vs energy, in calorimeter acceptance |
192 | TH1F * fhPtPrimMCAcc[7]; //! Number of generated photon vs pT, in calorimeter acceptance | |
193 | TH2F * fhPhiPrimMCAcc[7]; //! Phi of generted photon, in calorimeter acceptance | |
194 | TH2F * fhYPrimMCAcc[7]; //! Rapidity of generated photon, in calorimeter acceptance | |
f66d95af | 195 | |
521636d2 | 196 | // Shower Shape MC |
197 | ||
5812a064 | 198 | TH2F * fhMCELambda0[6] ; //! E vs Lambda0 from MC particle |
199 | TH2F * fhMCELambda1[6] ; //! E vs Lambda1 from MC particle | |
200 | TH2F * fhMCEDispersion[6] ; //! E vs Dispersion from MC particle | |
f66d95af | 201 | |
5812a064 | 202 | TH2F * fhMCPhotonELambda0NoOverlap ; //! E vs Lambda0 from MC photons, no overlap |
203 | TH2F * fhMCPhotonELambda0TwoOverlap ; //! E vs Lambda0 from MC photons, 2 particles overlap | |
204 | TH2F * fhMCPhotonELambda0NOverlap ; //! E vs Lambda0 from MC photons, N particles overlap | |
f66d95af | 205 | |
206 | TH2F * fhMCLambda0vsClusterMaxCellDiffE0[6]; //! Lambda0 vs fraction of energy of max cell for E < 2 GeV | |
207 | TH2F * fhMCLambda0vsClusterMaxCellDiffE2[6]; //! Lambda0 vs fraction of energy of max cell for 2< E < 6 GeV | |
208 | TH2F * fhMCLambda0vsClusterMaxCellDiffE6[6]; //! Lambda0 vs fraction of energy of max cell for E > 6 GeV | |
209 | TH2F * fhMCNCellsvsClusterMaxCellDiffE0[6]; //! NCells vs fraction of energy of max cell for E < 2 | |
210 | TH2F * fhMCNCellsvsClusterMaxCellDiffE2[6]; //! NCells vs fraction of energy of max cell for 2 < E < 6 GeV | |
211 | TH2F * fhMCNCellsvsClusterMaxCellDiffE6[6]; //! NCells vs fraction of energy of max cell for E > 6 | |
212 | TH2F * fhMCNCellsE[6]; //! NCells per cluster vs energy | |
213 | TH2F * fhMCMaxCellDiffClusterE[6]; //! Fraction of energy carried by cell with maximum energy | |
214 | ||
3d5d5078 | 215 | //Embedding |
5812a064 | 216 | TH2F * fhEmbeddedSignalFractionEnergy ; //! Fraction of photon energy of embedded signal vs cluster energy |
3d5d5078 | 217 | |
5812a064 | 218 | TH2F * fhEmbedPhotonELambda0FullSignal ; //! Lambda0 vs E for embedded photons with more than 90% of the cluster energy |
219 | TH2F * fhEmbedPhotonELambda0MostlySignal ; //! Lambda0 vs E for embedded photons with 90%<fraction<50% | |
220 | TH2F * fhEmbedPhotonELambda0MostlyBkg ; //! Lambda0 vs E for embedded photons with 50%<fraction<10% | |
221 | TH2F * fhEmbedPhotonELambda0FullBkg ; //! Lambda0 vs E for embedded photons with less than 10% of the cluster energy | |
3d5d5078 | 222 | |
5812a064 | 223 | TH2F * fhEmbedPi0ELambda0FullSignal ; //! Lambda0 vs E for embedded photons with more than 90% of the cluster energy |
224 | TH2F * fhEmbedPi0ELambda0MostlySignal ; //! Lambda0 vs E for embedded photons with 90%<fraction<50% | |
225 | TH2F * fhEmbedPi0ELambda0MostlyBkg ; //! Lambda0 vs E for embedded photons with 50%<fraction<10% | |
226 | TH2F * fhEmbedPi0ELambda0FullBkg ; //! Lambda0 vs E for embedded photons with less than 10% of the cluster energy | |
3d5d5078 | 227 | |
09273901 | 228 | // Track Matching |
4bfeae64 | 229 | TH2F * fhTrackMatchedDEta[2] ; //! Eta distance between track and cluster vs cluster E, after and before photon cuts |
230 | TH2F * fhTrackMatchedDPhi[2] ; //! Phi distance between track and cluster vs cluster E, after and before photon cuts | |
231 | TH2F * fhTrackMatchedDEtaDPhi[2] ; //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV, after and before photon cuts | |
232 | ||
233 | TH2F * fhTrackMatchedDEtaTRD[2] ; //! Eta distance between track and cluster vs cluster E, after and before photon cuts, behind TRD | |
234 | TH2F * fhTrackMatchedDPhiTRD[2] ; //! Phi distance between track and cluster vs cluster E, after and before photon cuts, behind TRD | |
235 | ||
236 | TH2F * fhTrackMatchedDEtaMCOverlap[2] ; //! Eta distance between track and cluster vs cluster E, several particle overlap, after and before photon cuts | |
237 | TH2F * fhTrackMatchedDPhiMCOverlap[2] ; //! Phi distance between track and cluster vs cluster E, several particle overlap, after and before photon cuts | |
238 | TH2F * fhTrackMatchedDEtaMCNoOverlap[2]; //! Eta distance between track and cluster vs cluster E, not other particle overlap, after and before photon cuts | |
239 | TH2F * fhTrackMatchedDPhiMCNoOverlap[2]; //! Phi distance between track and cluster vs cluster E, not other particle overlap, after and before photon cuts | |
240 | TH2F * fhTrackMatchedDEtaMCConversion[2]; //! Eta distance between track and cluster vs cluster E, originated in conversion, after and before photon cuts | |
241 | TH2F * fhTrackMatchedDPhiMCConversion[2]; //! Phi distance between track and cluster vs cluster E, originated in conversion, after and before photon cuts | |
242 | ||
243 | TH2F * fhTrackMatchedMCParticle[2]; //! Trace origin of matched particle | |
244 | TH2F * fhdEdx[2]; //! matched track dEdx vs cluster E, after and before photon cuts | |
245 | TH2F * fhEOverP[2]; //! matched track E cluster over P track vs cluster E, after dEdx cut, after and before photon cuts | |
246 | TH2F * fhEOverPTRD[2]; //! matched track E cluster over P track vs cluster E, after dEdx cut, after and before photon cuts, behind TRD | |
31ae6d59 | 247 | |
09273901 | 248 | AliAnaPhoton( const AliAnaPhoton & g) ; // cpy ctor |
c5693f62 | 249 | AliAnaPhoton & operator = (const AliAnaPhoton & g) ; // cpy assignment |
250 | ||
4bfeae64 | 251 | ClassDef(AliAnaPhoton,23) |
6639984f | 252 | |
1c5acb87 | 253 | } ; |
254 | ||
1c5acb87 | 255 | #endif//ALIANAPHOTON_H |
256 | ||
257 | ||
258 |