ATLAS Offline Software
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 
9 using CLHEP::GeV;
10 
11 namespace 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 
83  void ShowerShapesPlots::initializePlots(){
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 
154  void ShowerShapesPlots::fill(const xAOD::Egamma& egamma, const xAOD::EventInfo& eventInfo) {
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, weight);
316  E1->Fill(e1, weight);
317  E2->Fill(e2, weight);
318  E3->Fill(e3, weight);
319  }
320  }
321 
322 }
xAOD::EgammaParameters::e233
@ e233
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 3x3 (in cell units e...
Definition: EgammaEnums.h:68
Egamma::ShowerShapesPlots::wtots1
TH1 * wtots1
Definition: ShowerShapesPlots.h:33
Egamma::ShowerShapesPlots::f3vseta
TH2 * f3vseta
Definition: ShowerShapesPlots.h:81
Egamma::ShowerShapesPlots::Rphivseta
TH2 * Rphivseta
Definition: ShowerShapesPlots.h:88
GeV
#define GeV
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:17
xAOD::Iso::topoetcone20
@ topoetcone20
Topo-cluster ET-sum.
Definition: IsolationType.h:48
ParticleImpl::pt
virtual double pt() const
transverse momentum
Definition: ParticleImpl.h:554
ParticleGun_SamplingFraction.eta2
eta2
Definition: ParticleGun_SamplingFraction.py:96
Egamma::ShowerShapesPlots::fill
void fill(const xAOD::Egamma &egamma, const xAOD::EventInfo &eventInfo)
Definition: ShowerShapesPlots.cxx:154
PlotBase
Definition: PlotBase.h:34
Egamma::ShowerShapesPlots::fracs1vset
TH2 * fracs1vset
Definition: ShowerShapesPlots.h:57
Egamma::ShowerShapesPlots::Rhad1
TH1 * Rhad1
Definition: ShowerShapesPlots.h:44
Egamma::ShowerShapesPlots::f3vset
TH2 * f3vset
Definition: ShowerShapesPlots.h:60
xAOD::EgammaParameters::Reta
@ Reta
e237/e277
Definition: EgammaEnums.h:154
egammaEnergyPositionAllSamples::e1
double e1(const xAOD::CaloCluster &cluster)
return the uncorrected cluster energy in 1st sampling
Egamma::ShowerShapesPlots::Retavset
TH2 * Retavset
Definition: ShowerShapesPlots.h:66
Egamma
Definition: ClusMomentumPlots.cxx:11
python.copyTCTOutput.sDir
sDir
Definition: copyTCTOutput.py:60
Egamma::ShowerShapesPlots::Rhad1vseta
TH2 * Rhad1vseta
Definition: ShowerShapesPlots.h:86
Egamma::ShowerShapesPlots::fracs1vseta
TH2 * fracs1vseta
Definition: ShowerShapesPlots.h:78
Egamma::ShowerShapesPlots::m_sParticleType
std::string m_sParticleType
Definition: ShowerShapesPlots.h:21
Egamma::ShowerShapesPlots::devseta
TH2 * devseta
Definition: ShowerShapesPlots.h:77
Egamma::ShowerShapesPlots::demax1vset
TH2 * demax1vset
Definition: ShowerShapesPlots.h:54
Egamma::ShowerShapesPlots::middleevseta
TH2 * middleevseta
Definition: ShowerShapesPlots.h:70
PlotBase::Book2D
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
ShowerShapesPlots.h
Egamma::ShowerShapesPlots::hadleak
TH1 * hadleak
Definition: ShowerShapesPlots.h:23
xAOD::EgammaParameters::ethad1
@ ethad1
transverse energy in the first sampling of the hadronic calorimeters behind the cluster calculated fr...
Definition: EgammaEnums.h:42
Egamma::ShowerShapesPlots::pos7
TH1 * pos7
Definition: ShowerShapesPlots.h:36
Egamma::ShowerShapesPlots::E2
TH1 * E2
Definition: ShowerShapesPlots.h:41
xAOD::Egamma_v1
Definition: Egamma_v1.h:56
Egamma::ShowerShapesPlots::weta2vset
TH2 * weta2vset
Definition: ShowerShapesPlots.h:53
Egamma::ShowerShapesPlots::wtots1vset
TH2 * wtots1vset
Definition: ShowerShapesPlots.h:58
Egamma::ShowerShapesPlots::middlee
TH1 * middlee
Definition: ShowerShapesPlots.h:24
Egamma::ShowerShapesPlots::clusiso
TH1 * clusiso
Definition: ShowerShapesPlots.h:37
xAOD::EgammaParameters::Rphi
@ Rphi
e233/e237
Definition: EgammaEnums.h:156
Egamma::ShowerShapesPlots::weta2vseta
TH2 * weta2vseta
Definition: ShowerShapesPlots.h:74
Egamma::ShowerShapesPlots::de
TH1 * de
Definition: ShowerShapesPlots.h:31
Egamma::ShowerShapesPlots::wtots1vseta
TH2 * wtots1vseta
Definition: ShowerShapesPlots.h:79
Egamma::ShowerShapesPlots::pos7vseta
TH2 * pos7vseta
Definition: ShowerShapesPlots.h:82
Egamma::ShowerShapesPlots::Rhadvseta
TH2 * Rhadvseta
Definition: ShowerShapesPlots.h:85
xAOD::EgammaParameters::wtots1
@ wtots1
shower width is determined in a window detaxdphi = 0,0625 ×~0,2, corresponding typically to 20 strips...
Definition: EgammaEnums.h:140
Egamma::ShowerShapesPlots::weta2
TH1 * weta2
Definition: ShowerShapesPlots.h:28
egamma
Definition: egamma.h:58
Egamma::ShowerShapesPlots::initializePlots
virtual void initializePlots()
Definition: ShowerShapesPlots.cxx:83
Egamma::ShowerShapesPlots::reta33vseta
TH2 * reta33vseta
Definition: ShowerShapesPlots.h:72
dqt_zlumi_pandas.weight
int weight
Definition: dqt_zlumi_pandas.py:189
PlotBase::Book1D
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
Egamma::ShowerShapesPlots::weta1vseta
TH2 * weta1vseta
Definition: ShowerShapesPlots.h:73
Egamma::ShowerShapesPlots::Rhad1vset
TH2 * Rhad1vset
Definition: ShowerShapesPlots.h:65
Egamma::ShowerShapesPlots::Rphi
TH1 * Rphi
Definition: ShowerShapesPlots.h:46
Egamma::ShowerShapesPlots::fracs1
TH1 * fracs1
Definition: ShowerShapesPlots.h:32
ParticleImpl::e
virtual double e() const
energy
Definition: ParticleImpl.h:534
Egamma::ShowerShapesPlots::weta1
TH1 * weta1
Definition: ShowerShapesPlots.h:27
Egamma::ShowerShapesPlots::Rphivset
TH2 * Rphivset
Definition: ShowerShapesPlots.h:67
xAOD::EgammaParameters::f3
@ f3
fraction of energy reconstructed in 3rd sampling
Definition: EgammaEnums.h:54
Egamma::ShowerShapesPlots::f1vseta
TH2 * f1vseta
Definition: ShowerShapesPlots.h:80
xAOD::EgammaParameters::ethad
@ ethad
ET leakage into hadronic calorimeter with exclusion of energy in CaloSampling::TileGap3.
Definition: EgammaEnums.h:45
Egamma::ShowerShapesPlots::f3
TH1 * f3
Definition: ShowerShapesPlots.h:35
Egamma::ShowerShapesPlots::demax2vseta
TH2 * demax2vseta
Definition: ShowerShapesPlots.h:76
CheckAppliedSFs.e3
e3
Definition: CheckAppliedSFs.py:264
xAOD::EgammaParameters::f1
@ f1
E1/E = fraction of energy reconstructed in the first sampling, where E1 is energy in all strips belon...
Definition: EgammaEnums.h:52
Egamma::ShowerShapesPlots::Eratio
TH1 * Eratio
Definition: ShowerShapesPlots.h:38
Egamma::ShowerShapesPlots::Retavseta
TH2 * Retavseta
Definition: ShowerShapesPlots.h:87
Egamma::ShowerShapesPlots::reta37vset
TH2 * reta37vset
Definition: ShowerShapesPlots.h:50
xAOD::EgammaParameters::pos7
@ 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=...
Definition: EgammaEnums.h:133
xAOD::EgammaParameters::emins1
@ emins1
energy reconstructed in the strip with the minimal value between the first and second maximum
Definition: EgammaEnums.h:143
Egamma::ShowerShapesPlots::f1
TH1 * f1
Definition: ShowerShapesPlots.h:34
Egamma::ShowerShapesPlots::E3
TH1 * E3
Definition: ShowerShapesPlots.h:42
Egamma::ShowerShapesPlots::demax1vseta
TH2 * demax1vseta
Definition: ShowerShapesPlots.h:75
Egamma::ShowerShapesPlots::hadleakvseta
TH2 * hadleakvseta
Definition: ShowerShapesPlots.h:69
Egamma::ShowerShapesPlots::reta37vseta
TH2 * reta37vseta
Definition: ShowerShapesPlots.h:71
Egamma::ShowerShapesPlots::pos7vset
TH2 * pos7vset
Definition: ShowerShapesPlots.h:61
Egamma::ShowerShapesPlots::weta1vset
TH2 * weta1vset
Definition: ShowerShapesPlots.h:52
ParticleImpl::eta
virtual double eta() const
pseudo rapidity
Definition: ParticleImpl.h:514
Egamma::ShowerShapesPlots::demax2vset
TH2 * demax2vset
Definition: ShowerShapesPlots.h:55
xAOD::EgammaParameters::Rhad1
@ Rhad1
ethad1/et
Definition: EgammaEnums.h:162
Egamma::ShowerShapesPlots::E0
TH1 * E0
Definition: ShowerShapesPlots.h:39
Egamma::ShowerShapesPlots::Eratiovseta
TH2 * Eratiovseta
Definition: ShowerShapesPlots.h:84
Egamma::ShowerShapesPlots::reta33
TH1 * reta33
Definition: ShowerShapesPlots.h:26
Egamma::ShowerShapesPlots::ShowerShapesPlots
ShowerShapesPlots(PlotBase *pParent, const std::string &sDir, std::string sParticleType)
Definition: ShowerShapesPlots.cxx:14
Egamma::ShowerShapesPlots::Reta
TH1 * Reta
Definition: ShowerShapesPlots.h:45
Egamma::ShowerShapesPlots::demax2
TH1 * demax2
Definition: ShowerShapesPlots.h:30
xAOD::EventInfo_v1
Class describing the basic event information.
Definition: EventInfo_v1.h:43
Egamma::ShowerShapesPlots::f1vset
TH2 * f1vset
Definition: ShowerShapesPlots.h:59
Egamma::ShowerShapesPlots::clusisovset
TH2 * clusisovset
Definition: ShowerShapesPlots.h:62
Egamma::ShowerShapesPlots::clusisovseta
TH2 * clusisovseta
Definition: ShowerShapesPlots.h:83
egammaEnergyPositionAllSamples::e2
double e2(const xAOD::CaloCluster &cluster)
return the uncorrected cluster energy in 2nd sampling
Egamma::ShowerShapesPlots::hadleakvset
TH2 * hadleakvset
Definition: ShowerShapesPlots.h:48
Egamma::ShowerShapesPlots::Rhadvset
TH2 * Rhadvset
Definition: ShowerShapesPlots.h:64
Egamma::ShowerShapesPlots::E1
TH1 * E1
Definition: ShowerShapesPlots.h:40
xAOD::EgammaParameters::e277
@ e277
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 7x7
Definition: EgammaEnums.h:80
Egamma::ShowerShapesPlots::reta33vset
TH2 * reta33vset
Definition: ShowerShapesPlots.h:51
Egamma::ShowerShapesPlots::Rhad
TH1 * Rhad
Definition: ShowerShapesPlots.h:43
xAOD::EgammaParameters::weta1
@ weta1
shower width using +/-3 strips around the one with the maximal energy deposit: w3 strips = sqrt{sum(E...
Definition: EgammaEnums.h:97
xAOD::EgammaParameters::Eratio
@ Eratio
(emaxs1-e2tsts1)/(emaxs1+e2tsts1)
Definition: EgammaEnums.h:158
xAOD::EgammaParameters::e237
@ e237
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 3x7
Definition: EgammaEnums.h:77
Egamma::ShowerShapesPlots::demax1
TH1 * demax1
Definition: ShowerShapesPlots.h:29
egammaEnergyPositionAllSamples::e0
double e0(const xAOD::CaloCluster &cluster)
return the uncorrected cluster energy in pre-sampler
Egamma::ShowerShapesPlots::middleevset
TH2 * middleevset
Definition: ShowerShapesPlots.h:49
xAOD::EgammaParameters::e2tsts1
@ e2tsts1
energy of the cell corresponding to second energy maximum in the first sampling
Definition: EgammaEnums.h:108
xAOD::EgammaParameters::Rhad
@ Rhad
ethad/et
Definition: EgammaEnums.h:160
Egamma::ShowerShapesPlots::reta37
TH1 * reta37
Definition: ShowerShapesPlots.h:25
Egamma::ShowerShapesPlots::Eratiovset
TH2 * Eratiovset
Definition: ShowerShapesPlots.h:63
Egamma::ShowerShapesPlots::devset
TH2 * devset
Definition: ShowerShapesPlots.h:56
xAOD::EgammaParameters::emaxs1
@ emaxs1
energy of strip with maximal energy deposit
Definition: EgammaEnums.h:145
xAOD::EventInfo_v1::beamSpotWeight
float beamSpotWeight() const
Weight for beam spot size reweighting.
Definition: EventInfo_v1.cxx:970
xAOD::EgammaParameters::fracs1
@ fracs1
shower shape in the shower core : [E(+/-3)-E(+/-1)]/E(+/-1), where E(+/-n) is the energy in ± n strip...
Definition: EgammaEnums.h:111
read_hist_ntuple.f1
f1
Definition: read_hist_ntuple.py:4
xAOD::EgammaParameters::weta2
@ weta2
the lateral width is calculated with a window of 3x5 cells using the energy weighted sum over all cel...
Definition: EgammaEnums.h:103
Generated on Sun Dec 22 2024 21:18:06 for ATLAS Offline Software by   doxygen 1.8.18