]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG4/PartCorrDep/AliAnaPi0.h
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Reader: Add option to remove or not event with primary vertex not reconstructed
[u/mrichter/AliRoot.git] / PWG4 / PartCorrDep / AliAnaPi0.h
CommitLineData
1c5acb87 1#ifndef ALIANAPI0_H
2#define ALIANAPI0_H
3/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
5/* $Id: $ */
6
7//_________________________________________________________________________
6175da48 8// Class to fill two-photon invariant mass histograms
1c5acb87 9// to be used to extract pi0 raw yield.
6175da48 10// Input is produced by AliAnaPhoton (or any other analysis producing output AliAODPWG4Particles),
11// it will do nothing if executed alone
1c5acb87 12//
13//-- Author: Dmitri Peressounko (RRC "KI")
14//-- Adapted to PartCorr frame by Lamia Benhabib (SUBATECH)
15//-- and Gustavo Conesa (INFN-Frascati)
16
17//Root
18class TList;
19class TH3D ;
50f39b97 20class TH2D ;
0c1383b5 21class TObjString;
1c5acb87 22
23//Analysis
5025c139 24#include "AliAnaPartCorrBaseClass.h"
1c5acb87 25class AliAODEvent ;
26class AliESDEvent ;
c8fe2783 27class AliAODPWG4Particle ;
1c5acb87 28
1c5acb87 29class AliAnaPi0 : public AliAnaPartCorrBaseClass {
6639984f 30
78219bac 31 public:
6639984f 32 AliAnaPi0() ; // default ctor
6639984f 33 virtual ~AliAnaPi0() ;//virtual dtor
78219bac 34 private:
35 AliAnaPi0(const AliAnaPi0 & g) ; // cpy ctor
614701c6 36 AliAnaPi0 & operator = (const AliAnaPi0 & api0) ;//cpy assignment
6639984f 37
78219bac 38 public:
39
6175da48 40 //-------------------------------
41 // General analysis frame methods
42 //-------------------------------
43
0c1383b5 44 TObjString * GetAnalysisCuts();
45 TList * GetCreateOutputObjects();
6175da48 46 void Terminate(TList* outputList);
47 void ReadHistograms(TList * outputList); //Fill histograms with histograms in ouput list, needed in Terminate.
48 void Print(const Option_t * opt) const;
6639984f 49 //void MakeAnalysisFillAOD() {;} //Not needed
50 void MakeAnalysisFillHistograms();
6175da48 51 //void Init();
52 void InitParameters();
53
54 //Calorimeter options
156549ae 55 TString GetCalorimeter() const { return fCalorimeter; }
56 void SetCalorimeter(TString & det) { fCalorimeter = det ; }
57 void SetNumberOfModules(Int_t nmod) { fNModules = nmod; }
c8fe2783 58
6175da48 59 //-------------------------------
60 // EVENT Bin Methods
61 //-------------------------------
62
c8fe2783 63 virtual Int_t GetEventIndex(AliAODPWG4Particle * part, Double_t * vert) ;
64
20218aea 65 void CountAndGetAverages(Int_t &nClus,Int_t &nCell, Float_t &eClusTot,Float_t &eCellTot, Float_t &eDenClus,Float_t &eDenCell) ;
66
6175da48 67 //Setters for parameters of event buffers
68 void SetNCentrBin(Int_t n=5) {fNCentrBin=n ;} //number of bins in centrality
69//void SetNZvertBin(Int_t n=5) {fNZvertBin=n ;} //number of bins for vertex position
70//void SetNRPBin(Int_t n=6) {fNrpBin=n ;} //number of bins in reaction plain
71 void SetNMaxEvMix(Int_t n=20) {fNmaxMixEv=n ;} //Maximal number of events for mixing
72
73 //Switchs for event multiplicity bin option, by default, centrality
74 void SwitchOnTrackMultBins() {fUseTrackMultBins = kTRUE ; }
75 void SwitchOffTrackMultBins() {fUseTrackMultBins = kFALSE ; }
76
77 void SwitchOnPhotonMultBins() {fUsePhotonMultBins = kTRUE ; }
78 void SwitchOffPhotonMultBins() {fUsePhotonMultBins = kFALSE ; }
79
80 void SwitchOnClusterEBins() {fUseAverClusterEBins = kTRUE ; }
81 void SwitchOffClusterEBins() {fUseAverClusterEBins = kFALSE ; }
398c93cc 82
6175da48 83 void SwitchOnCellEBins() {fUseAverCellEBins = kTRUE ; }
84 void SwitchOffCellEBins() {fUseAverCellEBins = kFALSE ; }
7e7694bb 85
20218aea 86 void SwitchOnClusterEDenBins() {fUseAverClusterEDenBins = kTRUE ; }
87 void SwitchOffClusterEDenBins() {fUseAverClusterEDenBins = kFALSE ; }
156549ae 88
89// void SwitchOnClusterPairRBins() {fUseAverClusterPairRBins = kTRUE ; }
90// void SwitchOffClusterPairRBins() {fUseAverClusterPairRBins = kFALSE ; }
91//
92// void SwitchOnClusterPairRWeightBins() {fUseAverClusterPairRWeightBins = kTRUE ; }
93// void SwitchOffClusterPairRWeightBins(){fUseAverClusterPairRWeightBins = kFALSE ; }
94
95// void SwitchOnEMaxBins() {fUseEMaxBins = kTRUE ; }
96// void SwitchOffEMaxBins() {fUseEMaxBins = kFALSE ; }
97
6175da48 98 //-------------------------------
99 //Opening angle pair selection
100 //-------------------------------
101 void SwitchOnAngleSelection() {fUseAngleCut = kTRUE ; }
102 void SwitchOffAngleSelection() {fUseAngleCut = kFALSE ; }
103 void SwitchOnAngleEDepSelection() {fUseAngleEDepCut = kTRUE ; }
104 void SwitchOffAngleEDepSelection() {fUseAngleEDepCut = kFALSE ; }
105 void SetAngleCut(Float_t a) {fAngleCut = a ; }
106 void SetAngleMaxCut(Float_t a) {fAngleMaxCut = a ; }
107
108 //-------------------------------
109 // Use mixing code of this class
110 //-------------------------------
111 void SwitchOnOwnMix() {fDoOwnMix = kTRUE ; }
112 void SwitchOffOwnMix() {fDoOwnMix = kFALSE ; }
113
114 //------------------------------------------
20218aea 115 //Do analysis only with clusters in same SM or different combinations of SM
6175da48 116 //------------------------------------------
117 void SwitchOnSameSM() {fSameSM = kTRUE ; }
118 void SwitchOffSameSM() {fSameSM = kFALSE ; }
119
20218aea 120 void SwitchOnSMCombinations() {fFillSMCombinations = kTRUE ; }
121 void SwitchOffSMCombinations() {fFillSMCombinations = kFALSE ; }
122
6175da48 123 //-------------------------------
124 //Histogram filling options off by default
125 //-------------------------------
126 void SwitchOnInvPtWeight() {fMakeInvPtPlots = kTRUE ; }
127 void SwitchOffInvPtWeight() {fMakeInvPtPlots = kFALSE ; }
af7b3903 128
6175da48 129 void SwitchOnFillBadDistHisto() {fFillBadDistHisto = kTRUE;}
130 void SwitchOffFillBadDistHisto() {fFillBadDistHisto = kFALSE;}
131
132 //-------------------------------------------
133 //Cuts for multiple analysis, off by default
134 //-------------------------------------------
156549ae 135 void SwitchOnMultipleCutAnalysis() {fMultiCutAna = kTRUE ;}
6175da48 136 void SwitchOffMultipleCutAnalysis() {fMultiCutAna = kFALSE;}
5ae09196 137
6175da48 138 void SetNPtCuts (Int_t size) {if(size <= 10)fNPtCuts = size; }
139 void SetNAsymCuts (Int_t size) {if(size <= 10)fNAsymCuts = size; }
140 void SetNNCellCuts(Int_t size) {if(size <= 10)fNCellNCuts = size; }
141 void SetNPIDBits (Int_t size) {if(size <= 10)fNPIDBits = size; }
d7c10d78 142
af7b3903 143 void SetPtCutsAt (Int_t pos,Float_t val) {if(pos < 10)fPtCuts[pos] = val;}
144 void SetAsymCutsAt (Int_t pos,Float_t val) {if(pos < 10)fAsymCuts[pos] = val;}
145 void SetNCellCutsAt(Int_t pos,Int_t val) {if(pos < 10)fCellNCuts[pos] = val;}
146 void SetPIDBitsAt (Int_t pos,Int_t val) {if(pos < 10)fPIDBits[pos] = val;}
d7c10d78 147
6175da48 148 //MC analysis related methods
149 void FillAcceptanceHistograms();
150 void FillMCVersusRecDataHistograms(const Int_t index1, const Int_t index2,
151 const Float_t pt1, const Float_t pt2,
152 const Int_t ncells1, const Int_t ncells2,
153 const Double_t mass, const Double_t pt, const Double_t asym,
154 const Double_t deta, const Double_t dphi);
155
156 void SwitchOnMultipleCutAnalysisInSimulation() {fMultiCutAnaSim = kTRUE;}
157 void SwitchOffMultipleCutAnalysisInSimulation() {fMultiCutAnaSim = kFALSE;}
158
20218aea 159 void SwitchOnConversionChecker() { fCheckConversion = kTRUE ; }
160 void SwitchOffConversionChecker() { fCheckConversion = kFALSE ; }
161
6639984f 162
163 private:
6175da48 164 Bool_t fDoOwnMix; // Do combinatorial background not the one provided by the frame
165 Int_t fNCentrBin ; // Number of bins in event container for centrality
166//Int_t fNZvertBin ; // Number of bins in event container for vertex position
167//Int_t fNrpBin ; // Number of bins in event container for reaction plain
168 Int_t fNmaxMixEv ; // Maximal number of events stored in buffer for mixing
169 TString fCalorimeter ; // Select Calorimeter for IM
170 Int_t fNModules ; // Number of EMCAL/PHOS modules, set as many histogras as modules
171 Bool_t fUseAngleCut ; // Select pairs depending on their opening angle
172 Bool_t fUseAngleEDepCut ; // Select pairs depending on their opening angle
173 Float_t fAngleCut ; // Select pairs with opening angle larger than a threshold
174 Float_t fAngleMaxCut ; // Select pairs with opening angle smaller than a threshold
175 TList ** fEventsList ; //![fNCentrBin*GetNZvertBin()*GetNRPBin()] Containers for photons in stored events
176
177 //Multiple cuts analysis
178 Bool_t fMultiCutAna; // Do analysis with several or fixed cut
179 Bool_t fMultiCutAnaSim; // Do analysis with several or fixed cut, in the simulation related part
180 Int_t fNPtCuts; // Number of pt cuts
181 Float_t fPtCuts[10]; // Array with different pt cuts
182 Int_t fNAsymCuts; // Number of assymmetry cuts
183 Float_t fAsymCuts[10]; // Array with different assymetry cuts
184 Int_t fNCellNCuts; // Number of cuts with number of cells in cluster
185 Int_t fCellNCuts[10]; // Array with different cell number cluster cuts
186 Int_t fNPIDBits ; // Number of possible PID bit combinations
187 Int_t fPIDBits[10]; // Array with different PID bits
188
189 //Switchs of different analysis options
190 Bool_t fMakeInvPtPlots; // D plots with inverse pt weight
191 Bool_t fSameSM; // Select only pairs in same SM;
20218aea 192 Bool_t fFillSMCombinations; // Fill histograms with different cluster pairs in SM combinations
193 Bool_t fCheckConversion; // Fill histograms with tagged photons as conversion
6175da48 194 Bool_t fUseTrackMultBins; // Use track multiplicity and not centrality bins
195 Bool_t fUsePhotonMultBins; // Use photon multiplicity and not centrality bins
6175da48 196 Bool_t fUseAverCellEBins; // Use cell average energy and not centrality bins
156549ae 197 Bool_t fUseAverClusterEBins; // Use cluster average energy and not centrality bins
198 Bool_t fUseAverClusterEDenBins; // Use cluster average energy density and not centrality bins
199// Bool_t fUseAverClusterPairRBins; // Use cluster average energy and not centrality bins
200// Bool_t fUseAverClusterPairRWeightBins; // Use cluster average energy and not centrality bins
201// Bool_t fUseEMaxBins; // Use Emax bins
6175da48 202 Bool_t fFillBadDistHisto; // Do plots for different distances to bad channels
af7b3903 203
6639984f 204 //Histograms
6175da48 205
156549ae 206 //Event characterization
6175da48 207 TH1F* fhAverTotECluster; //! Average number of clusters in SM
208 TH1F* fhAverTotECell; //! Average number of cells in SM
156549ae 209 TH2F* fhAverTotECellvsCluster; //! Average number of cells in SM
210 TH1F* fhEDensityCluster; //! Deposited energy in event per cluster
211 TH1F* fhEDensityCell; //! Deposited energy in event per cell vs cluster
212 TH2F* fhEDensityCellvsCluster; //! Deposited energy in event per cell vs cluster
213// TH1F* fhClusterPairDist; //! Distance between clusters
214// TH1F* fhClusterPairDistWeight; //! Distance between clusters weighted by pair energy
215// TH1F* fhAverClusterPairDist; //! Average distance between cluster pairs
216// TH1F* fhAverClusterPairDistWeight;//! Average distance between cluster pairs weighted with pair energy
217// TH2F* fhAverClusterPairDistvsAverE; //! Average distance between cluster pairs vs average cluster energy
218// TH2F* fhAverClusterPairDistWeightvsAverE; //! Average distance between cluster pairs weighted with pair energy vs average cluster energy
219// TH2F* fhAverClusterPairDistvsN; //! Average distance between cluster pairs vs number of clusters
220// TH2F* fhAverClusterPairDistWeightvsN; //! Average distance between cluster pairs weighted with pair energy vs number of clusters
221// TH2F* fhMaxEvsClustMult; //!
222// TH2F* fhMaxEvsClustEDen; //!
223
6175da48 224
225 TH2D ** fhReMod ; //![fNModules] REAL two-photon invariant mass distribution for different calorimeter modules.
8d230fa8 226 TH2D ** fhReSameSideEMCALMod ; //![fNModules-2] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules.
227 TH2D ** fhReSameSectorEMCALMod ; //![fNModules/2] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules.
228 TH2D ** fhReDiffPHOSMod ; //![fNModules] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules.
6175da48 229 TH2D ** fhMiMod ; //![fNModules] MIXED two-photon invariant mass distribution for different calorimeter modules.
8d230fa8 230 TH2D ** fhMiSameSideEMCALMod ; //![fNModules-2] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules.
231 TH2D ** fhMiSameSectorEMCALMod ; //![fNModules/2] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules.
232 TH2D ** fhMiDiffPHOSMod ; //![fNModules-1] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules.
6175da48 233
234 // Pairs with at least one cluster tagged as conversion
235 TH2D * fhReConv ; //! REAL two-photon invariant mass distribution one of the pair was 2 clusters with small mass
236 TH2D * fhMiConv ; //! MIXED two-photon invariant mass distribution one of the pair was 2 clusters with small mass
237 TH2D * fhReConv2 ; //! REAL two-photon invariant mass distribution both pair photons recombined from 2 clusters with small mass
238 TH2D * fhMiConv2 ; //! MIXED two-photon invariant mass distribution both pair photons recombined from 2 clusters with small mass
239
240 TH2D ** fhRe1 ; //![fNCentrBin*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry
241 TH2D ** fhMi1 ; //![fNCentrBin*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry
242 TH2D ** fhRe2 ; //![fNCentrBin*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry
243 TH2D ** fhMi2 ; //![fNCentrBin*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry
244 TH2D ** fhRe3 ; //![fNCentrBin*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry
245 TH2D ** fhMi3 ; //![fNCentrBin*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry
246
247 //Histograms weighted by inverse pT
248 TH2D ** fhReInvPt1 ; //![fNCentrBin*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT
249 TH2D ** fhMiInvPt1 ; //![fNCentrBin*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT
250 TH2D ** fhReInvPt2 ; //![fNCentrBin*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT
251 TH2D ** fhMiInvPt2 ; //![fNCentrBin*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT
252 TH2D ** fhReInvPt3 ; //![fNCentrBin*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT
253 TH2D ** fhMiInvPt3 ; //![fNCentrBin*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT
254
255 //Multiple cuts: Assymmetry, pt, n cells, PID
256 TH2D ** fhRePtNCellAsymCuts ; //![fNPtCuts*fNAsymCuts*fNCellNCuts] REAL two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry
257 TH2D ** fhRePtNCellAsymCutsSM0 ; //![fNPtCuts*fNAsymCuts*fNCellNCuts] REAL two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry
258 TH2D ** fhRePtNCellAsymCutsSM1 ; //![fNPtCuts*fNAsymCuts*fNCellNCuts] REAL two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry
259 TH2D ** fhRePtNCellAsymCutsSM2 ; //![fNPtCuts*fNAsymCuts*fNCellNCuts] REAL two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry
260 TH2D ** fhRePtNCellAsymCutsSM3 ; //![fNPtCuts*fNAsymCuts*fNCellNCuts] REAL two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry
261 TH2D ** fhMiPtNCellAsymCuts ; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Mixed two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry
262 TH2D ** fhRePIDBits ; //![fNPIDBits] REAL two-photon invariant mass distribution for different PID bits
263 TH3D ** fhRePtMult ; //![fNAsymCuts] REAL two-photon invariant mass distribution for different track multiplicity and assymetry cuts
264
265 // Asymmetry vs pt, in pi0/eta regions
266 TH2D * fhRePtAsym ; //! REAL two-photon pt vs asymmetry
267 TH2D * fhRePtAsymPi0 ; //! REAL two-photon pt vs asymmetry, close to pi0 mass
268 TH2D * fhRePtAsymEta ; //! REAL two-photon pt vs asymmetry, close to eta mass
269
270 TH3D * fhEvents; //! Number of events per centrality, RP, zbin
20218aea 271 TH1D * fhCentrality; //! Histogram with centrality bins with at least one pare
272 TH1D * fhCentralityNoPair; //! Histogram with centrality bins with no pair
273
6175da48 274 // Pair opening angle
275 TH2D * fhRealOpeningAngle ; //! Opening angle of pair versus pair energy
276 TH2D * fhRealCosOpeningAngle ; //! Cosinus of opening angle of pair version pair energy
277 TH2D * fhMixedOpeningAngle ; //! Opening angle of pair versus pair energy
278 TH2D * fhMixedCosOpeningAngle ; //! Cosinus of opening angle of pair version pair energy
279
280 //MC analysis histograms
281 //Pi0 Acceptance
282 TH1D * fhPrimPi0Pt ; //! Spectrum of Primary
283 TH1D * fhPrimPi0AccPt ; //! Spectrum of primary with accepted daughters
156549ae 284 TH2D * fhPrimPi0Y ; //! Rapidity distribution of primary particles vs pT
285 TH2D * fhPrimPi0AccY ; //! Rapidity distribution of primary with accepted daughters vs pT
286 TH2D * fhPrimPi0Phi ; //! Azimutal distribution of primary particles vs pT
287 TH2D * fhPrimPi0AccPhi; //! Azimutal distribution of primary with accepted daughters vs pT
6175da48 288 TH2D * fhPrimPi0OpeningAngle ; //! Opening angle of pair versus pair energy, primaries
289 TH2D * fhPrimPi0CosOpeningAngle ; //! Cosinus of opening angle of pair version pair energy, primaries
290 //Eta acceptance
291 TH1D * fhPrimEtaPt ; //! Spectrum of Primary
292 TH1D * fhPrimEtaAccPt ; //! Spectrum of primary with accepted daughters
156549ae 293 TH2D * fhPrimEtaY ; //! Rapidity distribution of primary particles vs pT
294 TH2D * fhPrimEtaAccY ; //! Rapidity distribution of primary with accepted daughters vs pT
295 TH2D * fhPrimEtaPhi ; //! Azimutal distribution of primary particles vs pT
296 TH2D * fhPrimEtaAccPhi; //! Azimutal distribution of primary with accepted daughters vs pT
6175da48 297
08a56f5f 298 // Primaries origin
299 TH2D * fhPrimPi0PtOrigin ; //! Spectrum of generated pi0 vs mother
300 TH2D * fhPrimEtaPtOrigin ; //! Spectrum of generated eta vs mother
301
6175da48 302 //Pair origin
303 //Array of histograms ordered as follows: 0-Photon, 1-electron, 2-pi0, 3-eta, 4-a-proton, 5-a-neutron, 6-stable particles,
304 // 7-other decays, 8-string, 9-final parton, 10-initial parton, intermediate, 11-colliding proton, 12-unrelated
305 TH2D * fhMCOrgMass[13]; //! Mass vs pt of real pairs, check common origin of pair
306 TH2D * fhMCOrgAsym[13]; //! Asymmetry vs pt of real pairs, check common origin of pair
307 TH2D * fhMCOrgDeltaEta[13]; //! Delta Eta vs pt of real pairs, check common origin of pair
308 TH2D * fhMCOrgDeltaPhi[13]; //! Delta Phi vs pt of real pairs, check common origin of pair
309
310 //Multiple cuts in simulation, origin pi0 or eta
311 TH2D ** fhMCPi0MassPtRec; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real pi0 pairs, reconstructed mass vs reconstructed pt of original pair
312 TH2D ** fhMCPi0MassPtTrue; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real pi0 pairs, reconstructed mass vs generated pt of original pair
313 TH2D ** fhMCPi0PtTruePtRec; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real pi0 pairs, reconstructed pt vs generated pt of pair
314 TH2D ** fhMCEtaMassPtRec; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real eta pairs, reconstructed mass vs reconstructed pt of original pair
315 TH2D ** fhMCEtaMassPtTrue; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real eta pairs, reconstructed mass vs generated pt of original pair
316 TH2D ** fhMCEtaPtTruePtRec; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real eta pairs, reconstructed pt vs generated pt of pair
08a56f5f 317
318 TH2D * fhMCPi0PtOrigin ; //! Mass of reoconstructed pi0 pairs in calorimeter vs mother
319 TH2D * fhMCEtaPtOrigin ; //! Mass of reoconstructed pi0 pairs in calorimeter vs mother
6175da48 320
20218aea 321 ClassDef(AliAnaPi0,18)
1c5acb87 322} ;
323
324
325#endif //ALIANAPI0_H
326
327
328