From 76b05365860521da321310a0b928c7e000455b76 Mon Sep 17 00:00:00 2001
From: Maria Zurek <zurek@anl.gov>
Date: Thu, 16 Sep 2021 22:25:08 +0000
Subject: [PATCH] Save the fsam mean error in the json file
---
.../scripts/emcal_barrel_particles_analysis.cxx | 11 ++++++++++-
1 file changed, 10 insertions(+), 1 deletion(-)
diff --git a/benchmarks/barrel_ecal/scripts/emcal_barrel_particles_analysis.cxx b/benchmarks/barrel_ecal/scripts/emcal_barrel_particles_analysis.cxx
index 61edb15d..c92f3402 100644
--- a/benchmarks/barrel_ecal/scripts/emcal_barrel_particles_analysis.cxx
+++ b/benchmarks/barrel_ecal/scripts/emcal_barrel_particles_analysis.cxx
@@ -55,6 +55,9 @@ void emcal_barrel_particles_analysis(std::string particle_name = "electron", boo
double fSam_mean;
double fSam_img_mean;
double fSam_scfi_mean;
+ double fSam_mean_err;
+ double fSam_img_mean_err;
+ double fSam_scfi_mean_err;
ROOT::EnableImplicitMT();
std::string input_fname = fmt::format("sim_output/sim_emcal_barrel_{}.root", particle_name);
@@ -181,6 +184,7 @@ void emcal_barrel_particles_analysis(std::string particle_name = "electron", boo
h->GetXaxis()->SetRangeUser(down_fit,up_fit);
TF1 *gaus = h->GetFunction("gaus");
fSam_mean = gaus->GetParameter(1);
+ fSam_mean_err = gaus->GetParError(1);
gaus->SetLineWidth(2);
gaus->SetLineColor(kRed);
save_canvas(c4,"fsam",particle_name);
@@ -197,6 +201,7 @@ void emcal_barrel_particles_analysis(std::string particle_name = "electron", boo
h->GetXaxis()->SetRangeUser(down_fit,up_fit);
TF1 *gaus = h->GetFunction("gaus");
fSam_img_mean = gaus->GetParameter(1);
+ fSam_img_mean_err = gaus->GetParError(1);
gaus->SetLineWidth(2);
gaus->SetLineColor(kRed);
save_canvas(c5,"fsamImg",particle_name);
@@ -213,6 +218,7 @@ void emcal_barrel_particles_analysis(std::string particle_name = "electron", boo
h->GetXaxis()->SetRangeUser(down_fit,up_fit);
TF1 *gaus = h->GetFunction("gaus");
fSam_scfi_mean = gaus->GetParameter(1);
+ fSam_scfi_mean_err = gaus->GetParError(1);
gaus->SetLineWidth(2);
gaus->SetLineColor(kRed);
save_canvas(c6,"fsamScFi",particle_name);
@@ -222,8 +228,11 @@ void emcal_barrel_particles_analysis(std::string particle_name = "electron", boo
{"particle_name", particle_name},
{"thrown_energy", Ethr_mean},
{"sampling_fraction", fSam_mean},
+ {"sampling_fraction_error", fSam_mean_err},
{"sampling_fraction_img", fSam_img_mean},
- {"sampling_fraction_scfi", fSam_scfi_mean}
+ {"sampling_fraction_error_img", fSam_img_mean_err},
+ {"sampling_fraction_scfi", fSam_scfi_mean},
+ {"sampling_fraction_error_scfi", fSam_scfi_mean_err}
};
if (save_calib) {
std::string calib_output_path = "results/emcal_barrel_calibration.json";
--
GitLab