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