diff --git a/benchmarks/clustering/options/full_cal_reco.py b/benchmarks/clustering/options/full_cal_reco.py
index 519bd49a232f04ab85cd43d82275d51c3be426bc..299a78c3b69f5e76c12c0319242e578141928e2b 100644
--- a/benchmarks/clustering/options/full_cal_reco.py
+++ b/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",
diff --git a/benchmarks/ecal/options/barrel.py b/benchmarks/ecal/options/barrel.py
index dfee80836d73d35a4309dd248a8e87a747e5f12c..6bfe72fc96091af1f816c93770b1008d3ef433fa 100644
--- a/benchmarks/ecal/options/barrel.py
+++ b/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)
 
diff --git a/benchmarks/imaging_ecal/options/hybrid_cluster.py b/benchmarks/imaging_ecal/options/hybrid_cluster.py
index 5c50009c948fde883b46c30c63eede6122389fa4..a457fe69891275964e13046a532f15472f9c2441 100644
--- a/benchmarks/imaging_ecal/options/hybrid_cluster.py
+++ b/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",