ATLAS Offline Software
Loading...
Searching...
No Matches
PhysicsAnalysis/TauID/TauDQA/src/EfficiencyPlots.cxx
Go to the documentation of this file.
1/*
2 Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
3*/
4
5#include "EfficiencyPlots.h"
6#include "AthenaKernel/Units.h"
7
8namespace Tau{
9
10 EfficiencyPlots::EfficiencyPlots(PlotBase* pParent, const std::string& sDir, std::string sTauJetContainerName):
11 PlotBase(pParent, sDir),
12 m_sTauJetContainerName(std::move(sTauJetContainerName))
13 {
14 }
15
16 void EfficiencyPlots::initializePlots()
17 {
18 m_eff_pt_jetGNTauloose = BookTProfile("Eff_Pt_jetGNTauloose"," Matched Tau loose GNTau eff in pt; pt; eff", 20, 0., 150.0);
19 m_eff_pt_jetGNTaumed = BookTProfile("Eff_Pt_jetGNTaumed","Matched Tau med GNTau eff in pt; pt; eff", 20, 0.0, 150.0);
20 m_eff_pt_jetGNTautight = BookTProfile("Eff_Pt_jetGNTautight","Matched Tau tight GNTau eff in pt; pt; eff", 20, 0.0, 150.0);
21 m_eff_pt_jetGNTaulooseHighPt = BookTProfile("Eff_Pt_jetGNTaulooseHightPt"," Matched Tau loose GNTau eff in pt; pt; eff", 20, 0.0, 1500.0);
22 m_eff_pt_jetGNTaumedHighPt = BookTProfile("Eff_Pt_jetGNTaumedHightPt","Matched Tau med GNTau eff in pt; pt; eff", 20, 0.0, 1500.0);
23 m_eff_pt_jetGNTautightHighPt = BookTProfile("Eff_Pt_jetGNTautightHightPt","Matched Tau tight GNTau eff in pt; pt; eff", 20, 0.0, 1500.0);
24 m_eff_mu_jetGNTauloose = BookTProfile("Eff_Mu_jetGNTauloose"," Matched Tau loose GNTau eff in <mu>; <mu>; eff", 16, 0., 80.0);
25 m_eff_mu_jetGNTaumed = BookTProfile("Eff_Mu_jetGNTaumed","Matched Tau med GNTau eff in <mu>; <mu>; eff", 16, 0.0, 80.0);
26 m_eff_mu_jetGNTautight = BookTProfile("Eff_Mu_jetGNTautight","Matched Tau tight GNTau eff in <mu>; <mu>; eff", 16, 0.0, 80.0);
27 m_eff_jetGNTauloose = BookTProfile("Eff_jetGNTauloose"," Matched Tau loose GNTau eff total; bin; eff",3,-1.5,1.5);
28 m_eff_jetGNTaumed = BookTProfile("Eff_jetGNTaumed","Matched Tau med GNTau eff total; bin; eff",3,-1.5,1.5);
29 m_eff_jetGNTautight = BookTProfile("Eff_jetGNTautight","Matched Tau tight GNTau eff total; bin; eff",3,-1.5,1.5);
30
31 m_eff_pt_eVetoloose = BookTProfile("Eff_Pt_eVetoloose"," loose eVeto eff in pt; pt; eff", 7, 0.0, 140.0);
32 m_eff_pt_eVetomed = BookTProfile("Eff_Pt_eVetomed"," med eVeto eff in pt; pt; eff", 7, 0.0, 140.0);
33 m_eff_pt_eVetotight = BookTProfile("Eff_Pt_eVetotight","tight eVeto eff in pt; pt; eff", 7, 0.0, 140.0);
34 m_eff_pt_eVetolooseHighPt = BookTProfile("Eff_Pt_eVetolooseHightPt"," loose eVeto eff in pt; pt; eff", 20, 0.0, 1500.0);
35 m_eff_pt_eVetomedHighPt = BookTProfile("Eff_Pt_eVetomedHightPt","med eVeto eff in pt; pt; eff", 20, 0.0, 1500.0);
36 m_eff_pt_eVetotightHighPt = BookTProfile("Eff_Pt_eVetotightHightPt","tight eVeto eff in pt; pt; eff", 20, 0.0, 1500.0);
37 m_eff_eta_eVetoloose = BookTProfile("Eff_Eta_eVetoloose"," loose eVeto eff in eta; eta; eff", 15, -3.0,3.0);
38 m_eff_eta_eVetomed = BookTProfile("Eff_Eta_eVetomed"," med eVeto eff in eta; eta; eff", 15, -3.0,3.0);
39 m_eff_eta_eVetotight = BookTProfile("Eff_Eta_eVetotight","tight eVeto eff in eta; eta; eff", 15, -3.0,3.0);
40 m_eff_eVetoloose = BookTProfile("Eff_eVetoloose"," loose eVeto eff total; bin; eff",3,-1.5,1.5);
41 m_eff_eVetomed = BookTProfile("Eff_eVetomed","med eVeto eff total; bin; eff",3,-1.5,1.5);
42 m_eff_eVetotight = BookTProfile("Eff_eVetotight","tight eVeto eff total; bin; eff",3,-1.5,1.5);
43
44 }
45
46 void EfficiencyPlots::fill(const xAOD::TauJet& tau, float weight, float avg_mu)
47 {
48 static const SG::ConstAccessor<char> acc_GNTauL("GNTauL_v0prune");
49 double pass_loose = acc_GNTauL.withDefault(tau,false);
50 m_eff_pt_jetGNTauloose->Fill(tau.pt()/Athena::Units::GeV, pass_loose, weight);
51 m_eff_pt_jetGNTaulooseHighPt->Fill(tau.pt()/Athena::Units::GeV, pass_loose, weight);
52 m_eff_jetGNTauloose->Fill(0.,pass_loose , weight);
53 m_eff_mu_jetGNTauloose->Fill(avg_mu,pass_loose, weight);
54
55 static const SG::ConstAccessor<char> acc_GNTauM("GNTauM_v0prune");
56 double pass_medium = acc_GNTauM.withDefault(tau,false);
57 m_eff_pt_jetGNTaumed->Fill(tau.pt()/Athena::Units::GeV, pass_medium, weight);
58 m_eff_pt_jetGNTaumedHighPt->Fill(tau.pt()/Athena::Units::GeV, pass_medium, weight);
59 m_eff_jetGNTaumed->Fill(0., pass_medium, weight);
60 m_eff_mu_jetGNTaumed->Fill(avg_mu,pass_medium, weight);
61
62 static const SG::ConstAccessor<char> acc_GNTauT("GNTauT_v0prune");
63 double pass_tight = acc_GNTauT.withDefault(tau,false);
64 m_eff_pt_jetGNTautight->Fill(tau.pt()/Athena::Units::GeV, pass_tight, weight);
65 m_eff_pt_jetGNTautightHighPt->Fill(tau.pt()/Athena::Units::GeV, pass_tight, weight);
66 m_eff_jetGNTautight->Fill(0., pass_tight, weight);
67 m_eff_mu_jetGNTautight->Fill(avg_mu,pass_tight, weight);
68
69 static const SG::ConstAccessor<char> acc_RNNEleLoose("EleRNNLoose_v1");
70 double pass_eVeto_loose = acc_RNNEleLoose.withDefault(tau,false);
71 m_eff_pt_eVetoloose->Fill(tau.pt()/Athena::Units::GeV, pass_eVeto_loose, weight);
72 m_eff_pt_eVetolooseHighPt->Fill(tau.pt()/Athena::Units::GeV, pass_eVeto_loose, weight);
73 m_eff_eVetoloose->Fill(0., pass_eVeto_loose, weight);
74 m_eff_eta_eVetoloose->Fill(tau.eta(), pass_eVeto_loose, weight);
75
76
77 static const SG::ConstAccessor<char> acc_RNNEleMedium("EleRNNMedium_v1");
78 double pass_eVeto_medium = acc_RNNEleMedium.withDefault(tau,false);
79 m_eff_pt_eVetomed->Fill(tau.pt()/Athena::Units::GeV, pass_eVeto_medium, weight);
80 m_eff_pt_eVetomedHighPt->Fill(tau.pt()/Athena::Units::GeV, pass_eVeto_medium, weight);
81 m_eff_eVetomed->Fill(0., pass_eVeto_medium, weight);
82 m_eff_eta_eVetomed->Fill(tau.eta(), pass_eVeto_medium, weight);
83
84 static const SG::ConstAccessor<char> acc_RNNEleTight("EleRNNTight_v1");
85 double pass_eVeto_tight = acc_RNNEleTight.withDefault(tau,false);
86 m_eff_pt_eVetotight->Fill(tau.pt()/Athena::Units::GeV, pass_eVeto_tight, weight);
87 m_eff_pt_eVetotightHighPt->Fill(tau.pt()/Athena::Units::GeV, pass_eVeto_tight, weight);
88 m_eff_eVetotight->Fill(0., pass_eVeto_tight, weight);
89 m_eff_eta_eVetotight->Fill(tau.eta(), pass_eVeto_tight, weight);
90
91 }
92
93
94}
Units.h
Wrapper to avoid constant divisions when using units.
PlotBase::BookTProfile
TProfile * BookTProfile(std::string_view name, std::string_view labels, int nBinsX, float startX, float endX, float startY=-1, float endY=-1, bool prependDir=true, bool useRMS=false)
Book a TProfile histogram.
Definition PlotBase.cxx:186
PlotBase::PlotBase
PlotBase(PlotBase *parent, std::string_view sDir)
Definition PlotBase.cxx:29
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.
Tau::EfficiencyPlots::fill
void fill(const xAOD::TauJet &tau, float weight, float avg_mu)
Definition PhysicsAnalysis/TauID/TauDQA/src/EfficiencyPlots.cxx:46
Tau::EfficiencyPlots::EfficiencyPlots
EfficiencyPlots(PlotBase *pParent, const std::string &sDir, std::string sTauJetContainerName)
Definition PhysicsAnalysis/TauID/TauDQA/src/EfficiencyPlots.cxx:10
xAOD::TauJet_v3::pt
virtual double pt() const
The transverse momentum ( ) of the particle.
xAOD::TauJet_v3::eta
virtual double eta() const
The pseudorapidity ( ) of the particle.
std
STL namespace.
xAOD::TauJet
TauJet_v3 TauJet
Definition of the current "tau version".
Definition Event/xAOD/xAODTau/xAODTau/TauJet.h:17
Generated on for ATLAS Offline Software by doxygen 1.16.1