diff --git a/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx b/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx index de2a14ca2e0e490674c61da6ca053fdf904892f0..eef7ac6b4e877a3e3a85309b722a1b92805e0bf3 100644 --- a/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx +++ b/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx @@ -56,7 +56,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal ROOT::RDataFrame d0("events", input_fname); // Sampling Fraction - double samp_frac = 1.0; + double samp_frac = 0.0136; // Thrown Energy [GeV] auto Ethr = [](std::vector<dd4pod::Geant4ParticleData> const& input) { @@ -132,6 +132,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal auto hpid = d1.Histo1D({"hpid", "PID; PID; Count", 100, -220, 220}, "pid"); // Set sampling Fraction, ideally this will be taken from a json file + cout << "The sampling fraction mean is " << samp_frac << endl; samp_frac = hfsam -> GetMean(); cout << "The sampling fraction mean is " << samp_frac << endl;