From 141d9788de5607879da718ecc7149453454ed7ed Mon Sep 17 00:00:00 2001
From: Wouter Deconinck <wdconinc@gmail.com>
Date: Thu, 22 Dec 2022 00:22:59 +0000
Subject: [PATCH] fix: emcal_barrel_pi0 fallback to stats if fit fails
---
.../barrel_ecal/scripts/emcal_barrel_pi0_analysis.cxx | 10 +++++++---
1 file changed, 7 insertions(+), 3 deletions(-)
diff --git a/benchmarks/barrel_ecal/scripts/emcal_barrel_pi0_analysis.cxx b/benchmarks/barrel_ecal/scripts/emcal_barrel_pi0_analysis.cxx
index dd917e3d..6d4a51ef 100644
--- a/benchmarks/barrel_ecal/scripts/emcal_barrel_pi0_analysis.cxx
+++ b/benchmarks/barrel_ecal/scripts/emcal_barrel_pi0_analysis.cxx
@@ -190,9 +190,13 @@ void emcal_barrel_pi0_analysis(
h3->SetLineWidth(2);
h3->SetLineColor(kBlue);
auto fit = h3->Fit("gaus","","", fitRange[1][0], fitRange[1][1]);
- double* res = h3->GetFunction("gaus")->GetParameters();
- sigmaOverE = res[2] / meanE;
-
+ if (fit == 0) {
+ double* res = h3->GetFunction("gaus")->GetParameters();
+ sigmaOverE = res[2] / meanE;
+ } else {
+ std::printf("Fit failed\n");
+ sigmaOverE = h3->GetStdDev() / h3->GetMean();
+ }
c1->SaveAs((fmt::format("results/emcal_barrel_pi0_{}.png", col[1])).c_str());
c1->SaveAs((fmt::format("results/emcal_barrel_pi0_{}.pdf", col[1])).c_str());
std::printf("Resolution %d\n", 1);
--
GitLab