]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.h
correcting cone exess
[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//
12//-- Author: Gustavo Conesa (LPSC-Grenoble)
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
33 TObjString * GetAnalysisCuts();
34
35 TList * GetCreateOutputObjects();
dbe09c26 36
19391b8c 37 void FillSSWeightHistograms(AliVCluster *clus, const Int_t nlm, const Int_t absId1, const Int_t absId2);
dbe09c26 38
992b14a7 39 void Init();
40
41 void InitParameters();
42
43 void MakeAnalysisFillHistograms() ;
992b14a7 44
71e3889f 45 void Print(const Option_t * opt) const;
992b14a7 46
71e3889f 47 void SetCalorimeter(TString & det) { fCalorimeter = det ; }
48
49 void SetM02Cut(Float_t min=0, Float_t max=10) { fM02MinCut = min ; fM02MaxCut = max ; }
992b14a7 50
71e3889f 51 void SetMinNCells(Int_t cut) { fMinNCells = cut ; }
2cb134fb 52
53 void SetMinBadChannelDistance(Float_t cut) { fMinBadDist = cut ; }
992b14a7 54
8edbd100 55 void SetWCorrectionParameter(Float_t p = 0.07) { fWSimu = p ; }
56
883411b2 57 void SwitchOnFillAngleHistograms() { fFillAngleHisto = kTRUE ; }
58 void SwitchOffFillAngleHistograms() { fFillAngleHisto = kFALSE ; }
59
60 void SwitchOnFillExtraSSHistograms() { fFillSSExtraHisto = kTRUE ; }
61 void SwitchOffFillExtraSSHistograms() { fFillSSExtraHisto = kFALSE ; }
8e81c2cf 62
8edbd100 63 void SwitchOnFillTMHistograms() { fFillTMHisto = kTRUE ; }
64 void SwitchOffFillTMHistograms() { fFillTMHisto = kFALSE ; }
65
8e81c2cf 66 void SwitchOnFillTMResidualHistograms() { fFillTMResidualHisto = kTRUE ; }
67 void SwitchOffFillTMResidualHistograms() { fFillTMResidualHisto = kFALSE ; }
68
2a77f6f4 69 void SwitchOnMCFractionHistograms() { fFillMCFractionHisto = kTRUE ; }
70 void SwitchOffMCFractionHistograms() { fFillMCFractionHisto = kFALSE ; }
71
19391b8c 72 void SwitchOnFillSSWeightHistograms() { fFillSSWeightHisto = kTRUE ; }
73 void SwitchOffFillSSWeightHistograms() { fFillSSWeightHisto = kFALSE ; }
74
75 void SwitchOnFillEbinHistograms() { fFillEbinHisto = kTRUE ; }
76 void SwitchOffFillEbinHistograms() { fFillEbinHisto = kFALSE ; }
dbe09c26 77
19391b8c 78 void SetNWeightForShowerShape(Int_t n) { fSSWeightN = n ; }
dbe09c26 79 void SetWeightForShowerShape(Int_t i, Float_t v) { if (i < 10) fSSWeight[i] = v ; }
80
19391b8c 81 void SetNECellCutForShowerShape(Int_t n) { fSSECellCutN = n ; }
82 void SetECellCutForShowerShape(Int_t i, Float_t v) { if (i < 10) fSSECellCut[i] = v ; }
83
84
85 void RecalculateClusterShowerShapeParametersWithCellCut(const AliEMCALGeometry * geom, AliVCaloCells* cells, AliVCluster * cluster,
86 Float_t & l0, Float_t & l1,
87 Float_t & disp, Float_t & dEta, Float_t & dPhi,
88 Float_t & sEta, Float_t & sPhi, Float_t & sEtaPhi,Float_t eCellMin = 0.);
89
2a77f6f4 90
992b14a7 91 //For histograms
c5693f62 92 enum mcTypes { kmcPhoton = 1, kmcConversion = 2, kmcPi0 = 3,
de7d73e6 93 kmcEta = 4, kmcElectron = 5, kmcHadron = 6, kmcPi0Conv = 7 };
992b14a7 94
95 private:
96
5c46c992 97 TString fCalorimeter ; // Calorimeter where the gamma is searched
3c1d9afb 98 Float_t fM02MaxCut ; // Study clusters with l0 smaller than cut
99 Float_t fM02MinCut ; // Study clusters with l0 larger than cut
100 Int_t fMinNCells ; // Study clusters with ncells larger than cut
101 Float_t fMinBadDist ; // Minimal distance to bad channel to accept cluster
243c2909 102
883411b2 103 Bool_t fFillAngleHisto; // Fill splitted clusters angle histograms
8edbd100 104 Bool_t fFillTMHisto ; // Fill track matching histos,
8e81c2cf 105 Bool_t fFillTMResidualHisto ; // Fill track matching histos, residuals
106 Bool_t fFillSSExtraHisto ; // Fill shower shape extra histos
2a77f6f4 107 Bool_t fFillMCFractionHisto ; // Fill MC energy fraction histos
dbe09c26 108 Bool_t fFillSSWeightHisto ; // Fill weigth histograms
19391b8c 109 Bool_t fFillEbinHisto ; // Fill E bin histograms
8edbd100 110
dbe09c26 111 Float_t fSSWeight[10]; // List of weights to test
19391b8c 112 Int_t fSSWeightN; // Total number of weights to test
113
114 Float_t fSSECellCut[10]; // List of cell min energy cuts to test
115 Int_t fSSECellCutN; // Total number of cell min energy cuts to test
dbe09c26 116
8edbd100 117 Float_t fWSimu; // Slope of the linear correction factor for the shower
118 // shape weight in simulation, about 0.07
119
992b14a7 120 //Histograms
121
de7d73e6 122 TH2F * fhMassNLocMax1[8][2] ; //! Mass of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types
123 TH2F * fhMassNLocMax2[8][2] ; //! Mass of 2 cells local maxima vs E, 1-6 for different MC particle types
124 TH2F * fhMassNLocMaxN[8][2] ; //! Mass of >2 cells local maxima vs E, 1-6 for different MC particle types
5c46c992 125
de7d73e6 126 TH2F * fhAsymNLocMax1[8][2] ; //! Asymmetry of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types
127 TH2F * fhAsymNLocMax2[8][2] ; //! Asymmetry of 2 cells local maxima vs E, 1-6 for different MC particle types
128 TH2F * fhAsymNLocMaxN[8][2] ; //! Asymmetry of >2 cells local maxima vs E, 1-6 for different MC particle types
2a77f6f4 129
e671adc2 130 TH2F * fhSplitEFractionvsAsyNLocMax1[2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1 vs |A|
131 TH2F * fhSplitEFractionvsAsyNLocMax2[2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2 vs |A|
132 TH2F * fhSplitEFractionvsAsyNLocMaxN[2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2 vs |A|
133
77cadd95 134 TH2F * fhMassM02CutNLocMax1 ; //! M02(E) selection, not matched, Mass of split clusters, NLM = 1
135 TH2F * fhMassM02CutNLocMax2 ; //! M02(E) selection, not matched, Mass of split clusters, NLM = 1
136 TH2F * fhMassM02CutNLocMaxN ; //! M02(E) selection, not matched, Mass of split clusters, NLM > 2
2a77f6f4 137
77cadd95 138 TH2F * fhAsymM02CutNLocMax1 ; //! M02(E) selection, not matched, energy asymmetry of split clusters, NLM = 1
139 TH2F * fhAsymM02CutNLocMax2 ; //! M02(E) selection, not matched, energy asymmetry of split clusters, NLM = 2
140 TH2F * fhAsymM02CutNLocMaxN ; //! M02(E) selection, not matched, energy asymmetry of split clusters, NLM > 2
e671adc2 141
77cadd95 142 TH2F * fhMassSplitECutNLocMax1 ; //! 85% of split energy, not matched, Mass of split clusters, NLM = 1
143 TH2F * fhMassSplitECutNLocMax2 ; //! 85% of split energy, not matched, Mass of split clusters, NLM = 1
144 TH2F * fhMassSplitECutNLocMaxN ; //! 85% of split energy, not matched, Mass of split clusters, NLM > 2
667432ef 145
de7d73e6 146 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
147 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
148 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 149
b8eb40fc 150 TH2F * fhMassM02NLocMax1Ebin[4] ; //! Mass of splitted clusters when 1 local max vs M02, 4 E bins, neutral clusters
151 TH2F * fhMassM02NLocMax2Ebin[4] ; //! Mass of splitted clusters when 2 local max vs M02, 4 E bins, neutral clusters
152 TH2F * fhMassM02NLocMaxNEbin[4] ; //! Mass of splitted clusters when >2 local max vs M02, 4 E bins, neutral clusters
153
154 TH2F * fhMassAsyNLocMax1Ebin[4] ; //! Mass of Mass of splitted clusters when 1 local max vs asymmetry, 4 E bins, neutral clusters
155 TH2F * fhMassAsyNLocMax2Ebin[4] ; //! Mass of Mass of splitted clusters when 2 local max vs asymmetry, 4 E bins, neutral clusters
156 TH2F * fhMassAsyNLocMaxNEbin[4] ; //! Mass of Mass of splitted clusters when >2 local max vs asymmetry, 4 E bins, neutral clusters
157
158 TH2F * fhAsyMCGenRecoNLocMax1EbinPi0[4] ; //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when 1 local max, 4 E bins, neutral clusters
159 TH2F * fhAsyMCGenRecoNLocMax2EbinPi0[4] ; //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when 2 local max, 4 E bins, neutral clusters
160 TH2F * fhAsyMCGenRecoNLocMaxNEbinPi0[4] ; //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when >2 local max, 4 E bins, neutral clusters
7b686344 161
de7d73e6 162 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
163 TH2F * fhMassDispEtaNLocMax2[8][2] ; //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
164 TH2F * fhMassDispEtaNLocMaxN[8][2] ; //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
d2655d46 165
e671adc2 166 TH2F * fhMassDispEtaNLocMax1Ebin[4] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
167 TH2F * fhMassDispEtaNLocMax2Ebin[4] ; //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
168 TH2F * fhMassDispEtaNLocMaxNEbin[4] ; //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
d2655d46 169
de7d73e6 170 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
171 TH2F * fhMassDispPhiNLocMax2[8][2] ; //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
172 TH2F * fhMassDispPhiNLocMaxN[8][2] ; //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
d2655d46 173
e671adc2 174 TH2F * fhMassDispPhiNLocMax1Ebin[4] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
175 TH2F * fhMassDispPhiNLocMax2Ebin[4] ; //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
176 TH2F * fhMassDispPhiNLocMaxNEbin[4] ; //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
d2655d46 177
de7d73e6 178 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
179 TH2F * fhMassDispAsyNLocMax2[8][2] ; //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
180 TH2F * fhMassDispAsyNLocMaxN[8][2] ; //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
d2655d46 181
e671adc2 182 TH2F * fhMassDispAsyNLocMax1Ebin[4] ; //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
183 TH2F * fhMassDispAsyNLocMax2Ebin[4] ; //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
184 TH2F * fhMassDispAsyNLocMaxNEbin[4] ; //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
d2655d46 185
de7d73e6 186 TH2F * fhNLocMax [8][2] ; //! Number of maxima in cluster vs E, 1-6 for different MC particle types
187 TH2F * fhNLocMaxM02Cut[8][2] ; //! Number of maxima in cluster vs E, 1-6 for different MC particle types, after SS cut
5c46c992 188
de7d73e6 189 TH2F * fhM02NLocMax1 [8][2] ; //! M02 vs E for N max in cluster = 1, 1-6 for different MC particle types
190 TH2F * fhM02NLocMax2 [8][2] ; //! M02 vs E for N max in cluster = 2, 1-6 for different MC particle types
191 TH2F * fhM02NLocMaxN [8][2] ; //! M02 vs E for N max in cluster > 2, 1-6 for different MC particle types
8e81c2cf 192
e671adc2 193 TH2F * fhMCAsymM02NLocMax1MCPi0Ebin[4] ; //! M02 vs decay asymmetry for N max in cluster = 1, for 4 energy bins
194 TH2F * fhMCAsymM02NLocMax2MCPi0Ebin[4] ; //! M02 vs decay asymmetry for N max in cluster = 2, for 4 energy bins
195 TH2F * fhMCAsymM02NLocMaxNMCPi0Ebin[4] ; //! M02 vs decay asymmetry for N max in cluster > 2, for 4 energy bins
8e81c2cf 196
de7d73e6 197 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
198 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
199 TH2F * fhMCGenFracNLocMaxN[8][2] ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types
667432ef 200
b8eb40fc 201 TH2F * fhMCGenFracAfterCutsNLocMax1MCPi0 ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, MCPi0 after M02 and asymmetry cut
202 TH2F * fhMCGenFracAfterCutsNLocMax2MCPi0 ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, MCPi0, after M02 and asymmetry cut
203 TH2F * fhMCGenFracAfterCutsNLocMaxNMCPi0 ; //! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, MCPi0, after M02 and asymmetry cut
667432ef 204
de7d73e6 205 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
206 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
207 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 208
209 TH2F * fhMCGenSplitEFracAfterCutsNLocMax1MCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types
210 TH2F * fhMCGenSplitEFracAfterCutsNLocMax2MCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types
211 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 212
de7d73e6 213 TH2F * fhMCGenEFracvsSplitEFracNLocMax1[8][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster = 1, MC pi0
214 TH2F * fhMCGenEFracvsSplitEFracNLocMax2[8][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster = 2, MC pi0
215 TH2F * fhMCGenEFracvsSplitEFracNLocMaxN[8][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster > 2, MC pi0
bb2d339b 216
de7d73e6 217 TH2F * fhMCGenEvsSplitENLocMax1[8][2] ; //! E generated particle vs E1+E2 for N max in cluster = 1, 1-6 for different MC particle types
218 TH2F * fhMCGenEvsSplitENLocMax2[8][2] ; //! E generated particle vs E1+E2 for N max in cluster = 2, 1-6 for different MC particle types
219 TH2F * fhMCGenEvsSplitENLocMaxN[8][2] ; //! E generated particle vs E1+E2 for N max in cluster > 2, 1-6 for different MC particle types
bb2d339b 220
de7d73e6 221 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
222 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 223
de7d73e6 224 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
225 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
226 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 227
de7d73e6 228 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
229 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
230 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 231
de7d73e6 232 TH2F * fhNCellNLocMax1[8][2] ; //! n cells in cluster vs E for N max in cluster = 1, 1-6 for different MC particle types
233 TH2F * fhNCellNLocMax2[8][2] ; //! n cells in cluster vs E for N max in cluster = 2, 1-6 for different MC particle types
234 TH2F * fhNCellNLocMaxN[8][2] ; //! n cells in cluster vs E for N max in cluster > 2, 1-6 for different MC particle types
992b14a7 235
c8710850 236 TH2F * fhM02Pi0NLocMax1[8][2] ; //! M02 for Mass around pi0, N Local Maxima = 1
237 TH2F * fhM02EtaNLocMax1[8][2] ; //! M02 for Mass around eta, N Local Maxima = 1
238 TH2F * fhM02ConNLocMax1[8][2] ; //! M02 for Mass around close to 0, N Local Maxima = 1
239
240 TH2F * fhM02Pi0NLocMax2[8][2] ; //! M02 for Mass around pi0, N Local Maxima = 2
241 TH2F * fhM02EtaNLocMax2[8][2] ; //! M02 for Mass around eta, N Local Maxima = 2
242 TH2F * fhM02ConNLocMax2[8][2] ; //! M02 for Mass around close to 0, N Local Maxima = 2
5c46c992 243
c8710850 244 TH2F * fhM02Pi0NLocMaxN[8][2] ; //! M02 for Mass around pi0, N Local Maxima > 2
245 TH2F * fhM02EtaNLocMaxN[8][2] ; //! M02 for Mass around eta, N Local Maxima > 2
246 TH2F * fhM02ConNLocMaxN[8][2] ; //! M02 for Mass around close to 0, N Local Maxima > 2
fc01318e 247
c8710850 248 TH2F * fhMassPi0NLocMax1[8][2] ; //! Mass for selected pi0, N Local Maxima = 1
249 TH2F * fhMassEtaNLocMax1[8][2] ; //! Mass for selected around eta, N Local Maxima = 1
250 TH2F * fhMassConNLocMax1[8][2] ; //! Mass for selected around close to 0, N Local Maxima = 1
e671adc2 251
c8710850 252 TH2F * fhMassPi0NLocMax2[8][2] ; //! Mass for selected around pi0, N Local Maxima = 2
253 TH2F * fhMassEtaNLocMax2[8][2] ; //! Mass for selected around eta, N Local Maxima = 2
254 TH2F * fhMassConNLocMax2[8][2] ; //! Mass for selected around close to 0, N Local Maxima = 2
bb2d339b 255
c8710850 256 TH2F * fhMassPi0NLocMaxN[8][2] ; //! Mass for selected around pi0, N Local Maxima > 2
257 TH2F * fhMassEtaNLocMaxN[8][2] ; //! Mass for selected around eta, N Local Maxima > 2
258 TH2F * fhMassConNLocMaxN[8][2] ; //! Mass for selected around close to 0, N Local Maxima > 2
bb2d339b 259
74858845 260 TH2F * fhMassAfterCutsNLocMax1[8][2] ; //! Mass after M02, asymmetry cuts for pi0, N Local Maxima = 1
261 TH2F * fhMassAfterCutsNLocMax2[8][2] ; //! Mass after M02, asymmetry cuts for pi0, N Local Maxima = 2
262 TH2F * fhMassAfterCutsNLocMaxN[8][2] ; //! Mass after M02, asymmetry cuts for pi0, N Local Maxima > 2
167f2534 263
c8710850 264 TH2F * fhAsyPi0NLocMax1[8][2] ; //! Asy for Mass around pi0, N Local Maxima = 1
265 TH2F * fhAsyEtaNLocMax1[8][2] ; //! Asy for Mass around eta, N Local Maxima = 1
266 TH2F * fhAsyConNLocMax1[8][2] ; //! Asy for Mass around close to 0, N Local Maxima = 1
bb2d339b 267
c8710850 268 TH2F * fhAsyPi0NLocMax2[8][2] ; //! Asy for Mass around pi0, N Local Maxima = 2
269 TH2F * fhAsyEtaNLocMax2[8][2] ; //! Asy for Mass around eta, N Local Maxima = 2
270 TH2F * fhAsyConNLocMax2[8][2] ; //! Asy for Mass around close to 0, N Local Maxima = 2
fc01318e 271
c8710850 272 TH2F * fhAsyPi0NLocMaxN[8][2] ; //! Asy for Mass around pi0, N Local Maxima > 2
273 TH2F * fhAsyEtaNLocMaxN[8][2] ; //! Asy for Mass around eta, N Local Maxima > 2
274 TH2F * fhAsyConNLocMaxN[8][2] ; //! Asy for Mass around close to 0, N Local Maxima > 2
e671adc2 275
de7d73e6 276 TH2F * fhSplitEFractionNLocMax1[8][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1
277 TH2F * fhSplitEFractionNLocMax2[8][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2
278 TH2F * fhSplitEFractionNLocMaxN[8][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2
667432ef 279
74858845 280 TH2F * fhSplitEFractionAfterCutsNLocMax1[8][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1, after M02 and asymmetry cut
281 TH2F * fhSplitEFractionAfterCutsNLocMax2[8][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2, after M02 and asymmetry cut
282 TH2F * fhSplitEFractionAfterCutsNLocMaxN[8][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2, after M02 and asymmetry cut
667432ef 283
de7d73e6 284 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
285 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
286 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 287
c8710850 288 TH2F * fhAnglePairNLocMax1[2] ; //! pair opening angle vs E
289 TH2F * fhAnglePairNLocMax2[2] ; //! pair opening angle vs E
290 TH2F * fhAnglePairNLocMaxN[2] ; //! pair opening angle vs E
291
292 TH2F * fhAnglePairMassNLocMax1[2] ; //! pair opening angle vs Mass for E > 7 GeV
293 TH2F * fhAnglePairMassNLocMax2[2] ; //! pair opening angle vs Mass for E > 7 GeV
294 TH2F * fhAnglePairMassNLocMaxN[2] ; //! pair opening angle vs Mass for E > 7 GeV
295
296 TH2F * fhTrackMatchedDEtaNLocMax1[8] ; //! Eta distance between track and cluster vs cluster E, 1 local maximum
297 TH2F * fhTrackMatchedDPhiNLocMax1[8] ; //! Phi distance between track and cluster vs cluster E, 1 local maximum
298 TH2F * fhTrackMatchedDEtaNLocMax2[8] ; //! Eta distance between track and cluster vs cluster E, 2 local maximum
299 TH2F * fhTrackMatchedDPhiNLocMax2[8] ; //! Phi distance between track and cluster vs cluster E, 2 local maximum
300 TH2F * fhTrackMatchedDEtaNLocMaxN[8] ; //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
301 TH2F * fhTrackMatchedDPhiNLocMaxN[8] ; //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
302
17f5b4b6 303 TH2F * fhCentralityPi0NLocMax1[8][2] ; //! Centrality for selected pi0, N Local Maxima = 1
304 TH2F * fhCentralityEtaNLocMax1[8][2] ; //! Centrality for selected eta, N Local Maxima = 1
305 TH2F * fhCentralityPi0NLocMax2[8][2] ; //! Centrality for selected pi0, N Local Maxima = 2
306 TH2F * fhCentralityEtaNLocMax2[8][2] ; //! Centrality for selected eta, N Local Maxima = 2
307 TH2F * fhCentralityPi0NLocMaxN[8][2] ; //! Centrality for selected pi0, N Local Maxima > 2
308 TH2F * fhCentralityEtaNLocMaxN[8][2] ; //! Centrality for selected eta, N Local Maxima > 2
c8710850 309
310 TH2F * fhEventPlanePi0NLocMax1 ; //! Event plane for selected pi0, N Local Maxima = 1
311 TH2F * fhEventPlaneEtaNLocMax1 ; //! Event plane for selected eta, N Local Maxima = 1
312 TH2F * fhEventPlanePi0NLocMax2 ; //! Event plane for selected pi0, N Local Maxima = 2
313 TH2F * fhEventPlaneEtaNLocMax2 ; //! Event plane for selected eta, N Local Maxima = 2
314 TH2F * fhEventPlanePi0NLocMaxN ; //! Event plane for selected pi0, N Local Maxima > 2
315 TH2F * fhEventPlaneEtaNLocMaxN ; //! Event plane for selected eta, N Local Maxima > 2
316
9554fc65 317 TH2F * fhClusterEtaPhiNLocMax1 ; //! Eta vs Phi of clusters with N Local Maxima = 1, E > 8 GeV
318 TH2F * fhClusterEtaPhiNLocMax2 ; //! Eta vs Phi of clusters with N Local Maxima = 2, E > 8 GeV
319 TH2F * fhClusterEtaPhiNLocMaxN ; //! Eta vs Phi of clusters with N Local Maxima > 2, E > 8 GeV
320 TH2F * fhPi0EtaPhiNLocMax1 ; //! Eta vs Phi of pi0's with N Local Maxima = 1, E > 8 GeV
321 TH2F * fhPi0EtaPhiNLocMax2 ; //! Eta vs Phi of pi0's with N Local Maxima = 2, E > 8 GeV
322 TH2F * fhPi0EtaPhiNLocMaxN ; //! Eta vs Phi of pi0's with N Local Maxima > N, E > 8 GeV
323 TH2F * fhEtaEtaPhiNLocMax1 ; //! Eta vs Phi of eta's with N Local Maxima = 1, E > 8 GeV
324 TH2F * fhEtaEtaPhiNLocMax2 ; //! Eta vs Phi of eta's with N Local Maxima = 2, E > 8 GeV
325 TH2F * fhEtaEtaPhiNLocMaxN ; //! Eta vs Phi of eta's with N Local Maxima > N, E > 8 GeV
326
dbe09c26 327 TH2F * fhPi0CellE[3] ; //! pi0's energy vs cluster cell energy with NLM = 1, = 2, > 2
328 TH2F * fhPi0CellEFrac[3] ; //! pi0's energy vs cluster cell energy fraction with NLM = 1, = 2, > 2
329 TH2F * fhPi0CellLogEFrac[3] ; //! pi0's energy vs cluster log cell energy fraction with NLM = 1, = 2, > 2
19391b8c 330 TH2F * fhPi0CellEMaxEMax2Frac [3]; //! pi0's energy vs fraction of 2 main maxima energy with NLM = 1, = 2, > 2
331 TH2F * fhPi0CellEMaxClusterFrac [3]; //! pi0's energy vs energy fraction of main LM and cluster energy with NLM = 1, = 2, > 2
332 TH2F * fhPi0CellEMax2ClusterFrac[3]; //! pi0's energy vs energy fraction of second LM and cluster energy with NLM = 1, = 2, > 2
333 TH2F * fhPi0CellEMaxFrac [3]; //! pi0's energy vs energy fraction of main LM and cluster cell energy with NLM = 1, = 2, > 2
334 TH2F * fhPi0CellEMax2Frac [3]; //! pi0's energy vs energy fraction of second LM and cluster cell energy with NLM = 1, = 2, > 2
335
dbe09c26 336 TH2F * fhM02WeightPi0[3][10] ; //! M02 for selected pi0 with different weight, with NLM = 1, = 2, > 2
19391b8c 337 TH2F * fhM02ECellCutPi0[3][10] ; //! M02 for selected pi0 with different cut on cell energy, with NLM = 1, = 2, > 2
dbe09c26 338
339
2a77f6f4 340 AliAnaInsideClusterInvariantMass( const AliAnaInsideClusterInvariantMass & split) ; // cpy ctor
341 AliAnaInsideClusterInvariantMass & operator = (const AliAnaInsideClusterInvariantMass & split) ; // cpy assignment
992b14a7 342
9554fc65 343 ClassDef(AliAnaInsideClusterInvariantMass,20)
992b14a7 344
345} ;
346
347#endif //ALIANAINSIDECLUSTERINVARIANTMASS_H
348
349
350