ATLAS Offline Software
HIMonitoringPhotonsTool.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
3  */
4 
5 
6 #include "PATCore/AcceptData.h"
7 
8 #include "AthenaMonitoring/AthenaMonManager.h"
9 #include "HIMonitoringPhotonsTool.h"
10 
11 #include "ElectronPhotonSelectorTools/AsgPhotonIsEMSelector.h"
12 #include "ElectronPhotonSelectorTools/egammaPIDdefs.h"
13 
14 #include "xAODHIEvent/HIEventShapeContainer.h"
15 #include "xAODEgamma/PhotonContainer.h"
16 #include "AthContainers/ConstAccessor.h"
17 
18 #include "LWHists/TH1D_LW.h"
19 #include "LWHists/TH2D_LW.h"
20 
21 HIMonitoringPhotonsTool::
22  HIMonitoringPhotonsTool(const std::string& type, const std::string& name,
23  const IInterface* parent) : ManagedMonitorToolBase(type, name, parent) {
24  m_FCalEt = 0;
25 
26  //From Dominik Derendarz - manual cherry-picking of 345b289e
27  declareProperty( "FCalEt_nbins", m_FCalEt_nbins=95);
28  declareProperty( "lowFCalEt", m_FCalEt_low=-0.15);
29  declareProperty( "highFCalEt", m_FCalEt_high=0.8);
30  declareProperty( "FCalEt_peripheral", m_FCalEt_peripheral=20);
31  declareProperty( "FCalEt_central", m_FCalEt_central=50);
32  //In master, before the cherry picking, was:
33  //m_FCalEt_nbins = 50;
34  //m_FCalEt_low = -50;
35  //m_FCalEt_high = 200;
36 
37  m_PhotonPt_nbins = 50;
38  m_PhotonPt_low = 0;
39  m_PhotonPt_high = 500;
40 
41  m_PhotonEta_nbins = 60;
42  m_PhotonEta_low = -3;
43  m_PhotonEta_high = +3;
44 
45  m_PhotonPhi_nbins = 64;
46  m_PhotonPhi_low = -3.14159;
47  m_PhotonPhi_high = +3.14159;
48 
49  m_IsoWide_nbins = 100;
50  m_IsoWide_low = -50.0;
51  m_IsoWide_high = 50.0;
52 
53  m_IsoNarrow_nbins = 60;
54  m_IsoNarrow_low = -30.0;
55  m_IsoNarrow_high = 30.0;
56 
57  m_reta_nbins = 50;
58  m_reta_low = 0.70;
59  m_reta_high = 1.1;
60 
61  m_rphi_nbins = 50;
62  m_rphi_low = 0.60;
63  m_rphi_high = 1.1;
64 
65  m_weta2_nbins = 50;
66  m_weta2_low = 0.005;
67  m_weta2_high = 0.015;
68 
69  m_rhad_nbins = 50;
70  m_rhad_low = -0.3;
71  m_rhad_high = +0.3;
72 
73  m_wtots1_nbins = 50;
74  m_wtots1_low = 0;
75  m_wtots1_high = 6;
76 
77  m_fracs1_nbins = 50;
78  m_fracs1_low = -0.1;
79  m_fracs1_high = 1.2;
80 
81  m_deltae_nbins = 50;
82  m_deltae_low = -0.1;
83  m_deltae_high = 1.5;
84 
85  m_eratio_nbins = 50;
86  m_eratio_low = -0.1;
87  m_eratio_high = 1.1;
88 
89  m_f1_nbins = 50;
90  m_f1_low = -0.05;
91  m_f1_high = 0.60;
92 
93  StatusCode sc;
94  m_photonLooseIsEMSelector = new AsgPhotonIsEMSelector("PhotonLooseIsEMSelector");
95  sc &= m_photonLooseIsEMSelector->setProperty("isEMMask", egammaPID::PhotonLoose);
96  sc &= m_photonLooseIsEMSelector->setProperty("ConfigFile",
97  "ElectronPhotonSelectorTools/offline/mc15_20150712/PhotonIsEMLooseSelectorCutDefs.conf");
98  if (!m_photonLooseIsEMSelector->initialize().isSuccess()) {
99  Fatal("MyFunction", "Failed to initialize PhotonLooseIsEMSelector ");
100  }
101 
102  m_photonTightIsEMSelector = new AsgPhotonIsEMSelector("PhotonTightIsEMSelector");
103  sc &= m_photonTightIsEMSelector->setProperty("isEMMask", egammaPID::PhotonTight);
104  sc &= m_photonTightIsEMSelector->setProperty("ConfigFile",
105  "ElectronPhotonSelectorTools/offline/mc15_20150712/PhotonIsEMTightSelectorCutDefs.conf");
106  if (!m_photonTightIsEMSelector->initialize().isSuccess()) {
107  Fatal("MyFunction", "Failed to initialize PhotonTightIsEMSelector ");
108  }
109  if (sc.isFailure()) {
110  Warning("MyFunction", "Failed to set some properties");
111  }
112 }
113 
114 HIMonitoringPhotonsTool::~HIMonitoringPhotonsTool() {
115 }
116 
117 // Description: Used for rebooking unmanaged histograms
118 StatusCode HIMonitoringPhotonsTool::bookHistogramsRecurrent( ) {
119  return StatusCode::SUCCESS;
120 }
121 
122 // Description: Used for re-booking managed histograms
123 StatusCode HIMonitoringPhotonsTool::bookHistograms( ) {
124  if (m_environment == AthenaMonManager::online) {
125  // book histograms that are only made in the online environment...
126  }
127 
128  if (m_dataType == AthenaMonManager::cosmics) {
129  // book histograms that are only relevant for cosmics data...
130  }
131 
132  book_hist();
133 
134 
135  return StatusCode::SUCCESS;
136 }
137 
138 StatusCode HIMonitoringPhotonsTool::fillHistograms() {
139  StatusCode sc;
140 
142  const xAOD::HIEventShapeContainer* evtShape;
143 
144  sc = evtStore()->retrieve(evtShape, "HIEventShape");
145  if (sc.isFailure()) {
146  ATH_MSG_ERROR("Could not find HIEventShape");
147  return StatusCode::FAILURE;
148  } else {
149  ATH_MSG_INFO("HIEventShape retrieved from StoreGate");
150  }
151 
152  int size = evtShape->size();
153  for (int i = 0; i < size; i++) {
154  const xAOD::HIEventShape* sh = evtShape->at(i);
155  if (sh->layer() == 21 || sh->layer() == 22 || sh->layer() == 23) {
156  m_FCalEt += sh->et();
157  }
158  }
159 
160  m_FCalEt = m_FCalEt * 1e-3;
161 
162 
163  const xAOD::PhotonContainer* photons = 0;
164  sc = evtStore()->retrieve(photons, "Photons");
165  if (sc.isFailure()) {
166  ATH_MSG_ERROR("Could not find Photons");
167  return StatusCode::FAILURE;
168  } else {
169  ATH_MSG_INFO("Photons retrieved from StoreGate");
170  }
171 
172  xAOD::PhotonContainer::const_iterator photons_itr = photons->begin();
173  xAOD::PhotonContainer::const_iterator photons_end = photons->end();
174 
175  for (; photons_itr != photons_end; ++photons_itr) {
176  float pt = (*photons_itr)->pt() * 1e-3;
177  float eta = (*photons_itr)->caloCluster()->etaBE(2);
178  float phi = (*photons_itr)->phi();
179 
180  m_h_photon_pt_eta->Fill(pt, eta);
181  m_h_photon_pt_phi->Fill(pt, phi);
182  m_h_photon_pt_fcal->Fill(pt, phi);
183 
184  bool loose_MC15 = (bool) m_photonLooseIsEMSelector->accept(*photons_itr);
185  bool tight_MC15 = (bool) m_photonTightIsEMSelector->accept(*photons_itr);
186 
187  if (tight_MC15) {
188  m_h_photon_pt_eta_tight->Fill(pt, eta);
189  m_h_photon_pt_phi_tight->Fill(pt, phi);
190  m_h_photon_pt_fcal_tight->Fill(pt, phi);
191  }
192 
193  if (pt < 20) continue;
194 
195  m_h_photon_eta_phi_ptCut->Fill(eta, phi);
196 
197  if (tight_MC15) {
198  m_h_photon_eta_phi_ptCut_tight->Fill(eta, phi);
199  }
200 
201  static const SG::ConstAccessor< float > etcone20Acc("etcone20");
202  static const SG::ConstAccessor< float > etcone30Acc("etcone30");
203  static const SG::ConstAccessor< float > etcone40Acc("etcone40");
204  float etcone20 = etcone20Acc(**photons_itr) * 1e-3;
205  float etcone30 = etcone30Acc(**photons_itr) * 1e-3;
206  float etcone40 = etcone40Acc(**photons_itr) * 1e-3;
207 
208  m_h_photon_fcal_etcone20_ptCut->Fill(m_FCalEt, etcone20);
209  m_h_photon_fcal_etcone30_ptCut->Fill(m_FCalEt, etcone30);
210  m_h_photon_fcal_etcone40_ptCut->Fill(m_FCalEt, etcone40);
211 
212  //In master, before the cherry picking, was:
213  //if (m_FCalEt < 20) {
214  if (m_FCalEt < m_FCalEt_peripheral) {
215  m_h_photon_pt_etcone40_fcal0->Fill(pt, etcone40);
216 
217  m_h_photon_etcone20_ptCut_fcal0->Fill(etcone20);
218  m_h_photon_etcone30_ptCut_fcal0->Fill(etcone30);
219  m_h_photon_etcone40_ptCut_fcal0->Fill(etcone40);
220 
221  m_h_photon_etcone30_ptCut_fcal0_all->Fill(etcone30);
222  if (loose_MC15) m_h_photon_etcone30_ptCut_fcal0_loose->Fill(etcone30);
223  if (tight_MC15) m_h_photon_etcone30_ptCut_fcal0_tight->Fill(etcone30);
224  }
225  //if (m_FCalEt > 20 && m_FCalEt < 50) {
226  if (m_FCalEt > m_FCalEt_peripheral && m_FCalEt < m_FCalEt_central) {
227  m_h_photon_pt_etcone40_fcal1->Fill(pt, etcone40);
228 
229  m_h_photon_etcone20_ptCut_fcal1->Fill(etcone20);
230  m_h_photon_etcone30_ptCut_fcal1->Fill(etcone30);
231  m_h_photon_etcone40_ptCut_fcal1->Fill(etcone40);
232 
233  m_h_photon_etcone30_ptCut_fcal1_all->Fill(etcone30);
234  if (loose_MC15) m_h_photon_etcone30_ptCut_fcal1_loose->Fill(etcone30);
235  if (tight_MC15) m_h_photon_etcone30_ptCut_fcal1_tight->Fill(etcone30);
236  }
237  //if (m_FCalEt > 50) {
238  if (m_FCalEt > m_FCalEt_central) {
239  m_h_photon_pt_etcone40_fcal2->Fill(pt, etcone40);
240 
241  m_h_photon_etcone20_ptCut_fcal2->Fill(etcone20);
242  m_h_photon_etcone30_ptCut_fcal2->Fill(etcone30);
243  m_h_photon_etcone40_ptCut_fcal2->Fill(etcone40);
244 
245  m_h_photon_etcone30_ptCut_fcal2_all->Fill(etcone30);
246  if (loose_MC15) m_h_photon_etcone30_ptCut_fcal2_loose->Fill(etcone30);
247  if (tight_MC15) m_h_photon_etcone30_ptCut_fcal2_tight->Fill(etcone30);
248  }
249 
250  static const SG::ConstAccessor<float> RetaAcc("Reta");
251  static const SG::ConstAccessor<float> RphiAcc("Rphi");
252  static const SG::ConstAccessor< float > weta2Acc("weta2");
253  float reta = RetaAcc(**photons_itr);
254  float rphi = RphiAcc(**photons_itr);
255  float weta2 = weta2Acc(**photons_itr);
256 
257  static const SG::ConstAccessor< float > RhadAcc("Rhad");
258  static const SG::ConstAccessor< float > Rhad1Acc("Rhad1");
259  static const SG::ConstAccessor< float > wtots1Acc("wtots1");
260  static const SG::ConstAccessor< float > fracs1Acc("fracs1");
261  float rhad = RhadAcc(**photons_itr);
262  float rhad1 = Rhad1Acc(**photons_itr);
263  float wtots1 = wtots1Acc(**photons_itr);
264  float fracs1 = fracs1Acc(**photons_itr);
265 
266  static const SG::ConstAccessor< float > DeltaEAcc("DeltaE");
267  static const SG::ConstAccessor< float > EratioAcc("Eratio");
268  static const SG::ConstAccessor< float > f1Acc("f1");
269  float deltae = DeltaEAcc(**photons_itr) * 1e-3;
270  float eratio = EratioAcc(**photons_itr);
271  float f1 = f1Acc(**photons_itr);
272 
273  if (fabs(eta) < 1.37) {
274  m_h_photon_reta_fcal_ptCut_eta0->Fill(m_FCalEt, reta);
275  m_h_photon_rphi_fcal_ptCut_eta0->Fill(m_FCalEt, rphi);
276  m_h_photon_weta2_fcal_ptCut_eta0->Fill(m_FCalEt, weta2);
277  m_h_photon_rhad_fcal_ptCut_eta0->Fill(m_FCalEt, rhad);
278  m_h_photon_wtots1_fcal_ptCut_eta0->Fill(m_FCalEt, wtots1);
279  m_h_photon_fracs1_fcal_ptCut_eta0->Fill(m_FCalEt, fracs1);
280  m_h_photon_deltae_fcal_ptCut_eta0->Fill(m_FCalEt, deltae);
281  m_h_photon_eratio_fcal_ptCut_eta0->Fill(m_FCalEt, eratio);
282  m_h_photon_f1_fcal_ptCut_eta0->Fill(m_FCalEt, f1);
283  }
284  if (fabs(eta) > 1.52) {
285  m_h_photon_reta_fcal_ptCut_eta1->Fill(m_FCalEt, reta);
286  m_h_photon_rphi_fcal_ptCut_eta1->Fill(m_FCalEt, rphi);
287  m_h_photon_weta2_fcal_ptCut_eta1->Fill(m_FCalEt, weta2);
288  m_h_photon_rhad_fcal_ptCut_eta1->Fill(m_FCalEt, rhad1);
289  m_h_photon_wtots1_fcal_ptCut_eta1->Fill(m_FCalEt, wtots1);
290  m_h_photon_fracs1_fcal_ptCut_eta1->Fill(m_FCalEt, fracs1);
291  m_h_photon_deltae_fcal_ptCut_eta1->Fill(m_FCalEt, deltae);
292  m_h_photon_eratio_fcal_ptCut_eta1->Fill(m_FCalEt, eratio);
293  m_h_photon_f1_fcal_ptCut_eta1->Fill(m_FCalEt, f1);
294  }
295 
296  m_h_photon_reta_ptCut_all->Fill(reta);
297  m_h_photon_rphi_ptCut_all->Fill(rphi);
298  m_h_photon_weta2_ptCut_all->Fill(weta2);
299  m_h_photon_rhad_ptCut_all->Fill(rhad);
300  m_h_photon_wtots1_ptCut_all->Fill(wtots1);
301  m_h_photon_fracs1_ptCut_all->Fill(fracs1);
302  m_h_photon_deltae_ptCut_all->Fill(deltae);
303  m_h_photon_eratio_ptCut_all->Fill(eratio);
304  m_h_photon_f1_ptCut_all->Fill(f1);
305  if (loose_MC15) {
306  m_h_photon_reta_ptCut_loose->Fill(reta);
307  m_h_photon_rphi_ptCut_loose->Fill(rphi);
308  m_h_photon_weta2_ptCut_loose->Fill(weta2);
309  m_h_photon_rhad_ptCut_loose->Fill(rhad);
310  m_h_photon_wtots1_ptCut_loose->Fill(wtots1);
311  m_h_photon_fracs1_ptCut_loose->Fill(fracs1);
312  m_h_photon_deltae_ptCut_loose->Fill(deltae);
313  m_h_photon_eratio_ptCut_loose->Fill(eratio);
314  m_h_photon_f1_ptCut_loose->Fill(f1);
315  }
316  if (tight_MC15) {
317  m_h_photon_reta_ptCut_tight->Fill(reta);
318  m_h_photon_rphi_ptCut_tight->Fill(rphi);
319  m_h_photon_weta2_ptCut_tight->Fill(weta2);
320  m_h_photon_rhad_ptCut_tight->Fill(rhad);
321  m_h_photon_wtots1_ptCut_tight->Fill(wtots1);
322  m_h_photon_fracs1_ptCut_tight->Fill(fracs1);
323  m_h_photon_deltae_ptCut_tight->Fill(deltae);
324  m_h_photon_eratio_ptCut_tight->Fill(eratio);
325  m_h_photon_f1_ptCut_tight->Fill(f1);
326  }
327 
328  m_h_photon_reta_eta_ptCut->Fill(eta, reta);
329  m_h_photon_rphi_eta_ptCut->Fill(eta, rphi);
330  m_h_photon_weta2_eta_ptCut->Fill(eta, weta2);
331 
332  if (fabs(eta) < 1.52) m_h_photon_rhad_eta_ptCut->Fill(eta, rhad);
333  else m_h_photon_rhad_eta_ptCut->Fill(eta, rhad1);
334  m_h_photon_wtots1_eta_ptCut->Fill(eta, wtots1);
335  m_h_photon_fracs1_eta_ptCut->Fill(eta, fracs1);
336 
337  m_h_photon_deltae_eta_ptCut->Fill(eta, deltae);
338  m_h_photon_eratio_eta_ptCut->Fill(eta, eratio);
339  m_h_photon_f1_eta_ptCut->Fill(eta, f1);
340  }
341 
342  return StatusCode::SUCCESS;
343 }
344 
345 StatusCode HIMonitoringPhotonsTool::procHistograms( ) {
346  //if( endOfLowStatFlag() || endOfLumiBlockFlag() ) { }
347 
348 
349  if (endOfRunFlag()) {
350  if (m_h_photon_reta_ptCut_all->GetEntries() > 0) m_h_photon_reta_ptCut_all->Scale(
351  1. / m_h_photon_reta_ptCut_all->GetEntries());
352  if (m_h_photon_reta_ptCut_loose->GetEntries() > 0) m_h_photon_reta_ptCut_loose->Scale(
353  1. / m_h_photon_reta_ptCut_loose->GetEntries());
354  if (m_h_photon_reta_ptCut_tight->GetEntries() > 0) m_h_photon_reta_ptCut_tight->Scale(
355  1. / m_h_photon_reta_ptCut_tight->GetEntries());
356  if (m_h_photon_rphi_ptCut_all->GetEntries() > 0) m_h_photon_rphi_ptCut_all->Scale(1. / m_h_photon_rphi_ptCut_all->GetEntries());
357  if (m_h_photon_rphi_ptCut_loose->GetEntries() > 0) m_h_photon_rphi_ptCut_loose->Scale(1. / m_h_photon_rphi_ptCut_loose->GetEntries());
358  if (m_h_photon_rphi_ptCut_tight->GetEntries() > 0) m_h_photon_rphi_ptCut_tight->Scale(1. / m_h_photon_rphi_ptCut_tight->GetEntries());
359  if (m_h_photon_weta2_ptCut_all->GetEntries() > 0) m_h_photon_weta2_ptCut_all->Scale(1. / m_h_photon_weta2_ptCut_all->GetEntries());
360  if (m_h_photon_weta2_ptCut_loose->GetEntries() > 0) m_h_photon_weta2_ptCut_loose->Scale(1. / m_h_photon_weta2_ptCut_loose->GetEntries());
361  if (m_h_photon_weta2_ptCut_tight->GetEntries() > 0) m_h_photon_weta2_ptCut_tight->Scale(1. / m_h_photon_weta2_ptCut_tight->GetEntries());
362  if (m_h_photon_rhad_ptCut_all->GetEntries() > 0) m_h_photon_rhad_ptCut_all->Scale(1. / m_h_photon_rhad_ptCut_all->GetEntries());
363  if (m_h_photon_rhad_ptCut_loose->GetEntries() > 0) m_h_photon_rhad_ptCut_loose->Scale(1. / m_h_photon_rhad_ptCut_loose->GetEntries());
364  if (m_h_photon_rhad_ptCut_tight->GetEntries() > 0) m_h_photon_rhad_ptCut_tight->Scale(1. / m_h_photon_rhad_ptCut_tight->GetEntries());
365  if (m_h_photon_wtots1_ptCut_all->GetEntries() > 0) m_h_photon_wtots1_ptCut_all->Scale(1. / m_h_photon_wtots1_ptCut_all->GetEntries());
366  if (m_h_photon_wtots1_ptCut_loose->GetEntries() > 0) m_h_photon_wtots1_ptCut_loose->Scale(1. / m_h_photon_wtots1_ptCut_loose->GetEntries());
367  if (m_h_photon_wtots1_ptCut_tight->GetEntries() > 0) m_h_photon_wtots1_ptCut_tight->Scale(1. / m_h_photon_wtots1_ptCut_tight->GetEntries());
368  if (m_h_photon_fracs1_ptCut_all->GetEntries() > 0) m_h_photon_fracs1_ptCut_all->Scale(1. / m_h_photon_fracs1_ptCut_all->GetEntries());
369  if (m_h_photon_fracs1_ptCut_loose->GetEntries() > 0) m_h_photon_fracs1_ptCut_loose->Scale(1. / m_h_photon_fracs1_ptCut_loose->GetEntries());
370  if (m_h_photon_fracs1_ptCut_tight->GetEntries() > 0) m_h_photon_fracs1_ptCut_tight->Scale(1. / m_h_photon_fracs1_ptCut_tight->GetEntries());
371  if (m_h_photon_deltae_ptCut_all->GetEntries() > 0) m_h_photon_deltae_ptCut_all->Scale(1. / m_h_photon_deltae_ptCut_all->GetEntries());
372  if (m_h_photon_deltae_ptCut_loose->GetEntries() > 0) m_h_photon_deltae_ptCut_loose->Scale(1. / m_h_photon_deltae_ptCut_loose->GetEntries());
373  if (m_h_photon_deltae_ptCut_tight->GetEntries() > 0) m_h_photon_deltae_ptCut_tight->Scale(1. / m_h_photon_deltae_ptCut_tight->GetEntries());
374  if (m_h_photon_eratio_ptCut_all->GetEntries() > 0) m_h_photon_eratio_ptCut_all->Scale(1. / m_h_photon_eratio_ptCut_all->GetEntries());
375  if (m_h_photon_eratio_ptCut_loose->GetEntries() > 0) m_h_photon_eratio_ptCut_loose->Scale(1. / m_h_photon_eratio_ptCut_loose->GetEntries());
376  if (m_h_photon_eratio_ptCut_tight->GetEntries() > 0) m_h_photon_eratio_ptCut_tight->Scale(1. / m_h_photon_eratio_ptCut_tight->GetEntries());
377  if (m_h_photon_f1_ptCut_all->GetEntries() > 0) m_h_photon_f1_ptCut_all->Scale(1. / m_h_photon_f1_ptCut_all->GetEntries());
378  if (m_h_photon_f1_ptCut_loose->GetEntries() > 0) m_h_photon_f1_ptCut_loose->Scale(1. / m_h_photon_f1_ptCut_loose->GetEntries());
379  if (m_h_photon_f1_ptCut_tight->GetEntries() > 0) m_h_photon_f1_ptCut_tight->Scale(1. / m_h_photon_f1_ptCut_tight->GetEntries());
380  if (m_h_photon_etcone20_ptCut_fcal0->GetEntries() > 0) m_h_photon_etcone20_ptCut_fcal0->Scale(1. / m_h_photon_etcone20_ptCut_fcal0->GetEntries());
381  if (m_h_photon_etcone20_ptCut_fcal1->GetEntries() > 0) m_h_photon_etcone20_ptCut_fcal1->Scale(1. / m_h_photon_etcone20_ptCut_fcal1->GetEntries());
382  if (m_h_photon_etcone20_ptCut_fcal2->GetEntries() > 0) m_h_photon_etcone20_ptCut_fcal2->Scale(1. / m_h_photon_etcone20_ptCut_fcal2->GetEntries());
383  if (m_h_photon_etcone30_ptCut_fcal0->GetEntries() > 0) m_h_photon_etcone30_ptCut_fcal0->Scale(1. / m_h_photon_etcone30_ptCut_fcal0->GetEntries());
384  if (m_h_photon_etcone30_ptCut_fcal1->GetEntries() > 0) m_h_photon_etcone30_ptCut_fcal1->Scale(1. / m_h_photon_etcone30_ptCut_fcal1->GetEntries());
385  if (m_h_photon_etcone30_ptCut_fcal2->GetEntries() > 0) m_h_photon_etcone30_ptCut_fcal2->Scale(1. / m_h_photon_etcone30_ptCut_fcal2->GetEntries());
386  if (m_h_photon_etcone40_ptCut_fcal0->GetEntries() > 0) m_h_photon_etcone40_ptCut_fcal0->Scale(1. / m_h_photon_etcone40_ptCut_fcal0->GetEntries());
387  if (m_h_photon_etcone40_ptCut_fcal1->GetEntries() > 0) m_h_photon_etcone40_ptCut_fcal1->Scale(1. / m_h_photon_etcone40_ptCut_fcal1->GetEntries());
388  if (m_h_photon_etcone40_ptCut_fcal2->GetEntries() > 0) m_h_photon_etcone40_ptCut_fcal2->Scale(1. / m_h_photon_etcone40_ptCut_fcal2->GetEntries());
389  if (m_h_photon_etcone30_ptCut_fcal0_all->GetEntries() > 0) m_h_photon_etcone30_ptCut_fcal0_all->Scale(1. / m_h_photon_etcone30_ptCut_fcal0_all->GetEntries());
390  if (m_h_photon_etcone30_ptCut_fcal0_loose->GetEntries() > 0) m_h_photon_etcone30_ptCut_fcal0_loose->Scale(1. / m_h_photon_etcone30_ptCut_fcal0_loose->GetEntries());
391  if (m_h_photon_etcone30_ptCut_fcal0_tight->GetEntries() > 0) m_h_photon_etcone30_ptCut_fcal0_tight->Scale(1. / m_h_photon_etcone30_ptCut_fcal0_tight->GetEntries());
392  if (m_h_photon_etcone30_ptCut_fcal1_all->GetEntries() > 0) m_h_photon_etcone30_ptCut_fcal1_all->Scale(1. / m_h_photon_etcone30_ptCut_fcal1_all->GetEntries());
393  if (m_h_photon_etcone30_ptCut_fcal1_loose->GetEntries() > 0) m_h_photon_etcone30_ptCut_fcal1_loose->Scale(1. / m_h_photon_etcone30_ptCut_fcal1_loose->GetEntries());
394  if (m_h_photon_etcone30_ptCut_fcal1_tight->GetEntries() > 0) m_h_photon_etcone30_ptCut_fcal1_tight->Scale(1. / m_h_photon_etcone30_ptCut_fcal1_tight->GetEntries());
395  if (m_h_photon_etcone30_ptCut_fcal2_all->GetEntries() > 0) m_h_photon_etcone30_ptCut_fcal2_all->Scale(1. / m_h_photon_etcone30_ptCut_fcal2_all->GetEntries());
396  if (m_h_photon_etcone30_ptCut_fcal2_loose->GetEntries() > 0) m_h_photon_etcone30_ptCut_fcal2_loose->Scale(1. / m_h_photon_etcone30_ptCut_fcal2_loose->GetEntries());
397  if (m_h_photon_etcone30_ptCut_fcal2_tight->GetEntries() > 0) m_h_photon_etcone30_ptCut_fcal2_tight->Scale(1. / m_h_photon_etcone30_ptCut_fcal2_tight->GetEntries());
398  }
399 
400  return StatusCode::SUCCESS;
401 }
402 
404 void HIMonitoringPhotonsTool::book_hist() {
405  std::string path0 = "HeavyIon/Photons/Basic";
406 
407  m_h_photon_pt_eta = TH2D_LW::create("h_photon_pt_eta", "; photon pT; photon eta", m_PhotonPt_nbins, m_PhotonPt_low, m_PhotonPt_high, m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high);
408  regHist(m_h_photon_pt_eta, path0, run).ignore();
409 
410  m_h_photon_pt_phi = TH2D_LW::create("h_photon_pt_phi", "; photon pT; photon phi", m_PhotonPt_nbins, m_PhotonPt_low, m_PhotonPt_high, m_PhotonPhi_nbins, m_PhotonPhi_low, m_PhotonPhi_high);
411  regHist(m_h_photon_pt_phi, path0, run).ignore();
412 
413  m_h_photon_pt_fcal = TH2D_LW::create("h_photon_pt_fcal", "; photon pT; FCal SumET", m_PhotonPt_nbins, m_PhotonPt_low, m_PhotonPt_high, m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high);
414  regHist(m_h_photon_pt_fcal, path0, run).ignore();
415 
416  m_h_photon_eta_phi_ptCut = TH2D_LW::create("h_photon_eta_phi_ptCut", "; photon eta; photon phi", m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high, m_PhotonPhi_nbins, m_PhotonPhi_low, m_PhotonPhi_high);
417  regHist(m_h_photon_eta_phi_ptCut, path0, run).ignore();
418 
419  // tight versions
420  m_h_photon_pt_eta_tight = TH2D_LW::create("h_photon_pt_eta_tight", "; photon pT; photon eta", m_PhotonPt_nbins, m_PhotonPt_low, m_PhotonPt_high, m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high);
421  regHist(m_h_photon_pt_eta_tight, path0, run).ignore();
422 
423  m_h_photon_pt_phi_tight = TH2D_LW::create("h_photon_pt_phi_tight", "; photon pT; photon phi", m_PhotonPt_nbins, m_PhotonPt_low, m_PhotonPt_high, m_PhotonPhi_nbins, m_PhotonPhi_low, m_PhotonPhi_high);
424  regHist(m_h_photon_pt_phi_tight, path0, run).ignore();
425 
426  m_h_photon_pt_fcal_tight = TH2D_LW::create("h_photon_pt_fcal_tight", "; photon pT; FCal SumET", m_PhotonPt_nbins, m_PhotonPt_low, m_PhotonPt_high, m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high);
427  regHist(m_h_photon_pt_fcal_tight, path0, run).ignore();
428 
429  m_h_photon_eta_phi_ptCut_tight = TH2D_LW::create("h_photon_eta_phi_ptCut_tight", "; photon eta; photon phi", m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high, m_PhotonPhi_nbins, m_PhotonPhi_low, m_PhotonPhi_high);
430  regHist(m_h_photon_eta_phi_ptCut_tight, path0, run).ignore();
431 
432  std::string path1 = "HeavyIon/Photons/ShowerShapes";
433 
434  // 1-D vs. ID cuts
435 
436  m_h_photon_reta_ptCut_all = new TH1D("h_photon_reta_ptCut_all", "; photon reta ; ", m_reta_nbins, m_reta_low, m_reta_high);
437  regHist(m_h_photon_reta_ptCut_all, path1, run).ignore();
438  m_h_photon_reta_ptCut_loose = new TH1D("h_photon_reta_ptCut_loose", "; photon reta ; ", m_reta_nbins, m_reta_low, m_reta_high);
439  regHist(m_h_photon_reta_ptCut_loose, path1, run).ignore();
440  m_h_photon_reta_ptCut_tight = new TH1D("h_photon_reta_ptCut_tight", "; photon reta ; ", m_reta_nbins, m_reta_low, m_reta_high);
441  regHist(m_h_photon_reta_ptCut_tight, path1, run).ignore();
442 
443  m_h_photon_rphi_ptCut_all = new TH1D("h_photon_rphi_ptCut_all", "; photon rphi ; ", m_rphi_nbins, m_rphi_low, m_rphi_high);
444  regHist(m_h_photon_rphi_ptCut_all, path1, run).ignore();
445  m_h_photon_rphi_ptCut_loose = new TH1D("h_photon_rphi_ptCut_loose", "; photon rphi ; ", m_rphi_nbins, m_rphi_low, m_rphi_high);
446  regHist(m_h_photon_rphi_ptCut_loose, path1, run).ignore();
447  m_h_photon_rphi_ptCut_tight = new TH1D("h_photon_rphi_ptCut_tight", "; photon rphi ; ", m_rphi_nbins, m_rphi_low, m_rphi_high);
448  regHist(m_h_photon_rphi_ptCut_tight, path1, run).ignore();
449 
450  m_h_photon_weta2_ptCut_all = new TH1D("h_photon_weta2_ptCut_all", "; photon weta2 ; ", m_weta2_nbins, m_weta2_low, m_weta2_high);
451  regHist(m_h_photon_weta2_ptCut_all, path1, run).ignore();
452  m_h_photon_weta2_ptCut_loose = new TH1D("h_photon_weta2_ptCut_loose", "; photon weta2 ; ", m_weta2_nbins, m_weta2_low, m_weta2_high);
453  regHist(m_h_photon_weta2_ptCut_loose, path1, run).ignore();
454  m_h_photon_weta2_ptCut_tight = new TH1D("h_photon_weta2_ptCut_tight", "; photon weta2 ; ", m_weta2_nbins, m_weta2_low, m_weta2_high);
455  regHist(m_h_photon_weta2_ptCut_tight, path1, run).ignore();
456 
457  m_h_photon_rhad_ptCut_all = new TH1D("h_photon_rhad_ptCut_all", "; photon rhad ; ", m_rhad_nbins, m_rhad_low, m_rhad_high);
458  regHist(m_h_photon_rhad_ptCut_all, path1, run).ignore();
459  m_h_photon_rhad_ptCut_loose = new TH1D("h_photon_rhad_ptCut_loose", "; photon rhad ; ", m_rhad_nbins, m_rhad_low, m_rhad_high);
460  regHist(m_h_photon_rhad_ptCut_loose, path1, run).ignore();
461  m_h_photon_rhad_ptCut_tight = new TH1D("h_photon_rhad_ptCut_tight", "; photon rhad ; ", m_rhad_nbins, m_rhad_low, m_rhad_high);
462  regHist(m_h_photon_rhad_ptCut_tight, path1, run).ignore();
463 
464  m_h_photon_wtots1_ptCut_all = new TH1D("h_photon_wtots1_ptCut_all", "; photon wtots1 ; ", m_wtots1_nbins, m_wtots1_low, m_wtots1_high);
465  regHist(m_h_photon_wtots1_ptCut_all, path1, run).ignore();
466  m_h_photon_wtots1_ptCut_loose = new TH1D("h_photon_wtots1_ptCut_loose", "; photon wtots1 ; ", m_wtots1_nbins, m_wtots1_low, m_wtots1_high);
467  regHist(m_h_photon_wtots1_ptCut_loose, path1, run).ignore();
468  m_h_photon_wtots1_ptCut_tight = new TH1D("h_photon_wtots1_ptCut_tight", "; photon wtots1 ; ", m_wtots1_nbins, m_wtots1_low, m_wtots1_high);
469  regHist(m_h_photon_wtots1_ptCut_tight, path1, run).ignore();
470 
471  m_h_photon_fracs1_ptCut_all = new TH1D("h_photon_fracs1_ptCut_all", "; photon fracs1 ; ", m_fracs1_nbins, m_fracs1_low, m_fracs1_high);
472  regHist(m_h_photon_fracs1_ptCut_all, path1, run).ignore();
473  m_h_photon_fracs1_ptCut_loose = new TH1D("h_photon_fracs1_ptCut_loose", "; photon fracs1 ; ", m_fracs1_nbins, m_fracs1_low, m_fracs1_high);
474  regHist(m_h_photon_fracs1_ptCut_loose, path1, run).ignore();
475  m_h_photon_fracs1_ptCut_tight = new TH1D("h_photon_fracs1_ptCut_tight", "; photon fracs1 ; ", m_fracs1_nbins, m_fracs1_low, m_fracs1_high);
476  regHist(m_h_photon_fracs1_ptCut_tight, path1, run).ignore();
477 
478  m_h_photon_deltae_ptCut_all = new TH1D("h_photon_deltae_ptCut_all", "; photon deltae ; ", m_deltae_nbins, m_deltae_low, m_deltae_high);
479  regHist(m_h_photon_deltae_ptCut_all, path1, run).ignore();
480  m_h_photon_deltae_ptCut_loose = new TH1D("h_photon_deltae_ptCut_loose", "; photon deltae ; ", m_deltae_nbins, m_deltae_low, m_deltae_high);
481  regHist(m_h_photon_deltae_ptCut_loose, path1, run).ignore();
482  m_h_photon_deltae_ptCut_tight = new TH1D("h_photon_deltae_ptCut_tight", "; photon deltae ; ", m_deltae_nbins, m_deltae_low, m_deltae_high);
483  regHist(m_h_photon_deltae_ptCut_tight, path1, run).ignore();
484 
485  m_h_photon_eratio_ptCut_all = new TH1D("h_photon_eratio_ptCut_all", "; photon eratio ; ", m_eratio_nbins, m_eratio_low, m_eratio_high);
486  regHist(m_h_photon_eratio_ptCut_all, path1, run).ignore();
487  m_h_photon_eratio_ptCut_loose = new TH1D("h_photon_eratio_ptCut_loose", "; photon eratio ; ", m_eratio_nbins, m_eratio_low, m_eratio_high);
488  regHist(m_h_photon_eratio_ptCut_loose, path1, run).ignore();
489  m_h_photon_eratio_ptCut_tight = new TH1D("h_photon_eratio_ptCut_tight", "; photon eratio ; ", m_eratio_nbins, m_eratio_low, m_eratio_high);
490  regHist(m_h_photon_eratio_ptCut_tight, path1, run).ignore();
491 
492  m_h_photon_f1_ptCut_all = new TH1D("h_photon_f1_ptCut_all", "; photon f1 ; ", m_f1_nbins, m_f1_low, m_f1_high);
493  regHist(m_h_photon_f1_ptCut_all, path1, run).ignore();
494  m_h_photon_f1_ptCut_loose = new TH1D("h_photon_f1_ptCut_loose", "; photon f1 ; ", m_f1_nbins, m_f1_low, m_f1_high);
495  regHist(m_h_photon_f1_ptCut_loose, path1, run).ignore();
496  m_h_photon_f1_ptCut_tight = new TH1D("h_photon_f1_ptCut_tight", "; photon f1 ; ", m_f1_nbins, m_f1_low, m_f1_high);
497  regHist(m_h_photon_f1_ptCut_tight, path1, run).ignore();
498 
499  // vs. eta
500  m_h_photon_reta_eta_ptCut = TH2D_LW::create("h_photon_reta_eta_ptCut", "; photon eta; photon reta", m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high, m_reta_nbins, m_reta_low, m_reta_high);
501  regHist(m_h_photon_reta_eta_ptCut, path1, run).ignore();
502  m_h_photon_rphi_eta_ptCut = TH2D_LW::create("h_photon_rphi_eta_ptCut", "; photon eta; photon rphi", m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high, m_rphi_nbins, m_rphi_low, m_rphi_high);
503  regHist(m_h_photon_rphi_eta_ptCut, path1, run).ignore();
504  m_h_photon_weta2_eta_ptCut = TH2D_LW::create("h_photon_weta2_eta_ptCut", "; photon eta; photon weta2", m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high, m_weta2_nbins, m_weta2_low, m_weta2_high);
505  regHist(m_h_photon_weta2_eta_ptCut, path1, run).ignore();
506 
507  m_h_photon_rhad_eta_ptCut = TH2D_LW::create("h_photon_rhad_eta_ptCut", "; photon eta; photon rhad", m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high, m_rhad_nbins, m_rhad_low, m_rhad_high);
508  regHist(m_h_photon_rhad_eta_ptCut, path1, run).ignore();
509  m_h_photon_wtots1_eta_ptCut = TH2D_LW::create("h_photon_wtots1_eta_ptCut", "; photon eta; photon wtots1", m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high, m_wtots1_nbins, m_wtots1_low, m_wtots1_high);
510  regHist(m_h_photon_wtots1_eta_ptCut, path1, run).ignore();
511  m_h_photon_fracs1_eta_ptCut = TH2D_LW::create("h_photon_fracs1_eta_ptCut", "; photon eta; photon fracs1", m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high, m_fracs1_nbins, m_fracs1_low, m_fracs1_high);
512  regHist(m_h_photon_fracs1_eta_ptCut, path1, run).ignore();
513 
514  m_h_photon_deltae_eta_ptCut = TH2D_LW::create("h_photon_deltae_eta_ptCut", "; photon eta; photon deltae", m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high, m_deltae_nbins, m_deltae_low, m_deltae_high);
515  regHist(m_h_photon_deltae_eta_ptCut, path1, run).ignore();
516  m_h_photon_eratio_eta_ptCut = TH2D_LW::create("h_photon_eratio_eta_ptCut", "; photon eta; photon eratio", m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high, m_eratio_nbins, m_eratio_low, m_eratio_high);
517  regHist(m_h_photon_eratio_eta_ptCut, path1, run).ignore();
518  m_h_photon_f1_eta_ptCut = TH2D_LW::create("h_photon_f1_eta_ptCut", "; photon eta; photon f1", m_PhotonEta_nbins, m_PhotonEta_low, m_PhotonEta_high, m_f1_nbins, m_f1_low, m_f1_high);
519  regHist(m_h_photon_f1_eta_ptCut, path1, run).ignore();
520 
521 
522  // vs. FCal ET
523 
524  m_h_photon_reta_fcal_ptCut_eta0 = TH2D_LW::create("h_photon_reta_fcal_ptCut_eta0", "; FCal Sum ET; photon reta", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_reta_nbins, m_reta_low, m_reta_high);
525  regHist(m_h_photon_reta_fcal_ptCut_eta0, path1, run).ignore();
526  m_h_photon_reta_fcal_ptCut_eta1 = TH2D_LW::create("h_photon_reta_fcal_ptCut_eta1", "; FCal Sum ET; photon reta", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_reta_nbins, m_reta_low, m_reta_high);
527  regHist(m_h_photon_reta_fcal_ptCut_eta1, path1, run).ignore();
528 
529  m_h_photon_rphi_fcal_ptCut_eta0 = TH2D_LW::create("h_photon_rphi_fcal_ptCut_eta0", "; FCal Sum ET; photon rphi", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_rphi_nbins, m_rphi_low, m_rphi_high);
530  regHist(m_h_photon_rphi_fcal_ptCut_eta0, path1, run).ignore();
531  m_h_photon_rphi_fcal_ptCut_eta1 = TH2D_LW::create("h_photon_rphi_fcal_ptCut_eta1", "; FCal Sum ET; photon rphi", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_rphi_nbins, m_rphi_low, m_rphi_high);
532  regHist(m_h_photon_rphi_fcal_ptCut_eta1, path1, run).ignore();
533 
534  m_h_photon_weta2_fcal_ptCut_eta0 = TH2D_LW::create("h_photon_weta2_fcal_ptCut_eta0", "; FCal Sum ET; photon weta2", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_weta2_nbins, m_weta2_low, m_weta2_high);
535  regHist(m_h_photon_weta2_fcal_ptCut_eta0, path1, run).ignore();
536  m_h_photon_weta2_fcal_ptCut_eta1 = TH2D_LW::create("h_photon_weta2_fcal_ptCut_eta1", "; FCal Sum ET; photon weta2", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_weta2_nbins, m_weta2_low, m_weta2_high);
537  regHist(m_h_photon_weta2_fcal_ptCut_eta1, path1, run).ignore();
538 
539 
540  m_h_photon_rhad_fcal_ptCut_eta0 = TH2D_LW::create("h_photon_rhad_fcal_ptCut_eta0", "; FCal Sum ET; photon rhad", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_rhad_nbins, m_rhad_low, m_rhad_high);
541  regHist(m_h_photon_rhad_fcal_ptCut_eta0, path1, run).ignore();
542  m_h_photon_rhad_fcal_ptCut_eta1 = TH2D_LW::create("h_photon_rhad_fcal_ptCut_eta1", "; FCal Sum ET; photon rhad", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_rhad_nbins, m_rhad_low, m_rhad_high);
543  regHist(m_h_photon_rhad_fcal_ptCut_eta1, path1, run).ignore();
544 
545  m_h_photon_wtots1_fcal_ptCut_eta0 = TH2D_LW::create("h_photon_wtots1_fcal_ptCut_eta0", "; FCal Sum ET; photon wtots1", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_wtots1_nbins, m_wtots1_low, m_wtots1_high);
546  regHist(m_h_photon_wtots1_fcal_ptCut_eta0, path1, run).ignore();
547  m_h_photon_wtots1_fcal_ptCut_eta1 = TH2D_LW::create("h_photon_wtots1_fcal_ptCut_eta1", "; FCal Sum ET; photon wtots1", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_wtots1_nbins, m_wtots1_low, m_wtots1_high);
548  regHist(m_h_photon_wtots1_fcal_ptCut_eta1, path1, run).ignore();
549 
550  m_h_photon_fracs1_fcal_ptCut_eta0 = TH2D_LW::create("h_photon_fracs1_fcal_ptCut_eta0", "; FCal Sum ET; photon fracs1", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_fracs1_nbins, m_fracs1_low, m_fracs1_high);
551  regHist(m_h_photon_fracs1_fcal_ptCut_eta0, path1, run).ignore();
552  m_h_photon_fracs1_fcal_ptCut_eta1 = TH2D_LW::create("h_photon_fracs1_fcal_ptCut_eta1", "; FCal Sum ET; photon fracs1", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_fracs1_nbins, m_fracs1_low, m_fracs1_high);
553  regHist(m_h_photon_fracs1_fcal_ptCut_eta1, path1, run).ignore();
554 
555 
556  m_h_photon_deltae_fcal_ptCut_eta0 = TH2D_LW::create("h_photon_deltae_fcal_ptCut_eta0", "; FCal Sum ET; photon deltae", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_deltae_nbins, m_deltae_low, m_deltae_high);
557  regHist(m_h_photon_deltae_fcal_ptCut_eta0, path1, run).ignore();
558  m_h_photon_deltae_fcal_ptCut_eta1 = TH2D_LW::create("h_photon_deltae_fcal_ptCut_eta1", "; FCal Sum ET; photon deltae", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_deltae_nbins, m_deltae_low, m_deltae_high);
559  regHist(m_h_photon_deltae_fcal_ptCut_eta1, path1, run).ignore();
560 
561  m_h_photon_eratio_fcal_ptCut_eta0 = TH2D_LW::create("h_photon_eratio_fcal_ptCut_eta0", "; FCal Sum ET; photon eratio", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_eratio_nbins, m_eratio_low, m_eratio_high);
562  regHist(m_h_photon_eratio_fcal_ptCut_eta0, path1, run).ignore();
563  m_h_photon_eratio_fcal_ptCut_eta1 = TH2D_LW::create("h_photon_eratio_fcal_ptCut_eta1", "; FCal Sum ET; photon eratio", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_eratio_nbins, m_eratio_low, m_eratio_high);
564  regHist(m_h_photon_eratio_fcal_ptCut_eta1, path1, run).ignore();
565 
566  m_h_photon_f1_fcal_ptCut_eta0 = TH2D_LW::create("h_photon_f1_fcal_ptCut_eta0", "; FCal Sum ET; photon f1", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_f1_nbins, m_f1_low, m_f1_high);
567  regHist(m_h_photon_f1_fcal_ptCut_eta0, path1, run).ignore();
568  m_h_photon_f1_fcal_ptCut_eta1 = TH2D_LW::create("h_photon_f1_fcal_ptCut_eta1", "; FCal Sum ET; photon f1", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_f1_nbins, m_f1_low, m_f1_high);
569  regHist(m_h_photon_f1_fcal_ptCut_eta1, path1, run).ignore();
570 
571  std::string path2 = "HeavyIon/Photons/Isolation";
572 
573  m_h_photon_fcal_etcone20_ptCut = TH2D_LW::create("h_photon_fcal_etcone20_ptCut", "; FCal Sum ET; photon etcone20", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_IsoWide_nbins, m_IsoWide_low, m_IsoWide_high);
574  regHist(m_h_photon_fcal_etcone20_ptCut, path2, run).ignore();
575 
576  m_h_photon_fcal_etcone30_ptCut = TH2D_LW::create("h_photon_fcal_etcone30_ptCut", "; FCal Sum ET; photon etcone30", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_IsoWide_nbins, m_IsoWide_low, m_IsoWide_high);
577  regHist(m_h_photon_fcal_etcone30_ptCut, path2, run).ignore();
578 
579  m_h_photon_fcal_etcone40_ptCut = TH2D_LW::create("h_photon_fcal_etcone40_ptCut", "; FCal Sum ET; photon etcone40", m_FCalEt_nbins, m_FCalEt_low, m_FCalEt_high, m_IsoWide_nbins, m_IsoWide_low, m_IsoWide_high);
580  regHist(m_h_photon_fcal_etcone40_ptCut, path2, run).ignore();
581 
582  m_h_photon_pt_etcone40_fcal0 = TH2D_LW::create("h_photon_pt_etcone40_fcal0", "; photon pT; photon etcone40", m_PhotonPt_nbins, m_PhotonPt_low, m_PhotonPt_high, m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
583  regHist(m_h_photon_pt_etcone40_fcal0, path2, run).ignore();
584 
585  m_h_photon_pt_etcone40_fcal1 = TH2D_LW::create("h_photon_pt_etcone40_fcal1", "; photon pT; photon etcone40", m_PhotonPt_nbins, m_PhotonPt_low, m_PhotonPt_high, m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
586  regHist(m_h_photon_pt_etcone40_fcal1, path2, run).ignore();
587 
588  m_h_photon_pt_etcone40_fcal2 = TH2D_LW::create("h_photon_pt_etcone40_fcal2", "; photon pT; photon etcone40", m_PhotonPt_nbins, m_PhotonPt_low, m_PhotonPt_high, m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
589  regHist(m_h_photon_pt_etcone40_fcal2, path2, run).ignore();
590 
591  // 1-D vs. FCal ET
592  m_h_photon_etcone20_ptCut_fcal0 = new TH1D("h_photon_etcone20_ptCut_fcal0", "; photon etcone20", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
593  regHist(m_h_photon_etcone20_ptCut_fcal0, path2, run).ignore();
594 
595  m_h_photon_etcone20_ptCut_fcal1 = new TH1D("h_photon_etcone20_ptCut_fcal1", "; photon etcone20", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
596  regHist(m_h_photon_etcone20_ptCut_fcal1, path2, run).ignore();
597 
598  m_h_photon_etcone20_ptCut_fcal2 = new TH1D("h_photon_etcone20_ptCut_fcal2", "; photon etcone20", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
599  regHist(m_h_photon_etcone20_ptCut_fcal2, path2, run).ignore();
600 
601  m_h_photon_etcone30_ptCut_fcal0 = new TH1D("h_photon_etcone30_ptCut_fcal0", "; photon etcone30", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
602  regHist(m_h_photon_etcone30_ptCut_fcal0, path2, run).ignore();
603 
604  m_h_photon_etcone30_ptCut_fcal1 = new TH1D("h_photon_etcone30_ptCut_fcal1", "; photon etcone30", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
605  regHist(m_h_photon_etcone30_ptCut_fcal1, path2, run).ignore();
606 
607  m_h_photon_etcone30_ptCut_fcal2 = new TH1D("h_photon_etcone30_ptCut_fcal2", "; photon etcone30", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
608  regHist(m_h_photon_etcone30_ptCut_fcal2, path2, run).ignore();
609 
610  m_h_photon_etcone40_ptCut_fcal0 = new TH1D("h_photon_etcone40_ptCut_fcal0", "; photon etcone40", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
611  regHist(m_h_photon_etcone40_ptCut_fcal0, path2, run).ignore();
612 
613  m_h_photon_etcone40_ptCut_fcal1 = new TH1D("h_photon_etcone40_ptCut_fcal1", "; photon etcone40", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
614  regHist(m_h_photon_etcone40_ptCut_fcal1, path2, run).ignore();
615 
616  m_h_photon_etcone40_ptCut_fcal2 = new TH1D("h_photon_etcone40_ptCut_fcal2", "; photon etcone40", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
617  regHist(m_h_photon_etcone40_ptCut_fcal2, path2, run).ignore();
618 
619  // 1-D vs. ID
620 
621  m_h_photon_etcone30_ptCut_fcal0_all = new TH1D("h_photon_etcone30_ptCut_fcal0_all", "; photon etcone30", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
622  regHist(m_h_photon_etcone30_ptCut_fcal0_all, path2, run).ignore();
623  m_h_photon_etcone30_ptCut_fcal0_loose = new TH1D("h_photon_etcone30_ptCut_fcal0_loose", "; photon etcone30", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
624  regHist(m_h_photon_etcone30_ptCut_fcal0_loose, path2, run).ignore();
625  m_h_photon_etcone30_ptCut_fcal0_tight = new TH1D("h_photon_etcone30_ptCut_fcal0_tight", "; photon etcone30", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
626  regHist(m_h_photon_etcone30_ptCut_fcal0_tight, path2, run).ignore();
627 
628  m_h_photon_etcone30_ptCut_fcal1_all = new TH1D("h_photon_etcone30_ptCut_fcal1_all", "; photon etcone30", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
629  regHist(m_h_photon_etcone30_ptCut_fcal1_all, path2, run).ignore();
630  m_h_photon_etcone30_ptCut_fcal1_loose = new TH1D("h_photon_etcone30_ptCut_fcal1_loose", "; photon etcone30", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
631  regHist(m_h_photon_etcone30_ptCut_fcal1_loose, path2, run).ignore();
632  m_h_photon_etcone30_ptCut_fcal1_tight = new TH1D("h_photon_etcone30_ptCut_fcal1_tight", "; photon etcone30", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
633  regHist(m_h_photon_etcone30_ptCut_fcal1_tight, path2, run).ignore();
634 
635  m_h_photon_etcone30_ptCut_fcal2_all = new TH1D("h_photon_etcone30_ptCut_fcal2_all", "; photon etcone30", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
636  regHist(m_h_photon_etcone30_ptCut_fcal2_all, path2, run).ignore();
637  m_h_photon_etcone30_ptCut_fcal2_loose = new TH1D("h_photon_etcone30_ptCut_fcal2_loose", "; photon etcone30", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
638  regHist(m_h_photon_etcone30_ptCut_fcal2_loose, path2, run).ignore();
639  m_h_photon_etcone30_ptCut_fcal2_tight = new TH1D("h_photon_etcone30_ptCut_fcal2_tight", "; photon etcone30", m_IsoNarrow_nbins, m_IsoNarrow_low, m_IsoNarrow_high);
640  regHist(m_h_photon_etcone30_ptCut_fcal2_tight, path2, run).ignore();
641 }
HIMonitoringPhotonsTool::m_h_photon_rphi_fcal_ptCut_eta0
TH2D_LW * m_h_photon_rphi_fcal_ptCut_eta0
Definition: HIMonitoringPhotonsTool.h:168
HIMonitoringPhotonsTool::m_h_photon_etcone30_ptCut_fcal2_all
TH1D * m_h_photon_etcone30_ptCut_fcal2_all
Definition: HIMonitoringPhotonsTool.h:215
HIMonitoringPhotonsTool::m_PhotonEta_nbins
int m_PhotonEta_nbins
Definition: HIMonitoringPhotonsTool.h:46
HIMonitoringPhotonsTool::m_h_photon_deltae_fcal_ptCut_eta1
TH2D_LW * m_h_photon_deltae_fcal_ptCut_eta1
Definition: HIMonitoringPhotonsTool.h:179
AsgPhotonIsEMSelector::initialize
virtual StatusCode initialize()
AlgTool initialize method.
Definition: AsgPhotonIsEMSelector.cxx:89
xAOD::eratio
setCharge setNTRTHiThresholdHits eratio
Definition: TrigElectron_v1.cxx:96
HIMonitoringPhotonsTool::m_h_photon_f1_fcal_ptCut_eta1
TH2D_LW * m_h_photon_f1_fcal_ptCut_eta1
Definition: HIMonitoringPhotonsTool.h:183
HIMonitoringPhotonsTool::m_wtots1_nbins
int m_wtots1_nbins
Definition: HIMonitoringPhotonsTool.h:70
HIMonitoringPhotonsTool::m_FCalEt
double m_FCalEt
Definition: HIMonitoringPhotonsTool.h:33
egammaPID::PhotonTight
const unsigned int PhotonTight
Tight photon selection.
Definition: egammaPIDdefs.h:602
DataModel_detail::const_iterator
Const iterator class for DataVector/DataList.
Definition: DVLIterator.h:82
egammaPIDdefs.h
HIMonitoringPhotonsTool::m_h_photon_f1_ptCut_loose
TH1D * m_h_photon_f1_ptCut_loose
Definition: HIMonitoringPhotonsTool.h:137
HIMonitoringPhotonsTool::m_h_photon_fracs1_ptCut_tight
TH1D * m_h_photon_fracs1_ptCut_tight
Definition: HIMonitoringPhotonsTool.h:146
HIMonitoringPhotonsTool::m_deltae_nbins
int m_deltae_nbins
Definition: HIMonitoringPhotonsTool.h:78
phi
Scalar phi() const
phi method
Definition: AmgMatrixBasePlugin.h:64
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
HIMonitoringPhotonsTool::m_h_photon_pt_eta
TH2D_LW * m_h_photon_pt_eta
histograms
Definition: HIMonitoringPhotonsTool.h:101
HIMonitoringPhotonsTool::m_h_photon_weta2_fcal_ptCut_eta1
TH2D_LW * m_h_photon_weta2_fcal_ptCut_eta1
Definition: HIMonitoringPhotonsTool.h:171
HIMonitoringPhotonsTool::m_rphi_nbins
int m_rphi_nbins
Definition: HIMonitoringPhotonsTool.h:58
HIMonitoringPhotonsTool::m_photonLooseIsEMSelector
AsgPhotonIsEMSelector * m_photonLooseIsEMSelector
Definition: HIMonitoringPhotonsTool.h:31
egammaPID::PhotonLoose
const unsigned int PhotonLoose
Loose photon selection.
Definition: egammaPIDdefs.h:586
ManagedMonitorToolBase
Provides functionality for users to implement and save histograms, ntuples, and summary data,...
Definition: ManagedMonitorToolBase.h:74
HIMonitoringPhotonsTool::m_rphi_high
float m_rphi_high
Definition: HIMonitoringPhotonsTool.h:60
HIMonitoringPhotonsTool::m_rhad_nbins
int m_rhad_nbins
Definition: HIMonitoringPhotonsTool.h:66
HIMonitoringPhotonsTool::m_h_photon_wtots1_eta_ptCut
TH2D_LW * m_h_photon_wtots1_eta_ptCut
Definition: HIMonitoringPhotonsTool.h:159
eta
Scalar eta() const
pseudorapidity method
Definition: AmgMatrixBasePlugin.h:79
HIMonitoringPhotonsTool::m_h_photon_etcone40_ptCut_fcal1
TH1D * m_h_photon_etcone40_ptCut_fcal1
Definition: HIMonitoringPhotonsTool.h:204
HIMonitoringPhotonsTool::m_eratio_high
float m_eratio_high
Definition: HIMonitoringPhotonsTool.h:84
TH1D_LW.h
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
TH2D_LW.h
HIMonitoringPhotonsTool::m_h_photon_etcone30_ptCut_fcal0_tight
TH1D * m_h_photon_etcone30_ptCut_fcal0_tight
Definition: HIMonitoringPhotonsTool.h:209
HIMonitoringPhotonsTool::m_FCalEt_low
float m_FCalEt_low
Definition: HIMonitoringPhotonsTool.h:36
HIMonitoringPhotonsTool::m_h_photon_deltae_fcal_ptCut_eta0
TH2D_LW * m_h_photon_deltae_fcal_ptCut_eta0
Definition: HIMonitoringPhotonsTool.h:178
AthenaMonManager::cosmics
@ cosmics
Definition: AthenaMonManager.h:58
HIMonitoringPhotonsTool::m_IsoWide_low
float m_IsoWide_low
Definition: HIMonitoringPhotonsTool.h:91
HIMonitoringPhotonsTool::m_h_photon_pt_fcal
TH2D_LW * m_h_photon_pt_fcal
Definition: HIMonitoringPhotonsTool.h:103
HIMonitoringPhotonsTool::m_h_photon_etcone20_ptCut_fcal0
TH1D * m_h_photon_etcone20_ptCut_fcal0
Definition: HIMonitoringPhotonsTool.h:195
HIMonitoringPhotonsTool::m_PhotonPt_nbins
int m_PhotonPt_nbins
Definition: HIMonitoringPhotonsTool.h:42
HIMonitoringPhotonsTool::m_h_photon_pt_etcone40_fcal1
TH2D_LW * m_h_photon_pt_etcone40_fcal1
Definition: HIMonitoringPhotonsTool.h:192
TH1D
Definition: rootspy.cxx:342
test_pyathena.pt
pt
Definition: test_pyathena.py:11
HIMonitoringPhotonsTool::m_PhotonPhi_low
float m_PhotonPhi_low
Definition: HIMonitoringPhotonsTool.h:51
xAOD::Iso::etcone40
@ etcone40
Definition: IsolationType.h:34
HIMonitoringPhotonsTool::m_fracs1_high
float m_fracs1_high
Definition: HIMonitoringPhotonsTool.h:76
HIMonitoringPhotonsTool::m_h_photon_eratio_ptCut_loose
TH1D * m_h_photon_eratio_ptCut_loose
Definition: HIMonitoringPhotonsTool.h:136
ManagedMonitorToolBase::m_dataType
AthenaMonManager::DataType_t m_dataType
Definition: ManagedMonitorToolBase.h:901
HIMonitoringPhotonsTool::m_h_photon_f1_fcal_ptCut_eta0
TH2D_LW * m_h_photon_f1_fcal_ptCut_eta0
Definition: HIMonitoringPhotonsTool.h:182
xAOD::Iso::etcone30
@ etcone30
Definition: IsolationType.h:33
HIMonitoringPhotonsTool::m_h_photon_reta_ptCut_tight
TH1D * m_h_photon_reta_ptCut_tight
Definition: HIMonitoringPhotonsTool.h:140
HIMonitoringPhotonsTool::m_h_photon_f1_ptCut_tight
TH1D * m_h_photon_f1_ptCut_tight
Definition: HIMonitoringPhotonsTool.h:150
HIMonitoringPhotonsTool::m_h_photon_eratio_ptCut_all
TH1D * m_h_photon_eratio_ptCut_all
Definition: HIMonitoringPhotonsTool.h:123
SG::ConstAccessor< float >
HIMonitoringPhotonsTool::m_IsoNarrow_low
float m_IsoNarrow_low
Definition: HIMonitoringPhotonsTool.h:95
HIMonitoringPhotonsTool::m_reta_nbins
int m_reta_nbins
Definition: HIMonitoringPhotonsTool.h:54
HIMonitoringPhotonsTool.h
HIMonitoringPhotonsTool::m_weta2_high
float m_weta2_high
Definition: HIMonitoringPhotonsTool.h:64
HIMonitoringPhotonsTool::m_h_photon_reta_ptCut_loose
TH1D * m_h_photon_reta_ptCut_loose
Definition: HIMonitoringPhotonsTool.h:127
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
HIMonitoringPhotonsTool::m_h_photon_weta2_eta_ptCut
TH2D_LW * m_h_photon_weta2_eta_ptCut
Definition: HIMonitoringPhotonsTool.h:156
HIMonitoringPhotonsTool::m_rphi_low
float m_rphi_low
Definition: HIMonitoringPhotonsTool.h:59
HIMonitoringPhotonsTool::m_reta_low
float m_reta_low
Definition: HIMonitoringPhotonsTool.h:55
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
HIMonitoringPhotonsTool::m_h_photon_reta_ptCut_all
TH1D * m_h_photon_reta_ptCut_all
Definition: HIMonitoringPhotonsTool.h:114
HIMonitoringPhotonsTool::m_deltae_high
float m_deltae_high
Definition: HIMonitoringPhotonsTool.h:80
HIMonitoringPhotonsTool::m_PhotonPhi_high
float m_PhotonPhi_high
Definition: HIMonitoringPhotonsTool.h:52
HIMonitoringPhotonsTool::m_f1_nbins
int m_f1_nbins
Definition: HIMonitoringPhotonsTool.h:86
HIMonitoringPhotonsTool::m_h_photon_wtots1_ptCut_loose
TH1D * m_h_photon_wtots1_ptCut_loose
Definition: HIMonitoringPhotonsTool.h:132
HIMonitoringPhotonsTool::m_h_photon_fracs1_ptCut_all
TH1D * m_h_photon_fracs1_ptCut_all
Definition: HIMonitoringPhotonsTool.h:120
HIMonitoringPhotonsTool::m_h_photon_rhad_ptCut_all
TH1D * m_h_photon_rhad_ptCut_all
Definition: HIMonitoringPhotonsTool.h:118
python.setupRTTAlg.size
int size
Definition: setupRTTAlg.py:39
HIMonitoringPhotonsTool::m_FCalEt_nbins
int m_FCalEt_nbins
Definition: HIMonitoringPhotonsTool.h:35
xAOD::HIEventShape_v2
Interface class for the HI reconstruction EDM.
Definition: HIEventShape_v2.h:31
AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:85
HIMonitoringPhotonsTool::bookHistograms
virtual StatusCode bookHistograms()
An inheriting class should either override this function or bookHists().
Definition: HIMonitoringPhotonsTool.cxx:123
HIMonitoringPhotonsTool::m_h_photon_etcone30_ptCut_fcal0_loose
TH1D * m_h_photon_etcone30_ptCut_fcal0_loose
Definition: HIMonitoringPhotonsTool.h:208
HIMonitoringPhotonsTool::m_h_photon_fcal_etcone20_ptCut
TH2D_LW * m_h_photon_fcal_etcone20_ptCut
Definition: HIMonitoringPhotonsTool.h:187
HIMonitoringPhotonsTool::m_h_photon_etcone30_ptCut_fcal1_all
TH1D * m_h_photon_etcone30_ptCut_fcal1_all
Definition: HIMonitoringPhotonsTool.h:211
HIMonitoringPhotonsTool::m_photonTightIsEMSelector
AsgPhotonIsEMSelector * m_photonTightIsEMSelector
Definition: HIMonitoringPhotonsTool.h:30
HIMonitoringPhotonsTool::m_h_photon_pt_etcone40_fcal0
TH2D_LW * m_h_photon_pt_etcone40_fcal0
Definition: HIMonitoringPhotonsTool.h:191
HIMonitoringPhotonsTool::m_h_photon_etcone30_ptCut_fcal0
TH1D * m_h_photon_etcone30_ptCut_fcal0
Definition: HIMonitoringPhotonsTool.h:199
HIMonitoringPhotonsTool::m_h_photon_f1_eta_ptCut
TH2D_LW * m_h_photon_f1_eta_ptCut
Definition: HIMonitoringPhotonsTool.h:164
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
HIMonitoringPhotonsTool::m_h_photon_pt_eta_tight
TH2D_LW * m_h_photon_pt_eta_tight
Definition: HIMonitoringPhotonsTool.h:106
HIMonitoringPhotonsTool::m_h_photon_etcone20_ptCut_fcal1
TH1D * m_h_photon_etcone20_ptCut_fcal1
Definition: HIMonitoringPhotonsTool.h:196
HIMonitoringPhotonsTool::m_h_photon_wtots1_ptCut_all
TH1D * m_h_photon_wtots1_ptCut_all
Definition: HIMonitoringPhotonsTool.h:119
HIMonitoringPhotonsTool::m_h_photon_wtots1_ptCut_tight
TH1D * m_h_photon_wtots1_ptCut_tight
Definition: HIMonitoringPhotonsTool.h:145
HIMonitoringPhotonsTool::m_h_photon_weta2_ptCut_loose
TH1D * m_h_photon_weta2_ptCut_loose
Definition: HIMonitoringPhotonsTool.h:129
HIMonitoringPhotonsTool::m_h_photon_deltae_ptCut_all
TH1D * m_h_photon_deltae_ptCut_all
Definition: HIMonitoringPhotonsTool.h:122
lumiFormat.i
int i
Definition: lumiFormat.py:92
HIMonitoringPhotonsTool::m_h_photon_wtots1_fcal_ptCut_eta0
TH2D_LW * m_h_photon_wtots1_fcal_ptCut_eta0
Definition: HIMonitoringPhotonsTool.h:174
AsgPhotonIsEMSelector::accept
virtual asg::AcceptData accept(const xAOD::IParticle *part) const
Accept with generic interface.
Definition: AsgPhotonIsEMSelector.cxx:227
HIMonitoringPhotonsTool::m_fracs1_nbins
int m_fracs1_nbins
Definition: HIMonitoringPhotonsTool.h:74
HIMonitoringPhotonsTool::m_h_photon_etcone30_ptCut_fcal2
TH1D * m_h_photon_etcone30_ptCut_fcal2
Definition: HIMonitoringPhotonsTool.h:201
xAOD::Iso::etcone20
@ etcone20
Calorimeter isolation.
Definition: IsolationType.h:32
TH2D_LW::Fill
virtual void Fill(const double &x, const double &y) override
Definition: TH2D_LW.cxx:329
RCU::Shell
Definition: ShellExec.cxx:28
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
HIMonitoringPhotonsTool::m_h_photon_etcone20_ptCut_fcal2
TH1D * m_h_photon_etcone20_ptCut_fcal2
Definition: HIMonitoringPhotonsTool.h:197
AthenaMonManager::online
@ online
Definition: AthenaMonManager.h:49
HIMonitoringPhotonsTool::m_h_photon_weta2_fcal_ptCut_eta0
TH2D_LW * m_h_photon_weta2_fcal_ptCut_eta0
Definition: HIMonitoringPhotonsTool.h:170
HIMonitoringPhotonsTool::m_h_photon_pt_phi
TH2D_LW * m_h_photon_pt_phi
Definition: HIMonitoringPhotonsTool.h:102
HIMonitoringPhotonsTool::m_h_photon_pt_fcal_tight
TH2D_LW * m_h_photon_pt_fcal_tight
Definition: HIMonitoringPhotonsTool.h:108
HIEventShapeContainer.h
HIMonitoringPhotonsTool::m_h_photon_etcone30_ptCut_fcal1_loose
TH1D * m_h_photon_etcone30_ptCut_fcal1_loose
Definition: HIMonitoringPhotonsTool.h:212
test_pyathena.parent
parent
Definition: test_pyathena.py:15
AthenaMonManager.h
HIMonitoringPhotonsTool::m_h_photon_fracs1_ptCut_loose
TH1D * m_h_photon_fracs1_ptCut_loose
Definition: HIMonitoringPhotonsTool.h:133
run
Definition: run.py:1
HIMonitoringPhotonsTool::m_f1_low
float m_f1_low
Definition: HIMonitoringPhotonsTool.h:87
HIMonitoringPhotonsTool::m_rhad_low
float m_rhad_low
Definition: HIMonitoringPhotonsTool.h:67
DataVector
Derived DataVector<T>.
Definition: DataVector.h:581
AcceptData.h
HIMonitoringPhotonsTool::m_h_photon_deltae_eta_ptCut
TH2D_LW * m_h_photon_deltae_eta_ptCut
Definition: HIMonitoringPhotonsTool.h:162
ManagedMonitorToolBase::m_environment
AthenaMonManager::Environment_t m_environment
Definition: ManagedMonitorToolBase.h:902
HIMonitoringPhotonsTool::m_h_photon_deltae_ptCut_tight
TH1D * m_h_photon_deltae_ptCut_tight
Definition: HIMonitoringPhotonsTool.h:148
HIMonitoringPhotonsTool::m_h_photon_reta_fcal_ptCut_eta1
TH2D_LW * m_h_photon_reta_fcal_ptCut_eta1
Definition: HIMonitoringPhotonsTool.h:167
HIMonitoringPhotonsTool::m_h_photon_etcone30_ptCut_fcal0_all
TH1D * m_h_photon_etcone30_ptCut_fcal0_all
Definition: HIMonitoringPhotonsTool.h:207
HIMonitoringPhotonsTool::HIMonitoringPhotonsTool
HIMonitoringPhotonsTool(const std::string &type, const std::string &name, const IInterface *parent)
Definition: HIMonitoringPhotonsTool.cxx:22
HIMonitoringPhotonsTool::m_h_photon_etcone30_ptCut_fcal1
TH1D * m_h_photon_etcone30_ptCut_fcal1
Definition: HIMonitoringPhotonsTool.h:200
HIMonitoringPhotonsTool::m_h_photon_pt_phi_tight
TH2D_LW * m_h_photon_pt_phi_tight
Definition: HIMonitoringPhotonsTool.h:107
HIMonitoringPhotonsTool::bookHistogramsRecurrent
virtual StatusCode bookHistogramsRecurrent()
An inheriting class should either override this function, bookHists() or bookHistograms().
Definition: HIMonitoringPhotonsTool.cxx:118
HIMonitoringPhotonsTool::m_fracs1_low
float m_fracs1_low
Definition: HIMonitoringPhotonsTool.h:75
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:195
HIMonitoringPhotonsTool::m_h_photon_rphi_eta_ptCut
TH2D_LW * m_h_photon_rphi_eta_ptCut
Definition: HIMonitoringPhotonsTool.h:155
HIMonitoringPhotonsTool::m_deltae_low
float m_deltae_low
Definition: HIMonitoringPhotonsTool.h:79
HIMonitoringPhotonsTool::m_IsoNarrow_nbins
int m_IsoNarrow_nbins
Definition: HIMonitoringPhotonsTool.h:94
HIMonitoringPhotonsTool::~HIMonitoringPhotonsTool
virtual ~HIMonitoringPhotonsTool()
Definition: HIMonitoringPhotonsTool.cxx:114
HIMonitoringPhotonsTool::m_h_photon_weta2_ptCut_tight
TH1D * m_h_photon_weta2_ptCut_tight
Definition: HIMonitoringPhotonsTool.h:142
HIMonitoringPhotonsTool::m_h_photon_etcone30_ptCut_fcal1_tight
TH1D * m_h_photon_etcone30_ptCut_fcal1_tight
Definition: HIMonitoringPhotonsTool.h:213
HIMonitoringPhotonsTool::m_h_photon_rphi_ptCut_tight
TH1D * m_h_photon_rphi_ptCut_tight
Definition: HIMonitoringPhotonsTool.h:141
HIMonitoringPhotonsTool::m_h_photon_reta_eta_ptCut
TH2D_LW * m_h_photon_reta_eta_ptCut
Definition: HIMonitoringPhotonsTool.h:154
HIMonitoringPhotonsTool::m_h_photon_weta2_ptCut_all
TH1D * m_h_photon_weta2_ptCut_all
Definition: HIMonitoringPhotonsTool.h:116
HIMonitoringPhotonsTool::m_IsoWide_high
float m_IsoWide_high
Definition: HIMonitoringPhotonsTool.h:92
HIMonitoringPhotonsTool::m_h_photon_etcone40_ptCut_fcal2
TH1D * m_h_photon_etcone40_ptCut_fcal2
Definition: HIMonitoringPhotonsTool.h:205
HIMonitoringPhotonsTool::m_weta2_nbins
int m_weta2_nbins
Definition: HIMonitoringPhotonsTool.h:62
HIMonitoringPhotonsTool::m_PhotonEta_low
float m_PhotonEta_low
Definition: HIMonitoringPhotonsTool.h:47
HIMonitoringPhotonsTool::m_PhotonPt_high
float m_PhotonPt_high
Definition: HIMonitoringPhotonsTool.h:44
HIMonitoringPhotonsTool::m_h_photon_etcone30_ptCut_fcal2_loose
TH1D * m_h_photon_etcone30_ptCut_fcal2_loose
Definition: HIMonitoringPhotonsTool.h:216
HIMonitoringPhotonsTool::m_PhotonPhi_nbins
int m_PhotonPhi_nbins
Definition: HIMonitoringPhotonsTool.h:50
HIMonitoringPhotonsTool::m_h_photon_fracs1_eta_ptCut
TH2D_LW * m_h_photon_fracs1_eta_ptCut
Definition: HIMonitoringPhotonsTool.h:160
HIMonitoringPhotonsTool::m_IsoWide_nbins
int m_IsoWide_nbins
Definition: HIMonitoringPhotonsTool.h:90
HIMonitoringPhotonsTool::m_h_photon_rhad_ptCut_tight
TH1D * m_h_photon_rhad_ptCut_tight
Definition: HIMonitoringPhotonsTool.h:144
HIMonitoringPhotonsTool::m_eratio_nbins
int m_eratio_nbins
Definition: HIMonitoringPhotonsTool.h:82
DataVector::end
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
AsgPhotonIsEMSelector.h
HIMonitoringPhotonsTool::m_h_photon_eratio_ptCut_tight
TH1D * m_h_photon_eratio_ptCut_tight
Definition: HIMonitoringPhotonsTool.h:149
DiTauMassTools::MaxHistStrategyV2::e
e
Definition: PhysicsAnalysis/TauID/DiTauMassTools/DiTauMassTools/HelperFunctions.h:26
HIMonitoringPhotonsTool::m_h_photon_rhad_fcal_ptCut_eta0
TH2D_LW * m_h_photon_rhad_fcal_ptCut_eta0
Definition: HIMonitoringPhotonsTool.h:172
HIMonitoringPhotonsTool::m_h_photon_fcal_etcone40_ptCut
TH2D_LW * m_h_photon_fcal_etcone40_ptCut
Definition: HIMonitoringPhotonsTool.h:189
ManagedMonitorToolBase::endOfRunFlag
bool endOfRunFlag() const
Definition: ManagedMonitorToolBase.h:859
HIMonitoringPhotonsTool::m_reta_high
float m_reta_high
Definition: HIMonitoringPhotonsTool.h:56
HIMonitoringPhotonsTool::m_h_photon_eta_phi_ptCut
TH2D_LW * m_h_photon_eta_phi_ptCut
Definition: HIMonitoringPhotonsTool.h:104
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
HIMonitoringPhotonsTool::m_h_photon_reta_fcal_ptCut_eta0
TH2D_LW * m_h_photon_reta_fcal_ptCut_eta0
Definition: HIMonitoringPhotonsTool.h:166
HIMonitoringPhotonsTool::m_h_photon_eratio_fcal_ptCut_eta1
TH2D_LW * m_h_photon_eratio_fcal_ptCut_eta1
Definition: HIMonitoringPhotonsTool.h:181
HIMonitoringPhotonsTool::procHistograms
virtual StatusCode procHistograms()
An inheriting class should either override this function or finalHists().
Definition: HIMonitoringPhotonsTool.cxx:345
HIMonitoringPhotonsTool::m_f1_high
float m_f1_high
Definition: HIMonitoringPhotonsTool.h:88
HIMonitoringPhotonsTool::m_h_photon_fracs1_fcal_ptCut_eta1
TH2D_LW * m_h_photon_fracs1_fcal_ptCut_eta1
Definition: HIMonitoringPhotonsTool.h:177
HIMonitoringPhotonsTool::m_h_photon_fracs1_fcal_ptCut_eta0
TH2D_LW * m_h_photon_fracs1_fcal_ptCut_eta0
Definition: HIMonitoringPhotonsTool.h:176
HIMonitoringPhotonsTool::m_FCalEt_central
float m_FCalEt_central
Definition: HIMonitoringPhotonsTool.h:40
HIMonitoringPhotonsTool::m_h_photon_wtots1_fcal_ptCut_eta1
TH2D_LW * m_h_photon_wtots1_fcal_ptCut_eta1
Definition: HIMonitoringPhotonsTool.h:175
HIMonitoringPhotonsTool::m_h_photon_rhad_ptCut_loose
TH1D * m_h_photon_rhad_ptCut_loose
Definition: HIMonitoringPhotonsTool.h:131
HIMonitoringPhotonsTool::m_rhad_high
float m_rhad_high
Definition: HIMonitoringPhotonsTool.h:68
HIMonitoringPhotonsTool::m_eratio_low
float m_eratio_low
Definition: HIMonitoringPhotonsTool.h:83
HIMonitoringPhotonsTool::fillHistograms
virtual StatusCode fillHistograms()
An inheriting class should either override this function or fillHists().
Definition: HIMonitoringPhotonsTool.cxx:138
HIMonitoringPhotonsTool::m_h_photon_eratio_fcal_ptCut_eta0
TH2D_LW * m_h_photon_eratio_fcal_ptCut_eta0
Definition: HIMonitoringPhotonsTool.h:180
HIMonitoringPhotonsTool::m_h_photon_eratio_eta_ptCut
TH2D_LW * m_h_photon_eratio_eta_ptCut
Definition: HIMonitoringPhotonsTool.h:163
DataVector::at
const T * at(size_type n) const
Access an element, as an rvalue.
ConstAccessor.h
Helper class to provide constant type-safe access to aux data.
HIMonitoringPhotonsTool::m_h_photon_rhad_fcal_ptCut_eta1
TH2D_LW * m_h_photon_rhad_fcal_ptCut_eta1
Definition: HIMonitoringPhotonsTool.h:173
HIMonitoringPhotonsTool::m_wtots1_low
float m_wtots1_low
Definition: HIMonitoringPhotonsTool.h:71
HIMonitoringPhotonsTool::m_h_photon_etcone40_ptCut_fcal0
TH1D * m_h_photon_etcone40_ptCut_fcal0
Definition: HIMonitoringPhotonsTool.h:203
HIMonitoringPhotonsTool::m_h_photon_pt_etcone40_fcal2
TH2D_LW * m_h_photon_pt_etcone40_fcal2
Definition: HIMonitoringPhotonsTool.h:193
HIMonitoringPhotonsTool::m_h_photon_rhad_eta_ptCut
TH2D_LW * m_h_photon_rhad_eta_ptCut
Definition: HIMonitoringPhotonsTool.h:158
HIMonitoringPhotonsTool::m_h_photon_deltae_ptCut_loose
TH1D * m_h_photon_deltae_ptCut_loose
Definition: HIMonitoringPhotonsTool.h:135
xAOD::bool
setBGCode setTAP setLVL2ErrorBits bool
Definition: TrigDecision_v1.cxx:60
TH2D_LW::create
static TH2D_LW * create(const char *name, const char *title, unsigned nbinsx, const double &xmin, const double &xmax, unsigned nbinsy, const double &ymin, const double &ymax)
Definition: TH2D_LW.cxx:33
HIMonitoringPhotonsTool::m_h_photon_f1_ptCut_all
TH1D * m_h_photon_f1_ptCut_all
Definition: HIMonitoringPhotonsTool.h:124
HIMonitoringPhotonsTool::m_FCalEt_peripheral
float m_FCalEt_peripheral
Definition: HIMonitoringPhotonsTool.h:39
DataVector::size
size_type size() const noexcept
Returns the number of elements in the collection.
HIMonitoringPhotonsTool::m_wtots1_high
float m_wtots1_high
Definition: HIMonitoringPhotonsTool.h:72
PhotonContainer.h
HIMonitoringPhotonsTool::m_PhotonEta_high
float m_PhotonEta_high
Definition: HIMonitoringPhotonsTool.h:48
HIMonitoringPhotonsTool::m_FCalEt_high
float m_FCalEt_high
Definition: HIMonitoringPhotonsTool.h:37
ManagedMonitorToolBase::regHist
virtual StatusCode regHist(TH1 *h, const std::string &system, Interval_t interval, MgmtAttr_t histo_mgmt=ATTRIB_MANAGED, const std::string &chain="", const std::string &merge="")
Registers a TH1 (including TH2, TH3, and TProfile) to be included in the output stream using logical ...
Definition: ManagedMonitorToolBase.cxx:1454
HIMonitoringPhotonsTool::m_h_photon_rphi_fcal_ptCut_eta1
TH2D_LW * m_h_photon_rphi_fcal_ptCut_eta1
Definition: HIMonitoringPhotonsTool.h:169
HIMonitoringPhotonsTool::m_IsoNarrow_high
float m_IsoNarrow_high
Definition: HIMonitoringPhotonsTool.h:96
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
DataVector::begin
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
HIMonitoringPhotonsTool::book_hist
void book_hist()
*** FCal sum Et *** ///
Definition: HIMonitoringPhotonsTool.cxx:404
HIMonitoringPhotonsTool::m_h_photon_eta_phi_ptCut_tight
TH2D_LW * m_h_photon_eta_phi_ptCut_tight
Definition: HIMonitoringPhotonsTool.h:109
HIMonitoringPhotonsTool::m_h_photon_fcal_etcone30_ptCut
TH2D_LW * m_h_photon_fcal_etcone30_ptCut
Definition: HIMonitoringPhotonsTool.h:188
HIMonitoringPhotonsTool::m_h_photon_rphi_ptCut_loose
TH1D * m_h_photon_rphi_ptCut_loose
Definition: HIMonitoringPhotonsTool.h:128
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
HIMonitoringPhotonsTool::m_PhotonPt_low
float m_PhotonPt_low
Definition: HIMonitoringPhotonsTool.h:43
HIMonitoringPhotonsTool::m_h_photon_rphi_ptCut_all
TH1D * m_h_photon_rphi_ptCut_all
Definition: HIMonitoringPhotonsTool.h:115
HIMonitoringPhotonsTool::m_h_photon_etcone30_ptCut_fcal2_tight
TH1D * m_h_photon_etcone30_ptCut_fcal2_tight
Definition: HIMonitoringPhotonsTool.h:217
AsgPhotonIsEMSelector
Photon isEM selector.
Definition: AsgPhotonIsEMSelector.h:41
HIMonitoringPhotonsTool::m_weta2_low
float m_weta2_low
Definition: HIMonitoringPhotonsTool.h:63
Generated on Tue Jul 2 2024 21:16:44 for ATLAS Offline Software by   doxygen 1.8.18