Resolve "Fix the cluster benchmark to use sampling fraction"

Closes #24 (closed)

Added sampling fraction as input in option at the last step of reconstruction EIC/juggler!111 (merged)

benchmarks/clustering/options/calorimeter_clustering.py

@@ -12,6 +12,11 @@ if "JUGGLER_DETECTOR" in os.environ :
 if "JUGGLER_DETECTOR_PATH" in os.environ :
   detector_name = str(os.environ["JUGGLER_DETECTOR_PATH"]) + "/" + detector_name
 
+# get sampling fraction from system environment variable, 1.0 by default
+sf = 1.0
+if "CB_EMCAL_SAMP_FRAC" in os.environ :
+  sf = str(os.environ["CB_EMCAL_SAMP_FRAC"])
+
 # todo add checks
 input_sim_file  = str(os.environ["JUGGLER_SIM_FILE"])
 output_rec_file = str(os.environ["JUGGLER_REC_FILE"])
@@ -70,6 +75,7 @@ ec_barrel_cluster = IslandCluster("ec_barrel_cluster",
         minClusterCenterEdep=1*units.MeV, 
         groupRange=2.0,
         OutputLevel=DEBUG)
+
 crystal_ec_cluster = IslandCluster("crystal_ec_cluster", 
         inputHitCollection="RecoEcalHits", 
         outputClusterCollection="EcalClusters",
@@ -82,8 +88,17 @@ simple_cluster = SimpleClustering("simple_cluster",
         outputClusters="SimpleClusters",
         OutputLevel=DEBUG)
 
+ec_barrel_clusterreco = RecoCoG("ec_barrel_clusterreco",
+        clusterCollection="EcalBarrelClusters", 
+        logWeightBase=6.2,
+        samplingFraction=sf) 
 
-clusterreco = RecoCoG("cluster_reco", clusterCollection="EcalClusters", logWeightBase=4.2, moduleDimZName="CrystalBox_z_length",OutputLevel=DEBUG)
+clusterreco = RecoCoG("cluster_reco", 
+        clusterCollection="EcalClusters", 
+        logWeightBase=4.2, 
+        moduleDimZName="CrystalBox_z_length",
+        samplingFraction=sf, 
+        OutputLevel=DEBUG)
 
 out = PodioOutput("out", filename=output_rec_file)
 
@@ -93,7 +108,7 @@ ApplicationMgr(
     TopAlg = [podioinput, copier, calcopier,
               ecdigi, emcaldigi, 
               crystal_ec_reco, ecal_reco, 
-              ec_barrel_cluster, crystal_ec_cluster, clusterreco,
+              ec_barrel_cluster, crystal_ec_cluster, ec_barrel_clusterreco, clusterreco,
               simple_cluster,
               out],
     EvtSel = 'NONE',