ATLAS Offline Software
Loading...
Searching...
No Matches
ShowerShapesPlots.cxx
Go to the documentation of this file.
1/*
2 Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
3*/
4
5#include <utility>
6
7#include "ShowerShapesPlots.h"
8
9using CLHEP::GeV;
10
11namespace Egamma{
12
13
14 ShowerShapesPlots::ShowerShapesPlots(PlotBase* pParent, const std::string& sDir, std::string sParticleType):PlotBase(pParent, sDir),
15 m_sParticleType(std::move(sParticleType)),
16 hadleak(nullptr),
17 middlee(nullptr),
18 reta37 (nullptr),
19 reta33 (nullptr),
20 weta1 (nullptr),
21 weta2 (nullptr),
22 demax1 (nullptr),
23 demax2 (nullptr),
24 de (nullptr),
25 fracs1 (nullptr),
26 wtots1 (nullptr),
27 f1 (nullptr),
28 f3 (nullptr),
29 pos7 (nullptr),
30 clusiso(nullptr),
31 Eratio(nullptr),
32 E0(nullptr),
33 E1(nullptr),
34 E2(nullptr),
35 E3(nullptr),
36 Rhad(nullptr),
37 Rhad1(nullptr),
38 Reta(nullptr),
39 Rphi(nullptr),
40 hadleakvset(nullptr),
41 middleevset(nullptr),
42 reta37vset (nullptr),
43 reta33vset (nullptr),
44 weta1vset (nullptr),
45 weta2vset (nullptr),
46 demax1vset (nullptr),
47 demax2vset (nullptr),
48 devset (nullptr),
49 fracs1vset (nullptr),
50 wtots1vset (nullptr),
51 f1vset (nullptr),
52 f3vset (nullptr),
53 pos7vset (nullptr),
54 clusisovset(nullptr),
55 Eratiovset(nullptr),
56 Rhadvset(nullptr),
57 Rhad1vset(nullptr),
58 Retavset(nullptr),
59 Rphivset(nullptr),
60 hadleakvseta(nullptr),
61 middleevseta(nullptr),
62 reta37vseta(nullptr),
63 reta33vseta(nullptr),
64 weta1vseta(nullptr),
65 weta2vseta(nullptr),
66 demax1vseta(nullptr),
67 demax2vseta(nullptr),
68 devseta(nullptr),
69 fracs1vseta(nullptr),
70 wtots1vseta(nullptr),
71 f1vseta(nullptr),
72 f3vseta(nullptr),
73 pos7vseta(nullptr),
74 clusisovseta(nullptr),
75 Eratiovseta(nullptr),
76 Rhadvseta(nullptr),
77 Rhad1vseta(nullptr),
78 Retavseta(nullptr),
79 Rphivseta(nullptr)
80
81 {}
82
84
85 hadleak = Book1D("hadleak", "Hadronic leakage of " + m_sParticleType+"; E_{hadleak} (GeV);Entries", 100, -0.07, 0.13);
86 middlee = Book1D("middlee", "E_{middle} of "+ m_sParticleType+"; E_{277} (GeV); Entries" , 210,-10.,200.);
87 reta37 = Book1D("reta37", "R_{#eta37} of "+ m_sParticleType+";R_{#eta37};Entries", 355, 0., 1.1005);
88 reta33 = Book1D("reta33", "R_{#eta33} of "+ m_sParticleType+";R_{#eta33};Entries", 355, 0., 1.1005);
89 weta1 = Book1D("weta1", "W_{#etas1} of "+ m_sParticleType+";W_{#etas1};Entries", 100, 0., 1.);
90 weta2 = Book1D("weta2", "W_{#etas2} of "+ m_sParticleType+";W_{#etas2};Entries", 100, 0., 0.03);
91 demax1 = Book1D("demax1", "#DeltaE_{maxs1} of "+ m_sParticleType+";#DeltaE_{maxs1};Entries" , 100, 0., 1.);
92 demax2 = Book1D("demax2", "#DeltaE_{maxs2} of "+ m_sParticleType+";#DeltaE_{maxs2};Entries" , 80, 0., 0.8);
93 de = Book1D("de", "#DeltaE of "+ m_sParticleType+";#DeltaE (GeV);Entries", 250, 0., 0.5);
94 fracs1 = Book1D("fracs1", "Fracs1 of "+ m_sParticleType+";Fracs1;Entries", 350, 0., 3.5);
95 wtots1 = Book1D("wtots1", "W_{tots1} of "+ m_sParticleType+";W_{tots1};Entries", 100, 0., 10.);
96 f1 = Book1D("f1", "f1 of "+ m_sParticleType+";f1;Entries" , 120, -0.2, 1.0);
97 f3 = Book1D("f3", "f3 of "+ m_sParticleType+";f3;Entries" , 120, -0.2, 1.0);
98 Eratio = Book1D("Eratio", "Eratio of "+ m_sParticleType+";Eratio;Entries" , 100, 0., 1.0);
99 pos7 = Book1D("pos7", "pos7 of "+ m_sParticleType+";pos7;Entries", 80, -4., 4.);
100 clusiso = Book1D("clusiso","Cluster Isolation of "+ m_sParticleType+"; E_{T}^{cone20}/E_{T};Entries", 100, 0., 1.);
101 E0 = Book1D("E0", "E0 of "+ m_sParticleType+"; E0 (GeV); Entries" , 210,-10.,200.);
102 E1 = Book1D("E1", "E1 of "+ m_sParticleType+"; E1 (GeV); Entries" , 210,-10.,200.);
103 E2 = Book1D("E2", "E2 of "+ m_sParticleType+"; E2 (GeV); Entries" , 210,-10.,200.);
104 E3 = Book1D("E3", "E3 of "+ m_sParticleType+"; E3 (GeV); Entries" , 210,-10.,200.);
105 Rhad = Book1D("Rhad", "Rhad of "+ m_sParticleType+";Rhad;Entries", 100, -0.5, 2.);
106 Rhad1 = Book1D("Rhad1", "Rhad1 of "+ m_sParticleType+";Rhad1;Entries", 100, -0.5, 1.5);
107 Reta = Book1D("Reta", "Reta of "+ m_sParticleType+";Reta;Entries", 100, -0., 3.);
108 Rphi = Book1D("Rphi", "Rphi of "+ m_sParticleType+";Rphi;Entries", 100, -1., 1.);
109
110 hadleakvset = Book2D("hadleakvset", "Hadronic leakage vs E_{T} of " + m_sParticleType+"; E_{hadleak} (GeV) ; E_{T} (GeV) ", 100, -0.07, 0.13, 200, 0., 200);
111 middleevset = Book2D("middleevset", "E_{middle} vs E_{T} of "+ m_sParticleType+"; E_{277} (GeV); E_{T} (GeV) " , 210,-10.,200., 200, 0., 200);
112 reta37vset = Book2D("reta37vset", "R_{#eta37} vs E_{T} of "+ m_sParticleType+";R_{#eta37}; E_{T} (GeV) ", 355, 0., 1.1005, 200, 0., 200);
113 reta33vset = Book2D("reta33vset", "R_{#eta33} vs E_{T} of "+ m_sParticleType+";R_{#eta33}; E_{T} (GeV) ", 355, 0., 1.1005, 200, 0., 200);
114 weta1vset = Book2D("weta1vset", "W_{#etas1} vs E_{T} of "+ m_sParticleType+";W_{#etas1}; E_{T} (GeV) ", 100, 0., 1., 200, 0., 200);
115 weta2vset = Book2D("weta2vset", "W_{#etas2} vs E_{T} of "+ m_sParticleType+";W_{#etas2}; E_{T} (GeV) ", 100, 0., 0.03, 200, 0., 200);
116 demax1vset = Book2D("demax1vset", "#DeltaE_{maxs1} vs E_{T} of "+ m_sParticleType+";#DeltaE_{maxs1}; E_{T} (GeV) " , 100, 0., 1., 200, 0., 200);
117 demax2vset = Book2D("demax2vset", "#DeltaE_{maxs2} vs E_{T} of "+ m_sParticleType+";#DeltaE_{maxs2}; E_{T} (GeV) " , 80, 0., 0.8, 200, 0., 200);
118 devset = Book2D("devset", "#DeltaE vs E_{T} of "+ m_sParticleType+";#DeltaE (GeV); E_{T} (GeV) ", 250, 0., 0.5, 200, 0., 200);
119 fracs1vset = Book2D("fracs1vset", "Fracs1 vs E_{T} of "+ m_sParticleType+";Fracs1; E_{T} (GeV) ", 350, 0., 3.5, 200, 0., 200);
120 wtots1vset = Book2D("wtots1vset", "W_{tots1} vsE_{T} of "+ m_sParticleType+";W_{tots1}; E_{T} (GeV) ", 100, 0., 10., 200, 0., 200);
121 f1vset = Book2D("f1vset", "f1 vs E_{T} of "+ m_sParticleType+";f1; E_{T} (GeV) " , 100, 0., 1.0, 200, 0., 200);
122 f3vset = Book2D("f3vset", "f3 vs E_{T} of "+ m_sParticleType+";f3; E_{T} (GeV) " , 120, -0.2, 1.0, 200, 0., 200);
123 Eratiovset = Book2D("Eratiovset", "Eratio vs E_{T} of "+ m_sParticleType+";Eratio; E_{T} (GeV) " , 100, 0., 1.0, 200, 0., 200);
124 pos7vset = Book2D("pos7vset", "pos7 vs E_{T} of "+ m_sParticleType+";pos7; E_{T} (GeV) ", 80, -4., 4., 200, 0., 200);
125 clusisovset = Book2D("clusisovset","Cluster Isolation vs E_{T} of "+ m_sParticleType+"; E_{T}^{cone20}/E_{T}; E_{T} (GeV) ", 100, 0., 1., 200, 0., 200);
126 Rhadvset = Book2D("Rhadvset", "Rhad vs E_{T} of "+ m_sParticleType+";Rhad;E_{T} (GeV)", 100, -0.5, 2., 200, 0., 200);
127 Rhad1vset = Book2D("Rhad1vset", "Rhad1 vs E_{T} of "+ m_sParticleType+";Rhad1;E_{T} (GeV)", 100, -0.5, 1.5, 200, 0., 200);
128 Retavset = Book2D("Retavset", "Reta vs E_{T} of "+ m_sParticleType+";Reta;E_{T} (GeV)", 100, -0., 3., 200, 0., 200);
129 Rphivset = Book2D("Rphivset", "Rphi vs E_{T} of "+ m_sParticleType+";Rphi;E_{T} (GeV)", 100, -1., 1., 200, 0., 200);
130
131 hadleakvseta = Book2D("hadleakvseta", "Hadronic leakage vs E_{T} of " + m_sParticleType+"; E_{hadleak} (GeV) ; #eta ", 100, -0.07, 0.13, 1000,-5.,5.);
132 middleevseta = Book2D("middleevseta", "E_{middle} vs E_{T} of "+ m_sParticleType+"; E_{277} (GeV); #eta " , 210,-10.,200., 1000,-5.,5.);
133 reta37vseta = Book2D("reta37vseta", "R_{#eta37} vs E_{T} of "+ m_sParticleType+";R_{#eta37}; #eta ", 355, 0., 1.1005, 1000,-5.,5.);
134 reta33vseta = Book2D("reta33vseta", "R_{#eta33} vs E_{T} of "+ m_sParticleType+";R_{#eta33}; #eta ", 355, 0., 1.1005, 1000,-5.,5.);
135 weta1vseta = Book2D("weta1vseta", "W_{#etas1} vs E_{T} of "+ m_sParticleType+";W_{#etas1}; #eta ", 100, 0., 1., 1000,-5.,5.);
136 weta2vseta = Book2D("weta2vseta", "W_{#etas2} vs E_{T} of "+ m_sParticleType+";W_{#etas2}; #eta ", 100, 0., 0.03, 1000,-5.,5.);
137 demax1vseta = Book2D("demax1vseta", "#DeltaE_{maxs1} vs E_{T} of "+ m_sParticleType+";#DeltaE_{maxs1}; #eta " , 100, 0., 1., 1000,-5.,5.);
138 demax2vseta = Book2D("demax2vseta", "#DeltaE_{maxs2} vs E_{T} of "+ m_sParticleType+";#DeltaE_{maxs2}; #eta " , 80, 0., 0.8, 1000,-5.,5.);
139 devseta = Book2D("devseta", "#DeltaE vs E_{T} of "+ m_sParticleType+";#DeltaE (GeV); #eta ", 250, 0., 0.5, 1000,-5.,5.);
140 fracs1vseta = Book2D("fracs1vseta", "Fracs1 vs E_{T} of "+ m_sParticleType+";Fracs1; #eta ", 350, 0., 3.5, 1000,-5.,5.);
141 wtots1vseta = Book2D("wtots1vseta", "W_{tots1} vsE_{T} of "+ m_sParticleType+";W_{tots1}; #eta ", 100, 0., 10., 1000,-5.,5.);
142 f1vseta = Book2D("f1vseta", "f1 vs E_{T} of "+ m_sParticleType+";f1; #eta " , 100, 0., 1.0, 1000,-5.,5.);
143 f3vseta = Book2D("f3vseta", "f3 vs #eta of "+ m_sParticleType+";f3; #eta " , 120, -0.2, 1.0, 1000,-5.,5.);
144 Eratiovseta = Book2D("Eratiovseta", "Eratio vs #eta of "+ m_sParticleType+";Eratio; #eta " , 100, 0., 1.0, 1000,-5.,5.);
145 pos7vseta = Book2D("pos7vseta", "pos7 vs #eta of "+ m_sParticleType+";pos7; #eta ", 80, -4., 4., 1000,-5.,5.);
146 clusisovseta = Book2D("clusisovseta","Cluster Isolation vs #eta of "+ m_sParticleType+"; E_{T}^{topoEtcone20}/E_{T}; #eta ", 100, 0., 1., 1000,-5.,5.);
147 Rhadvseta = Book2D("Rhadvseta", "Rhad vs #eta of "+ m_sParticleType+";Rhad;#eta", 100, -0.5, 2., 1000,-5.,5.);
148 Rhad1vseta = Book2D("Rhad1vseta", "Rhad1 vs #eta of "+ m_sParticleType+";Rhad1;#eta", 100, -0.5, 1.5, 1000,-5.,5.);
149 Retavseta = Book2D("Retavseta", "Reta vs #eta of "+ m_sParticleType+";Reta;#eta", 100, -0., 3., 1000,-5.,5.);
150 Rphivseta = Book2D("Rphivseta", "Rphi vs #eta of "+ m_sParticleType+";Rphi;#eta", 100, -1., 1., 1000,-5.,5.);
151
152 }
153
155
156 float weight = 1.;
157 weight = eventInfo.beamSpotWeight();
158
159 float eta2 = fabs (egamma.caloCluster()->etaBE (2));
160 float et37 = egamma.caloCluster()->e() / cosh (eta2);
161 float ethad(0);
162 float ethad1(0);
163 float raphad(0);
164 float raphad1(0);
165 if(egamma.showerShapeValue(ethad, xAOD::EgammaParameters::ethad )){
166 raphad = et37 > 0. ? ethad / et37 : 0.;
167 }
168 if(egamma.showerShapeValue(ethad1, xAOD::EgammaParameters::ethad1 )){
169 raphad1 = et37 > 0. ? ethad1 / et37 : 0.;
170 }
171 float hadrleak = (eta2 >= 0.8 && eta2 < 1.37) ? raphad : raphad1;
172 hadleak->Fill(hadrleak, weight);
173 hadleakvset->Fill(hadrleak, egamma.pt()/GeV, weight);
174 hadleakvseta->Fill(hadrleak, egamma.eta(), weight);
175
176 float e277(0);
177 float e237(0);
178 float e233(0);
179
180 if(egamma.showerShapeValue(e277, xAOD::EgammaParameters::e277 ) &&
181 egamma.showerShapeValue(e237, xAOD::EgammaParameters::e237 ) &&
182 egamma.showerShapeValue(e233, xAOD::EgammaParameters::e233 ) ){
183
184 middlee->Fill(e277/GeV, weight);
185 middleevset->Fill(e277/GeV, egamma.pt()/GeV, weight);
186 middleevseta->Fill(e277/GeV, egamma.eta(), weight);
187 float rateta37 = e277!= 0. ? e237/e277:0.;
188 float rateta33 = e277!= 0. ? e233/e277:0.;
189 reta37->Fill(rateta37, weight);
190 reta33->Fill(rateta33, weight);
191 reta37vset->Fill(rateta37, egamma.pt()/GeV, weight);
192 reta33vset->Fill(rateta33, egamma.pt()/GeV, weight);
193 reta37vseta->Fill(rateta37, egamma.eta(), weight);
194 reta33vseta->Fill(rateta33, egamma.eta(), weight);
195 }
196
197 float shweta1(0);
198 float shweta2(0);
199 if(egamma.showerShapeValue(shweta1, xAOD::EgammaParameters::weta1 )&&
200 egamma.showerShapeValue(shweta2, xAOD::EgammaParameters::weta2 )){
201 weta1->Fill(shweta1, weight);
202 weta2->Fill(shweta2, weight);
203 weta1vset->Fill(shweta1, egamma.pt()/GeV, weight);
204 weta2vset->Fill(shweta2, egamma.pt()/GeV, weight);
205 weta1vseta->Fill(shweta1, egamma.eta(), weight);
206 weta2vseta->Fill(shweta2, egamma.eta(), weight);
207 }
208
209 float emax(0);
210 float emin(0);
211 float emax2(0);
212 if(egamma.showerShapeValue(emax, xAOD::EgammaParameters::emaxs1 )&&
213 egamma.showerShapeValue(emin, xAOD::EgammaParameters::emins1 )&&
214 egamma.showerShapeValue(emax2, xAOD::EgammaParameters::e2tsts1 )){
215
216 float deltemax1 = fabs (emax + emax2) > 0. ? (emax - emax2) / (emax + emax2) : 0.;
217 demax1->Fill(deltemax1, weight);
218 demax2->Fill(emax2 / (1. + 0.009 * et37/GeV)/GeV, weight);
219 de->Fill( (emax2 - emin)/GeV, weight);
220
221 demax1vset->Fill(deltemax1, egamma.pt()/GeV, weight);
222 demax2vset->Fill(emax2 / (1. + 0.009 * et37/GeV)/GeV, egamma.pt()/GeV, weight);
223 devset->Fill( (emax2 - emin)/GeV, egamma.pt()/GeV, weight);
224 demax1vseta->Fill(deltemax1, egamma.eta(), weight);
225 demax2vseta->Fill(emax2 / (1. + 0.009 * et37/GeV)/GeV, egamma.eta(), weight);
226 devseta->Fill( (emax2 - emin)/GeV, egamma.eta(), weight);
227 }
228
229 float shfracs1(0);
230 float shwtots1(0);
231 float fracf1(0);
232 float fracf3(0);
233 float eRatio(0);
234 float difpos7(0);
235 float rhad(0);
236 float rhad1(0);
237 float reta(0);
238 float rphi(0);
239
240 if(egamma.showerShapeValue(shfracs1, xAOD::EgammaParameters::fracs1 )){
241 fracs1->Fill(shfracs1, weight);
242 fracs1vset->Fill(shfracs1, egamma.pt()/GeV, weight);
243 fracs1vseta->Fill(shfracs1, egamma.eta(), weight);
244 }
245
246 if(egamma.showerShapeValue(shwtots1, xAOD::EgammaParameters::wtots1 )){
247 wtots1->Fill(shwtots1, weight);
248 wtots1vset->Fill(shwtots1, egamma.pt()/GeV, weight);
249 wtots1vseta->Fill(shwtots1, egamma.eta(), weight);
250 }
251
252 if(egamma.showerShapeValue(fracf1, xAOD::EgammaParameters::f1 )){
253 f1->Fill(fracf1, weight);
254 f1vset->Fill(fracf1, egamma.pt()/GeV, weight);
255 f1vseta->Fill(fracf1, egamma.eta(), weight);
256 }
257
258 if(egamma.showerShapeValue(fracf3, xAOD::EgammaParameters::f3 )){
259 f3->Fill(fracf3, weight);
260 f3vset->Fill(fracf3, egamma.pt()/GeV, weight);
261 f3vseta->Fill(fracf3, egamma.eta(), weight);
262 }
263 if(egamma.showerShapeValue(eRatio, xAOD::EgammaParameters::Eratio )){
264 Eratio->Fill(eRatio, weight);
265 Eratiovset->Fill(eRatio, egamma.pt()/GeV, weight);
266 Eratiovseta->Fill(eRatio, egamma.eta(), weight);
267 }
268
269 if(egamma.showerShapeValue(difpos7, xAOD::EgammaParameters:: pos7)){
270 pos7->Fill(difpos7, weight);
271 pos7vset->Fill(difpos7, egamma.pt()/GeV, weight);
272 pos7vseta->Fill(difpos7, egamma.eta(), weight);
273 }
274
275 if(egamma.showerShapeValue(rhad, xAOD::EgammaParameters::Rhad )){
276 Rhad->Fill(rhad, weight);
277 Rhadvset->Fill(rhad, egamma.pt()/GeV, weight);
278 Rhadvseta->Fill(rhad, egamma.eta(), weight);
279 }
280
281 if(egamma.showerShapeValue(rhad1, xAOD::EgammaParameters::Rhad1 )){
282 Rhad1->Fill(rhad1, weight);
283 Rhad1vset->Fill(rhad1, egamma.pt()/GeV, weight);
284 Rhad1vseta->Fill(rhad1, egamma.eta(), weight);
285 }
286
287 if(egamma.showerShapeValue(reta, xAOD::EgammaParameters::Reta )){
288 Reta->Fill(reta, weight);
289 Retavset->Fill(reta, egamma.pt()/GeV, weight);
290 Retavseta->Fill(reta, egamma.eta(), weight);
291 }
292
293 if(egamma.showerShapeValue(rphi, xAOD::EgammaParameters::Rphi )){
294 Rphi->Fill(rphi, weight);
295 Rphivset->Fill(rphi, egamma.pt()/GeV, weight);
296 Rphivseta->Fill(rphi, egamma.eta(), weight);
297 }
298
299 float topoEtconeIso(0);
300 float shetcone20(0);
301 if(egamma.isolationValue(shetcone20, xAOD::Iso::topoetcone20)){ //rel20
302// if(egamma.isolationValue(shetcone20, xAOD::EgammaParameters::topoetcone20)){//rel19
303 topoEtconeIso = et37!=0. ? shetcone20/et37 : 0.;
304 clusiso->Fill(topoEtconeIso, weight);
305 clusisovset->Fill(topoEtconeIso, egamma.pt()/GeV, weight);
306 clusisovseta->Fill(topoEtconeIso, egamma.eta(), weight);
307 }
308
309 if(egamma.caloCluster()){
310 const double e0 = egamma.caloCluster()->energyBE(0);
311 const double e1 = egamma.caloCluster()->energyBE(1);
312 const double e2 = egamma.caloCluster()->energyBE(2);
313 const double e3 = egamma.caloCluster()->energyBE(3);
314
315 E0->Fill(e0 / GeV, weight);
316 E1->Fill(e1 / GeV, weight);
317 E2->Fill(e2 / GeV, weight);
318 E3->Fill(e3 / GeV, weight);
319 }
320 }
321
322}
ShowerShapesPlots(PlotBase *pParent, const std::string &sDir, std::string sParticleType)
void fill(const xAOD::Egamma &egamma, const xAOD::EventInfo &eventInfo)
virtual double e() const
energy
virtual double pt() const
transverse momentum
virtual double eta() const
pseudo rapidity
TH1D * Book1D(const std::string &name, const std::string &labels, int nBins, float start, float end, bool prependDir=true)
Book a TH1D histogram.
Definition PlotBase.cxx:94
PlotBase(PlotBase *parent, const std::string &sDir)
Definition PlotBase.cxx:29
TH2F * Book2D(const std::string &name, const std::string &labels, int nBinsX, float startX, float endX, int nBinsY, float startY, float endY, bool prependDir=true)
Book a TH2F histogram.
Definition PlotBase.cxx:123
elec/gamma data class.
Definition egamma.h:58
float beamSpotWeight() const
Weight for beam spot size reweighting.
STL namespace.
@ wtots1
shower width is determined in a window detaxdphi = 0,0625 ×~0,2, corresponding typically to 20 strips...
@ pos7
Difference between the track and the shower positions: sum_{i=i_m-7}^{i=i_m+7}E_i x (i-i_m) / sum_{i=...
@ ethad1
transverse energy in the first sampling of the hadronic calorimeters behind the cluster calculated fr...
Definition EgammaEnums.h:43
@ e277
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 7x7
Definition EgammaEnums.h:81
@ e237
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 3x7
Definition EgammaEnums.h:78
@ ethad
ET leakage into hadronic calorimeter with exclusion of energy in CaloSampling::TileGap3.
Definition EgammaEnums.h:46
@ f3
fraction of energy reconstructed in 3rd sampling
Definition EgammaEnums.h:55
@ f1
E1/E = fraction of energy reconstructed in the first sampling, where E1 is energy in all strips belon...
Definition EgammaEnums.h:53
@ Eratio
(emaxs1-e2tsts1)/(emaxs1+e2tsts1)
@ emaxs1
energy of strip with maximal energy deposit
@ e233
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 3x3 (in cell units e...
Definition EgammaEnums.h:69
@ fracs1
shower shape in the shower core : [E(+/-3)-E(+/-1)]/E(+/-1), where E(+/-n) is the energy in ± n strip...
@ weta2
the lateral width is calculated with a window of 3x5 cells using the energy weighted sum over all cel...
@ weta1
shower width using +/-3 strips around the one with the maximal energy deposit: w3 strips = sqrt{sum(E...
Definition EgammaEnums.h:98
@ e2tsts1
energy of the cell corresponding to second energy maximum in the first sampling
@ emins1
energy reconstructed in the strip with the minimal value between the first and second maximum
@ topoetcone20
Topo-cluster ET-sum.
EventInfo_v1 EventInfo
Definition of the latest event info version.
Egamma_v1 Egamma
Definition of the current "egamma version".
Definition Egamma.h:17