]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.h
add option to study variation of NLM depedending on cuts
[u/mrichter/AliRoot.git] / PWGGA / CaloTrackCorrelations / AliAnaInsideClusterInvariantMass.h
CommitLineData
992b14a7 1#ifndef ALIANAINSIDECLUSTERINVARIANTMASS_H
2#define ALIANAINSIDECLUSTERINVARIANTMASS_H
3/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
992b14a7 5
6//_________________________________________________________________________
7//
8// Split clusters with some criteria and calculate invariant mass
9// to identify them as pi0 or conversion
10//
11//
83351853 12//-- Author: Gustavo Conesa (LPSC-Grenoble)
992b14a7 13//_________________________________________________________________________
14
15
16// --- ROOT system ---
17class TList ;
18class TObjString;
19class TLorentzVector;
20
21// --- ANALYSIS system ---
5c46c992 22class AliAODCaloCluster;
23
745913ae 24#include "AliAnaCaloTrackCorrBaseClass.h"
992b14a7 25
745913ae 26class AliAnaInsideClusterInvariantMass : public AliAnaCaloTrackCorrBaseClass {
992b14a7 27
28 public:
3c1d9afb 29
992b14a7 30 AliAnaInsideClusterInvariantMass() ; // default ctor
31 virtual ~AliAnaInsideClusterInvariantMass() { ; } //virtual dtor
992b14a7 32
b94e038e 33 void CheckLocalMaximaMCOrigin(AliVCluster* cluster, Int_t mcindex, Int_t noverlaps,
34 Float_t e1, Float_t e2, Float_t mass);
b583134f 35 //, Float_t m02, TLorentzVector l1, TLorentzVector l2);
b2e375c7 36
992b14a7 37 TObjString * GetAnalysisCuts();
38
39 TList * GetCreateOutputObjects();
dbe09c26 40
4914e781 41 void GetMCIndex(AliVCluster * cluster, Int_t & mcindex, Int_t & mcTag);
b2e375c7 42
b94e038e 43 void GetMCPrimaryKine(AliVCluster* cluster, Int_t mcindex, Int_t mcTag, Bool_t matched,
0186b6a2 44 Float_t & eprim, Float_t & asymGen, Float_t & angleGen, Int_t & noverlaps );
b2e375c7 45
b94e038e 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);
b2e375c7 50
b94e038e 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);
cb99fa55 54
b94e038e 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);
29555e96 58
b94e038e 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);
61
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);
68
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);
74
75 void FillSSWeightHistograms(AliVCluster *cluster, Int_t nlm, Int_t absId1, Int_t absId2);
76
77 void FillSSExtraHistograms(AliVCluster *cluster, Int_t nMax,
78 Bool_t matched, Int_t mcindex,
79 Float_t mass , Int_t ebin);
4914e781 80
ff6aa4ce 81 void FillNLMDiffCutHistograms(AliVCluster *cluster, AliVCaloCells *cells, Bool_t matched);
82
b94e038e 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);
b2e375c7 86
b94e038e 87 void FillTrackMatchingHistograms(AliVCluster * cluster,Int_t nMax, Int_t mcindex);
dbe09c26 88
b94e038e 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);
1253480f 93
94
b94e038e 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);
99
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);
105
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);
111
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);
1253480f 115
992b14a7 116 void Init();
117
118 void InitParameters();
1253480f 119
120 void MakeAnalysisFillHistograms() ;
121
71e3889f 122 void Print(const Option_t * opt) const;
992b14a7 123
71e3889f 124 void SetCalorimeter(TString & det) { fCalorimeter = det ; }
1253480f 125
71e3889f 126 void SetMinNCells(Int_t cut) { fMinNCells = cut ; }
2cb134fb 127
128 void SetMinBadChannelDistance(Float_t cut) { fMinBadDist = cut ; }
992b14a7 129
8edbd100 130 void SetWCorrectionParameter(Float_t p = 0.07) { fWSimu = p ; }
131
883411b2 132 void SwitchOnFillAngleHistograms() { fFillAngleHisto = kTRUE ; }
133 void SwitchOffFillAngleHistograms() { fFillAngleHisto = kFALSE ; }
cb99fa55 134
135 void SwitchOnFillArmenterosHistograms() { fFillArmenterosHisto = kTRUE ; }
136 void SwitchOffFillArmenterosHistograms() { fFillArmenterosHisto = kFALSE ; }
137
29555e96 138 void SwitchOnFillThetaStarHistograms() { fFillThetaStarHisto = kTRUE ; }
139 void SwitchOffFillThetaStarHistograms() { fFillThetaStarHisto = kFALSE ; }
140
883411b2 141 void SwitchOnFillExtraSSHistograms() { fFillSSExtraHisto = kTRUE ; }
142 void SwitchOffFillExtraSSHistograms() { fFillSSExtraHisto = kFALSE ; }
1253480f 143
144 void SwitchOnFillHighMultHistograms() { fFillHighMultHisto = kTRUE ; }
145 void SwitchOffFillHighMultHistograms() { fFillHighMultHisto = kFALSE ; }
146
147 void SwitchOnFillIdConvHistograms() { fFillIdConvHisto = kTRUE ; }
148 void SwitchOffFillIdConvHistograms() { fFillIdConvHisto = kFALSE ; }
8e81c2cf 149
1253480f 150 void SwitchOnFillIdEtaHistograms() { fFillIdEtaHisto = kTRUE ; }
151 void SwitchOffFillIdEtaHistograms() { fFillIdEtaHisto = kFALSE ; }
152
8edbd100 153 void SwitchOnFillTMHistograms() { fFillTMHisto = kTRUE ; }
154 void SwitchOffFillTMHistograms() { fFillTMHisto = kFALSE ; }
155
8e81c2cf 156 void SwitchOnFillTMResidualHistograms() { fFillTMResidualHisto = kTRUE ; }
157 void SwitchOffFillTMResidualHistograms() { fFillTMResidualHisto = kFALSE ; }
158
b2e375c7 159 void SwitchOnFillMCPrimaryHistograms() { fFillMCHisto = kTRUE ; }
160 void SwitchOffFillMCPrimaryHistograms() { fFillMCHisto = kFALSE ; }
2a77f6f4 161
19391b8c 162 void SwitchOnFillSSWeightHistograms() { fFillSSWeightHisto = kTRUE ; }
163 void SwitchOffFillSSWeightHistograms() { fFillSSWeightHisto = kFALSE ; }
164
ff6aa4ce 165 void SwitchOnFillNLMDiffCutsHistograms() { fFillNLMDiffCutHisto = kTRUE ; }
166 void SwitchOffFillNLMDiffCutsHistograms() { fFillNLMDiffCutHisto = kFALSE ; }
167
19391b8c 168 void SwitchOnFillEbinHistograms() { fFillEbinHisto = kTRUE ; }
169 void SwitchOffFillEbinHistograms() { fFillEbinHisto = kFALSE ; }
dbe09c26 170
b2e375c7 171 void SwitchOnFillMCOverlapHistograms() { fFillMCOverlapHisto = kTRUE ; }
172 void SwitchOffFillMCOverlapHistograms() { fFillMCOverlapHisto = kFALSE ; }
4914e781 173
174 void SwitchOnFillNCellHistograms() { fFillNCellHisto = kTRUE ; }
175 void SwitchOffFillNCellHistograms() { fFillNCellHisto = kFALSE ; }
b2e375c7 176
1253480f 177 void SwitchOnSplitClusterDistToBad() { fCheckSplitDistToBad = kTRUE ; }
178 void SwitchOffSplitClusterDistToBad() { fCheckSplitDistToBad = kFALSE ; }
179
ff6aa4ce 180 void SetNWeightForShowerShape(Int_t n) { fSSWeightN = n ; }
181 void SetWeightForShowerShape(Int_t i, Float_t v)
182 { if (i < 10) fSSWeight[i] = v ; }
dbe09c26 183
ff6aa4ce 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 ; }
187
188
189 void SetNECellCutForShowerShape(Int_t n) { fSSECellCutN = n ; }
190 void SetECellCutForShowerShape(Int_t i, Float_t v)
191 { if (i < 10) fSSECellCut[i] = v ; }
19391b8c 192
1253480f 193
19391b8c 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.);
198
2a77f6f4 199
992b14a7 200 //For histograms
cb99fa55 201 enum mcTypes { kmcPhoton = 1, kmcConversion = 2, kmcPi0 = 3, kmcPi0Conv = 4,
202 kmcEta = 5, kmcHadron = 6 };
992b14a7 203
204 private:
205
5c46c992 206 TString fCalorimeter ; // Calorimeter where the gamma is searched
3c1d9afb 207 Int_t fMinNCells ; // Study clusters with ncells larger than cut
208 Float_t fMinBadDist ; // Minimal distance to bad channel to accept cluster
b2e375c7 209 Float_t fHistoECut ; // Fixed E cut for some histograms
1253480f 210 Bool_t fCheckSplitDistToBad; // Check the distance to bad channel and to EMCal borders of split clusters
243c2909 211
883411b2 212 Bool_t fFillAngleHisto; // Fill splitted clusters angle histograms
8edbd100 213 Bool_t fFillTMHisto ; // Fill track matching histos,
8e81c2cf 214 Bool_t fFillTMResidualHisto ; // Fill track matching histos, residuals
215 Bool_t fFillSSExtraHisto ; // Fill shower shape extra histos
74e3eb22 216 Bool_t fFillMCHisto ; // Fill MC energy fraction histos
dbe09c26 217 Bool_t fFillSSWeightHisto ; // Fill weigth histograms
ff6aa4ce 218 Bool_t fFillNLMDiffCutHisto ; // Fill NLM histograms for different settings
19391b8c 219 Bool_t fFillEbinHisto ; // Fill E bin histograms
b2e375c7 220 Bool_t fFillMCOverlapHisto ; // Fill MC particles overlap histograms
4914e781 221 Bool_t fFillNCellHisto ; // Fill n cells in cluster dependent histograms
1253480f 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
cb99fa55 225 Bool_t fFillArmenterosHisto; // Fill armenteros type histo
29555e96 226 Bool_t fFillThetaStarHisto; // Fill cosThetaStar histos
1253480f 227
dbe09c26 228 Float_t fSSWeight[10]; // List of weights to test
19391b8c 229 Int_t fSSWeightN; // Total number of weights to test
230
231 Float_t fSSECellCut[10]; // List of cell min energy cuts to test
232 Int_t fSSECellCutN; // Total number of cell min energy cuts to test
ff6aa4ce 233
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
dbe09c26 237
8edbd100 238 Float_t fWSimu; // Slope of the linear correction factor for the shower
239 // shape weight in simulation, about 0.07
240
992b14a7 241 //Histograms
242
c7e9a4e9 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
5c46c992 246
c7e9a4e9 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
250
cb99fa55 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
254
e671adc2 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|
cb99fa55 258
2c36e041 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
cb99fa55 262
2c36e041 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
e671adc2 266
2c36e041 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
2a77f6f4 270
2c36e041 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
274
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
278
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
282
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
a6d3b0a8 286
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
e671adc2 290
77cadd95 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
667432ef 294
cb99fa55 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
0137016b 298
b8eb40fc 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
302
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
306
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
7b686344 310
b583134f 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
313
cb99fa55 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
d2655d46 317
e671adc2 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
d2655d46 321
cb99fa55 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
d2655d46 325
e671adc2 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
d2655d46 329
cb99fa55 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
d2655d46 333
e671adc2 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
d2655d46 337
cb99fa55 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
29555e96 340 TH2F * fhNLocMaxIdPi0 [7][2] ; //! Number of maxima in cluster vs E, 1-6 for different MC particle types, after pi0 selection
5c46c992 341
cb99fa55 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
29555e96 344
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
348
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
cc909e6f 352
cb99fa55 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
8e81c2cf 356
e671adc2 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
8e81c2cf 360
cb99fa55 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
4914e781 364
cb99fa55 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
4914e781 368
b8eb40fc 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
667432ef 372
cb99fa55 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
667432ef 376
cb99fa55 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
4914e781 380
667432ef 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
bb2d339b 384
cb99fa55 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
bb2d339b 388
cb99fa55 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
bb2d339b 392
cb99fa55 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
53f2c382 395
cb99fa55 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
8e81c2cf 399
cb99fa55 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
8e81c2cf 403
cb99fa55 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
992b14a7 407
4914e781 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
412
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
417
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
422
cb99fa55 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
c8710850 426
cb99fa55 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
5c46c992 430
cb99fa55 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
fc01318e 434
cb99fa55 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
e671adc2 438
cb99fa55 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
bb2d339b 442
cb99fa55 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
bb2d339b 446
cb99fa55 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
453
c7e9a4e9 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
457
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
461
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
167f2534 465
cb99fa55 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
bb2d339b 469
cb99fa55 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
fc01318e 473
cb99fa55 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
e671adc2 477
cb99fa55 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
667432ef 481
cb99fa55 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
667432ef 485
cb99fa55 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
5094c724 489
0186b6a2 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
c8710850 493
29555e96 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
497
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
501
0186b6a2 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
d261ae67 505
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
0186b6a2 509
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
513
29555e96 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
517
4e5a94c5 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
521
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
525
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
529
29555e96 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
533
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
537
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
541
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
545
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
549
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
c8710850 553
cb99fa55 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
560
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
567
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
b2e375c7 574
1253480f 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
c8710850 581
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
588
9554fc65 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
598
dbe09c26 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
19391b8c 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
607
dbe09c26 608 TH2F * fhM02WeightPi0[3][10] ; //! M02 for selected pi0 with different weight, with NLM = 1, = 2, > 2
19391b8c 609 TH2F * fhM02ECellCutPi0[3][10] ; //! M02 for selected pi0 with different cut on cell energy, with NLM = 1, = 2, > 2
dbe09c26 610
a1fd1b69 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
617
cb99fa55 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
624
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
631
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
638
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
645
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
652
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
659
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
cc909e6f 663
b2e375c7 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
670
cb99fa55 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
b2e375c7 673
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
b2e375c7 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
36769d30 683
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
693
83351853 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
36769d30 698 TH2F * fhMCPi0DecayPhotonAdjacent; //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM
83351853 699 TH2F * fhMCPi0DecayPhotonHitNoLM; //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas
b2e375c7 700
36769d30 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
707
b583134f 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
712
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
717
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
722
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
29555e96 727
728
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
733
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
738
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
743
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
748
b583134f 749
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
756
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
763
764
b2e375c7 765 TH2F * fhMCEOverlapType; //! what particles overlap with pi0, neutral clusters
766 TH2F * fhMCEOverlapTypeMatch; //! what particles overlap with pi0, charged clusters
dbe09c26 767
ce49dd72 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
ce49dd72 772
ff6aa4ce 773
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
777
2a77f6f4 778 AliAnaInsideClusterInvariantMass( const AliAnaInsideClusterInvariantMass & split) ; // cpy ctor
779 AliAnaInsideClusterInvariantMass & operator = (const AliAnaInsideClusterInvariantMass & split) ; // cpy assignment
992b14a7 780
ff6aa4ce 781 ClassDef(AliAnaInsideClusterInvariantMass,29)
992b14a7 782
783} ;
784
785#endif //ALIANAINSIDECLUSTERINVARIANTMASS_H
786
787
788