ATLAS Offline Software
Loading...
Searching...
No Matches
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
19HIJetUEMonitoring::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
26StatusCode HIJetUEMonitoring::initialize() {
27 return StatusCode::SUCCESS;
28}
29
30int 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
138int 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
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition AthMsgStreamMacros.h:31
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition AthMsgStreamMacros.h:29
ConstAccessor.h
Helper class to provide constant type-safe access to aux data.
CHECK
#define CHECK(...)
Evaluate an expression and check for errors.
Definition Control/AthenaKernel/AthenaKernel/errorcheck.h:422
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition AthCommonDataStore.h:145
AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
Definition AthCommonDataStore.h:85
HIJetUEMonitoring::m_SubtractedET_2Dphi_60_100
TProfile * m_SubtractedET_2Dphi_60_100
Definition HIJetUEMonitoring.h:135
HIJetUEMonitoring::m_SubtractedET_2Dphi
TProfile * m_SubtractedET_2Dphi
Definition HIJetUEMonitoring.h:58
HIJetUEMonitoring::m_2dSubtractedET_Expected_Centrality
TH2F * m_2dSubtractedET_Expected_Centrality
Definition HIJetUEMonitoring.h:50
HIJetUEMonitoring::m_SubtractedET_eta_20_40
TProfile * m_SubtractedET_eta_20_40
Definition HIJetUEMonitoring.h:114
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_m_60_100
TH1F * m_JetSubtractedScaleMomentum_m_60_100
Definition HIJetUEMonitoring.h:126
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_phi_20_40
TH1F * m_JetSubtractedScaleMomentum_phi_20_40
Definition HIJetUEMonitoring.h:106
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_m_60_100
TH1F * m_JetUnsubtractedScaleMomentum_m_60_100
Definition HIJetUEMonitoring.h:122
HIJetUEMonitoring::HIJetUEMonitoring
HIJetUEMonitoring(const std::string &name)
Definition HIJetUEMonitoring.cxx:19
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_m_20_40
TH1F * m_JetUnsubtractedScaleMomentum_m_20_40
Definition HIJetUEMonitoring.h:103
HIJetUEMonitoring::m_SubtractedET_2Dphi_20_40
TProfile * m_SubtractedET_2Dphi_20_40
Definition HIJetUEMonitoring.h:116
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_phi_60_100
TH1F * m_JetUnsubtractedScaleMomentum_phi_60_100
Definition HIJetUEMonitoring.h:121
HIJetUEMonitoring::m_SubtractedET_Centrality
TProfile * m_SubtractedET_Centrality
Definition HIJetUEMonitoring.h:48
HIJetUEMonitoring::m_2dSubtractedET_2Dphi_10_20
TH2F * m_2dSubtractedET_2Dphi_10_20
Definition HIJetUEMonitoring.h:93
HIJetUEMonitoring::m_SubtractedET_pt_20_40
TProfile * m_SubtractedET_pt_20_40
Definition HIJetUEMonitoring.h:115
HIJetUEMonitoring::m_SubtractedET_Expected_eta_60_100
TProfile * m_SubtractedET_Expected_eta_60_100
Definition HIJetUEMonitoring.h:129
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_phi_60_100
TH1F * m_JetSubtractedScaleMomentum_phi_60_100
Definition HIJetUEMonitoring.h:125
HIJetUEMonitoring::m_2dSubtractedET_Expected_eta_20_40
TH2F * m_2dSubtractedET_Expected_eta_20_40
Definition HIJetUEMonitoring.h:109
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_m_20_40
TH1F * m_JetSubtractedScaleMomentum_m_20_40
Definition HIJetUEMonitoring.h:107
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_eta_0_10
TH1F * m_JetUnsubtractedScaleMomentum_eta_0_10
Definition HIJetUEMonitoring.h:62
HIJetUEMonitoring::m_histoDef
ToolHandleArray< HistoDefinitionTool > m_histoDef
Definition HIJetUEMonitoring.h:42
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_eta_20_40
TH1F * m_JetUnsubtractedScaleMomentum_eta_20_40
Definition HIJetUEMonitoring.h:101
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_pt_60_100
TH1F * m_JetUnsubtractedScaleMomentum_pt_60_100
Definition HIJetUEMonitoring.h:119
HIJetUEMonitoring::m_SubtractedET_Expected_eta_20_40
TProfile * m_SubtractedET_Expected_eta_20_40
Definition HIJetUEMonitoring.h:110
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_phi_0_10
TH1F * m_JetSubtractedScaleMomentum_phi_0_10
Definition HIJetUEMonitoring.h:67
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_eta_10_20
TH1F * m_JetSubtractedScaleMomentum_eta_10_20
Definition HIJetUEMonitoring.h:86
HIJetUEMonitoring::m_SubtractedET_eta
TProfile * m_SubtractedET_eta
Definition HIJetUEMonitoring.h:56
HIJetUEMonitoring::m_SubtractedET_eta_10_20
TProfile * m_SubtractedET_eta_10_20
Definition HIJetUEMonitoring.h:95
HIJetUEMonitoring::m_container_key
std::string m_container_key
Definition HIJetUEMonitoring.h:30
HIJetUEMonitoring::m_SubtractedET_Expected_eta
TProfile * m_SubtractedET_Expected_eta
Definition HIJetUEMonitoring.h:54
HIJetUEMonitoring::m_eventShape
const xAOD::HIEventShapeContainer * m_eventShape
Definition HIJetUEMonitoring.h:31
HIJetUEMonitoring::m_2dSubtractedET_Expected_eta_60_100
TH2F * m_2dSubtractedET_Expected_eta_60_100
Definition HIJetUEMonitoring.h:128
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_pt_0_10
TH1F * m_JetSubtractedScaleMomentum_pt_0_10
Definition HIJetUEMonitoring.h:65
HIJetUEMonitoring::m_FCalET
float m_FCalET
Centrality.
Definition HIJetUEMonitoring.h:34
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_phi_20_40
TH1F * m_JetUnsubtractedScaleMomentum_phi_20_40
Definition HIJetUEMonitoring.h:102
HIJetUEMonitoring::m_SubtractedET_pt
TProfile * m_SubtractedET_pt
Definition HIJetUEMonitoring.h:57
HIJetUEMonitoring::m_SubtractedET_eta_0_10
TProfile * m_SubtractedET_eta_0_10
Definition HIJetUEMonitoring.h:75
HIJetUEMonitoring::m_SubtractedET_pt_10_20
TProfile * m_SubtractedET_pt_10_20
Definition HIJetUEMonitoring.h:96
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_eta_60_100
TH1F * m_JetSubtractedScaleMomentum_eta_60_100
Definition HIJetUEMonitoring.h:124
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_phi_10_20
TH1F * m_JetUnsubtractedScaleMomentum_phi_10_20
Definition HIJetUEMonitoring.h:83
HIJetUEMonitoring::m_SubtractedET_Expected_eta_10_20
TProfile * m_SubtractedET_Expected_eta_10_20
Definition HIJetUEMonitoring.h:91
HIJetUEMonitoring::m_SubtractedET_2Dphi_10_20
TProfile * m_SubtractedET_2Dphi_10_20
Definition HIJetUEMonitoring.h:97
HIJetUEMonitoring::m_2dSubtractedET_Expected_eta_0_10
TH2F * m_2dSubtractedET_Expected_eta_0_10
Definition HIJetUEMonitoring.h:70
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_eta_60_100
TH1F * m_JetUnsubtractedScaleMomentum_eta_60_100
Definition HIJetUEMonitoring.h:120
HIJetUEMonitoring::m_SubtractedET_pt_60_100
TProfile * m_SubtractedET_pt_60_100
Definition HIJetUEMonitoring.h:134
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_eta_10_20
TH1F * m_JetUnsubtractedScaleMomentum_eta_10_20
Definition HIJetUEMonitoring.h:82
HIJetUEMonitoring::m_refContainerName
std::string m_refContainerName
Definition HIJetUEMonitoring.h:43
HIJetUEMonitoring::initialize
virtual StatusCode initialize()
Dummy implementation of the initialisation function.
Definition HIJetUEMonitoring.cxx:26
HIJetUEMonitoring::m_SubtractedET_eta_60_100
TProfile * m_SubtractedET_eta_60_100
Definition HIJetUEMonitoring.h:133
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_pt_0_10
TH1F * m_JetUnsubtractedScaleMomentum_pt_0_10
Definition HIJetUEMonitoring.h:61
HIJetUEMonitoring::buildHistos
virtual int buildHistos()
Definition HIJetUEMonitoring.cxx:30
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_eta_20_40
TH1F * m_JetSubtractedScaleMomentum_eta_20_40
Definition HIJetUEMonitoring.h:105
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_eta_0_10
TH1F * m_JetSubtractedScaleMomentum_eta_0_10
Definition HIJetUEMonitoring.h:66
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_pt_20_40
TH1F * m_JetUnsubtractedScaleMomentum_pt_20_40
Definition HIJetUEMonitoring.h:100
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_m_0_10
TH1F * m_JetUnsubtractedScaleMomentum_m_0_10
Definition HIJetUEMonitoring.h:64
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_m_0_10
TH1F * m_JetSubtractedScaleMomentum_m_0_10
Definition HIJetUEMonitoring.h:68
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_m_10_20
TH1F * m_JetUnsubtractedScaleMomentum_m_10_20
Definition HIJetUEMonitoring.h:84
HIJetUEMonitoring::m_2dSubtractedET_Expected_eta_10_20
TH2F * m_2dSubtractedET_Expected_eta_10_20
Definition HIJetUEMonitoring.h:90
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_pt_20_40
TH1F * m_JetSubtractedScaleMomentum_pt_20_40
Definition HIJetUEMonitoring.h:104
HIJetUEMonitoring::m_SubtractedET_2Dphi_0_10
TProfile * m_SubtractedET_2Dphi_0_10
Definition HIJetUEMonitoring.h:78
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_pt_10_20
TH1F * m_JetSubtractedScaleMomentum_pt_10_20
Definition HIJetUEMonitoring.h:85
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_pt_60_100
TH1F * m_JetSubtractedScaleMomentum_pt_60_100
Definition HIJetUEMonitoring.h:123
HIJetUEMonitoring::m_SubtractedET_Expected_eta_0_10
TProfile * m_SubtractedET_Expected_eta_0_10
Definition HIJetUEMonitoring.h:71
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_phi_10_20
TH1F * m_JetSubtractedScaleMomentum_phi_10_20
Definition HIJetUEMonitoring.h:87
HIJetUEMonitoring::m_2dSubtractedET_Expected_eta
TH2F * m_2dSubtractedET_Expected_eta
Definition HIJetUEMonitoring.h:53
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_pt_10_20
TH1F * m_JetUnsubtractedScaleMomentum_pt_10_20
Definition HIJetUEMonitoring.h:81
HIJetUEMonitoring::m_JetSubtractedScaleMomentum_m_10_20
TH1F * m_JetSubtractedScaleMomentum_m_10_20
Definition HIJetUEMonitoring.h:88
HIJetUEMonitoring::m_SubtractedET_pt_0_10
TProfile * m_SubtractedET_pt_0_10
Definition HIJetUEMonitoring.h:76
HIJetUEMonitoring::m_harmonic
unsigned int m_harmonic
Definition HIJetUEMonitoring.h:38
HIJetUEMonitoring::m_JetUnsubtractedScaleMomentum_phi_0_10
TH1F * m_JetUnsubtractedScaleMomentum_phi_0_10
Definition HIJetUEMonitoring.h:63
HIJetUEMonitoring::fillHistosFromJet
virtual int fillHistosFromJet(const xAOD::Jet &j, float weight)
Definition HIJetUEMonitoring.cxx:138
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
JetHistoBase::JetHistoBase
JetHistoBase(const std::string &t)
Definition JetHistoBase.cxx:10
SG::ConstAccessor
Helper class to provide constant type-safe access to aux data.
Definition ConstAccessor.h:55
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.
xAOD::EventInfo_v1::LAr
@ LAr
The LAr calorimeter.
Definition EventInfo_v1.h:335
xAOD::EventInfo_v1::Error
@ Error
The sub-detector issued an error.
Definition EventInfo_v1.h:349
xAOD::EventInfo_v1::errorState
EventFlagErrorState errorState(EventFlagSubDet subDet) const
Get the error state for a particular sub-detector.
Definition EventInfo_v1.cxx:817
xAOD::Jet_v1::getAttribute
bool getAttribute(AttributeID type, T &value) const
Retrieve attribute moment by enum.
xAOD
ICaloAffectedTool is abstract interface for tools checking if 4 mom is in calo affected region.
Definition ICaloAffectedTool.h:24
xAOD::Jet
Jet_v1 Jet
Definition of the current "jet version".
Definition Event/xAOD/xAODJet/xAODJet/Jet.h:17
xAOD::EventInfo
EventInfo_v1 EventInfo
Definition of the latest event info version.
Definition IEventInfoCnvTool.h:18
xAOD::HIEventShape
HIEventShape_v2 HIEventShape
Definition of the latest event info version.
Definition HIEventShape.h:16
Generated on for ATLAS Offline Software by doxygen 1.14.0