]>
Commit | Line | Data |
---|---|---|
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 | |
17 | class TList; | |
0333ede6 | 18 | class TH3F ; |
19 | class TH2F ; | |
0c1383b5 | 20 | class TObjString; |
1c5acb87 | 21 | |
22 | //Analysis | |
745913ae | 23 | #include "AliAnaCaloTrackCorrBaseClass.h" |
1c5acb87 | 24 | class AliAODEvent ; |
25 | class AliESDEvent ; | |
c8fe2783 | 26 | class AliAODPWG4Particle ; |
1c5acb87 | 27 | |
745913ae | 28 | class 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 | ||
133 | void FillArmenterosThetaStar(Int_t pdg, TLorentzVector meson, | |
134 | TLorentzVector daugh1, TLorentzVector daugh2); | |
be894c1d | 135 | |
136 | ||
6639984f | 137 | private: |
0333ede6 | 138 | |
0333ede6 | 139 | TList ** fEventsList ; //![GetNCentrBin()*GetNZvertBin()*GetNRPBin()] Containers for photons in stored events |
140 | ||
0333ede6 | 141 | Int_t fNModules ; // Number of EMCAL/PHOS modules, set as many histogras as modules |
142 | ||
143 | Bool_t fUseAngleCut ; // Select pairs depending on their opening angle | |
144 | Bool_t fUseAngleEDepCut ; // Select pairs depending on their opening angle | |
145 | Float_t fAngleCut ; // Select pairs with opening angle larger than a threshold | |
146 | Float_t fAngleMaxCut ; // Select pairs with opening angle smaller than a threshold | |
6175da48 | 147 | |
148 | //Multiple cuts analysis | |
0333ede6 | 149 | Bool_t fMultiCutAna; // Do analysis with several or fixed cut |
150 | Bool_t fMultiCutAnaSim; // Do analysis with several or fixed cut, in the simulation related part | |
151 | Int_t fNPtCuts; // Number of pt cuts | |
152 | Float_t fPtCuts[10]; // Array with different pt cuts | |
153 | Int_t fNAsymCuts; // Number of assymmetry cuts | |
154 | Float_t fAsymCuts[10]; // Array with different assymetry cuts | |
155 | Int_t fNCellNCuts; // Number of cuts with number of cells in cluster | |
156 | Int_t fCellNCuts[10]; // Array with different cell number cluster cuts | |
157 | Int_t fNPIDBits ; // Number of possible PID bit combinations | |
158 | Int_t fPIDBits[10]; // Array with different PID bits | |
6175da48 | 159 | |
160 | //Switchs of different analysis options | |
0333ede6 | 161 | Bool_t fMakeInvPtPlots; // D plots with inverse pt weight |
162 | Bool_t fSameSM; // Select only pairs in same SM; | |
163 | Bool_t fFillSMCombinations; // Fill histograms with different cluster pairs in SM combinations | |
164 | Bool_t fCheckConversion; // Fill histograms with tagged photons as conversion | |
0333ede6 | 165 | Bool_t fFillBadDistHisto; // Do plots for different distances to bad channels |
7a972c0c | 166 | Bool_t fFillSSCombinations; // Do invariant mass for different combination of shower shape clusters |
167 | Bool_t fFillAngleHisto; // Fill histograms with pair opening angle | |
168 | Bool_t fFillAsymmetryHisto; // Fill histograms with asymmetry vs pt | |
020e681b | 169 | Bool_t fFillOriginHisto; // Fill histograms depending on their origin |
29555e96 | 170 | Bool_t fFillArmenterosThetaStar; // Fill armenteros histograms |
be894c1d | 171 | |
07fc077e | 172 | Bool_t fCheckAccInSector; // Check that the decay pi0 falls in the same SM or sector |
173 | ||
6639984f | 174 | //Histograms |
6175da48 | 175 | |
156549ae | 176 | //Event characterization |
0333ede6 | 177 | TH1F * fhAverTotECluster; //! Average number of clusters in SM |
178 | TH1F * fhAverTotECell; //! Average number of cells in SM | |
179 | TH2F * fhAverTotECellvsCluster; //! Average number of cells in SM | |
180 | TH1F * fhEDensityCluster; //! Deposited energy in event per cluster | |
181 | TH1F * fhEDensityCell; //! Deposited energy in event per cell vs cluster | |
182 | TH2F * fhEDensityCellvsCluster; //! Deposited energy in event per cell vs cluster | |
156549ae | 183 | |
7a972c0c | 184 | TH2F ** fhReMod ; //![fNModules] REAL two-photon invariant mass distribution for different calorimeter modules. |
185 | TH2F ** fhReSameSideEMCALMod ; //![fNModules-2] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules. | |
186 | TH2F ** fhReSameSectorEMCALMod ; //![fNModules/2] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules. | |
187 | TH2F ** fhReDiffPHOSMod ; //![fNModules] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules. | |
188 | TH2F ** fhMiMod ; //![fNModules] MIXED two-photon invariant mass distribution for different calorimeter modules. | |
189 | TH2F ** fhMiSameSideEMCALMod ; //![fNModules-2] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules. | |
190 | TH2F ** fhMiSameSectorEMCALMod ; //![fNModules/2] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules. | |
191 | TH2F ** fhMiDiffPHOSMod ; //![fNModules-1] REAL two-photon invariant mass distribution for different clusters in different calorimeter modules. | |
6175da48 | 192 | |
193 | // Pairs with at least one cluster tagged as conversion | |
7a972c0c | 194 | TH2F * fhReConv ; //! REAL two-photon invariant mass distribution one of the pair was 2 clusters with small mass |
195 | TH2F * fhMiConv ; //! MIXED two-photon invariant mass distribution one of the pair was 2 clusters with small mass | |
196 | TH2F * fhReConv2 ; //! REAL two-photon invariant mass distribution both pair photons recombined from 2 clusters with small mass | |
197 | TH2F * fhMiConv2 ; //! MIXED two-photon invariant mass distribution both pair photons recombined from 2 clusters with small mass | |
0333ede6 | 198 | |
7a972c0c | 199 | TH2F ** fhRe1 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry |
200 | TH2F ** fhMi1 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry | |
201 | TH2F ** fhRe2 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry | |
202 | TH2F ** fhMi2 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry | |
203 | TH2F ** fhRe3 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry | |
204 | TH2F ** fhMi3 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry | |
6175da48 | 205 | |
6175da48 | 206 | //Histograms weighted by inverse pT |
0333ede6 | 207 | TH2F ** fhReInvPt1 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT |
208 | TH2F ** fhMiInvPt1 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT | |
209 | TH2F ** fhReInvPt2 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT | |
210 | TH2F ** fhMiInvPt2 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT | |
211 | TH2F ** fhReInvPt3 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] REAL two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT | |
212 | TH2F ** fhMiInvPt3 ; //![GetNCentrBin()*fNPIDBits*fNAsymCuts] MIXED two-photon invariant mass distribution for different centralities and Asymmetry, inverse pT | |
6175da48 | 213 | |
214 | //Multiple cuts: Assymmetry, pt, n cells, PID | |
0333ede6 | 215 | TH2F ** fhRePtNCellAsymCuts ; //![fNPtCuts*fNAsymCuts*fNCellNCuts*] REAL two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry |
216 | TH2F ** fhMiPtNCellAsymCuts ; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Mixed two-photon invariant mass distribution for different pt cut, n cell cuts and assymetry | |
99b8e903 | 217 | 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 | 218 | |
219 | TH2F ** fhRePIDBits ; //![fNPIDBits] REAL two-photon invariant mass distribution for different PID bits | |
220 | TH3F ** fhRePtMult ; //![fNAsymCuts] REAL two-photon invariant mass distribution for different track multiplicity and assymetry cuts | |
221 | TH2F * fhReSS[3] ; //! Combine clusters with 3 different cuts on shower shape | |
6175da48 | 222 | |
223 | // Asymmetry vs pt, in pi0/eta regions | |
0333ede6 | 224 | TH2F * fhRePtAsym ; //! REAL two-photon pt vs asymmetry |
225 | TH2F * fhRePtAsymPi0 ; //! REAL two-photon pt vs asymmetry, close to pi0 mass | |
226 | TH2F * fhRePtAsymEta ; //! REAL two-photon pt vs asymmetry, close to eta mass | |
c4a7d28a | 227 | |
72542aba | 228 | //Centrality, Event plane bins |
2e876d85 | 229 | TH1I * fhEventBin; //! Number of real pairs in a particular bin (cen,vz,rp) |
230 | TH1I * fhEventMixBin; //! Number of mixed pairs in a particular bin (cen,vz,rp) | |
0333ede6 | 231 | TH1F * fhCentrality; //! Histogram with centrality bins with at least one pare |
232 | TH1F * fhCentralityNoPair; //! Histogram with centrality bins with no pair | |
20218aea | 233 | |
0333ede6 | 234 | TH2F * fhEventPlaneResolution; //! Histogram with Event plane resolution vs centrality |
72542aba | 235 | |
6175da48 | 236 | // Pair opening angle |
0333ede6 | 237 | TH2F * fhRealOpeningAngle ; //! Opening angle of pair versus pair energy |
238 | TH2F * fhRealCosOpeningAngle ; //! Cosinus of opening angle of pair version pair energy | |
239 | TH2F * fhMixedOpeningAngle ; //! Opening angle of pair versus pair energy | |
240 | TH2F * fhMixedCosOpeningAngle ; //! Cosinus of opening angle of pair version pair energy | |
6175da48 | 241 | |
242 | //MC analysis histograms | |
243 | //Pi0 Acceptance | |
29250849 | 244 | TH1F * fhPrimPi0E ; //! Spectrum of Primary |
245 | TH1F * fhPrimPi0Pt ; //! Spectrum of Primary | |
246 | TH1F * fhPrimPi0AccE ; //! Spectrum of primary with accepted daughters | |
247 | TH1F * fhPrimPi0AccPt ; //! Spectrum of primary with accepted daughters | |
0333ede6 | 248 | TH2F * fhPrimPi0Y ; //! Rapidity distribution of primary particles vs pT |
249 | TH2F * fhPrimPi0AccY ; //! Rapidity distribution of primary with accepted daughters vs pT | |
fba46696 | 250 | TH2F * fhPrimPi0Yeta ; //! PseudoRapidity distribution of primary particles vs pT |
251 | TH2F * fhPrimPi0YetaYcut ; //! PseudoRapidity distribution of primary particles vs pT, Y<1 | |
252 | TH2F * fhPrimPi0AccYeta ; //! PseudoRapidity distribution of primary with accepted daughters vs pT | |
0333ede6 | 253 | TH2F * fhPrimPi0Phi ; //! Azimutal distribution of primary particles vs pT |
254 | TH2F * fhPrimPi0AccPhi; //! Azimutal distribution of primary with accepted daughters vs pT | |
255 | TH2F * fhPrimPi0OpeningAngle ; //! Opening angle of pair versus pair energy, primaries | |
3eb6ab95 | 256 | TH2F * fhPrimPi0OpeningAngleAsym ; //! Opening angle of pair versus pair E asymmetry, pi0 primaries |
257 | TH2F * fhPrimPi0CosOpeningAngle ; //! Cosinus of opening angle of pair version pair energy, pi0 primaries | |
c8710850 | 258 | TH2F * fhPrimPi0PtCentrality ; //! primary pi0 reconstructed centrality vs pT |
259 | TH2F * fhPrimPi0PtEventPlane ; //! primary pi0 reconstructed event plane vs pT | |
260 | TH2F * fhPrimPi0AccPtCentrality ; //! primary pi0 with accepted daughters reconstructed centrality vs pT | |
261 | TH2F * fhPrimPi0AccPtEventPlane ; //! primary pi0 with accepted daughters reconstructed event plane vs pT | |
262 | ||
6175da48 | 263 | //Eta acceptance |
29250849 | 264 | TH1F * fhPrimEtaE ; //! Spectrum of Primary |
265 | TH1F * fhPrimEtaPt ; //! Spectrum of Primary | |
266 | TH1F * fhPrimEtaAccE ; //! Spectrum of primary with accepted daughters | |
267 | TH1F * fhPrimEtaAccPt ; //! Spectrum of primary with accepted daughters | |
0333ede6 | 268 | TH2F * fhPrimEtaY ; //! Rapidity distribution of primary particles vs pT |
269 | TH2F * fhPrimEtaAccY ; //! Rapidity distribution of primary with accepted daughters vs pT | |
fba46696 | 270 | TH2F * fhPrimEtaYeta ; //! PseudoRapidity distribution of primary particles vs pT |
271 | TH2F * fhPrimEtaYetaYcut ; //! PseudoRapidity distribution of primary particles vs pT, Y<1 | |
272 | TH2F * fhPrimEtaAccYeta ; //! PseudoRapidity distribution of primary with accepted daughters vs pT | |
0333ede6 | 273 | TH2F * fhPrimEtaPhi ; //! Azimutal distribution of primary particles vs pT |
274 | TH2F * fhPrimEtaAccPhi; //! Azimutal distribution of primary with accepted daughters vs pT | |
c8710850 | 275 | TH2F * fhPrimEtaOpeningAngle ; //! Opening angle of pair versus pair energy, eta primaries |
276 | TH2F * fhPrimEtaOpeningAngleAsym ; //! Opening angle of pair versus pair E asymmetry, eta primaries | |
277 | TH2F * fhPrimEtaCosOpeningAngle ; //! Cosinus of opening angle of pair version pair energy, eta primaries | |
278 | TH2F * fhPrimEtaPtCentrality ; //! primary eta reconstructed centrality vs pT | |
279 | TH2F * fhPrimEtaPtEventPlane ; //! primary eta reconstructed event plane vs pT | |
280 | TH2F * fhPrimEtaAccPtCentrality ; //! primary eta with accepted daughters reconstructed centrality vs pT | |
281 | TH2F * fhPrimEtaAccPtEventPlane ; //! primary eta with accepted daughters reconstructed event plane vs pT | |
6175da48 | 282 | |
08a56f5f | 283 | // Primaries origin |
0333ede6 | 284 | TH2F * fhPrimPi0PtOrigin ; //! Spectrum of generated pi0 vs mother |
285 | TH2F * fhPrimEtaPtOrigin ; //! Spectrum of generated eta vs mother | |
08a56f5f | 286 | |
6175da48 | 287 | //Pair origin |
288 | //Array of histograms ordered as follows: 0-Photon, 1-electron, 2-pi0, 3-eta, 4-a-proton, 5-a-neutron, 6-stable particles, | |
289 | // 7-other decays, 8-string, 9-final parton, 10-initial parton, intermediate, 11-colliding proton, 12-unrelated | |
0333ede6 | 290 | TH2F * fhMCOrgMass[13]; //! Mass vs pt of real pairs, check common origin of pair |
291 | TH2F * fhMCOrgAsym[13]; //! Asymmetry vs pt of real pairs, check common origin of pair | |
292 | TH2F * fhMCOrgDeltaEta[13]; //! Delta Eta vs pt of real pairs, check common origin of pair | |
293 | TH2F * fhMCOrgDeltaPhi[13]; //! Delta Phi vs pt of real pairs, check common origin of pair | |
6175da48 | 294 | |
295 | //Multiple cuts in simulation, origin pi0 or eta | |
0333ede6 | 296 | TH2F ** fhMCPi0MassPtRec; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real pi0 pairs, reconstructed mass vs reconstructed pt of original pair |
297 | TH2F ** fhMCPi0MassPtTrue; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real pi0 pairs, reconstructed mass vs generated pt of original pair | |
298 | TH2F ** fhMCPi0PtTruePtRec; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real pi0 pairs, reconstructed pt vs generated pt of pair | |
299 | TH2F ** fhMCEtaMassPtRec; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real eta pairs, reconstructed mass vs reconstructed pt of original pair | |
300 | TH2F ** fhMCEtaMassPtTrue; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real eta pairs, reconstructed mass vs generated pt of original pair | |
301 | TH2F ** fhMCEtaPtTruePtRec; //![fNPtCuts*fNAsymCuts*fNCellNCuts] Real eta pairs, reconstructed pt vs generated pt of pair | |
302 | ||
303 | TH2F * fhMCPi0PtOrigin ; //! Mass of reoconstructed pi0 pairs in calorimeter vs mother | |
304 | TH2F * fhMCEtaPtOrigin ; //! Mass of reoconstructed pi0 pairs in calorimeter vs mother | |
6175da48 | 305 | |
59b85683 | 306 | TH2F * fhMCPi0ProdVertex; //! Spectrum of selected pi0 vs production vertex |
307 | TH2F * fhMCEtaProdVertex; //! Spectrum of selected eta vs production vertex | |
00ae36c1 | 308 | TH2F * fhPrimPi0ProdVertex; //! Spectrum of primary pi0 vs production vertex |
309 | TH2F * fhPrimEtaProdVertex; //! Spectrum of primary eta vs production vertex | |
310 | ||
99b8e903 | 311 | TH2F * fhReMCFromConversion ; //! Invariant mass of 2 clusters originated in conversions |
312 | TH2F * fhReMCFromNotConversion ; //! Invariant mass of 2 clusters not originated in conversions | |
313 | TH2F * fhReMCFromMixConversion ; //! Invariant mass of 2 clusters one from conversion and the other not | |
314 | ||
fba46696 | 315 | TH2F * fhArmPrimPi0[4]; //! Armenteros plots for primary pi0 in 6 energy bins |
316 | TH2F * fhArmPrimEta[4]; //! Armenteros plots for primary eta in 6 energy bins | |
317 | TH2F * fhCosThStarPrimPi0; //! cos(theta*) plots vs E for primary pi0, same as asymmetry ... | |
318 | TH2F * fhCosThStarPrimEta; //! cos(theta*) plots vs E for primary eta, same as asymmetry ... | |
be894c1d | 319 | |
745913ae | 320 | AliAnaPi0( const AliAnaPi0 & api0) ; // cpy ctor |
321 | AliAnaPi0 & operator = (const AliAnaPi0 & api0) ; // cpy assignment | |
99b8e903 | 322 | |
59b85683 | 323 | ClassDef(AliAnaPi0,28) |
1c5acb87 | 324 | } ; |
325 | ||
326 | ||
327 | #endif //ALIANAPI0_H | |
328 | ||
329 | ||
330 |