]>
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 | ||
22ad7981 | 54 | Bool_t ClusterSelected(AliVCluster* cl, TLorentzVector mom, Int_t nlm) ; |
1c5acb87 | 55 | |
3d5d5078 | 56 | void FillAcceptanceHistograms(); |
57 | ||
22ad7981 | 58 | void FillShowerShapeHistograms( AliVCluster* cluster, Int_t mcTag) ; |
3d5d5078 | 59 | |
c2a62a94 | 60 | void SwitchOnFillShowerShapeHistograms() { fFillSSHistograms = kTRUE ; } |
61 | void SwitchOffFillShowerShapeHistograms() { fFillSSHistograms = kFALSE ; } | |
3d5d5078 | 62 | |
764ab1f4 | 63 | void SwitchOnOnlySimpleSSHistoFill() { fFillOnlySimpleSSHisto = kTRUE ; } |
64 | void SwitchOffOnlySimpleHistoFill() { fFillOnlySimpleSSHisto = kFALSE ; } | |
65 | ||
22ad7981 | 66 | void FillTrackMatchingResidualHistograms(AliVCluster* calo, Int_t cut); |
4bfeae64 | 67 | |
c2a62a94 | 68 | void SwitchOnTMHistoFill() { fFillTMHisto = kTRUE ; } |
69 | void SwitchOffTMHistoFill() { fFillTMHisto = kFALSE ; } | |
4bfeae64 | 70 | |
22ad7981 | 71 | void FillPileUpHistograms(Float_t energy, Float_t pt, Float_t time) ; |
acd56ca4 | 72 | void FillPileUpHistogramsPerEvent(TObjArray * clusters) ; |
73 | ||
c2a62a94 | 74 | void SwitchOnFillPileUpHistograms() { fFillPileUpHistograms = kTRUE ; } |
75 | void SwitchOffFillPileUpHistograms() { fFillPileUpHistograms = kFALSE ; } | |
76 | ||
77 | void SwitchOnFillEMCALBCHistograms() { fFillEMCALBCHistograms = kTRUE ; } | |
78 | void SwitchOffFillEMCALBCHistograms() { fFillEMCALBCHistograms = kFALSE ; } | |
79 | ||
3d5d5078 | 80 | |
6175da48 | 81 | // Analysis parameters setters getters |
c4a7d28a | 82 | |
521636d2 | 83 | TString GetCalorimeter() const { return fCalorimeter ; } |
84 | void SetCalorimeter(TString & det) { fCalorimeter = det ; } | |
85 | ||
6175da48 | 86 | // ** Cluster selection methods ** |
87 | ||
c4a7d28a | 88 | void SetMinDistanceToBadChannel(Float_t m1, Float_t m2, Float_t m3) { |
521636d2 | 89 | fMinDist = m1; fMinDist2 = m2; fMinDist3 = m3; } |
6175da48 | 90 | |
c4a7d28a | 91 | void SetTimeCut(Double_t min, Double_t max) { fTimeCutMin = min; |
521636d2 | 92 | fTimeCutMax = max ; } |
93 | Double_t GetTimeCutMin() const { return fTimeCutMin ; } | |
94 | Double_t GetTimeCutMax() const { return fTimeCutMax ; } | |
1e86c71e | 95 | |
521636d2 | 96 | void SetNCellCut(Int_t n) { fNCellsCut = n ; } |
97 | Double_t GetNCellCut() const { return fNCellsCut ; } | |
c4a7d28a | 98 | |
c2a62a94 | 99 | void SetNLMCut(Int_t min, Int_t max) { fNLMCutMin = min; |
9e51e29a | 100 | fNLMCutMax = max ; } |
c2a62a94 | 101 | Int_t GetNLMCutMin() const { return fNLMCutMin ; } |
102 | Int_t GetNLMCutMax() const { return fNLMCutMax ; } | |
9e51e29a | 103 | |
104 | ||
c4a7d28a | 105 | Bool_t IsTrackMatchRejectionOn() const { return fRejectTrackMatch ; } |
106 | void SwitchOnTrackMatchRejection() { fRejectTrackMatch = kTRUE ; } | |
107 | void SwitchOffTrackMatchRejection() { fRejectTrackMatch = kFALSE ; } | |
09273901 | 108 | |
f66d95af | 109 | void FillNOriginHistograms(Int_t n) { fNOriginHistograms = n ; |
110 | if(n > 14) fNOriginHistograms = 14; } | |
111 | void FillNPrimaryHistograms(Int_t n) { fNPrimaryHistograms= n ; | |
112 | if(n > 7) fNPrimaryHistograms = 7; } | |
113 | ||
3d5d5078 | 114 | // For histograms in arrays, index in the array, corresponding to a particle |
c5693f62 | 115 | enum mcTypes { kmcPhoton = 0, kmcPi0Decay = 1, kmcOtherDecay = 2, |
116 | kmcPi0 = 3, kmcEta = 4, kmcElectron = 5, | |
117 | kmcConversion = 6, kmcOther = 7, kmcAntiNeutron = 8, | |
118 | kmcAntiProton = 9, kmcPrompt = 10, kmcFragmentation = 11, | |
119 | kmcISR = 12, kmcString = 13 }; | |
41121cfe | 120 | |
c5693f62 | 121 | enum mcPTypes { kmcPPhoton = 0, kmcPPi0Decay = 1, kmcPOtherDecay = 2, kmcPOther = 3, |
122 | kmcPPrompt = 4, kmcPFragmentation = 5, kmcPISR = 6 }; | |
f66d95af | 123 | |
c5693f62 | 124 | enum mcssTypes { kmcssPhoton = 0, kmcssOther = 1, kmcssPi0 = 2, |
125 | kmcssEta = 3, kmcssConversion = 4, kmcssElectron = 5 }; | |
3d5d5078 | 126 | |
1c5acb87 | 127 | private: |
128 | ||
6175da48 | 129 | TString fCalorimeter ; // Calorimeter where the gamma is searched; |
130 | Float_t fMinDist ; // Minimal distance to bad channel to accept cluster | |
131 | Float_t fMinDist2; // Cuts on Minimal distance to study acceptance evaluation | |
132 | Float_t fMinDist3; // One more cut on distance used for acceptance-efficiency study | |
133 | Bool_t fRejectTrackMatch ; // If PID on, reject clusters which have an associated TPC track | |
09273901 | 134 | Bool_t fFillTMHisto; // Fill track matching plots |
6175da48 | 135 | Double_t fTimeCutMin ; // Remove clusters/cells with time smaller than this value, in ns |
136 | Double_t fTimeCutMax ; // Remove clusters/cells with time larger than this value, in ns | |
137 | Int_t fNCellsCut ; // Accept for the analysis clusters with more than fNCellsCut cells | |
9e51e29a | 138 | Int_t fNLMCutMin ; // Remove clusters/cells with number of local maxima smaller than this value |
139 | Int_t fNLMCutMax ; // Remove clusters/cells with number of local maxima larger than this value | |
c4a7d28a | 140 | Bool_t fFillSSHistograms ; // Fill shower shape histograms |
764ab1f4 | 141 | Bool_t fFillOnlySimpleSSHisto; // Fill selected cluster histograms, selected SS histograms |
f66d95af | 142 | Int_t fNOriginHistograms; // Fill only NOriginHistograms of the 14 defined types |
143 | Int_t fNPrimaryHistograms; // Fill only NPrimaryHistograms of the 7 defined types | |
2ad19c3d | 144 | Bool_t fFillPileUpHistograms; // Fill pile-up related histograms |
c2a62a94 | 145 | Bool_t fFillEMCALBCHistograms; // Fill eta-phi BC dependent histograms |
2ad19c3d | 146 | |
2244659d | 147 | //Histograms |
9e51e29a | 148 | TH1F * fhClusterCuts[10]; //! control histogram on the different photon selection cuts |
c4a7d28a | 149 | TH2F * fhNCellsE; //! number of cells in cluster vs E |
5c46c992 | 150 | TH2F * fhCellsE; //! energy of cells in cluster vs E of cluster |
f66d95af | 151 | TH2F * fhMaxCellDiffClusterE; //! Fraction of energy carried by cell with maximum energy |
c2a62a94 | 152 | TH2F * fhTimeE; //! time of cluster vs E |
153 | ||
154 | TH2F * fhEtaPhi ; //! Pseudorapidity vs Phi of clusters for E > 0.5 | |
155 | TH2F * fhEtaPhiEMCALBC0 ; //! Pseudorapidity vs Phi of clusters for E > 0.5 | |
156 | TH2F * fhEtaPhiEMCALBC1 ; //! Pseudorapidity vs Phi of clusters for E > 0.5 | |
157 | TH2F * fhEtaPhiEMCALBCN ; //! Pseudorapidity vs Phi of clusters for E > 0.5 | |
158 | ||
afb3af8a | 159 | TH2F * fhEtaPhiTriggerEMCALBC[11] ; //! Pseudorapidity vs Phi of clusters for E > 2 |
160 | TH2F * fhTimeTriggerEMCALBC [11] ; //! Time distribution of clusters, when trigger is in a given BC | |
161 | TH2F * fhTimeTriggerEMCALBCPileUpSPD[11]; //! Time distribution of clusters, when trigger is in a given BC, tagged as pile-up SPD | |
162 | ||
163 | TH2F * fhEtaPhiTriggerEMCALBCUM[11] ; //! Pseudorapidity vs Phi of clusters for E > 2, not matched to trigger | |
164 | TH2F * fhTimeTriggerEMCALBCUM [11] ; //! Time distribution of clusters, when trigger is in a given BC, not matched to trigger | |
165 | ||
166 | TH2F * fhEtaPhiTriggerEMCALBCCluster [11] ; //! Pseudorapidity vs Phi of trigger clusters | |
167 | TH2F * fhTimeTriggerEMCALBCCluster [11] ; //! Time distribution of clusters, when trigger cluster is in a given BC | |
168 | TH2F * fhEtaPhiTriggerEMCALBCUMCluster[11] ; //! Pseudorapidity vs Phi of highest E cluster in event, not matched to trigger | |
169 | TH2F * fhTimeTriggerEMCALBCUMCluster [11] ; //! Time distribution of highest energy cluster in event, when trigger is in a given BC, not | |
170 | ||
171 | TH2F * fhEtaPhiTriggerEMCALBCClusterOverTh ; //! Pseudorapidity vs Phi of trigger clusters, over nominal threshold | |
172 | TH2F * fhEtaPhiTriggerEMCALBCUMClusterOverTh ; //! Pseudorapidity vs Phi of highest E cluster in event, not matched to trigger, over nominal threshold | |
173 | TH2F * fhEtaPhiTriggerEMCALBCClusterBelowTh1 ; //! Pseudorapidity vs Phi of trigger clusters, 1 GeV below nominal threshold | |
174 | TH2F * fhEtaPhiTriggerEMCALBCUMClusterBelowTh1 ; //! Pseudorapidity vs Phi of highest E cluster in event, not matched to trigger, 2 GeV below nominal threshold | |
175 | TH2F * fhEtaPhiTriggerEMCALBCClusterBelowTh2 ; //! Pseudorapidity vs Phi of trigger clusters, 1 GeV below nominal threshold | |
176 | TH2F * fhEtaPhiTriggerEMCALBCUMClusterBelowTh2 ; //! Pseudorapidity vs Phi of highest E cluster in event, not matched to trigger, 2 GeV below nominal threshold | |
177 | ||
178 | TH2F * fhEtaPhiTriggerEMCALBCExotic ; //! Pseudorapidity vs Phi of trigger exotic clusters | |
179 | TH2F * fhTimeTriggerEMCALBCExotic ; //! Time distribution of clusters, when trigger exotic cluster is in a given BC | |
180 | TH2F * fhEtaPhiTriggerEMCALBCUMExotic ; //! Pseudorapidity vs Phi of highest E exotic cluster in event, not matched to trigger | |
181 | TH2F * fhTimeTriggerEMCALBCUMExotic ; //! Time distribution of highest energy exotic cluster in event, when trigger is in a given BC, not | |
182 | ||
183 | TH2F * fhEtaPhiTriggerEMCALBCBad ; //! Pseudorapidity vs Phi of trigger exotic clusters | |
184 | TH2F * fhTimeTriggerEMCALBCBad ; //! Time distribution of clusters, when trigger exotic cluster is in a given BC | |
185 | TH2F * fhEtaPhiTriggerEMCALBCUMBad ; //! Pseudorapidity vs Phi of highest E exotic cluster in event, not matched to trigger | |
186 | TH2F * fhTimeTriggerEMCALBCUMBad ; //! Time distribution of highest energy exotic cluster in event, when trigger is in a given BC, not | |
187 | ||
188 | TH2F * fhEtaPhiTriggerEMCALBCExoticCluster ; //! Pseudorapidity vs Phi of trigger exotic clusters | |
189 | TH2F * fhTimeTriggerEMCALBCExoticCluster ; //! Time distribution of clusters, when trigger exotic cluster is in a given BC | |
190 | TH2F * fhEtaPhiTriggerEMCALBCUMExoticCluster ; //! Pseudorapidity vs Phi of highest E exotic cluster in event, not matched to trigger | |
191 | TH2F * fhTimeTriggerEMCALBCUMExoticCluster ; //! Time distribution of highest energy exotic cluster in event, when trigger is in a given BC, not | |
192 | ||
193 | TH2F * fhEtaPhiTriggerEMCALBCBadCluster ; //! Pseudorapidity vs Phi of trigger exotic clusters | |
194 | TH2F * fhTimeTriggerEMCALBCBadCluster ; //! Time distribution of clusters, when trigger exotic cluster is in a given BC | |
195 | TH2F * fhEtaPhiTriggerEMCALBCUMBadCluster ; //! Pseudorapidity vs Phi of highest E exotic cluster in event, not matched to trigger | |
196 | TH2F * fhTimeTriggerEMCALBCUMBadCluster ; //! Time distribution of highest energy exotic cluster in event, when trigger is in a given BC, not | |
197 | ||
198 | TH2F * fhEtaPhiNoTrigger ; //! Pseudorapidity vs Phi of highest E exotic cluster in event, no trigger at all | |
199 | TH2F * fhTimeNoTrigger ; //! Time distribution of highest energy exotic cluster in event, no trigger at all | |
200 | ||
20218aea | 201 | TH1F * fhEPhoton ; //! Number of identified photon vs energy |
6175da48 | 202 | TH1F * fhPtPhoton ; //! Number of identified photon vs transerse momentum |
203 | TH2F * fhPhiPhoton ; //! Azimuthal angle of identified photon vs transerse momentum | |
204 | TH2F * fhEtaPhoton ; //! Pseudorapidity of identified photon vs transerse momentum | |
c2a62a94 | 205 | TH2F * fhEtaPhiPhoton ; //! Pseudorapidity vs Phi of identified photon for E > 0.5 |
206 | TH2F * fhEtaPhi05Photon ; //! Pseudorapidity vs Phi of identified photon for E < 0.5 | |
207 | TH2F * fhEtaPhiPhotonEMCALBC0 ; //! Pseudorapidity vs Phi of identified photon for E > 0.5 | |
208 | TH2F * fhEtaPhiPhotonEMCALBC1 ; //! Pseudorapidity vs Phi of identified photon for E > 0.5 | |
209 | TH2F * fhEtaPhiPhotonEMCALBCN ; //! Pseudorapidity vs Phi of identified photon for E > 0.5 | |
afb3af8a | 210 | TH2F * fhEtaPhiPhotonTriggerEMCALBC[11]; //! Pseudorapidity vs Phi of photons for E > 0.5 |
211 | TH2F * fhTimePhotonTriggerEMCALBC [11]; //! Time distribution of photons, when trigger is in a given BC | |
212 | TH2F * fhTimePhotonTriggerEMCALBCPileUpSPD[11] ; //! Time distribution of photons, when trigger is in a given BC, tagged as pile-up SPD | |
213 | TH2F * fhEtaPhiPhotonTriggerEMCALBCUM[11]; //! Pseudorapidity vs Phi of photons for E > 2, not matched to trigger | |
214 | TH2F * fhTimePhotonTriggerEMCALBCUM [11]; //! Time distribution of photons, when trigger is in a given BC, not matched to trigger | |
215 | ||
c8710850 | 216 | TH2F * fhPtCentralityPhoton ; //! centrality vs photon pT |
217 | TH2F * fhPtEventPlanePhoton ; //! event plane vs photon pT | |
fedea415 | 218 | |
521636d2 | 219 | //Shower shape |
9e51e29a | 220 | TH2F * fhNLocMax; //! number of maxima in selected clusters |
221 | ||
521636d2 | 222 | TH2F * fhDispE; //! cluster dispersion vs E |
223 | TH2F * fhLam0E; //! cluster lambda0 vs E | |
224 | TH2F * fhLam1E; //! cluster lambda1 vs E | |
7c65ad18 | 225 | |
521636d2 | 226 | TH2F * fhDispETRD; //! cluster dispersion vs E, SM covered by TRD |
227 | TH2F * fhLam0ETRD; //! cluster lambda0 vs E, SM covered by TRD | |
228 | TH2F * fhLam1ETRD; //! cluster lambda1 vs E, SM covered by TRD | |
7c65ad18 | 229 | |
b5dbb99b | 230 | TH2F * fhDispETM; //! cluster dispersion vs E, cut on Track Matching residual |
231 | TH2F * fhLam0ETM; //! cluster lambda0 vs E, cut on Track Matching residual | |
232 | TH2F * fhLam1ETM; //! cluster lambda1 vs E, cut on Track Matching residual | |
233 | ||
234 | TH2F * fhDispETMTRD; //! cluster dispersion vs E, SM covered by TRD, cut on Track Matching residual | |
235 | TH2F * fhLam0ETMTRD; //! cluster lambda0 vs E, SM covered by TRD, cut on Track Matching residual | |
236 | TH2F * fhLam1ETMTRD; //! cluster lambda1 vs E, SM covered by TRD, cut on Track Matching residual | |
237 | ||
521636d2 | 238 | TH2F * fhNCellsLam0LowE; //! number of cells in cluster vs lambda0 |
239 | TH2F * fhNCellsLam1LowE; //! number of cells in cluster vs lambda1 | |
240 | TH2F * fhNCellsDispLowE; //! number of cells in cluster vs dispersion | |
241 | TH2F * fhNCellsLam0HighE; //! number of cells in cluster vs lambda0, E>2 | |
242 | TH2F * fhNCellsLam1HighE; //! number of cells in cluster vs lambda1, E>2 | |
243 | TH2F * fhNCellsDispHighE; //! number of cells in cluster vs dispersion, E>2 | |
244 | ||
521636d2 | 245 | TH2F * fhEtaLam0LowE; //! cluster eta vs lambda0, E<2 |
246 | TH2F * fhPhiLam0LowE; //! cluster phi vs lambda0, E<2 | |
247 | TH2F * fhEtaLam0HighE; //! cluster eta vs lambda0, E>2 | |
248 | TH2F * fhPhiLam0HighE; //! cluster phi vs lambda0, E>2 | |
249 | TH2F * fhLam0DispLowE; //! cluster lambda0 vs dispersion, E<2 | |
250 | TH2F * fhLam0DispHighE; //! cluster lambda0 vs dispersion, E>2 | |
251 | TH2F * fhLam1Lam0LowE; //! cluster lambda1 vs lambda0, E<2 | |
252 | TH2F * fhLam1Lam0HighE; //! cluster lambda1 vs lambda0, E>2 | |
253 | TH2F * fhDispLam1LowE; //! cluster disp vs lambda1, E<2 | |
254 | TH2F * fhDispLam1HighE; //! cluster disp vs lambda1, E>2 | |
7c65ad18 | 255 | |
34c16486 | 256 | TH2F * fhDispEtaE ; //! shower dispersion in eta direction |
257 | TH2F * fhDispPhiE ; //! shower dispersion in phi direction | |
258 | TH2F * fhSumEtaE ; //! shower dispersion in eta direction | |
259 | TH2F * fhSumPhiE ; //! shower dispersion in phi direction | |
260 | TH2F * fhSumEtaPhiE ; //! shower dispersion in eta and phi direction | |
261 | TH2F * fhDispEtaPhiDiffE ; //! shower dispersion eta - phi | |
262 | TH2F * fhSphericityE ; //! shower sphericity in eta vs phi | |
263 | TH2F * fhDispSumEtaDiffE ; //! difference of 2 eta dispersions | |
264 | TH2F * fhDispSumPhiDiffE ; //! difference of 2 phi dispersions | |
d2655d46 | 265 | TH2F * fhDispEtaDispPhi[7] ; //! shower dispersion in eta direction vs phi direction for 5 E bins [0-2],[2-4],[4-6],[6-10],[> 10] |
266 | TH2F * fhLambda0DispEta[7] ; //! shower shape correlation l0 vs disp eta | |
267 | TH2F * fhLambda0DispPhi[7] ; //! shower shape correlation l0 vs disp phi | |
bfdcf7fb | 268 | |
4c8f7c2e | 269 | //Fill MC dependent histograms, Origin of this cluster is ... |
270 | ||
5812a064 | 271 | TH2F * fhMCDeltaE[14] ; //! MC-Reco E distribution coming from MC particle |
272 | TH2F * fhMCDeltaPt[14] ; //! MC-Reco pT distribution coming from MC particle | |
273 | TH2F * fhMC2E[14] ; //! E distribution, Reco vs MC coming from MC particle | |
274 | TH2F * fhMC2Pt[14] ; //! pT distribution, Reco vs MC coming from MC particle | |
4c8f7c2e | 275 | |
5812a064 | 276 | TH1F * fhMCE[14]; //! Number of identified photon vs cluster energy coming from MC particle |
277 | TH1F * fhMCPt[14]; //! Number of identified photon vs cluster pT coming from MC particle | |
278 | TH2F * fhMCPhi[14]; //! Phi of identified photon coming from MC particle | |
279 | TH2F * fhMCEta[14]; //! eta of identified photon coming from MC particle | |
3d5d5078 | 280 | |
5812a064 | 281 | TH1F * fhEPrimMC[7]; //! Number of generated photon vs energy |
282 | TH1F * fhPtPrimMC[7]; //! Number of generated photon vs pT | |
283 | TH2F * fhPhiPrimMC[7]; //! Phi of generted photon | |
284 | TH2F * fhYPrimMC[7]; //! Rapidity of generated photon | |
3d5d5078 | 285 | |
5812a064 | 286 | TH1F * fhEPrimMCAcc[7]; //! Number of generated photon vs energy, in calorimeter acceptance |
287 | TH1F * fhPtPrimMCAcc[7]; //! Number of generated photon vs pT, in calorimeter acceptance | |
288 | TH2F * fhPhiPrimMCAcc[7]; //! Phi of generted photon, in calorimeter acceptance | |
289 | TH2F * fhYPrimMCAcc[7]; //! Rapidity of generated photon, in calorimeter acceptance | |
f66d95af | 290 | |
521636d2 | 291 | // Shower Shape MC |
292 | ||
5812a064 | 293 | TH2F * fhMCELambda0[6] ; //! E vs Lambda0 from MC particle |
294 | TH2F * fhMCELambda1[6] ; //! E vs Lambda1 from MC particle | |
295 | TH2F * fhMCEDispersion[6] ; //! E vs Dispersion from MC particle | |
f66d95af | 296 | |
5812a064 | 297 | TH2F * fhMCPhotonELambda0NoOverlap ; //! E vs Lambda0 from MC photons, no overlap |
298 | TH2F * fhMCPhotonELambda0TwoOverlap ; //! E vs Lambda0 from MC photons, 2 particles overlap | |
299 | TH2F * fhMCPhotonELambda0NOverlap ; //! E vs Lambda0 from MC photons, N particles overlap | |
f66d95af | 300 | |
301 | TH2F * fhMCLambda0vsClusterMaxCellDiffE0[6]; //! Lambda0 vs fraction of energy of max cell for E < 2 GeV | |
302 | TH2F * fhMCLambda0vsClusterMaxCellDiffE2[6]; //! Lambda0 vs fraction of energy of max cell for 2< E < 6 GeV | |
303 | TH2F * fhMCLambda0vsClusterMaxCellDiffE6[6]; //! Lambda0 vs fraction of energy of max cell for E > 6 GeV | |
304 | TH2F * fhMCNCellsvsClusterMaxCellDiffE0[6]; //! NCells vs fraction of energy of max cell for E < 2 | |
305 | TH2F * fhMCNCellsvsClusterMaxCellDiffE2[6]; //! NCells vs fraction of energy of max cell for 2 < E < 6 GeV | |
306 | TH2F * fhMCNCellsvsClusterMaxCellDiffE6[6]; //! NCells vs fraction of energy of max cell for E > 6 | |
307 | TH2F * fhMCNCellsE[6]; //! NCells per cluster vs energy | |
308 | TH2F * fhMCMaxCellDiffClusterE[6]; //! Fraction of energy carried by cell with maximum energy | |
309 | ||
34c16486 | 310 | TH2F * fhMCEDispEta[6] ; //! shower dispersion in eta direction |
311 | TH2F * fhMCEDispPhi[6] ; //! shower dispersion in phi direction | |
312 | TH2F * fhMCESumEtaPhi[6] ; //! shower dispersion in eta vs phi direction | |
313 | TH2F * fhMCEDispEtaPhiDiff[6] ; //! shower dispersion in eta -phi direction | |
314 | TH2F * fhMCESphericity[6] ; //! shower sphericity, eta vs phi | |
d2655d46 | 315 | TH2F * fhMCDispEtaDispPhi[7][6] ; //! shower dispersion in eta direction vs phi direction for 5 E bins [0-2],[2-4],[4-6],[6-10],[> 10] |
316 | TH2F * fhMCLambda0DispEta[7][6] ; //! shower shape correlation l0 vs disp eta | |
317 | TH2F * fhMCLambda0DispPhi[7][6] ; //! shower shape correlation l0 vs disp phi | |
34c16486 | 318 | |
3d5d5078 | 319 | //Embedding |
5812a064 | 320 | TH2F * fhEmbeddedSignalFractionEnergy ; //! Fraction of photon energy of embedded signal vs cluster energy |
3d5d5078 | 321 | |
5812a064 | 322 | TH2F * fhEmbedPhotonELambda0FullSignal ; //! Lambda0 vs E for embedded photons with more than 90% of the cluster energy |
323 | TH2F * fhEmbedPhotonELambda0MostlySignal ; //! Lambda0 vs E for embedded photons with 90%<fraction<50% | |
324 | TH2F * fhEmbedPhotonELambda0MostlyBkg ; //! Lambda0 vs E for embedded photons with 50%<fraction<10% | |
325 | TH2F * fhEmbedPhotonELambda0FullBkg ; //! Lambda0 vs E for embedded photons with less than 10% of the cluster energy | |
3d5d5078 | 326 | |
5812a064 | 327 | TH2F * fhEmbedPi0ELambda0FullSignal ; //! Lambda0 vs E for embedded photons with more than 90% of the cluster energy |
328 | TH2F * fhEmbedPi0ELambda0MostlySignal ; //! Lambda0 vs E for embedded photons with 90%<fraction<50% | |
329 | TH2F * fhEmbedPi0ELambda0MostlyBkg ; //! Lambda0 vs E for embedded photons with 50%<fraction<10% | |
330 | TH2F * fhEmbedPi0ELambda0FullBkg ; //! Lambda0 vs E for embedded photons with less than 10% of the cluster energy | |
3d5d5078 | 331 | |
09273901 | 332 | // Track Matching |
4bfeae64 | 333 | TH2F * fhTrackMatchedDEta[2] ; //! Eta distance between track and cluster vs cluster E, after and before photon cuts |
334 | TH2F * fhTrackMatchedDPhi[2] ; //! Phi distance between track and cluster vs cluster E, after and before photon cuts | |
335 | TH2F * fhTrackMatchedDEtaDPhi[2] ; //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV, after and before photon cuts | |
336 | ||
337 | TH2F * fhTrackMatchedDEtaTRD[2] ; //! Eta distance between track and cluster vs cluster E, after and before photon cuts, behind TRD | |
338 | TH2F * fhTrackMatchedDPhiTRD[2] ; //! Phi distance between track and cluster vs cluster E, after and before photon cuts, behind TRD | |
339 | ||
340 | TH2F * fhTrackMatchedDEtaMCOverlap[2] ; //! Eta distance between track and cluster vs cluster E, several particle overlap, after and before photon cuts | |
341 | TH2F * fhTrackMatchedDPhiMCOverlap[2] ; //! Phi distance between track and cluster vs cluster E, several particle overlap, after and before photon cuts | |
342 | TH2F * fhTrackMatchedDEtaMCNoOverlap[2]; //! Eta distance between track and cluster vs cluster E, not other particle overlap, after and before photon cuts | |
343 | TH2F * fhTrackMatchedDPhiMCNoOverlap[2]; //! Phi distance between track and cluster vs cluster E, not other particle overlap, after and before photon cuts | |
344 | TH2F * fhTrackMatchedDEtaMCConversion[2]; //! Eta distance between track and cluster vs cluster E, originated in conversion, after and before photon cuts | |
345 | TH2F * fhTrackMatchedDPhiMCConversion[2]; //! Phi distance between track and cluster vs cluster E, originated in conversion, after and before photon cuts | |
346 | ||
347 | TH2F * fhTrackMatchedMCParticle[2]; //! Trace origin of matched particle | |
348 | TH2F * fhdEdx[2]; //! matched track dEdx vs cluster E, after and before photon cuts | |
349 | TH2F * fhEOverP[2]; //! matched track E cluster over P track vs cluster E, after dEdx cut, after and before photon cuts | |
350 | TH2F * fhEOverPTRD[2]; //! matched track E cluster over P track vs cluster E, after dEdx cut, after and before photon cuts, behind TRD | |
31ae6d59 | 351 | |
2ad19c3d | 352 | // Pile-up |
fad96885 | 353 | TH1F * fhPtPileUp[7]; //! pT distribution of clusters before any selection |
354 | TH1F * fhPtChargedPileUp[7]; //! pT distribution of track matched clusters | |
5e5e056f | 355 | TH1F * fhPtPhotonPileUp[7]; //! pT distribution of selected photons |
fad96885 | 356 | TH2F * fhLambda0PileUp[7]; //! E vs M02 distribution of clusters, before any selection |
357 | TH2F * fhLambda0ChargedPileUp[7]; //! E vs M02 distribution of clusters, track matched clusters | |
358 | TH2F * fhClusterTimeDiffPileUp[7]; //! E vs Time difference inside cluster, before any selection | |
359 | TH2F * fhClusterTimeDiffChargedPileUp[7]; //! E vs Time difference inside cluster for track matched clusters | |
360 | TH2F * fhClusterTimeDiffPhotonPileUp[7]; //! E vs Time difference inside cluster for selected photons | |
650d1938 | 361 | TH2F * fhClusterEFracLongTimePileUp[7]; //! E vs fraction of cluster energy from cells with large time |
5e5e056f | 362 | TH2F * fhTimeENoCut; //! time of cluster vs E, no cut |
2ad19c3d | 363 | TH2F * fhTimeESPD; //! time of cluster vs E, IsSPDPileUp |
364 | TH2F * fhTimeESPDMulti; //! time of cluster vs E, IsSPDPileUpMulti | |
365 | TH2F * fhTimeNPileUpVertSPD; //! time of cluster vs n pile-up vertices from SPD | |
366 | TH2F * fhTimeNPileUpVertTrack; //! time of cluster vs n pile-up vertices from Tracks | |
367 | TH2F * fhTimeNPileUpVertContributors; //! time of cluster vs n pile-up vertex from SPD contributors | |
368 | TH2F * fhTimePileUpMainVertexZDistance; //! time of cluster vs difference of z main vertex and pile-up vertex | |
369 | TH2F * fhTimePileUpMainVertexZDiamond; //! time of cluster vs difference of z diamond and pile-up vertex | |
acd56ca4 | 370 | TH2F * fhClusterMultSPDPileUp[4]; //! E max cluster vs event cluster multiplicity, for tmax-tdiff cuts, pile up event |
371 | TH2F * fhClusterMultNoPileUp[4]; //! E max cluster vs event cluster multiplicity, for tmax-tdiff cuts, not pile up event | |
fedea415 | 372 | TH2F * fhEtaPhiBC0; //! eta/phi of clusters in BC=0 |
373 | TH2F * fhEtaPhiBCPlus; //! eta/phi of clusters in BC>0 | |
374 | TH2F * fhEtaPhiBCMinus; //! eta/phi of clusters in BC<0 | |
375 | TH2F * fhEtaPhiBC0PileUpSPD; //! eta/phi of clusters in BC=0, SPD pile-up | |
376 | TH2F * fhEtaPhiBCPlusPileUpSPD; //! eta/phi of clusters in BC>0, SPD pile-up | |
377 | TH2F * fhEtaPhiBCMinusPileUpSPD; //! eta/phi of clusters in BC<0, SPD pile-up | |
378 | ||
09273901 | 379 | AliAnaPhoton( const AliAnaPhoton & g) ; // cpy ctor |
c5693f62 | 380 | AliAnaPhoton & operator = (const AliAnaPhoton & g) ; // cpy assignment |
381 | ||
afb3af8a | 382 | ClassDef(AliAnaPhoton,31) |
6639984f | 383 | |
1c5acb87 | 384 | } ; |
385 | ||
1c5acb87 | 386 | #endif//ALIANAPHOTON_H |
387 | ||
388 | ||
389 |