Resolve "pi0 resolution for ECal barrel"

Closes #19 (closed)

benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx

@@ -8,11 +8,11 @@
 
 #include "dd4pod/Geant4ParticleCollection.h"
 #include "dd4pod/CalorimeterHitCollection.h"
-
+/*
 #include "benchmark.h"
 #include "mt.h"
 #include "util.h"
-
+*/
 #include "TCanvas.h"
 #include "TStyle.h"
 #include "TMath.h"
@@ -42,7 +42,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
   //TODO: Change test_tag to something else
   std:string detector = "Barrel_emcal";
   double resolutionTarget = 0.1;
-  /*
+
   eic::util::Test pion0_Energy_resolution{
       {{"name", fmt::format("{}_energy_resolution", test_tag)},
        {"title", "Pion0 Energy resolution"},
@@ -50,7 +50,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
         fmt::format("Pion0 energy resolution with {}, estimated using a Gaussian fit.", detector)},
        {"quantity", "resolution (in %)"},
        {"target", std::to_string(resolutionTarget)}}};
-*/
+
 
   ROOT::EnableImplicitMT();
   ROOT::RDataFrame d0("events", input_fname);
@@ -105,7 +105,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
     return result;
   };
 
-  // Relative Energy Resolution = (Esampling - Ethrown)/Ethrownnnn
+  // Relative Energy Resolution = (Esampling - Ethrown)/Ethrown
   auto eResol_rel = [](const std::vector<double>& sampled, const std::vector<double>& thrown) {
     std::vector<double> result;
     auto it_sam = sampled.cbegin();