for debug hcal_clustering.py with XY merger

benchmarks/clustering/barrel_clusters.sh

@@ -15,7 +15,7 @@ print_env.sh
 export JUGGLER_DETECTOR_PATH=${DETECTOR_PATH}
 
 if [[ ! -n  "${JUGGLER_N_EVENTS}" ]] ; then 
-  export JUGGLER_N_EVENTS=100
+  export JUGGLER_N_EVENTS=1
 fi
 
 if [[ ! -n "${CB_EMCAL_SAMP_FRAC}" ]] ; then
@@ -31,7 +31,7 @@ export JUGGLER_REC_FILE="rec_${JUGGLER_FILE_NAME_TAG}.root"
 echo "JUGGLER_N_EVENTS = ${JUGGLER_N_EVENTS}"
 echo "JUGGLER_DETECTOR = ${JUGGLER_DETECTOR}"
 
-root -b -q "benchmarks/clustering/scripts/gen_central_electrons.cxx(${JUGGLER_N_EVENTS}, \"${JUGGLER_FILE_NAME_TAG}.hepmc\")"
+root -b -q "benchmarks/clustering/scripts/gen_central_pions.cxx(${JUGGLER_N_EVENTS}, \"${JUGGLER_FILE_NAME_TAG}.hepmc\")"
 
 ### run geant4 simulations
 npsim --runType batch \

benchmarks/clustering/options/hcal_clustering.py

@@ -87,13 +87,16 @@ hcal_reco = HCalReconstruction("hcal_reco",
         OutputLevel=DEBUG)
 
 #clusters
-xymerger = CalorimeterHitsMerger("hcal_xy_merger",                             
+xymerger = CalorimeterHitsMerger("hcal_xy_merger",                            # 
           fields=["layer", "slice"],
+          fieldRefNumbers=[0,0],
           inputHitCollection="RecHcalBarrelHits",
-          outputHitCollection="RecHcalBarrelHitsXY")
+          outputHitCollection="RecHcalBarrelHitsXY",
+          OutputLevel=DEBUG)                       
+          
 
 hcal_barrel_cluster = IslandCluster("hcal_barrel_cluster", 
-        inputHitCollection="RecHcalBarrelHitsXY", 
+        inputHitCollection="RecHcalBarrelHits", 
         outputClusterCollection="HcalBarrelClusters",
         minClusterCenterEdep=30*units.MeV, 
         groupRange=2.0,

benchmarks/clustering/scripts/gen_central_pions.cxx

@@ -39,7 +39,7 @@ void gen_central_pions(int n_events = 100,
         FourVector(0.0, 0.0, 0.0, 0.938), 2212, 4);
 
     // Define momentum
-    Double_t p     = r1->Uniform(100.0, 100.1);
+    Double_t p     = r1->Uniform(20.0, 20.1);
     Double_t phi   = r1->Uniform(0.0, 0.25 * M_PI);
     Double_t costh = r1->Uniform(cos_theta_min, cos_theta_max);
     Double_t th    = std::acos(costh);