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	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, 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
ITHistSvc *	m_rootHistSvc = nullptr

Private Attributes
bool	m_reducedHistSet
bool	m_hasFwd = false

Detailed Description

Definition at line 16 of file TruthElectronHistograms.h.

Constructor & Destructor Documentation

~TruthElectronHistograms()

virtual egammaMonitoring::TruthElectronHistograms::~TruthElectronHistograms ( )

inlinevirtual

Definition at line 29 of file TruthElectronHistograms.h.

{};

Member Function Documentation

fill() [1/5]

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

virtualinherited

Implements egammaMonitoring::IHistograms.

Reimplemented in egammaMonitoring::TruthPhotonHistograms.

Definition at line 41 of file ParticleHistograms.cxx.

                                                           {
    ParticleHistograms::fill(egamma,0.);
  }

fill() [2/5]

void egammaMonitoring::ParticleHistograms::fill

Definition at line 23 of file ParticleHistograms.cxx.

                                                           {
    ParticleHistograms::fill(egamma,0.);
  }

fill() [3/5]

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

virtualinherited

Implements egammaMonitoring::IHistograms.

Reimplemented in egammaMonitoring::TruthPhotonHistograms.

Definition at line 45 of file ParticleHistograms.cxx.

                                                                    {
 
    double pT = egamma.pt()*1e-3;
    if (pT > 0)
      histoMap["pT"]->Fill(pT);
    histoMap["eta"]->Fill(egamma.eta());
    histoMap["phi"]->Fill(egamma.phi());
 
    if (pT > 15) {
      histoMap["eta_15GeV"]->Fill(egamma.eta());
      histoMap["phi_15GeV"]->Fill(egamma.phi());
    }
 
    if (m_hasFwd && std::abs(egamma.eta()) < 4.) {
      histoMap["pT_eta4"]->Fill(pT);
      histoMap["phi_eta4"]->Fill(egamma.phi());
    }
 
  }

fill() [4/5]

void egammaMonitoring::ParticleHistograms::fill

Definition at line 24 of file ParticleHistograms.cxx.

                                                                    {
 
    double pT = egamma.pt()*1e-3;
    if (pT > 0)
      histoMap["pT"]->Fill(pT);
    histoMap["eta"]->Fill(egamma.eta());
    histoMap["phi"]->Fill(egamma.phi());
 
    if (pT > 15) {
      histoMap["eta_15GeV"]->Fill(egamma.eta());
      histoMap["phi_15GeV"]->Fill(egamma.phi());
    }
 
    if (m_hasFwd && std::abs(egamma.eta()) < 4.) {
      histoMap["pT_eta4"]->Fill(pT);
      histoMap["phi_eta4"]->Fill(egamma.phi());
    }
 
  }

fill() [5/5]

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

Definition at line 67 of file TruthElectronHistograms.cxx.

                                                                                               {
 
  ParticleHistograms::fill(*truth);
 
  if (!electron || m_reducedHistSet) return;
 
  const xAOD::TrackParticle* track = electron->trackParticle();
 
  // This can happen if we use it for forwardElectron
  if (!track) return;
 
  bool has2DHis = !histoMap2D.empty();
 
  float dphires2(0.);
  float dphi2(0.);
  float deta2(0);
 
  if (electron->trackCaloMatchValue(dphires2, xAOD::EgammaParameters::deltaPhiRescaled2)) {
    histoMap["deltaPhiRescaled2"]->Fill(dphires2);
    if (has2DHis) histoMap2D["eta_deltaPhiRescaled2"]->Fill(electron->eta(),dphires2);
  }
  if (electron->trackCaloMatchValue(dphi2, xAOD::EgammaParameters::deltaPhi2)) {
    histoMap["deltaPhi2"]->Fill(dphi2);
    if (has2DHis) histoMap2D["eta_deltaPhi2"]->Fill(electron->eta(),dphi2);
  }
  if (electron->trackCaloMatchValue(deta2, xAOD::EgammaParameters::deltaEta2)) {
    histoMap["deltaEta2"]->Fill(deta2);
    if (has2DHis) histoMap2D["eta_deltaEta2"]->Fill(electron->eta(),deta2);
  }
 
  float d0 = track->d0(); 
  float reco_qp = track->qOverP();
  float truth_qp = truth->charge()/(truth->pt()*cosh(truth->eta()));
  float vard0 = track->definingParametersCovMatrix()(0, 0);
 
  if (vard0 > 0) {
    histoMap["d0Oversigmad0"]->Fill(d0/sqrtf(vard0));
    if (has2DHis) histoMap2D["eta_d0Oversigmad0"]->Fill(electron->eta(),d0/sqrtf(vard0));
  }
  
  if (truth_qp > 0) {
    histoMap["qOverp_resolution"]->Fill((reco_qp-truth_qp)/truth_qp);
    if (has2DHis) histoMap2D["eta_qOverp_resolution"]->Fill(electron->eta(),(reco_qp-truth_qp)/truth_qp);
  }
  
 
}

hasFwd()

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

inlineinherited

Definition at line 26 of file ParticleHistograms.h.

{ m_hasFwd = a; }

IHistograms() [1/2]

egammaMonitoring::IHistograms::IHistograms

inherited

IHistograms() [2/2]

egammaMonitoring::IHistograms::IHistograms

inlineinherited

Definition at line 23 of file IHistograms.h.

          :
      m_name(std::move(std::move(name))),
      m_title(std::move(std::move(title))),
      m_folder(std::move(std::move(folder))),
      m_rootHistSvc(rootHistSvc) {}

initializePlots() [1/2]

StatusCode TruthElectronHistograms::initializePlots ( )

overridevirtual

Reimplemented from egammaMonitoring::ParticleHistograms.

Definition at line 17 of file TruthElectronHistograms.cxx.

                                                    {
  return initializePlots (false);
}

initializePlots() [2/2]

StatusCode TruthElectronHistograms::initializePlots

(

bool

reducedHistSet

)

Definition at line 21 of file TruthElectronHistograms.cxx.

                                                                       {
 
  const char* fN = m_name.c_str();
 
  if (!reducedHistSet) {
    histoMap["deltaPhi2"] = new TH1D(Form("%s_deltaPhi2",fN), ";deltaPhi2; Events", 40, -0.06, 0.06);
    histoMap["deltaEta2"] = new TH1D(Form("%s_deltaEta2",fN), ";deltaEta2; Events", 40, -0.04, 0.04);
    histoMap["deltaPhiRescaled2"] = new TH1D(Form("%s_deltaPhiRescaled2",fN), ";deltaPhiRescaled2; Events", 40, -0.04, 0.04);
 
    histoMap["d0Oversigmad0"] = new TH1D(Form("%s_d0Oversigmad0",fN), "; d0Oversigmad0; Events", 40, -10, 10);
    histoMap["qOverp_resolution"] = new TH1D(Form("%s_qOverp_resolution",fN), ";(q/P reco - q/P truth)/ q/p truth; Events", 60, -1, 1.5);
 
    ATH_CHECK(m_rootHistSvc->regHist(m_folder+"deltaPhi2", histoMap["deltaPhi2"]));
    ATH_CHECK(m_rootHistSvc->regHist(m_folder+"deltaEta2", histoMap["deltaEta2"]));
    ATH_CHECK(m_rootHistSvc->regHist(m_folder+"deltaPhiRescaled2", histoMap["deltaPhiRescaled2"]));
    ATH_CHECK(m_rootHistSvc->regHist(m_folder+"d0Oversigmad0", histoMap["d0Oversigmad0"]));
    ATH_CHECK(m_rootHistSvc->regHist(m_folder+"qOverp_resolution", histoMap["qOverp_resolution"]));
 
    // 2D only for truthPromptElectronWithRecoTrack (temporary)
    if (m_name == "truthPromptElectronWithRecoTrack") {
      histoMap2D["eta_deltaPhi2"] = new TH2D(Form("%s_eta_deltaPhi2",fN),
    ";#eta;#Delta#phi_{2}; Events", 90, -4.5, 4.5, 40, -0.06, 0.06);
      histoMap2D["eta_deltaEta2"] = new TH2D(Form("%s_eta_deltaEta2",fN),
    ";#eta;#Delta#eta_{2}; Events", 90, -4.5, 4.5, 40, -0.06, 0.06);
      histoMap2D["eta_deltaPhiRescaled2"] = new TH2D(Form("%s_eta_deltaPhiRescaled2",fN),
    ";#eta;#Delta#phi_{2}^{Rescaled}; Events", 90, -4.5, 4.5, 40, -0.06, 0.06);
      histoMap2D["eta_d0Oversigmad0"] = new TH2D(Form("%s_eta_d0Oversigmad0",fN),
    ";#eta;d_{0}/#sigma_{d_{0}}; Events", 90, -4.5, 4.5, 40, -10, 10);
      histoMap2D["eta_qOverp_resolution"] = new TH2D(Form("%s_eta_qOverp_resolution",fN),
    ";#eta;(q/P_{reco})/(q/P_{truth}) -1; Events", 90, -4.5, 4.5, 60, -1, 1.5);
 
      for (const auto& e : histoMap2D) {
    ATH_CHECK(m_rootHistSvc->regHist(m_folder+e.first, e.second));
      }
    }
  }
  
  ATH_CHECK(ParticleHistograms::initializePlots());
 
  m_reducedHistSet = reducedHistSet;
 
  return StatusCode::SUCCESS;
 
} 

Member Data Documentation

histoMap

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

inherited

Definition at line 39 of file IHistograms.h.

histoMap2D

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

Definition at line 31 of file TruthElectronHistograms.h.

m_folder

std::string egammaMonitoring::IHistograms::m_folder

protectedinherited

Definition at line 44 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 42 of file IHistograms.h.

m_reducedHistSet

bool egammaMonitoring::TruthElectronHistograms::m_reducedHistSet

private

Definition at line 35 of file TruthElectronHistograms.h.

m_rootHistSvc

ITHistSvc* egammaMonitoring::IHistograms::m_rootHistSvc = nullptr

protectedinherited

Definition at line 45 of file IHistograms.h.

m_title

std::string egammaMonitoring::IHistograms::m_title

protectedinherited

Definition at line 43 of file IHistograms.h.

---

The documentation for this class was generated from the following files:

• TruthElectronHistograms.h
• TruthElectronHistograms.cxx