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 */
6 //_________________________________________________________________________
8 // Split clusters with some criteria and calculate invariant mass
9 // to identify them as pi0 or conversion
12 //-- Author: Gustavo Conesa (LPSC-Grenoble)
13 //_________________________________________________________________________
16 // --- ROOT system ---
21 // --- ANALYSIS system ---
22 class AliAODCaloCluster;
24 #include "AliAnaCaloTrackCorrBaseClass.h"
26 class AliAnaInsideClusterInvariantMass : public AliAnaCaloTrackCorrBaseClass {
30 AliAnaInsideClusterInvariantMass() ; // default ctor
31 virtual ~AliAnaInsideClusterInvariantMass() { ; } //virtual dtor
33 void CheckLocalMaximaMCOrigin(AliVCluster* cluster, Int_t mcindex, Int_t noverlaps,
34 Float_t e1, Float_t e2, Float_t mass);
35 //, Float_t m02, TLorentzVector l1, TLorentzVector l2);
37 TObjString * GetAnalysisCuts();
39 TList * GetCreateOutputObjects();
41 void GetMCIndex(AliVCluster * cluster, Int_t & mcindex, Int_t & mcTag);
43 void GetMCPrimaryKine(AliVCluster* cluster, Int_t mcindex, Int_t mcTag, Bool_t matched,
44 Float_t & eprim, Float_t & asymGen, Float_t & angleGen, Int_t & noverlaps );
46 void FillAngleHistograms(Int_t nMax, Bool_t matched, Int_t mcindex,
47 Float_t en, Float_t e1 , Float_t e2,
48 Float_t angle, Float_t mass, Float_t anglePrim,
49 Float_t m02, Float_t asym, Int_t pid, Int_t noverlaps);
51 void FillArmenterosHistograms(Int_t nMax, Int_t ebin, Int_t mcindex,
52 Float_t pi0E, TLorentzVector g1, TLorentzVector g2,
53 Float_t m02, Int_t pid);
55 void FillThetaStarHistograms(Int_t nMax, Bool_t matched, Int_t mcindex,
56 Float_t pi0E, TLorentzVector g1, TLorentzVector g2,
57 Float_t m02, Int_t pid);
59 void FillEBinHistograms(Int_t ebin, Int_t nMax, Int_t mcindex, Float_t splitFrac,
60 Float_t mass, Float_t asym, Float_t l0);
62 void FillMCHistograms(Float_t en, Float_t e1 , Float_t e2,
63 Int_t ebin, Int_t mcindex,Int_t noverlaps,
64 Float_t l0, Float_t mass,
65 Int_t nMax, Bool_t matched,
66 Float_t splitFrac, Float_t asym,
67 Float_t eprim, Float_t asymGen);
69 void FillMCOverlapHistograms(Float_t en, Float_t enprim,
70 Int_t nc, Float_t mass, Float_t l0,
71 Float_t asym, Float_t splitFrac,
72 Int_t nlm, Int_t ebin, Bool_t matched,
73 Int_t mcindex, Int_t noverlaps);
75 void FillSSWeightHistograms(AliVCluster *cluster, Int_t nlm, Int_t absId1, Int_t absId2);
77 void FillSSExtraHistograms(AliVCluster *cluster, Int_t nMax,
78 Bool_t matched, Int_t mcindex,
79 Float_t mass , Int_t ebin);
81 void FillNCellHistograms(Int_t ncells, Float_t energy, Int_t nMax,
82 Bool_t matched, Int_t mcindex,
83 Float_t mass , Float_t l0);
85 void FillTrackMatchingHistograms(AliVCluster * cluster,Int_t nMax, Int_t mcindex);
87 void FillHistograms1(Float_t en, Float_t e1, Float_t e2,
88 Int_t nMax, Float_t mass, Float_t l0,
89 Float_t eta, Float_t phi,
90 Bool_t matched, Int_t mcindex);
93 void FillHistograms2(Float_t en, Float_t eprim,
94 Float_t e1, Float_t e2, Int_t nMax,
95 Float_t mass, Float_t l0,
96 Bool_t matched, Int_t mcindex);
98 void FillIdPi0Histograms(Float_t en, Float_t e1, Float_t e2,
99 Int_t nc, Int_t nMax, Float_t t12diff,
100 Float_t mass, Float_t l0,
101 Float_t eta, Float_t phi,
102 Bool_t matched, Int_t mcindex);
104 void FillIdEtaHistograms(Float_t en, Float_t e1, Float_t e2,
105 Int_t nc, Int_t nMax, Float_t t12diff,
106 Float_t mass, Float_t l0,
107 Float_t eta, Float_t phi,
108 Bool_t matched, Int_t mcindex);
110 void FillIdConvHistograms(Float_t en, Int_t nMax, Float_t asym,
111 Float_t mass, Float_t l0,
112 Bool_t matched, Int_t mcindex);
116 void InitParameters();
118 void MakeAnalysisFillHistograms() ;
120 void Print(const Option_t * opt) const;
122 void SetCalorimeter(TString & det) { fCalorimeter = det ; }
124 void SetMinNCells(Int_t cut) { fMinNCells = cut ; }
126 void SetMinBadChannelDistance(Float_t cut) { fMinBadDist = cut ; }
128 void SetWCorrectionParameter(Float_t p = 0.07) { fWSimu = p ; }
130 void SwitchOnFillAngleHistograms() { fFillAngleHisto = kTRUE ; }
131 void SwitchOffFillAngleHistograms() { fFillAngleHisto = kFALSE ; }
133 void SwitchOnFillArmenterosHistograms() { fFillArmenterosHisto = kTRUE ; }
134 void SwitchOffFillArmenterosHistograms() { fFillArmenterosHisto = kFALSE ; }
136 void SwitchOnFillThetaStarHistograms() { fFillThetaStarHisto = kTRUE ; }
137 void SwitchOffFillThetaStarHistograms() { fFillThetaStarHisto = kFALSE ; }
139 void SwitchOnFillExtraSSHistograms() { fFillSSExtraHisto = kTRUE ; }
140 void SwitchOffFillExtraSSHistograms() { fFillSSExtraHisto = kFALSE ; }
142 void SwitchOnFillHighMultHistograms() { fFillHighMultHisto = kTRUE ; }
143 void SwitchOffFillHighMultHistograms() { fFillHighMultHisto = kFALSE ; }
145 void SwitchOnFillIdConvHistograms() { fFillIdConvHisto = kTRUE ; }
146 void SwitchOffFillIdConvHistograms() { fFillIdConvHisto = kFALSE ; }
148 void SwitchOnFillIdEtaHistograms() { fFillIdEtaHisto = kTRUE ; }
149 void SwitchOffFillIdEtaHistograms() { fFillIdEtaHisto = kFALSE ; }
151 void SwitchOnFillTMHistograms() { fFillTMHisto = kTRUE ; }
152 void SwitchOffFillTMHistograms() { fFillTMHisto = kFALSE ; }
154 void SwitchOnFillTMResidualHistograms() { fFillTMResidualHisto = kTRUE ; }
155 void SwitchOffFillTMResidualHistograms() { fFillTMResidualHisto = kFALSE ; }
157 void SwitchOnFillMCPrimaryHistograms() { fFillMCHisto = kTRUE ; }
158 void SwitchOffFillMCPrimaryHistograms() { fFillMCHisto = kFALSE ; }
160 void SwitchOnFillSSWeightHistograms() { fFillSSWeightHisto = kTRUE ; }
161 void SwitchOffFillSSWeightHistograms() { fFillSSWeightHisto = kFALSE ; }
163 void SwitchOnFillEbinHistograms() { fFillEbinHisto = kTRUE ; }
164 void SwitchOffFillEbinHistograms() { fFillEbinHisto = kFALSE ; }
166 void SwitchOnFillMCOverlapHistograms() { fFillMCOverlapHisto = kTRUE ; }
167 void SwitchOffFillMCOverlapHistograms() { fFillMCOverlapHisto = kFALSE ; }
169 void SwitchOnFillNCellHistograms() { fFillNCellHisto = kTRUE ; }
170 void SwitchOffFillNCellHistograms() { fFillNCellHisto = kFALSE ; }
172 void SwitchOnSplitClusterDistToBad() { fCheckSplitDistToBad = kTRUE ; }
173 void SwitchOffSplitClusterDistToBad() { fCheckSplitDistToBad = kFALSE ; }
175 void SetNWeightForShowerShape(Int_t n) { fSSWeightN = n ; }
176 void SetWeightForShowerShape(Int_t i, Float_t v) { if (i < 10) fSSWeight[i] = v ; }
178 void SetNECellCutForShowerShape(Int_t n) { fSSECellCutN = n ; }
179 void SetECellCutForShowerShape(Int_t i, Float_t v) { if (i < 10) fSSECellCut[i] = v ; }
182 void RecalculateClusterShowerShapeParametersWithCellCut(const AliEMCALGeometry * geom, AliVCaloCells* cells, AliVCluster * cluster,
183 Float_t & l0, Float_t & l1,
184 Float_t & disp, Float_t & dEta, Float_t & dPhi,
185 Float_t & sEta, Float_t & sPhi, Float_t & sEtaPhi,Float_t eCellMin = 0.);
189 enum mcTypes { kmcPhoton = 1, kmcConversion = 2, kmcPi0 = 3, kmcPi0Conv = 4,
190 kmcEta = 5, kmcHadron = 6 };
194 TString fCalorimeter ; // Calorimeter where the gamma is searched
195 Int_t fMinNCells ; // Study clusters with ncells larger than cut
196 Float_t fMinBadDist ; // Minimal distance to bad channel to accept cluster
197 Float_t fHistoECut ; // Fixed E cut for some histograms
198 Bool_t fCheckSplitDistToBad; // Check the distance to bad channel and to EMCal borders of split clusters
200 Bool_t fFillAngleHisto; // Fill splitted clusters angle histograms
201 Bool_t fFillTMHisto ; // Fill track matching histos,
202 Bool_t fFillTMResidualHisto ; // Fill track matching histos, residuals
203 Bool_t fFillSSExtraHisto ; // Fill shower shape extra histos
204 Bool_t fFillMCHisto ; // Fill MC energy fraction histos
205 Bool_t fFillSSWeightHisto ; // Fill weigth histograms
206 Bool_t fFillEbinHisto ; // Fill E bin histograms
207 Bool_t fFillMCOverlapHisto ; // Fill MC particles overlap histograms
208 Bool_t fFillNCellHisto ; // Fill n cells in cluster dependent histograms
209 Bool_t fFillIdConvHisto ; // Fill histograms for clusters identified as conversion
210 Bool_t fFillIdEtaHisto ; // Fill histograms for clusters identified as Eta
211 Bool_t fFillHighMultHisto; // Fill centrality/event plane histograms
212 Bool_t fFillArmenterosHisto; // Fill armenteros type histo
213 Bool_t fFillThetaStarHisto; // Fill cosThetaStar histos
215 Float_t fSSWeight[10]; // List of weights to test
216 Int_t fSSWeightN; // Total number of weights to test
218 Float_t fSSECellCut[10]; // List of cell min energy cuts to test
219 Int_t fSSECellCutN; // Total number of cell min energy cuts to test
221 Float_t fWSimu; // Slope of the linear correction factor for the shower
222 // shape weight in simulation, about 0.07
226 TH2F * fhMassNLocMax1[7][2] ; //! Split Inv Mass vs cluster E, NLM=1, different MC particle types, track matching on/off
227 TH2F * fhMassNLocMax2[7][2] ; //! Split Inv Mass vs cluster E, NLM=2, different MC particle types, track matching on/off
228 TH2F * fhMassNLocMaxN[7][2] ; //! Split Inv Mass vs cluster E, NLM>2, different MC particle types, track matching on/off
230 TH2F * fhMassSplitENLocMax1[7][2] ; //! Split Inv Mass vs E1+E2, NLM=1, different MC particle types, track matching on/off
231 TH2F * fhMassSplitENLocMax2[7][2] ; //! Split Inv Mass vs E1+E2, NLM=2, different MC particle types, track matching on/off
232 TH2F * fhMassSplitENLocMaxN[7][2] ; //! Split Inv Mass vs E1+E2, NLM>2, different MC particle types, track matching on/off
234 TH2F * fhAsymNLocMax1[7][2] ; //! Asymmetry of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types
235 TH2F * fhAsymNLocMax2[7][2] ; //! Asymmetry of 2 cells local maxima vs E, 1-6 for different MC particle types
236 TH2F * fhAsymNLocMaxN[7][2] ; //! Asymmetry of >2 cells local maxima vs E, 1-6 for different MC particle types
238 TH2F * fhSplitEFractionvsAsyNLocMax1[2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1 vs |A|
239 TH2F * fhSplitEFractionvsAsyNLocMax2[2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2 vs |A|
240 TH2F * fhSplitEFractionvsAsyNLocMaxN[2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2 vs |A|
242 TH2F * fhMassAsyCutNLocMax1 ; //! Mass(E) asym selection, not matched, Mass of split clusters, NLM = 1
243 TH2F * fhMassAsyCutNLocMax2 ; //! Mass(E) asym selection, not matched, Mass of split clusters, NLM = 1
244 TH2F * fhMassAsyCutNLocMaxN ; //! Mass(E) asym selection, not matched, Mass of split clusters, NLM > 2
246 TH2F * fhM02AsyCutNLocMax1 ; //! M02(E) asym selection, not matched, M02, NLM = 1
247 TH2F * fhM02AsyCutNLocMax2 ; //! M02(E) asym selection, not matched, M02, NLM = 2
248 TH2F * fhM02AsyCutNLocMaxN ; //! M02(E) asym selection, not matched, M02, NLM > 2
250 TH2F * fhMassM02CutNLocMax1 ; //! Mass(E) M02 selection, not matched, Mass of split clusters, NLM = 1
251 TH2F * fhMassM02CutNLocMax2 ; //! Mass(E) M02 selection, not matched, Mass of split clusters, NLM = 1
252 TH2F * fhMassM02CutNLocMaxN ; //! Mass(E) M02 selection, not matched, Mass of split clusters, NLM > 2
254 TH2F * fhAsymM02CutNLocMax1 ; //! Asym(E) M02 selection, not matched, energy asymmetry of split clusters, NLM = 1
255 TH2F * fhAsymM02CutNLocMax2 ; //! Asym(E) M02 selection, not matched, energy asymmetry of split clusters, NLM = 2
256 TH2F * fhAsymM02CutNLocMaxN ; //! Asym(E) M02 selection, not matched, energy asymmetry of split clusters, NLM > 2
258 TH2F * fhMassEnCutNLocMax1 ; //! Mass(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
259 TH2F * fhMassEnCutNLocMax2 ; //! Mass(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
260 TH2F * fhMassEnCutNLocMaxN ; //! Mass(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
262 TH2F * fhM02EnCutNLocMax1 ; //! M02(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
263 TH2F * fhM02EnCutNLocMax2 ; //! M02(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
264 TH2F * fhM02EnCutNLocMaxN ; //! M02(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
266 TH2F * fhAsymEnCutNLocMax1 ; //! Asym(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
267 TH2F * fhAsymEnCutNLocMax2 ; //! Asym(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
268 TH2F * fhAsymEnCutNLocMaxN ; //! Asym(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
270 TH2F * fhSplitEFracEnCutNLocMax1 ; //! Split E fraction (E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
271 TH2F * fhSplitEFracEnCutNLocMax2 ; //! Split E fraction (E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
272 TH2F * fhSplitEFracEnCutNLocMaxN ; //! Split E fraction (E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
274 TH2F * fhMassSplitECutNLocMax1 ; //! 85% of split energy, not matched, Mass of split clusters, NLM = 1
275 TH2F * fhMassSplitECutNLocMax2 ; //! 85% of split energy, not matched, Mass of split clusters, NLM = 1
276 TH2F * fhMassSplitECutNLocMaxN ; //! 85% of split energy, not matched, Mass of split clusters, NLM > 2
278 TH2F * fhMassM02NLocMax1[7][2] ; //! Mass of splitted clusters when 1 local max vs M02, for E > 8 GeV, 1-6 for different MC particle types
279 TH2F * fhMassM02NLocMax2[7][2] ; //! Mass of splitted clusters when 2 local max vs M02, for E > 8 GeV, 1-6 for different MC particle types
280 TH2F * fhMassM02NLocMaxN[7][2] ; //! Mass of splitted clusters when >2 local max vs M02, for E > 8 GeV, 1-6 for different MC particle types
282 TH2F * fhMassM02NLocMax1Ebin[4] ; //! Mass of splitted clusters when 1 local max vs M02, 4 E bins, neutral clusters
283 TH2F * fhMassM02NLocMax2Ebin[4] ; //! Mass of splitted clusters when 2 local max vs M02, 4 E bins, neutral clusters
284 TH2F * fhMassM02NLocMaxNEbin[4] ; //! Mass of splitted clusters when >2 local max vs M02, 4 E bins, neutral clusters
286 TH2F * fhMassAsyNLocMax1Ebin[4] ; //! Mass of Mass of splitted clusters when 1 local max vs asymmetry, 4 E bins, neutral clusters
287 TH2F * fhMassAsyNLocMax2Ebin[4] ; //! Mass of Mass of splitted clusters when 2 local max vs asymmetry, 4 E bins, neutral clusters
288 TH2F * fhMassAsyNLocMaxNEbin[4] ; //! Mass of Mass of splitted clusters when >2 local max vs asymmetry, 4 E bins, neutral clusters
290 TH2F * fhAsyMCGenRecoNLocMax1EbinPi0[4] ; //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when 1 local max, 4 E bins, neutral clusters
291 TH2F * fhAsyMCGenRecoNLocMax2EbinPi0[4] ; //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when 2 local max, 4 E bins, neutral clusters
292 TH2F * fhAsyMCGenRecoNLocMaxNEbinPi0[4] ; //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when >2 local max, 4 E bins, neutral clusters
294 TH2F * fhAsyMCGenRecoDiffMCPi0[3]; //! reconstructed-generated asymmetry of splitted clusters vs E from pi0, for 3 NLM cases
295 TH2F * fhAsyMCGenRecoDiffMCPi0Conv[3]; //! reconstructed-generated asymmetry of splitted clusters vs E from converted pi0, for 3 NLM cases
297 TH2F * fhMassDispEtaNLocMax1[7][2] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, for E > 8 GeV, 1-6 for different MC particle types
298 TH2F * fhMassDispEtaNLocMax2[7][2] ; //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
299 TH2F * fhMassDispEtaNLocMaxN[7][2] ; //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
301 TH2F * fhMassDispEtaNLocMax1Ebin[4] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
302 TH2F * fhMassDispEtaNLocMax2Ebin[4] ; //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
303 TH2F * fhMassDispEtaNLocMaxNEbin[4] ; //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
305 TH2F * fhMassDispPhiNLocMax1[7][2] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, for E > 8 GeV, 1-6 for different MC particle types
306 TH2F * fhMassDispPhiNLocMax2[7][2] ; //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
307 TH2F * fhMassDispPhiNLocMaxN[7][2] ; //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
309 TH2F * fhMassDispPhiNLocMax1Ebin[4] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
310 TH2F * fhMassDispPhiNLocMax2Ebin[4] ; //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
311 TH2F * fhMassDispPhiNLocMaxNEbin[4] ; //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
313 TH2F * fhMassDispAsyNLocMax1[7][2] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, for E > 8 GeV, 1-6 for different MC particle types
314 TH2F * fhMassDispAsyNLocMax2[7][2] ; //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
315 TH2F * fhMassDispAsyNLocMaxN[7][2] ; //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
317 TH2F * fhMassDispAsyNLocMax1Ebin[4] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
318 TH2F * fhMassDispAsyNLocMax2Ebin[4] ; //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
319 TH2F * fhMassDispAsyNLocMaxNEbin[4] ; //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
321 TH2F * fhNLocMax [7][2] ; //! Number of maxima in cluster vs E, 1-6 for different MC particle types
322 TH2F * fhNLocMaxM02Cut[7][2] ; //! Number of maxima in cluster vs E, 1-6 for different MC particle types, after SS cut
323 TH2F * fhNLocMaxIdPi0 [7][2] ; //! Number of maxima in cluster vs E, 1-6 for different MC particle types, after pi0 selection
325 TH2F * fhSplitClusterENLocMax[7][2] ; //! Number of maxima in cluster vs E of splitted clusters, 1-6 for different MC particle types
326 TH2F * fhSplitClusterEPi0NLocMax[7][2] ; //! Number of maxima in cluster vs E of splitted clusters when cluster id as pi0, 1-6 for different MC particle types
328 TH2F * fhLM1NLocMax [7][2] ; //! Split cluster 1 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types
329 TH2F * fhLM1NLocMaxM02Cut[7][2] ; //! Split cluster 1 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types, after SS cut
330 TH2F * fhLM1NLocMaxIdPi0 [7][2] ; //! Split cluster 1 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types, pi0 selection
332 TH2F * fhLM2NLocMax [7][2] ; //! Split cluster 2 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types
333 TH2F * fhLM2NLocMaxM02Cut[7][2] ; //! Split cluster 2 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types, after SS cut
334 TH2F * fhLM2NLocMaxIdPi0 [7][2] ; //! Split cluster 2 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types, pi0 selection
336 TH2F * fhM02NLocMax1 [7][2] ; //! M02 vs E for N max in cluster = 1, 1-6 for different MC particle types
337 TH2F * fhM02NLocMax2 [7][2] ; //! M02 vs E for N max in cluster = 2, 1-6 for different MC particle types
338 TH2F * fhM02NLocMaxN [7][2] ; //! M02 vs E for N max in cluster > 2, 1-6 for different MC particle types
340 TH2F * fhMCAsymM02NLocMax1MCPi0Ebin[4] ; //! M02 vs decay asymmetry for N max in cluster = 1, for 4 energy bins
341 TH2F * fhMCAsymM02NLocMax2MCPi0Ebin[4] ; //! M02 vs decay asymmetry for N max in cluster = 2, for 4 energy bins
342 TH2F * fhMCAsymM02NLocMaxNMCPi0Ebin[4] ; //! M02 vs decay asymmetry for N max in cluster > 2, for 4 energy bins
344 TH2F * fhMCGenFracNLocMax1[7][2] ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types
345 TH2F * fhMCGenFracNLocMax2[7][2] ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types
346 TH2F * fhMCGenFracNLocMaxN[7][2] ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types
348 TH2F * fhMCGenFracNLocMax1NoOverlap[7][2] ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types, no overlap found
349 TH2F * fhMCGenFracNLocMax2NoOverlap[7][2] ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types, no overlap found
350 TH2F * fhMCGenFracNLocMaxNNoOverlap[7][2] ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types, no overlap found
352 TH2F * fhMCGenFracAfterCutsNLocMax1MCPi0 ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, MCPi0 after M02 and asymmetry cut
353 TH2F * fhMCGenFracAfterCutsNLocMax2MCPi0 ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, MCPi0, after M02 and asymmetry cut
354 TH2F * fhMCGenFracAfterCutsNLocMaxNMCPi0 ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, MCPi0, after M02 and asymmetry cut
356 TH2F * fhMCGenSplitEFracNLocMax1[7][2] ; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types
357 TH2F * fhMCGenSplitEFracNLocMax2[7][2] ; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types
358 TH2F * fhMCGenSplitEFracNLocMaxN[7][2] ; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types
360 TH2F * fhMCGenSplitEFracNLocMax1NoOverlap[7][2];//! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types, no overlap
361 TH2F * fhMCGenSplitEFracNLocMax2NoOverlap[7][2];//! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types, no overlap
362 TH2F * fhMCGenSplitEFracNLocMaxNNoOverlap[7][2];//! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types, no overlap
364 TH2F * fhMCGenSplitEFracAfterCutsNLocMax1MCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types
365 TH2F * fhMCGenSplitEFracAfterCutsNLocMax2MCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types
366 TH2F * fhMCGenSplitEFracAfterCutsNLocMaxNMCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types
368 TH2F * fhMCGenEFracvsSplitEFracNLocMax1[7][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster = 1, MC pi0
369 TH2F * fhMCGenEFracvsSplitEFracNLocMax2[7][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster = 2, MC pi0
370 TH2F * fhMCGenEFracvsSplitEFracNLocMaxN[7][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster > 2, MC pi0
372 TH2F * fhMCGenEvsSplitENLocMax1[7][2] ; //! E generated particle vs E1+E2 for N max in cluster = 1, 1-6 for different MC particle types
373 TH2F * fhMCGenEvsSplitENLocMax2[7][2] ; //! E generated particle vs E1+E2 for N max in cluster = 2, 1-6 for different MC particle types
374 TH2F * fhMCGenEvsSplitENLocMaxN[7][2] ; //! E generated particle vs E1+E2 for N max in cluster > 2, 1-6 for different MC particle types
376 TH2F * fhMCGenFracNLocMaxEbin[7][4] ; //! NLM vs E generated particle / E reconstructed vs E reconstructed 1-6 for different MC particle types, not matched to track
377 TH2F * fhMCGenFracNLocMaxEbinMatched[7][4] ; //! NLM vs E generated particle / E reconstructed vs E reconstructed 1-6 for different MC particle types, matched to track
379 TH2F * fhM02MCGenFracNLocMax1Ebin[7][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
380 TH2F * fhM02MCGenFracNLocMax2Ebin[7][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
381 TH2F * fhM02MCGenFracNLocMaxNEbin[7][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
383 TH2F * fhMassMCGenFracNLocMax1Ebin[7][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
384 TH2F * fhMassMCGenFracNLocMax2Ebin[7][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
385 TH2F * fhMassMCGenFracNLocMaxNEbin[7][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
387 TH2F * fhNCellNLocMax1[7][2] ; //! n cells in cluster vs E for N max in cluster = 1, 1-6 for different MC particle types
388 TH2F * fhNCellNLocMax2[7][2] ; //! n cells in cluster vs E for N max in cluster = 2, 1-6 for different MC particle types
389 TH2F * fhNCellNLocMaxN[7][2] ; //! n cells in cluster vs E for N max in cluster > 2, 1-6 for different MC particle types
391 TH2F * fhNCellMassEHighNLocMax1MCPi0 ; //! n cells in cluster vs mass for high energy clusters, for N max in cluster = 1, for MC pi0
392 TH2F * fhNCellM02EHighNLocMax1MCPi0 ; //! n cells in cluster vs m02 for high energy clusters, for N max in cluster = 1, for MC pi0
393 TH2F * fhNCellMassELowNLocMax1MCPi0 ; //! n cells in cluster vs mass for low energy clusters, for N max in cluster = 1, for MC pi0
394 TH2F * fhNCellM02ELowNLocMax1MCPi0 ; //! n cells in cluster vs m02 for low energy clusters, for N max in cluster = 1, for MC pi0
396 TH2F * fhNCellMassEHighNLocMax2MCPi0 ; //! n cells in cluster vs mass for high energy clusters, for N max in cluster = 2, for MC pi0
397 TH2F * fhNCellM02EHighNLocMax2MCPi0 ; //! n cells in cluster vs m02 for high energy clusters, for N max in cluster = 2, for MC pi0
398 TH2F * fhNCellMassELowNLocMax2MCPi0 ; //! n cells in cluster vs mass for low energy clusters, for N max in cluster = 2, for MC pi0
399 TH2F * fhNCellM02ELowNLocMax2MCPi0 ; //! n cells in cluster vs m02 for low energy clusters, for N max in cluster = 2, for MC pi0
401 TH2F * fhNCellMassEHighNLocMaxNMCPi0 ; //! n cells in cluster vs mass for high energy clusters, for N max in cluster > 2, for MC pi0
402 TH2F * fhNCellM02EHighNLocMaxNMCPi0 ; //! n cells in cluster vs m02 for high energy clusters, for N max in cluster > 2, for MC pi0
403 TH2F * fhNCellMassELowNLocMaxNMCPi0 ; //! n cells in cluster vs mass for low energy clusters, for N max in cluster > 2, for MC pi0
404 TH2F * fhNCellM02ELowNLocMaxNMCPi0 ; //! n cells in cluster vs m02 for low energy clusters, for N max in cluster > 2, for MC pi0
406 TH2F * fhM02Pi0NLocMax1[7][2] ; //! M02 for Mass around pi0, N Local Maxima = 1
407 TH2F * fhM02EtaNLocMax1[7][2] ; //! M02 for Mass around eta, N Local Maxima = 1
408 TH2F * fhM02ConNLocMax1[7][2] ; //! M02 for Mass around close to 0, N Local Maxima = 1
410 TH2F * fhM02Pi0NLocMax2[7][2] ; //! M02 for Mass around pi0, N Local Maxima = 2
411 TH2F * fhM02EtaNLocMax2[7][2] ; //! M02 for Mass around eta, N Local Maxima = 2
412 TH2F * fhM02ConNLocMax2[7][2] ; //! M02 for Mass around close to 0, N Local Maxima = 2
414 TH2F * fhM02Pi0NLocMaxN[7][2] ; //! M02 for Mass around pi0, N Local Maxima > 2
415 TH2F * fhM02EtaNLocMaxN[7][2] ; //! M02 for Mass around eta, N Local Maxima > 2
416 TH2F * fhM02ConNLocMaxN[7][2] ; //! M02 for Mass around close to 0, N Local Maxima > 2
418 TH2F * fhMassPi0NLocMax1[7][2] ; //! Mass for selected pi0, N Local Maxima = 1
419 TH2F * fhMassEtaNLocMax1[7][2] ; //! Mass for selected around eta, N Local Maxima = 1
420 TH2F * fhMassConNLocMax1[7][2] ; //! Mass for selected around close to 0, N Local Maxima = 1
422 TH2F * fhMassPi0NLocMax2[7][2] ; //! Mass for selected around pi0, N Local Maxima = 2
423 TH2F * fhMassEtaNLocMax2[7][2] ; //! Mass for selected around eta, N Local Maxima = 2
424 TH2F * fhMassConNLocMax2[7][2] ; //! Mass for selected around close to 0, N Local Maxima = 2
426 TH2F * fhMassPi0NLocMaxN[7][2] ; //! Mass for selected around pi0, N Local Maxima > 2
427 TH2F * fhMassEtaNLocMaxN[7][2] ; //! Mass for selected around eta, N Local Maxima > 2
428 TH2F * fhMassConNLocMaxN[7][2] ; //! Mass for selected around close to 0, N Local Maxima > 2
430 TH2F * fhNCellPi0NLocMax1[7][2] ; //! n cells for selected around pi0, N Local Maxima = 1
431 TH2F * fhNCellEtaNLocMax1[7][2] ; //! n cells for selected around eta, N Local Maxima = 1
432 TH2F * fhNCellPi0NLocMax2[7][2] ; //! n cells for selected around pi0, N Local Maxima = 2
433 TH2F * fhNCellEtaNLocMax2[7][2] ; //! n cells for selected around eta, N Local Maxima = 2
434 TH2F * fhNCellPi0NLocMaxN[7][2] ; //! n cells for selected around pi0, N Local Maxima > 2
435 TH2F * fhNCellEtaNLocMaxN[7][2] ; //! n cells for selected around eta, N Local Maxima > 2
437 TH2F * fhMassAfterCutsNLocMax1[7][2] ; //! Mass after M02, asymmetry cuts for MC part, N Local Maxima = 1
438 TH2F * fhMassAfterCutsNLocMax2[7][2] ; //! Mass after M02, asymmetry cuts for MC part, N Local Maxima = 2
439 TH2F * fhMassAfterCutsNLocMaxN[7][2] ; //! Mass after M02, asymmetry cuts for MC part, N Local Maxima > 2
441 TH2F * fhMassSplitEAfterCutsNLocMax1[7][2] ; //! Split Inv Mass vs E1+E2, NLM=1, after M02, asymmetry cuts, different MC particle types, track matching on/off
442 TH2F * fhMassSplitEAfterCutsNLocMax2[7][2] ; //! Split Inv Mass vs E1+E2, NLM=2, after M02, asymmetry cuts, different MC particle types, track matching on/off
443 TH2F * fhMassSplitEAfterCutsNLocMaxN[7][2] ; //! Split Inv Mass vs E1+E2, NLM>2, after M02, asymmetry cuts, different MC particle types, track matching on/off
445 TH2F * fhMassSplitEPi0NLocMax1[7][2] ; //! Split Inv Mass vs E1+E2, NLM=1, after pi0 selection, different MC particle types, track matching on/off
446 TH2F * fhMassSplitEPi0NLocMax2[7][2] ; //! Split Inv Mass vs E1+E2, NLM=2, after pi0 selection, different MC particle types, track matching on/off
447 TH2F * fhMassSplitEPi0NLocMaxN[7][2] ; //! Split Inv Mass vs E1+E2, NLM>2, after pi0 selection, different MC particle types, track matching on/off
449 TH2F * fhAsyPi0NLocMax1[7][2] ; //! Asy for Mass around pi0, N Local Maxima = 1
450 TH2F * fhAsyEtaNLocMax1[7][2] ; //! Asy for Mass around eta, N Local Maxima = 1
451 TH2F * fhAsyConNLocMax1[7][2] ; //! Asy for Mass around close to 0, N Local Maxima = 1
453 TH2F * fhAsyPi0NLocMax2[7][2] ; //! Asy for Mass around pi0, N Local Maxima = 2
454 TH2F * fhAsyEtaNLocMax2[7][2] ; //! Asy for Mass around eta, N Local Maxima = 2
455 TH2F * fhAsyConNLocMax2[7][2] ; //! Asy for Mass around close to 0, N Local Maxima = 2
457 TH2F * fhAsyPi0NLocMaxN[7][2] ; //! Asy for Mass around pi0, N Local Maxima > 2
458 TH2F * fhAsyEtaNLocMaxN[7][2] ; //! Asy for Mass around eta, N Local Maxima > 2
459 TH2F * fhAsyConNLocMaxN[7][2] ; //! Asy for Mass around close to 0, N Local Maxima > 2
461 TH2F * fhSplitEFractionNLocMax1[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1
462 TH2F * fhSplitEFractionNLocMax2[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2
463 TH2F * fhSplitEFractionNLocMaxN[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2
465 TH2F * fhSplitEFractionAfterCutsNLocMax1[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1, after M02 and asymmetry cut
466 TH2F * fhSplitEFractionAfterCutsNLocMax2[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2, after M02 and asymmetry cut
467 TH2F * fhSplitEFractionAfterCutsNLocMaxN[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2, after M02 and asymmetry cut
469 TH2F * fhMassSplitEFractionNLocMax1Ebin[7][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
470 TH2F * fhMassSplitEFractionNLocMax2Ebin[7][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
471 TH2F * fhMassSplitEFractionNLocMaxNEbin[7][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
473 TH2F * fhAnglePairNLocMax1[7][2] ; //! pair opening angle vs E
474 TH2F * fhAnglePairNLocMax2[7][2] ; //! pair opening angle vs E
475 TH2F * fhAnglePairNLocMaxN[7][2] ; //! pair opening angle vs E
477 TH2F * fhAnglePairAfterCutsNLocMax1[7][2] ; //! pair opening angle vs E
478 TH2F * fhAnglePairAfterCutsNLocMax2[7][2] ; //! pair opening angle vs E
479 TH2F * fhAnglePairAfterCutsNLocMaxN[7][2] ; //! pair opening angle vs E
481 TH2F * fhAnglePairPi0NLocMax1[7][2] ; //! pair opening angle vs E
482 TH2F * fhAnglePairPi0NLocMax2[7][2] ; //! pair opening angle vs E
483 TH2F * fhAnglePairPi0NLocMaxN[7][2] ; //! pair opening angle vs E
485 TH2F * fhAnglePairMassNLocMax1[7][2] ; //! pair opening angle vs Mass for E > 7 GeV
486 TH2F * fhAnglePairMassNLocMax2[7][2] ; //! pair opening angle vs Mass for E > 7 GeV
487 TH2F * fhAnglePairMassNLocMaxN[7][2] ; //! pair opening angle vs Mass for E > 7 GeV
489 TH2F * fhAnglePairM02NLocMax1[7][2] ; //! pair opening angle vs M02 for E > 7 GeV
490 TH2F * fhAnglePairM02NLocMax2[7][2] ; //! pair opening angle vs M02 for E > 7 GeV
491 TH2F * fhAnglePairM02NLocMaxN[7][2] ; //! pair opening angle vs M02 for E > 7 GeV
493 TH2F * fhAnglePairPrimPi0RecoNLocMax1; //! pair opening angle pi0 generated/reconstructed vs E
494 TH2F * fhAnglePairPrimPi0RecoNLocMax2; //! pair opening angle pi0 generated/reconstructed vs E
495 TH2F * fhAnglePairPrimPi0RecoNLocMaxN; //! pair opening angle pi0 generated/reconstructed vs E
497 TH2F * fhAnglePairPrimPi0vsRecoNLocMax1; //! pair opening angle pi0 generated vs reconstructed
498 TH2F * fhAnglePairPrimPi0vsRecoNLocMax2; //! pair opening angle pi0 generated vs reconstructed
499 TH2F * fhAnglePairPrimPi0vsRecoNLocMaxN; //! pair opening angle pi0 generated vs reconstructed
501 TH2F * fhAnglePairOverM02NLocMax1[7][2]; //! pair opening angle / m02 vs E, NLM=1
502 TH2F * fhAnglePairOverM02NLocMax2[7][2]; //! pair opening angle / m02 vs E, NLM=2
503 TH2F * fhAnglePairOverM02NLocMaxN[7][2]; //! pair opening angle / m02 vs E, NLM=N
505 TH2F * fhAnglePairOverM02NLocMax1Overlap0[7][2];//! pair opening angle / m02 vs E, NLM=1
506 TH2F * fhAnglePairOverM02NLocMax2Overlap0[7][2];//! pair opening angle / m02 vs E, NLM=2
507 TH2F * fhAnglePairOverM02NLocMaxNOverlap0[7][2];//! pair opening angle / m02 vs E, NLM=N
509 TH2F * fhAnglePairPrimPi0OverM02NLocMax1; //! pair opening angle / m02 vs E, NLM=1, prim pi0
510 TH2F * fhAnglePairPrimPi0OverM02NLocMax2; //! pair opening angle / m02 vs E, NLM=2, prim pi0
511 TH2F * fhAnglePairPrimPi0OverM02NLocMaxN; //! pair opening angle / m02 vs E, NLM=N, prim pi0
513 TH2F * fhArmNLocMax1[7][4] ; //! Armenteros of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types
514 TH2F * fhArmNLocMax2[7][4] ; //! Armenteros of 2 cells local maxima vs E, 1-6 for different MC particle types
515 TH2F * fhArmNLocMaxN[7][4] ; //! Armenteros of >2 cells local maxima vs E, 1-6 for different MC particle types
517 TH2F * fhArmAfterCutsNLocMax1[7][4] ; //! Armenteros after M02, asymmetry cuts for pi0, N Local Maxima = 1
518 TH2F * fhArmAfterCutsNLocMax2[7][4] ; //! Armenteros after M02, asymmetry cuts for pi0, N Local Maxima = 2
519 TH2F * fhArmAfterCutsNLocMaxN[7][4] ; //! Armenteros after M02, asymmetry cuts for pi0, N Local Maxima > 2
521 TH2F * fhArmPi0NLocMax1[7][4] ; //! Armenteros for selected pi0, N Local Maxima = 1
522 TH2F * fhArmPi0NLocMax2[7][4] ; //! Armenteros for selected pi0, N Local Maxima = 2
523 TH2F * fhArmPi0NLocMaxN[7][4] ; //! Armenteros for selected pi0, N Local Maxima > 2
525 TH2F * fhCosThStarNLocMax1[7][2] ; //! cos(theta^star) vs E, NLM=1
526 TH2F * fhCosThStarNLocMax2[7][2] ; //! cos(theta^star) vs E, NLM=2
527 TH2F * fhCosThStarNLocMaxN[7][2] ; //! cos(theta^star) vs E, NLM>2
529 TH2F * fhCosThStarAfterCutsNLocMax1[7][2] ; //! cos(theta^star) vs E, after M02, asymmetry cuts, NLM=1
530 TH2F * fhCosThStarAfterCutsNLocMax2[7][2] ; //! cos(theta^star) vs E, after M02, asymmetry cuts, NLM=2
531 TH2F * fhCosThStarAfterCutsNLocMaxN[7][2] ; //! cos(theta^star) vs E, after M02, asymmetry cuts, NLM>2
533 TH2F * fhCosThStarPi0NLocMax1[7][2] ; //! cos(theta^star) vs E, after M02, asymmetry and pi0 mass cuts, NLM=1
534 TH2F * fhCosThStarPi0NLocMax2[7][2] ; //! cos(theta^star) vs E, after M02, asymmetry and pi0 mass cuts, NLM=2
535 TH2F * fhCosThStarPi0NLocMaxN[7][2] ; //! cos(theta^star) vs E, after M02, asymmetry and pi0 mass cuts, NLM>2
537 TH2F * fhTrackMatchedDEtaNLocMax1[7] ; //! Eta distance between track and cluster vs cluster E, 1 local maximum
538 TH2F * fhTrackMatchedDPhiNLocMax1[7] ; //! Phi distance between track and cluster vs cluster E, 1 local maximum
539 TH2F * fhTrackMatchedDEtaNLocMax2[7] ; //! Eta distance between track and cluster vs cluster E, 2 local maximum
540 TH2F * fhTrackMatchedDPhiNLocMax2[7] ; //! Phi distance between track and cluster vs cluster E, 2 local maximum
541 TH2F * fhTrackMatchedDEtaNLocMaxN[7] ; //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
542 TH2F * fhTrackMatchedDPhiNLocMaxN[7] ; //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
544 TH2F * fhTrackMatchedDEtaNLocMax1Pos[7] ; //! Eta distance between track and cluster vs cluster E, 1 local maximum
545 TH2F * fhTrackMatchedDPhiNLocMax1Pos[7] ; //! Phi distance between track and cluster vs cluster E, 1 local maximum
546 TH2F * fhTrackMatchedDEtaNLocMax2Pos[7] ; //! Eta distance between track and cluster vs cluster E, 2 local maximum
547 TH2F * fhTrackMatchedDPhiNLocMax2Pos[7] ; //! Phi distance between track and cluster vs cluster E, 2 local maximum
548 TH2F * fhTrackMatchedDEtaNLocMaxNPos[7] ; //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
549 TH2F * fhTrackMatchedDPhiNLocMaxNPos[7] ; //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
551 TH2F * fhTrackMatchedDEtaNLocMax1Neg[7] ; //! Eta distance between track and cluster vs cluster E, 1 local maximum
552 TH2F * fhTrackMatchedDPhiNLocMax1Neg[7] ; //! Phi distance between track and cluster vs cluster E, 1 local maximum
553 TH2F * fhTrackMatchedDEtaNLocMax2Neg[7] ; //! Eta distance between track and cluster vs cluster E, 2 local maximum
554 TH2F * fhTrackMatchedDPhiNLocMax2Neg[7] ; //! Phi distance between track and cluster vs cluster E, 2 local maximum
555 TH2F * fhTrackMatchedDEtaNLocMaxNNeg[7] ; //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
556 TH2F * fhTrackMatchedDPhiNLocMaxNNeg[7] ; //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
558 TH2F * fhCentralityPi0NLocMax1 ; //! Centrality for selected pi0, N Local Maxima = 1
559 TH2F * fhCentralityEtaNLocMax1 ; //! Centrality for selected eta, N Local Maxima = 1
560 TH2F * fhCentralityPi0NLocMax2 ; //! Centrality for selected pi0, N Local Maxima = 2
561 TH2F * fhCentralityEtaNLocMax2 ; //! Centrality for selected eta, N Local Maxima = 2
562 TH2F * fhCentralityPi0NLocMaxN ; //! Centrality for selected pi0, N Local Maxima > 2
563 TH2F * fhCentralityEtaNLocMaxN ; //! Centrality for selected eta, N Local Maxima > 2
565 TH2F * fhEventPlanePi0NLocMax1 ; //! Event plane for selected pi0, N Local Maxima = 1
566 TH2F * fhEventPlaneEtaNLocMax1 ; //! Event plane for selected eta, N Local Maxima = 1
567 TH2F * fhEventPlanePi0NLocMax2 ; //! Event plane for selected pi0, N Local Maxima = 2
568 TH2F * fhEventPlaneEtaNLocMax2 ; //! Event plane for selected eta, N Local Maxima = 2
569 TH2F * fhEventPlanePi0NLocMaxN ; //! Event plane for selected pi0, N Local Maxima > 2
570 TH2F * fhEventPlaneEtaNLocMaxN ; //! Event plane for selected eta, N Local Maxima > 2
572 TH2F * fhClusterEtaPhiNLocMax1 ; //! Eta vs Phi of clusters with N Local Maxima = 1, E > 8 GeV
573 TH2F * fhClusterEtaPhiNLocMax2 ; //! Eta vs Phi of clusters with N Local Maxima = 2, E > 8 GeV
574 TH2F * fhClusterEtaPhiNLocMaxN ; //! Eta vs Phi of clusters with N Local Maxima > 2, E > 8 GeV
575 TH2F * fhPi0EtaPhiNLocMax1 ; //! Eta vs Phi of pi0's with N Local Maxima = 1, E > 8 GeV
576 TH2F * fhPi0EtaPhiNLocMax2 ; //! Eta vs Phi of pi0's with N Local Maxima = 2, E > 8 GeV
577 TH2F * fhPi0EtaPhiNLocMaxN ; //! Eta vs Phi of pi0's with N Local Maxima > N, E > 8 GeV
578 TH2F * fhEtaEtaPhiNLocMax1 ; //! Eta vs Phi of eta's with N Local Maxima = 1, E > 8 GeV
579 TH2F * fhEtaEtaPhiNLocMax2 ; //! Eta vs Phi of eta's with N Local Maxima = 2, E > 8 GeV
580 TH2F * fhEtaEtaPhiNLocMaxN ; //! Eta vs Phi of eta's with N Local Maxima > N, E > 8 GeV
582 TH2F * fhPi0CellE[3] ; //! pi0's energy vs cluster cell energy with NLM = 1, = 2, > 2
583 TH2F * fhPi0CellEFrac[3] ; //! pi0's energy vs cluster cell energy fraction with NLM = 1, = 2, > 2
584 TH2F * fhPi0CellLogEFrac[3] ; //! pi0's energy vs cluster log cell energy fraction with NLM = 1, = 2, > 2
585 TH2F * fhPi0CellEMaxEMax2Frac [3]; //! pi0's energy vs fraction of 2 main maxima energy with NLM = 1, = 2, > 2
586 TH2F * fhPi0CellEMaxClusterFrac [3]; //! pi0's energy vs energy fraction of main LM and cluster energy with NLM = 1, = 2, > 2
587 TH2F * fhPi0CellEMax2ClusterFrac[3]; //! pi0's energy vs energy fraction of second LM and cluster energy with NLM = 1, = 2, > 2
588 TH2F * fhPi0CellEMaxFrac [3]; //! pi0's energy vs energy fraction of main LM and cluster cell energy with NLM = 1, = 2, > 2
589 TH2F * fhPi0CellEMax2Frac [3]; //! pi0's energy vs energy fraction of second LM and cluster cell energy with NLM = 1, = 2, > 2
591 TH2F * fhM02WeightPi0[3][10] ; //! M02 for selected pi0 with different weight, with NLM = 1, = 2, > 2
592 TH2F * fhM02ECellCutPi0[3][10] ; //! M02 for selected pi0 with different cut on cell energy, with NLM = 1, = 2, > 2
594 TH2F * fhPi0EPairDiffTimeNLM1; //! E vs Pair of clusters time difference vs E, for selected pi0, NLM=1
595 TH2F * fhPi0EPairDiffTimeNLM2; //! E vs Pair of clusters time difference vs E, for selected pi0, NLM=2
596 TH2F * fhPi0EPairDiffTimeNLMN; //! E vs Pair of clusters time difference vs E, for selected pi0, NLM>2
597 TH2F * fhEtaEPairDiffTimeNLM1; //! E vs Pair of clusters time difference vs E, for selected eta, NLM=1
598 TH2F * fhEtaEPairDiffTimeNLM2; //! E vs Pair of clusters time difference vs E, for selected eta, NLM=2
599 TH2F * fhEtaEPairDiffTimeNLMN; //! E vs Pair of clusters time difference vs E, for selected eta, NLM>2
601 TH2F * fhMCEM02Overlap0[3][7]; //! E vs M02 for different MC origin, no other MC particles contributes, neutral cluster
602 TH2F * fhMCEM02Overlap1[3][7]; //! E vs M02 for different MC origin, 1 other MC particles contributes, neutral cluster
603 TH2F * fhMCEM02OverlapN[3][7]; //! E vs M02 for different MC origin, N other MC particles contributes, neutral cluster
604 TH2F * fhMCEM02Overlap0Match[3][7]; //! E vs M02 for different MC origin, no other MC particles contributes, charged cluster
605 TH2F * fhMCEM02Overlap1Match[3][7]; //! E vs M02 for different MC origin, 1 other MC particles contributes, charged cluster
606 TH2F * fhMCEM02OverlapNMatch[3][7]; //! E vs M02 for different MC origin, N other MC particles contributes, charged cluster
608 TH2F * fhMCEMassOverlap0[3][7]; //! E vs Mass for different MC origin, no other MC particles contributes, neutral cluster
609 TH2F * fhMCEMassOverlap1[3][7]; //! E vs Mass for different MC origin, 1 other MC particles contributes, neutral cluster
610 TH2F * fhMCEMassOverlapN[3][7]; //! E vs Mass for different MC origin, N other MC particles contributes, neutral cluster
611 TH2F * fhMCEMassOverlap0Match[3][7]; //! E vs Mass for different MC origin, no other MC particles contributes, charged cluster
612 TH2F * fhMCEMassOverlap1Match[3][7]; //! E vs Mass for different MC origin, 1 other MC particles contributes, charged cluster
613 TH2F * fhMCEMassOverlapNMatch[3][7]; //! E vs Mass for different MC origin, N other MC particles contributes, charged cluster
615 TH2F * fhMCESplitEFracOverlap0[3][7]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, no other MC particles contributes, neutral cluster
616 TH2F * fhMCESplitEFracOverlap1[3][7]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, 1 other MC particles contributes, neutral cluster
617 TH2F * fhMCESplitEFracOverlapN[3][7]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, N other MC particles contributes, neutral cluster
618 TH2F * fhMCESplitEFracOverlap0Match[3][7]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, no other MC particles contributes, charged cluster
619 TH2F * fhMCESplitEFracOverlap1Match[3][7]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, 1 other MC particles contributes, charged cluster
620 TH2F * fhMCESplitEFracOverlapNMatch[3][7]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, N other MC particles contributes, charged cluster
622 TH2F * fhMCEAsymOverlap0[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contributes, neutral cluster
623 TH2F * fhMCEAsymOverlap1[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, 1 other MC particles contributes, neutral cluster
624 TH2F * fhMCEAsymOverlapN[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, N other MC particles contributes, neutral cluster
625 TH2F * fhMCEAsymOverlap0Match[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contributes, charged cluster
626 TH2F * fhMCEAsymOverlap1Match[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, 1 other MC particles contributes, charged cluster
627 TH2F * fhMCEAsymOverlapNMatch[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, N other MC particles contributes, charged cluster
629 TH2F * fhMCENCellOverlap0[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contributes, neutral cluster
630 TH2F * fhMCENCellOverlap1[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, 1 other MC particles contributes, neutral cluster
631 TH2F * fhMCENCellOverlapN[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, N other MC particles contributes, neutral cluster
632 TH2F * fhMCENCellOverlap0Match[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contributes, charged cluster
633 TH2F * fhMCENCellOverlap1Match[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, 1 other MC particles contributes, charged cluster
634 TH2F * fhMCENCellOverlapNMatch[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, N other MC particles contributes, charged cluster
636 TH2F * fhMCEEpriOverlap0[3][7]; //! E reco vs primary for different MC origin, no other MC particles contributes, neutral cluster
637 TH2F * fhMCEEpriOverlap1[3][7]; //! E reco vs primary for different MC origin, 1 other MC particles contributes, neutral cluster
638 TH2F * fhMCEEpriOverlapN[3][7]; //! E reco vs primary for different MC origin, N other MC particles contributes, neutral cluster
639 TH2F * fhMCEEpriOverlap0Match[3][7]; //! E reco vs primary for different MC origin, no other MC particles contributes, charged cluster
640 TH2F * fhMCEEpriOverlap1Match[3][7]; //! E reco vs primary for different MC origin, 1 other MC particles contributes, charged cluster
641 TH2F * fhMCEEpriOverlapNMatch[3][7]; //! E reco vs primary for different MC origin, N other MC particles contributes, charged cluster
643 TH2F * fhMCEEpriOverlap0IdPi0[3][7]; //! E reco vs primary for different MC origin, no other MC particles contributes, neutral cluster, neutral clusters id as pi0
644 TH2F * fhMCEEpriOverlap1IdPi0[3][7]; //! E reco vs primary for different MC origin, 1 other MC particles contributes, neutral cluster, neutral clusters id as pi0
645 TH2F * fhMCEEpriOverlapNIdPi0[3][7]; //! E reco vs primary for different MC origin, 1 other MC particles contributes, neutral cluster, neutral clusters is as pi0
647 TH2F * fhMCPi0MassM02Overlap0[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, no other MC particles contributes, neutral cluster, 4 E bins
648 TH2F * fhMCPi0MassM02Overlap1[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, 1 other MC particles contributes, neutral cluster, 4 E bins
649 TH2F * fhMCPi0MassM02OverlapN[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, N other MC particles contributes, neutral cluster, 4 E bins
650 TH2F * fhMCPi0MassM02Overlap0Match[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, no other MC particles contributes, charged cluster, 4 E bins
651 TH2F * fhMCPi0MassM02Overlap1Match[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, 1 other MC particles contributes, charged cluster, 4 E bins
652 TH2F * fhMCPi0MassM02OverlapNMatch[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, N other MC particles contributes, charged cluster, 4 E bins
654 TH2F * fhMCENOverlaps[3][7]; //! E vs number of Overlaps in MC, neutral cluster
655 TH2F * fhMCENOverlapsMatch[3][7]; //! E vs number of Overlaps in MC, charged cluster
657 TH2F * fhMCPi0HighNLMPair; //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons
658 TH2F * fhMCPi0LowNLMPair; //! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy correspond to 2 photons
659 TH2F * fhMCPi0AnyNLMPair; //! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs correspond to 2 photons
660 TH2F * fhMCPi0NoneNLMPair; //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons
661 // No match between highest energy local maxima and highest energy MC particle
662 TH2F * fhMCPi0HighNLMPairNoMCMatch; //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons
663 TH2F * fhMCPi0LowNLMPairNoMCMatch; //! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy correspond to 2 photons
664 TH2F * fhMCPi0AnyNLMPairNoMCMatch; //! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs correspond to 2 photons
665 TH2F * fhMCPi0NoneNLMPairNoMCMatch; //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons
667 TH2F * fhMCPi0HighNLMPairOverlap; //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons, overlap
668 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
669 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
670 TH2F * fhMCPi0NoneNLMPairOverlap; //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons, overlap
671 // No match between highest energy local maxima and highest energy MC particle
672 TH2F * fhMCPi0HighNLMPairNoMCMatchOverlap; //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons, overlap
673 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
674 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
675 TH2F * fhMCPi0NoneNLMPairNoMCMatchOverlap; //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons, overlap
677 TH2F * fhMCPi0DecayPhotonHitHighLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
678 TH2F * fhMCPi0DecayPhotonAdjHighLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
679 TH2F * fhMCPi0DecayPhotonHitOtherLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
680 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
681 TH2F * fhMCPi0DecayPhotonAdjacent; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM
682 TH2F * fhMCPi0DecayPhotonHitNoLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas
684 TH2F * fhMCPi0DecayPhotonHitHighLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maxima, overlap
685 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
686 TH2F * fhMCPi0DecayPhotonHitOtherLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
687 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
688 TH2F * fhMCPi0DecayPhotonAdjacentOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM, overlap
689 TH2F * fhMCPi0DecayPhotonHitNoLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas, overlap
691 TH2F * fhMCPi0DecayPhotonHitHighLMDiffELM1[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
692 TH2F * fhMCPi0DecayPhotonAdjHighLMDiffELM1[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
693 TH2F * fhMCPi0DecayPhotonHitOtherLMDiffELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
694 TH2F * fhMCPi0DecayPhotonAdjOtherLMDiffELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high
696 TH2F * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM1[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
697 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
698 TH2F * fhMCPi0DecayPhotonHitOtherLMOverlapDiffELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
699 TH2F * fhMCPi0DecayPhotonAdjOtherLMOverlapDiffELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high, overlap
701 TH2F * fhMCPi0DecayPhotonHitHighLMDiffELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
702 TH2F * fhMCPi0DecayPhotonAdjHighLMDiffELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
703 TH2F * fhMCPi0DecayPhotonHitOtherLMDiffELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
704 TH2F * fhMCPi0DecayPhotonAdjOtherLMDiffELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high
706 TH2F * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM2[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
707 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
708 TH2F * fhMCPi0DecayPhotonHitOtherLMOverlapDiffELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
709 TH2F * fhMCPi0DecayPhotonAdjOtherLMOverlapDiffELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high, overlap
712 TH2F * fhMCPi0DecayPhotonHitHighLMDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
713 TH2F * fhMCPi0DecayPhotonAdjHighLMDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
714 TH2F * fhMCPi0DecayPhotonHitOtherLMDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
715 TH2F * fhMCPi0DecayPhotonAdjOtherLMDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high
717 TH2F * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
718 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
719 TH2F * fhMCPi0DecayPhotonHitOtherLMOverlapDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
720 TH2F * fhMCPi0DecayPhotonAdjOtherLMOverlapDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high, overlap
722 TH2F * fhMCPi0DecayPhotonHitHighLMDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
723 TH2F * fhMCPi0DecayPhotonAdjHighLMDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
724 TH2F * fhMCPi0DecayPhotonHitOtherLMDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
725 TH2F * fhMCPi0DecayPhotonAdjOtherLMDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high
727 TH2F * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
728 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
729 TH2F * fhMCPi0DecayPhotonHitOtherLMOverlapDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
730 TH2F * fhMCPi0DecayPhotonAdjOtherLMOverlapDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high, overlap
733 TH2F * fhMCPi0DecayPhotonHitHighLMMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
734 TH2F * fhMCPi0DecayPhotonAdjHighLMMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
735 TH2F * fhMCPi0DecayPhotonHitOtherLMMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
736 TH2F * fhMCPi0DecayPhotonAdjOtherLMMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high
737 TH2F * fhMCPi0DecayPhotonAdjacentMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM
738 TH2F * fhMCPi0DecayPhotonHitNoLMMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas
740 TH2F * fhMCPi0DecayPhotonHitHighLMOverlapMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maxima, overlap
741 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
742 TH2F * fhMCPi0DecayPhotonHitOtherLMOverlapMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
743 TH2F * fhMCPi0DecayPhotonAdjOtherLMOverlapMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high, overlap
744 TH2F * fhMCPi0DecayPhotonAdjacentOverlapMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM, overlap
745 TH2F * fhMCPi0DecayPhotonHitNoLMOverlapMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas, overlap
748 TH2F * fhMCEOverlapType; //! what particles overlap with pi0, neutral clusters
749 TH2F * fhMCEOverlapTypeMatch; //! what particles overlap with pi0, charged clusters
751 TH2F * fhMassBadDistClose[3]; //! split mass of clusters with second LM close to bad channel
752 TH2F * fhM02BadDistClose[3]; //! m02 of clusters with second LM close to bad channel
753 TH2F * fhMassOnBorder[3]; //! split mass of clusters with second LM on EMCAL border
754 TH2F * fhM02OnBorder[3]; //! m02 of clusters with second LM close to EMCAL border
756 AliAnaInsideClusterInvariantMass( const AliAnaInsideClusterInvariantMass & split) ; // cpy ctor
757 AliAnaInsideClusterInvariantMass & operator = (const AliAnaInsideClusterInvariantMass & split) ; // cpy assignment
759 ClassDef(AliAnaInsideClusterInvariantMass,28)
763 #endif //ALIANAINSIDECLUSTERINVARIANTMASS_H