]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.h
add optional debug prints
[u/mrichter/AliRoot.git] / PWGGA / CaloTrackCorrelations / AliAnaInsideClusterInvariantMass.h
CommitLineData
992b14a7 1#ifndef ALIANAINSIDECLUSTERINVARIANTMASS_H
2#define ALIANAINSIDECLUSTERINVARIANTMASS_H
3/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
992b14a7 5
6//_________________________________________________________________________
7//
8// Split clusters with some criteria and calculate invariant mass
9// to identify them as pi0 or conversion
10//
11//
83351853 12//-- Author: Gustavo Conesa (LPSC-Grenoble)
992b14a7 13//_________________________________________________________________________
14
15
16// --- ROOT system ---
17class TList ;
18class TObjString;
19class TLorentzVector;
20
21// --- ANALYSIS system ---
5c46c992 22class AliAODCaloCluster;
23
745913ae 24#include "AliAnaCaloTrackCorrBaseClass.h"
992b14a7 25
745913ae 26class AliAnaInsideClusterInvariantMass : public AliAnaCaloTrackCorrBaseClass {
992b14a7 27
28 public:
3c1d9afb 29
992b14a7 30 AliAnaInsideClusterInvariantMass() ; // default ctor
31 virtual ~AliAnaInsideClusterInvariantMass() { ; } //virtual dtor
992b14a7 32
36769d30 33 void CheckLocalMaximaMCOrigin(AliVCluster* cluster, const Int_t mcindex, const Int_t noverlaps);
b2e375c7 34
992b14a7 35 TObjString * GetAnalysisCuts();
36
37 TList * GetCreateOutputObjects();
dbe09c26 38
b2e375c7 39 void GetMCIndex(AliVCluster * cluster, Int_t & mcindex);
40
41 void GetMCPrimaryKine(AliVCluster* cluster, const Int_t mcindex, const Bool_t matched,
42 Float_t & eprim, Float_t & asymGen, Int_t & noverlaps );
43
44 void FillAngleHistograms(const Int_t nMax, const Bool_t matched,
45 const Float_t en, const Float_t angle, const Float_t mass);
46
47 void FillEBinHistograms(const Int_t ebin, const Int_t nMax, const Int_t mcindex, const Float_t splitFrac,
48 const Float_t mass, const Float_t asym, const Float_t l0);
49
50 void FillMCHistograms(const Float_t en, const Float_t e1 , const Float_t e2,
51 const Int_t ebin, const Int_t mcindex,
52 const Float_t l0, const Float_t mass,
53 const Int_t nMax, const Bool_t matched,
54 const Float_t splitFrac, const Float_t asym,
55 const Float_t eprim, const Float_t asymGen);
56
83351853 57 void FillMCOverlapHistograms(const Float_t en, const Float_t mass, const Float_t l0, const Float_t splitFrac,
b2e375c7 58 const Int_t nlm, const Int_t ebin, const Bool_t matched,
59 const Int_t mcindex, const Int_t noverlaps);
60
61 void FillSSWeightHistograms(AliVCluster *cluster, const Int_t nlm, const Int_t absId1, const Int_t absId2);
62
63 void FillSSExtraHistograms(AliVCluster *cluster, const Int_t nMax,
64 const Bool_t matched, const Int_t mcindex,
65 const Float_t mass , const Int_t ebin);
66
67 void FillTrackMatchingHistograms(AliVCluster * cluster,const Int_t nMax, const Int_t mcindex);
dbe09c26 68
992b14a7 69 void Init();
70
71 void InitParameters();
72
73 void MakeAnalysisFillHistograms() ;
992b14a7 74
71e3889f 75 void Print(const Option_t * opt) const;
992b14a7 76
71e3889f 77 void SetCalorimeter(TString & det) { fCalorimeter = det ; }
78
79 void SetM02Cut(Float_t min=0, Float_t max=10) { fM02MinCut = min ; fM02MaxCut = max ; }
992b14a7 80
71e3889f 81 void SetMinNCells(Int_t cut) { fMinNCells = cut ; }
2cb134fb 82
83 void SetMinBadChannelDistance(Float_t cut) { fMinBadDist = cut ; }
992b14a7 84
8edbd100 85 void SetWCorrectionParameter(Float_t p = 0.07) { fWSimu = p ; }
86
883411b2 87 void SwitchOnFillAngleHistograms() { fFillAngleHisto = kTRUE ; }
88 void SwitchOffFillAngleHistograms() { fFillAngleHisto = kFALSE ; }
89
90 void SwitchOnFillExtraSSHistograms() { fFillSSExtraHisto = kTRUE ; }
91 void SwitchOffFillExtraSSHistograms() { fFillSSExtraHisto = kFALSE ; }
8e81c2cf 92
8edbd100 93 void SwitchOnFillTMHistograms() { fFillTMHisto = kTRUE ; }
94 void SwitchOffFillTMHistograms() { fFillTMHisto = kFALSE ; }
95
8e81c2cf 96 void SwitchOnFillTMResidualHistograms() { fFillTMResidualHisto = kTRUE ; }
97 void SwitchOffFillTMResidualHistograms() { fFillTMResidualHisto = kFALSE ; }
98
b2e375c7 99 void SwitchOnFillMCPrimaryHistograms() { fFillMCHisto = kTRUE ; }
100 void SwitchOffFillMCPrimaryHistograms() { fFillMCHisto = kFALSE ; }
2a77f6f4 101
19391b8c 102 void SwitchOnFillSSWeightHistograms() { fFillSSWeightHisto = kTRUE ; }
103 void SwitchOffFillSSWeightHistograms() { fFillSSWeightHisto = kFALSE ; }
104
105 void SwitchOnFillEbinHistograms() { fFillEbinHisto = kTRUE ; }
106 void SwitchOffFillEbinHistograms() { fFillEbinHisto = kFALSE ; }
dbe09c26 107
b2e375c7 108 void SwitchOnFillMCOverlapHistograms() { fFillMCOverlapHisto = kTRUE ; }
109 void SwitchOffFillMCOverlapHistograms() { fFillMCOverlapHisto = kFALSE ; }
110
19391b8c 111 void SetNWeightForShowerShape(Int_t n) { fSSWeightN = n ; }
dbe09c26 112 void SetWeightForShowerShape(Int_t i, Float_t v) { if (i < 10) fSSWeight[i] = v ; }
113
19391b8c 114 void SetNECellCutForShowerShape(Int_t n) { fSSECellCutN = n ; }
115 void SetECellCutForShowerShape(Int_t i, Float_t v) { if (i < 10) fSSECellCut[i] = v ; }
116
117
118 void RecalculateClusterShowerShapeParametersWithCellCut(const AliEMCALGeometry * geom, AliVCaloCells* cells, AliVCluster * cluster,
119 Float_t & l0, Float_t & l1,
120 Float_t & disp, Float_t & dEta, Float_t & dPhi,
121 Float_t & sEta, Float_t & sPhi, Float_t & sEtaPhi,Float_t eCellMin = 0.);
122
2a77f6f4 123
992b14a7 124 //For histograms
c5693f62 125 enum mcTypes { kmcPhoton = 1, kmcConversion = 2, kmcPi0 = 3,
de7d73e6 126 kmcEta = 4, kmcElectron = 5, kmcHadron = 6, kmcPi0Conv = 7 };
992b14a7 127
128 private:
129
5c46c992 130 TString fCalorimeter ; // Calorimeter where the gamma is searched
3c1d9afb 131 Float_t fM02MaxCut ; // Study clusters with l0 smaller than cut
132 Float_t fM02MinCut ; // Study clusters with l0 larger than cut
133 Int_t fMinNCells ; // Study clusters with ncells larger than cut
134 Float_t fMinBadDist ; // Minimal distance to bad channel to accept cluster
b2e375c7 135 Float_t fHistoECut ; // Fixed E cut for some histograms
243c2909 136
883411b2 137 Bool_t fFillAngleHisto; // Fill splitted clusters angle histograms
8edbd100 138 Bool_t fFillTMHisto ; // Fill track matching histos,
8e81c2cf 139 Bool_t fFillTMResidualHisto ; // Fill track matching histos, residuals
140 Bool_t fFillSSExtraHisto ; // Fill shower shape extra histos
b2e375c7 141 Bool_t fFillMCHisto ; // Fill MC energy fraction histos
dbe09c26 142 Bool_t fFillSSWeightHisto ; // Fill weigth histograms
19391b8c 143 Bool_t fFillEbinHisto ; // Fill E bin histograms
b2e375c7 144 Bool_t fFillMCOverlapHisto ; // Fill MC particles overlap histograms
8edbd100 145
dbe09c26 146 Float_t fSSWeight[10]; // List of weights to test
19391b8c 147 Int_t fSSWeightN; // Total number of weights to test
148
149 Float_t fSSECellCut[10]; // List of cell min energy cuts to test
150 Int_t fSSECellCutN; // Total number of cell min energy cuts to test
dbe09c26 151
8edbd100 152 Float_t fWSimu; // Slope of the linear correction factor for the shower
153 // shape weight in simulation, about 0.07
154
992b14a7 155 //Histograms
156
de7d73e6 157 TH2F * fhMassNLocMax1[8][2] ; //! Mass of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types
158 TH2F * fhMassNLocMax2[8][2] ; //! Mass of 2 cells local maxima vs E, 1-6 for different MC particle types
159 TH2F * fhMassNLocMaxN[8][2] ; //! Mass of >2 cells local maxima vs E, 1-6 for different MC particle types
5c46c992 160
de7d73e6 161 TH2F * fhAsymNLocMax1[8][2] ; //! Asymmetry of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types
162 TH2F * fhAsymNLocMax2[8][2] ; //! Asymmetry of 2 cells local maxima vs E, 1-6 for different MC particle types
163 TH2F * fhAsymNLocMaxN[8][2] ; //! Asymmetry of >2 cells local maxima vs E, 1-6 for different MC particle types
2a77f6f4 164
e671adc2 165 TH2F * fhSplitEFractionvsAsyNLocMax1[2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1 vs |A|
166 TH2F * fhSplitEFractionvsAsyNLocMax2[2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2 vs |A|
167 TH2F * fhSplitEFractionvsAsyNLocMaxN[2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2 vs |A|
168
77cadd95 169 TH2F * fhMassM02CutNLocMax1 ; //! M02(E) selection, not matched, Mass of split clusters, NLM = 1
170 TH2F * fhMassM02CutNLocMax2 ; //! M02(E) selection, not matched, Mass of split clusters, NLM = 1
171 TH2F * fhMassM02CutNLocMaxN ; //! M02(E) selection, not matched, Mass of split clusters, NLM > 2
2a77f6f4 172
77cadd95 173 TH2F * fhAsymM02CutNLocMax1 ; //! M02(E) selection, not matched, energy asymmetry of split clusters, NLM = 1
174 TH2F * fhAsymM02CutNLocMax2 ; //! M02(E) selection, not matched, energy asymmetry of split clusters, NLM = 2
175 TH2F * fhAsymM02CutNLocMaxN ; //! M02(E) selection, not matched, energy asymmetry of split clusters, NLM > 2
e671adc2 176
77cadd95 177 TH2F * fhMassSplitECutNLocMax1 ; //! 85% of split energy, not matched, Mass of split clusters, NLM = 1
178 TH2F * fhMassSplitECutNLocMax2 ; //! 85% of split energy, not matched, Mass of split clusters, NLM = 1
179 TH2F * fhMassSplitECutNLocMaxN ; //! 85% of split energy, not matched, Mass of split clusters, NLM > 2
667432ef 180
de7d73e6 181 TH2F * fhMassM02NLocMax1[8][2] ; //! Mass of splitted clusters when 1 local max vs M02, for E > 8 GeV, 1-6 for different MC particle types
182 TH2F * fhMassM02NLocMax2[8][2] ; //! Mass of splitted clusters when 2 local max vs M02, for E > 8 GeV, 1-6 for different MC particle types
183 TH2F * fhMassM02NLocMaxN[8][2] ; //! Mass of splitted clusters when >2 local max vs M02, for E > 8 GeV, 1-6 for different MC particle types
0137016b 184
b8eb40fc 185 TH2F * fhMassM02NLocMax1Ebin[4] ; //! Mass of splitted clusters when 1 local max vs M02, 4 E bins, neutral clusters
186 TH2F * fhMassM02NLocMax2Ebin[4] ; //! Mass of splitted clusters when 2 local max vs M02, 4 E bins, neutral clusters
187 TH2F * fhMassM02NLocMaxNEbin[4] ; //! Mass of splitted clusters when >2 local max vs M02, 4 E bins, neutral clusters
188
189 TH2F * fhMassAsyNLocMax1Ebin[4] ; //! Mass of Mass of splitted clusters when 1 local max vs asymmetry, 4 E bins, neutral clusters
190 TH2F * fhMassAsyNLocMax2Ebin[4] ; //! Mass of Mass of splitted clusters when 2 local max vs asymmetry, 4 E bins, neutral clusters
191 TH2F * fhMassAsyNLocMaxNEbin[4] ; //! Mass of Mass of splitted clusters when >2 local max vs asymmetry, 4 E bins, neutral clusters
192
193 TH2F * fhAsyMCGenRecoNLocMax1EbinPi0[4] ; //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when 1 local max, 4 E bins, neutral clusters
194 TH2F * fhAsyMCGenRecoNLocMax2EbinPi0[4] ; //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when 2 local max, 4 E bins, neutral clusters
195 TH2F * fhAsyMCGenRecoNLocMaxNEbinPi0[4] ; //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when >2 local max, 4 E bins, neutral clusters
7b686344 196
de7d73e6 197 TH2F * fhMassDispEtaNLocMax1[8][2] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, for E > 8 GeV, 1-6 for different MC particle types
198 TH2F * fhMassDispEtaNLocMax2[8][2] ; //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
199 TH2F * fhMassDispEtaNLocMaxN[8][2] ; //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
d2655d46 200
e671adc2 201 TH2F * fhMassDispEtaNLocMax1Ebin[4] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
202 TH2F * fhMassDispEtaNLocMax2Ebin[4] ; //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
203 TH2F * fhMassDispEtaNLocMaxNEbin[4] ; //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
d2655d46 204
de7d73e6 205 TH2F * fhMassDispPhiNLocMax1[8][2] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, for E > 8 GeV, 1-6 for different MC particle types
206 TH2F * fhMassDispPhiNLocMax2[8][2] ; //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
207 TH2F * fhMassDispPhiNLocMaxN[8][2] ; //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
d2655d46 208
e671adc2 209 TH2F * fhMassDispPhiNLocMax1Ebin[4] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
210 TH2F * fhMassDispPhiNLocMax2Ebin[4] ; //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
211 TH2F * fhMassDispPhiNLocMaxNEbin[4] ; //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
d2655d46 212
de7d73e6 213 TH2F * fhMassDispAsyNLocMax1[8][2] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, for E > 8 GeV, 1-6 for different MC particle types
214 TH2F * fhMassDispAsyNLocMax2[8][2] ; //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
215 TH2F * fhMassDispAsyNLocMaxN[8][2] ; //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
d2655d46 216
e671adc2 217 TH2F * fhMassDispAsyNLocMax1Ebin[4] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
218 TH2F * fhMassDispAsyNLocMax2Ebin[4] ; //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
219 TH2F * fhMassDispAsyNLocMaxNEbin[4] ; //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
d2655d46 220
de7d73e6 221 TH2F * fhNLocMax [8][2] ; //! Number of maxima in cluster vs E, 1-6 for different MC particle types
222 TH2F * fhNLocMaxM02Cut[8][2] ; //! Number of maxima in cluster vs E, 1-6 for different MC particle types, after SS cut
5c46c992 223
de7d73e6 224 TH2F * fhM02NLocMax1 [8][2] ; //! M02 vs E for N max in cluster = 1, 1-6 for different MC particle types
225 TH2F * fhM02NLocMax2 [8][2] ; //! M02 vs E for N max in cluster = 2, 1-6 for different MC particle types
226 TH2F * fhM02NLocMaxN [8][2] ; //! M02 vs E for N max in cluster > 2, 1-6 for different MC particle types
8e81c2cf 227
e671adc2 228 TH2F * fhMCAsymM02NLocMax1MCPi0Ebin[4] ; //! M02 vs decay asymmetry for N max in cluster = 1, for 4 energy bins
229 TH2F * fhMCAsymM02NLocMax2MCPi0Ebin[4] ; //! M02 vs decay asymmetry for N max in cluster = 2, for 4 energy bins
230 TH2F * fhMCAsymM02NLocMaxNMCPi0Ebin[4] ; //! M02 vs decay asymmetry for N max in cluster > 2, for 4 energy bins
8e81c2cf 231
de7d73e6 232 TH2F * fhMCGenFracNLocMax1[8][2] ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types
233 TH2F * fhMCGenFracNLocMax2[8][2] ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types
234 TH2F * fhMCGenFracNLocMaxN[8][2] ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types
667432ef 235
b8eb40fc 236 TH2F * fhMCGenFracAfterCutsNLocMax1MCPi0 ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, MCPi0 after M02 and asymmetry cut
237 TH2F * fhMCGenFracAfterCutsNLocMax2MCPi0 ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, MCPi0, after M02 and asymmetry cut
238 TH2F * fhMCGenFracAfterCutsNLocMaxNMCPi0 ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, MCPi0, after M02 and asymmetry cut
667432ef 239
de7d73e6 240 TH2F * fhMCGenSplitEFracNLocMax1[8][2] ; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types
241 TH2F * fhMCGenSplitEFracNLocMax2[8][2] ; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types
242 TH2F * fhMCGenSplitEFracNLocMaxN[8][2] ; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types
667432ef 243
244 TH2F * fhMCGenSplitEFracAfterCutsNLocMax1MCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types
245 TH2F * fhMCGenSplitEFracAfterCutsNLocMax2MCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types
246 TH2F * fhMCGenSplitEFracAfterCutsNLocMaxNMCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types
bb2d339b 247
de7d73e6 248 TH2F * fhMCGenEFracvsSplitEFracNLocMax1[8][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster = 1, MC pi0
249 TH2F * fhMCGenEFracvsSplitEFracNLocMax2[8][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster = 2, MC pi0
250 TH2F * fhMCGenEFracvsSplitEFracNLocMaxN[8][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster > 2, MC pi0
bb2d339b 251
de7d73e6 252 TH2F * fhMCGenEvsSplitENLocMax1[8][2] ; //! E generated particle vs E1+E2 for N max in cluster = 1, 1-6 for different MC particle types
253 TH2F * fhMCGenEvsSplitENLocMax2[8][2] ; //! E generated particle vs E1+E2 for N max in cluster = 2, 1-6 for different MC particle types
254 TH2F * fhMCGenEvsSplitENLocMaxN[8][2] ; //! E generated particle vs E1+E2 for N max in cluster > 2, 1-6 for different MC particle types
bb2d339b 255
de7d73e6 256 TH2F * fhMCGenFracNLocMaxEbin[8][4] ; //! NLM vs E generated particle / E reconstructed vs E reconstructed 1-6 for different MC particle types, not matched to track
257 TH2F * fhMCGenFracNLocMaxEbinMatched[8][4] ; //! NLM vs E generated particle / E reconstructed vs E reconstructed 1-6 for different MC particle types, matched to track
53f2c382 258
de7d73e6 259 TH2F * fhM02MCGenFracNLocMax1Ebin[8][4] ; //! M02 vs E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types, not track matched
260 TH2F * fhM02MCGenFracNLocMax2Ebin[8][4] ; //! M02 vs E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types, not track matched
261 TH2F * fhM02MCGenFracNLocMaxNEbin[8][4] ; //! M02 vs E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types, not track matched
8e81c2cf 262
de7d73e6 263 TH2F * fhMassMCGenFracNLocMax1Ebin[8][4] ; //! Mass vs E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types, not track matched
264 TH2F * fhMassMCGenFracNLocMax2Ebin[8][4] ; //! Mass vs E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types, not track matched
265 TH2F * fhMassMCGenFracNLocMaxNEbin[8][4] ; //! Mass vs E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types, not track matched
8e81c2cf 266
de7d73e6 267 TH2F * fhNCellNLocMax1[8][2] ; //! n cells in cluster vs E for N max in cluster = 1, 1-6 for different MC particle types
268 TH2F * fhNCellNLocMax2[8][2] ; //! n cells in cluster vs E for N max in cluster = 2, 1-6 for different MC particle types
269 TH2F * fhNCellNLocMaxN[8][2] ; //! n cells in cluster vs E for N max in cluster > 2, 1-6 for different MC particle types
992b14a7 270
c8710850 271 TH2F * fhM02Pi0NLocMax1[8][2] ; //! M02 for Mass around pi0, N Local Maxima = 1
272 TH2F * fhM02EtaNLocMax1[8][2] ; //! M02 for Mass around eta, N Local Maxima = 1
273 TH2F * fhM02ConNLocMax1[8][2] ; //! M02 for Mass around close to 0, N Local Maxima = 1
274
275 TH2F * fhM02Pi0NLocMax2[8][2] ; //! M02 for Mass around pi0, N Local Maxima = 2
276 TH2F * fhM02EtaNLocMax2[8][2] ; //! M02 for Mass around eta, N Local Maxima = 2
277 TH2F * fhM02ConNLocMax2[8][2] ; //! M02 for Mass around close to 0, N Local Maxima = 2
5c46c992 278
c8710850 279 TH2F * fhM02Pi0NLocMaxN[8][2] ; //! M02 for Mass around pi0, N Local Maxima > 2
280 TH2F * fhM02EtaNLocMaxN[8][2] ; //! M02 for Mass around eta, N Local Maxima > 2
281 TH2F * fhM02ConNLocMaxN[8][2] ; //! M02 for Mass around close to 0, N Local Maxima > 2
fc01318e 282
c8710850 283 TH2F * fhMassPi0NLocMax1[8][2] ; //! Mass for selected pi0, N Local Maxima = 1
284 TH2F * fhMassEtaNLocMax1[8][2] ; //! Mass for selected around eta, N Local Maxima = 1
285 TH2F * fhMassConNLocMax1[8][2] ; //! Mass for selected around close to 0, N Local Maxima = 1
e671adc2 286
c8710850 287 TH2F * fhMassPi0NLocMax2[8][2] ; //! Mass for selected around pi0, N Local Maxima = 2
288 TH2F * fhMassEtaNLocMax2[8][2] ; //! Mass for selected around eta, N Local Maxima = 2
289 TH2F * fhMassConNLocMax2[8][2] ; //! Mass for selected around close to 0, N Local Maxima = 2
bb2d339b 290
c8710850 291 TH2F * fhMassPi0NLocMaxN[8][2] ; //! Mass for selected around pi0, N Local Maxima > 2
292 TH2F * fhMassEtaNLocMaxN[8][2] ; //! Mass for selected around eta, N Local Maxima > 2
293 TH2F * fhMassConNLocMaxN[8][2] ; //! Mass for selected around close to 0, N Local Maxima > 2
bb2d339b 294
74858845 295 TH2F * fhMassAfterCutsNLocMax1[8][2] ; //! Mass after M02, asymmetry cuts for pi0, N Local Maxima = 1
296 TH2F * fhMassAfterCutsNLocMax2[8][2] ; //! Mass after M02, asymmetry cuts for pi0, N Local Maxima = 2
297 TH2F * fhMassAfterCutsNLocMaxN[8][2] ; //! Mass after M02, asymmetry cuts for pi0, N Local Maxima > 2
167f2534 298
c8710850 299 TH2F * fhAsyPi0NLocMax1[8][2] ; //! Asy for Mass around pi0, N Local Maxima = 1
300 TH2F * fhAsyEtaNLocMax1[8][2] ; //! Asy for Mass around eta, N Local Maxima = 1
301 TH2F * fhAsyConNLocMax1[8][2] ; //! Asy for Mass around close to 0, N Local Maxima = 1
bb2d339b 302
c8710850 303 TH2F * fhAsyPi0NLocMax2[8][2] ; //! Asy for Mass around pi0, N Local Maxima = 2
304 TH2F * fhAsyEtaNLocMax2[8][2] ; //! Asy for Mass around eta, N Local Maxima = 2
305 TH2F * fhAsyConNLocMax2[8][2] ; //! Asy for Mass around close to 0, N Local Maxima = 2
fc01318e 306
c8710850 307 TH2F * fhAsyPi0NLocMaxN[8][2] ; //! Asy for Mass around pi0, N Local Maxima > 2
308 TH2F * fhAsyEtaNLocMaxN[8][2] ; //! Asy for Mass around eta, N Local Maxima > 2
309 TH2F * fhAsyConNLocMaxN[8][2] ; //! Asy for Mass around close to 0, N Local Maxima > 2
e671adc2 310
de7d73e6 311 TH2F * fhSplitEFractionNLocMax1[8][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1
312 TH2F * fhSplitEFractionNLocMax2[8][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2
313 TH2F * fhSplitEFractionNLocMaxN[8][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2
667432ef 314
74858845 315 TH2F * fhSplitEFractionAfterCutsNLocMax1[8][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1, after M02 and asymmetry cut
316 TH2F * fhSplitEFractionAfterCutsNLocMax2[8][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2, after M02 and asymmetry cut
317 TH2F * fhSplitEFractionAfterCutsNLocMaxN[8][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2, after M02 and asymmetry cut
667432ef 318
de7d73e6 319 TH2F * fhMassSplitEFractionNLocMax1Ebin[8][4] ; //! Mass vs sum of splitted cluster energy / cluster energy for N max in cluster = 1, 1-6 for different MC particle types, not track matched
320 TH2F * fhMassSplitEFractionNLocMax2Ebin[8][4] ; //! Mass vs sum of splitted cluster energy / cluster energy for N max in cluster = 2, 1-6 for different MC particle types, not track matched
321 TH2F * fhMassSplitEFractionNLocMaxNEbin[8][4] ; //! Mass vs sum of splitted cluster energy / cluster energy for N max in cluster > 2, 1-6 for different MC particle types, not track matched
5094c724 322
c8710850 323 TH2F * fhAnglePairNLocMax1[2] ; //! pair opening angle vs E
324 TH2F * fhAnglePairNLocMax2[2] ; //! pair opening angle vs E
325 TH2F * fhAnglePairNLocMaxN[2] ; //! pair opening angle vs E
326
327 TH2F * fhAnglePairMassNLocMax1[2] ; //! pair opening angle vs Mass for E > 7 GeV
328 TH2F * fhAnglePairMassNLocMax2[2] ; //! pair opening angle vs Mass for E > 7 GeV
329 TH2F * fhAnglePairMassNLocMaxN[2] ; //! pair opening angle vs Mass for E > 7 GeV
330
331 TH2F * fhTrackMatchedDEtaNLocMax1[8] ; //! Eta distance between track and cluster vs cluster E, 1 local maximum
332 TH2F * fhTrackMatchedDPhiNLocMax1[8] ; //! Phi distance between track and cluster vs cluster E, 1 local maximum
333 TH2F * fhTrackMatchedDEtaNLocMax2[8] ; //! Eta distance between track and cluster vs cluster E, 2 local maximum
334 TH2F * fhTrackMatchedDPhiNLocMax2[8] ; //! Phi distance between track and cluster vs cluster E, 2 local maximum
335 TH2F * fhTrackMatchedDEtaNLocMaxN[8] ; //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
336 TH2F * fhTrackMatchedDPhiNLocMaxN[8] ; //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
b2e375c7 337
338 TH2F * fhTrackMatchedDEtaNLocMax1Pos[8] ; //! Eta distance between track and cluster vs cluster E, 1 local maximum
339 TH2F * fhTrackMatchedDPhiNLocMax1Pos[8] ; //! Phi distance between track and cluster vs cluster E, 1 local maximum
340 TH2F * fhTrackMatchedDEtaNLocMax2Pos[8] ; //! Eta distance between track and cluster vs cluster E, 2 local maximum
341 TH2F * fhTrackMatchedDPhiNLocMax2Pos[8] ; //! Phi distance between track and cluster vs cluster E, 2 local maximum
342 TH2F * fhTrackMatchedDEtaNLocMaxNPos[8] ; //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
343 TH2F * fhTrackMatchedDPhiNLocMaxNPos[8] ; //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
344
345 TH2F * fhTrackMatchedDEtaNLocMax1Neg[8] ; //! Eta distance between track and cluster vs cluster E, 1 local maximum
346 TH2F * fhTrackMatchedDPhiNLocMax1Neg[8] ; //! Phi distance between track and cluster vs cluster E, 1 local maximum
347 TH2F * fhTrackMatchedDEtaNLocMax2Neg[8] ; //! Eta distance between track and cluster vs cluster E, 2 local maximum
348 TH2F * fhTrackMatchedDPhiNLocMax2Neg[8] ; //! Phi distance between track and cluster vs cluster E, 2 local maximum
349 TH2F * fhTrackMatchedDEtaNLocMaxNNeg[8] ; //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
350 TH2F * fhTrackMatchedDPhiNLocMaxNNeg[8] ; //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
351
17f5b4b6 352 TH2F * fhCentralityPi0NLocMax1[8][2] ; //! Centrality for selected pi0, N Local Maxima = 1
353 TH2F * fhCentralityEtaNLocMax1[8][2] ; //! Centrality for selected eta, N Local Maxima = 1
354 TH2F * fhCentralityPi0NLocMax2[8][2] ; //! Centrality for selected pi0, N Local Maxima = 2
355 TH2F * fhCentralityEtaNLocMax2[8][2] ; //! Centrality for selected eta, N Local Maxima = 2
356 TH2F * fhCentralityPi0NLocMaxN[8][2] ; //! Centrality for selected pi0, N Local Maxima > 2
357 TH2F * fhCentralityEtaNLocMaxN[8][2] ; //! Centrality for selected eta, N Local Maxima > 2
c8710850 358
359 TH2F * fhEventPlanePi0NLocMax1 ; //! Event plane for selected pi0, N Local Maxima = 1
360 TH2F * fhEventPlaneEtaNLocMax1 ; //! Event plane for selected eta, N Local Maxima = 1
361 TH2F * fhEventPlanePi0NLocMax2 ; //! Event plane for selected pi0, N Local Maxima = 2
362 TH2F * fhEventPlaneEtaNLocMax2 ; //! Event plane for selected eta, N Local Maxima = 2
363 TH2F * fhEventPlanePi0NLocMaxN ; //! Event plane for selected pi0, N Local Maxima > 2
364 TH2F * fhEventPlaneEtaNLocMaxN ; //! Event plane for selected eta, N Local Maxima > 2
365
9554fc65 366 TH2F * fhClusterEtaPhiNLocMax1 ; //! Eta vs Phi of clusters with N Local Maxima = 1, E > 8 GeV
367 TH2F * fhClusterEtaPhiNLocMax2 ; //! Eta vs Phi of clusters with N Local Maxima = 2, E > 8 GeV
368 TH2F * fhClusterEtaPhiNLocMaxN ; //! Eta vs Phi of clusters with N Local Maxima > 2, E > 8 GeV
369 TH2F * fhPi0EtaPhiNLocMax1 ; //! Eta vs Phi of pi0's with N Local Maxima = 1, E > 8 GeV
370 TH2F * fhPi0EtaPhiNLocMax2 ; //! Eta vs Phi of pi0's with N Local Maxima = 2, E > 8 GeV
371 TH2F * fhPi0EtaPhiNLocMaxN ; //! Eta vs Phi of pi0's with N Local Maxima > N, E > 8 GeV
372 TH2F * fhEtaEtaPhiNLocMax1 ; //! Eta vs Phi of eta's with N Local Maxima = 1, E > 8 GeV
373 TH2F * fhEtaEtaPhiNLocMax2 ; //! Eta vs Phi of eta's with N Local Maxima = 2, E > 8 GeV
374 TH2F * fhEtaEtaPhiNLocMaxN ; //! Eta vs Phi of eta's with N Local Maxima > N, E > 8 GeV
375
dbe09c26 376 TH2F * fhPi0CellE[3] ; //! pi0's energy vs cluster cell energy with NLM = 1, = 2, > 2
377 TH2F * fhPi0CellEFrac[3] ; //! pi0's energy vs cluster cell energy fraction with NLM = 1, = 2, > 2
378 TH2F * fhPi0CellLogEFrac[3] ; //! pi0's energy vs cluster log cell energy fraction with NLM = 1, = 2, > 2
19391b8c 379 TH2F * fhPi0CellEMaxEMax2Frac [3]; //! pi0's energy vs fraction of 2 main maxima energy with NLM = 1, = 2, > 2
380 TH2F * fhPi0CellEMaxClusterFrac [3]; //! pi0's energy vs energy fraction of main LM and cluster energy with NLM = 1, = 2, > 2
381 TH2F * fhPi0CellEMax2ClusterFrac[3]; //! pi0's energy vs energy fraction of second LM and cluster energy with NLM = 1, = 2, > 2
382 TH2F * fhPi0CellEMaxFrac [3]; //! pi0's energy vs energy fraction of main LM and cluster cell energy with NLM = 1, = 2, > 2
383 TH2F * fhPi0CellEMax2Frac [3]; //! pi0's energy vs energy fraction of second LM and cluster cell energy with NLM = 1, = 2, > 2
384
dbe09c26 385 TH2F * fhM02WeightPi0[3][10] ; //! M02 for selected pi0 with different weight, with NLM = 1, = 2, > 2
19391b8c 386 TH2F * fhM02ECellCutPi0[3][10] ; //! M02 for selected pi0 with different cut on cell energy, with NLM = 1, = 2, > 2
dbe09c26 387
a1fd1b69 388 TH2F * fhPi0EPairDiffTimeNLM1; //! E vs Pair of clusters time difference vs E, for selected pi0, NLM=1
389 TH2F * fhPi0EPairDiffTimeNLM2; //! E vs Pair of clusters time difference vs E, for selected pi0, NLM=2
390 TH2F * fhPi0EPairDiffTimeNLMN; //! E vs Pair of clusters time difference vs E, for selected pi0, NLM>2
391 TH2F * fhEtaEPairDiffTimeNLM1; //! E vs Pair of clusters time difference vs E, for selected eta, NLM=1
392 TH2F * fhEtaEPairDiffTimeNLM2; //! E vs Pair of clusters time difference vs E, for selected eta, NLM=2
393 TH2F * fhEtaEPairDiffTimeNLMN; //! E vs Pair of clusters time difference vs E, for selected eta, NLM>2
394
b2e375c7 395 TH2F * fhMCEM02Overlap0[3][8]; //! E vs M02 for different MC origin, no other MC particles contributes, neutral cluster
396 TH2F * fhMCEM02Overlap1[3][8]; //! E vs M02 for different MC origin, 1 other MC particles contributes, neutral cluster
397 TH2F * fhMCEM02OverlapN[3][8]; //! E vs M02 for different MC origin, N other MC particles contributes, neutral cluster
398 TH2F * fhMCEM02Overlap0Match[3][8]; //! E vs M02 for different MC origin, no other MC particles contributes, charged cluster
399 TH2F * fhMCEM02Overlap1Match[3][8]; //! E vs M02 for different MC origin, 1 other MC particles contributes, charged cluster
400 TH2F * fhMCEM02OverlapNMatch[3][8]; //! E vs M02 for different MC origin, N other MC particles contributes, charged cluster
401
402 TH2F * fhMCEMassOverlap0[3][8]; //! E vs Mass for different MC origin, no other MC particles contributes, neutral cluster
403 TH2F * fhMCEMassOverlap1[3][8]; //! E vs Mass for different MC origin, 1 other MC particles contributes, neutral cluster
404 TH2F * fhMCEMassOverlapN[3][8]; //! E vs Mass for different MC origin, N other MC particles contributes, neutral cluster
405 TH2F * fhMCEMassOverlap0Match[3][8]; //! E vs Mass for different MC origin, no other MC particles contributes, charged cluster
406 TH2F * fhMCEMassOverlap1Match[3][8]; //! E vs Mass for different MC origin, 1 other MC particles contributes, charged cluster
407 TH2F * fhMCEMassOverlapNMatch[3][8]; //! E vs Mass for different MC origin, N other MC particles contributes, charged cluster
408
83351853 409 TH2F * fhMCESplitEFracOverlap0[3][8]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, no other MC particles contributes, neutral cluster
410 TH2F * fhMCESplitEFracOverlap1[3][8]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, 1 other MC particles contributes, neutral cluster
411 TH2F * fhMCESplitEFracOverlapN[3][8]; //! E vs sum of splitted cluster energy / cluster energyfor different MC origin, N other MC particles contributes, neutral cluster
412 TH2F * fhMCESplitEFracOverlap0Match[3][8]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, no other MC particles contributes, charged cluster
413 TH2F * fhMCESplitEFracOverlap1Match[3][8]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, 1 other MC particles contributes, charged cluster
414 TH2F * fhMCESplitEFracOverlapNMatch[3][8]; //! E vs sum of splitted cluster energy / cluster energyc for different MC origin, N other MC particles contributes, charged cluster
415
b2e375c7 416 TH2F * fhMCPi0MassM02Overlap0[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, no other MC particles contributes, neutral cluster, 4 E bins
417 TH2F * fhMCPi0MassM02Overlap1[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, 1 other MC particles contributes, neutral cluster, 4 E bins
418 TH2F * fhMCPi0MassM02OverlapN[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, N other MC particles contributes, neutral cluster, 4 E bins
419 TH2F * fhMCPi0MassM02Overlap0Match[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, no other MC particles contributes, charged cluster, 4 E bins
420 TH2F * fhMCPi0MassM02Overlap1Match[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, 1 other MC particles contributes, charged cluster, 4 E bins
421 TH2F * fhMCPi0MassM02OverlapNMatch[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, N other MC particles contributes, charged cluster, 4 E bins
422
423 TH2F * fhMCENOverlaps[3][8]; //! E vs number of Overlaps in MC, neutral cluster
424 TH2F * fhMCENOverlapsMatch[3][8]; //! E vs number of Overlaps in MC, charged cluster
425
426 TH2F * fhMCPi0HighNLMPair; //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons
427 TH2F * fhMCPi0LowNLMPair; //! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy correspond to 2 photons
428 TH2F * fhMCPi0AnyNLMPair; //! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs correspond to 2 photons
429 TH2F * fhMCPi0NoneNLMPair; //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons
b2e375c7 430 // No match between highest energy local maxima and highest energy MC particle
431 TH2F * fhMCPi0HighNLMPairNoMCMatch; //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons
432 TH2F * fhMCPi0LowNLMPairNoMCMatch; //! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy correspond to 2 photons
433 TH2F * fhMCPi0AnyNLMPairNoMCMatch; //! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs correspond to 2 photons
434 TH2F * fhMCPi0NoneNLMPairNoMCMatch; //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons
36769d30 435
436 TH2F * fhMCPi0HighNLMPairOverlap; //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons, overlap
437 TH2F * fhMCPi0LowNLMPairOverlap; //! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy correspond to 2 photons, overlap
438 TH2F * fhMCPi0AnyNLMPairOverlap; //! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs correspond to 2 photons, overlap
439 TH2F * fhMCPi0NoneNLMPairOverlap; //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons, overlap
440 // No match between highest energy local maxima and highest energy MC particle
441 TH2F * fhMCPi0HighNLMPairNoMCMatchOverlap; //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons, overlap
442 TH2F * fhMCPi0LowNLMPairNoMCMatchOverlap; //! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy correspond to 2 photons, overlap
443 TH2F * fhMCPi0AnyNLMPairNoMCMatchOverlap; //! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs correspond to 2 photons, overlap
444 TH2F * fhMCPi0NoneNLMPairNoMCMatchOverlap; //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons, overlap
445
83351853 446 TH2F * fhMCPi0DecayPhotonHitHighLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
447 TH2F * fhMCPi0DecayPhotonAdjHighLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
448 TH2F * fhMCPi0DecayPhotonHitOtherLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
449 TH2F * fhMCPi0DecayPhotonAdjOtherLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high
36769d30 450 TH2F * fhMCPi0DecayPhotonAdjacent; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM
83351853 451 TH2F * fhMCPi0DecayPhotonHitNoLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas
b2e375c7 452
36769d30 453 TH2F * fhMCPi0DecayPhotonHitHighLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maxima, overlap
454 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
455 TH2F * fhMCPi0DecayPhotonHitOtherLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
456 TH2F * fhMCPi0DecayPhotonAdjOtherLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high, overlap
457 TH2F * fhMCPi0DecayPhotonAdjacentOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM, overlap
458 TH2F * fhMCPi0DecayPhotonHitNoLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas, overlap
459
b2e375c7 460 TH2F * fhMCEOverlapType; //! what particles overlap with pi0, neutral clusters
461 TH2F * fhMCEOverlapTypeMatch; //! what particles overlap with pi0, charged clusters
dbe09c26 462
2a77f6f4 463 AliAnaInsideClusterInvariantMass( const AliAnaInsideClusterInvariantMass & split) ; // cpy ctor
464 AliAnaInsideClusterInvariantMass & operator = (const AliAnaInsideClusterInvariantMass & split) ; // cpy assignment
992b14a7 465
b2e375c7 466 ClassDef(AliAnaInsideClusterInvariantMass,21)
992b14a7 467
468} ;
469
470#endif //ALIANAINSIDECLUSTERINVARIANTMASS_H
471
472
473