From c88e7ac2fce2d6b48b9a1c0bb2b133a84ef2542d Mon Sep 17 00:00:00 2001
From: Marshall Scott <mbscott@anl.gov>
Date: Wed, 5 May 2021 20:53:45 -0400
Subject: [PATCH] Changed binning
---
.../barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx | 7 +++----
1 file changed, 3 insertions(+), 4 deletions(-)
diff --git a/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx b/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx
index 486755fc..52594492 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 {
--
GitLab