diff --git a/benchmarks/barrel_ecal/config.yml b/benchmarks/barrel_ecal/config.yml index ff84a74dcf09e06ebc7e21f58b3d4c957ae3c92e..dcf41f273946ad191130f3d4eb0440c9cfc9cc7e 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 e845ea7a0f36f4720f10cce59568acc5357cd0b0..a56ececde8c905c1fa9811d20195c50ade9aa0f7 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";