Merge branch 'ARTplots_egamma' into 'master'

egamma ART plots

See merge request atlas/athena!39556

Reconstruction/egamma/egammaValidation/src/ClusterHistograms.cxx

@@ -17,6 +17,14 @@ StatusCode ClusterHistograms::initializePlots() {
   histo2DMap["number_cells_vs_et_in_layer_2_2D"] = (new TH2D(Form("%s_%s",m_name.c_str(),"number_cells_vs_et_in_layer_2_2D"), "Number of cells;E_{T}", 60,0,300,50,   0,100.));
   histo2DMap["number_cells_vs_et_in_layer_3_2D"] = (new TH2D(Form("%s_%s",m_name.c_str(),"number_cells_vs_et_in_layer_3_2D"), "Number of cells;E_{T}", 60,0,300,50,   0,100.));
 
+  histo2DMap["number_cells_vs_eta_in_layer_0_2D"] = (new TH2D(Form("%s_%s",m_name.c_str(),"number_cells_vs_eta_in_layer_0_2D"), "Number of cells;truth #eta", 40,-3,3, 50,0,100));
+  histo2DMap["number_cells_vs_eta_in_layer_1_2D"] = (new TH2D(Form("%s_%s",m_name.c_str(),"number_cells_vs_eta_in_layer_1_2D"), "Number of cells;truth #eta", 40,-3,3, 50,0,100));
+  histo2DMap["number_cells_vs_eta_in_layer_2_2D"] = (new TH2D(Form("%s_%s",m_name.c_str(),"number_cells_vs_eta_in_layer_2_2D"), "Number of cells;truth #eta", 40,-3,3, 50,0,100));
+  histo2DMap["number_cells_vs_eta_in_layer_3_2D"] = (new TH2D(Form("%s_%s",m_name.c_str(),"number_cells_vs_eta_in_layer_3_2D"), "Number of cells;truth #eta", 40,-3,3, 50,0,100));
+
+  histo2DMap["number_topocluster_vs_et_2D"] = (new TH2D(Form("%s_%s",m_name.c_str(),"number_topocluster_vs_et_2D"), "Number of topocluster;E_{T}", 60,0,300, 15,0,15));
+  histo2DMap["number_topocluster_vs_eta_2D"] = (new TH2D(Form("%s_%s",m_name.c_str(),"number_topocluster_vs_eta_2D"), "Number of topocluster;truth #eta", 40,-3,3, 15,0,15));
+
   histo2DMap["number_cell_in_layer"] = (new TH2D(Form("%s_%s",m_name.c_str(),"number_cell_in_layer"), ";Number of cells;Layer",100,0,200, 4,0,4));  
   histo2DMap["mu_energy_resolution_2D"] = (new TH2D(Form("%s_%s",m_name.c_str(),"mu_energy_resolution_2D"), ";<#mu>; Energy Resolution", 5,0,80,20,-1,1));
 
@@ -28,6 +36,14 @@ StatusCode ClusterHistograms::initializePlots() {
   ATH_CHECK(m_rootHistSvc->regHist(m_folder+"number_cells_vs_et_in_layer_2_2D", histo2DMap["number_cells_vs_et_in_layer_2_2D"]));
   ATH_CHECK(m_rootHistSvc->regHist(m_folder+"number_cells_vs_et_in_layer_3_2D", histo2DMap["number_cells_vs_et_in_layer_3_2D"]));
  
+  ATH_CHECK(m_rootHistSvc->regHist(m_folder+"number_cells_vs_eta_in_layer_0_2D", histo2DMap["number_cells_vs_eta_in_layer_0_2D"]));
+  ATH_CHECK(m_rootHistSvc->regHist(m_folder+"number_cells_vs_eta_in_layer_1_2D", histo2DMap["number_cells_vs_eta_in_layer_1_2D"]));
+  ATH_CHECK(m_rootHistSvc->regHist(m_folder+"number_cells_vs_eta_in_layer_2_2D", histo2DMap["number_cells_vs_eta_in_layer_2_2D"]));
+  ATH_CHECK(m_rootHistSvc->regHist(m_folder+"number_cells_vs_eta_in_layer_3_2D", histo2DMap["number_cells_vs_eta_in_layer_3_2D"]));
+
+  ATH_CHECK(m_rootHistSvc->regHist(m_folder+"number_topocluster_vs_et_2D", histo2DMap["number_topocluster_vs_et_2D"]));
+  ATH_CHECK(m_rootHistSvc->regHist(m_folder+"number_topocluster_vs_eta_2D", histo2DMap["number_topocluster_vs_eta_2D"]));
+
  
   ATH_CHECK(m_rootHistSvc->regHist(m_folder+"number_cell_in_layer", histo2DMap["number_cell_in_layer"]));
   ATH_CHECK(m_rootHistSvc->regHist(m_folder+"mu_energy_resolution_2D", histo2DMap["mu_energy_resolution_2D"]));
@@ -80,9 +96,14 @@ void ClusterHistograms::fill(const xAOD::Egamma& egamma, float mu = 0) {
 
   }  
 
+  auto associatedTopoCluster = xAOD::EgammaHelpers::getAssociatedTopoClusters(cluster);
+  histo2DMap["number_topocluster_vs_et_2D"]->Fill(truth_egamma->e()/1000, associatedTopoCluster.size());
+  histo2DMap["number_topocluster_vs_eta_2D"]->Fill(truth_egamma->eta(), associatedTopoCluster.size());
+
   for (auto const& x : cells_per_layer) {
    histo2DMap["number_cell_in_layer"]->Fill(x.second,x.first);
    histo2DMap["number_cells_vs_et_in_layer_"+std::to_string(x.first)+"_2D"]->Fill(truth_egamma->e()/1000, x.second);
+   histo2DMap["number_cells_vs_eta_in_layer_"+std::to_string(x.first)+"_2D"]->Fill(truth_egamma->eta(), x.second);
 
   }