Merge pull request #49158 from rdelliga/occupancy_ETL

ETL occupancy study

Author: cmsbuild

Geometry/MTDGeometryBuilder/test/MTDDigiGeometryAnalyzer.cc

@@ -56,10 +56,24 @@ class MTDDigiGeometryAnalyzer : public edm::one::EDAnalyzer<> {
 
   void checkRectangularMTDTopology(const RectangularMTDTopology&);
   void checkPixelsAcceptance(const GeomDetUnit& det);
+  void CheckETLstructure(const MTDGeometry&);
 
   std::stringstream sunitt_;
 
   edm::ESGetToken<MTDGeometry, MTDDigiGeometryRecord> mtdgeoToken_;
+
+  // Constants to define the bins for Eta
+  static constexpr int n_bin_Eta = 3;
+  static constexpr double eta_bins_edges_neg[n_bin_Eta + 1] = {-3.0, -2.5, -2.1, -1.5};
+  static constexpr double eta_bins_edges_pos[n_bin_Eta + 1] = {1.5, 2.1, 2.5, 3.0};
+
+  // LGAD counter per Disc, DiscSide, and Sector: [disk][discSide][sector]
+  static constexpr int n_discSide = 2;
+  static constexpr int n_sector = 3;  // Use size 3 to allow 1-based indexing (1 to 2)
+  uint32_t LGADsPerDiscSideSector_[4][n_discSide][n_sector] = {};
+
+  // Counter for total LGADs per disk per eta bin: [disk][eta_bin]
+  uint32_t LGADsPerDiskperEtaBin_[4][n_bin_Eta] = {{0}};
 };
 
 MTDDigiGeometryAnalyzer::MTDDigiGeometryAnalyzer(const edm::ParameterSet& iConfig) {
@@ -128,6 +142,10 @@ void MTDDigiGeometryAnalyzer::analyze(const edm::Event& iEvent, const edm::Event
   auto const& etldet = *(dynamic_cast<const MTDGeomDetUnit*>(pDD->detsETL().front()));
   checkPixelsAcceptance(etldet);
 
+  // ETL structure prints
+  edm::LogVerbatim("MTDDigiGeometryAnalyzer") << "\n";
+  CheckETLstructure(*pDD);
+
   edm::LogVerbatim("MTDUnitTest") << sunitt_.str();
 }
 
@@ -178,5 +196,106 @@ void MTDDigiGeometryAnalyzer::checkPixelsAcceptance(const GeomDetUnit& det) {
   sunitt_ << " Acceptance: " << fround(acc, 3) << " +/- " << fround(accerr, 3);
 }
 
+void MTDDigiGeometryAnalyzer::CheckETLstructure(const MTDGeometry& geom) {
+  edm::LogVerbatim("MTDDigiGeometryAnalyzer") << "\n--- ETL Structure Validation ---";
+
+  // Reset counters
+  for (int d = 0; d < 4; ++d) {
+    for (int eta = 0; eta < n_bin_Eta; ++eta) {
+      LGADsPerDiskperEtaBin_[d][eta] = 0;
+    }
+    for (int k = 0; k < n_discSide; ++k) {
+      for (int l = 0; l < n_sector; ++l) {
+        LGADsPerDiscSideSector_[d][k][l] = 0;
+      }
+    }
+  }
+
+  uint32_t totalETLdets = 0;
+  for (const auto& det : geom.detsETL()) {
+    const GeomDet* thedet = det;
+    ETLDetId detId(thedet->geographicalId());
+
+    // Get the global position of the detector center
+    const GlobalPoint& global_point = thedet->position();
+    double eta = global_point.eta();
+
+    int discSide = detId.discSide();  // 0 to 1
+    int sector = detId.sector();      // 1 to 2
+
+    int idet = 999;
+    if ((detId.zside() == -1) && (detId.nDisc() == 1)) {
+      idet = 0;
+    } else if ((detId.zside() == -1) && (detId.nDisc() == 2)) {
+      idet = 1;
+    } else if ((detId.zside() == 1) && (detId.nDisc() == 1)) {
+      idet = 2;
+    } else if ((detId.zside() == 1) && (detId.nDisc() == 2)) {
+      idet = 3;
+    } else {
+      edm::LogWarning("EtlDigiHitsValidation") << "Unknown ETL DetId configuration: " << detId;
+      continue;
+    }
+
+    totalETLdets++;
+
+    // Count LGADs per Disc, Side, Sector
+    LGADsPerDiscSideSector_[idet][discSide][sector]++;
+
+    // Count LGADs per disk per eta bin
+    const double* eta_edges = (idet < 2) ? eta_bins_edges_neg : eta_bins_edges_pos;
+
+    for (int j = 0; j < n_bin_Eta; j++) {
+      double lower_edge = eta_edges[j];
+      double upper_edge = eta_edges[j + 1];
+
+      // Check if the center of the LGAD is within the bin
+      if ((eta >= lower_edge && eta < upper_edge) || (idet < 2 && j == n_bin_Eta - 1 && eta <= upper_edge)) {
+        LGADsPerDiskperEtaBin_[idet][j]++;
+        break;  // Found the bin
+      }
+    }
+  }
+
+  // --- Print Summary ---
+
+  edm::LogVerbatim("MTDDigiGeometryAnalyzer") << " Total ETL Detectors (LGADs): " << totalETLdets << "\n";
+  const char* diskNames[4] = {"Disc 1 (-Z)", "Disc 2 (-Z)", "Disc 1 (+Z)", "Disc 2 (+Z)"};
+
+  edm::LogVerbatim("MTDDigiGeometryAnalyzer") << "\n--- LGADs per Eta Bin and per Disk, DiscSide, Sector ---\n";
+  for (int d = 0; d < 4; ++d) {  // Physical Disk loop (0-3)
+    std::string disk_name = diskNames[d];
+    uint32_t total_disk = 0;
+    for (int k = 0; k < n_discSide; ++k) {
+      for (int l = 1; l < n_sector; ++l) {
+        total_disk += LGADsPerDiscSideSector_[d][k][l];
+      }
+    }
+    edm::LogVerbatim("MTDDigiGeometryAnalyzer") << "Region: " << disk_name << " | Total LGADs: " << total_disk << "\n";
+    edm::LogVerbatim("MTDDigiGeometryAnalyzer") << "  - LGADs per Eta Bin:\n";
+    const double* eta_edges = (d < 2) ? eta_bins_edges_neg : eta_bins_edges_pos;
+    for (int j = 0; j < n_bin_Eta; ++j) {
+      edm::LogVerbatim("MTDDigiGeometryAnalyzer")
+          << "    Eta [" << std::setprecision(1) << std::fixed << eta_edges[j] << ", " << eta_edges[j + 1]
+          << "): " << LGADsPerDiskperEtaBin_[d][j] << "\n";
+    }
+    for (int k = 0; k < n_discSide; ++k) {
+      uint32_t total_discside = 0;
+      for (int l = 1; l < n_sector; ++l) {
+        total_discside += LGADsPerDiscSideSector_[d][k][l];
+      }
+      edm::LogVerbatim("MTDDigiGeometryAnalyzer") << "  - Side: " << k << " | Total LGADs: " << total_discside << "\n";
+      edm::LogVerbatim("MTDDigiGeometryAnalyzer") << "    - Sectors: ";
+      for (int l = 1; l < n_sector; ++l) {
+        if (LGADsPerDiscSideSector_[d][k][l] > 0) {
+          edm::LogVerbatim("MTDDigiGeometryAnalyzer")
+              << "Sec " << l << ": " << LGADsPerDiscSideSector_[d][k][l] << " | ";
+        }
+      }
+      edm::LogVerbatim("MTDDigiGeometryAnalyzer") << "\n";
+    }
+  }
+}
+
 //define this as a plug-in
 DEFINE_FWK_MODULE(MTDDigiGeometryAnalyzer);