TgcRawDataMonitorTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
#include "TgcRawDataMonitorTool.h"
#include "TVector2.h"
#include "MuonReadoutGeometry/TgcReadoutElement.h"
#include "TrkSurfaces/TrapezoidBounds.h"
 
TgcRawDataMonitorTool::TgcRawDataMonitorTool( const std::string& type,
                          const std::string& name,
                          const IInterface* parent)
  : base_class(type,name,parent)
{
}
 
TgcRawDataMonitorTool::~TgcRawDataMonitorTool(){}
StatusCode TgcRawDataMonitorTool::initialize() {
  ATH_MSG_DEBUG("TgcRawDataMonitorTool::initialize(): Start");
  ATH_CHECK(m_DetectorManagerKey.initialize());
  ATH_MSG_DEBUG("TgcRawDataMonitorTool::initialize(): End");
  return StatusCode::SUCCESS;
}
std::set<std::string> TgcRawDataMonitorTool::getPassedChambers(const Amg::Vector3D& extVec) const {
  ATH_MSG_DEBUG("TgcRawDataMonitorTool::getPassedChambers(): Start");  
  double ext_x = extVec[0];
  double ext_y = extVec[1];
  double ext_z = extVec[2];
  std::set<std::string> passed_chambers;
  SG::ReadCondHandle<MuonGM::MuonDetectorManager> DetectorManagerHandle{m_DetectorManagerKey,Gaudi::Hive::currentContext()};
  const MuonGM::MuonDetectorManager* muonMgr = DetectorManagerHandle.cptr();
  if(muonMgr==nullptr){
    ATH_MSG_ERROR("Null pointer to the read MuonDetectorManager conditions object.");
    return passed_chambers;
  }     
  for(int iStationName = 41 ; iStationName <= 48 ; iStationName++){
    if ( (iStationName==41 || iStationName==42) && // M1
     ( std::abs(ext_z)<m_M1Mod1Pl1L1_offsetZ-10. || std::abs(ext_z)>m_M1Mod2Pl2L3_offsetZ+10. ) )continue;
    if ( (iStationName==43 || iStationName==44) && // M2
     ( std::abs(ext_z)<m_M2Mod1Pl1L1_offsetZ-10. || std::abs(ext_z)>m_M2Mod2Pl2L2_offsetZ+10. ) )continue;
    if ((iStationName==45 || iStationName==46) && // M3
    ( std::abs(ext_z)<m_M3Mod1Pl1L1_offsetZ-10. || std::abs(ext_z)>m_M3Mod2Pl2L2_offsetZ+10. ) )continue;
    if ( iStationName==48 && // EI
     ( std::abs(ext_z)<m_EIMod1Pl1L1_offsetZ-10. || std::abs(ext_z)>m_EIMod2Pl1L2_offsetZ+10. ) )continue;
    if ( iStationName==47 && // FI
     ( std::abs(ext_z)<m_FIMod1Pl1L1_offsetZ-10. || std::abs(ext_z)>m_FIMod2Pl1L2_offsetZ+10. ) )continue;
    bool endcap = (iStationName%2==0);
    int maxEta = 1;
    if(endcap && iStationName==42) maxEta = 4;
    else if(endcap && (iStationName==44 || iStationName==46)) maxEta = 5;
    for(int iStationEta = -maxEta ; iStationEta <= maxEta ; iStationEta++){
      if(iStationEta==0)continue;
      if(ext_z > 0 && iStationEta < 0)continue;
      if(ext_z < 0 && iStationEta > 0)continue;
      for(int iStationPhi = muonMgr->tgcIdHelper()->stationPhiMin(endcap)  ;
      iStationPhi <= muonMgr->tgcIdHelper()->stationPhiMax(endcap) ; iStationPhi++){
    if(iStationName==48 && iStationPhi>21)continue;
    TGC::TgcChamber tgcChamber; tgcChamber.initChamber(iStationEta, iStationPhi, iStationName);
     bool isValid{false};
     const Identifier id = muonMgr->tgcIdHelper()->elementID(iStationName, iStationEta, iStationPhi, isValid);
     if (!isValid) continue;
    
    const MuonGM::TgcReadoutElement* detEle = muonMgr->getTgcReadoutElement(id);
    const Trk::TrapezoidBounds* tb = static_cast<const Trk::TrapezoidBounds*>(&detEle->bounds());
    Amg::Vector3D trkLocVec3D = detEle->transform().inverse() * Amg::Vector3D(ext_x,ext_y,ext_z);
    if(std::abs(trkLocVec3D.x())>1000.)continue;
    if(std::abs(trkLocVec3D.y())>1500.)continue;
    Amg::Vector2D trkvec(trkLocVec3D.x(),trkLocVec3D.y());
    if( tb->minDistance(trkvec) < m_fiduciFrame ) passed_chambers.insert(tgcChamber.cham_name());
      }
    }
  }
 
  ATH_MSG_DEBUG("TgcRawDataMonitorTool::getPassedChambers(): End");
  ATH_MSG_DEBUG("Number of passed TGC chambers: " << passed_chambers.size());
  return passed_chambers;
}
 
bool TgcRawDataMonitorTool::getMapIndex(const TGC::TgcHit& tgcHit,
                    int& etamap_index, int& phimap_index, int& phimap_global_index ) const {
 
  ATH_MSG_DEBUG("TgcRawDataMonitorTool::getMapIndex(const TGC::TgcHit& tgcHit,: Start");  
  return getMapIndex(tgcHit,tgcHit.gasGap(),etamap_index,phimap_index,phimap_global_index);
  ATH_MSG_DEBUG("TgcRawDataMonitorTool::getMapIndex(const TGC::TgcHit& tgcHit,: End");  
}
bool TgcRawDataMonitorTool::getMapIndex(const TGC::TgcChamber& tgcCham, int iLay,
                    int& etamap_index, int& phimap_index, int& phimap_global_index ) const {
  ATH_MSG_DEBUG("TgcRawDataMonitorTool::getMapIndex(const TGC::TgcChamber& tgcCham,: Start");  
  if (tgcCham.iM() > 1 && iLay > 2){
    ATH_MSG_DEBUG("tgcCham.iM() > 1 && iLay > 2");  
    return false;
  }
  int nLay = (tgcCham.iM()==1)?(3):(2);
  if (tgcCham.iM() == 4) { // EI/FI
    etamap_index = tgcCham.iEta() * nLay + iLay; // 1..4
    if(tgcCham.iEta() == 0){ //FI
      phimap_index = tgcCham.iPhi(); // 1..24
    }else{ // EI 
      phimap_index = tgcCham.iPhi() + int(tgcCham.iSec() / 2) + tgcCham.iSec(); // 1..24
    }
    phimap_global_index = tgcCham.iEta() * 24*nLay + (phimap_index - 1) * nLay + iLay; // 1..96
  }else{ // BigWheel
    etamap_index = tgcCham.iEta() * nLay + iLay; // 1..15 for M1, 1..12 for M2M3
    phimap_index = (tgcCham.iSec() - 1) * 4 + tgcCham.iPhi() + 1; // 1..48
    phimap_global_index = tgcCham.iEta() * 48*nLay + (phimap_index - 1) * nLay + iLay; // 1..720 for M1, 1..576 for M2M3
  }
  ATH_MSG_DEBUG("TgcRawDataMonitorTool::getMapIndex(const TGC::TgcChamber& tgcCham,: End");
  return true;
}
 
bool TgcRawDataMonitorTool::getMapIndexOnline(const std::string& chamber_type_name, // e.g. A01M01f00E01L03W
                          int& etamap_index, int& phimap_index, int& phimap_global_index) const {
  etamap_index=-1;
  phimap_index=-1;
  phimap_global_index=-1;
  if( chamber_type_name.size() != 16 ) return false; // e.g. A01M01f01E01L01S
 
  const int offset = (chamber_type_name[5]=='1')?(0):((chamber_type_name[5]=='2')?(3):((chamber_type_name[5]=='3')?(5):(7)));
  const int layerNumber = std::atoi(chamber_type_name.substr(14,1).c_str())-1 + offset;/*0..6..8*/
  const bool isEndcap = chamber_type_name[9] == 'E';
  const int phiId =// 0..47 for E, 0..23 for F, A/M1/sector01/phi0/E1 is 0,  A/M1/sector01/phi2/F is also 1
    (chamber_type_name[5]=='4')
    ?( std::atoi(chamber_type_name.substr(8,1).c_str()) + (std::atoi(chamber_type_name.substr(1,2).c_str())/2)*3 ) /*0..23*/
    :( (isEndcap)
       ?( std::atoi(chamber_type_name.substr(8,1).c_str()) + (std::atoi(chamber_type_name.substr(1,2).c_str())-1)*4 ) /*0..47*/
       :( std::atoi(chamber_type_name.substr(8,1).c_str())/2 + (std::atoi(chamber_type_name.substr(1,2).c_str())-1)*2 ) /*0..23*/ );
  const int eNumber = std::atoi(chamber_type_name.substr(11,1).c_str())-1;/*0..4*/
  const bool isStrip = chamber_type_name[15] == 'S';
 
  int sector = isEndcap ? phiId/4+1 : phiId/2+1;
  int phi = isEndcap ? phiId%4 : (phiId%2)*2;
  int efNumber = (isEndcap) ? eNumber : ( (layerNumber<3)?(4):(5)) ; // 0..3 or 4, plus 1 for forward
 
  phimap_index = (!isEndcap || chamber_type_name[5]=='4') ? (phiId*2+1) : (phiId+1);
 
  int station = 0;
  int iLay = 0;
  if(layerNumber<3){ // M1
    station=1;
    iLay = layerNumber+1;
    etamap_index = ((isStrip)?(2):(3)) * efNumber + ((isStrip)?(iLay/2+1):(iLay));
  }else if(layerNumber<5){ // M2
    station=2;
    iLay = layerNumber+1-3;
    etamap_index = 2 * efNumber + iLay + ((isStrip)?(10):(15));
  }else if(layerNumber<7){ // M3
    station=3;
    iLay = layerNumber+1-5;
    etamap_index = 2 * efNumber + iLay + ((isStrip)?(22):(27));
  }else{ // M4
    station=4;
    iLay = layerNumber+1-7;
    etamap_index = iLay + ((isStrip)?(34):(39));
  }
 
  if(layerNumber<7){
    phimap_global_index = iLay + efNumber*3 + phi*18 + (station-1)*72 + (sector-1)*216;
  }else{
    phimap_global_index = 2592 + iLay + phiId*2;
  }
 
  return true;
}