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 FillNLMDiffCutHistograms(AliVCluster *cluster, AliVCaloCells *cells, Bool_t matched);
83 void FillNCellHistograms(Int_t ncells, Float_t energy, Int_t nMax,
84 Bool_t matched, Int_t mcindex,
85 Float_t mass , Float_t l0);
87 void FillTrackMatchingHistograms(AliVCluster * cluster,Int_t nMax, Int_t mcindex);
89 void FillHistograms1(Float_t en, Float_t e1, Float_t e2,
90 Int_t nMax, Float_t mass, Float_t l0,
91 Float_t eta, Float_t phi,
92 Bool_t matched, Int_t mcindex);
95 void FillHistograms2(Float_t en, Float_t eprim,
96 Float_t e1, Float_t e2, Int_t nMax,
97 Float_t mass, Float_t l0,
98 Bool_t matched, Int_t mcindex);
100 void FillIdPi0Histograms(Float_t en, Float_t e1, Float_t e2,
101 Int_t nc, Int_t nMax, Float_t t12diff,
102 Float_t mass, Float_t l0,
103 Float_t eta, Float_t phi,
104 Bool_t matched, Int_t mcindex);
106 void FillIdEtaHistograms(Float_t en, Float_t e1, Float_t e2,
107 Int_t nc, Int_t nMax, Float_t t12diff,
108 Float_t mass, Float_t l0,
109 Float_t eta, Float_t phi,
110 Bool_t matched, Int_t mcindex);
112 void FillIdConvHistograms(Float_t en, Int_t nMax, Float_t asym,
113 Float_t mass, Float_t l0,
114 Bool_t matched, Int_t mcindex);
118 void InitParameters();
120 void MakeAnalysisFillHistograms() ;
122 void Print(const Option_t * opt) const;
124 void SetCalorimeter(TString & det) { fCalorimeter = det ; }
126 void SetMinNCells(Int_t cut) { fMinNCells = cut ; }
128 void SetMinBadChannelDistance(Float_t cut) { fMinBadDist = cut ; }
130 void SetWCorrectionParameter(Int_t i, Float_t p = 0.07) { if( i<2 ) fWSimu[i] = p; }
132 void SwitchOnFillAngleHistograms() { fFillAngleHisto = kTRUE ; }
133 void SwitchOffFillAngleHistograms() { fFillAngleHisto = kFALSE ; }
135 void SwitchOnFillArmenterosHistograms() { fFillArmenterosHisto = kTRUE ; }
136 void SwitchOffFillArmenterosHistograms() { fFillArmenterosHisto = kFALSE ; }
138 void SwitchOnFillThetaStarHistograms() { fFillThetaStarHisto = kTRUE ; }
139 void SwitchOffFillThetaStarHistograms() { fFillThetaStarHisto = kFALSE ; }
141 void SwitchOnFillExtraSSHistograms() { fFillSSExtraHisto = kTRUE ; }
142 void SwitchOffFillExtraSSHistograms() { fFillSSExtraHisto = kFALSE ; }
144 void SwitchOnFillHighMultHistograms() { fFillHighMultHisto = kTRUE ; }
145 void SwitchOffFillHighMultHistograms() { fFillHighMultHisto = kFALSE ; }
147 void SwitchOnFillIdConvHistograms() { fFillIdConvHisto = kTRUE ; }
148 void SwitchOffFillIdConvHistograms() { fFillIdConvHisto = kFALSE ; }
150 void SwitchOnFillIdEtaHistograms() { fFillIdEtaHisto = kTRUE ; }
151 void SwitchOffFillIdEtaHistograms() { fFillIdEtaHisto = kFALSE ; }
153 void SwitchOnFillTMHistograms() { fFillTMHisto = kTRUE ; }
154 void SwitchOffFillTMHistograms() { fFillTMHisto = kFALSE ; }
156 void SwitchOnFillTMResidualHistograms() { fFillTMResidualHisto = kTRUE ; }
157 void SwitchOffFillTMResidualHistograms() { fFillTMResidualHisto = kFALSE ; }
159 void SwitchOnFillMCPrimaryHistograms() { fFillMCHisto = kTRUE ; }
160 void SwitchOffFillMCPrimaryHistograms() { fFillMCHisto = kFALSE ; }
162 void SwitchOnFillSSWeightHistograms() { fFillSSWeightHisto = kTRUE ; }
163 void SwitchOffFillSSWeightHistograms() { fFillSSWeightHisto = kFALSE ; }
165 void SwitchOnFillNLMDiffCutsHistograms() { fFillNLMDiffCutHisto = kTRUE ; }
166 void SwitchOffFillNLMDiffCutsHistograms() { fFillNLMDiffCutHisto = kFALSE ; }
168 void SwitchOnFillEbinHistograms() { fFillEbinHisto = kTRUE ; }
169 void SwitchOffFillEbinHistograms() { fFillEbinHisto = kFALSE ; }
171 void SwitchOnFillMCOverlapHistograms() { fFillMCOverlapHisto = kTRUE ; }
172 void SwitchOffFillMCOverlapHistograms() { fFillMCOverlapHisto = kFALSE ; }
174 void SwitchOnFillNCellHistograms() { fFillNCellHisto = kTRUE ; }
175 void SwitchOffFillNCellHistograms() { fFillNCellHisto = kFALSE ; }
177 void SwitchOnSplitClusterDistToBad() { fCheckSplitDistToBad = kTRUE ; }
178 void SwitchOffSplitClusterDistToBad() { fCheckSplitDistToBad = kFALSE ; }
180 void SetNWeightForShowerShape(Int_t n) { fSSWeightN = n ; }
181 void SetWeightForShowerShape(Int_t i, Float_t v)
182 { if (i < 20) fSSWeight[i] = v ; }
184 void SetNumberOfNLocMaxSettings(Int_t n) { fNLMSettingN = n ; }
185 void SetNLocMaxMinE (Int_t i, Float_t v) { if (i < 5) fNLMMinE [i] = v ; }
186 void SetNLocMaxMinDiff(Int_t i, Float_t v) { if (i < 5) fNLMMinDiff[i] = v ; }
189 void SetNECellCutForShowerShape(Int_t n) { fSSECellCutN = n ; }
190 void SetECellCutForShowerShape(Int_t i, Float_t v)
191 { if (i < 20) fSSECellCut[i] = v ; }
194 void RecalculateClusterShowerShapeParametersWithCellCut(const AliEMCALGeometry * geom, AliVCaloCells* cells, AliVCluster * cluster,
195 Float_t & l0, Float_t & l1,
196 Float_t & disp, Float_t & dEta, Float_t & dPhi,
197 Float_t & sEta, Float_t & sPhi, Float_t & sEtaPhi,Float_t eCellMin = 0.);
201 enum mcTypes { kmcPhoton = 1, kmcConversion = 2, kmcPi0 = 3, kmcPi0Conv = 4,
202 kmcEta = 5, kmcHadron = 6 };
206 TString fCalorimeter ; // Calorimeter where the gamma is searched
207 Int_t fMinNCells ; // Study clusters with ncells larger than cut
208 Float_t fMinBadDist ; // Minimal distance to bad channel to accept cluster
209 Float_t fHistoECut ; // Fixed E cut for some histograms
210 Bool_t fCheckSplitDistToBad; // Check the distance to bad channel and to EMCal borders of split clusters
212 Bool_t fFillAngleHisto; // Fill splitted clusters angle histograms
213 Bool_t fFillTMHisto ; // Fill track matching histos,
214 Bool_t fFillTMResidualHisto ; // Fill track matching histos, residuals
215 Bool_t fFillSSExtraHisto ; // Fill shower shape extra histos
216 Bool_t fFillMCHisto ; // Fill MC energy fraction histos
217 Bool_t fFillSSWeightHisto ; // Fill weigth histograms
218 Bool_t fFillNLMDiffCutHisto ; // Fill NLM histograms for different settings
219 Bool_t fFillEbinHisto ; // Fill E bin histograms
220 Bool_t fFillMCOverlapHisto ; // Fill MC particles overlap histograms
221 Bool_t fFillNCellHisto ; // Fill n cells in cluster dependent histograms
222 Bool_t fFillIdConvHisto ; // Fill histograms for clusters identified as conversion
223 Bool_t fFillIdEtaHisto ; // Fill histograms for clusters identified as Eta
224 Bool_t fFillHighMultHisto; // Fill centrality/event plane histograms
225 Bool_t fFillArmenterosHisto; // Fill armenteros type histo
226 Bool_t fFillThetaStarHisto; // Fill cosThetaStar histos
228 Float_t fSSWeight[20]; // List of weights to test
229 Int_t fSSWeightN; // Total number of weights to test
231 Float_t fSSECellCut[20]; // List of cell min energy cuts to test
232 Int_t fSSECellCutN; // Total number of cell min energy cuts to test
234 Float_t fNLMMinE [5]; // List of local maxima min energy
235 Float_t fNLMMinDiff[5]; // List of local maxima min difference cell energy
236 Int_t fNLMSettingN; // Total number of NLM settings to test
238 Float_t fWSimu[2]; // Constant and slope of the linear correction factor for the shower
239 // shape weight in simulation, about 1-0.07*w
243 TH2F * fhMassNLocMax1[7][2] ; //! Split Inv Mass vs cluster E, NLM=1, different MC particle types, track matching on/off
244 TH2F * fhMassNLocMax2[7][2] ; //! Split Inv Mass vs cluster E, NLM=2, different MC particle types, track matching on/off
245 TH2F * fhMassNLocMaxN[7][2] ; //! Split Inv Mass vs cluster E, NLM>2, different MC particle types, track matching on/off
247 TH2F * fhMassSplitENLocMax1[7][2] ; //! Split Inv Mass vs E1+E2, NLM=1, different MC particle types, track matching on/off
248 TH2F * fhMassSplitENLocMax2[7][2] ; //! Split Inv Mass vs E1+E2, NLM=2, different MC particle types, track matching on/off
249 TH2F * fhMassSplitENLocMaxN[7][2] ; //! Split Inv Mass vs E1+E2, NLM>2, different MC particle types, track matching on/off
251 TH2F * fhAsymNLocMax1[7][2] ; //! Asymmetry of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types
252 TH2F * fhAsymNLocMax2[7][2] ; //! Asymmetry of 2 cells local maxima vs E, 1-6 for different MC particle types
253 TH2F * fhAsymNLocMaxN[7][2] ; //! Asymmetry of >2 cells local maxima vs E, 1-6 for different MC particle types
255 TH2F * fhSplitEFractionvsAsyNLocMax1[2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1 vs |A|
256 TH2F * fhSplitEFractionvsAsyNLocMax2[2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2 vs |A|
257 TH2F * fhSplitEFractionvsAsyNLocMaxN[2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2 vs |A|
259 TH2F * fhMassAsyCutNLocMax1 ; //! Mass(E) asym selection, not matched, Mass of split clusters, NLM = 1
260 TH2F * fhMassAsyCutNLocMax2 ; //! Mass(E) asym selection, not matched, Mass of split clusters, NLM = 1
261 TH2F * fhMassAsyCutNLocMaxN ; //! Mass(E) asym selection, not matched, Mass of split clusters, NLM > 2
263 TH2F * fhM02AsyCutNLocMax1 ; //! M02(E) asym selection, not matched, M02, NLM = 1
264 TH2F * fhM02AsyCutNLocMax2 ; //! M02(E) asym selection, not matched, M02, NLM = 2
265 TH2F * fhM02AsyCutNLocMaxN ; //! M02(E) asym selection, not matched, M02, NLM > 2
267 TH2F * fhMassM02CutNLocMax1 ; //! Mass(E) M02 selection, not matched, Mass of split clusters, NLM = 1
268 TH2F * fhMassM02CutNLocMax2 ; //! Mass(E) M02 selection, not matched, Mass of split clusters, NLM = 1
269 TH2F * fhMassM02CutNLocMaxN ; //! Mass(E) M02 selection, not matched, Mass of split clusters, NLM > 2
271 TH2F * fhAsymM02CutNLocMax1 ; //! Asym(E) M02 selection, not matched, energy asymmetry of split clusters, NLM = 1
272 TH2F * fhAsymM02CutNLocMax2 ; //! Asym(E) M02 selection, not matched, energy asymmetry of split clusters, NLM = 2
273 TH2F * fhAsymM02CutNLocMaxN ; //! Asym(E) M02 selection, not matched, energy asymmetry of split clusters, NLM > 2
275 TH2F * fhMassEnCutNLocMax1 ; //! Mass(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
276 TH2F * fhMassEnCutNLocMax2 ; //! Mass(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
277 TH2F * fhMassEnCutNLocMaxN ; //! Mass(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
279 TH2F * fhM02EnCutNLocMax1 ; //! M02(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
280 TH2F * fhM02EnCutNLocMax2 ; //! M02(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
281 TH2F * fhM02EnCutNLocMaxN ; //! M02(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
283 TH2F * fhAsymEnCutNLocMax1 ; //! Asym(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
284 TH2F * fhAsymEnCutNLocMax2 ; //! Asym(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
285 TH2F * fhAsymEnCutNLocMaxN ; //! Asym(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
287 TH2F * fhSplitEFracEnCutNLocMax1 ; //! Split E fraction (E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
288 TH2F * fhSplitEFracEnCutNLocMax2 ; //! Split E fraction (E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
289 TH2F * fhSplitEFracEnCutNLocMaxN ; //! Split E fraction (E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
291 TH2F * fhMassSplitECutNLocMax1 ; //! 85% of split energy, not matched, Mass of split clusters, NLM = 1
292 TH2F * fhMassSplitECutNLocMax2 ; //! 85% of split energy, not matched, Mass of split clusters, NLM = 1
293 TH2F * fhMassSplitECutNLocMaxN ; //! 85% of split energy, not matched, Mass of split clusters, NLM > 2
295 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
296 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
297 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
299 TH2F * fhMassM02NLocMax1Ebin[4] ; //! Mass of splitted clusters when 1 local max vs M02, 4 E bins, neutral clusters
300 TH2F * fhMassM02NLocMax2Ebin[4] ; //! Mass of splitted clusters when 2 local max vs M02, 4 E bins, neutral clusters
301 TH2F * fhMassM02NLocMaxNEbin[4] ; //! Mass of splitted clusters when >2 local max vs M02, 4 E bins, neutral clusters
303 TH2F * fhMassAsyNLocMax1Ebin[4] ; //! Mass of Mass of splitted clusters when 1 local max vs asymmetry, 4 E bins, neutral clusters
304 TH2F * fhMassAsyNLocMax2Ebin[4] ; //! Mass of Mass of splitted clusters when 2 local max vs asymmetry, 4 E bins, neutral clusters
305 TH2F * fhMassAsyNLocMaxNEbin[4] ; //! Mass of Mass of splitted clusters when >2 local max vs asymmetry, 4 E bins, neutral clusters
307 TH2F * fhAsyMCGenRecoNLocMax1EbinPi0[4] ; //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when 1 local max, 4 E bins, neutral clusters
308 TH2F * fhAsyMCGenRecoNLocMax2EbinPi0[4] ; //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when 2 local max, 4 E bins, neutral clusters
309 TH2F * fhAsyMCGenRecoNLocMaxNEbinPi0[4] ; //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when >2 local max, 4 E bins, neutral clusters
311 TH2F * fhAsyMCGenRecoDiffMCPi0[3]; //! reconstructed-generated asymmetry of splitted clusters vs E from pi0, for 3 NLM cases
312 TH2F * fhAsyMCGenRecoDiffMCPi0Conv[3]; //! reconstructed-generated asymmetry of splitted clusters vs E from converted pi0, for 3 NLM cases
314 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
315 TH2F * fhMassDispEtaNLocMax2[7][2] ; //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
316 TH2F * fhMassDispEtaNLocMaxN[7][2] ; //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
318 TH2F * fhMassDispEtaNLocMax1Ebin[4] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
319 TH2F * fhMassDispEtaNLocMax2Ebin[4] ; //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
320 TH2F * fhMassDispEtaNLocMaxNEbin[4] ; //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
322 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
323 TH2F * fhMassDispPhiNLocMax2[7][2] ; //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
324 TH2F * fhMassDispPhiNLocMaxN[7][2] ; //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
326 TH2F * fhMassDispPhiNLocMax1Ebin[4] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
327 TH2F * fhMassDispPhiNLocMax2Ebin[4] ; //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
328 TH2F * fhMassDispPhiNLocMaxNEbin[4] ; //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
330 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
331 TH2F * fhMassDispAsyNLocMax2[7][2] ; //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
332 TH2F * fhMassDispAsyNLocMaxN[7][2] ; //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
334 TH2F * fhMassDispAsyNLocMax1Ebin[4] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
335 TH2F * fhMassDispAsyNLocMax2Ebin[4] ; //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
336 TH2F * fhMassDispAsyNLocMaxNEbin[4] ; //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
338 TH2F * fhNLocMax [7][2] ; //! Number of maxima in cluster vs E, 1-6 for different MC particle types
339 TH2F * fhNLocMaxM02Cut[7][2] ; //! Number of maxima in cluster vs E, 1-6 for different MC particle types, after SS cut
340 TH2F * fhNLocMaxIdPi0 [7][2] ; //! Number of maxima in cluster vs E, 1-6 for different MC particle types, after pi0 selection
342 TH2F * fhSplitClusterENLocMax[7][2] ; //! Number of maxima in cluster vs E of splitted clusters, 1-6 for different MC particle types
343 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
345 TH2F * fhLM1NLocMax [7][2] ; //! Split cluster 1 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types
346 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
347 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
349 TH2F * fhLM2NLocMax [7][2] ; //! Split cluster 2 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types
350 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
351 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
353 TH2F * fhM02NLocMax1 [7][2] ; //! M02 vs E for N max in cluster = 1, 1-6 for different MC particle types
354 TH2F * fhM02NLocMax2 [7][2] ; //! M02 vs E for N max in cluster = 2, 1-6 for different MC particle types
355 TH2F * fhM02NLocMaxN [7][2] ; //! M02 vs E for N max in cluster > 2, 1-6 for different MC particle types
357 TH2F * fhMCAsymM02NLocMax1MCPi0Ebin[4] ; //! M02 vs decay asymmetry for N max in cluster = 1, for 4 energy bins
358 TH2F * fhMCAsymM02NLocMax2MCPi0Ebin[4] ; //! M02 vs decay asymmetry for N max in cluster = 2, for 4 energy bins
359 TH2F * fhMCAsymM02NLocMaxNMCPi0Ebin[4] ; //! M02 vs decay asymmetry for N max in cluster > 2, for 4 energy bins
361 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
362 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
363 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
365 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
366 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
367 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
369 TH2F * fhMCGenFracAfterCutsNLocMax1MCPi0 ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, MCPi0 after M02 and asymmetry cut
370 TH2F * fhMCGenFracAfterCutsNLocMax2MCPi0 ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, MCPi0, after M02 and asymmetry cut
371 TH2F * fhMCGenFracAfterCutsNLocMaxNMCPi0 ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, MCPi0, after M02 and asymmetry cut
373 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
374 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
375 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
377 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
378 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
379 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
381 TH2F * fhMCGenSplitEFracAfterCutsNLocMax1MCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types
382 TH2F * fhMCGenSplitEFracAfterCutsNLocMax2MCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types
383 TH2F * fhMCGenSplitEFracAfterCutsNLocMaxNMCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types
385 TH2F * fhMCGenEFracvsSplitEFracNLocMax1[7][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster = 1, MC pi0
386 TH2F * fhMCGenEFracvsSplitEFracNLocMax2[7][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster = 2, MC pi0
387 TH2F * fhMCGenEFracvsSplitEFracNLocMaxN[7][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster > 2, MC pi0
389 TH2F * fhMCGenEvsSplitENLocMax1[7][2] ; //! E generated particle vs E1+E2 for N max in cluster = 1, 1-6 for different MC particle types
390 TH2F * fhMCGenEvsSplitENLocMax2[7][2] ; //! E generated particle vs E1+E2 for N max in cluster = 2, 1-6 for different MC particle types
391 TH2F * fhMCGenEvsSplitENLocMaxN[7][2] ; //! E generated particle vs E1+E2 for N max in cluster > 2, 1-6 for different MC particle types
393 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
394 TH2F * fhMCGenFracNLocMaxEbinMatched[7][4] ; //! NLM vs E generated particle / E reconstructed vs E reconstructed 1-6 for different MC particle types, matched to track
396 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
397 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
398 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
400 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
401 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
402 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
404 TH2F * fhNCellNLocMax1[7][2] ; //! n cells in cluster vs E for N max in cluster = 1, 1-6 for different MC particle types
405 TH2F * fhNCellNLocMax2[7][2] ; //! n cells in cluster vs E for N max in cluster = 2, 1-6 for different MC particle types
406 TH2F * fhNCellNLocMaxN[7][2] ; //! n cells in cluster vs E for N max in cluster > 2, 1-6 for different MC particle types
408 TH2F * fhNCellMassEHighNLocMax1MCPi0 ; //! n cells in cluster vs mass for high energy clusters, for N max in cluster = 1, for MC pi0
409 TH2F * fhNCellM02EHighNLocMax1MCPi0 ; //! n cells in cluster vs m02 for high energy clusters, for N max in cluster = 1, for MC pi0
410 TH2F * fhNCellMassELowNLocMax1MCPi0 ; //! n cells in cluster vs mass for low energy clusters, for N max in cluster = 1, for MC pi0
411 TH2F * fhNCellM02ELowNLocMax1MCPi0 ; //! n cells in cluster vs m02 for low energy clusters, for N max in cluster = 1, for MC pi0
413 TH2F * fhNCellMassEHighNLocMax2MCPi0 ; //! n cells in cluster vs mass for high energy clusters, for N max in cluster = 2, for MC pi0
414 TH2F * fhNCellM02EHighNLocMax2MCPi0 ; //! n cells in cluster vs m02 for high energy clusters, for N max in cluster = 2, for MC pi0
415 TH2F * fhNCellMassELowNLocMax2MCPi0 ; //! n cells in cluster vs mass for low energy clusters, for N max in cluster = 2, for MC pi0
416 TH2F * fhNCellM02ELowNLocMax2MCPi0 ; //! n cells in cluster vs m02 for low energy clusters, for N max in cluster = 2, for MC pi0
418 TH2F * fhNCellMassEHighNLocMaxNMCPi0 ; //! n cells in cluster vs mass for high energy clusters, for N max in cluster > 2, for MC pi0
419 TH2F * fhNCellM02EHighNLocMaxNMCPi0 ; //! n cells in cluster vs m02 for high energy clusters, for N max in cluster > 2, for MC pi0
420 TH2F * fhNCellMassELowNLocMaxNMCPi0 ; //! n cells in cluster vs mass for low energy clusters, for N max in cluster > 2, for MC pi0
421 TH2F * fhNCellM02ELowNLocMaxNMCPi0 ; //! n cells in cluster vs m02 for low energy clusters, for N max in cluster > 2, for MC pi0
423 TH2F * fhM02Pi0NLocMax1[7][2] ; //! M02 for Mass around pi0, N Local Maxima = 1
424 TH2F * fhM02EtaNLocMax1[7][2] ; //! M02 for Mass around eta, N Local Maxima = 1
425 TH2F * fhM02ConNLocMax1[7][2] ; //! M02 for Mass around close to 0, N Local Maxima = 1
427 TH2F * fhM02Pi0NLocMax2[7][2] ; //! M02 for Mass around pi0, N Local Maxima = 2
428 TH2F * fhM02EtaNLocMax2[7][2] ; //! M02 for Mass around eta, N Local Maxima = 2
429 TH2F * fhM02ConNLocMax2[7][2] ; //! M02 for Mass around close to 0, N Local Maxima = 2
431 TH2F * fhM02Pi0NLocMaxN[7][2] ; //! M02 for Mass around pi0, N Local Maxima > 2
432 TH2F * fhM02EtaNLocMaxN[7][2] ; //! M02 for Mass around eta, N Local Maxima > 2
433 TH2F * fhM02ConNLocMaxN[7][2] ; //! M02 for Mass around close to 0, N Local Maxima > 2
435 TH2F * fhMassPi0NLocMax1[7][2] ; //! Mass for selected pi0, N Local Maxima = 1
436 TH2F * fhMassEtaNLocMax1[7][2] ; //! Mass for selected around eta, N Local Maxima = 1
437 TH2F * fhMassConNLocMax1[7][2] ; //! Mass for selected around close to 0, N Local Maxima = 1
439 TH2F * fhMassPi0NLocMax2[7][2] ; //! Mass for selected around pi0, N Local Maxima = 2
440 TH2F * fhMassEtaNLocMax2[7][2] ; //! Mass for selected around eta, N Local Maxima = 2
441 TH2F * fhMassConNLocMax2[7][2] ; //! Mass for selected around close to 0, N Local Maxima = 2
443 TH2F * fhMassPi0NLocMaxN[7][2] ; //! Mass for selected around pi0, N Local Maxima > 2
444 TH2F * fhMassEtaNLocMaxN[7][2] ; //! Mass for selected around eta, N Local Maxima > 2
445 TH2F * fhMassConNLocMaxN[7][2] ; //! Mass for selected around close to 0, N Local Maxima > 2
447 TH2F * fhNCellPi0NLocMax1[7][2] ; //! n cells for selected around pi0, N Local Maxima = 1
448 TH2F * fhNCellEtaNLocMax1[7][2] ; //! n cells for selected around eta, N Local Maxima = 1
449 TH2F * fhNCellPi0NLocMax2[7][2] ; //! n cells for selected around pi0, N Local Maxima = 2
450 TH2F * fhNCellEtaNLocMax2[7][2] ; //! n cells for selected around eta, N Local Maxima = 2
451 TH2F * fhNCellPi0NLocMaxN[7][2] ; //! n cells for selected around pi0, N Local Maxima > 2
452 TH2F * fhNCellEtaNLocMaxN[7][2] ; //! n cells for selected around eta, N Local Maxima > 2
454 TH2F * fhMassAfterCutsNLocMax1[7][2] ; //! Mass after M02, asymmetry cuts for MC part, N Local Maxima = 1
455 TH2F * fhMassAfterCutsNLocMax2[7][2] ; //! Mass after M02, asymmetry cuts for MC part, N Local Maxima = 2
456 TH2F * fhMassAfterCutsNLocMaxN[7][2] ; //! Mass after M02, asymmetry cuts for MC part, N Local Maxima > 2
458 TH2F * fhMassSplitEAfterCutsNLocMax1[7][2] ; //! Split Inv Mass vs E1+E2, NLM=1, after M02, asymmetry cuts, different MC particle types, track matching on/off
459 TH2F * fhMassSplitEAfterCutsNLocMax2[7][2] ; //! Split Inv Mass vs E1+E2, NLM=2, after M02, asymmetry cuts, different MC particle types, track matching on/off
460 TH2F * fhMassSplitEAfterCutsNLocMaxN[7][2] ; //! Split Inv Mass vs E1+E2, NLM>2, after M02, asymmetry cuts, different MC particle types, track matching on/off
462 TH2F * fhMassSplitEPi0NLocMax1[7][2] ; //! Split Inv Mass vs E1+E2, NLM=1, after pi0 selection, different MC particle types, track matching on/off
463 TH2F * fhMassSplitEPi0NLocMax2[7][2] ; //! Split Inv Mass vs E1+E2, NLM=2, after pi0 selection, different MC particle types, track matching on/off
464 TH2F * fhMassSplitEPi0NLocMaxN[7][2] ; //! Split Inv Mass vs E1+E2, NLM>2, after pi0 selection, different MC particle types, track matching on/off
466 TH2F * fhAsyPi0NLocMax1[7][2] ; //! Asy for Mass around pi0, N Local Maxima = 1
467 TH2F * fhAsyEtaNLocMax1[7][2] ; //! Asy for Mass around eta, N Local Maxima = 1
468 TH2F * fhAsyConNLocMax1[7][2] ; //! Asy for Mass around close to 0, N Local Maxima = 1
470 TH2F * fhAsyPi0NLocMax2[7][2] ; //! Asy for Mass around pi0, N Local Maxima = 2
471 TH2F * fhAsyEtaNLocMax2[7][2] ; //! Asy for Mass around eta, N Local Maxima = 2
472 TH2F * fhAsyConNLocMax2[7][2] ; //! Asy for Mass around close to 0, N Local Maxima = 2
474 TH2F * fhAsyPi0NLocMaxN[7][2] ; //! Asy for Mass around pi0, N Local Maxima > 2
475 TH2F * fhAsyEtaNLocMaxN[7][2] ; //! Asy for Mass around eta, N Local Maxima > 2
476 TH2F * fhAsyConNLocMaxN[7][2] ; //! Asy for Mass around close to 0, N Local Maxima > 2
478 TH2F * fhSplitEFractionNLocMax1[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1
479 TH2F * fhSplitEFractionNLocMax2[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2
480 TH2F * fhSplitEFractionNLocMaxN[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2
482 TH2F * fhSplitEFractionAfterCutsNLocMax1[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1, after M02 and asymmetry cut
483 TH2F * fhSplitEFractionAfterCutsNLocMax2[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2, after M02 and asymmetry cut
484 TH2F * fhSplitEFractionAfterCutsNLocMaxN[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2, after M02 and asymmetry cut
486 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
487 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
488 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
490 TH2F * fhAnglePairNLocMax1[7][2] ; //! pair opening angle vs E
491 TH2F * fhAnglePairNLocMax2[7][2] ; //! pair opening angle vs E
492 TH2F * fhAnglePairNLocMaxN[7][2] ; //! pair opening angle vs E
494 TH2F * fhAnglePairAfterCutsNLocMax1[7][2] ; //! pair opening angle vs E
495 TH2F * fhAnglePairAfterCutsNLocMax2[7][2] ; //! pair opening angle vs E
496 TH2F * fhAnglePairAfterCutsNLocMaxN[7][2] ; //! pair opening angle vs E
498 TH2F * fhAnglePairPi0NLocMax1[7][2] ; //! pair opening angle vs E
499 TH2F * fhAnglePairPi0NLocMax2[7][2] ; //! pair opening angle vs E
500 TH2F * fhAnglePairPi0NLocMaxN[7][2] ; //! pair opening angle vs E
502 TH2F * fhAnglePairMassNLocMax1[7][2] ; //! pair opening angle vs Mass for E > 7 GeV
503 TH2F * fhAnglePairMassNLocMax2[7][2] ; //! pair opening angle vs Mass for E > 7 GeV
504 TH2F * fhAnglePairMassNLocMaxN[7][2] ; //! pair opening angle vs Mass for E > 7 GeV
506 TH2F * fhAnglePairM02NLocMax1[7][2] ; //! pair opening angle vs M02 for E > 7 GeV
507 TH2F * fhAnglePairM02NLocMax2[7][2] ; //! pair opening angle vs M02 for E > 7 GeV
508 TH2F * fhAnglePairM02NLocMaxN[7][2] ; //! pair opening angle vs M02 for E > 7 GeV
510 TH2F * fhAnglePairPrimPi0RecoNLocMax1; //! pair opening angle pi0 generated/reconstructed vs E
511 TH2F * fhAnglePairPrimPi0RecoNLocMax2; //! pair opening angle pi0 generated/reconstructed vs E
512 TH2F * fhAnglePairPrimPi0RecoNLocMaxN; //! pair opening angle pi0 generated/reconstructed vs E
514 TH2F * fhAnglePairPrimPi0vsRecoNLocMax1; //! pair opening angle pi0 generated vs reconstructed
515 TH2F * fhAnglePairPrimPi0vsRecoNLocMax2; //! pair opening angle pi0 generated vs reconstructed
516 TH2F * fhAnglePairPrimPi0vsRecoNLocMaxN; //! pair opening angle pi0 generated vs reconstructed
518 TH2F * fhAnglePairOverM02NLocMax1[7][2]; //! pair opening angle / m02 vs E, NLM=1
519 TH2F * fhAnglePairOverM02NLocMax2[7][2]; //! pair opening angle / m02 vs E, NLM=2
520 TH2F * fhAnglePairOverM02NLocMaxN[7][2]; //! pair opening angle / m02 vs E, NLM=N
522 TH2F * fhAnglePairOverM02NLocMax1Overlap0[7][2];//! pair opening angle / m02 vs E, NLM=1
523 TH2F * fhAnglePairOverM02NLocMax2Overlap0[7][2];//! pair opening angle / m02 vs E, NLM=2
524 TH2F * fhAnglePairOverM02NLocMaxNOverlap0[7][2];//! pair opening angle / m02 vs E, NLM=N
526 TH2F * fhAnglePairPrimPi0OverM02NLocMax1; //! pair opening angle / m02 vs E, NLM=1, prim pi0
527 TH2F * fhAnglePairPrimPi0OverM02NLocMax2; //! pair opening angle / m02 vs E, NLM=2, prim pi0
528 TH2F * fhAnglePairPrimPi0OverM02NLocMaxN; //! pair opening angle / m02 vs E, NLM=N, prim pi0
530 TH2F * fhArmNLocMax1[7][4] ; //! Armenteros of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types
531 TH2F * fhArmNLocMax2[7][4] ; //! Armenteros of 2 cells local maxima vs E, 1-6 for different MC particle types
532 TH2F * fhArmNLocMaxN[7][4] ; //! Armenteros of >2 cells local maxima vs E, 1-6 for different MC particle types
534 TH2F * fhArmAfterCutsNLocMax1[7][4] ; //! Armenteros after M02, asymmetry cuts for pi0, N Local Maxima = 1
535 TH2F * fhArmAfterCutsNLocMax2[7][4] ; //! Armenteros after M02, asymmetry cuts for pi0, N Local Maxima = 2
536 TH2F * fhArmAfterCutsNLocMaxN[7][4] ; //! Armenteros after M02, asymmetry cuts for pi0, N Local Maxima > 2
538 TH2F * fhArmPi0NLocMax1[7][4] ; //! Armenteros for selected pi0, N Local Maxima = 1
539 TH2F * fhArmPi0NLocMax2[7][4] ; //! Armenteros for selected pi0, N Local Maxima = 2
540 TH2F * fhArmPi0NLocMaxN[7][4] ; //! Armenteros for selected pi0, N Local Maxima > 2
542 TH2F * fhCosThStarNLocMax1[7][2] ; //! cos(theta^star) vs E, NLM=1
543 TH2F * fhCosThStarNLocMax2[7][2] ; //! cos(theta^star) vs E, NLM=2
544 TH2F * fhCosThStarNLocMaxN[7][2] ; //! cos(theta^star) vs E, NLM>2
546 TH2F * fhCosThStarAfterCutsNLocMax1[7][2] ; //! cos(theta^star) vs E, after M02, asymmetry cuts, NLM=1
547 TH2F * fhCosThStarAfterCutsNLocMax2[7][2] ; //! cos(theta^star) vs E, after M02, asymmetry cuts, NLM=2
548 TH2F * fhCosThStarAfterCutsNLocMaxN[7][2] ; //! cos(theta^star) vs E, after M02, asymmetry cuts, NLM>2
550 TH2F * fhCosThStarPi0NLocMax1[7][2] ; //! cos(theta^star) vs E, after M02, asymmetry and pi0 mass cuts, NLM=1
551 TH2F * fhCosThStarPi0NLocMax2[7][2] ; //! cos(theta^star) vs E, after M02, asymmetry and pi0 mass cuts, NLM=2
552 TH2F * fhCosThStarPi0NLocMaxN[7][2] ; //! cos(theta^star) vs E, after M02, asymmetry and pi0 mass cuts, NLM>2
554 TH2F * fhTrackMatchedDEtaNLocMax1[7] ; //! Eta distance between track and cluster vs cluster E, 1 local maximum
555 TH2F * fhTrackMatchedDPhiNLocMax1[7] ; //! Phi distance between track and cluster vs cluster E, 1 local maximum
556 TH2F * fhTrackMatchedDEtaNLocMax2[7] ; //! Eta distance between track and cluster vs cluster E, 2 local maximum
557 TH2F * fhTrackMatchedDPhiNLocMax2[7] ; //! Phi distance between track and cluster vs cluster E, 2 local maximum
558 TH2F * fhTrackMatchedDEtaNLocMaxN[7] ; //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
559 TH2F * fhTrackMatchedDPhiNLocMaxN[7] ; //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
561 TH2F * fhTrackMatchedDEtaNLocMax1Pos[7] ; //! Eta distance between track and cluster vs cluster E, 1 local maximum
562 TH2F * fhTrackMatchedDPhiNLocMax1Pos[7] ; //! Phi distance between track and cluster vs cluster E, 1 local maximum
563 TH2F * fhTrackMatchedDEtaNLocMax2Pos[7] ; //! Eta distance between track and cluster vs cluster E, 2 local maximum
564 TH2F * fhTrackMatchedDPhiNLocMax2Pos[7] ; //! Phi distance between track and cluster vs cluster E, 2 local maximum
565 TH2F * fhTrackMatchedDEtaNLocMaxNPos[7] ; //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
566 TH2F * fhTrackMatchedDPhiNLocMaxNPos[7] ; //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
568 TH2F * fhTrackMatchedDEtaNLocMax1Neg[7] ; //! Eta distance between track and cluster vs cluster E, 1 local maximum
569 TH2F * fhTrackMatchedDPhiNLocMax1Neg[7] ; //! Phi distance between track and cluster vs cluster E, 1 local maximum
570 TH2F * fhTrackMatchedDEtaNLocMax2Neg[7] ; //! Eta distance between track and cluster vs cluster E, 2 local maximum
571 TH2F * fhTrackMatchedDPhiNLocMax2Neg[7] ; //! Phi distance between track and cluster vs cluster E, 2 local maximum
572 TH2F * fhTrackMatchedDEtaNLocMaxNNeg[7] ; //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
573 TH2F * fhTrackMatchedDPhiNLocMaxNNeg[7] ; //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
575 TH2F * fhCentralityPi0NLocMax1 ; //! Centrality for selected pi0, N Local Maxima = 1
576 TH2F * fhCentralityEtaNLocMax1 ; //! Centrality for selected eta, N Local Maxima = 1
577 TH2F * fhCentralityPi0NLocMax2 ; //! Centrality for selected pi0, N Local Maxima = 2
578 TH2F * fhCentralityEtaNLocMax2 ; //! Centrality for selected eta, N Local Maxima = 2
579 TH2F * fhCentralityPi0NLocMaxN ; //! Centrality for selected pi0, N Local Maxima > 2
580 TH2F * fhCentralityEtaNLocMaxN ; //! Centrality for selected eta, N Local Maxima > 2
582 TH2F * fhEventPlanePi0NLocMax1 ; //! Event plane for selected pi0, N Local Maxima = 1
583 TH2F * fhEventPlaneEtaNLocMax1 ; //! Event plane for selected eta, N Local Maxima = 1
584 TH2F * fhEventPlanePi0NLocMax2 ; //! Event plane for selected pi0, N Local Maxima = 2
585 TH2F * fhEventPlaneEtaNLocMax2 ; //! Event plane for selected eta, N Local Maxima = 2
586 TH2F * fhEventPlanePi0NLocMaxN ; //! Event plane for selected pi0, N Local Maxima > 2
587 TH2F * fhEventPlaneEtaNLocMaxN ; //! Event plane for selected eta, N Local Maxima > 2
589 TH2F * fhClusterEtaPhiNLocMax1 ; //! Eta vs Phi of clusters with N Local Maxima = 1, E > 8 GeV
590 TH2F * fhClusterEtaPhiNLocMax2 ; //! Eta vs Phi of clusters with N Local Maxima = 2, E > 8 GeV
591 TH2F * fhClusterEtaPhiNLocMaxN ; //! Eta vs Phi of clusters with N Local Maxima > 2, E > 8 GeV
592 TH2F * fhPi0EtaPhiNLocMax1 ; //! Eta vs Phi of pi0's with N Local Maxima = 1, E > 8 GeV
593 TH2F * fhPi0EtaPhiNLocMax2 ; //! Eta vs Phi of pi0's with N Local Maxima = 2, E > 8 GeV
594 TH2F * fhPi0EtaPhiNLocMaxN ; //! Eta vs Phi of pi0's with N Local Maxima > N, E > 8 GeV
595 TH2F * fhEtaEtaPhiNLocMax1 ; //! Eta vs Phi of eta's with N Local Maxima = 1, E > 8 GeV
596 TH2F * fhEtaEtaPhiNLocMax2 ; //! Eta vs Phi of eta's with N Local Maxima = 2, E > 8 GeV
597 TH2F * fhEtaEtaPhiNLocMaxN ; //! Eta vs Phi of eta's with N Local Maxima > N, E > 8 GeV
599 TH2F * fhPi0CellE[3] ; //! pi0's energy vs cluster cell energy with NLM = 1, = 2, > 2
600 TH2F * fhPi0CellEFrac[3] ; //! pi0's energy vs cluster cell energy fraction with NLM = 1, = 2, > 2
601 TH2F * fhPi0CellLogEFrac[3] ; //! pi0's energy vs cluster log cell energy fraction with NLM = 1, = 2, > 2
602 TH2F * fhPi0CellEMaxEMax2Frac [3]; //! pi0's energy vs fraction of 2 main maxima energy with NLM = 1, = 2, > 2
603 TH2F * fhPi0CellEMaxClusterFrac [3]; //! pi0's energy vs energy fraction of main LM and cluster energy with NLM = 1, = 2, > 2
604 TH2F * fhPi0CellEMax2ClusterFrac[3]; //! pi0's energy vs energy fraction of second LM and cluster energy with NLM = 1, = 2, > 2
605 TH2F * fhPi0CellEMaxFrac [3]; //! pi0's energy vs energy fraction of main LM and cluster cell energy with NLM = 1, = 2, > 2
606 TH2F * fhPi0CellEMax2Frac [3]; //! pi0's energy vs energy fraction of second LM and cluster cell energy with NLM = 1, = 2, > 2
608 TH2F * fhM02WeightPi0[3][20] ; //! M02 for selected pi0 with different weight, with NLM = 1, = 2, > 2
609 TH2F * fhM02ECellCutPi0[3][20] ; //! M02 for selected pi0 with different cut on cell energy, with NLM = 1, = 2, > 2
611 TH2F * fhPi0EPairDiffTimeNLM1; //! E vs Pair of clusters time difference vs E, for selected pi0, NLM=1
612 TH2F * fhPi0EPairDiffTimeNLM2; //! E vs Pair of clusters time difference vs E, for selected pi0, NLM=2
613 TH2F * fhPi0EPairDiffTimeNLMN; //! E vs Pair of clusters time difference vs E, for selected pi0, NLM>2
614 TH2F * fhEtaEPairDiffTimeNLM1; //! E vs Pair of clusters time difference vs E, for selected eta, NLM=1
615 TH2F * fhEtaEPairDiffTimeNLM2; //! E vs Pair of clusters time difference vs E, for selected eta, NLM=2
616 TH2F * fhEtaEPairDiffTimeNLMN; //! E vs Pair of clusters time difference vs E, for selected eta, NLM>2
618 TH2F * fhMCEM02Overlap0[3][7]; //! E vs M02 for different MC origin, no other MC particles contributes, neutral cluster
619 TH2F * fhMCEM02Overlap1[3][7]; //! E vs M02 for different MC origin, 1 other MC particles contributes, neutral cluster
620 TH2F * fhMCEM02OverlapN[3][7]; //! E vs M02 for different MC origin, N other MC particles contributes, neutral cluster
621 TH2F * fhMCEM02Overlap0Match[3][7]; //! E vs M02 for different MC origin, no other MC particles contributes, charged cluster
622 TH2F * fhMCEM02Overlap1Match[3][7]; //! E vs M02 for different MC origin, 1 other MC particles contributes, charged cluster
623 TH2F * fhMCEM02OverlapNMatch[3][7]; //! E vs M02 for different MC origin, N other MC particles contributes, charged cluster
625 TH2F * fhMCEMassOverlap0[3][7]; //! E vs Mass for different MC origin, no other MC particles contributes, neutral cluster
626 TH2F * fhMCEMassOverlap1[3][7]; //! E vs Mass for different MC origin, 1 other MC particles contributes, neutral cluster
627 TH2F * fhMCEMassOverlapN[3][7]; //! E vs Mass for different MC origin, N other MC particles contributes, neutral cluster
628 TH2F * fhMCEMassOverlap0Match[3][7]; //! E vs Mass for different MC origin, no other MC particles contributes, charged cluster
629 TH2F * fhMCEMassOverlap1Match[3][7]; //! E vs Mass for different MC origin, 1 other MC particles contributes, charged cluster
630 TH2F * fhMCEMassOverlapNMatch[3][7]; //! E vs Mass for different MC origin, N other MC particles contributes, charged cluster
632 TH2F * fhMCESplitEFracOverlap0[3][7]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, no other MC particles contributes, neutral cluster
633 TH2F * fhMCESplitEFracOverlap1[3][7]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, 1 other MC particles contributes, neutral cluster
634 TH2F * fhMCESplitEFracOverlapN[3][7]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, N other MC particles contributes, neutral cluster
635 TH2F * fhMCESplitEFracOverlap0Match[3][7]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, no other MC particles contributes, charged cluster
636 TH2F * fhMCESplitEFracOverlap1Match[3][7]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, 1 other MC particles contributes, charged cluster
637 TH2F * fhMCESplitEFracOverlapNMatch[3][7]; //! E vs sum of splitted cluster energy / cluster energy for different MC origin, N other MC particles contributes, charged cluster
639 TH2F * fhMCEAsymOverlap0[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contributes, neutral cluster
640 TH2F * fhMCEAsymOverlap1[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, 1 other MC particles contributes, neutral cluster
641 TH2F * fhMCEAsymOverlapN[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, N other MC particles contributes, neutral cluster
642 TH2F * fhMCEAsymOverlap0Match[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contributes, charged cluster
643 TH2F * fhMCEAsymOverlap1Match[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, 1 other MC particles contributes, charged cluster
644 TH2F * fhMCEAsymOverlapNMatch[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, N other MC particles contributes, charged cluster
646 TH2F * fhMCENCellOverlap0[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contributes, neutral cluster
647 TH2F * fhMCENCellOverlap1[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, 1 other MC particles contributes, neutral cluster
648 TH2F * fhMCENCellOverlapN[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, N other MC particles contributes, neutral cluster
649 TH2F * fhMCENCellOverlap0Match[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contributes, charged cluster
650 TH2F * fhMCENCellOverlap1Match[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, 1 other MC particles contributes, charged cluster
651 TH2F * fhMCENCellOverlapNMatch[3][7]; //! E vs sum of splitted cluster energy asymmetry for different MC origin, N other MC particles contributes, charged cluster
653 TH2F * fhMCEEpriOverlap0[3][7]; //! E reco vs primary for different MC origin, no other MC particles contributes, neutral cluster
654 TH2F * fhMCEEpriOverlap1[3][7]; //! E reco vs primary for different MC origin, 1 other MC particles contributes, neutral cluster
655 TH2F * fhMCEEpriOverlapN[3][7]; //! E reco vs primary for different MC origin, N other MC particles contributes, neutral cluster
656 TH2F * fhMCEEpriOverlap0Match[3][7]; //! E reco vs primary for different MC origin, no other MC particles contributes, charged cluster
657 TH2F * fhMCEEpriOverlap1Match[3][7]; //! E reco vs primary for different MC origin, 1 other MC particles contributes, charged cluster
658 TH2F * fhMCEEpriOverlapNMatch[3][7]; //! E reco vs primary for different MC origin, N other MC particles contributes, charged cluster
660 TH2F * fhMCEEpriOverlap0IdPi0[3][7]; //! E reco vs primary for different MC origin, no other MC particles contributes, neutral cluster, neutral clusters id as pi0
661 TH2F * fhMCEEpriOverlap1IdPi0[3][7]; //! E reco vs primary for different MC origin, 1 other MC particles contributes, neutral cluster, neutral clusters id as pi0
662 TH2F * fhMCEEpriOverlapNIdPi0[3][7]; //! E reco vs primary for different MC origin, 1 other MC particles contributes, neutral cluster, neutral clusters is as pi0
664 TH2F * fhMCPi0MassM02Overlap0[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, no other MC particles contributes, neutral cluster, 4 E bins
665 TH2F * fhMCPi0MassM02Overlap1[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, 1 other MC particles contributes, neutral cluster, 4 E bins
666 TH2F * fhMCPi0MassM02OverlapN[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, N other MC particles contributes, neutral cluster, 4 E bins
667 TH2F * fhMCPi0MassM02Overlap0Match[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, no other MC particles contributes, charged cluster, 4 E bins
668 TH2F * fhMCPi0MassM02Overlap1Match[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, 1 other MC particles contributes, charged cluster, 4 E bins
669 TH2F * fhMCPi0MassM02OverlapNMatch[3][4]; //! MC Pi0 M02 vs Mass for different MC origin, N other MC particles contributes, charged cluster, 4 E bins
671 TH2F * fhMCENOverlaps[3][7]; //! E vs number of Overlaps in MC, neutral cluster
672 TH2F * fhMCENOverlapsMatch[3][7]; //! E vs number of Overlaps in MC, charged cluster
674 TH2F * fhMCPi0HighNLMPair; //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons
675 TH2F * fhMCPi0LowNLMPair; //! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy correspond to 2 photons
676 TH2F * fhMCPi0AnyNLMPair; //! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs correspond to 2 photons
677 TH2F * fhMCPi0NoneNLMPair; //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons
678 // No match between highest energy local maxima and highest energy MC particle
679 TH2F * fhMCPi0HighNLMPairNoMCMatch; //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons
680 TH2F * fhMCPi0LowNLMPairNoMCMatch; //! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy correspond to 2 photons
681 TH2F * fhMCPi0AnyNLMPairNoMCMatch; //! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs correspond to 2 photons
682 TH2F * fhMCPi0NoneNLMPairNoMCMatch; //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons
684 TH2F * fhMCPi0HighNLMPairOverlap; //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons, overlap
685 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
686 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
687 TH2F * fhMCPi0NoneNLMPairOverlap; //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons, overlap
688 // No match between highest energy local maxima and highest energy MC particle
689 TH2F * fhMCPi0HighNLMPairNoMCMatchOverlap; //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons, overlap
690 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
691 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
692 TH2F * fhMCPi0NoneNLMPairNoMCMatchOverlap; //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons, overlap
694 TH2F * fhMCPi0DecayPhotonHitHighLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
695 TH2F * fhMCPi0DecayPhotonAdjHighLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
696 TH2F * fhMCPi0DecayPhotonHitOtherLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
697 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
698 TH2F * fhMCPi0DecayPhotonAdjacent; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM
699 TH2F * fhMCPi0DecayPhotonHitNoLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas
701 TH2F * fhMCPi0DecayPhotonHitHighLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maxima, overlap
702 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
703 TH2F * fhMCPi0DecayPhotonHitOtherLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
704 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
705 TH2F * fhMCPi0DecayPhotonAdjacentOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM, overlap
706 TH2F * fhMCPi0DecayPhotonHitNoLMOverlap; //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas, overlap
708 TH2F * fhMCPi0DecayPhotonHitHighLMDiffELM1[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
709 TH2F * fhMCPi0DecayPhotonAdjHighLMDiffELM1[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
710 TH2F * fhMCPi0DecayPhotonHitOtherLMDiffELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
711 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
713 TH2F * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM1[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
714 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
715 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
716 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
718 TH2F * fhMCPi0DecayPhotonHitHighLMDiffELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
719 TH2F * fhMCPi0DecayPhotonAdjHighLMDiffELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
720 TH2F * fhMCPi0DecayPhotonHitOtherLMDiffELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
721 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
723 TH2F * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM2[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
724 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
725 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
726 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
729 TH2F * fhMCPi0DecayPhotonHitHighLMDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
730 TH2F * fhMCPi0DecayPhotonAdjHighLMDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
731 TH2F * fhMCPi0DecayPhotonHitOtherLMDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
732 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
734 TH2F * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
735 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM1vsELM1[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
736 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
737 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
739 TH2F * fhMCPi0DecayPhotonHitHighLMDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
740 TH2F * fhMCPi0DecayPhotonAdjHighLMDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
741 TH2F * fhMCPi0DecayPhotonHitOtherLMDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
742 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
744 TH2F * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
745 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM2vsELM2[3]; //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
746 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
747 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
750 TH2F * fhMCPi0DecayPhotonHitHighLMMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
751 TH2F * fhMCPi0DecayPhotonAdjHighLMMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
752 TH2F * fhMCPi0DecayPhotonHitOtherLMMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
753 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
754 TH2F * fhMCPi0DecayPhotonAdjacentMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM
755 TH2F * fhMCPi0DecayPhotonHitNoLMMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas
757 TH2F * fhMCPi0DecayPhotonHitHighLMOverlapMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maxima, overlap
758 TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
759 TH2F * fhMCPi0DecayPhotonHitOtherLMOverlapMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
760 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
761 TH2F * fhMCPi0DecayPhotonAdjacentOverlapMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM, overlap
762 TH2F * fhMCPi0DecayPhotonHitNoLMOverlapMass[3]; //! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas, overlap
765 TH2F * fhMCEOverlapType; //! what particles overlap with pi0, neutral clusters
766 TH2F * fhMCEOverlapTypeMatch; //! what particles overlap with pi0, charged clusters
768 TH2F * fhMassBadDistClose[3]; //! split mass of clusters with second LM close to bad channel
769 TH2F * fhM02BadDistClose[3]; //! m02 of clusters with second LM close to bad channel
770 TH2F * fhMassOnBorder[3]; //! split mass of clusters with second LM on EMCAL border
771 TH2F * fhM02OnBorder[3]; //! m02 of clusters with second LM close to EMCAL border
774 TH2F * fhNLocMaxDiffCut [5][5] [2] ; //! Number of maxima for different values of min Loc Max value and min difference between cells, matched/unmatched with tracks
775 TH2F * fhM02NLocMaxDiffCut[5][5][3][2] ; //! M02 for 3 kinds of number of maxima for different values of min Loc Max value and min difference between cells, matched/unmatched with tracks
776 TH2F * fhMassNLocMaxDiffCut[5][5][3][2] ; //! Mass for 3 kinds of number of maxima for different values of min Loc Max value and min difference between cells, matched/unmatched with tracks
778 TH2F * fhNLocMaxDiffCutPi0 [5][5] [2] ; //! Number of maxima for different values of min Loc Max value and min difference between cells, matched/unmatched with tracks, cluster selected as pi0
779 TH2F * fhM02NLocMaxDiffCutPi0[5][5][3][2] ; //! M02 for 3 kinds of number of maxima for different values of min Loc Max value and min difference between cells, matched/unmatched with tracks, cluster selected as pi0
780 TH2F * fhMassNLocMaxDiffCutPi0[5][5][3][2] ; //! M02 for 3 kinds of number of maxima for different values of min Loc Max value and min difference between cells, matched/unmatched with tracks
782 AliAnaInsideClusterInvariantMass( const AliAnaInsideClusterInvariantMass & split) ; // cpy ctor
783 AliAnaInsideClusterInvariantMass & operator = (const AliAnaInsideClusterInvariantMass & split) ; // cpy assignment
785 ClassDef(AliAnaInsideClusterInvariantMass,30)
789 #endif //ALIANAINSIDECLUSTERINVARIANTMASS_H