Use only two energies

benchmarks/imaging_ecal/config.yml

@@ -51,7 +51,7 @@ imaging_ecal_energy_scan_ph:
     - bash benchmarks/imaging_ecal/run_emcal_barrel.sh -t emcal_barrel_${PARTICLE}_${ENERGY} -n 100 -p "${PARTICLE}" -e "${ENERGY}" && echo "${ENERGY}" >> "$E_file"
   parallel:
     matrix:
-      - ENERGY: ["0.5", "1", "2", "5", "10"]
+      - ENERGY: ["0.5", "1"]
         PARTICLE: [photon]
 
 imaging_ecal_energy_scan_analysis_electrons:

benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis.cxx

@@ -225,8 +225,8 @@ std::tuple <double, double, double, double, double, double, double, double> extr
   hfsamImg->SetLineColor(kBlue);
   set_histo_range(hfsamImg);
   hfsamImg->DrawClone();
-  double up_fit = hfsamImg->GetMean() + 5*hfsamImg->GetStdDev();
-  double down_fit = hfsamImg->GetMean() - 5*hfsamImg->GetStdDev();
+  auto up_fit = hfsamImg->GetMean() + 5*hfsamImg->GetStdDev();
+  auto down_fit = hfsamImg->GetMean() - 5*hfsamImg->GetStdDev();
   if(down_fit <=0 ) down_fit = hfsamImg->GetXaxis()->GetBinUpEdge(1);
   hfsamImg->Fit("gaus", "", "", down_fit, up_fit);
   TF1 *gausImg = hfsamImg->GetFunction("gaus");
@@ -246,8 +246,8 @@ std::tuple <double, double, double, double, double, double, double, double> extr
   set_histo_range(hfsamScFi);
   hfsamScFi->DrawClone();
   hfsamScFi->Fit("gaus", "", "", 0.01, 0.1);
-  up_fit = hfsamScFi->GetMean() + 5*hfsamScFi->GetStdDev();
-  down_fit = hfsamScFi->GetMean() - 5*hfsamScFi->GetStdDev();
+  auto up_fit = hfsamScFi->GetMean() + 5*hfsamScFi->GetStdDev();
+  auto down_fit = hfsamScFi->GetMean() - 5*hfsamScFi->GetStdDev();
   if(down_fit <=0 ) down_fit = hfsamScFi->GetXaxis()->GetBinUpEdge(1);
   hfsamScFi->Fit("gaus", "", "", down_fit, up_fit);
   TF1 *gaus = hfsamScFi->GetFunction("gaus");