From 9c945ef73e9654370b1afe0cc3a1c79ec7d7c319 Mon Sep 17 00:00:00 2001
From: Maria <zurek@anl.gov>
Date: Mon, 20 Sep 2021 10:54:03 -0500
Subject: [PATCH] Calculate nb of thrown events
---
benchmarks/imaging_ecal/config.yml | 2 +-
.../scripts/emcal_barrel_energy_scan_analysis_pi_e.cxx | 6 +++---
2 files changed, 4 insertions(+), 4 deletions(-)
diff --git a/benchmarks/imaging_ecal/config.yml b/benchmarks/imaging_ecal/config.yml
index 4162525f..6de9dee5 100644
--- a/benchmarks/imaging_ecal/config.yml
+++ b/benchmarks/imaging_ecal/config.yml
@@ -37,7 +37,7 @@ imaging_ecal_energy_scan_e:
- bash benchmarks/imaging_ecal/run_emcal_barrel.sh -t emcal_barrel_${PARTICLE}_${ENERGY} -n 200 -p "${PARTICLE}" -e "${ENERGY}" && echo "${ENERGY}" >> "$E_file"
parallel:
matrix:
- - ENERGY: ["2", "5", "10"]
+ - ENERGY: ["2"]
PARTICLE: [electron, pion-]
imaging_ecal_energy_scan_ph:
diff --git a/benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis_pi_e.cxx b/benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis_pi_e.cxx
index 444a6615..08530a8d 100644
--- a/benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis_pi_e.cxx
+++ b/benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis_pi_e.cxx
@@ -191,18 +191,18 @@ std::tuple <double, double> extract_sampling_fraction_parameters(std::string E_l
auto nbins = helectron->GetXaxis()->GetNbins();
auto integral_helectron = helectron->Integral();
auto ngen_helectron = delectrons.Count();
- std::cout << infoHists[col].name << ": Thrown electron events: " << ngen_helectron << " Histo integral: " << integral_helectron << std::endl;
+ std::cout << infoHists[col].name << ": Thrown electron events: " << (*ngen_helectron) << " Histo integral: " << integral_helectron << std::endl;
auto efficiency_cut_bin = 0;
for(int bin=1; bin<=nbins;bin++){
auto inegral_helectron_bin = helectron->Integral(1,bin);
- auto efficiency = 1. - inegral_helectron_bin/ngen_helectron;
+ auto efficiency = 1. - inegral_helectron_bin/(*ngen_helectron);
if(efficiency > requested_efficiency) continue;
else{
efficiency_cut_bin = bin-1;
break;
}
}
- std::cout << "efficiency cut bin: " << efficiency_cut_bin << " efficiency: " << 1. - (helectron->Integral(1,efficiency_cut_bin))/ngen_helectron << std::endl;
+ std::cout << "efficiency cut bin: " << efficiency_cut_bin << " efficiency: " << 1. - (helectron->Integral(1,efficiency_cut_bin))/(*ngen_helectron) << std::endl;
auto ngen_hpion = dpions.Count();
auto pion_suppresion = hpion->Integral(efficiency_cut_bin, nbins)/ngen_hpion;
std::cout << "pion supression: " << pion_suppresion << std::end;
--
GitLab