]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGGA/CaloTrackCorrelations/AliAnaPi0.h
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Merge branch 'feature-movesplit'
[u/mrichter/AliRoot.git] / PWGGA / CaloTrackCorrelations / 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 */
1c5acb87 5
6//_________________________________________________________________________
6175da48 7// Class to fill two-photon invariant mass histograms
1c5acb87 8// to be used to extract pi0 raw yield.
6175da48 9// Input is produced by AliAnaPhoton (or any other analysis producing output AliAODPWG4Particles),
10// it will do nothing if executed alone
1c5acb87 11//
12//-- Author: Dmitri Peressounko (RRC "KI")
745913ae 13//-- Adapted to CaloTrackCorr frame by Lamia Benhabib (SUBATECH)
1c5acb87 14//-- and Gustavo Conesa (INFN-Frascati)
15
16//Root
17class TList;
0333ede6 18class TH3F ;
19class TH2F ;
0c1383b5 20class TObjString;
1c5acb87 21
22//Analysis
745913ae 23#include "AliAnaCaloTrackCorrBaseClass.h"
1c5acb87 24class AliAODEvent ;
25class AliESDEvent ;
c8fe2783 26class AliAODPWG4Particle ;
1c5acb87 27
745913ae 28class AliAnaPi0 : public AliAnaCaloTrackCorrBaseClass {
6639984f 29
78219bac 30 public:
6639984f 31 AliAnaPi0() ; // default ctor
6639984f 32 virtual ~AliAnaPi0() ;//virtual dtor
33
6175da48 34 //-------------------------------
35 // General analysis frame methods
36 //-------------------------------
37
0c1383b5 38 TObjString * GetAnalysisCuts();
0333ede6 39
0c1383b5 40 TList * GetCreateOutputObjects();
0333ede6 41
42 void Print(const Option_t * opt) const;
43
44 void MakeAnalysisFillHistograms();
45
46 void InitParameters();
c8fe2783 47
6175da48 48 //-------------------------------
49 // EVENT Bin Methods
50 //-------------------------------
51
2e876d85 52 Int_t GetEventIndex(AliAODPWG4Particle * part, Double_t * vert) ;
156549ae 53
6175da48 54 //-------------------------------
55 //Opening angle pair selection
56 //-------------------------------
0333ede6 57 void SwitchOnAngleSelection() { fUseAngleCut = kTRUE ; }
58 void SwitchOffAngleSelection() { fUseAngleCut = kFALSE ; }
59
60 void SwitchOnAngleEDepSelection() { fUseAngleEDepCut = kTRUE ; }
61 void SwitchOffAngleEDepSelection() { fUseAngleEDepCut = kFALSE ; }
62
63 void SetAngleCut(Float_t a) { fAngleCut = a ; }
64 void SetAngleMaxCut(Float_t a) { fAngleMaxCut = a ; }
6175da48 65
7a972c0c 66 void SwitchOnFillAngleHisto() { fFillAngleHisto = kTRUE ; }
67 void SwitchOffFillAngleHisto() { fFillAngleHisto = kFALSE ; }
68
6175da48 69 //------------------------------------------
20218aea 70 //Do analysis only with clusters in same SM or different combinations of SM
6175da48 71 //------------------------------------------
0333ede6 72 void SwitchOnSameSM() { fSameSM = kTRUE ; }
73 void SwitchOffSameSM() { fSameSM = kFALSE ; }
6175da48 74
0333ede6 75 void SwitchOnSMCombinations() { fFillSMCombinations = kTRUE ; }
76 void SwitchOffSMCombinations() { fFillSMCombinations = kFALSE ; }
20218aea 77
6175da48 78 //-------------------------------
79 //Histogram filling options off by default
80 //-------------------------------
0333ede6 81 void SwitchOnInvPtWeight() { fMakeInvPtPlots = kTRUE ; }
82 void SwitchOffInvPtWeight() { fMakeInvPtPlots = kFALSE ; }
af7b3903 83
0333ede6 84 void SwitchOnFillBadDistHisto() { fFillBadDistHisto = kTRUE ; }
85 void SwitchOffFillBadDistHisto() { fFillBadDistHisto = kFALSE ; }
6175da48 86
87 //-------------------------------------------
88 //Cuts for multiple analysis, off by default
89 //-------------------------------------------
0333ede6 90 void SwitchOnMultipleCutAnalysis() { fMultiCutAna = kTRUE ; }
91 void SwitchOffMultipleCutAnalysis() { fMultiCutAna = kFALSE ; }
5ae09196 92
0333ede6 93 void SetNPtCuts (Int_t s) { if(s <= 10)fNPtCuts = s ; }
94 void SetNAsymCuts (Int_t s) { if(s <= 10)fNAsymCuts = s ; }
95 void SetNNCellCuts(Int_t s) { if(s <= 10)fNCellNCuts = s ; }
96 void SetNPIDBits (Int_t s) { if(s <= 10)fNPIDBits = s ; }
d7c10d78 97
0333ede6 98 void SetPtCutsAt (Int_t p,Float_t v) { if(p < 10)fPtCuts[p] = v ; }
99 void SetAsymCutsAt(Int_t p,Float_t v) { if(p < 10)fAsymCuts[p] = v ; }
100 void SetNCellCutsAt(Int_t p,Int_t v) { if(p < 10)fCellNCuts[p]= v ; }
101 void SetPIDBitsAt (Int_t p,Int_t v) { if(p < 10)fPIDBits[p] = v ; }
d7c10d78 102
7a972c0c 103 void SwitchOnFillSSCombinations() { fFillSSCombinations = kTRUE ; }
104 void SwitchOffFillSSCombinations() { fFillSSCombinations = kFALSE ; }
105
106 void SwitchOnFillAsymmetryHisto() { fFillAsymmetryHisto = kTRUE ; }
107 void SwitchOffFillAsymmetryHisto() { fFillAsymmetryHisto = kFALSE ; }
108
020e681b 109 void SwitchOnFillOriginHisto() { fFillOriginHisto = kTRUE ; }
110 void SwitchOffFillOriginHisto() { fFillOriginHisto = kFALSE ; }
be894c1d 111
29555e96 112 void SwitchOnFillArmenterosThetaStarHisto() { fFillArmenterosThetaStar = kTRUE ; }
113 void SwitchOffFillArmenterosThetaStarHisto() { fFillArmenterosThetaStar = kFALSE ; }
7a972c0c 114
6175da48 115 //MC analysis related methods
0333ede6 116
117 void SwitchOnConversionChecker() { fCheckConversion = kTRUE ; }
118 void SwitchOffConversionChecker() { fCheckConversion = kFALSE ; }
119
120 void SwitchOnMultipleCutAnalysisInSimulation() { fMultiCutAnaSim = kTRUE ; }
121 void SwitchOffMultipleCutAnalysisInSimulation() { fMultiCutAnaSim = kFALSE ; }
20218aea 122
07fc077e 123 void SwitchOnCheckAcceptanceInSector() { fCheckAccInSector = kTRUE ; }
124 void SwitchOffCheckAcceptanceInSector(){ fCheckAccInSector = kFALSE ; }
125
0333ede6 126 void FillAcceptanceHistograms();
b94e038e 127 void FillMCVersusRecDataHistograms(Int_t index1, Int_t index2,
128 Float_t pt1, Float_t pt2,
129 Int_t ncells1, Int_t ncells2,
130 Double_t mass, Double_t pt, Double_t asym,
131 Double_t deta, Double_t dphi);
132
1a8c88c1 133 void FillArmenterosThetaStar(Int_t pdg);
be894c1d 134
135
6639984f 136 private:
0333ede6 137
0333ede6 138 TList ** fEventsList ; //![GetNCentrBin()*GetNZvertBin()*GetNRPBin()] Containers for photons in stored events
139
0333ede6 140 Int_t fNModules ; // Number of EMCAL/PHOS modules, set as many histogras as modules
141
142 Bool_t fUseAngleCut ; // Select pairs depending on their opening angle
143 Bool_t fUseAngleEDepCut ; // Select pairs depending on their opening angle
144 Float_t fAngleCut ; // Select pairs with opening angle larger than a threshold
145 Float_t fAngleMaxCut ; // Select pairs with opening angle smaller than a threshold
6175da48 146
147 //Multiple cuts analysis
0333ede6 148 Bool_t fMultiCutAna; // Do analysis with several or fixed cut
149 Bool_t fMultiCutAnaSim; // Do analysis with several or fixed cut, in the simulation related part
150 Int_t fNPtCuts; // Number of pt cuts
151 Float_t fPtCuts[10]; // Array with different pt cuts
152 Int_t fNAsymCuts; // Number of assymmetry cuts
153 Float_t fAsymCuts[10]; // Array with different assymetry cuts
154 Int_t fNCellNCuts; // Number of cuts with number of cells in cluster
155 Int_t fCellNCuts[10]; // Array with different cell number cluster cuts
156 Int_t fNPIDBits ; // Number of possible PID bit combinations
157 Int_t fPIDBits[10]; // Array with different PID bits
6175da48 158
159 //Switchs of different analysis options
0333ede6 160 Bool_t fMakeInvPtPlots; // D plots with inverse pt weight
161 Bool_t fSameSM; // Select only pairs in same SM;
162 Bool_t fFillSMCombinations; // Fill histograms with different cluster pairs in SM combinations
163 Bool_t fCheckConversion; // Fill histograms with tagged photons as conversion
0333ede6 164 Bool_t fFillBadDistHisto; // Do plots for different distances to bad channels
7a972c0c 165 Bool_t fFillSSCombinations; // Do invariant mass for different combination of shower shape clusters
166 Bool_t fFillAngleHisto; // Fill histograms with pair opening angle
167 Bool_t fFillAsymmetryHisto; // Fill histograms with asymmetry vs pt
020e681b 168 Bool_t fFillOriginHisto; // Fill histograms depending on their origin
29555e96 169 Bool_t fFillArmenterosThetaStar; // Fill armenteros histograms
be894c1d 170
07fc077e 171 Bool_t fCheckAccInSector; // Check that the decay pi0 falls in the same SM or sector
172
1a8c88c1 173 TLorentzVector fPhotonMom1; //! photon cluster momentum
174 TLorentzVector fPhotonMom1Boost; //! photon cluster momentum
175 TLorentzVector fPhotonMom2; //! photon cluster momentum
176 TLorentzVector fPi0Mom; //! pi0 cluster momentum
177 TVector3 fProdVertex; //! production vertex
178
6639984f 179 //Histograms
6175da48 180
156549ae 181 //Event characterization
0333ede6 182 TH1F * fhAverTotECluster; //! Average number of clusters in SM
183 TH1F * fhAverTotECell; //! Average number of cells in SM
184 TH2F * fhAverTotECellvsCluster; //! Average number of cells in SM
185 TH1F * fhEDensityCluster; //! Deposited energy in event per cluster
186 TH1F * fhEDensityCell; //! Deposited energy in event per cell vs cluster
187 TH2F * fhEDensityCellvsCluster; //! Deposited energy in event per cell vs cluster
156549ae 188
7a972c0c 189 TH2F ** fhReMod ; //![fNModules] REAL two-photon invariant mass distribution for different calorimeter modules.
190 TH2F ** fhReSameSideEMCALMod ; //![fNModules-2] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules.
191 TH2F ** fhReSameSectorEMCALMod ; //![fNModules/2] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules.
192 TH2F ** fhReDiffPHOSMod ; //![fNModules] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules.
193 TH2F ** fhMiMod ; //![fNModules] MIXED two-photon invariant mass distribution for different calorimeter modules.
194 TH2F ** fhMiSameSideEMCALMod ; //![fNModules-2] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules.
195 TH2F ** fhMiSameSectorEMCALMod ; //![fNModules/2] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules.
196 TH2F ** fhMiDiffPHOSMod ; //![fNModules-1] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules.
6175da48 197
198 // Pairs with at least one cluster tagged as conversion
7a972c0c 199 TH2F * fhReConv ; //! REAL two-photon invariant mass distribution one of the pair was 2 clusters with small mass
200 TH2F * fhMiConv ; //! MIXED two-photon invariant mass distribution one of the pair was 2 clusters with small mass
201 TH2F * fhReConv2 ; //! REAL two-photon invariant mass distribution both pair photons recombined from 2 clusters with small mass
202 TH2F * fhMiConv2 ; //! MIXED two-photon invariant mass distribution both pair photons recombined from 2 clusters with small mass
0333ede6 203
7a972c0c 204 TH2F ** fhRe1 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry
205 TH2F ** fhMi1 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry
206 TH2F ** fhRe2 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry
207 TH2F ** fhMi2 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry
208 TH2F ** fhRe3 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry
209 TH2F ** fhMi3 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry
6175da48 210
6175da48 211 //Histograms weighted by inverse pT
0333ede6 212 TH2F ** fhReInvPt1 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT
213 TH2F ** fhMiInvPt1 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT
214 TH2F ** fhReInvPt2 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT
215 TH2F ** fhMiInvPt2 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT
216 TH2F ** fhReInvPt3 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT
217 TH2F ** fhMiInvPt3 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT
6175da48 218
219 //Multiple cuts: Assymmetry, pt, n cells, PID
0333ede6 220 TH2F ** fhRePtNCellAsymCuts ; //![fNPtCuts*fNAsymCuts*fNCellNCuts*] REAL two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry
221 TH2F ** fhMiPtNCellAsymCuts ; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Mixed two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry
99b8e903 222 TH2F ** fhRePtNCellAsymCutsSM[12] ; //![fNPtCuts*fNAsymCuts*fNCellNCutsfNModules] REAL two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry for each module
0333ede6 223
224 TH2F ** fhRePIDBits ; //![fNPIDBits] REAL two-photon invariant mass distribution for different PID bits
225 TH3F ** fhRePtMult ; //![fNAsymCuts] REAL two-photon invariant mass distribution for different track multiplicity and assymetry cuts
226 TH2F * fhReSS[3] ; //! Combine clusters with 3 different cuts on shower shape
6175da48 227
228 // Asymmetry vs pt, in pi0/eta regions
0333ede6 229 TH2F * fhRePtAsym ; //! REAL two-photon pt vs asymmetry
230 TH2F * fhRePtAsymPi0 ; //! REAL two-photon pt vs asymmetry, close to pi0 mass
231 TH2F * fhRePtAsymEta ; //! REAL two-photon pt vs asymmetry, close to eta mass
c4a7d28a 232
72542aba 233 //Centrality, Event plane bins
2e876d85 234 TH1I * fhEventBin; //! Number of real pairs in a particular bin (cen,vz,rp)
235 TH1I * fhEventMixBin; //! Number of mixed pairs in a particular bin (cen,vz,rp)
0333ede6 236 TH1F * fhCentrality; //! Histogram with centrality bins with at least one pare
237 TH1F * fhCentralityNoPair; //! Histogram with centrality bins with no pair
20218aea 238
0333ede6 239 TH2F * fhEventPlaneResolution; //! Histogram with Event plane resolution vs centrality
72542aba 240
6175da48 241 // Pair opening angle
0333ede6 242 TH2F * fhRealOpeningAngle ; //! Opening angle of pair versus pair energy
243 TH2F * fhRealCosOpeningAngle ; //! Cosinus of opening angle of pair version pair energy
244 TH2F * fhMixedOpeningAngle ; //! Opening angle of pair versus pair energy
245 TH2F * fhMixedCosOpeningAngle ; //! Cosinus of opening angle of pair version pair energy
6175da48 246
247 //MC analysis histograms
248 //Pi0 Acceptance
29250849 249 TH1F * fhPrimPi0E ; //! Spectrum of Primary
250 TH1F * fhPrimPi0Pt ; //! Spectrum of Primary
251 TH1F * fhPrimPi0AccE ; //! Spectrum of primary with accepted daughters
252 TH1F * fhPrimPi0AccPt ; //! Spectrum of primary with accepted daughters
0333ede6 253 TH2F * fhPrimPi0Y ; //! Rapidity distribution of primary particles vs pT
254 TH2F * fhPrimPi0AccY ; //! Rapidity distribution of primary with accepted daughters vs pT
fba46696 255 TH2F * fhPrimPi0Yeta ; //! PseudoRapidity distribution of primary particles vs pT
256 TH2F * fhPrimPi0YetaYcut ; //! PseudoRapidity distribution of primary particles vs pT, Y<1
257 TH2F * fhPrimPi0AccYeta ; //! PseudoRapidity distribution of primary with accepted daughters vs pT
0333ede6 258 TH2F * fhPrimPi0Phi ; //! Azimutal distribution of primary particles vs pT
259 TH2F * fhPrimPi0AccPhi; //! Azimutal distribution of primary with accepted daughters vs pT
260 TH2F * fhPrimPi0OpeningAngle ; //! Opening angle of pair versus pair energy, primaries
3eb6ab95 261 TH2F * fhPrimPi0OpeningAngleAsym ; //! Opening angle of pair versus pair E asymmetry, pi0 primaries
262 TH2F * fhPrimPi0CosOpeningAngle ; //! Cosinus of opening angle of pair version pair energy, pi0 primaries
c8710850 263 TH2F * fhPrimPi0PtCentrality ; //! primary pi0 reconstructed centrality vs pT
264 TH2F * fhPrimPi0PtEventPlane ; //! primary pi0 reconstructed event plane vs pT
265 TH2F * fhPrimPi0AccPtCentrality ; //! primary pi0 with accepted daughters reconstructed centrality vs pT
266 TH2F * fhPrimPi0AccPtEventPlane ; //! primary pi0 with accepted daughters reconstructed event plane vs pT
267
6175da48 268 //Eta acceptance
29250849 269 TH1F * fhPrimEtaE ; //! Spectrum of Primary
270 TH1F * fhPrimEtaPt ; //! Spectrum of Primary
271 TH1F * fhPrimEtaAccE ; //! Spectrum of primary with accepted daughters
272 TH1F * fhPrimEtaAccPt ; //! Spectrum of primary with accepted daughters
0333ede6 273 TH2F * fhPrimEtaY ; //! Rapidity distribution of primary particles vs pT
274 TH2F * fhPrimEtaAccY ; //! Rapidity distribution of primary with accepted daughters vs pT
fba46696 275 TH2F * fhPrimEtaYeta ; //! PseudoRapidity distribution of primary particles vs pT
276 TH2F * fhPrimEtaYetaYcut ; //! PseudoRapidity distribution of primary particles vs pT, Y<1
277 TH2F * fhPrimEtaAccYeta ; //! PseudoRapidity distribution of primary with accepted daughters vs pT
0333ede6 278 TH2F * fhPrimEtaPhi ; //! Azimutal distribution of primary particles vs pT
279 TH2F * fhPrimEtaAccPhi; //! Azimutal distribution of primary with accepted daughters vs pT
c8710850 280 TH2F * fhPrimEtaOpeningAngle ; //! Opening angle of pair versus pair energy, eta primaries
281 TH2F * fhPrimEtaOpeningAngleAsym ; //! Opening angle of pair versus pair E asymmetry, eta primaries
282 TH2F * fhPrimEtaCosOpeningAngle ; //! Cosinus of opening angle of pair version pair energy, eta primaries
283 TH2F * fhPrimEtaPtCentrality ; //! primary eta reconstructed centrality vs pT
284 TH2F * fhPrimEtaPtEventPlane ; //! primary eta reconstructed event plane vs pT
285 TH2F * fhPrimEtaAccPtCentrality ; //! primary eta with accepted daughters reconstructed centrality vs pT
286 TH2F * fhPrimEtaAccPtEventPlane ; //! primary eta with accepted daughters reconstructed event plane vs pT
6175da48 287
08a56f5f 288 // Primaries origin
0333ede6 289 TH2F * fhPrimPi0PtOrigin ; //! Spectrum of generated pi0 vs mother
290 TH2F * fhPrimEtaPtOrigin ; //! Spectrum of generated eta vs mother
08a56f5f 291
6175da48 292 //Pair origin
293 //Array of histograms ordered as follows: 0-Photon, 1-electron, 2-pi0, 3-eta, 4-a-proton, 5-a-neutron, 6-stable particles,
294 // 7-other decays, 8-string, 9-final parton, 10-initial parton, intermediate, 11-colliding proton, 12-unrelated
0333ede6 295 TH2F * fhMCOrgMass[13]; //! Mass vs pt of real pairs, check common origin of pair
296 TH2F * fhMCOrgAsym[13]; //! Asymmetry vs pt of real pairs, check common origin of pair
297 TH2F * fhMCOrgDeltaEta[13]; //! Delta Eta vs pt of real pairs, check common origin of pair
298 TH2F * fhMCOrgDeltaPhi[13]; //! Delta Phi vs pt of real pairs, check common origin of pair
6175da48 299
300 //Multiple cuts in simulation, origin pi0 or eta
0333ede6 301 TH2F ** fhMCPi0MassPtRec; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real pi0 pairs, reconstructed mass vs reconstructed pt of original pair
302 TH2F ** fhMCPi0MassPtTrue; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real pi0 pairs, reconstructed mass vs generated pt of original pair
303 TH2F ** fhMCPi0PtTruePtRec; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real pi0 pairs, reconstructed pt vs generated pt of pair
304 TH2F ** fhMCEtaMassPtRec; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real eta pairs, reconstructed mass vs reconstructed pt of original pair
305 TH2F ** fhMCEtaMassPtTrue; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real eta pairs, reconstructed mass vs generated pt of original pair
306 TH2F ** fhMCEtaPtTruePtRec; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real eta pairs, reconstructed pt vs generated pt of pair
307
308 TH2F * fhMCPi0PtOrigin ; //! Mass of reoconstructed pi0 pairs in calorimeter vs mother
309 TH2F * fhMCEtaPtOrigin ; //! Mass of reoconstructed pi0 pairs in calorimeter vs mother
6175da48 310
59b85683 311 TH2F * fhMCPi0ProdVertex; //! Spectrum of selected pi0 vs production vertex
312 TH2F * fhMCEtaProdVertex; //! Spectrum of selected eta vs production vertex
00ae36c1 313 TH2F * fhPrimPi0ProdVertex; //! Spectrum of primary pi0 vs production vertex
314 TH2F * fhPrimEtaProdVertex; //! Spectrum of primary eta vs production vertex
315
99b8e903 316 TH2F * fhReMCFromConversion ; //! Invariant mass of 2 clusters originated in conversions
317 TH2F * fhReMCFromNotConversion ; //! Invariant mass of 2 clusters not originated in conversions
318 TH2F * fhReMCFromMixConversion ; //! Invariant mass of 2 clusters one from conversion and the other not
319
fba46696 320 TH2F * fhArmPrimPi0[4]; //! Armenteros plots for primary pi0 in 6 energy bins
321 TH2F * fhArmPrimEta[4]; //! Armenteros plots for primary eta in 6 energy bins
322 TH2F * fhCosThStarPrimPi0; //! cos(theta*) plots vs E for primary pi0, same as asymmetry ...
323 TH2F * fhCosThStarPrimEta; //! cos(theta*) plots vs E for primary eta, same as asymmetry ...
be894c1d 324
5cd6b1a6 325 TH2F * fhEPairDiffTime; //! E pair vs Pair of clusters time difference vs E
326
745913ae 327 AliAnaPi0( const AliAnaPi0 & api0) ; // cpy ctor
328 AliAnaPi0 & operator = (const AliAnaPi0 & api0) ; // cpy assignment
99b8e903 329
1a8c88c1 330 ClassDef(AliAnaPi0,29)
1c5acb87 331} ;
332
333
334#endif //ALIANAPI0_H
335
336
337