From df2472ab8b3944ab44e9c25f8b9db1ac8f0b3ca0 Mon Sep 17 00:00:00 2001
From: Maria Zurek <zurek@anl.gov>
Date: Sun, 8 Aug 2021 17:29:34 +0000
Subject: [PATCH] Resolve "Add sampling fractions from BEMC to the calibration
 file."

---
 benchmarks/barrel_ecal/config.yml                         | 1 +
 .../scripts/emcal_barrel_particles_analysis.cxx           | 8 +++++++-
 2 files changed, 8 insertions(+), 1 deletion(-)

diff --git a/benchmarks/barrel_ecal/config.yml b/benchmarks/barrel_ecal/config.yml
index ff84a74d..dcf41f27 100644
--- a/benchmarks/barrel_ecal/config.yml
+++ b/benchmarks/barrel_ecal/config.yml
@@ -16,6 +16,7 @@ sim:emcal_barrel_electrons:
   stage: simulate
   script:
     - if [[ "$RUN_EXTENDED_BENCHMARK" == "true" ]] ; then bash benchmarks/barrel_ecal/run_emcal_barrel_energy_scan.sh electron ; fi
+    - export JUGGLER_N_EVENTS=400
     - bash benchmarks/barrel_ecal/run_emcal_barrel_particles.sh electron
 
 sim:emcal_barrel_photons:
diff --git a/benchmarks/barrel_ecal/scripts/emcal_barrel_particles_analysis.cxx b/benchmarks/barrel_ecal/scripts/emcal_barrel_particles_analysis.cxx
index e845ea7a..a56ececd 100644
--- a/benchmarks/barrel_ecal/scripts/emcal_barrel_particles_analysis.cxx
+++ b/benchmarks/barrel_ecal/scripts/emcal_barrel_particles_analysis.cxx
@@ -53,6 +53,8 @@ void emcal_barrel_particles_analysis(std::string particle_name = "electron", boo
   // variables that will be saved in the JSON file
   double Ethr_mean;
   double fSam_mean;
+  double fSam_img_mean;
+  double fSam_scfi_mean;
 
   ROOT::EnableImplicitMT();
   std::string input_fname = fmt::format("sim_output/sim_emcal_barrel_{}.root", particle_name);
@@ -194,6 +196,7 @@ void emcal_barrel_particles_analysis(std::string particle_name = "electron", boo
     h->Fit("gaus", "", "", down_fit, up_fit);
     h->GetXaxis()->SetRangeUser(down_fit,up_fit);
     TF1 *gaus = h->GetFunction("gaus");
+    fSam_img_mean = gaus->GetParameter(1);
     gaus->SetLineWidth(2);
     gaus->SetLineColor(kRed); 
     save_canvas(c5,"fsamImg",particle_name);
@@ -209,6 +212,7 @@ void emcal_barrel_particles_analysis(std::string particle_name = "electron", boo
     h->Fit("gaus", "", "", down_fit, up_fit);
     h->GetXaxis()->SetRangeUser(down_fit,up_fit);
     TF1 *gaus = h->GetFunction("gaus");
+    fSam_scfi_mean = gaus->GetParameter(1);
     gaus->SetLineWidth(2);
     gaus->SetLineColor(kRed); 
     save_canvas(c6,"fsamScFi",particle_name);
@@ -217,7 +221,9 @@ void emcal_barrel_particles_analysis(std::string particle_name = "electron", boo
   j[particle_name] = {
     {"particle_name", particle_name},
     {"thrown_energy", Ethr_mean},
-    {"sampling_fraction", fSam_mean}
+    {"sampling_fraction", fSam_mean},
+    {"sampling_fraction_img", fSam_img_mean},
+    {"sampling_fraction_scfi", fSam_scfi_mean}      
   };
   if (save_calib) {
     std::string calib_output_path = "results/emcal_barrel_calibration.json";
-- 
GitLab