diff --git a/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx b/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx index 486755fc123679e543ef35096d77fd3899e8b0d5..52594492407cce8215ff6197c0ab4c1d42be313a 100644 --- a/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx +++ b/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx @@ -213,14 +213,13 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal c6->SaveAs("results/emcal_barrel_pions_dau.pdf"); //Energy Resolution Calculation - auto hdE = d2.Histo1D({"hdE", "dE; dE[GeV]; Events", 20, -2.5, 2.5}, "dE");//changed from 100 - auto hdE_rel = d2.Histo1D({"hdE_rel", "dE Relative; dE Relative; Events", 20, -2.5, 2.5}, "dE_rel");//changed from 100 - hdE->Fit("gaus", "", "", -2.5, 2.5); + auto hdE = d2.Histo1D({"hdE", "dE; dE[GeV]; Events", 100, -3.0, 3.0}, "dE"); + auto hdE_rel = d2.Histo1D({"hdE_rel", "dE Relative; dE Relative; Events", 100, -3.0, 3.0}, "dE_rel"); + hdE->Fit("gaus", "", "", -3.0, 3.0); double* res = hdE->GetFunction("gaus")->GetParameters(); double sigmaOverE = res[2] / thrown_energy; //Pass/Fail - if (sigmaOverE <= resolutionTarget) { pi0_energy_resolution.pass(sigmaOverE); } else {