StgcRawDataMonAlg.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/

// Package : sTgcRawDataMonAlg
// Author: Sebastian Fuenzalida Garrido
// Local supervisor: Edson Carquin Lopez
// Technical supervisor: Gerardo Vasquez
//
// DESCRIPTION:
// Subject: sTgc --> sTgc raw data monitoring
 
#include "StgcRawDataMonitoring/StgcRawDataMonAlg.h"

// *********************************************************************
// Public Methods
// *********************************************************************
 
sTgcRawDataMonAlg::sTgcRawDataMonAlg(const std::string& name, ISvcLocator* pSvcLocator) : AthMonitorAlgorithm(name, pSvcLocator) {
  //Declare the property
}
 
StatusCode sTgcRawDataMonAlg::initialize() {
  ATH_CHECK(AthMonitorAlgorithm::initialize());
  ATH_CHECK(m_idHelperSvc.retrieve());
  ATH_CHECK(m_sTgcContainerKey.initialize());
  ATH_CHECK(m_detectorManagerKey.initialize());
  ATH_CHECK(m_meTrkKey.initialize());
  ATH_CHECK(m_residualPullCalculator.retrieve());
  ATH_CHECK(m_muonKey.initialize());
  ATH_CHECK(m_rdoKey.initialize(SG::AllowEmpty));
 
  return StatusCode::SUCCESS;
}
 
StatusCode sTgcRawDataMonAlg::fillHistograms(const EventContext& ctx) const {
  const int lumiblock = GetEventInfo(ctx) -> lumiBlock();
 
  SG::ReadHandle<Muon::sTgcPrepDataContainer> sTgcContainer(m_sTgcContainerKey, ctx);
  SG::ReadCondHandle<MuonGM::MuonDetectorManager> detectorManagerKey(m_detectorManagerKey, ctx);
  SG::ReadHandle<xAOD::TrackParticleContainer> meTPContainer(m_meTrkKey, ctx);
  SG::ReadHandle<xAOD::MuonContainer> muonContainer(m_muonKey, ctx);
 
  if (!meTPContainer.isValid()) {
    ATH_MSG_FATAL("Could not get track particle container: " << m_meTrkKey.fullKey());
    return StatusCode::FAILURE;
  }
 
  if (!muonContainer.isValid()) {
    ATH_MSG_FATAL("Could not get muon container: " << m_muonKey.fullKey());
    return StatusCode::FAILURE;
  }
 
  if(!m_rdoKey.key().empty()){
    SG::ReadHandle<Muon::NSW_PadTriggerDataContainer> NSWpadTriggerContainer(m_rdoKey, ctx);
    if (!NSWpadTriggerContainer.isValid()) {
      ATH_MSG_FATAL("Could not get pad trigger data container: " << m_rdoKey.fullKey());
      return StatusCode::FAILURE;
    }
    fillsTgcPadTriggerDataHistograms(muonContainer.cptr(), NSWpadTriggerContainer.cptr(), lumiblock);
    fillsTgcPadTriggerEfficiencyHistograms(muonContainer.cptr(), NSWpadTriggerContainer.cptr(), detectorManagerKey.cptr());
  }
 
  ATH_MSG_DEBUG("Lumiblock: " << lumiblock);  
 
  fillsTgcClusterFromTrackHistograms(meTPContainer.cptr());
  fillsTgcEfficiencyHistograms(muonContainer.cptr(), detectorManagerKey.cptr());
  fillsTgcOccupancyHistograms(sTgcContainer.cptr(), detectorManagerKey.cptr(), lumiblock);
  fillsTgcLumiblockHistograms(sTgcContainer.cptr(), lumiblock);
 
  return StatusCode::SUCCESS;
}
 
void sTgcRawDataMonAlg::fillsTgcOccupancyHistograms(const Muon::sTgcPrepDataContainer* sTgcContainer, const MuonGM::MuonDetectorManager* muonDetectorManagerObject, int lb) const {
  ATH_MSG_DEBUG("ON fillsTgcOccupancyHistograms");
  for(const Muon::sTgcPrepDataCollection* coll : *sTgcContainer) {
    for (const Muon::sTgcPrepData* prd : *coll) {
      Identifier id = prd -> identify();
 
      if(!id.is_valid()) {
    ATH_MSG_WARNING("Invalid identifier found in Muon::sTgcPrepData");
    return;
      }
 
      std::string stationName = m_idHelperSvc -> stgcIdHelper().stationNameString(m_idHelperSvc -> stgcIdHelper().stationName(id));
      int stationEta          = m_idHelperSvc -> stgcIdHelper().stationEta(id);
      int stationPhi          = m_idHelperSvc -> stgcIdHelper().stationPhi(id);
      int multiplet           = m_idHelperSvc -> stgcIdHelper().multilayer(id);
      int gasGap              = m_idHelperSvc -> stgcIdHelper().gasGap(id);
      int channelType         = m_idHelperSvc -> stgcIdHelper().channelType(id);
      int sector              = m_idHelperSvc -> sector(id);
      std::string sectorStr   = std::to_string(sector);
      int sectorsTotal        = getSectors(id);
      int layer               = getLayer(multiplet, gasGap);
      int stationEtaAbs       = std::abs(stationEta);
      std::string layerStr    = std::to_string(layer);
      int iside               = (stationEta > 0) ? 1 : 0;      
      std::string side        = GeometricSectors::sTgcSide[iside];
      std::string channelName = "";
 
      if (channelType == sTgcIdHelper::sTgcChannelTypes::Pad) {
        channelName = "pad";
    int padNumber = m_idHelperSvc -> stgcIdHelper().channel(id);
    Identifier idPadQ1 = m_idHelperSvc -> stgcIdHelper().channelID(stationName, 1, stationPhi, multiplet, gasGap, channelType, 1);
    Identifier idPadQ2 = m_idHelperSvc -> stgcIdHelper().channelID(stationName, 2, stationPhi, multiplet, gasGap, channelType, 1);
    const MuonGM::sTgcReadoutElement* sTgcReadoutObjectPadQ1 = muonDetectorManagerObject -> getsTgcReadoutElement(idPadQ1);
    const MuonGM::sTgcReadoutElement* sTgcReadoutObjectPadQ2 = muonDetectorManagerObject -> getsTgcReadoutElement(idPadQ2);
    int maxPadNumberQ1 = sTgcReadoutObjectPadQ1 -> maxPadNumber(idPadQ1);
    int maxPadNumberQ2 = sTgcReadoutObjectPadQ2 -> maxPadNumber(idPadQ2);
 
    if (stationEtaAbs == 2) padNumber = padNumber + maxPadNumberQ1;
    else if (stationEtaAbs == 3) padNumber =  padNumber + maxPadNumberQ1 + maxPadNumberQ2;
        
    auto sectorMon    = Monitored::Scalar<int>("sector_layer_" + layerStr, sectorsTotal);
    auto padNumberMon = Monitored::Scalar<int>("padNumber_layer_" + layerStr, padNumber);
  auto padLbMon     = Monitored::Scalar<int>(channelName + "Lb", lb);
  auto padsectMon   = Monitored::Scalar<int>(channelName + "Sector", sectorsTotal);
 
    fill("Occupancy", sectorMon, padNumberMon);
 
  auto layerMon  = Monitored::Scalar<int>(channelName+"layer_" + side + sectorStr, layer);
  auto quadMon   = Monitored::Scalar<int>(channelName+"quad_" + side + sectorStr, stationEtaAbs);
  fill("Occupancy", layerMon, quadMon);
      }
 
      else if (channelType == sTgcIdHelper::sTgcChannelTypes::Strip) {
        channelName = "strip";        
    int stripNumber      = m_idHelperSvc -> stgcIdHelper().channel(id);
    Identifier idStripQ1 = m_idHelperSvc -> stgcIdHelper().channelID(stationName, 1, stationPhi, multiplet, gasGap, channelType, 1);
    Identifier idStripQ2 = m_idHelperSvc -> stgcIdHelper().channelID(stationName, 2, stationPhi, multiplet, gasGap, channelType, 1);
    const MuonGM::sTgcReadoutElement* sTgcReadoutObjectStripQ1 = muonDetectorManagerObject -> getsTgcReadoutElement(idStripQ1);
    const MuonGM::sTgcReadoutElement* sTgcReadoutObjectStripQ2 = muonDetectorManagerObject -> getsTgcReadoutElement(idStripQ2);
    int maxStripNumberQ1 = sTgcReadoutObjectStripQ1 -> numberOfStrips(idStripQ1);
    int maxStripNumberQ2 = sTgcReadoutObjectStripQ2 -> numberOfStrips(idStripQ2);
 
    if (stationEtaAbs == 2) stripNumber = stripNumber + maxStripNumberQ1 + 1;
    else if (stationEtaAbs == 3) stripNumber =  stripNumber + maxStripNumberQ1 + maxStripNumberQ2 + 1;
 
    auto sectorMon      = Monitored::Scalar<int>("sector_layer_" + layerStr, sectorsTotal);
    auto stripNumberMon = Monitored::Scalar<int>("stripNumber_layer_" + layerStr, stripNumber);
  auto stripLbMon     = Monitored::Scalar<int>(channelName + "Lb", lb);
  auto stripsectMon   = Monitored::Scalar<int>(channelName + "Sector", sectorsTotal);
 
    fill("Occupancy", sectorMon, stripNumberMon);
 
  auto layerMon  = Monitored::Scalar<int>(channelName+"layer_" + side + sectorStr, layer);
  auto quadMon   = Monitored::Scalar<int>(channelName+"quad_" + side + sectorStr, stationEtaAbs);
  fill("Occupancy", layerMon, quadMon);
      }
 
      else if (channelType == sTgcIdHelper::sTgcChannelTypes::Wire) {
        channelName = "wire";        
    int wireGroupNumber      = m_idHelperSvc -> stgcIdHelper().channel(id);
    Identifier idWireGroupQ3 = m_idHelperSvc -> stgcIdHelper().channelID("STL", 3, stationPhi, 1, 3, channelType, 1);
    const MuonGM::sTgcReadoutElement* sTgcReadoutObjectWireGroupQ3 = muonDetectorManagerObject -> getsTgcReadoutElement(idWireGroupQ3);
    int maxWireGroupNumberQ3 = sTgcReadoutObjectWireGroupQ3 -> numberOfStrips(idWireGroupQ3);
 
    auto stationEtaMon      = Monitored::Scalar<int>("stationEta_layer_" + layerStr, stationEta);
    auto wireGroupNumberMon = Monitored::Scalar<int>("wireGroupNumber_layer_" + layerStr, wireGroupNumber + (sector - 1)*maxWireGroupNumberQ3);
    auto wireLbMon          = Monitored::Scalar<int>(channelName + "Lb", lb);
  auto wiresectMon        = Monitored::Scalar<int>(channelName + "Sector", sectorsTotal);
  fill("Occupancy", stationEtaMon, wireGroupNumberMon);
 
  auto layerMon  = Monitored::Scalar<int>(channelName+"layer_" + side + sectorStr, layer);
  auto quadMon   = Monitored::Scalar<int>(channelName+"quad_" + side + sectorStr, stationEtaAbs);
  fill("Occupancy", layerMon, quadMon);
      }
      auto sectorMon = Monitored::Scalar<int>(channelName+"Sector", sectorsTotal);      
      auto febMon    = Monitored::Scalar<int>(channelName+"Feb", getFEBs(stationEta,layer));
      auto lbMon     = Monitored::Scalar<int>(channelName+"Lb", lb);
 
      fill("Occupancy", sectorMon, febMon);
      fill("LBShifterGroup", lbMon, sectorMon);    
    }
  }
  ATH_MSG_DEBUG("Off fillsTgcOccupancyHistograms");
}
 
 
void sTgcRawDataMonAlg::fillsTgcLumiblockHistograms(const Muon::sTgcPrepDataContainer* sTgcContainer, int lb) const {
  ATH_MSG_DEBUG("On fillsTgcLumiblockHistograms");
  for(const Muon::sTgcPrepDataCollection* coll : *sTgcContainer) {    
    for (const Muon::sTgcPrepData* prd : *coll) {
      Identifier id = prd -> identify();
 
 
      std::string stationName = m_idHelperSvc -> stgcIdHelper().stationNameString(m_idHelperSvc -> stgcIdHelper().stationName(id));
      int stationEta          = m_idHelperSvc -> stgcIdHelper().stationEta(id);
      int multiplet           = m_idHelperSvc -> stgcIdHelper().multilayer(id);
      int gasGap              = m_idHelperSvc -> stgcIdHelper().gasGap(id);
      int channelType         = m_idHelperSvc -> stgcIdHelper().channelType(id);
      int sector              = m_idHelperSvc -> sector(id);
      std::string sectorStr   = std::to_string(sector);
      int sectorsTotal        = getSectors(id);
      int numberOfHits        = coll->size();
      int layer               = getLayer(multiplet, gasGap);
      std::string layerStr    = std::to_string(layer);
      int iside               = (stationEta > 0) ? 1 : 0;      
      std::string side        = GeometricSectors::sTgcSide[iside];
      int feb                 = getFEBs(stationEta,layer);   
      std::string channelName = "";
 
      if (channelType == sTgcIdHelper::sTgcChannelTypes::Pad) channelName = "pad";
      else if (channelType == sTgcIdHelper::sTgcChannelTypes::Strip) channelName = "strip";
      else if (channelType == sTgcIdHelper::sTgcChannelTypes::Wire) channelName = "wire";
 
      auto febMon = Monitored::Scalar<int>(channelName+"FEBsector"+side+sectorStr, feb);
      auto lbMon  = Monitored::Scalar<int>(channelName+"LBsector"+side+sectorStr, lb);
      auto sectorMon = Monitored::Scalar<int>(channelName+"Sector", sectorsTotal); 
      auto lbbMon  = Monitored::Scalar<int>(channelName+"Lb", lb);
      auto numberOfHitsMon = Monitored::Scalar<int>(channelName+"NumberOfHits", numberOfHits);   
      
      fill("LBExpertGroup", lbMon, febMon);
    } 
  }
  ATH_MSG_DEBUG("Off fillsTgcLumiblockHistograms");  
}
 
 
void sTgcRawDataMonAlg::fillsTgcClusterFromTrackHistograms(const xAOD::TrackParticleContainer*  trkPartCont) const {
  ATH_MSG_DEBUG("On fillsTgcClusterFromTrackHistograms");  
  for (const xAOD::TrackParticle* meTP : *trkPartCont) {
    const Trk::Track* meTrack = meTP -> track();
    if(!meTrack) continue;
    for (const Trk::TrackStateOnSurface* trkState : *meTrack -> trackStateOnSurfaces()) {
      std::optional<std::tuple<Identifier, const Trk::RIO_OnTrack*>> status = getRotIdAndRotObject(trkState);
      if (!status.has_value()) continue;
      std::tuple<Identifier, const Trk::RIO_OnTrack*> rotIDtuple = status.value();
 
      Identifier rot_id = std::get<Identifier>(rotIDtuple);
 
      const Muon::sTgcClusterOnTrack* cluster = dynamic_cast<const Muon::sTgcClusterOnTrack*>(std::get<const Trk::RIO_OnTrack*>(rotIDtuple));
      if(!cluster) continue;
 
      const Muon::sTgcPrepData* prd = cluster -> prepRawData();
      if (!prd) continue;
 
      const MuonGM::sTgcReadoutElement* rdoEl = prd->detectorElement();
      if (!rdoEl) continue;
 
      const Amg::Vector2D& localPos = prd->localPosition();
 
      int channelType = m_idHelperSvc  -> stgcIdHelper().channelType(rot_id);
      int stEta       = m_idHelperSvc -> stgcIdHelper().stationEta(rot_id);
      int multi       = m_idHelperSvc  -> stgcIdHelper().multilayer(rot_id);
      int gap         = m_idHelperSvc  -> stgcIdHelper().gasGap(rot_id);
      int sector      = m_idHelperSvc -> sector(rot_id);
      int sectorsTotal        = getSectors(rot_id);
      int layer               = getLayer(multi, gap);
      int iside               = (stEta > 0) ? 1 : 0;
      std::string side        = GeometricSectors::sTgcSide[iside];
      std::string channelName = "";      
      std::string sectorStr   = std::to_string(sector);
 
      std::string quadLabel = "";
      bool isInner = false;
      bool isQ1 = (std::abs(stEta) == 1);
      if (isQ1) {
        isInner = rdoEl->isEtaZero(rot_id, localPos);
        quadLabel = isInner ? "Q1Inner" : "Q1Outer";
      } else {
        quadLabel = "Q" + std::to_string(std::abs(stEta));
      }
 
      if (channelType == sTgcIdHelper::sTgcChannelTypes::Pad) {
    float padCharge     = prd -> charge();
  auto padChargeMon = Monitored::Scalar<float>("padTrackCharge_" + side + "_quad_" + quadLabel + "_sector_" + std::to_string(sector) + "_layer_" + std::to_string(layer), padCharge);
  fill("padCharge_" + side + std::to_string(sector) + "_quad_" + quadLabel, padChargeMon);
 
    short int padTiming = prd -> time();
    auto padSectorSidedMon = Monitored::Scalar<int>("padTrackSectorSided_layer_" + std::to_string(layer), sectorsTotal);
    auto padTimingMon      = Monitored::Scalar<float>("padTrackTiming_layer_" + std::to_string(layer), padTiming);
    fill("sTgcTiming", padSectorSidedMon, padTimingMon);
 
    auto padSectorSidedExpertMon = Monitored::Scalar<int>("padTrackSectorSided_quad_" + std::to_string(std::abs(stEta)) + "_layer_" + std::to_string(layer), sectorsTotal);
    auto padTimingExpertMon      = Monitored::Scalar<float>("padTrackTiming_quad_" + std::to_string(std::abs(stEta)) + "_layer_" + std::to_string(layer), padTiming);
    fill("padTiming_quad_" + std::to_string(std::abs(stEta)), padSectorSidedExpertMon, padTimingExpertMon);
 
        ATH_MSG_DEBUG("Pad Timing: " << padTiming);
 
        auto febMon = Monitored::Scalar<int>("padFEB"+ side + sectorStr, getFEBs(stEta,layer));
        auto timeMon  = Monitored::Scalar<int>("padTiming"+ side + sectorStr, padTiming);
        fill("sTgcTiming", timeMon, febMon);              
      }
 
      else if (channelType == sTgcIdHelper::sTgcChannelTypes::Strip) {
        channelName = "strip";        
    const std::vector<Identifier>& stripIds = prd->rdoList();
    unsigned int csize = stripIds.size();
 
    std::vector<short int> stripTimesVec = prd -> stripTimes();
    std::vector<int> stripChargesVec = prd -> stripCharges();
 
    float stripClusterTimes = 0;
    float stripClusterCharges = 0;
 
    for (unsigned int sIdx = 0; sIdx < csize; ++sIdx) {
      stripClusterTimes += stripTimesVec.at(sIdx);
      stripClusterCharges += stripChargesVec.at(sIdx);
    }
 
    stripClusterTimes /= stripTimesVec.size();
 
  auto stripClusterChargesPerSideQuadMon = Monitored::Scalar<float>("stripTrackCharge_" + side  + "_quad_" + quadLabel + "_sector_" + std::to_string(sector)  +  "_layer_" + std::to_string(layer), stripClusterCharges);
  fill("stripCharge_" + side + std::to_string(sector) + "_quad_" + quadLabel, stripClusterChargesPerSideQuadMon);
 
 
    auto stripClusterSectorSidedMon = Monitored::Scalar<int>("stripTrackSectorSided_layer_" + std::to_string(layer), sectorsTotal);
    auto stripClusterTimesMon       = Monitored::Scalar<float>("stripTrackTiming_layer_" + std::to_string(layer), stripClusterTimes);
    auto stripClusterSizeMon        = Monitored::Scalar<unsigned int>("stripTrackClusterSize_layer_" + std::to_string(layer), csize);
    fill("sTgcTiming", stripClusterSectorSidedMon, stripClusterTimesMon);
    fill("padTriggerExpert", stripClusterSectorSidedMon, stripClusterSizeMon);
 
        ATH_MSG_DEBUG("Strip Timing: " << stripClusterTimes);        
        
        auto febMon = Monitored::Scalar<int>("stripFEB"+ side + sectorStr, getFEBs(stEta,layer));
        auto timeMon  = Monitored::Scalar<int>("stripTiming"+ side + sectorStr, stripClusterTimes);
        fill("sTgcTiming", timeMon, febMon);      
        
 
    auto stripSectorSidedExpertMon = Monitored::Scalar<int>("stripTrackSectorSided_quad_" + std::to_string(std::abs(stEta)) + "_layer_" + std::to_string(layer), sectorsTotal);
    auto stripTimingExpertMon      = Monitored::Scalar<float>("stripTrackTiming_quad_" + std::to_string(std::abs(stEta)) + "_layer_" + std::to_string(layer), stripClusterTimes);
    fill("stripTiming_quad_" + std::to_string(std::abs(stEta)), stripSectorSidedExpertMon, stripTimingExpertMon);
 
 
    std::optional<Trk::ResidualPull> resPull(m_residualPullCalculator -> residualPull(trkState -> measurementOnTrack(), trkState -> trackParameters(), Trk::ResidualPull::ResidualType::Biased));
 
    if (resPull) {
      float residual = resPull -> residual()[Trk::locX];
      auto residualMon  = Monitored::Scalar<float>("residual_" + side + "_quad_" + std::to_string(std::abs(stEta)) + "_sector_" + std::to_string(sector) + "_layer_" + std::to_string(layer), residual);
      fill("sTgcResiduals_" + side + std::to_string(sector) + "_quad_" + std::to_string(std::abs(stEta)), residualMon);
    }        
      }
 
      else if (channelType == sTgcIdHelper::sTgcChannelTypes::Wire) {
    float wireGroupCharge     = prd -> charge();
    auto wireGroupChargeMon = Monitored::Scalar<float>("wireGroupTrackCharge_" + side + "_quad_" + std::to_string(std::abs(stEta)) + "_sector_" + std::to_string(sector)  + "_layer_" + std::to_string(layer), wireGroupCharge);
    fill("wireGroupCharge_" + side + std::to_string(sector) + "_quad_" + std::to_string(std::abs(stEta)), wireGroupChargeMon);
 
    short int wireGroupTiming = prd -> time();
    auto wireGroupSectorSidedMon = Monitored::Scalar<int>("wireGroupTrackSectorSided_layer_" + std::to_string(layer), sectorsTotal);
    auto wireGroupTimingMon      = Monitored::Scalar<float>("wireGroupTrackTiming_layer_" + std::to_string(layer), wireGroupTiming);
    fill("sTgcTiming", wireGroupSectorSidedMon, wireGroupTimingMon);
 
    auto wireSectorSidedExpertMon = Monitored::Scalar<int>("wireTrackSectorSided_quad_" + std::to_string(std::abs(stEta)) + "_layer_" + std::to_string(layer), sectorsTotal);
    auto wireTimingExpertMon      = Monitored::Scalar<float>("wireTrackTiming_quad_" + std::to_string(std::abs(stEta)) + "_layer_" + std::to_string(layer), wireGroupTiming);
    fill("wireTiming_quad_" + std::to_string(std::abs(stEta)), wireSectorSidedExpertMon, wireTimingExpertMon);
 
        ATH_MSG_DEBUG("Wire Timing: " << wireGroupTiming);        
 
        auto febMon = Monitored::Scalar<int>("wireFEB"+ side + sectorStr, getFEBs(stEta,layer));
        auto timeMon  = Monitored::Scalar<int>("wireTiming"+ side + sectorStr, wireGroupTiming);
        fill("sTgcTiming", timeMon, febMon);              
      }      
    }
  }
  ATH_MSG_DEBUG("Off fillsTgcClusterFromTrackHistograms");  
}
 
void sTgcRawDataMonAlg::fillsTgcPadTriggerDataHistograms(const xAOD::MuonContainer*  muonContainer, const Muon::NSW_PadTriggerDataContainer* NSWpadTriggerObject, int lb) const {
 ATH_MSG_DEBUG("On fillsTgcPadTriggerDataHistograms");  
 for (const Muon::NSW_PadTriggerData* rdo : *NSWpadTriggerObject) {
    bool sideA = rdo -> sideA();
    bool largeSector = rdo -> largeSector();
 
    int iside = (sideA) ? 1 : 0;
    int isize = (largeSector) ? 1 : 0;
 
    std::string side = GeometricSectors::sTgcSide[iside];
    std::string size = GeometricSectors::sTgcSize[isize];
 
    size_t numberOfTriggers = rdo -> getNumberOfTriggers();
    size_t numberOfHits     = rdo -> getNumberOfHits();
 
   for (size_t trigger = 0; trigger < numberOfTriggers; ++trigger) {
      int triggerPhiIdsUnsigned = rdo -> getTriggerPhiIds().at(trigger);
      int triggerBandIds        = rdo -> getTriggerBandIds().at(trigger);
      int triggerRelBCID        = rdo -> getTriggerRelBcids().at(trigger);
      int sourceId     = rdo -> getSourceid();
      int sectorNumber = sourceidToSector(sourceId, sideA);
 
      if (triggerPhiIdsUnsigned == m_cutTriggerPhiId || triggerBandIds == m_cutTriggerBandId) continue;
 
      int triggerPhiIds = getSignedPhiId(triggerPhiIdsUnsigned);
 
      auto phiIdsPerSideSizeMon  = Monitored::Scalar<int>("phiIds_"  + side + "_" + size, triggerPhiIds);
      auto bandIdsPerSideSizeMon = Monitored::Scalar<int>("bandIds_" + side + "_" + size, triggerBandIds);
      fill("padTriggerExpert", phiIdsPerSideSizeMon, bandIdsPerSideSizeMon);
 
      auto lbMon = Monitored::Scalar<int>("lb", lb);
      auto relBCIDMon = Monitored::Scalar<int>("relBCID", triggerRelBCID);
      auto sectorMon = Monitored::Scalar<int>("sector", sectorNumber);
      auto numberOfTriggersMon = Monitored::Scalar<int>("numberOfTriggers", numberOfTriggers);
      fill("Overview", lbMon, sectorMon, numberOfTriggersMon);
      fill("padTriggerExpert", lbMon, relBCIDMon);
 
      auto numberOfTriggersPerSectorMon = Monitored::Scalar<int>("numberOfTriggers_" + side + "_sector_" + std::to_string(std::abs(sectorNumber)), numberOfTriggers);
      auto phiIdsSidedSizedPerSectorMon  = Monitored::Scalar<int>("phiIds_"  + side + "_sector_" + std::to_string(std::abs(sectorNumber)), triggerPhiIds);
      auto bandIdsSidedSizedPerSectorMon = Monitored::Scalar<int>("bandIds_" + side + "_sector_" + std::to_string(std::abs(sectorNumber)), triggerBandIds);
      auto lbPerSectorMon = Monitored::Scalar<int>("lb_" + side + "_sector_" + std::to_string(std::abs(sectorNumber)), lb);
      auto relBCIDperSectorMon = Monitored::Scalar<int>("relBCID_"  + side + "_sector_" + std::to_string(std::abs(sectorNumber)), triggerRelBCID);
      fill("padTriggerExpert", numberOfTriggersPerSectorMon, phiIdsSidedSizedPerSectorMon, bandIdsSidedSizedPerSectorMon, lbPerSectorMon, relBCIDperSectorMon);
 
      auto RelBCIDPerSectorMon = Monitored::Scalar<int>("relBCID_"+ side + "_sector_" + std::to_string(std::abs(sectorNumber)), triggerRelBCID);
      auto PhiIDPerSectorMon   = Monitored::Scalar<int>("phiIds_"  + side + "_sector_" + std::to_string(std::abs(sectorNumber)), triggerPhiIds);
      auto BandIDPerSectorMon  = Monitored::Scalar<int>("bandID_" + side + "_sector_" + std::to_string(std::abs(sectorNumber)), triggerBandIds);
      fill("padTriggerExpert", RelBCIDPerSectorMon, PhiIDPerSectorMon, BandIDPerSectorMon);
   } 
 
   for (size_t hits = 0; hits < numberOfHits; ++hits){
      std::optional<Identifier> status = getPadId(rdo->getSourceid(), rdo->getHitPfebs().at(hits), rdo->getHitTdsChannels().at(hits));
      if (!status.has_value()) continue;
      Identifier pad_id = status.value();
 
      for (const xAOD::Muon* mu : *muonContainer) {
         if(mu -> pt() < m_cutPt) continue;
         if(!(mu -> author() == xAOD::Muon::Author::MuidCo || mu -> author() == xAOD::Muon::Author::MuidSA)) continue;
        
         const xAOD::TrackParticle* meTP = mu -> trackParticle(xAOD::Muon::TrackParticleType::ExtrapolatedMuonSpectrometerTrackParticle);
         if(meTP == nullptr) continue;
 
         const Trk::Track* meTrack = meTP -> track();
         if(!meTrack) continue;
 
         for(const Trk::TrackStateOnSurface* trkState : *meTrack->trackStateOnSurfaces()) {
            std::optional<Identifier> status = getRotId(trkState);
            if (!status.has_value()) continue;
            Identifier rot_id = status.value();
 
            int channelType = m_idHelperSvc -> stgcIdHelper().channelType(rot_id);
            if (channelType != sTgcIdHelper::sTgcChannelTypes::Pad) continue;
          
 
        if (rot_id != pad_id) continue;
 
        int sourceIds = rdo -> getSourceid();
        int sectorNumbers = sourceidToSector(sourceIds, sideA);
        int hitRelBCID = rdo -> getHitRelBcids().at(hits);
        int hitpfebs = rdo -> getHitPfebs().at(hits);
        int hitTdsChannels = rdo->getHitTdsChannels().at(hits);
 
        std::optional<std::tuple<int, int, std::string, std::string, int>> statusPadEtaPhi = getPadEtaPhiTuple(sourceIds, hitpfebs, hitTdsChannels);
        if (!statusPadEtaPhi.has_value()) continue;
        std::tuple<int, int, std::string, std::string, int> padEtaPhiTuple = statusPadEtaPhi.value();
 
        int padPhi = std::get<0>(padEtaPhiTuple);
        int padEta = std::get<1>(padEtaPhiTuple);
        std::string sideName = std::get<2>(padEtaPhiTuple);
        std::string sizeName = std::get<3>(padEtaPhiTuple);
        int layer = std::get<4>(padEtaPhiTuple);
                
        auto padPhiMon = Monitored::Scalar<int>("padPhi_" + sideName + "_" + sizeName + "_layer_" + std::to_string(layer), padPhi);
        auto padEtaMon = Monitored::Scalar<int>("padEta_" + sideName + "_" + sizeName + "_layer_" + std::to_string(layer), padEta);
        fill("padTriggerOccupancy", padPhiMon, padEtaMon);
 
        auto hitRelBCIDmon = Monitored::Scalar<int>("hitRelBCID", hitRelBCID);
        auto hitPfebsMon   = Monitored::Scalar<int>("hitPfebs", hitpfebs);
        auto sectorMon     = Monitored::Scalar<int>("sector", sectorNumbers);
        fill("padTriggerExpert", hitRelBCIDmon, hitPfebsMon, sectorMon);
 
        auto hitRelBCIDPerSectorMon = Monitored::Scalar<int>("hitRelBCID_"+side+"_sector_"+std::to_string(std::abs(sectorNumbers)), hitRelBCID);
        auto hitpfebsPerSectorMon = Monitored::Scalar<int>("hitPfebs_"+side+"_sector_"+std::to_string(std::abs(sectorNumbers)), hitpfebs);
        fill("padTriggerExpert", hitRelBCIDPerSectorMon, hitpfebsPerSectorMon);
         } // end TrackStateOnSurface loop
      } // end Muon loop
    } // end NSW_PadTriggerData loop
  }// end Number of Hits loop
  ATH_MSG_DEBUG("Off fillsTgcPadTriggerDataHistograms");   
} 
 
void sTgcRawDataMonAlg::fillsTgcEfficiencyHistograms(const xAOD::MuonContainer*  muonContainer, const MuonGM::MuonDetectorManager* muonDetectorManagerObject) const {
  ATH_MSG_DEBUG("On fillsTgcEfficiencyHistograms");    
  for (const xAOD::Muon* mu : *muonContainer) {
    if (mu -> pt() < m_cutPt) continue;
    if (!(mu -> author() == xAOD::Muon::Author::MuidCo || mu -> author() == xAOD::Muon::Author::MuidSA)) continue;
   
    struct sTGCeff {
      std::array<int, 8> quadMultiplet{};
      std::array<int, 8> layerMultiplet{};
      std::array<float, 8> xPosMultiplet{};
      std::array<float, 8> yPosMultiplet{};
      std::array<float, 8> zPosMultiplet{};
    };
 
    std::array<std::array<sTGCeff, 16>, 2> effPlots; // Store active layers per side (2) and sectors (16) 
        
    const xAOD::TrackParticle* meTP = mu -> trackParticle(xAOD::Muon::TrackParticleType::ExtrapolatedMuonSpectrometerTrackParticle);
    if(meTP == nullptr) continue;
 
    const Trk::Track* meTrack = meTP -> track();
    if(!meTrack) continue;
 
    for(const Trk::TrackStateOnSurface* trkState : *meTrack->trackStateOnSurfaces()) {
      std::optional<Identifier> status = getRotId(trkState);
      if (!status.has_value()) continue;
 
      Identifier rot_id = status.value();
 
      int channelType = m_idHelperSvc -> stgcIdHelper().channelType(rot_id);
      if (channelType != sTgcIdHelper::sTgcChannelTypes::Strip) continue;
 
      int stEta  = m_idHelperSvc -> stgcIdHelper().stationEta(rot_id);
      int iside  = (stEta > 0) ? 1 : 0;
      int sector = m_idHelperSvc -> sector(rot_id);
      int multi  = m_idHelperSvc -> stgcIdHelper().multilayer(rot_id);
      int gap    = m_idHelperSvc -> stgcIdHelper().gasGap(rot_id);
      int layer  = getLayer(multi, gap);
            
      const Amg::Vector2D& positionsMultiplet = (trkState) -> trackParameters() -> localPosition();
      float xPosStripInMultipletLocal = positionsMultiplet.x();
      float yPosStripInMultipletLocal = positionsMultiplet.y();
 
      Amg::Vector2D localPos(xPosStripInMultipletLocal, yPosStripInMultipletLocal);
      Amg::Vector3D globalPos(Amg::Vector3D::Zero());
      const MuonGM::sTgcReadoutElement* sTgcReadoutObjectStrip = muonDetectorManagerObject -> getsTgcReadoutElement(rot_id);
      sTgcReadoutObjectStrip -> surface(rot_id).localToGlobal(localPos, Amg::Vector3D::Zero(), globalPos);
      float xPosStripInMultiplet = globalPos.x();
      float yPosStripInMultiplet = globalPos.y();
      
      Amg::Vector3D posStripGlobal{Amg::Vector3D::Zero()};
      (muonDetectorManagerObject -> getsTgcReadoutElement(rot_id)) -> stripGlobalPosition(rot_id, posStripGlobal);
      float zPosStripInMultiplet = posStripGlobal.z();
 
      auto& sTGCelements = effPlots[iside][sector - 1];
 
      sTGCelements.quadMultiplet.at(layer - 1) = stEta;
      sTGCelements.layerMultiplet.at(layer - 1) = layer;
      sTGCelements.xPosMultiplet.at(layer - 1) = xPosStripInMultiplet;
      sTGCelements.yPosMultiplet.at(layer - 1) = yPosStripInMultiplet;
      sTGCelements.zPosMultiplet.at(layer - 1) = zPosStripInMultiplet;      
    } // end track state loop
 
    for (unsigned int isideIndex = 0; isideIndex <= 1; ++isideIndex) {
      for (unsigned int sectorIndex = 1; sectorIndex <= 16; ++sectorIndex) {
        auto& sTGCelements = effPlots[isideIndex][sectorIndex - 1];
        bool fourOutEight = std::count_if(sTGCelements.layerMultiplet.begin(), sTGCelements.layerMultiplet.end(), [](int i) { return i != 0; }) >= 4;
        bool oneRefLayer = std::count_if(sTGCelements.layerMultiplet.begin(), sTGCelements.layerMultiplet.end(), [](int i) { return i != 0; }) >= 1;
        
        if (fourOutEight) {
          for (auto layerIndex = static_cast<std::array<int, 8>::size_type>(1); layerIndex <= sTGCelements.layerMultiplet.size(); ++layerIndex) {
            if (sTGCelements.layerMultiplet.at(layerIndex - 1) == 0) continue;
            
            float xPos = sTGCelements.xPosMultiplet.at(layerIndex - 1);
            float yPos = sTGCelements.yPosMultiplet.at(layerIndex - 1);
            float rPos = std::hypot(xPos, yPos);
            
            std::string side = GeometricSectors::sTgcSide[isideIndex];
            
            auto effQuestionMon = Monitored::Scalar<bool>("hitLayer", true);
            auto rPosStripMon = Monitored::Scalar<float>("rPosStrip_" + side + "_sector_" + std::to_string(sectorIndex)  + "_layer_" + std::to_string(layerIndex), rPos);
            fill("rPosStrip_" + side + std::to_string(sectorIndex), rPosStripMon, effQuestionMon);
            //GlobalRGroup
            
            auto xPosStripmon = Monitored::Scalar<float>("xPosStrip_" + side + "_layer_" + std::to_string(layerIndex), xPos);
            auto yPosStripmon = Monitored::Scalar<float>("yPosStrip_" + side + "_layer_" + std::to_string(layerIndex), yPos);
            fill("padTriggerShifter", xPosStripmon, yPosStripmon, effQuestionMon);
            //StripEfficiency
          } // End of loop over efficient layers
        } // End of efficient case
        
        else if (!fourOutEight && oneRefLayer) {
          auto refLayerIndex = std::distance(sTGCelements.layerMultiplet.begin(), std::find_if(sTGCelements.layerMultiplet.begin(), sTGCelements.layerMultiplet.end(), [](int i) {return i != 0;}));
 
          for (auto layerIndex = static_cast<std::array<int, 8>::size_type>(1); layerIndex <= sTGCelements.layerMultiplet.size(); ++layerIndex) {
            if (sTGCelements.layerMultiplet.at(layerIndex - 1) != 0) continue;
            int quad = sTGCelements.quadMultiplet.at(refLayerIndex);
            int multi = (layerIndex <= static_cast<std::array<int, 8>::size_type>(m_idHelperSvc->stgcIdHelper().gasGapMax())) ? m_idHelperSvc->stgcIdHelper().multilayerMin() : m_idHelperSvc->stgcIdHelper().multilayerMax();
            int gap = (layerIndex <= static_cast<std::array<int, 8>::size_type>(m_idHelperSvc->stgcIdHelper().gasGapMax())) ? layerIndex : layerIndex - static_cast<std::array<int, 8>::size_type>(m_idHelperSvc->stgcIdHelper().gasGapMax()); 
                        
            bool isValid = false;
            const Identifier idProbe = m_idHelperSvc -> stgcIdHelper().channelID((sectorIndex % 2 == 0) ? "STS" : "STL", quad, (sectorIndex % 2 == 0) ? sectorIndex/2 : (sectorIndex + 1)/2, multi, gap, sTgcIdHelper::sTgcChannelTypes::Strip, 1, isValid);
            
            if (!isValid) {
              ATH_MSG_WARNING("Identifier of probe layer is invalid");
              continue;
            }
            
            Amg::Vector3D posProbe{Amg::Vector3D::Zero()};
            (muonDetectorManagerObject -> getsTgcReadoutElement(idProbe)) -> stripGlobalPosition(idProbe, posProbe);
            float posZprobe = posProbe.z();
              
            float xSlope = sTGCelements.xPosMultiplet.at(refLayerIndex)/sTGCelements.zPosMultiplet.at(refLayerIndex);
            float ySlope = sTGCelements.yPosMultiplet.at(refLayerIndex)/sTGCelements.zPosMultiplet.at(refLayerIndex);
            
            float xPos = sTGCelements.xPosMultiplet.at(refLayerIndex) + xSlope*(posZprobe - sTGCelements.zPosMultiplet.at(refLayerIndex));
            float yPos = sTGCelements.yPosMultiplet.at(refLayerIndex) + ySlope*(posZprobe - sTGCelements.zPosMultiplet.at(refLayerIndex));
            float rPos = std::hypot(xPos, yPos);
 
            std::string side = GeometricSectors::sTgcSide[isideIndex];
           
            auto effQuestionMon = Monitored::Scalar<bool>("hitLayer", false);
            
            auto rPosStripProbemon = Monitored::Scalar<float>("rPosStrip_" + side + "_sector_" + std::to_string(sectorIndex)  + "_layer_" + std::to_string(layerIndex), rPos);
            fill("rPosStrip_" + side + std::to_string(sectorIndex), rPosStripProbemon, effQuestionMon);
              
            auto xPosStripProbemon = Monitored::Scalar<float>("xPosStrip_" + side + "_layer_" + std::to_string(layerIndex), xPos);
            auto yPosStripProbemon = Monitored::Scalar<float>("yPosStrip_" + side + "_layer_" + std::to_string(layerIndex), yPos);
            fill("padTriggerShifter", xPosStripProbemon, yPosStripProbemon, effQuestionMon);
          } // End of loop over probe layers
        } // End of non-efficient case
      } // End of sector loop
    } // End of iside loop
  } // End muon container loop
  ATH_MSG_DEBUG("Off fillsTgcClusterFromTrackHistograms");
} // end stgc strip function
 
void sTgcRawDataMonAlg::fillsTgcPadTriggerEfficiencyHistograms(const xAOD::MuonContainer*  muonContainer, const Muon::NSW_PadTriggerDataContainer* NSWpadTriggerObject, const MuonGM::MuonDetectorManager* muonDetectorManagerObject) const {
  ATH_MSG_DEBUG("On fillsTgcPadTriggerEfficiencyHistograms");  
  for (const xAOD::Muon* mu : *muonContainer) {
    if (!(mu -> author() == xAOD::Muon::Author::MuidCo || mu -> author() == xAOD::Muon::Author::MuidSA)) continue;
    if (mu -> pt() < m_cutPt) continue;
    if (std::abs(mu -> eta()) < m_cutEtaDown || std::abs(mu -> eta()) > m_cutEtaUp) continue; 
 
    double recoMuonEta = mu -> eta();
    double recoMuonPhi = mu -> phi();
 
    std::string sideRecoMuon = GeometricSectors::sTgcSide[recoMuonEta > 0];
 
    std::string minSideRecoMuon = "", minSideTrigger = "";
    
    double minTriggerEta = 999., minTriggerPhi = 999.,
      minRecoEta = 999., minRecoPhi = 999., minDeltaR = 999.;
    
    auto minDeltaRtrigIt = -1;
  
    for (const Muon::NSW_PadTriggerData* rdo : *NSWpadTriggerObject) {
      bool sideA = rdo -> sideA();
      bool largeSector = rdo -> largeSector();     
      std::string sideTrigger = GeometricSectors::sTgcSide[sideA];
      size_t numberOfTriggers = rdo -> getNumberOfTriggers();
      
      for (size_t trigger = 0; trigger < numberOfTriggers; ++trigger) {
        int triggerPhiIdsUnsigned = rdo -> getTriggerPhiIds().at(trigger);    
        int triggerBandIds = rdo -> getTriggerBandIds().at(trigger);
        int sourceId = rdo -> getSourceid();
        
        if (triggerPhiIdsUnsigned == m_cutTriggerPhiId || triggerBandIds == m_cutTriggerBandId) continue;
        
        int triggerPhiIds = getSignedPhiId(triggerPhiIdsUnsigned);
        std::optional<double> status = bandId2eta(triggerBandIds, largeSector, sideA, muonDetectorManagerObject);
        if (!status.has_value()) continue;
        double triggerBandIdToEta = status.value();
        double triggerPhiIDtoPhi = triggersectorphiid2phi(sourceId, triggerPhiIds);
        double deltaR =  xAOD::P4Helpers::deltaR(recoMuonEta, recoMuonPhi, triggerBandIdToEta, triggerPhiIDtoPhi); 
 
        if (std::abs(triggerBandIdToEta) < m_cutEtaDown || std::abs(triggerBandIdToEta) > m_cutEtaUp) continue;
    
        if (sideRecoMuon == sideTrigger) { 
          if (deltaR < minDeltaR) {
            minSideRecoMuon = sideRecoMuon;
            minSideTrigger = sideTrigger;
            minTriggerEta = triggerBandIdToEta;
            minTriggerPhi = triggerPhiIDtoPhi;
            minRecoEta = recoMuonEta;
            minRecoPhi = recoMuonPhi;
            minDeltaR = deltaR;
            minDeltaRtrigIt = trigger;
          }
        }
      } // end number of triggers loop
    } // end pad trigger data container loop
    
    bool muonRecoTriggerMatch = false;
    
    if (minDeltaRtrigIt != -1) {
      if (minDeltaR < m_minDeltaR) {
        muonRecoTriggerMatch = true;
      }
    }
        
    auto deltaRmon = Monitored::Scalar<double>("deltaR", minDeltaR);
    fill("Overview", deltaRmon);
    
    auto etaRecoMuonMon = Monitored::Scalar<double>("etaRecoMuon", minRecoEta);
    auto phiRecoMuonMon = Monitored::Scalar<double>("phiRecoMuon", minRecoPhi);
    fill("Overview", etaRecoMuonMon, phiRecoMuonMon);
    
    auto etaPadTriggerMon = Monitored::Scalar<double>("etaPadTrigger", minTriggerEta);
    auto phiPadTriggerMon = Monitored::Scalar<double>("phiPadTrigger", minTriggerPhi);
    fill("Overview", etaPadTriggerMon, phiPadTriggerMon);
 
    auto phiRecoMuonSidedMon = Monitored::Scalar<double>("phiRecoMuon_" + minSideRecoMuon, minRecoPhi);
    auto phiPadTriggerSidedMon = Monitored::Scalar<double>("phiPadTrigger_" + minSideTrigger, minTriggerPhi);
    fill("Overview", phiRecoMuonSidedMon, phiPadTriggerSidedMon);
    
    auto muonRecoTriggerMatchMon = Monitored::Scalar<bool>("muonRecoTriggerMatch", muonRecoTriggerMatch);
    auto etaRecoMuonEffMon = Monitored::Scalar<double>("etaRecoMuonEff", minRecoEta);
    auto phiRecoMuonEffMon = Monitored::Scalar<double>("phiRecoMuonEff", minRecoPhi);
    fill("Overview", muonRecoTriggerMatchMon, etaRecoMuonEffMon, phiRecoMuonEffMon);
 
  } // end muon container loop
  ATH_MSG_DEBUG("Off fillsTgcPadTriggerEfficiencyHistograms");    
}