ATLAS Offline Software
HIJetUEMonitoring.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 //HIJetUEMonitoring.cxx
6 
7 #include "JetMonitoring/HIJetUEMonitoring.h"
8 #include "JetMonitoring/ToolHandleHistoHelper.h"
9 #include "xAODEventInfo/EventInfo.h"
10 #include "AthContainers/ConstAccessor.h"
11 
12 #include "TH1.h"
13 #include "TH2.h"
14 #include "TProfile.h"
15 
16 #define toGeV 1/1000.
17 #define toTeV 1/1000000.
18 
19 HIJetUEMonitoring::HIJetUEMonitoring(const std::string& name)
20  : JetHistoBase(name), m_histoDef(this), m_jetScale(xAOD::JetAssignedScaleMomentum) {
21  declareProperty("HIEventShapeContainerName", m_container_key="CaloSums");
22  declareProperty("HistoDef", m_histoDef, "The list of HistoDefinitionTool defining the histos to be used in this tool");
23  declareProperty("RefContainer", m_refContainerName);
24 }
25 
26 StatusCode HIJetUEMonitoring::initialize() {
27  return StatusCode::SUCCESS;
28 }
29 
30 int HIJetUEMonitoring::buildHistos(){
31  ATH_MSG_INFO(" buildHistos num of histos : "<< m_histoDef.size() );
32 
33  ToolHandleHistoHelper::HistoRetriever hbuilder(m_histoDef);
34 
35  // Histos are defined in jobOptions !
36  // For each histo, ask hbuilder if a corresponding definition exists in the jobOption list.
37  // -> if so a valid histo is returned (and booked)
38  // -> else NULL is returned
39  // m_peripheral = bookHisto( hbuilder.build<TH1F>("peripheral") );
40 
41  m_FCALET = bookHisto( hbuilder.build<TH1F>("FCalET") );
42 
43  //inclusive
44  m_2dSubtractedET_pT = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_pT") );
45  m_SubtractedET_Centrality = bookHisto( hbuilder.build<TProfile>("SubtractedET_Centrality") );
46  m_2dSubtractedET_Centrality = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_Centrality") );
47  m_2dSubtractedET_Expected_Centrality = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_Expected_Centrality") );
48  m_2dSubtractedET_Expected_eta = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_Expected_eta") );
49  m_SubtractedET_Expected_eta = bookHisto( hbuilder.build<TProfile>("SubtractedET_Expected_eta") );
50  m_2dSubtractedET_2Dphi = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_2Dphi") );
51  m_SubtractedET_eta = bookHisto( hbuilder.build<TProfile>("SubtractedET_eta") );
52  m_SubtractedET_pt = bookHisto( hbuilder.build<TProfile>("SubtractedET_pt") );
53  m_SubtractedET_2Dphi = bookHisto( hbuilder.build<TProfile>("SubtractedET_2Dphi") );
54 
55  //0-10%
56  m_JetUnsubtractedScaleMomentum_pt_0_10 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_pt_0_10") );
57  m_JetUnsubtractedScaleMomentum_eta_0_10 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_eta_0_10") );
58  m_JetUnsubtractedScaleMomentum_phi_0_10 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_phi_0_10") );
59  m_JetUnsubtractedScaleMomentum_m_0_10 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_m_0_10") );
60  m_JetSubtractedScaleMomentum_pt_0_10 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_pt_0_10") );
61  m_JetSubtractedScaleMomentum_eta_0_10 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_eta_0_10") );
62  m_JetSubtractedScaleMomentum_phi_0_10 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_phi_0_10") );
63  m_JetSubtractedScaleMomentum_m_0_10 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_m_0_10") );
64 
65  m_2dSubtractedET_Expected_eta_0_10 = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_Expected_eta_0_10") );
66  m_SubtractedET_Expected_eta_0_10 = bookHisto( hbuilder.build<TProfile>("SubtractedET_Expected_eta_0_10") );
67 
68  m_2dSubtractedET_2Dphi_0_10 = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_2Dphi_0_10") );
69  m_SubtractedET_eta_0_10 = bookHisto( hbuilder.build<TProfile>("SubtractedET_eta_0_10") );
70  m_SubtractedET_pt_0_10 = bookHisto( hbuilder.build<TProfile>("SubtractedET_pt_0_10") );
71 
72  m_SubtractedET_2Dphi_0_10 = bookHisto( hbuilder.build<TProfile>("SubtractedET_2Dphi_0_10") );
73 
74  // // //10-20%
75  m_JetUnsubtractedScaleMomentum_pt_10_20 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_pt_10_20") );
76  m_JetUnsubtractedScaleMomentum_eta_10_20 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_eta_10_20") );
77  m_JetUnsubtractedScaleMomentum_phi_10_20 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_phi_10_20") );
78  m_JetUnsubtractedScaleMomentum_m_10_20 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_m_10_20") );
79  m_JetSubtractedScaleMomentum_pt_10_20 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_pt_10_20") );
80  m_JetSubtractedScaleMomentum_eta_10_20 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_eta_10_20") );
81  m_JetSubtractedScaleMomentum_phi_10_20 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_phi_10_20") );
82  m_JetSubtractedScaleMomentum_m_10_20 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_m_10_20") );
83 
84  m_2dSubtractedET_Expected_eta_10_20 = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_Expected_eta_10_20") );
85  m_SubtractedET_Expected_eta_10_20 = bookHisto( hbuilder.build<TProfile>("SubtractedET_Expected_eta_10_20") );
86 
87  m_2dSubtractedET_2Dphi_10_20 = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_2Dphi_10_20") );
88 
89  m_SubtractedET_eta_10_20 = bookHisto( hbuilder.build<TProfile>("SubtractedET_eta_10_20") );
90  m_SubtractedET_pt_10_20 = bookHisto( hbuilder.build<TProfile>("SubtractedET_pt_10_20") );
91 
92  m_SubtractedET_2Dphi_10_20 = bookHisto( hbuilder.build<TProfile>("SubtractedET_2Dphi_10_20") );
93 
94 
95  // //20-40%
96  m_JetUnsubtractedScaleMomentum_pt_20_40 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_pt_20_40") );
97  m_JetUnsubtractedScaleMomentum_eta_20_40 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_eta_20_40") );
98  m_JetUnsubtractedScaleMomentum_phi_20_40 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_phi_20_40") );
99  m_JetUnsubtractedScaleMomentum_m_20_40 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_m_20_40") );
100  m_JetSubtractedScaleMomentum_pt_20_40 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_pt_20_40") );
101  m_JetSubtractedScaleMomentum_eta_20_40 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_eta_20_40") );
102  m_JetSubtractedScaleMomentum_phi_20_40 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_phi_20_40") );
103  m_JetSubtractedScaleMomentum_m_20_40 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_m_20_40") );
104 
105  m_2dSubtractedET_Expected_eta_20_40 = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_Expected_eta_20_40") );
106  m_SubtractedET_Expected_eta_20_40 = bookHisto( hbuilder.build<TProfile>("SubtractedET_Expected_eta_20_40") );
107 
108  m_2dSubtractedET_2Dphi_20_40 = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_2Dphi_20_40") );
109 
110  m_SubtractedET_eta_20_40 = bookHisto( hbuilder.build<TProfile>("SubtractedET_eta_20_40") );
111  m_SubtractedET_pt_20_40 = bookHisto( hbuilder.build<TProfile>("SubtractedET_pt_20_40") );
112 
113  m_SubtractedET_2Dphi_20_40 = bookHisto( hbuilder.build<TProfile>("SubtractedET_2Dphi_20_40") );
114 
115  // //60-100%
116  m_JetUnsubtractedScaleMomentum_pt_60_100 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_pt_60_100") );
117  m_JetUnsubtractedScaleMomentum_eta_60_100 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_eta_60_100") );
118  m_JetUnsubtractedScaleMomentum_phi_60_100 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_phi_60_100") );
119  m_JetUnsubtractedScaleMomentum_m_60_100 = bookHisto( hbuilder.build<TH1F>("JetUnsubtractedScaleMomentum_m_60_100") );
120  m_JetSubtractedScaleMomentum_pt_60_100 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_pt_60_100") );
121  m_JetSubtractedScaleMomentum_eta_60_100 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_eta_60_100") );
122  m_JetSubtractedScaleMomentum_phi_60_100 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_phi_60_100") );
123  m_JetSubtractedScaleMomentum_m_60_100 = bookHisto( hbuilder.build<TH1F>("JetSubtractedScaleMomentum_m_60_100") );
124 
125  m_2dSubtractedET_Expected_eta_60_100 = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_Expected_eta_60_100") );
126  m_SubtractedET_Expected_eta_60_100 = bookHisto( hbuilder.build<TProfile>("SubtractedET_Expected_eta_60_100") );
127 
128  m_2dSubtractedET_2Dphi_60_100 = bookHisto( hbuilder.build<TH2F>("2dSubtractedET_2Dphi_60_100") );
129 
130  m_SubtractedET_eta_60_100 = bookHisto( hbuilder.build<TProfile>("SubtractedET_eta_60_100") );
131  m_SubtractedET_pt_60_100 = bookHisto( hbuilder.build<TProfile>("SubtractedET_pt_60_100") );
132 
133  m_SubtractedET_2Dphi_60_100 = bookHisto( hbuilder.build<TProfile>("SubtractedET_2Dphi_60_100") );
134 
135  return 0;
136 }
137 
138 int HIJetUEMonitoring::fillHistosFromJet(const xAOD::Jet &j, float weight){
139 
140  const xAOD::EventInfo* evtInfo = nullptr;
141  CHECK(evtStore()->retrieve( evtInfo ), 1);
142 
143 //LAr event veto: skip events rejected by LAr
144  if(evtInfo->errorState(xAOD::EventInfo::LAr)==xAOD::EventInfo::Error){
145  ATH_MSG_DEBUG("SKIP for LAR error");
146  return 1;
147  }
148 
149  m_n=2;
150  m_harmonic=m_n-1;
151  m_eventShape=nullptr;
152  CHECK( evtStore()->retrieve(m_eventShape,m_container_key), 1 );
153  m_FCalET=0;
154  m_psiN_FCal=0;
155  // m_vN_fcal=0;
156  for(const xAOD::HIEventShape* sh : *m_eventShape){
157  static const SG::ConstAccessor<std::string> SummaryAcc("Summary");
158  std::string summary = SummaryAcc.withDefault(*sh, "");
159  if(summary.compare("FCal")==0){
160  m_FCalET=sh->et()*toTeV;
161  float qx=sh->etCos().at(m_harmonic);
162  float qy=sh->etSin().at(m_harmonic);
163  m_psiN_FCal=std::atan2(qy,qx)/float(m_n); // Psi2 (m_n=2)
164  m_vN_fcal=-999.;
165  if (m_FCalET != 0.) m_vN_fcal = std::sqrt(qx+qx+qy*qy)/m_FCalET;
166  break;
167  }
168  }
169 
170  m_FCALET->Fill(m_FCalET, weight);
171 
172  // m_jetScale is a property of the base tool
173  //const xAOD::JetFourMom_t p4 = j.jetP4( (xAOD::JetScale) m_jetScale);
174 
175  float SubtractedET = j.getAttribute<float>("JetUnsubtractedScaleMomentum_pt")*toGeV - j.getAttribute<float>("JetSubtractedScaleMomentum_pt")*toGeV;
176 
177  float Acos = std::acos(std::cos(2*(j.getAttribute<float>("JetEtaJESScaleMomentum_phi") - m_psiN_FCal)));
178  // float diff = fabs(j.getAttribute<float>("JetSubtractedScaleMomentum_phi") - m_psiN_FCal);
179  // while (diff > TMath::Pi()/2. ) diff = TMath::Pi() - diff;
180 
181  if (fabs(j.getAttribute<float>("JetEtaJESScaleMomentum_eta"))<2.8){
182  if (m_FCalET > 2.7){//0-10%
183  m_SubtractedET_pt_0_10->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_pt")*toGeV,SubtractedET, weight );
184  }
185  if (m_FCalET < 2.7 && m_FCalET > 1.75 ){//10-20%
186  m_SubtractedET_pt_10_20->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_pt")*toGeV,SubtractedET, weight );
187  }
188  if (m_FCalET < 1.75 && m_FCalET > 0.65 ){//20-40%
189  m_SubtractedET_pt_20_40->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_pt")*toGeV,SubtractedET, weight );
190  }
191  if (m_FCalET < 0.20 ){//60-100%
192  m_SubtractedET_pt_60_100->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_pt")*toGeV,SubtractedET, weight );
193  }
194  }
195 
196  if (j.getAttribute<float>("JetEtaJESScaleMomentum_pt")*toGeV > m_ptcut){
197  // ATH_MSG_INFO(" FCal ET: "<< m_FCalET<<" SubtractedE:, "<<SubtractedET<<" m_psiN_FCal: "<<m_psiN_FCal<<" Acos: "<<Acos <<" m_vN_fcal: "<<m_vN_fcal<<" ptcut: "<<m_ptcut);
198  m_2dSubtractedET_pT->Fill( j.getAttribute<float>("JetEtaJESScaleMomentum_pt")*toGeV, SubtractedET, weight);
199 
200  double subtractedET_Expected = -999.;
201  if (m_FCalET != 0.) subtractedET_Expected = (SubtractedET/m_FCalET)*0.025;
202 
203  m_SubtractedET_Centrality->Fill( m_FCalET, SubtractedET, weight);
204  m_2dSubtractedET_Centrality->Fill( m_FCalET, SubtractedET, weight);
205  m_2dSubtractedET_Expected_Centrality->Fill( m_FCalET, subtractedET_Expected, weight);
206 
207  m_2dSubtractedET_Expected_eta->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta") , subtractedET_Expected, weight);
208  m_SubtractedET_Expected_eta->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta") , subtractedET_Expected, weight);
209  m_2dSubtractedET_2Dphi->Fill( Acos,SubtractedET, weight);
210  m_SubtractedET_eta->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta"),SubtractedET, weight);
211  m_SubtractedET_pt->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_pt")*toGeV,SubtractedET, weight );
212  m_SubtractedET_2Dphi->Fill( Acos,SubtractedET, weight);
213 
214  if (m_FCalET > 2.7){//0-10%
215  m_JetUnsubtractedScaleMomentum_pt_0_10->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_pt")*toGeV, weight );
216  m_JetUnsubtractedScaleMomentum_eta_0_10->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_eta"), weight );
217  m_JetUnsubtractedScaleMomentum_phi_0_10->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_phi"), weight );
218  m_JetUnsubtractedScaleMomentum_m_0_10->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_m")*toGeV, weight );
219  m_JetSubtractedScaleMomentum_pt_0_10->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_pt")*toGeV, weight );
220  m_JetSubtractedScaleMomentum_eta_0_10->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_eta"), weight );
221  m_JetSubtractedScaleMomentum_phi_0_10->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_phi"), weight );
222  m_JetSubtractedScaleMomentum_m_0_10->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_m")*toGeV, weight );
223  m_2dSubtractedET_Expected_eta_0_10->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta") , subtractedET_Expected, weight);
224  m_SubtractedET_Expected_eta_0_10->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta") , subtractedET_Expected, weight);
225 
226  m_2dSubtractedET_2Dphi_0_10->Fill( Acos,SubtractedET, weight);
227  m_SubtractedET_eta_0_10->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta"),SubtractedET, weight );
228  m_SubtractedET_2Dphi_0_10->Fill( Acos,SubtractedET, weight);
229  }
230  if (m_FCalET < 2.7 && m_FCalET > 1.75 ){//10-20%
231  m_JetUnsubtractedScaleMomentum_pt_10_20->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_pt")*toGeV, weight );
232  m_JetUnsubtractedScaleMomentum_eta_10_20->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_eta"), weight );
233  m_JetUnsubtractedScaleMomentum_phi_10_20->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_phi"), weight );
234  m_JetUnsubtractedScaleMomentum_m_10_20->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_m")*toGeV, weight );
235  m_JetSubtractedScaleMomentum_pt_10_20->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_pt")*toGeV, weight );
236  m_JetSubtractedScaleMomentum_eta_10_20->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_eta"), weight );
237  m_JetSubtractedScaleMomentum_phi_10_20->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_phi"), weight );
238  m_JetSubtractedScaleMomentum_m_10_20->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_m")*toGeV, weight );
239 
240  m_2dSubtractedET_Expected_eta_10_20->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta") , subtractedET_Expected, weight);
241  m_SubtractedET_Expected_eta_10_20->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta") , subtractedET_Expected, weight);
242 
243  m_2dSubtractedET_2Dphi_10_20->Fill( Acos,SubtractedET, weight);
244 
245  m_SubtractedET_eta_10_20->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta"),SubtractedET, weight );
246  m_SubtractedET_2Dphi_10_20->Fill( Acos,SubtractedET, weight);
247  }
248 
249  if (m_FCalET < 1.75 && m_FCalET > 0.65 ){//20-40%
250  m_JetUnsubtractedScaleMomentum_pt_20_40->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_pt")*toGeV, weight );
251  m_JetUnsubtractedScaleMomentum_eta_20_40->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_eta"), weight );
252  m_JetUnsubtractedScaleMomentum_phi_20_40->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_phi"), weight );
253  m_JetUnsubtractedScaleMomentum_m_20_40->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_m")*toGeV, weight );
254  m_JetSubtractedScaleMomentum_pt_20_40->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_pt")*toGeV, weight );
255  m_JetSubtractedScaleMomentum_eta_20_40->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_eta"), weight );
256  m_JetSubtractedScaleMomentum_phi_20_40->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_phi"), weight );
257  m_JetSubtractedScaleMomentum_m_20_40->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_m")*toGeV, weight );
258 
259  m_2dSubtractedET_Expected_eta_20_40->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta") , subtractedET_Expected, weight);
260  m_SubtractedET_Expected_eta_20_40->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta") , subtractedET_Expected, weight);
261 
262  m_2dSubtractedET_2Dphi_20_40->Fill( Acos,SubtractedET, weight);
263 
264  m_SubtractedET_eta_20_40->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta"),SubtractedET, weight );
265  m_SubtractedET_2Dphi_20_40->Fill( Acos,SubtractedET, weight);
266  }
267  if (m_FCalET < 0.20 ){//60-100%
268  m_JetUnsubtractedScaleMomentum_pt_60_100->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_pt")*toGeV, weight );
269  m_JetUnsubtractedScaleMomentum_eta_60_100->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_eta"), weight );
270  m_JetUnsubtractedScaleMomentum_phi_60_100->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_phi"), weight );
271  m_JetUnsubtractedScaleMomentum_m_60_100->Fill( j.getAttribute<float>("JetUnsubtractedScaleMomentum_m")*toGeV, weight );
272  m_JetSubtractedScaleMomentum_pt_60_100->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_pt")*toGeV, weight );
273  m_JetSubtractedScaleMomentum_eta_60_100->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_eta"), weight );
274  m_JetSubtractedScaleMomentum_phi_60_100->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_phi"), weight );
275  m_JetSubtractedScaleMomentum_m_60_100->Fill( j.getAttribute<float>("JetSubtractedScaleMomentum_m")*toGeV, weight );
276 
277  m_2dSubtractedET_Expected_eta_60_100->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta") , subtractedET_Expected, weight);
278  m_SubtractedET_Expected_eta_60_100->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta") , subtractedET_Expected, weight);
279 
280  m_2dSubtractedET_2Dphi_60_100->Fill( Acos,SubtractedET, weight);
281 
282  m_SubtractedET_eta_60_100->Fill(j.getAttribute<float>("JetEtaJESScaleMomentum_eta"),SubtractedET, weight );
283  m_SubtractedET_2Dphi_60_100->Fill( Acos,SubtractedET, weight);
284  }
285 
286  }
287 
288  // ATH_MSG_INFO(" FCal ET: "<< m_FCalET<<" SubtractedE:, "<<SubtractedET<<" m_psiN_FCal: "<<m_psiN_FCal<<" Acos: "<<Acos <<" m_vN_fcal: "<<m_vN_fcal<<" ptcut: "<<m_ptcut);
289  //ATH_MSG_INFO(" Teste : "<<p4.Eta()<<" "<<j.getAttribute<float>("JetSubtractedScaleMomentum_eta"));
290  return 0;
291 }
292 
293 
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
HIJetUEMonitoring::m_SubtractedET_pt
TProfile * m_SubtractedET_pt
Definition: HIJetUEMonitoring.h:59
HIJetUEMonitoring::m_SubtractedET_eta
TProfile * m_SubtractedET_eta
Definition: HIJetUEMonitoring.h:58
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_pt_20_40
TH1F * m_JetUnsubtractedScaleMomentum_pt_20_40
Definition: HIJetUEMonitoring.h:102
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_phi_60_100
TH1F * m_JetUnsubtractedScaleMomentum_phi_60_100
Definition: HIJetUEMonitoring.h:123
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_m_10_20
TH1F * m_JetSubtractedScaleMomentum_m_10_20
Definition: HIJetUEMonitoring.h:90
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_m_20_40
TH1F * m_JetUnsubtractedScaleMomentum_m_20_40
Definition: HIJetUEMonitoring.h:105
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
HIJetUEMonitoring::m_vN_fcal
float m_vN_fcal
Definition: HIJetUEMonitoring.h:38
TH2F
Definition: rootspy.cxx:420
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_pt_0_10
TH1F * m_JetUnsubtractedScaleMomentum_pt_0_10
Definition: HIJetUEMonitoring.h:63
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_eta_20_40
TH1F * m_JetSubtractedScaleMomentum_eta_20_40
Definition: HIJetUEMonitoring.h:107
HIJetUEMonitoring::m_SubtractedET_pt_10_20
TProfile * m_SubtractedET_pt_10_20
Definition: HIJetUEMonitoring.h:98
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_pt_10_20
TH1F * m_JetUnsubtractedScaleMomentum_pt_10_20
Definition: HIJetUEMonitoring.h:83
HIJetUEMonitoring::m_FCalET
float m_FCalET
Centrality.
Definition: HIJetUEMonitoring.h:36
HIJetUEMonitoring::m_SubtractedET_Expected_eta
TProfile * m_SubtractedET_Expected_eta
Definition: HIJetUEMonitoring.h:56
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_m_0_10
TH1F * m_JetSubtractedScaleMomentum_m_0_10
Definition: HIJetUEMonitoring.h:70
HIJetUEMonitoring::m_SubtractedET_2Dphi
TProfile * m_SubtractedET_2Dphi
Definition: HIJetUEMonitoring.h:60
HIJetUEMonitoring::m_histoDef
ToolHandleArray< HistoDefinitionTool > m_histoDef
Definition: HIJetUEMonitoring.h:44
xAOD
ICaloAffectedTool is abstract interface for tools checking if 4 mom is in calo affected region.
Definition: ICaloAffectedTool.h:24
HIJetUEMonitoring::fillHistosFromJet
virtual int fillHistosFromJet(const xAOD::Jet &j, float weight)
Definition: HIJetUEMonitoring.cxx:138
SG::ConstAccessor
Helper class to provide constant type-safe access to aux data.
Definition: ConstAccessor.h:54
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_pt_10_20
TH1F * m_JetSubtractedScaleMomentum_pt_10_20
Definition: HIJetUEMonitoring.h:87
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
HIJetUEMonitoring::m_2dSubtractedET_Expected_eta_0_10
TH2F * m_2dSubtractedET_Expected_eta_0_10
Definition: HIJetUEMonitoring.h:72
HIJetUEMonitoring::m_SubtractedET_2Dphi_10_20
TProfile * m_SubtractedET_2Dphi_10_20
Definition: HIJetUEMonitoring.h:99
xAOD::EventInfo_v1::LAr
@ LAr
The LAr calorimeter.
Definition: EventInfo_v1.h:335
HIJetUEMonitoring::m_SubtractedET_pt_20_40
TProfile * m_SubtractedET_pt_20_40
Definition: HIJetUEMonitoring.h:117
xAOD::Jet_v1::getAttribute
bool getAttribute(AttributeID type, T &value) const
Retrieve attribute moment by enum.
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_eta_0_10
TH1F * m_JetSubtractedScaleMomentum_eta_0_10
Definition: HIJetUEMonitoring.h:68
HIJetUEMonitoring::m_2dSubtractedET_2Dphi_10_20
TH2F * m_2dSubtractedET_2Dphi_10_20
Definition: HIJetUEMonitoring.h:95
HIJetUEMonitoring::m_2dSubtractedET_2Dphi_60_100
TH2F * m_2dSubtractedET_2Dphi_60_100
Definition: HIJetUEMonitoring.h:133
HIJetUEMonitoring::m_psiN_FCal
float m_psiN_FCal
Definition: HIJetUEMonitoring.h:37
xAOD::EventInfo_v1::Error
@ Error
The sub-detector issued an error.
Definition: EventInfo_v1.h:349
HIJetUEMonitoring::m_n
unsigned int m_n
Definition: HIJetUEMonitoring.h:39
dqt_zlumi_pandas.weight
int weight
Definition: dqt_zlumi_pandas.py:200
HIJetUEMonitoring::m_2dSubtractedET_pT
TH2F * m_2dSubtractedET_pT
Definition: HIJetUEMonitoring.h:53
HIJetUEMonitoring::initialize
virtual StatusCode initialize()
Dummy implementation of the initialisation function.
Definition: HIJetUEMonitoring.cxx:26
JetHistoBase
Definition: JetHistoBase.h:27
HIJetUEMonitoring::m_SubtractedET_Centrality
TProfile * m_SubtractedET_Centrality
Definition: HIJetUEMonitoring.h:50
HIJetUEMonitoring::m_FCALET
TH1F * m_FCALET
Definition: HIJetUEMonitoring.h:48
HIJetUEMonitoring::m_SubtractedET_Expected_eta_0_10
TProfile * m_SubtractedET_Expected_eta_0_10
Definition: HIJetUEMonitoring.h:73
HIJetUEMonitoring::m_refContainerName
std::string m_refContainerName
Definition: HIJetUEMonitoring.h:45
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
HIJetUEMonitoring::m_SubtractedET_Expected_eta_60_100
TProfile * m_SubtractedET_Expected_eta_60_100
Definition: HIJetUEMonitoring.h:131
HIJetUEMonitoring::m_2dSubtractedET_Expected_eta
TH2F * m_2dSubtractedET_Expected_eta
Definition: HIJetUEMonitoring.h:55
HIJetUEMonitoring::m_SubtractedET_eta_10_20
TProfile * m_SubtractedET_eta_10_20
Definition: HIJetUEMonitoring.h:97
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_phi_20_40
TH1F * m_JetUnsubtractedScaleMomentum_phi_20_40
Definition: HIJetUEMonitoring.h:104
HIJetUEMonitoring::buildHistos
virtual int buildHistos()
Definition: HIJetUEMonitoring.cxx:30
RCU::Shell
Definition: ShellExec.cxx:28
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
HIJetUEMonitoring::m_2dSubtractedET_2Dphi_0_10
TH2F * m_2dSubtractedET_2Dphi_0_10
Definition: HIJetUEMonitoring.h:75
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_m_60_100
TH1F * m_JetUnsubtractedScaleMomentum_m_60_100
Definition: HIJetUEMonitoring.h:124
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_eta_60_100
TH1F * m_JetUnsubtractedScaleMomentum_eta_60_100
Definition: HIJetUEMonitoring.h:122
HIJetUEMonitoring.h
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_pt_0_10
TH1F * m_JetSubtractedScaleMomentum_pt_0_10
Definition: HIJetUEMonitoring.h:67
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_eta_20_40
TH1F * m_JetUnsubtractedScaleMomentum_eta_20_40
Definition: HIJetUEMonitoring.h:103
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_phi_10_20
TH1F * m_JetUnsubtractedScaleMomentum_phi_10_20
Definition: HIJetUEMonitoring.h:85
HIJetUEMonitoring::m_SubtractedET_2Dphi_20_40
TProfile * m_SubtractedET_2Dphi_20_40
Definition: HIJetUEMonitoring.h:118
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_eta_10_20
TH1F * m_JetSubtractedScaleMomentum_eta_10_20
Definition: HIJetUEMonitoring.h:88
CHECK
#define CHECK(...)
Evaluate an expression and check for errors.
Definition: Control/AthenaKernel/AthenaKernel/errorcheck.h:422
HIJetUEMonitoring::m_ptcut
float m_ptcut
Definition: HIJetUEMonitoring.h:41
ToolHandleHistoHelper::HistoRetriever
Definition: ToolHandleHistoHelper.h:26
HIJetUEMonitoring::m_2dSubtractedET_2Dphi_20_40
TH2F * m_2dSubtractedET_2Dphi_20_40
Definition: HIJetUEMonitoring.h:114
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_m_20_40
TH1F * m_JetSubtractedScaleMomentum_m_20_40
Definition: HIJetUEMonitoring.h:109
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_m_10_20
TH1F * m_JetUnsubtractedScaleMomentum_m_10_20
Definition: HIJetUEMonitoring.h:86
HIJetUEMonitoring::m_SubtractedET_eta_0_10
TProfile * m_SubtractedET_eta_0_10
Definition: HIJetUEMonitoring.h:77
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_phi_60_100
TH1F * m_JetSubtractedScaleMomentum_phi_60_100
Definition: HIJetUEMonitoring.h:127
HIJetUEMonitoring::m_2dSubtractedET_Expected_eta_60_100
TH2F * m_2dSubtractedET_Expected_eta_60_100
Definition: HIJetUEMonitoring.h:130
HIJetUEMonitoring::m_harmonic
unsigned int m_harmonic
Definition: HIJetUEMonitoring.h:40
HIJetUEMonitoring::m_container_key
std::string m_container_key
Definition: HIJetUEMonitoring.h:32
HIJetUEMonitoring::m_2dSubtractedET_Expected_eta_20_40
TH2F * m_2dSubtractedET_Expected_eta_20_40
Definition: HIJetUEMonitoring.h:111
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:192
toTeV
#define toTeV
Definition: HIJetUEMonitoring.cxx:17
toGeV
#define toGeV
Definition: HIJetUEMonitoring.cxx:16
TProfile
Definition: rootspy.cxx:515
HIJetUEMonitoring::m_SubtractedET_Expected_eta_10_20
TProfile * m_SubtractedET_Expected_eta_10_20
Definition: HIJetUEMonitoring.h:93
HIJetUEMonitoring::m_2dSubtractedET_2Dphi
TH2F * m_2dSubtractedET_2Dphi
Definition: HIJetUEMonitoring.h:57
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_pt_60_100
TH1F * m_JetUnsubtractedScaleMomentum_pt_60_100
Definition: HIJetUEMonitoring.h:121
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_phi_0_10
TH1F * m_JetUnsubtractedScaleMomentum_phi_0_10
Definition: HIJetUEMonitoring.h:65
xAOD::JetAssignedScaleMomentum
@ JetAssignedScaleMomentum
Definition: JetTypes.h:30
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_m_60_100
TH1F * m_JetSubtractedScaleMomentum_m_60_100
Definition: HIJetUEMonitoring.h:128
HIJetUEMonitoring::m_SubtractedET_eta_20_40
TProfile * m_SubtractedET_eta_20_40
Definition: HIJetUEMonitoring.h:116
HIJetUEMonitoring::m_SubtractedET_pt_0_10
TProfile * m_SubtractedET_pt_0_10
Definition: HIJetUEMonitoring.h:78
EventInfo.h
xAOD::EventInfo_v1
Class describing the basic event information.
Definition: EventInfo_v1.h:43
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_eta_60_100
TH1F * m_JetSubtractedScaleMomentum_eta_60_100
Definition: HIJetUEMonitoring.h:126
ToolHandleHistoHelper::HistoRetriever::build
T * build(const std::string &n)
Definition: ToolHandleHistoHelper.h:44
xAOD::Jet_v1
Class describing a jet.
Definition: Jet_v1.h:57
HIJetUEMonitoring::m_SubtractedET_eta_60_100
TProfile * m_SubtractedET_eta_60_100
Definition: HIJetUEMonitoring.h:135
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_pt_60_100
TH1F * m_JetSubtractedScaleMomentum_pt_60_100
Definition: HIJetUEMonitoring.h:125
TProfile::Fill
int Fill(double, double)
Definition: rootspy.cxx:523
TH1F
Definition: rootspy.cxx:320
HIJetUEMonitoring::m_SubtractedET_2Dphi_0_10
TProfile * m_SubtractedET_2Dphi_0_10
Definition: HIJetUEMonitoring.h:80
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_phi_10_20
TH1F * m_JetSubtractedScaleMomentum_phi_10_20
Definition: HIJetUEMonitoring.h:89
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_m_0_10
TH1F * m_JetUnsubtractedScaleMomentum_m_0_10
Definition: HIJetUEMonitoring.h:66
HIJetUEMonitoring::m_eventShape
const xAOD::HIEventShapeContainer * m_eventShape
Definition: HIJetUEMonitoring.h:33
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_pt_20_40
TH1F * m_JetSubtractedScaleMomentum_pt_20_40
Definition: HIJetUEMonitoring.h:106
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_eta_0_10
TH1F * m_JetUnsubtractedScaleMomentum_eta_0_10
Definition: HIJetUEMonitoring.h:64
HIJetUEMonitoring::m_2dSubtractedET_Expected_Centrality
TH2F * m_2dSubtractedET_Expected_Centrality
Definition: HIJetUEMonitoring.h:52
xAOD::EventInfo_v1::errorState
EventFlagErrorState errorState(EventFlagSubDet subDet) const
Get the error state for a particular sub-detector.
Definition: EventInfo_v1.cxx:817
HistoGroupBase::bookHisto
T * bookHisto(T *h, Interval_t ityp=useToolInterval)
register the histo h in this group (if h!=NULL). The histo name is changed if m_prefixedHistoName==tr...
Definition: HistoGroupBase.h:69
HIJetUEMonitoring::m_SubtractedET_pt_60_100
TProfile * m_SubtractedET_pt_60_100
Definition: HIJetUEMonitoring.h:136
ConstAccessor.h
Helper class to provide constant type-safe access to aux data.
HIJetUEMonitoring::m_2dSubtractedET_Expected_eta_10_20
TH2F * m_2dSubtractedET_Expected_eta_10_20
Definition: HIJetUEMonitoring.h:92
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_eta_10_20
TH1F * m_JetUnsubtractedScaleMomentum_eta_10_20
Definition: HIJetUEMonitoring.h:84
HIJetUEMonitoring::HIJetUEMonitoring
HIJetUEMonitoring(const std::string &name)
Definition: HIJetUEMonitoring.cxx:19
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_phi_0_10
TH1F * m_JetSubtractedScaleMomentum_phi_0_10
Definition: HIJetUEMonitoring.h:69
HIJetUEMonitoring::m_SubtractedET_Expected_eta_20_40
TProfile * m_SubtractedET_Expected_eta_20_40
Definition: HIJetUEMonitoring.h:112
HIJetUEMonitoring::m_SubtractedET_2Dphi_60_100
TProfile * m_SubtractedET_2Dphi_60_100
Definition: HIJetUEMonitoring.h:137
readCCLHist.float
float
Definition: readCCLHist.py:83
HIJetUEMonitoring::m_2dSubtractedET_Centrality
TH2F * m_2dSubtractedET_Centrality
Definition: HIJetUEMonitoring.h:51
ToolHandleHistoHelper.h
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_phi_20_40
TH1F * m_JetSubtractedScaleMomentum_phi_20_40
Definition: HIJetUEMonitoring.h:108
SG::ConstAccessor::withDefault
const_reference_type withDefault(const ELT &e, const T &deflt) const
Fetch the variable for one element, as a const reference, with a default.
SCT_Monitoring::summary
@ summary
Definition: SCT_MonitoringNumbers.h:65
Generated on Sun Jun 30 2024 21:16:52 for ATLAS Offline Software by   doxygen 1.8.18