**/options: enable adjacencyMatrix in IslandCluster for SciGlass

benchmarks/clustering/options/full_cal_reco.py

@@ -258,7 +258,11 @@ else:
         splitCluster=False,
         minClusterHitEdep=1.0*MeV,  # discard low energy hits
         minClusterCenterEdep=30*MeV,
-        sectorDist=5.0*cm)
+        # Magic constants:
+        #  24 - number of sectors
+        #  5  - number of towers per sector
+        adjacencyMatrix = "(abs(tower_1 - tower_2) + (abs((sector_1 - sector_2) * 5 + row_1 - row_2) == 1) + (abs((sector_1 - sector_2) * 5 + row_1 - row_2) == (24 * 5 - 1))) == 1",
+        readoutClass = "EcalBarrelSciGlassHits")
     algs.append(cb_ecal_cl)
 
     cb_ecal_clreco = ImagingClusterReco("cb_ecal_clreco",

benchmarks/ecal/options/barrel.py

@@ -161,7 +161,11 @@ else:
         splitCluster=False,
         minClusterHitEdep=1.0 * MeV,  # discard low energy hits
         minClusterCenterEdep=30 * MeV,
-        sectorDist=5.0 * cm,
+        # Magic constants:
+        #  24 - number of sectors
+        #  5  - number of towers per sector
+        adjacencyMatrix = "(abs(tower_1 - tower_2) + (abs((sector_1 - sector_2) * 5 + row_1 - row_2) == 1) + (abs((sector_1 - sector_2) * 5 + row_1 - row_2) == (24 * 5 - 1))) == 1",
+        readoutClass = "EcalBarrelSciGlassHits",
     )
     algorithms.append(sciglass_ecal_cl)
 

benchmarks/imaging_ecal/options/hybrid_cluster.py

@@ -212,7 +212,11 @@ else:
         splitCluster=False,
         minClusterHitEdep=1.0*MeV,  # discard low energy hits
         minClusterCenterEdep=30*MeV,
-        sectorDist=5.0*cm)
+        # Magic constants:
+        #  24 - number of sectors
+        #  5  - number of towers per sector
+        adjacencyMatrix = "(abs(tower_1 - tower_2) + (abs((sector_1 - sector_2) * 5 + row_1 - row_2) == 1) + (abs((sector_1 - sector_2) * 5 + row_1 - row_2) == (24 * 5 - 1))) == 1",
+        readoutClass = "EcalBarrelSciGlassHits")
     algorithms.append(sciglass_ecal_cl)
 
     sciglass_ecal_clreco = ImagingClusterReco("sciglass_ecal_clreco",