ATLAS Offline Software
Loading...
Searching...
No Matches
EMShower.h
Go to the documentation of this file.
1/*
2 Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
3*/
4
5#ifndef EGAMMAEVENT_EMSHOWER_H
6#define EGAMMAEVENT_EMSHOWER_H
16
19#include <vector>
20//#include <pair>
21
22class EMShower : public egDetail
23
24{
25 public:
26
28 EMShower();
30 virtual ~EMShower() = default;
32 virtual const std::string& className() const;
34 virtual double parameter(egammaParameters::ParamDef) const ;
36 virtual bool hasParameter(egammaParameters::ParamDef) const ;
37
39 void set_parameter(egammaParameters::ParamDef, double, bool overwrite=false) ;
40
41
42 // return values of shower parameters
44 double e011() const;
46 double e033() const;
48 double e132() const;
50 double e1152() const;
53 double ethad1() const;
55 double ethad() const;
57 double ehad1() const;
59 double f1() const;
61 double f3() const;
63 double f1core() const;
65 double f3core() const;
67 double e233() const;
69 double e235() const;
71 double e255() const;
73 double e237() const;
75 double e277() const;
77 double e333() const;
79 double e335() const;
81 double e337() const;
83 double e377() const;
86 double weta1() const;
89 double weta2() const;
91 double e2ts1() const;
93 double e2tsts1() const;
95 double fracs1() const;
98 double widths1() const;
100 double widths2() const;
102 double poscs1() const;
104 double poscs2() const;
106 double asy1() const;
109 double pos() const;
116 double pos7() const;
118 double barys1() const;
120 double wtots1() const;
122 double emins1() const;
124 double emaxs1() const;
126 double etcone() const;
128 double etcone15() const;
130 double etcone20() const;
132 double etcone25() const;
134 double etcone30() const;
136 double etcone35() const;
138 double etcone40() const;
139
141 double etconoisedR02SigAbs2() const;
143 double etconoisedR02SigAbs3() const;
144
146 double etconoisedR03SigAbs2() const;
148 double etconoisedR03SigAbs3() const;
149
151 double etconoisedR04SigAbs2() const;
153 double etconoisedR04SigAbs3() const;
154
156 double etconoisedR04Sig2() const;
158 double etconoisedR04Sig3() const;
159
162 double r33over37allcalo() const;
165 double ecore() const;
167 double zvertex() const;
169 double errz() const;
171 double etap() const;
173 double depth() const;
174
181
188
195
196
199 double etringnoisedR03Sig2() const;
202 double etringnoisedR03Sig3() const;
205 double etringnoisedR03Sig4() const;
206
209 double etringnoisedR04Sig2() const;
212 double etringnoisedR04Sig3() const;
213
215 int linkIndex() const;
218 double ptcone20 () const;
220 double ptcone30 () const;
222 double ptcone40 () const;
223
225 double nucone20 () const;
227 double nucone30 () const;
229 double nucone40 () const;
230
232 double ptcone20_zpv05() const;
234 double ptcone30_zpv05() const;
236 double ptcone40_zpv05() const;
238 double nucone20_zpv05() const;
240 double nucone30_zpv05() const;
242 double nucone40_zpv05() const;
243
245 double etcone15_ptcorrected() const;
247 double etcone20_ptcorrected() const;
249 double etcone25_ptcorrected() const;
251 double etcone30_ptcorrected() const;
253 double etcone35_ptcorrected() const;
255 double etcone40_ptcorrected() const;
256
258 double etcone20_corrected() const;
260 double etcone30_corrected() const;
262 double etcone40_corrected() const;
263
264 //et cone for topo cluster
266 double topoetcone20() const;
268 double topoetcone30() const;
270 double topoetcone40() const;
275
276
277
278
279 // set the values of the different shower shape and isolation variables
280 void set_e011(double);
281 void set_e033(double);
282 void set_e132(double);
283 void set_e1152(double);
284 void set_ethad1(double);
285 void set_ethad(double);
286 void set_ehad1(double);
287 void set_f1(double);
288 void set_f3(double);
289 void set_f1core(double);
290 void set_f3core(double);
291 void set_e233(double);
292 void set_e235(double);
293 void set_e255(double);
294 void set_e237(double);
295 void set_e277(double);
296 void set_e333(double);
297 void set_e335(double);
298 void set_e337(double);
299 void set_e377(double);
300 void set_weta1(double);
301 void set_weta2(double);
302 void set_e2ts1(double);
303 void set_e2tsts1(double);
304 void set_fracs1(double);
305 void set_widths1(double);
306 void set_widths2(double);
307 void set_poscs1(double);
308 void set_poscs2(double);
309 void set_asy1(double);
310 void set_pos(double);
311 void set_pos7(double);
312 void set_barys1(double);
313 void set_wtots1(double);
314 void set_emins1(double);
315 void set_emaxs1(double);
317 void set_ecore(double);
318 void set_zvertex(double);
319 void set_errz(double);
320 void set_etap(double);
321 void set_depth(double);
322 void set_etcone(double);
323 void set_etcone15(double);
324 void set_etcone20(double);
325 void set_etcone25(double);
326 void set_etcone30(double);
327 void set_etcone35(double);
328 void set_etcone40(double);
331
338
342
345
352
353 void set_linkIndex(int);
354 void set_ptcone20(double);
355 void set_ptcone30(double);
356 void set_ptcone40(double);
357 void set_nucone20(double);
358 void set_nucone30(double);
359 void set_nucone40(double);
360 void set_ptcone20_zpv05(double);
361 void set_ptcone30_zpv05(double);
362 void set_ptcone40_zpv05(double);
363 void set_nucone20_zpv05(double);
364 void set_nucone30_zpv05(double);
365 void set_nucone40_zpv05(double);
372
376
377 void set_topoetcone20(double);
378 void set_topoetcone30(double);
379 void set_topoetcone40(double);
382
383 void print() const;
384
385//
386 private:
387
388
389 std::vector < std::pair<egammaParameters::ParamDef , double > > m_parameters;
390 static const std::string s_className;
391};
392
394
395#endif
396
397
398
399
400
401
402
virtual ~EMShower()=default
Destructor.
double e337() const
uncalibrated energy (sum of cells) of the third sampling in a rectangle of size 3x7
void set_etcone25(double)
void set_topoetcone20(double)
void set_etringnoisedR03SigAbs2(double)
void set_etringnoisedR04Sig2(double)
double e255() const
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 5x5
double ptcone30_zpv05() const
summed pt of tracks in a cone with half-opening angle 0.2, |zpv| < 5mm
double nucone40() const
number of tracks in a cone with half-opening angle 0.4
double f1core() const
E1(3x1)/E = fraction of the energy reconstructed in the first longitudinal compartment of the electro...
void set_e011(double)
double etconoisedR02SigAbs2() const
ET in a ring of dR<20, cell abs(ET) above 2 sigma from noise.
EMShower()
Default constructor.
Definition EMShower.cxx:58
void set_etcone40(double)
double weta1() const
shower width using +/-3 strips around the one with the maximal energy deposit: w3 strips = sqrt{sum(E...
double e335() const
uncalibrated energy (sum of cells) of the third sampling in a rectangle of size 3x5
double ehad1() const
E leakage into 1st sampling of had calo (CaloSampling::HEC0 + CaloSampling::TileBar0 + CaloSampling::...
void set_etringnoisedR02SigAbs3(double)
double etconoisedR03SigAbs2() const
ET in a ring of dR<30, cell abs(ET) above 2 sigma from noise.
double widths1() const
same as egammaParameters::weta1 but without corrections on particle impact point inside the cell
void set_barys1(double)
void set_ptcone20_zpv05(double)
void set_etcone30(double)
void set_etconoisedR04SigAbs3(double)
double e011() const
uncalibrated energy (sum of cells) in presampler in a 1x1 window in cells in eta X phi
double etcone40_corrected() const
pt + ED corrected ET in a cone with half-opening angle 0.4
void set_e2tsts1(double)
void set_etconoisedR04SigAbs2(double)
void set_f3(double)
double etcone25_ptcorrected() const
pt corrected ET in a cone with half-opening angle 0.25
void set_etcone20_corrected(double)
double ethad() const
ET leakage into hadronic calorimeter with exclusion of energy in CaloSampling::TileGap3.
void set_etcone35(double)
void set_ethad1(double)
double etringnoisedR04SigAbs3() const
all energy in a ring : (0.1<DR<0.4) with abs(energy) > 3 sigma above noise
double etcone30() const
ET in a cone with half-opening angle 0.3, with exclusion of a window of size 7x5 in electromagnetic c...
void set_ehad1(double)
double etcone40_ptcorrected() const
ptcorrected ET in a cone with half-opening angle 0.4
void set_ecore(double)
double asy1() const
uncorr asymmetry in 3 strips in the 1st sampling
void set_widths1(double)
void set_etringnoisedR04SigAbs3(double)
double ethad1() const
transverse energy in the first sampling of the hadronic calorimeters behind the cluster calculated fr...
void set_f3core(double)
double etcone20() const
ET in a cone with half-opening angle 0.2, with exclusion of a window of size 7x5 in electromagnetic c...
void set_etconoisedR03SigAbs3(double)
double pos7() const
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=...
void set_wtots1(double)
void set_fracs1(double)
double e237() const
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 3x7
double e235() const
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 3x5
double nucone30_zpv05() const
number of tracks in a cone with half-opening angle 0.2, |zpv| < 5mm
double ptcone20() const
summed pt of tracks in a cone with half-opening angle 0.2 (no zvx cut photons, 1mm electrons)
void set_etconoisedR02SigAbs3(double)
void set_etconoisedR04Sig2(double)
void set_asy1(double)
void set_etap(double)
double etringnoisedR03Sig2() const
all energy in a ring : (0.1<DR<0.3) with energy > 2 sigma above noise
void set_etcone30_corrected(double)
void set_e255(double)
double e277() const
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 7x7
double etconoisedR04Sig3() const
ET in a ring of DR<0.40 above 3 sigma of noise.
double etringnoisedR04SigAbs2() const
all energy in a ring : (0.1<DR<0.4) with abs(energy) > 2 sigma above noise
void set_etringnoisedR03Sig2(double)
void set_ptcone30(double)
virtual bool hasParameter(egammaParameters::ParamDef) const
Definition EMShower.cxx:227
static const std::string s_className
interface methods
Definition EMShower.h:390
double e333() const
uncalibrated energy (sum of cells) of the third sampling in a rectangle of size 3x3
double f3core() const
E3(3x3)/E fraction of the energy reconstructed in the third compartment of the electromagnetic calori...
void set_e333(double)
double nucone20() const
number of tracks in a cone with half-opening angle 0.2
void set_linkIndex(int)
double fracs1() const
shower shape in the shower core : [E(+/-3)-E(+/-1)]/E(+/-1), where E(+/-n) is the energy in +- n stri...
double weta2() const
the lateral width is calculated with a window of 3x5 cells using the energy weighted sum over all cel...
double ptcone40_zpv05() const
summed pt of tracks in a cone with half-opening angle 0.2, |zpv| < 5mm
void print() const
Definition EMShower.cxx:118
void set_etcone20_ptcorrected(double)
double etcone20_corrected() const
pt + ED corrected ET in a cone with half-opening angle 0.2
double etringnoisedR03Sig3() const
all energy in a ring : (0.1<DR<0.3) with energy > 3 sigma above noise
double e2tsts1() const
energy of the cell corresponding to second energy maximum in the first sampling
double etcone() const
ET in a cone of R=0.45 in calo, with exclusion of a window of size 7x5 in electromagnetic calorimeter...
void set_f1(double)
void set_etconoisedR02SigAbs2(double)
void set_nucone20(double)
void set_ptcone40_zpv05(double)
double topoetcone40() const
ET in a cone with half-opening angle 0.40.
double etap() const
pointing eta reconstructed from the cluster (first and second sampling)
int linkIndex() const
link index
void set_ethad(double)
double e132() const
uncalibrated energy (sum of cells) in strips in a 3x2 window in cells in eta X phi
void set_ptcone20(double)
double nucone20_zpv05() const
number of tracks in a cone with half-opening angle 0.2, |zpv| < 5mm
void set_r33over37allcalo(double)
void set_e335(double)
double nucone30() const
number of tracks in a cone with half-opening angle 0.3
void set_etringnoisedR03Sig4(double)
void set_emins1(double)
double zvertex() const
pointing z at vertex reconstructed from the cluster
void set_etcone30_ptcorrected(double)
double nucone40_zpv05() const
number of tracks in a cone with half-opening angle 0.2, |zpv| < 5mm
void set_depth(double)
double emins1() const
energy reconstructed in the strip with the minimal value between the first and second maximum
double etringnoisedR03Sig4() const
all energy in a ring : (0.1<DR<0.3) with energy > 4 sigma above noise
double e233() const
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 3x3 (in cell units e...
double topoetcone40_ptcorrected() const
pt corrected ET in a cone with half-opening angle 0.40
void set_weta2(double)
void set_etringnoisedR03SigAbs3(double)
double ptcone20_zpv05() const
summed pt of tracks in a cone with half-opening angle 0.2, |zpv| < 5mm
void set_e237(double)
double widths2() const
same as egammaParameters::weta2 but without corrections on particle impact point inside the cell
void set_etconoisedR03SigAbs2(double)
void set_etcone15(double)
void set_poscs2(double)
void set_zvertex(double)
double topoetcone30() const
ET in a cone with half-opening angle 0.30.
double etringnoisedR03SigAbs3() const
all energy in a ring : (0.1<DR<0.3) with abs(energy) > 3 sigma above noise
void set_topoetcone40_ptcorrected(double)
double etringnoisedR04Sig3() const
all energy in a ring : (0.1<DR<0.4) with energy > 3 sigma above noise
double topoetcone40_corrected() const
fully corrected ET in a cone with half-opening angle 0.40
void set_pos(double)
double etringnoisedR03SigAbs2() const
all energy in a ring : (0.1<DR<0.3) with abs(energy) > 2 sigma above noise
void set_topoetcone40_corrected(double)
double errz() const
error associated to zvertex
double poscs1() const
relative position in eta within cell in 1st sampling
void set_etringnoisedR02SigAbs2(double)
double e2ts1() const
2nd max in strips calc by summing 3 strips
double poscs2() const
relative position in eta within cell in 2nd sampling
void set_e377(double)
double etconoisedR03SigAbs3() const
ET in a ring of dR<30 cell abs (ET) above 3 sigma from noise.
void set_e337(double)
void set_etconoisedR04Sig3(double)
double ecore() const
core energy in em calo E(core) = E0(3x3) + E1(15x2) + E2(5x5) + E3(3x5)
double etcone20_ptcorrected() const
pt corrected ET in a cone with half-opening angle 0.2
void set_poscs1(double)
void set_etcone15_ptcorrected(double)
double etringnoisedR02SigAbs3() const
all energy in a ring : (0.1<DR<0.2) with abs(energy) > 3 sigma above noise
void set_etcone25_ptcorrected(double)
void set_etcone(double)
void set_etringnoisedR04Sig3(double)
void set_emaxs1(double)
void set_e033(double)
void set_nucone40(double)
virtual double parameter(egammaParameters::ParamDef) const
Definition EMShower.cxx:71
double etringnoisedR04Sig2() const
all energy in a ring : (0.1<DR<0.4) with energy > 2 sigma above noise
double etcone30_corrected() const
pt + ED corrected ET in a cone with half-opening angle 0.3
void set_ptcone30_zpv05(double)
void set_pos7(double)
double emaxs1() const
energy of strip with maximal energy deposit
void set_e235(double)
double barys1() const
barycentre in sampling 1 calculated in 3 strips
double r33over37allcalo() const
1-ratio of energy in 3x3 over 3x7 cells; E(3x3) = E0(1x1) + E1(3x1) + E2(3x3) + E3(3x3); E(3x7) = E0(...
double etconoisedR02SigAbs3() const
ET in a ring of dR<20 cell abs (ET) above 3 sigma from noise.
void set_parameter(egammaParameters::ParamDef, double, bool overwrite=false)
Definition EMShower.cxx:90
std::vector< std::pair< egammaParameters::ParamDef, double > > m_parameters
Definition EMShower.h:389
void set_e2ts1(double)
void set_widths2(double)
void set_e233(double)
void set_topoetcone30(double)
double etcone15() const
ET in a cone with half-opening angle 0.15, with exclusion of a window of size 7x5 in electromagnetic ...
double e1152() const
uncalibrated energy (sum of cells) in strips in a 15x2 window in cells in eta X phi
void set_errz(double)
double ptcone30() const
summed pt of tracks in a cone with half-opening angle 0.3
void set_topoetcone40(double)
double etcone15_ptcorrected() const
pt corrected ET in a cone with half-opening angle 0.15
void set_weta1(double)
void set_nucone20_zpv05(double)
void set_nucone30_zpv05(double)
void set_nucone40_zpv05(double)
double depth() const
pointing depth of the shower as calculated in egammaqgcld
double etcone35() const
ET in a cone with half-opening angle 0.35, with exclusion of a window of size 7x5 in electromagnetic ...
void set_nucone30(double)
double topoetcone20() const
ET in a cone with half-opening angle 0.20.
double etcone30_ptcorrected() const
pt corrected ET in a cone with half-opening angle 0.3
double etconoisedR04Sig2() const
ET in a ring of DR<0.40 above 2 sigma of noise.
void set_etcone40_corrected(double)
void set_etcone40_ptcorrected(double)
void set_etcone20(double)
double etringnoisedR02SigAbs2() const
all energy in a ring : (0.1<DR<0.2) with abs(energy) > 2 sigma above noise
double pos() const
difference between shower cell and predicted track in +/- 1 cells
double etconoisedR04SigAbs2() const
ET in a ring of dR<40, cell abs(ET) above 2 sigma from noise.
double etconoisedR04SigAbs3() const
ET in a ring of dR<40 cell abs (ET) above 3 sigma from noise.
double etcone40() const
ET in a cone with half-opening angle 0.4, with exclusion of a window of size 7x5 in electromagnetic c...
void set_e277(double)
double wtots1() const
shower width is determined in a window detaxdphi = 0,0625 x~0,2, corresponding typically to 20 strips...
void set_f1core(double)
void set_etringnoisedR03Sig3(double)
double e377() const
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 7x7
double etcone25() const
ET in a cone with half-opening angle 0.25, with exclusion of a window of size 7x5 in electromagnetic ...
void set_etcone35_ptcorrected(double)
void set_ptcone40(double)
double etcone35_ptcorrected() const
pt corrected ET in a cone with half-opening angle 0.35
void set_etringnoisedR04SigAbs2(double)
void set_e1152(double)
double e033() const
uncalibrated energy (sum of cells) in presampler in a 3x3 window in cells in eta X phi
virtual const std::string & className() const
Definition EMShower.cxx:65
void set_e132(double)
double f3() const
fraction of energy reconstructed in 3rd sampling
double f1() const
E1/E = fraction of energy reconstructed in the first sampling, where E1 is energy in all strips belon...
double ptcone40() const
summed pt of tracks in a cone with half-opening angle 0.4
egDetail()=default
Default constructor.