Resolve "Add energy scan benchmark for energy resolution"

Closes #36 (closed)

benchmarks/clustering/options/full_cal_reco.py

@@ -2,6 +2,7 @@
     An example option file to digitize/reconstruct/clustering calorimeter hits
 '''
 from Gaudi.Configuration import *
+import json
 import os
 import ROOT
 
@@ -12,13 +13,21 @@ detector_name = str(os.environ.get("JUGGLER_DETECTOR", "athena"))
 detector_path = str(os.environ.get("DETECTOR_PATH", "."))
 compact_path = os.path.join(detector_path, detector_name)
 
+# input arguments from calibration file
+with open('config/emcal_barrel_calibration.json') as f:
+    calib_data = json.load(f)['electron']
+
+print(calib_data)
+
+cb_ecal_sf = float(calib_data['sampling_fraction_img'])
+scifi_barrel_sf = float(calib_data['sampling_fraction_scfi'])
+
 # get sampling fractions from system environment variable, 1.0 by default
 ci_ecal_sf = float(os.environ.get("CI_ECAL_SAMP_FRAC", 0.253))
-cb_ecal_sf = float(os.environ.get("CB_ECAL_SAMP_FRAC", 0.01324))
 cb_hcal_sf = float(os.environ.get("CB_HCAL_SAMP_FRAC", 0.038))
 ci_hcal_sf = float(os.environ.get("CI_HCAL_SAMP_FRAC", 0.025))
 ce_hcal_sf = float(os.environ.get("CE_HCAL_SAMP_FRAC", 0.025))
-scifi_barrel_sf = float(os.environ.get("CB_EMCAL_SCFI_SAMP_FRAC",0.0938))
+
 # input and output
 input_sims = [f.strip() for f in str.split(os.environ["JUGGLER_SIM_FILE"], ",") if f.strip()]
 output_rec = str(os.environ["JUGGLER_REC_FILE"])