Calculate nb of thrown events

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:

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;