diff --git a/benchmarks/imaging_ecal/config.yml b/benchmarks/imaging_ecal/config.yml index 4162525fa700de08af94b664a3cf7d0cecf62deb..6de9dee5e8a807007695fb74d4b6cccd599b4b90 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 444a66151cda1d5c3980c75f0965b7b5c4bfd1e9..08530a8d6eb14fb164d33be52767bf96172ddac2 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;