egammaMonitoring::ClusterHistograms Class Reference

#include <ClusterHistograms.h>

Collaboration diagram for egammaMonitoring::ClusterHistograms:

Public Member Functions
	ClusterHistograms (std::string name, std::string title, std::string folder, SmartIF< ITHistSvc > rootHistSvc)
StatusCode	initializePlots ()
void	fill (const xAOD::Egamma &egamma)
void	fill (const xAOD::Egamma &egamma, float mu)

Public Attributes
std::map< std::string, TH2D * >	histo2DMap
std::map< std::string, TProfile * >	profileMap

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

Detailed Description

Definition at line 21 of file ClusterHistograms.h.

Constructor & Destructor Documentation

ClusterHistograms()

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

inline

Definition at line 25 of file ClusterHistograms.h.

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

Member Function Documentation

fill() [1/2]

void ClusterHistograms::fill

(

const xAOD::Egamma &

egamma

)

Definition at line 69 of file ClusterHistograms.cxx.

                                                     {
  fill(egamma,0.);
}

fill() [2/2]

void ClusterHistograms::fill	(	const xAOD::Egamma &	egamma,
		float	mu = 0
	)

Definition at line 73 of file ClusterHistograms.cxx.

                                                                   {
 
  const xAOD::CaloCluster *cluster = egamma.caloCluster();
  
  const xAOD::TruthParticle *truth_egamma = xAOD::TruthHelpers::getTruthParticle(egamma);
 
  if ( !truth_egamma ) return;
 
  const auto Ereco = cluster->rawE();
  const auto Etruth = truth_egamma->e();
  const auto Eres = (Ereco - Etruth)/Etruth;
 
  profileMap["Eraw_Etruth_vs_Etruth_profile"]->Fill(Etruth/1000,Ereco/Etruth);
  profileMap["Eraw_Etruth_vs_eta_profile"]->Fill(truth_egamma->eta(),Ereco/Etruth);
  histo2DMap["mu_energy_resolution_2D"]->Fill(mu,Eres);
 
  const CaloClusterCellLink* cellLinks = cluster->getCellLinks();
  
  if ( !cellLinks ) return;
 
  std::map<int, int > cells_per_layer;
  std::map<int, int > cells_per_sam;
 
  for (const CaloCell* cell : *cellLinks) {
    if (cell) {
      int layer = cell->caloDDE()->getLayer();
      cells_per_layer[layer]++;
      int sam = cell->caloDDE()->getSampling();
      cells_per_sam[sam]++;
    }
  }
 
  auto associatedTopoCluster = xAOD::EgammaHelpers::getAssociatedTopoClusters(cluster);
  profileMap["number_topocluster_vs_e_profile"]->Fill(Etruth/1000, associatedTopoCluster.size());
  profileMap["number_topocluster_vs_eta_profile"]->Fill(truth_egamma->eta(), associatedTopoCluster.size());
 
  for (auto const& x : cells_per_layer) {
   profileMap["number_cell_in_layer"]->Fill(x.second,x.first);
   profileMap["number_cells_vs_e_in_layer_"+std::to_string(x.first)+"_profile"]->Fill(Etruth/1000, x.second);
   profileMap["number_cells_vs_eta_in_layer_"+std::to_string(x.first)+"_profile"]->Fill(truth_egamma->eta(), x.second);
  }
  for (auto const& x : cells_per_sam) {
    profileMap[Form("number_cells_vs_eta_in_sampling_%i_profile",x.first)]->Fill(truth_egamma->eta(), x.second);
  }
  
}

initializePlots()

StatusCode ClusterHistograms::initializePlots

(

)

Definition at line 20 of file ClusterHistograms.cxx.

                                              {
 
  const char* fN = m_name.c_str();
  
  for (int is = 0; is < 4; is++) {
    TString pNm = Form("number_cells_vs_eta_in_layer_%i_profile",is);
    TString pN  = Form("%s_number_cells_vs_eta_in_layer_%i_profile",fN,is);
    profileMap[pNm.Data()] = new TProfile(pN.Data(), "Number of cells;truth #eta", 90,-4.5,4.5, 0, 100);
    pN = Form("%s%s",m_folder.c_str(),pNm.Data());
    ATH_CHECK(m_rootHistSvc->regHist(pN.Data(), profileMap[pNm.Data()]));
 
    pNm = Form("number_cells_vs_e_in_layer_%i_profile",is);
    pN  = Form("%s_number_cells_vs_e_in_layer_%i_profile",fN,is);
    profileMap[pNm.Data()] = new TProfile(pN.Data(), "Number of cells;truth E [GeV]", 60, 0, 300, 0, 100);
    pN = Form("%s%s",m_folder.c_str(),pNm.Data());
    ATH_CHECK(m_rootHistSvc->regHist(pN.Data(), profileMap[pNm.Data()]));
  }
 
  for (int is = 0; is < 24; is++) {
    TString pNm = Form("number_cells_vs_eta_in_sampling_%i_profile",is);
    TString pN = Form("%s_number_cells_vs_eta_in_sampling_%i_profile",fN,is);
    profileMap[pNm.Data()] = new TProfile(pN.Data(), "Number of cells;truth #eta", 90,-4.5,4.5, 0,100);
    pN = Form("%s%s",m_folder.c_str(),pNm.Data());
    ATH_CHECK(m_rootHistSvc->regHist(pN.Data(), profileMap[pNm.Data()]));
  }
 
  profileMap["Eraw_Etruth_vs_Etruth_profile"] = new TProfile(Form("%s_Eraw_Etruth_vs_Etruth_profile",fN), ";E^{truth};E^{raw}/E^{truth}", 100, 0., 200., 0.5, 1.5);
  profileMap["Eraw_Etruth_vs_eta_profile"]    = new TProfile(Form("%s_Eraw_Etruth_vs_eta_profile",fN),    ";truth #eta;E^{raw}/E^{truth}", 90, -4.5, 4.5, 0.5, 1.5);
 
  profileMap["number_topocluster_vs_e_profile"]   = new TProfile(Form("%s_number_topocluster_vs_e_profile",fN),  "Number of topocluster;truth E [GeV]", 60, 0, 300, -0.5, 14.5);
  profileMap["number_topocluster_vs_eta_profile"] = new TProfile(Form("%s_number_topocluster_vs_eta_profile",fN), "Number of topocluster;truth #eta", 90,-4.5,4.5, -0.5, 14.5);
 
  // do not undertstand this histo...
  profileMap["number_cell_in_layer"]    = new TProfile(Form("%s_number_cell_in_layer",fN), ";Number of cells;Layer", 100, 0, 200,0, 4);
  histo2DMap["mu_energy_resolution_2D"] = new TH2D(Form("%s_mu_energy_resolution_2D",fN), ";<#mu>; Energy Resolution", 5, 0, 80, 20, -1, 1);
 
  ATH_CHECK(m_rootHistSvc->regHist(m_folder+"Eraw_Etruth_vs_Etruth_profile", profileMap["Eraw_Etruth_vs_Etruth_profile"]));
  ATH_CHECK(m_rootHistSvc->regHist(m_folder+"Eraw_Etruth_vs_eta_profile", profileMap["Eraw_Etruth_vs_eta_profile"]));
 
  ATH_CHECK(m_rootHistSvc->regHist(m_folder+"number_topocluster_vs_e_profile", profileMap["number_topocluster_vs_e_profile"]));
  ATH_CHECK(m_rootHistSvc->regHist(m_folder+"number_topocluster_vs_eta_profile", profileMap["number_topocluster_vs_eta_profile"]));
 
  ATH_CHECK(m_rootHistSvc->regHist(m_folder+"number_cell_in_layer", profileMap["number_cell_in_layer"]));
  ATH_CHECK(m_rootHistSvc->regHist(m_folder+"mu_energy_resolution_2D", histo2DMap["mu_energy_resolution_2D"]));
 
  return StatusCode::SUCCESS;
  
} // initializePlots

Member Data Documentation

histo2DMap

std::map<std::string, TH2D* > egammaMonitoring::ClusterHistograms::histo2DMap

Definition at line 35 of file ClusterHistograms.h.

m_folder

std::string egammaMonitoring::ClusterHistograms::m_folder

protected

Definition at line 45 of file ClusterHistograms.h.

m_name

std::string egammaMonitoring::ClusterHistograms::m_name

protected

Definition at line 43 of file ClusterHistograms.h.

m_rootHistSvc

SmartIF<ITHistSvc> egammaMonitoring::ClusterHistograms::m_rootHistSvc

protected

Definition at line 46 of file ClusterHistograms.h.

m_title

std::string egammaMonitoring::ClusterHistograms::m_title

protected

Definition at line 44 of file ClusterHistograms.h.

profileMap

std::map<std::string, TProfile* > egammaMonitoring::ClusterHistograms::profileMap

Definition at line 36 of file ClusterHistograms.h.

---

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

• ClusterHistograms.h
• ClusterHistograms.cxx