FPGAOutputValidationAlg.cxx

Go to the documentation of this file.

/*
 *   Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
 */
 
#include "EFTrackingFPGAOutputValidation/FPGAOutputValidationAlg.h"
#include "InDetMeasurementUtilities/Helpers.h"
#include "xAODInDetMeasurement/Utilities.h"
 
namespace {
 
  template <typename T>
  struct ClusterPair {
    std::pair<const T*, const T*> clusters;
    size_t commonRDOs = 0;
 
    ClusterPair(const T* cluster1, const T* cluster2, size_t commonRDOsCount)
      : clusters(cluster1, cluster2), commonRDOs(commonRDOsCount) {};
  };

  template <typename T>
  size_t compareClusters(const T* cluster1, const T* cluster2, size_t* nonCommonRdo1 = nullptr, size_t* nonCommonRdo2 = nullptr) {
    const auto& rdoList1 = cluster1->rdoList();
    const auto& rdoList2 = cluster2->rdoList();
    std::unordered_set<Identifier> rdoSet1(rdoList1.begin(), rdoList1.end());
    std::unordered_set<Identifier> rdoSet2(rdoList2.begin(), rdoList2.end());
 
    size_t nCommonRdo = 0;
    for (const auto& rdo : rdoList1) {
      if (rdoSet2.count(rdo)) {
        ++nCommonRdo;
      }
    }
    if (nonCommonRdo1) *nonCommonRdo1 = rdoList1.size() - nCommonRdo;
    if (nonCommonRdo2) *nonCommonRdo2 = rdoList2.size() - nCommonRdo;
 
    return nCommonRdo;
  }

  template <typename T>
  std::vector<const T*> findMatchingCluster(const T* cluster0, 
                                            const std::unordered_multimap<xAOD::DetectorIdentType, const T*>& clusterMap,
                                            const std::unordered_multimap<xAOD::DetectorIDHashType, const T*>& clusterMapHashIdMap,
                                            bool matchByID, const size_t& allowedMisses) {
    std::vector<const T*> matchedClusters{};
 
    std::unordered_set<Identifier> rdoSet0{};
    for (const auto& rdo : cluster0->rdoList()) {
      rdoSet0.insert(rdo);
    }
    const size_t rdoToMiss = (allowedMisses < rdoSet0.size()) ? allowedMisses : (rdoSet0.size()-1);
    if (matchByID) {
      auto range = clusterMap.equal_range(cluster0->identifier());
      for (auto it = range.first; it != range.second; ++it) {
      size_t nCommonRdo = compareClusters(cluster0, it->second);
      if (nCommonRdo >= rdoSet0.size() - rdoToMiss) {
        matchedClusters.push_back(it->second);
      }
      }
    } else {
      auto range = clusterMapHashIdMap.equal_range(cluster0->identifierHash());
      for (auto it = range.first; it != range.second; ++it) {
      const auto& cluster1 = it->second;
      size_t nCommonRdo = compareClusters(cluster0, cluster1);
      if (nCommonRdo >= rdoSet0.size() - rdoToMiss) {
        matchedClusters.push_back(cluster1);
      }
      }
    }
    return matchedClusters;
  }
 

  template <typename T>
  std::vector<ClusterPair<T>> findNonMergedClusters(const std::unordered_multimap<xAOD::DetectorIDHashType, const T*>& clusterMap) {
    std::vector<ClusterPair<T>> clusterPairs;
    std::unordered_set<xAOD::DetectorIDHashType> uniqueKeys;
    for (const auto& pair : clusterMap) {
      uniqueKeys.insert(pair.first);
    }
    for (const auto& key : uniqueKeys) {
      auto range = clusterMap.equal_range(key);
      for (auto it1 = range.first; it1 != range.second; ++it1) {
        for (auto it2 = std::next(it1); it2 != range.second; ++it2) {
          size_t nonCommonRdo1 = 0, nonCommonRdo2 = 0;
          size_t nCommonRdo = compareClusters(it1->second, it2->second, &nonCommonRdo1, &nonCommonRdo2);
          if (nCommonRdo > 0) {
            clusterPairs.emplace_back(it1->second, it2->second, nCommonRdo);
          }
        }
      }
    }
    return clusterPairs;
  }
 
  template <typename T>
  std::string logMultipleClusterMatches(const std::vector<const T*>& matchedClusters) {
    std::string out;
 
    out += "Found " + std::to_string(matchedClusters.size()) + " FPGA ";
    out += (matchedClusters[0]->type() == xAOD::UncalibMeasType::PixelClusterType ? "Pixel" : "Strip");
    out += " clusters matching to the same offline cluster:\n";
    for (const auto& cluster : matchedClusters) {
      out += std::to_string(cluster->identifier()) + " x: " + std::to_string(cluster->globalPosition()[0]) +
             " y: " + std::to_string(cluster->globalPosition()[1]) + " z: " + std::to_string(cluster->globalPosition()[2]) + "\n";
      for (const auto& rdo : cluster->rdoList()) {
        out += "\t" + std::to_string(rdo.get_compact()) + "\n";
      }
      out += "\n";
    }
    return out;
  }
}
 
FPGAOutputValidationAlg::FPGAOutputValidationAlg(
  const std::string& name,
  ISvcLocator* pSvcLocator
) : AthReentrantAlgorithm(name, pSvcLocator)
{}
 
StatusCode FPGAOutputValidationAlg::initialize() {
  ATH_MSG_INFO("Initializing FPGAOutputValidationAlg");
 
  ATH_CHECK(m_pixelKeys.initialize(!m_pixelKeys.empty()));
  ATH_CHECK(m_stripKeys.initialize(!m_stripKeys.empty()));
 
  ATH_CHECK(m_monitoringTool.retrieve());
  
  ATH_CHECK(m_chrono.retrieve());
 
  ATH_CHECK(detStore()->retrieve(m_pixelid, "PixelID"));
  ATH_CHECK(detStore()->retrieve(m_stripid, "SCT_ID"));
  ATH_CHECK(detStore()->retrieve(m_SCT_mgr, "ITkStrip"));
  ATH_CHECK(detStore()->retrieve(m_PIX_mgr, "ITkPixel"));
 
  return StatusCode::SUCCESS;
}
 
StatusCode FPGAOutputValidationAlg::execute(const EventContext& ctx) const { 
  if (m_pixelKeys.size() == 2 && m_doDiffHistograms) { 
    m_chrono->chronoStart("FPGAOutputValidationAlg::pixel diff");
    const SG::ReadHandleKey<xAOD::PixelClusterContainer>& key0 = m_pixelKeys[0];
    const SG::ReadHandleKey<xAOD::PixelClusterContainer>& key1 = m_pixelKeys[1];
    SG::ReadHandle<xAOD::PixelClusterContainer> handle0{key0, ctx};
    ATH_CHECK(handle0.isValid());
    SG::ReadHandle<xAOD::PixelClusterContainer> handle1{key1, ctx};
    ATH_CHECK(handle1.isValid());
    
    if (m_checkClusterRdos) {
      const auto& pixelClusters0 = *handle0;
      std::unordered_multimap<xAOD::DetectorIDHashType, const xAOD::PixelCluster*> pixelClustersHashIdMap0; // assumes that the first key is the FPGA one
      for (const auto* cluster0 : pixelClusters0) {
        const xAOD::DetectorIDHashType hashId0 = cluster0->identifierHash();
        pixelClustersHashIdMap0.insert(std::make_pair(hashId0, cluster0));
      }
      const std::vector<ClusterPair<xAOD::PixelCluster>> pixelClusterPairsWithCommonRdos = findNonMergedClusters(pixelClustersHashIdMap0);
      if (pixelClusterPairsWithCommonRdos.size() > 0) {
        std::stringstream ss;
        for (const auto& pair : pixelClusterPairsWithCommonRdos) {
          ss << "Found " << pair.commonRDOs << " common RDOs between clusters with hash "
            << pair.clusters.first->identifierHash() << ": "
            << pair.clusters.first->identifier() << " and "
            << pair.clusters.second->identifier() << "\n";
        }
        ATH_MSG_ERROR("Pixel cluster pairs with common RDOs:\n" << ss.str());
      }
    }
 
    const xAOD::PixelClusterContainer pixelClusters1 = *handle1;
    std::unordered_multimap<xAOD::DetectorIdentType, const xAOD::PixelCluster*> pixelClustersMap1;
    std::unordered_multimap<xAOD::DetectorIDHashType, const xAOD::PixelCluster*> pixelClustersHashIdMap1;
    for (const auto* cluster1 : pixelClusters1) {
      const xAOD::DetectorIdentType id1 = cluster1->identifier();
      const xAOD::DetectorIDHashType hashId1 = cluster1->identifierHash();
      pixelClustersMap1.insert(std::make_pair(id1, cluster1));
      pixelClustersHashIdMap1.insert(std::make_pair(hashId1, cluster1));
    }
 
    for (auto cluster0 : *handle0) {
      const std::vector<const xAOD::PixelCluster*> matchedClusters = findMatchingCluster(cluster0, pixelClustersMap1, pixelClustersHashIdMap1,
                                                                                         m_matchByID,
                                                                                         m_allowedRdoMisses);
 
      if (matchedClusters.size() == 0) {
        std::vector<std::string> regions {"all"};
        if(m_pixelid->barrel_ec(cluster0->rdoList()[0]) == 0) regions.push_back("barrel");
        else regions.push_back("endcap");
        for (const auto& region : regions) {
          Monitored::Group(
            m_monitoringTool,
            Monitored::Scalar<float>(handle0.key() + "_UNMATCHED_GLOBALPOSITION_Z_" + region, cluster0->globalPosition()[2]),
            Monitored::Scalar<float>(handle0.key() + "_UNMATCHED_GLOBALPOSITION_R_" + region, sqrt(cluster0->globalPosition()[0]*cluster0->globalPosition()[0] + 
                                                                                                   cluster0->globalPosition()[1]*cluster0->globalPosition()[1])),
            Monitored::Scalar<float>("nmatched_pixel_clusters_"+region, matchedClusters.size() - 0.5)
          );
        }
        continue;
      }
      if (matchedClusters.size() > 1) {
        ATH_MSG_ERROR(logMultipleClusterMatches(matchedClusters));
        return StatusCode::FAILURE;
      }
 
      const xAOD::PixelCluster *cluster1 = matchedClusters[0];
      
      std::vector<std::string> regions {"all"};
      if(m_pixelid->barrel_ec(cluster0->rdoList()[0]) == 0) regions.push_back("barrel");
      else regions.push_back("endcap");
 
      for(auto const& region: regions)
      {
        Monitored::Group(
          m_monitoringTool,
          Monitored::Scalar<float>("diff_pixel_locx_" +region , cluster0->localPosition<2>()[0] - cluster1->localPosition<2>()[0]),
          Monitored::Scalar<float>("diff_pixel_locy_" +region , cluster0->localPosition<2>()[1] - cluster1->localPosition<2>()[1]),
          Monitored::Scalar<float>("diff_pixel_covxx_" +region , cluster0->localCovariance<2>()(0, 0) - cluster1->localCovariance<2>()(0, 0)),
          Monitored::Scalar<float>("diff_pixel_covyy_" +region , cluster0->localCovariance<2>()(1, 1) - cluster1->localCovariance<2>()(1, 1)),
          Monitored::Scalar<float>("diff_pixel_globalx_" +region , cluster0->globalPosition()[0] - cluster1->globalPosition()[0]),
          Monitored::Scalar<float>("diff_pixel_globaly_" +region , cluster0->globalPosition()[1] - cluster1->globalPosition()[1]),
          Monitored::Scalar<float>("diff_pixel_globalz_" +region , cluster0->globalPosition()[2] - cluster1->globalPosition()[2]),
          Monitored::Scalar<int>("diff_pixel_channelsphi_" +region , cluster0->channelsInPhi() - cluster1->channelsInPhi()),
          Monitored::Scalar<int>("diff_pixel_channelseta_" +region , cluster0->channelsInEta() - cluster1->channelsInEta()),
          Monitored::Scalar<float>("diff_pixel_widtheta_" +region , cluster0->widthInEta() - cluster1->widthInEta()),
          Monitored::Scalar<int>("diff_pixel_tot_" +region , xAOD::xAODInDetMeasurement::Utilities::computeTotalToT(*cluster0) - xAOD::xAODInDetMeasurement::Utilities::computeTotalToT(*cluster1)),
          Monitored::Scalar<int>("diff_pixel_rdos_" +region , cluster0->rdoList().size() - cluster1->rdoList().size()),
          Monitored::Scalar<float>("pixel_globalR_ref_" + region, sqrt(cluster1->globalPosition()[0]*cluster1->globalPosition()[0] + 
                                                                       cluster1->globalPosition()[1]*cluster1->globalPosition()[1])),
          Monitored::Scalar<float>("pixel_globalZ_ref_" + region, cluster1->globalPosition()[2]),
          Monitored::Scalar<float>("nmatched_pixel_clusters_"+region, matchedClusters.size() - 0.5)
        );
        if(region != "all")
        {
          Monitored::Group(
            m_monitoringTool,
              Monitored::Scalar<float>("diff_pixel_locx_"+region+"Layer" + std::to_string(m_pixelid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->localPosition<2>()[0] - cluster1->localPosition<2>()[0]),
              Monitored::Scalar<float>("diff_pixel_locy_"+region+"Layer" + std::to_string(m_pixelid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->localPosition<2>()[1] - cluster1->localPosition<2>()[1]),
              Monitored::Scalar<float>("diff_pixel_covxx_"+region+"Layer" + std::to_string(m_pixelid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->localCovariance<2>()(0, 0) - cluster1->localCovariance<2>()(0, 0)),
              Monitored::Scalar<float>("diff_pixel_covyy_"+region+"Layer" + std::to_string(m_pixelid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->localCovariance<2>()(1, 1) - cluster1->localCovariance<2>()(1, 1)),
              Monitored::Scalar<float>("diff_pixel_globalX_"+region+"Layer" + std::to_string(m_pixelid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->globalPosition()[0] - cluster1->globalPosition()[0]),
              Monitored::Scalar<float>("diff_pixel_globalY_"+region+"Layer" + std::to_string(m_pixelid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->globalPosition()[1] - cluster1->globalPosition()[1]),
              Monitored::Scalar<float>("diff_pixel_globalZ_"+region+"Layer" + std::to_string(m_pixelid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->globalPosition()[2] - cluster1->globalPosition()[2]),
              Monitored::Scalar<float>("nmatched_pixel_clusters_"+region, matchedClusters.size() - 0.5)
            );
        }
      }
    }
    m_chrono->chronoStop("FPGAOutputValidationAlg::pixel diff");
  }
 
  if (m_stripKeys.size() == 2 && m_doDiffHistograms) { 
    m_chrono->chronoStart("FPGAOutputValidationAlg::strip diff");
    const SG::ReadHandleKey<xAOD::StripClusterContainer>& key0 = m_stripKeys[0];
    const SG::ReadHandleKey<xAOD::StripClusterContainer>& key1 = m_stripKeys[1];
    SG::ReadHandle<xAOD::StripClusterContainer> handle0{key0, ctx};
    ATH_CHECK(handle0.isValid());
    SG::ReadHandle<xAOD::StripClusterContainer> handle1{key1, ctx};
    ATH_CHECK(handle1.isValid());
 
    if (m_checkClusterRdos) {
      const xAOD::StripClusterContainer stripClusters0 = *handle0;
      std::unordered_multimap<xAOD::DetectorIDHashType, const xAOD::StripCluster*> stripClustersHashIdMap0; // assumes that the first key is the FPGA one
      for (const auto* cluster0 : stripClusters0) {
        const xAOD::DetectorIDHashType hashId0 = cluster0->identifierHash();
        stripClustersHashIdMap0.insert(std::make_pair(hashId0, cluster0));
      }
 
      // Find all pairs of clusters within the same detector element that share at least one RDO
      const std::vector<ClusterPair<xAOD::StripCluster>> stripPairsWithCommonRdos = findNonMergedClusters(stripClustersHashIdMap0);
 
      // If any such pairs are found, log an error with details for debugging
      if (stripPairsWithCommonRdos.size() > 0) {
        std::stringstream ss;
        for (const auto& pair : stripPairsWithCommonRdos) {
          ss << "Found " << pair.commonRDOs << " common RDOs between clusters with hash "
            << pair.clusters.first->identifierHash() << ": "
            << pair.clusters.first->identifier() << " and "
            << pair.clusters.second->identifier() << "\n";
        }
        ATH_MSG_ERROR("Strip cluster pairs with common RDOs:\n" << ss.str());
      }
    }
    const xAOD::StripClusterContainer stripClusters1 = *handle1; 
    std::unordered_multimap<xAOD::DetectorIdentType, const xAOD::StripCluster*> stripClustersMap1;
    std::unordered_multimap<xAOD::DetectorIDHashType, const xAOD::StripCluster*> stripClustersHashIdMap1;
    for (const auto *cluster1 : stripClusters1) {
      const xAOD::DetectorIdentType id1 = cluster1->identifier();
      const xAOD::DetectorIDHashType hashId1 = cluster1->identifierHash();
      stripClustersMap1.insert(std::make_pair(id1, cluster1));
      stripClustersHashIdMap1.insert(std::make_pair(hashId1, cluster1));
    }
 
    
    for (auto cluster0 : *handle0) {
      const std::vector<const xAOD::StripCluster*> matchedClusters = findMatchingCluster(cluster0, stripClustersMap1, stripClustersHashIdMap1,
                                                                                         m_matchByID,
                                                                                         m_allowedRdoMisses);
 
      if (matchedClusters.size() == 0 && handle1->size() > 0) {
        std::vector<std::string> regions {"all"};
        if(m_stripid->barrel_ec(cluster0->rdoList()[0]) == 0) regions.push_back("barrel");
        else regions.push_back("endcap");
        for (const auto& region : regions) {
          Monitored::Group(
            m_monitoringTool,
            Monitored::Scalar<float>(handle0.key() + "_UNMATCHED_GLOBALPOSITION_Z_" + region, cluster0->globalPosition()[2]),
            Monitored::Scalar<float>(handle0.key() + "_UNMATCHED_GLOBALPOSITION_R_" + region, sqrt(cluster0->globalPosition()[0]*cluster0->globalPosition()[0] + 
                                                                                                   cluster0->globalPosition()[1]*cluster0->globalPosition()[1])),
            Monitored::Scalar<float>("nmatched_strip_clusters_"+region, matchedClusters.size() - 0.5)
          );
        }
        continue;
      }
      if (matchedClusters.size() > 1) {
        ATH_MSG_ERROR(logMultipleClusterMatches(matchedClusters));
        return StatusCode::FAILURE;
      }
 
      const xAOD::StripCluster *cluster1 = matchedClusters[0];
 
      std::vector<std::string> regions {"all"};
      if(m_stripid->barrel_ec(cluster0->rdoList()[0]) == 0) regions.push_back("barrel");
      else regions.push_back("endcap");
 
      for(auto const& region: regions)
      {
        Monitored::Group(
          m_monitoringTool,
          Monitored::Scalar<float>("diff_strip_locx_" + region, cluster0->localPosition<1>()[0] - cluster1->localPosition<1>()[0]),
          Monitored::Scalar<float>("diff_strip_locxZoom_" + region, cluster0->localPosition<1>()[0] - cluster1->localPosition<1>()[0]),
          Monitored::Scalar<float>("diff_strip_covxx_" + region, cluster0->localCovariance<1>()(0, 0) - cluster1->localCovariance<1>()(0, 0)),
          Monitored::Scalar<float>("diff_strip_globalx_" + region, cluster0->globalPosition()[0] - cluster1->globalPosition()[0]),
          Monitored::Scalar<float>("diff_strip_globaly_" + region, cluster0->globalPosition()[1] - cluster1->globalPosition()[1]),
          Monitored::Scalar<float>("diff_strip_globalz_" + region, cluster0->globalPosition()[2] - cluster1->globalPosition()[2]),
          Monitored::Scalar<float>("diff_strip_channelsphi_" + region, cluster0->channelsInPhi() - cluster1->channelsInPhi()),
          Monitored::Scalar<int>("diff_strip_rdos_" +region , cluster0->rdoList().size() - cluster1->rdoList().size()),
          Monitored::Scalar<float>("strip_globalR_ref_" + region, sqrt(cluster1->globalPosition()[0]*cluster1->globalPosition()[0] + 
                                                                       cluster1->globalPosition()[1]*cluster1->globalPosition()[1])),
          Monitored::Scalar<float>("strip_globalZ_ref_" + region, cluster1->globalPosition()[2]),
          Monitored::Scalar<float>("nmatched_strip_clusters_"+region, matchedClusters.size() - 0.5)
        );
        if(region != "all")
        {
          Monitored::Group(
            m_monitoringTool,
              Monitored::Scalar<float>("diff_strip_locx_"+region+"Layer" + std::to_string(m_stripid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->localPosition<1>()[0] - cluster1->localPosition<1>()[0]),
              Monitored::Scalar<float>("diff_strip_locxZoom_"+region+"Layer" + std::to_string(m_stripid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->localPosition<1>()[0] - cluster1->localPosition<1>()[0]),
              Monitored::Scalar<float>("diff_strip_covxx_"+region+"Layer" + std::to_string(m_stripid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->localCovariance<1>()(0, 0) - cluster1->localCovariance<1>()(0, 0)),
              Monitored::Scalar<float>("diff_strip_globalX_"+region+"Layer" + std::to_string(m_stripid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->globalPosition()[0] - cluster1->globalPosition()[0]),
              Monitored::Scalar<float>("diff_strip_globalY_"+region+"Layer" + std::to_string(m_stripid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->globalPosition()[1] - cluster1->globalPosition()[1]),
              Monitored::Scalar<float>("diff_strip_globalZ_"+region+"Layer" + std::to_string(m_stripid->layer_disk(m_pixelid->wafer_id(cluster0->identifierHash()))),
              cluster0->globalPosition()[2] - cluster1->globalPosition()[2]),
              Monitored::Scalar<float>("nmatched_strip_clusters_"+region, matchedClusters.size() - 0.5)
            );
        }
      }
    }
    m_chrono->chronoStop("FPGAOutputValidationAlg::strip diff");
  }
 
  for (std::size_t index = 0; index < m_pixelKeys.size(); index++) {
    const SG::ReadHandleKey<xAOD::PixelClusterContainer>& key = m_pixelKeys[index];
    SG::ReadHandle<xAOD::PixelClusterContainer> handle{key, ctx};
    ATH_CHECK(handle.isValid());
 
 
    for(auto cluster : *handle)
    {
      std::vector<std::string> regions {"all"};
      if(m_pixelid->barrel_ec(cluster->rdoList()[0]) == 0) regions.push_back("barrel");
      else regions.push_back("endcap");
 
      for(auto const& region: regions)
      {
        Monitored::Group(
          m_monitoringTool,
          Monitored::Scalar<float>(key.key() + "_LOCALPOSITION_X_" + region, cluster->localPosition<2>()[0]),
          Monitored::Scalar<float>(key.key() + "_LOCALPOSITION_Y_" + region, cluster->localPosition<2>()[1]),
          Monitored::Scalar<float>(key.key() + "_LOCALCOVARIANCE_XX_" + region, cluster->localCovariance<2>()(0, 0)),
          Monitored::Scalar<float>(key.key() + "_LOCALCOVARIANCE_YY_" + region, cluster->localCovariance<2>()(1, 1)),
          Monitored::Scalar<float>(key.key() + "_GLOBALPOSITION_X_" + region, cluster->globalPosition()[0]),
          Monitored::Scalar<float>(key.key() + "_GLOBALPOSITION_Y_" + region, cluster->globalPosition()[1]),
          Monitored::Scalar<float>(key.key() + "_GLOBALPOSITION_Z_" + region, cluster->globalPosition()[2]),
          Monitored::Scalar<int>(key.key() + "_CHANNELS_IN_PHI_" + region, cluster->channelsInPhi()),
          Monitored::Scalar<int>(key.key() + "_CHANNELS_IN_ETA_" + region, cluster->channelsInEta()),
          Monitored::Scalar<float>(key.key() + "_WIDTH_IN_ETA_" + region, cluster->widthInEta()),
          Monitored::Scalar<int>(key.key() + "_TOTAL_TOT_" + region, xAOD::xAODInDetMeasurement::Utilities::computeTotalToT(*cluster))
        );
      }
    }
  }
      
  for (std::size_t index = 0; index < m_stripKeys.size(); index++) {
    const SG::ReadHandleKey<xAOD::StripClusterContainer>& key = m_stripKeys[index];
    SG::ReadHandle<xAOD::StripClusterContainer> handle{key, ctx};
    ATH_CHECK(handle.isValid());
 
    for(auto cluster : *handle)
    {
      std::vector<std::string> regions {"all"};
      if(m_stripid->barrel_ec(cluster->rdoList()[0]) == 0) regions.push_back("barrel");
      else regions.push_back("endcap");
 
      for(auto const& region: regions)
      {
        Monitored::Group(
          m_monitoringTool,
          Monitored::Scalar<float>(key.key() + "_LOCALPOSITIONZOOM_X_" + region, cluster->localPosition<1>()(0,0)),
          Monitored::Scalar<float>(key.key() + "_LOCALPOSITION_X_" + region, cluster->localPosition<1>()(0,0)),
          Monitored::Scalar<float>(key.key() + "_LOCALCOVARIANCE_XX_" + region, cluster->localCovariance<1>()(0, 0)),
          Monitored::Scalar<float>(key.key() + "_GLOBALPOSITION_X_" + region, cluster->globalPosition()[0]),
          Monitored::Scalar<float>(key.key() + "_GLOBALPOSITION_Y_" + region, cluster->globalPosition()[1]),
          Monitored::Scalar<float>(key.key() + "_GLOBALPOSITION_Z_" + region, cluster->globalPosition()[2]),
          Monitored::Scalar<float>(key.key() + "_CHANNELS_IN_PHI_" + region, cluster->channelsInPhi())
        );
      }
    }
  }
 
  return StatusCode::SUCCESS;
}