ATLAS Offline Software
Loading...
Searching...
No Matches
BTaggingValidationPlots.h
Go to the documentation of this file.
1/*
2 Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
3*/
4
5#ifndef JETTAGDQA_BTagPLOTS_H
6#define JETTAGDQA_BTagPLOTS_H
7
8#include "xAODBase/IParticle.h"
9#include "TrkValHistUtils/PlotBase.h"
10#include "AthenaBaseComps/AthMessaging.h"
11#include "xAODBTagging/BTagging.h" //typedef
12#include "xAODJet/JetFwd.h" //lightweight typedef
13#include "xAODTracking/VertexFwd.h"//lightweight typedef
14#include "xAODEventInfo/EventInfo.h" //typedef
15
16#include <string>
17#include <vector>
18#include <map>
19
20class TH1;
21
22
23namespace JetTagDQA{
24
25 class BTaggingValidationPlots : public PlotBase, public AthMessaging {
26 public:
27 // constructor
28 BTaggingValidationPlots(PlotBase* pParent, const std::string& sDir, std::string sParticleType);
29
30 // fill methods
31 void fillJetKinVars(const xAOD::Jet* jet, const int& truth_label, const bool& onZprime, const xAOD::EventInfo* event);
32 void fillDiscriminantVariables(const xAOD::Jet* jet, const double& jet_Lxy, const int& truth_label, const bool& onZprime, std::map<std::string, int>& nJetsThatPassedWPCuts, const xAOD::EventInfo* event);
33 void fillDiscriminantVariables_for_largeRjet(const xAOD::Jet* jet, const int& truth_label, const bool& onZprime, std::map<std::string, int>& nJetsThatPassedWPCuts, const xAOD::EventInfo* event);
34 void fillMultiplicities(const unsigned int& nJets, const unsigned int& nTracks, const int& nPrimVtx, const unsigned int& nTracksPrimVtx, const unsigned int& nJetsWithMuon, const unsigned int& nJetsWithSV, std::map<std::string, int>& nJetsThatPassedWPCuts, const xAOD::EventInfo* event);
35 void fillPVVariables(const double& PV_x, const double& PV_y, const double& PV_z, const xAOD::EventInfo* event);
36 void fillOther(const xAOD::Jet* jet, bool& contains_muon, double& jet_Lxy, const int& truth_label, const xAOD::EventInfo* event);
37 void fillTrackVariables(const xAOD::Jet* jet, const xAOD::Vertex *myVertex, const std::map<const xAOD::TrackParticle*, int> & track_truth_associations, const bool& contains_muon, const int& truth_label, int& num_HF_tracks_in_jet, const xAOD::EventInfo* event);
38 void fillTrackVariables_for_largeRjet(const xAOD::Jet* jet, const xAOD::Vertex *myVertex, const int& truth_label, const xAOD::EventInfo* event);
39 void fillSVVariables(const xAOD::Jet* jet, const std::map<const xAOD::TrackParticle*, int> & track_truth_associations, const bool& contains_muon, const int& truth_label, const int& num_HF_tracks_in_jet, bool& contains_SV, const xAOD::EventInfo* event);
40
41 void bookNJetsThatPassedWPCutsHistos();
42 void initializeNJetsThatPassedWPCutsMap(std::map<std::string, int>& nJetsThatPassedWPCuts);
43 void updateNJetsThatPassedWPCutsMap(std::map<std::string, int>& nJetsThatPassedWPCuts, const double& GN2v01, const double& GN3XPV01);
44 void fillNJetsThatPassedWPCutsHistos(std::map<std::string, int>& nJetsThatPassedWPCuts, const xAOD::EventInfo* event);
45
46 void setTaggerInfos();
47 void bookEffHistos();
48
49 // Reco only information
50 std::string m_sParticleType;
51
52 // multiplicities
53 TH1* m_nJets = nullptr;
54 TH1* m_nTracks = nullptr;
55 TH1* m_nPrimVtx = nullptr;
56 TH1* m_nTracksPrimVtx = nullptr;
57 TH1* m_nJetsWithMuon = nullptr;
58 TH1* m_nJetsWithSV = nullptr;
59 TH1* m_fracJetsWithMuon = nullptr;
60 TH1* m_fracJetsWithSV = nullptr;
61
62 // PV variables
63 TH1* m_PV_x = nullptr;
64 TH1* m_PV_y = nullptr;
65 TH1* m_PV_z = nullptr;
66
67 // jet kinematic variables
68 TH1* m_jet_e = nullptr;
69 TH1* m_jet_e_Zprime = nullptr;
70 TH1* m_jet_pt = nullptr;
71 TH1* m_jet_pt_Zprime = nullptr;
72 TH1* m_jet_eta = nullptr;
73 TH1* m_jet_phi = nullptr;
74
75 // muon vars
76 TH1* m_leading_muon_pT_frac = nullptr;
77 TH1* m_subleading_muon_pT_frac = nullptr;
78
79 // truth info
80 TH1* m_truthLabel = nullptr;
81
82 TH1* m_jet_pt_b = nullptr;
83 TH1* m_jet_pt_c = nullptr;
84 TH1* m_jet_pt_l = nullptr;
85 TH1* m_jet_pt_top = nullptr;
86 TH1* m_jet_pt_Zprime_b = nullptr;
87 TH1* m_jet_pt_Zprime_c = nullptr;
88 TH1* m_jet_pt_Zprime_l = nullptr;
89 TH1* m_jet_eta_b = nullptr;
90 TH1* m_jet_eta_c = nullptr;
91 TH1* m_jet_eta_l = nullptr;
92 TH1* m_jet_eta_top = nullptr;
93
94
95 // SV1 related vars
96 TH1* m_SV1_numSVs_incl = nullptr;
97 TH1* m_SV1_masssvx_incl = nullptr;
98 TH1* m_SV1_N2Tpair_incl = nullptr;
99 TH1* m_SV1_efracsvx_incl = nullptr;
100 TH1* m_SV1_deltaR_incl = nullptr;
101 TH1* m_SV1_significance3d_incl = nullptr;
102 TH1* m_SV1_energyTrkInJet_incl = nullptr;
103 TH1* m_SV1_NGTinSvx_incl = nullptr;
104 TH1* m_SV1_Lxy_incl = nullptr;
105 TH1* m_SV1_purity_incl = nullptr;
106 TH1* m_SV1_numSVs_b = nullptr;
107 TH1* m_SV1_masssvx_b = nullptr;
108 TH1* m_SV1_N2Tpair_b = nullptr;
109 TH1* m_SV1_efracsvx_b = nullptr;
110 TH1* m_SV1_deltaR_b = nullptr;
111 TH1* m_SV1_significance3d_b = nullptr;
112 TH1* m_SV1_energyTrkInJet_b = nullptr;
113 TH1* m_SV1_NGTinSvx_b = nullptr;
114 TH1* m_SV1_Lxy_b = nullptr;
115 TH1* m_SV1_purity_b = nullptr;
116 TH1* m_SV1_numSVs_c = nullptr;
117 TH1* m_SV1_masssvx_c = nullptr;
118 TH1* m_SV1_N2Tpair_c = nullptr;
119 TH1* m_SV1_efracsvx_c = nullptr;
120 TH1* m_SV1_deltaR_c = nullptr;
121 TH1* m_SV1_significance3d_c = nullptr;
122 TH1* m_SV1_energyTrkInJet_c = nullptr;
123 TH1* m_SV1_NGTinSvx_c = nullptr;
124 TH1* m_SV1_Lxy_c = nullptr;
125 TH1* m_SV1_purity_c = nullptr;
126 TH1* m_SV1_numSVs_l = nullptr;
127 TH1* m_SV1_masssvx_l = nullptr;
128 TH1* m_SV1_N2Tpair_l = nullptr;
129 TH1* m_SV1_efracsvx_l = nullptr;
130 TH1* m_SV1_deltaR_l = nullptr;
131 TH1* m_SV1_significance3d_l = nullptr;
132 TH1* m_SV1_energyTrkInJet_l = nullptr;
133 TH1* m_SV1_NGTinSvx_l = nullptr;
134 TH1* m_SV1_Lxy_l = nullptr;
135 TH1* m_SV1_purity_l = nullptr;
136 TH1* m_SV1_numSVs_muon = nullptr;
137 TH1* m_SV1_masssvx_muon = nullptr;
138 TH1* m_SV1_N2Tpair_muon = nullptr;
139 TH1* m_SV1_efracsvx_muon = nullptr;
140 TH1* m_SV1_deltaR_muon = nullptr;
141 TH1* m_SV1_significance3d_muon = nullptr;
142 TH1* m_SV1_energyTrkInJet_muon = nullptr;
143 TH1* m_SV1_NGTinSvx_muon = nullptr;
144 TH1* m_SV1_Lxy_muon = nullptr;
145 TH1* m_SV1_purity_muon = nullptr;
146
147 TH1* m_SV1_fracHFTracksInJet_incl = nullptr;
148 TH1* m_SV1_fracHFTracksInJet_b = nullptr;
149 TH1* m_SV1_fracHFTracksInJet_c = nullptr;
150 TH1* m_SV1_fracHFTracksInJet_l = nullptr;
151 TH1* m_SV1_fracHFTracksInJet_muon = nullptr;
152
153 TH1* m_SV1_fracTracks_fromB_incl = nullptr;
154 TH1* m_SV1_fracTracks_fromB_b = nullptr;
155 TH1* m_SV1_fracTracks_fromB_c = nullptr;
156 TH1* m_SV1_fracTracks_fromB_l = nullptr;
157 TH1* m_SV1_fracTracks_fromB_muon = nullptr;
158 TH1* m_SV1_fracTracks_fromC_incl = nullptr;
159 TH1* m_SV1_fracTracks_fromC_b = nullptr;
160 TH1* m_SV1_fracTracks_fromC_c = nullptr;
161 TH1* m_SV1_fracTracks_fromC_l = nullptr;
162 TH1* m_SV1_fracTracks_fromC_muon = nullptr;
163 TH1* m_SV1_fracTracks_fromFragmentation_incl = nullptr;
164 TH1* m_SV1_fracTracks_fromFragmentation_b = nullptr;
165 TH1* m_SV1_fracTracks_fromFragmentation_c = nullptr;
166 TH1* m_SV1_fracTracks_fromFragmentation_l = nullptr;
167 TH1* m_SV1_fracTracks_fromFragmentation_muon = nullptr;
168 TH1* m_SV1_fracTracks_fromSecondaries_incl = nullptr;
169 TH1* m_SV1_fracTracks_fromSecondaries_b = nullptr;
170 TH1* m_SV1_fracTracks_fromSecondaries_c = nullptr;
171 TH1* m_SV1_fracTracks_fromSecondaries_l = nullptr;
172 TH1* m_SV1_fracTracks_fromSecondaries_muon = nullptr;
173 TH1* m_SV1_fracTracks_fromPileup_incl = nullptr;
174 TH1* m_SV1_fracTracks_fromPileup_b = nullptr;
175 TH1* m_SV1_fracTracks_fromPileup_c = nullptr;
176 TH1* m_SV1_fracTracks_fromPileup_l = nullptr;
177 TH1* m_SV1_fracTracks_fromPileup_muon = nullptr;
178 TH1* m_SV1_fracTracks_fromFake_incl = nullptr;
179 TH1* m_SV1_fracTracks_fromFake_b = nullptr;
180 TH1* m_SV1_fracTracks_fromFake_c = nullptr;
181 TH1* m_SV1_fracTracks_fromFake_l = nullptr;
182 TH1* m_SV1_fracTracks_fromFake_muon = nullptr;
183
184 // these are only filled for detail level above 10
185 TH1* m_SV1_fracTracks_Secondaries_KshortDecay_incl = nullptr;
186 TH1* m_SV1_fracTracks_Secondaries_KshortDecay_b = nullptr;
187 TH1* m_SV1_fracTracks_Secondaries_KshortDecay_c = nullptr;
188 TH1* m_SV1_fracTracks_Secondaries_KshortDecay_u = nullptr;
189 TH1* m_SV1_fracTracks_Secondaries_KshortDecay_muon = nullptr;
190 TH1* m_SV1_fracTracks_Secondaries_LambdaDecay_incl = nullptr;
191 TH1* m_SV1_fracTracks_Secondaries_LambdaDecay_b = nullptr;
192 TH1* m_SV1_fracTracks_Secondaries_LambdaDecay_c = nullptr;
193 TH1* m_SV1_fracTracks_Secondaries_LambdaDecay_u = nullptr;
194 TH1* m_SV1_fracTracks_Secondaries_LambdaDecay_muon = nullptr;
195 TH1* m_SV1_fracTracks_Secondaries_GammaConversion_incl = nullptr;
196 TH1* m_SV1_fracTracks_Secondaries_GammaConversion_b = nullptr;
197 TH1* m_SV1_fracTracks_Secondaries_GammaConversion_c = nullptr;
198 TH1* m_SV1_fracTracks_Secondaries_GammaConversion_u = nullptr;
199 TH1* m_SV1_fracTracks_Secondaries_GammaConversion_muon = nullptr;
200 TH1* m_SV1_fracTracks_Secondaries_OtherDecay_incl = nullptr;
201 TH1* m_SV1_fracTracks_Secondaries_OtherDecay_b = nullptr;
202 TH1* m_SV1_fracTracks_Secondaries_OtherDecay_c = nullptr;
203 TH1* m_SV1_fracTracks_Secondaries_OtherDecay_u = nullptr;
204 TH1* m_SV1_fracTracks_Secondaries_OtherDecay_muon = nullptr;
205 TH1* m_SV1_fracTracks_Secondaries_HadronicInteraction_incl = nullptr;
206 TH1* m_SV1_fracTracks_Secondaries_HadronicInteraction_b = nullptr;
207 TH1* m_SV1_fracTracks_Secondaries_HadronicInteraction_c = nullptr;
208 TH1* m_SV1_fracTracks_Secondaries_HadronicInteraction_u = nullptr;
209 TH1* m_SV1_fracTracks_Secondaries_HadronicInteraction_muon = nullptr;
210 TH1* m_SV1_fracTracks_Secondaries_OtherSecondary_incl = nullptr;
211 TH1* m_SV1_fracTracks_Secondaries_OtherSecondary_b = nullptr;
212 TH1* m_SV1_fracTracks_Secondaries_OtherSecondary_c = nullptr;
213 TH1* m_SV1_fracTracks_Secondaries_OtherSecondary_u = nullptr;
214 TH1* m_SV1_fracTracks_Secondaries_OtherSecondary_muon = nullptr;
215 TH1* m_SV1_fracTracks_OtherOrigin_incl = nullptr;
216 TH1* m_SV1_fracTracks_OtherOrigin_b = nullptr;
217 TH1* m_SV1_fracTracks_OtherOrigin_c = nullptr;
218 TH1* m_SV1_fracTracks_OtherOrigin_u = nullptr;
219 TH1* m_SV1_fracTracks_OtherOrigin_muon = nullptr;
220
221 // IPs and IP significances
222 TH1* m_track_d0_incl = nullptr;
223 TH1* m_track_z0_incl = nullptr;
224 TH1* m_track_sigd0_incl = nullptr;
225 TH1* m_track_sigz0_incl = nullptr;
226
227 TH1* m_track_d0_b = nullptr;
228 TH1* m_track_z0_b = nullptr;
229 TH1* m_track_sigd0_b = nullptr;
230 TH1* m_track_sigz0_b = nullptr;
231
232 TH1* m_track_d0_c = nullptr;
233 TH1* m_track_z0_c = nullptr;
234 TH1* m_track_sigd0_c = nullptr;
235 TH1* m_track_sigz0_c = nullptr;
236
237 TH1* m_track_d0_u = nullptr;
238 TH1* m_track_z0_u = nullptr;
239 TH1* m_track_sigd0_u = nullptr;
240 TH1* m_track_sigz0_u = nullptr;
241
242 TH1* m_track_d0_muon = nullptr;
243 TH1* m_track_z0_muon = nullptr;
244 TH1* m_track_sigd0_muon = nullptr;
245 TH1* m_track_sigz0_muon = nullptr;
246
247 // pT_frac
248 TH1* m_track_pT_frac_incl = nullptr;
249 TH1* m_track_pT_frac_b = nullptr;
250 TH1* m_track_pT_frac_c = nullptr;
251 TH1* m_track_pT_frac_u = nullptr;
252 TH1* m_track_pT_frac_muon = nullptr;
253
254 // DeltaR_jet_track
255 TH1* m_DeltaR_jet_track_incl = nullptr;
256 TH1* m_DeltaR_jet_track_b = nullptr;
257 TH1* m_DeltaR_jet_track_c = nullptr;
258 TH1* m_DeltaR_jet_track_u = nullptr;
259 TH1* m_DeltaR_jet_track_muon = nullptr;
260
261 // tracker hits
262 TH1* m_nInnHits_incl = nullptr;
263 TH1* m_nNextToInnHits_incl = nullptr;
264 TH1* m_nBLHits_incl = nullptr;
265 TH1* m_nsharedBLHits_incl = nullptr;
266 TH1* m_nsplitBLHits_incl = nullptr;
267 TH1* m_nPixHits_incl = nullptr;
268 TH1* m_nPixHoles_incl = nullptr;
269 TH1* m_nsharedPixHits_incl = nullptr;
270 TH1* m_nsplitPixHits_incl = nullptr;
271 TH1* m_nSCTHits_incl = nullptr;
272 TH1* m_nSCTHoles_incl = nullptr;
273 TH1* m_nsharedSCTHits_incl = nullptr;
274
275 TH1* m_nInnHits_b = nullptr;
276 TH1* m_nNextToInnHits_b = nullptr;
277 TH1* m_nBLHits_b = nullptr;
278 TH1* m_nsharedBLHits_b = nullptr;
279 TH1* m_nsplitBLHits_b = nullptr;
280 TH1* m_nPixHits_b = nullptr;
281 TH1* m_nPixHoles_b = nullptr;
282 TH1* m_nsharedPixHits_b = nullptr;
283 TH1* m_nsplitPixHits_b = nullptr;
284 TH1* m_nSCTHits_b = nullptr;
285 TH1* m_nSCTHoles_b = nullptr;
286 TH1* m_nsharedSCTHits_b = nullptr;
287
288 TH1* m_nInnHits_c = nullptr;
289 TH1* m_nNextToInnHits_c = nullptr;
290 TH1* m_nBLHits_c = nullptr;
291 TH1* m_nsharedBLHits_c = nullptr;
292 TH1* m_nsplitBLHits_c = nullptr;
293 TH1* m_nPixHits_c = nullptr;
294 TH1* m_nPixHoles_c = nullptr;
295 TH1* m_nsharedPixHits_c = nullptr;
296 TH1* m_nsplitPixHits_c = nullptr;
297 TH1* m_nSCTHits_c = nullptr;
298 TH1* m_nSCTHoles_c = nullptr;
299 TH1* m_nsharedSCTHits_c = nullptr;
300
301 TH1* m_nInnHits_u = nullptr;
302 TH1* m_nNextToInnHits_u = nullptr;
303 TH1* m_nBLHits_u = nullptr;
304 TH1* m_nsharedBLHits_u = nullptr;
305 TH1* m_nsplitBLHits_u = nullptr;
306 TH1* m_nPixHits_u = nullptr;
307 TH1* m_nPixHoles_u = nullptr;
308 TH1* m_nsharedPixHits_u = nullptr;
309 TH1* m_nsplitPixHits_u = nullptr;
310 TH1* m_nSCTHits_u = nullptr;
311 TH1* m_nSCTHoles_u = nullptr;
312 TH1* m_nsharedSCTHits_u = nullptr;
313
314 TH1* m_nInnHits_muon = nullptr;
315 TH1* m_nNextToInnHits_muon = nullptr;
316 TH1* m_nBLHits_muon = nullptr;
317 TH1* m_nsharedBLHits_muon = nullptr;
318 TH1* m_nsplitBLHits_muon = nullptr;
319 TH1* m_nPixHits_muon = nullptr;
320 TH1* m_nPixHoles_muon = nullptr;
321 TH1* m_nsharedPixHits_muon = nullptr;
322 TH1* m_nsplitPixHits_muon = nullptr;
323 TH1* m_nSCTHits_muon = nullptr;
324 TH1* m_nSCTHoles_muon = nullptr;
325 TH1* m_nsharedSCTHits_muon = nullptr;
326
327 // numTracks_perJet
328 TH1* m_numTracks_perJet_incl = nullptr;
329 TH1* m_numTracks_perJet_b = nullptr;
330 TH1* m_numTracks_perJet_c = nullptr;
331 TH1* m_numTracks_perJet_u = nullptr;
332 TH1* m_numTracks_perJet_top = nullptr;
333 TH1* m_numTracks_perJet_muon = nullptr;
334
335 // number of tracks variables
336 TH1* m_numTracks_B_incl = nullptr;
337 TH1* m_numTracks_C_incl = nullptr;
338 TH1* m_numTracks_Fragmentation_incl = nullptr;
339 TH1* m_numTracks_Secondaries_incl = nullptr;
340 TH1* m_numTracks_Pileup_incl = nullptr;
341 TH1* m_numTracks_Fake_incl = nullptr;
342
343 TH1* m_numTracks_B_b = nullptr;
344 TH1* m_numTracks_C_b = nullptr;
345 TH1* m_numTracks_Fragmentation_b = nullptr;
346 TH1* m_numTracks_Secondaries_b = nullptr;
347 TH1* m_numTracks_Pileup_b = nullptr;
348 TH1* m_numTracks_Fake_b = nullptr;
349
350 TH1* m_numTracks_B_c = nullptr;
351 TH1* m_numTracks_C_c = nullptr;
352 TH1* m_numTracks_Fragmentation_c = nullptr;
353 TH1* m_numTracks_Secondaries_c = nullptr;
354 TH1* m_numTracks_Pileup_c = nullptr;
355 TH1* m_numTracks_Fake_c = nullptr;
356
357 TH1* m_numTracks_B_u = nullptr;
358 TH1* m_numTracks_C_u = nullptr;
359 TH1* m_numTracks_Fragmentation_u = nullptr;
360 TH1* m_numTracks_Secondaries_u = nullptr;
361 TH1* m_numTracks_Pileup_u = nullptr;
362 TH1* m_numTracks_Fake_u = nullptr;
363
364 TH1* m_numTracks_B_muon = nullptr;
365 TH1* m_numTracks_C_muon = nullptr;
366 TH1* m_numTracks_Fragmentation_muon = nullptr;
367 TH1* m_numTracks_Secondaries_muon = nullptr;
368 TH1* m_numTracks_Pileup_muon = nullptr;
369 TH1* m_numTracks_Fake_muon = nullptr;
370
371 // these are only filled for detail level above 10
372 TH1* m_numTracks_Secondaries_KshortDecay_incl = nullptr;
373 TH1* m_numTracks_Secondaries_KshortDecay_b = nullptr;
374 TH1* m_numTracks_Secondaries_KshortDecay_c = nullptr;
375 TH1* m_numTracks_Secondaries_KshortDecay_u = nullptr;
376 TH1* m_numTracks_Secondaries_KshortDecay_muon = nullptr;
377
378 TH1* m_numTracks_Secondaries_LambdaDecay_incl = nullptr;
379 TH1* m_numTracks_Secondaries_LambdaDecay_b = nullptr;
380 TH1* m_numTracks_Secondaries_LambdaDecay_c = nullptr;
381 TH1* m_numTracks_Secondaries_LambdaDecay_u = nullptr;
382 TH1* m_numTracks_Secondaries_LambdaDecay_muon = nullptr;
383
384 TH1* m_numTracks_Secondaries_GammaConversion_incl = nullptr;
385 TH1* m_numTracks_Secondaries_GammaConversion_b = nullptr;
386 TH1* m_numTracks_Secondaries_GammaConversion_c = nullptr;
387 TH1* m_numTracks_Secondaries_GammaConversion_u = nullptr;
388 TH1* m_numTracks_Secondaries_GammaConversion_muon = nullptr;
389
390 TH1* m_numTracks_Secondaries_OtherDecay_incl = nullptr;
391 TH1* m_numTracks_Secondaries_OtherDecay_b = nullptr;
392 TH1* m_numTracks_Secondaries_OtherDecay_c = nullptr;
393 TH1* m_numTracks_Secondaries_OtherDecay_u = nullptr;
394 TH1* m_numTracks_Secondaries_OtherDecay_muon = nullptr;
395
396 TH1* m_numTracks_Secondaries_HadronicInteraction_incl = nullptr;
397 TH1* m_numTracks_Secondaries_HadronicInteraction_b = nullptr;
398 TH1* m_numTracks_Secondaries_HadronicInteraction_c = nullptr;
399 TH1* m_numTracks_Secondaries_HadronicInteraction_u = nullptr;
400 TH1* m_numTracks_Secondaries_HadronicInteraction_muon = nullptr;
401
402 TH1* m_numTracks_Secondaries_OtherSecondary_incl = nullptr;
403 TH1* m_numTracks_Secondaries_OtherSecondary_b = nullptr;
404 TH1* m_numTracks_Secondaries_OtherSecondary_c = nullptr;
405 TH1* m_numTracks_Secondaries_OtherSecondary_u = nullptr;
406 TH1* m_numTracks_Secondaries_OtherSecondary_muon = nullptr;
407
408 TH1* m_numTracks_OtherOrigin_incl = nullptr;
409 TH1* m_numTracks_OtherOrigin_b = nullptr;
410 TH1* m_numTracks_OtherOrigin_c = nullptr;
411 TH1* m_numTracks_OtherOrigin_u = nullptr;
412 TH1* m_numTracks_OtherOrigin_muon = nullptr;
413
414 // features for largeRjet tagger
415 TH1* m_track_d0_top = nullptr;
416 TH1* m_track_z0_top = nullptr;
417 TH1* m_track_sigd0_top = nullptr;
418 TH1* m_track_sigz0_top = nullptr;
419 // pT_frac
420 TH1* m_track_pT_frac_top = nullptr;
421 // DeltaR_jet_track
422 TH1* m_DeltaR_jet_track_top = nullptr;
423 // tracker hits
424 TH1* m_nInnHits_top = nullptr;
425 TH1* m_nNextToInnHits_top = nullptr;
426 TH1* m_nBLHits_top = nullptr;
427 TH1* m_nsharedBLHits_top = nullptr;
428 TH1* m_nsplitBLHits_top = nullptr;
429 TH1* m_nPixHits_top = nullptr;
430 TH1* m_nPixHoles_top = nullptr;
431 TH1* m_nsharedPixHits_top = nullptr;
432 TH1* m_nsplitPixHits_top = nullptr;
433 TH1* m_nSCTHits_top = nullptr;
434 TH1* m_nSCTHoles_top = nullptr;
435 TH1* m_nsharedSCTHits_top = nullptr;
436
437 // tagger
438 TH1* m_GN2v01_pb = nullptr;
439 TH1* m_GN2v01_pc = nullptr;
440 TH1* m_GN2v01_pu = nullptr;
441 TH1* m_GN2v01_ptau = nullptr;
442
443 TH1* m_GN3XPV01_phtautauhad = nullptr;
444 TH1* m_GN3XPV01_phbb = nullptr;
445 TH1* m_GN3XPV01_phcc = nullptr;
446 TH1* m_GN3XPV01_ptop = nullptr;
447 TH1* m_GN3XPV01_pqcdbb = nullptr;
448 TH1* m_GN3XPV01_pqcdbx = nullptr;
449 TH1* m_GN3XPV01_pqcdcx = nullptr;
450 TH1* m_GN3XPV01_pqcdll = nullptr;
451 TH1* m_GN3XPV01_pwqq = nullptr;
452
453
454 // B hadron Lxy
455 TH1* m_Truth_Lxy_b = nullptr;
456 TH1* m_Truth_Lxy_c = nullptr;
457
458 // B hadron deltaR wrt jet
459 TH1* m_deltaR_truthBHadron_jet_b = nullptr;
460 TH1* m_deltaR_truthCHadron_jet_c = nullptr;
461
462 std::vector<std::string> m_taggers;
463 std::map<std::string, int> m_truthLabels;
464 std::map<std::string, double> m_GN2v01_workingPoints;
465 std::map<std::string, double> m_GN3XPV01_workingPoints; // TODO: Change this in the future since GN3 has WPs pT and mass dependent
466
467 double m_GN2v01_fc = 0.0;
468 double m_GN2v01_ftau = 0.0;
469 double m_GN3XPV01_hcc_fc = 0.0;
470 double m_GN3XPV01_top_fc = 0.0;
471 std::map<std::string, TH1*> m_weight_histos;
472
473 std::map<std::string, TH1*> m_nJetsThatPassedWPCutsHistos;
474
475 // detail level
476 void setDetailLevel(const unsigned int& detailLevel);
477
478 // a setter for the HistogramDefinitions and the jvt and TMP cuts
479 void setHistogramDefinitions( std::map< std::string, std::vector< std::string > > HistogramDefinitions);
480 void setIsDataJVTCutsAndTMPCut(bool isData, float JVTCutAntiKt4EMTopoJets, float JVTCutLargerEtaAntiKt4EMTopoJets, float JVTCutAntiKt4EMPFlowJets, float truthMatchProbabilityCut);
481 void setTaggerNames(const std::string& GN2v01Name, const std::string& GN3XPV01Name);
482
483 // jvt variables
484 bool m_JVT_defined{};
485 float m_JVT_cut = 0.0F;
486 bool m_JVTLargerEta_defined;
487 float m_JVTLargerEta_cut = 0.0F;
488
489 private:
490 virtual void initializePlots();
491 virtual void finalizePlots();
492
493 // detail level
494 unsigned int m_detailLevel = 0U;
495
496 // map with histogram definitions and the corresponding enum
497 std::map< std::string, std::vector< std::string > > m_HistogramDefinitions;
498 enum position{histo_name, histo_title, histo_path, histo_xbins, histo_xmin, histo_xmax, histo_type, histo_ymin, histo_ymax};
499 float m_truthMatchProbabilityCut = 0.0F;
500 bool m_isData = false;
501 // some helper functions
502 TH1* bookHistogram(std::string histo_name, const std::string& var_name, const std::string& part = "", const std::string& prefix = "");
503 int getTrackHits(const xAOD::TrackParticle& part, xAOD::SummaryType info);
504 void fillDiscriminantHistograms(const std::string& tagger_name, const double& discriminant_value, const std::map<std::string, double>& working_points, const int& truth_label, std::map<std::string, TH1*>::const_iterator hist_iter, std::map<std::string, int>::const_iterator label_iter, const bool& pass_nTracksCut, const double& jet_pT, const double& jet_Lxy, const bool& onZprime, const xAOD::EventInfo* event);
505 void bookDiscriminantVsPTAndLxyHistograms(const std::string& tagger_name, const std::map<std::string, double>& workingPoints, const bool& isOldTagger, std::map<std::string, int>::const_iterator label_iter, const std::string& m_sParticleType);
506 template <class T>
507 void fillHistoWithTruthCases(T value, TH1* histo_incl, TH1* histo_b, TH1* histo_c, TH1* histo_l, TH1* histo_muon, const int& truth_label, const bool& has_muon, const xAOD::EventInfo* event);
508 template <class T>
509 void fillHistoWithTruthCases_for_largeRjet(T value, TH1* histo_incl, TH1* histo_bb, TH1* histo_cc, TH1* histo_uu, TH1* histo_top, const int& truth_label, const xAOD::EventInfo* event);
510
511 // tagger names
512 std::string m_GN2v01Name;
513 std::string m_GN3XPV01Name;
514
515 };
516
517}
518
519#endif
AthMessaging::AthMessaging
AthMessaging(IMessageSvc *msgSvc, const std::string &name)
Constructor.
Definition AthMessaging.cxx:13
JetTagDQA::BTaggingValidationPlots::fillTrackVariables_for_largeRjet
void fillTrackVariables_for_largeRjet(const xAOD::Jet *jet, const xAOD::Vertex *myVertex, const int &truth_label, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:945
JetTagDQA::BTaggingValidationPlots::updateNJetsThatPassedWPCutsMap
void updateNJetsThatPassedWPCutsMap(std::map< std::string, int > &nJetsThatPassedWPCuts, const double &GN2v01, const double &GN3XPV01)
Definition BTaggingValidationPlots.cxx:1341
JetTagDQA::BTaggingValidationPlots::bookDiscriminantVsPTAndLxyHistograms
void bookDiscriminantVsPTAndLxyHistograms(const std::string &tagger_name, const std::map< std::string, double > &workingPoints, const bool &isOldTagger, std::map< std::string, int >::const_iterator label_iter, const std::string &m_sParticleType)
Definition BTaggingValidationPlots.cxx:152
JetTagDQA::BTaggingValidationPlots::setDetailLevel
void setDetailLevel(const unsigned int &detailLevel)
Definition BTaggingValidationPlots.cxx:35
JetTagDQA::BTaggingValidationPlots::fillTrackVariables
void fillTrackVariables(const xAOD::Jet *jet, const xAOD::Vertex *myVertex, const std::map< const xAOD::TrackParticle *, int > &track_truth_associations, const bool &contains_muon, const int &truth_label, int &num_HF_tracks_in_jet, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:798
JetTagDQA::BTaggingValidationPlots::fillHistoWithTruthCases
void fillHistoWithTruthCases(T value, TH1 *histo_incl, TH1 *histo_b, TH1 *histo_c, TH1 *histo_l, TH1 *histo_muon, const int &truth_label, const bool &has_muon, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:180
JetTagDQA::BTaggingValidationPlots::m_nJetsThatPassedWPCutsHistos
std::map< std::string, TH1 * > m_nJetsThatPassedWPCutsHistos
Definition BTaggingValidationPlots.h:473
JetTagDQA::BTaggingValidationPlots::fillDiscriminantVariables_for_largeRjet
void fillDiscriminantVariables_for_largeRjet(const xAOD::Jet *jet, const int &truth_label, const bool &onZprime, std::map< std::string, int > &nJetsThatPassedWPCuts, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:1245
JetTagDQA::BTaggingValidationPlots::setTaggerNames
void setTaggerNames(const std::string &GN2v01Name, const std::string &GN3XPV01Name)
Definition BTaggingValidationPlots.cxx:60
JetTagDQA::BTaggingValidationPlots::fillNJetsThatPassedWPCutsHistos
void fillNJetsThatPassedWPCutsHistos(std::map< std::string, int > &nJetsThatPassedWPCuts, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:1360
JetTagDQA::BTaggingValidationPlots::bookHistogram
TH1 * bookHistogram(std::string histo_name, const std::string &var_name, const std::string &part="", const std::string &prefix="")
Definition BTaggingValidationPlots.cxx:67
JetTagDQA::BTaggingValidationPlots::BTaggingValidationPlots
BTaggingValidationPlots(PlotBase *pParent, const std::string &sDir, std::string sParticleType)
Definition BTaggingValidationPlots.cxx:23
JetTagDQA::BTaggingValidationPlots::m_GN3XPV01_workingPoints
std::map< std::string, double > m_GN3XPV01_workingPoints
Definition BTaggingValidationPlots.h:465
JetTagDQA::BTaggingValidationPlots::fillHistoWithTruthCases_for_largeRjet
void fillHistoWithTruthCases_for_largeRjet(T value, TH1 *histo_incl, TH1 *histo_bb, TH1 *histo_cc, TH1 *histo_uu, TH1 *histo_top, const int &truth_label, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:203
JetTagDQA::BTaggingValidationPlots::fillPVVariables
void fillPVVariables(const double &PV_x, const double &PV_y, const double &PV_z, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:671
JetTagDQA::BTaggingValidationPlots::m_weight_histos
std::map< std::string, TH1 * > m_weight_histos
Definition BTaggingValidationPlots.h:471
JetTagDQA::BTaggingValidationPlots::setIsDataJVTCutsAndTMPCut
void setIsDataJVTCutsAndTMPCut(bool isData, float JVTCutAntiKt4EMTopoJets, float JVTCutLargerEtaAntiKt4EMTopoJets, float JVTCutAntiKt4EMPFlowJets, float truthMatchProbabilityCut)
Definition BTaggingValidationPlots.cxx:45
JetTagDQA::BTaggingValidationPlots::fillSVVariables
void fillSVVariables(const xAOD::Jet *jet, const std::map< const xAOD::TrackParticle *, int > &track_truth_associations, const bool &contains_muon, const int &truth_label, const int &num_HF_tracks_in_jet, bool &contains_SV, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:1032
JetTagDQA::BTaggingValidationPlots::m_truthLabels
std::map< std::string, int > m_truthLabels
Definition BTaggingValidationPlots.h:463
JetTagDQA::BTaggingValidationPlots::fillDiscriminantVariables
void fillDiscriminantVariables(const xAOD::Jet *jet, const double &jet_Lxy, const int &truth_label, const bool &onZprime, std::map< std::string, int > &nJetsThatPassedWPCuts, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:1207
JetTagDQA::BTaggingValidationPlots::setHistogramDefinitions
void setHistogramDefinitions(std::map< std::string, std::vector< std::string > > HistogramDefinitions)
Definition BTaggingValidationPlots.cxx:40
JetTagDQA::BTaggingValidationPlots::initializeNJetsThatPassedWPCutsMap
void initializeNJetsThatPassedWPCutsMap(std::map< std::string, int > &nJetsThatPassedWPCuts)
Definition BTaggingValidationPlots.cxx:1325
JetTagDQA::BTaggingValidationPlots::m_HistogramDefinitions
std::map< std::string, std::vector< std::string > > m_HistogramDefinitions
Definition BTaggingValidationPlots.h:497
JetTagDQA::BTaggingValidationPlots::fillOther
void fillOther(const xAOD::Jet *jet, bool &contains_muon, double &jet_Lxy, const int &truth_label, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:698
JetTagDQA::BTaggingValidationPlots::m_GN2v01_workingPoints
std::map< std::string, double > m_GN2v01_workingPoints
Definition BTaggingValidationPlots.h:464
JetTagDQA::BTaggingValidationPlots::fillJetKinVars
void fillJetKinVars(const xAOD::Jet *jet, const int &truth_label, const bool &onZprime, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:678
JetTagDQA::BTaggingValidationPlots::fillMultiplicities
void fillMultiplicities(const unsigned int &nJets, const unsigned int &nTracks, const int &nPrimVtx, const unsigned int &nTracksPrimVtx, const unsigned int &nJetsWithMuon, const unsigned int &nJetsWithSV, std::map< std::string, int > &nJetsThatPassedWPCuts, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:654
JetTagDQA::BTaggingValidationPlots::fillDiscriminantHistograms
void fillDiscriminantHistograms(const std::string &tagger_name, const double &discriminant_value, const std::map< std::string, double > &working_points, const int &truth_label, std::map< std::string, TH1 * >::const_iterator hist_iter, std::map< std::string, int >::const_iterator label_iter, const bool &pass_nTracksCut, const double &jet_pT, const double &jet_Lxy, const bool &onZprime, const xAOD::EventInfo *event)
Definition BTaggingValidationPlots.cxx:101
JetTagDQA::BTaggingValidationPlots::getTrackHits
int getTrackHits(const xAOD::TrackParticle &part, xAOD::SummaryType info)
Definition BTaggingValidationPlots.cxx:93
PlotBase::PlotBase
PlotBase(PlotBase *parent, const std::string &sDir)
Definition PlotBase.cxx:29
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::TrackParticle
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Definition Event/xAOD/xAODTracking/xAODTracking/TrackParticle.h:13
xAOD::Vertex
Vertex_v1 Vertex
Define the latest version of the vertex class.
Definition Event/xAOD/xAODTracking/xAODTracking/Vertex.h:16
xAOD::SummaryType
SummaryType
Enumerates the different types of information stored in Summary.
Definition TrackingPrimitives.h:229
Generated on for ATLAS Offline Software by doxygen 1.14.0