TGCCableSLBToSSW.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
#include "MuonTGC_Cabling/TGCCableSLBToSSW.h"
 
#include "MuonTGC_Cabling/TGCDatabaseSLBToROD.h" 
#include "MuonTGC_Cabling/TGCModuleSLB.h"
#include "MuonTGC_Cabling/TGCModuleSSW.h"
 
namespace MuonTGC_Cabling {
 
TGCCableSLBToSSW::TGCCableSLBToSSW(const std::string& filename)
  : TGCCable(TGCCable::SLBToSSW)
{
  m_database[TGCId::Endcap][TGCId::WT] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB EWT");
  m_database[TGCId::Endcap][TGCId::WD] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB EWD");
  m_database[TGCId::Endcap][TGCId::ST] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB EST");
  m_database[TGCId::Endcap][TGCId::SD] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB ESD");
  m_database[TGCId::Endcap][TGCId::WI] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB EWI");
  m_database[TGCId::Endcap][TGCId::SI] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB ESI");
  m_database[TGCId::Endcap][SL] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB ESL");
 
  m_database[TGCId::Forward][TGCId::WT] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB FWT");
  m_database[TGCId::Forward][TGCId::WD] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB FWD");
  m_database[TGCId::Forward][TGCId::ST] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB FST");
  m_database[TGCId::Forward][TGCId::SD] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB FSD");
  m_database[TGCId::Forward][TGCId::WI] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB FWI");
  m_database[TGCId::Forward][TGCId::SI] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB FSI");
  m_database[TGCId::Forward][SL] = std::make_unique<TGCDatabaseSLBToROD>(filename, "SLB FSL");
}
 
TGCModuleMap* TGCCableSLBToSSW::getModule(const TGCModuleId* moduleId) const {
  if(moduleId){
    if(moduleId->getModuleIdType()==TGCModuleId::SLB)
      return getModuleOut(moduleId);
    if(moduleId->getModuleIdType()==TGCModuleId::SSW)
      return getModuleIn(moduleId);
  }
  return nullptr;
}
 
TGCModuleMap* TGCCableSLBToSSW::getModuleIn(const TGCModuleId* ssw) const {
  if(ssw->isValid()==false) return nullptr;
 
  const int sswId = ssw->getId();
 
  TGCDatabase* databaseP[TGCId::MaxRegionType * MaxModuleType];
  TGCId::ModuleType module[TGCId::MaxRegionType * MaxModuleType];
  TGCId::RegionType region[TGCId::MaxRegionType * MaxModuleType];
  for(int i=0; i<TGCId::MaxRegionType; i++){
    for(int j=0; j< MaxModuleType; j++){
      databaseP[i*MaxModuleType + j] = m_database[i][j].get();
      region[i*MaxModuleType + j] = static_cast<TGCId::RegionType>(i);
      module[i*MaxModuleType + j] = static_cast<TGCId::ModuleType>(j);
    }
  }
 
  TGCModuleMap* mapId = nullptr;
  for(int type=0; type<TGCId::MaxRegionType* MaxModuleType; type++){
    const int MaxEntry = databaseP[type]->getMaxEntry();
    for(int i=0; i<MaxEntry; i++){
 
      if(databaseP[type]->getEntry(i,4) == sswId) {
    int sector = databaseP[type]->getEntry(i,0);
 
    // Inner
    if(module[type]==TGCId::WI) {
      // EI/FI 
          //  wire(TGCId::WI) and strip(TGCId::SI) of a chamber
          //  use the same SLB chip             
 
      int sswSector = ssw->getReadoutSector();
          // one ROD (rodId = 2,5,8,11) covers 6 Inner sectors
      if(sswSector %3 != 1) continue;
          sswSector -= 1;
 
      sector += sswSector*2;
      int id = databaseP[type]->getEntry(i,1);
      int block = databaseP[type]->getEntry(i,5);
      int sbLoc = databaseP[type]->getEntry(i,2);
          int slbAddr = databaseP[type]->getEntry(i,3);
      for(int ip= 0; ip<3 ; ip++){ 
        int psector  = sector  + ip*2;
        int pblock   = block   + ip*4;
        int psbLoc   = sbLoc   + ip*4;      
            int pslbAddr = slbAddr + ip*2;
        TGCModuleSLB* slb = new TGCModuleSLB(ssw->getSideType(),
                         module[type],
                         region[type],
                         psector,
                         id,
                         psbLoc,
                         pslbAddr);
      
        if(mapId==nullptr) mapId = new TGCModuleMap();
        mapId->insert(pblock,slb);
      }
 
    } else if(module[type]==TGCId::SI) {
      // EI/FI 
          //  wire(TGCId::WI) and strip(TGCId::SI) of a chamber
          //  use the same SLB chip             
 
      // do nothing
          
    } else {
      if(region[type]==TGCId::Endcap) {
        sector += ssw->getReadoutSector() * (TGCId::NUM_ENDCAP_SECTOR/TGCId::N_RODS);
      } else {
        sector += ssw->getReadoutSector() * (TGCId::NUM_FORWARD_SECTOR/TGCId::N_RODS);
      }
      int id = databaseP[type]->getEntry(i,1);
      int sbLoc = databaseP[type]->getEntry(i,2);
          int slbAddr = databaseP[type]->getEntry(i,3);
      int block = databaseP[type]->getEntry(i,5);
      TGCModuleSLB* slb = new TGCModuleSLB(ssw->getSideType(),
                           module[type],
                           region[type],
                           sector,
                           id,
                           sbLoc,
                           slbAddr);
      
      if(mapId==nullptr) mapId = new TGCModuleMap();
      mapId->insert(block,slb);
    } 
 
      }
    }
  }
 
  return mapId;
}
 
 
TGCModuleMap* TGCCableSLBToSSW::getModuleOut(const TGCModuleId* slb) const {
  if(slb->isValid()==false) return nullptr;
 
  const int slbId = slb->getId();
  const int sector = slb->getSectorInReadout();
  int readoutSector = slb->getReadoutSector();
  
  TGCModuleMap* mapId = nullptr;
 
  TGCDatabase* databaseP =m_database[slb->getRegionType()][slb->getModuleType()].get();
  const int MaxEntry = databaseP->getMaxEntry();
  if(slb->getModuleType() == TGCId::WI ||
     slb->getModuleType() == TGCId::SI) {
    // inner 
    for(int i=0; i<MaxEntry; i++){
      if(databaseP->getEntry(i,0)==sector &&
     databaseP->getEntry(i,1)==slbId) {
    
        int sswSectorRO =   readoutSector - (readoutSector%3) + 1;
    int id = databaseP->getEntry(i,4);
    int block = databaseP->getEntry(i,5) + 4*(readoutSector%3);
        
    TGCModuleSSW* ssw = new TGCModuleSSW(slb->getSideType(),
                         sswSectorRO,
                         id);
        mapId = new TGCModuleMap();
        mapId->insert(block,ssw);
    break;
      } 
    }
    
  } else { 
    for(int i=0; i<MaxEntry; i++){
      if(databaseP->getEntry(i,0)==sector&&
     databaseP->getEntry(i,1)==slbId) {
    
    int id = databaseP->getEntry(i,4);
    int block = databaseP->getEntry(i,5);
    TGCModuleSSW* ssw = new TGCModuleSSW(slb->getSideType(),
                         readoutSector,
                         id);
        mapId = new TGCModuleMap();
        mapId->insert(block,ssw);
    break;
      } 
    }
 
  }
  return mapId;
}
  
} //end of namespace