ATLAS Offline Software
Loading...
Searching...
No Matches
Public Member Functions | Public Attributes | Protected Attributes | Private Attributes | List of all members
egammaMonitoring::TruthElectronHistograms Class Reference

#include <TruthElectronHistograms.h>

Inheritance diagram for egammaMonitoring::TruthElectronHistograms:
Collaboration diagram for egammaMonitoring::TruthElectronHistograms:

Public Member Functions

virtual StatusCode initializePlots () override
StatusCode initializePlots (bool reducedHistSet)
void fill (const xAOD::TruthParticle *truth, const xAOD::Electron *el=nullptr)
virtual ~TruthElectronHistograms ()
void fill (const xAOD::IParticle &egamma)
void fill (const xAOD::IParticle &egamma, float mu)
void hasFwd (bool a=true)
 IHistograms ()
 IHistograms (std::string name, std::string title, std::string folder, SmartIF< ITHistSvc > rootHistSvc)

Public Attributes

std::map< std::string, TH2D * > histoMap2D
std::map< std::string, TH1D * > histoMap

Protected Attributes

std::string m_name
std::string m_title
std::string m_folder
SmartIF< ITHistSvc > m_rootHistSvc

Private Attributes

bool m_reducedHistSet {}
bool m_hasFwd = false

Detailed Description

Definition at line 18 of file TruthElectronHistograms.h.

Constructor & Destructor Documentation

◆ ~TruthElectronHistograms()

virtual egammaMonitoring::TruthElectronHistograms::~TruthElectronHistograms ( )
inlinevirtual

Definition at line 31 of file TruthElectronHistograms.h.

31{};

Member Function Documentation

◆ fill() [1/3]

void egammaMonitoring::ParticleHistograms::fill ( const xAOD::IParticle & egamma)
virtual

Reimplemented from egammaMonitoring::ParticleHistograms.

Definition at line 23 of file ParticleHistograms.cxx.

41 {
42 ParticleHistograms::fill(egamma,0.);
43 }
egammaMonitoring::ParticleHistograms::fill
void fill(const xAOD::IParticle &egamma)
Definition ParticleHistograms.cxx:41

◆ fill() [2/3]

void egammaMonitoring::ParticleHistograms::fill ( const xAOD::IParticle & egamma,
float mu )
virtual

Reimplemented from egammaMonitoring::ParticleHistograms.

Definition at line 24 of file ParticleHistograms.cxx.

45 {
46
47 double pT = egamma.pt()*1e-3;
48 if (pT > 0)
49 histoMap["pT"]->Fill(pT);
50 histoMap["eta"]->Fill(egamma.eta());
51 histoMap["phi"]->Fill(egamma.phi());
52
53 if (pT > 15) {
54 histoMap["eta_15GeV"]->Fill(egamma.eta());
55 histoMap["phi_15GeV"]->Fill(egamma.phi());
56 }
57
58 if (m_hasFwd && std::abs(egamma.eta()) < 4.) {
59 histoMap["pT_eta4"]->Fill(pT);
60 histoMap["phi_eta4"]->Fill(egamma.phi());
61 }
62
63 }
egammaMonitoring::IHistograms::histoMap
std::map< std::string, TH1D * > histoMap
Definition IHistograms.h:42

◆ fill() [3/3]

void TruthElectronHistograms::fill ( const xAOD::TruthParticle * truth,
const xAOD::Electron * el = nullptr )

Definition at line 68 of file TruthElectronHistograms.cxx.

68 {
69
70 ParticleHistograms::fill(*truth);
71
72 if (!electron || m_reducedHistSet) return;
73
74 const xAOD::TrackParticle* track = electron->trackParticle();
75
76 // This can happen if we use it for forwardElectron
77 if (!track) return;
78
79 bool has2DHis = !histoMap2D.empty();
80
81 float dphires2(0.);
82 float dphi2(0.);
83 float deta2(0);
84
85 if (electron->trackCaloMatchValue(dphires2, xAOD::EgammaParameters::deltaPhiRescaled2)) {
86 histoMap["deltaPhiRescaled2"]->Fill(dphires2);
87 if (has2DHis) histoMap2D["eta_deltaPhiRescaled2"]->Fill(electron->eta(),dphires2);
88 }
89 if (electron->trackCaloMatchValue(dphi2, xAOD::EgammaParameters::deltaPhi2)) {
90 histoMap["deltaPhi2"]->Fill(dphi2);
91 if (has2DHis) histoMap2D["eta_deltaPhi2"]->Fill(electron->eta(),dphi2);
92 }
93 if (electron->trackCaloMatchValue(deta2, xAOD::EgammaParameters::deltaEta2)) {
94 histoMap["deltaEta2"]->Fill(deta2);
95 if (has2DHis) histoMap2D["eta_deltaEta2"]->Fill(electron->eta(),deta2);
96 }
97
98 float d0 = track->d0();
99 float reco_qp = track->qOverP();
100 float truth_qp = truth->charge()/(truth->pt()*cosh(truth->eta()));
101 float vard0 = track->definingParametersCovMatrix()(0, 0);
102
103 if (vard0 > 0) {
104 histoMap["d0Oversigmad0"]->Fill(d0/sqrtf(vard0));
105 if (has2DHis) histoMap2D["eta_d0Oversigmad0"]->Fill(electron->eta(),d0/sqrtf(vard0));
106 }
107
108 if (truth_qp > 0) {
109 histoMap["qOverp_resolution"]->Fill((reco_qp-truth_qp)/truth_qp);
110 if (has2DHis) histoMap2D["eta_qOverp_resolution"]->Fill(electron->eta(),(reco_qp-truth_qp)/truth_qp);
111 }
112
113
114}
xAOD::TruthParticle_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition TruthParticle_v1.cxx:161
xAOD::TruthParticle_v1::charge
double charge() const
Physical charge.
xAOD::TruthParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition TruthParticle_v1.cxx:169
xAOD::EgammaParameters::deltaEta2
@ deltaEta2
difference between the cluster eta (second sampling) and the eta of the track extrapolated to the sec...
Definition EgammaEnums.h:188
xAOD::EgammaParameters::deltaPhiRescaled2
@ deltaPhiRescaled2
difference between the cluster phi (second sampling) and the phi of the track extrapolated to the sec...
Definition EgammaEnums.h:226
xAOD::EgammaParameters::deltaPhi2
@ deltaPhi2
difference between the cluster phi (second sampling) and the phi of the track extrapolated to the sec...
Definition EgammaEnums.h:205
xAOD::TrackParticle
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Definition Event/xAOD/xAODTracking/xAODTracking/TrackParticle.h:13

◆ hasFwd()

void egammaMonitoring::ParticleHistograms::hasFwd ( bool a = true)
inlineinherited

Definition at line 26 of file ParticleHistograms.h.

26{ m_hasFwd = a; }
a
static Double_t a
Definition LArPhysWaveHECTool.cxx:38

◆ IHistograms() [1/2]

egammaMonitoring::IHistograms::IHistograms ( )
inherited

◆ IHistograms() [2/2]

egammaMonitoring::IHistograms::IHistograms ( std::string name,
std::string title,
std::string folder,
SmartIF< ITHistSvc > rootHistSvc )
inlineinherited

Definition at line 26 of file IHistograms.h.

30 :
31 m_name(std::move(std::move(name))),
32 m_title(std::move(std::move(title))),
33 m_folder(std::move(std::move(folder))),
34 m_rootHistSvc(std::move(rootHistSvc)) {}
egammaMonitoring::IHistograms::m_rootHistSvc
SmartIF< ITHistSvc > m_rootHistSvc
Definition IHistograms.h:48

◆ initializePlots() [1/2]

StatusCode TruthElectronHistograms::initializePlots ( )
overridevirtual

Reimplemented from egammaMonitoring::ParticleHistograms.

Definition at line 18 of file TruthElectronHistograms.cxx.

18 {
19 return initializePlots (false);
20}

◆ initializePlots() [2/2]

StatusCode TruthElectronHistograms::initializePlots ( bool reducedHistSet)

Definition at line 22 of file TruthElectronHistograms.cxx.

22 {
23
24 const char* fN = m_name.c_str();
25
26 if (!reducedHistSet) {
27 histoMap["deltaPhi2"] = new TH1D(Form("%s_deltaPhi2",fN), ";deltaPhi2; Events", 40, -0.06, 0.06);
28 histoMap["deltaEta2"] = new TH1D(Form("%s_deltaEta2",fN), ";deltaEta2; Events", 40, -0.04, 0.04);
29 histoMap["deltaPhiRescaled2"] = new TH1D(Form("%s_deltaPhiRescaled2",fN), ";deltaPhiRescaled2; Events", 40, -0.04, 0.04);
30
31 histoMap["d0Oversigmad0"] = new TH1D(Form("%s_d0Oversigmad0",fN), "; d0Oversigmad0; Events", 40, -10, 10);
32 histoMap["qOverp_resolution"] = new TH1D(Form("%s_qOverp_resolution",fN), ";(q/P reco - q/P truth)/ q/p truth; Events", 60, -1, 1.5);
33
34 ATH_CHECK(m_rootHistSvc->regHist(m_folder+"deltaPhi2", histoMap["deltaPhi2"]));
35 ATH_CHECK(m_rootHistSvc->regHist(m_folder+"deltaEta2", histoMap["deltaEta2"]));
36 ATH_CHECK(m_rootHistSvc->regHist(m_folder+"deltaPhiRescaled2", histoMap["deltaPhiRescaled2"]));
37 ATH_CHECK(m_rootHistSvc->regHist(m_folder+"d0Oversigmad0", histoMap["d0Oversigmad0"]));
38 ATH_CHECK(m_rootHistSvc->regHist(m_folder+"qOverp_resolution", histoMap["qOverp_resolution"]));
39
40 // 2D only for truthPromptElectronWithRecoTrack (temporary)
41 if (m_name == "truthPromptElectronWithRecoTrack") {
42 histoMap2D["eta_deltaPhi2"] = new TH2D(Form("%s_eta_deltaPhi2",fN),
43 ";#eta;#Delta#phi_{2}; Events", 90, -4.5, 4.5, 40, -0.06, 0.06);
44 histoMap2D["eta_deltaEta2"] = new TH2D(Form("%s_eta_deltaEta2",fN),
45 ";#eta;#Delta#eta_{2}; Events", 90, -4.5, 4.5, 40, -0.06, 0.06);
46 histoMap2D["eta_deltaPhiRescaled2"] = new TH2D(Form("%s_eta_deltaPhiRescaled2",fN),
47 ";#eta;#Delta#phi_{2}^{Rescaled}; Events", 90, -4.5, 4.5, 40, -0.06, 0.06);
48 histoMap2D["eta_d0Oversigmad0"] = new TH2D(Form("%s_eta_d0Oversigmad0",fN),
49 ";#eta;d_{0}/#sigma_{d_{0}}; Events", 90, -4.5, 4.5, 40, -10, 10);
50 histoMap2D["eta_qOverp_resolution"] = new TH2D(Form("%s_eta_qOverp_resolution",fN),
51 ";#eta;(q/P_{reco})/(q/P_{truth}) -1; Events", 90, -4.5, 4.5, 60, -1, 1.5);
52
53 for (const auto& e : histoMap2D) {
54 ATH_CHECK(m_rootHistSvc->regHist(m_folder+e.first, e.second));
55 }
56 }
57 }
58
59 ATH_CHECK(ParticleHistograms::initializePlots());
60
61 m_reducedHistSet = reducedHistSet;
62
63 return StatusCode::SUCCESS;
64
65}
ATH_CHECK
#define ATH_CHECK
Evaluate an expression and check for errors.
Definition AthCheckMacros.h:40

Member Data Documentation

◆ histoMap

std::map<std::string, TH1D*> egammaMonitoring::IHistograms::histoMap
inherited

Definition at line 42 of file IHistograms.h.

◆ histoMap2D

std::map<std::string, TH2D*> egammaMonitoring::TruthElectronHistograms::histoMap2D

Definition at line 33 of file TruthElectronHistograms.h.

◆ m_folder

std::string egammaMonitoring::IHistograms::m_folder
protectedinherited

Definition at line 47 of file IHistograms.h.

◆ m_hasFwd

bool egammaMonitoring::ParticleHistograms::m_hasFwd = false
privateinherited

Definition at line 29 of file ParticleHistograms.h.

◆ m_name

std::string egammaMonitoring::IHistograms::m_name
protectedinherited

Definition at line 45 of file IHistograms.h.

◆ m_reducedHistSet

bool egammaMonitoring::TruthElectronHistograms::m_reducedHistSet {}
private

Definition at line 37 of file TruthElectronHistograms.h.

37{};

◆ m_rootHistSvc

SmartIF<ITHistSvc> egammaMonitoring::IHistograms::m_rootHistSvc
protectedinherited

Definition at line 48 of file IHistograms.h.

◆ m_title

std::string egammaMonitoring::IHistograms::m_title
protectedinherited

Definition at line 46 of file IHistograms.h.

The documentation for this class was generated from the following files:
Generated on for ATLAS Offline Software by doxygen 1.16.1