Trying with const double meanE

benchmarks/barrel_ecal/scripts/emcal_barrel_pi0_analysis.cxx

@@ -49,7 +49,7 @@ void emcal_barrel_pi0_analysis(const char* input_fname = "sim_output/sim_emcal_b
   double samp_frac = 1.0;
 
   // Thrown Energy [GeV]
-  double meanE = 5; // Calculated later
+  //double meanE = 5; // Calculated later
   auto Ethr = [](std::vector<dd4pod::Geant4ParticleData> const& input) {
     return TMath::Sqrt(input[2].psx*input[2].psx + input[2].psy*input[2].psy + input[2].psz*input[2].psz + input[2].mass*input[2].mass);
   };
@@ -112,7 +112,7 @@ void emcal_barrel_pi0_analysis(const char* input_fname = "sim_output/sim_emcal_b
 
   // Gather Sampling fraction and mean Energy, ideally this will be taken from a json file 
   samp_frac = hfsam->GetMean();
-  meanE     = hEthr->GetMean(); 
+  const double meanE     = hEthr->GetMean(); 
 
   // Event Counts
   auto nevents_thrown = d1.Count();