MdtCablingJsonDumpAlg Class Reference

#include <MdtCablingJsonDumpAlg.h>

Inheritance diagram for MdtCablingJsonDumpAlg:

Collaboration diagram for MdtCablingJsonDumpAlg:

Public Types
using	CablingData = MuonMDT_CablingMap::CablingData

Public Member Functions
virtual	~MdtCablingJsonDumpAlg ()=default
virtual StatusCode	initialize () override
virtual StatusCode	execute (const EventContext &ctx) override
	Execute method.
virtual unsigned int	cardinality () const override final
	AthAlgorithm (const std::string &name, ISvcLocator *pSvcLocator)
	Constructor with parameters:
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
bool	filterPassed (const EventContext &ctx) const
void	setFilterPassed (bool state, const EventContext &ctx) const
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const
const EventContext &	getContext () const
	Deprecated methods (use the ones with EventContext).
bool	filterPassed () const
void	setFilterPassed (bool state) const

Protected Member Functions
virtual bool	isReEntrant () const override final
	Legacy algorithms are not thread-safe.
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
SG::ReadCondHandleKey< MuonGM::MuonDetectorManager >	m_DetectorManagerKey
SG::ReadCondHandleKey< MuonMDT_CablingMap >	m_cablingKey {this, "ReadKey", "MuonMDT_CablingMap", "Key of input MDT cabling map"}
Gaudi::Property< std::string >	m_summaryTxt {this, "SummaryFile", "SummaryFile.txt", "Summary of the extracted mapping"}
Gaudi::Property< std::string >	m_mezzJSON {this, "OutMezzanineJSON", "MezzMapping.json", "Mezzanine JSON"}
Gaudi::Property< std::string >	m_cablingJSON {this, "OutCablingJSON", "MdtCabling.json", "Cabling JSON"}
Gaudi::Property< bool >	m_insertBISCabling {this, "insertBISCabling", false, "override the BIS cabling in the JSON file"}
DataObjIDColl	m_extendedExtraObjects
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default).
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 22 of file MdtCablingJsonDumpAlg.h.

Member Typedef Documentation

CablingData

using MdtCablingJsonDumpAlg::CablingData = MuonMDT_CablingMap::CablingData

Definition at line 30 of file MdtCablingJsonDumpAlg.h.

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

~MdtCablingJsonDumpAlg()

virtual MdtCablingJsonDumpAlg::~MdtCablingJsonDumpAlg ( )

virtualdefault

Member Function Documentation

AthAlgorithm()

AthAlgorithm::AthAlgorithm	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

Constructor with parameters:

Definition at line 57 of file AthAlgorithm.cxx.

                                                       :
  ::AthCommonDataStore<AthCommonMsg<Algorithm>>   ( name, pSvcLocator )
{
  // default cardinality for non-reentrant algorithms
  setProperty( "Cardinality", 1 ).orThrow("Unable to set property 'Cardinality'", name);
 
  // Set up to run AthAlgorithmDHUpdate in sysInitialize before
  // merging dependency lists.  This extends the output dependency
  // list with any symlinks implied by inheritance relations.
  m_updateDataHandles =
    std::make_unique<AthenaBaseComps::AthAlgorithmDHUpdate>
    (m_extendedExtraObjects,
     std::move (m_updateDataHandles));
}

cardinality()

virtual unsigned int MdtCablingJsonDumpAlg::cardinality ( ) const

inlinefinaloverridevirtual

Definition at line 28 of file MdtCablingJsonDumpAlg.h.

{return 1;}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl, const SG::VarHandleKeyType &  )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode MdtCablingJsonDumpAlg::execute ( const EventContext & ctx )

overridevirtual

Execute method.

Provides access to the EventContext if needed but is non-const as opposed to AthReentrantAlgorithm.

Struct to perform the mapping between online offline tube

Create the cabling object

Test if the online channel can be found

Next step is to filter all the mezzanine cards

Assign MROD / SubDet / CSM

Offline Identifier (stationIdx, eta,phi, ml)

Tube zero

The staggering of the hedgehog cards does not coincide with the multi-layer staggering i.e. tubeZero is 2,3,4,5 or 2,3,4,5,6,7.

Check whether the same layout is already used somewhere

Subdetector identifiers 97: barrel (A), 98: barrel (C), 99 endcap (A), 100, endcap (C)

Summary file

Write mezzanine file

Implements AthAlgorithm.

Definition at line 45 of file MdtCablingJsonDumpAlg.cxx.

                                                                 {
  ATH_MSG_INFO("Dump cabling & mezzanines into JSON file");
 
  const MuonGM::MuonDetectorManager* detectorMgr{nullptr};
  ATH_CHECK(SG::get(detectorMgr, m_DetectorManagerKey, ctx));
 
  const MuonMDT_CablingMap* cabling{nullptr};
  ATH_CHECK(SG::get(cabling, m_cablingKey,ctx));
 
  const MdtIdHelper& idHelper = m_idHelperSvc->mdtIdHelper();
 
  std::vector<MdtMezzanineCard> cached_cards{};
  std::set<MdtCablingData> cached_chnls{};
  constexpr auto unsetArray{make_array<uint8_t,24>(MdtMezzanineCard::NOTSET)};
  for (auto det_itr = idHelper.detectorElement_begin(); det_itr != idHelper.detectorElement_end(); ++det_itr){
    const MuonGM::MdtReadoutElement* readEle = detectorMgr->getMdtReadoutElement(*det_itr);
    if (!readEle) {
        ATH_MSG_DEBUG("Detector element does not exist. ");
        continue;
    }
    const Identifier station_id = idHelper.elementID(readEle->identify());
    if(m_insertBISCabling && m_idHelperSvc->stationNameString(station_id)=="BIS" && readEle->getStationEta() < 7){
      ATH_MSG_DEBUG("Will skip cabling for station "
                  <<m_idHelperSvc->toStringChamber(readEle->identify()) 
                  << " since it will be inserted manually");
      continue;
    }
  
    MdtMezzanineCard dummy_card(Mapping{}, readEle->getNLayers(), -1);
    std::map<TdcIdentifier, Mapping> chamber_mezz{};
    const int nLayers{readEle->getNLayers()};
    const int nTubes{readEle->getNtubesperlayer()};
    const int multiLayer{readEle->getMultilayer()};
 
    for (int layer = 1 ; layer <= nLayers; ++layer){
      for (int tubeInLayer = 1 ; tubeInLayer <= nTubes; ++tubeInLayer) {
          bool is_valid{false};
          const Identifier tube_id = idHelper.channelID(station_id, multiLayer, 
                                                        layer, tubeInLayer, is_valid);
              if (!is_valid) {
                ATH_MSG_VERBOSE("Invalid element");
                continue;
              }
              MdtCablingData cabling_data{};
              cabling->convert(tube_id,cabling_data);
              if (!cabling->getOnlineId(cabling_data, msgStream())) {
                    ATH_MSG_WARNING("Could no retrieve a valid online channel for "<<m_idHelperSvc->toString(tube_id));
                    continue;
              }
              const TdcIdentifier tdc_id{cabling_data};
              chamber_mezz.try_emplace(tdc_id,unsetArray);
              chamber_mezz[tdc_id][cabling_data.channelId] = dummy_card.tubeNumber(layer, tubeInLayer);
         }
    }
    for (auto&[tdc, mezz_mapping] : chamber_mezz) {
        MdtCablingData mezzCablingId{};
        
        static_cast<MdtCablingOnData&>(mezzCablingId) = tdc;
        mezzCablingId.tdcId = tdc.tdcId;
        
        MdtTdcOnlSorter chipInCab =  cabling->getOnlineConvMap().at(tdc).all_modules.at(tdc.tdcId);
        static_cast<MdtCablingOffData&>(mezzCablingId) = chipInCab->offId();
        mezzCablingId.tube = chipInCab->tubeZero();
        
        const uint8_t tubeOffSet = (mezzCablingId.tube-1)%dummy_card.numTubesPerLayer();
        if (tubeOffSet) {
           MdtMezzanineCard remap{mezz_mapping, dummy_card.numTubeLayers(), 0};
           for (size_t chan = 0 ; chan < mezz_mapping.size(); ++chan) {           
              const OfflineCh tube_lay = remap.offlineTube(chan, msgStream());
              if (!tube_lay.isValid) continue;
              uint8_t tubeNumber = tube_lay.tube + tubeOffSet + 1;
              mezz_mapping[chan] = remap.tubeNumber(tube_lay.layer, tubeNumber);
           }
        }
        Mapping& mtmp = mezz_mapping; // Work around clang15 compilation error.
        std::vector<MdtMezzanineCard>::const_iterator itr = std::find_if(cached_cards.begin(), cached_cards.end(), 
            [&dummy_card, &mtmp](const MdtMezzanineCard& card ){
               if (dummy_card.numTubeLayers() != card.numTubeLayers()) return false;
               for (size_t ch =0; ch < mtmp.size(); ++ch) {
                  if (mtmp[ch] != card.tdcToTubeMap()[ch]) return false;
               }
              return true;
        });
        if (itr != cached_cards.end()) {
          mezzCablingId.mezzanine_type = itr->id();
        } else {
          cached_cards.emplace_back(mezz_mapping, dummy_card.numTubeLayers(), cached_cards.size() + 10);
          if (!cached_cards.back().checkConsistency(msgStream())) {
            ATH_MSG_ERROR("Wrong assignment for "<<mezzCablingId);
            return StatusCode::FAILURE;
          }
          mezzCablingId.mezzanine_type = cached_cards.back().id();
        }
        cached_chnls.insert(std::move(mezzCablingId));
    }
  }
 
  if(m_insertBISCabling){
    ATH_MSG_INFO("Inserting BIS cabling manually");
    // defining the costum mapping to be inserted for BIS chambers
    
    /* The BIS chambers are build using two different tube stagering configurations, one for the A side and one for the C side
    For sectors 12 and 16 the chambers of the A and C side are swapped following space constrains in the detector,
    bringing the readout side of those chambers to the other side. This also swaps the tube staggering configuration.
    Following the staggering configuration two types of mezzanine cards are used to read out the BIS chambers. 
 
    Mezanine type 6 are used for the A-side type chambers. It has the following channel mapping:
    
    23  19  15  11  7   3
      22  18  14  10  6   2
    21  17  13  9   5   1 
      20  16  12  8   4   0
    
    For the C-side type chambers mezzanine type 7 is used, with the following channel mapping:
 
      1   5   9   13  17  21
    0   4   8   12  16  20
      3   7   11  15  19  23
    2   6   10  14  18  22
 
    For the mezz type 7 channel 0 is connected to layer 3  while for mezz type 6 channel 0 is connected to layer 4 so below the mapping will be rotated accordingly. Since the chambers have Nx6 - 2 tubes there are special mezzanine cards with a 4x5 configuration. Since they maintain the above design but missing the last row of tubes (the row on the opposite side of channel 0) so a 6b and 7b design will be introduced with 4 channels set to read tube 255 as proxy for not connected channels. 
    
    The CSM modules of each chamber can read up to 20 mezzanine cards. For BIS 2-7 this is enough to read out all the tube of the chamber.
    For the BIS1 chambers the 24 mezzanine cards are readout by two CSMs each reading out 12 mezzanine cards. In general the channel number of the CSM to which the mezzanine is connected to is even for mezzanine cards connected to multilayer 1 and odd for multilayer 2, starting to count from 0. The channel number is increased in the direction of the global z coordinate axis. 
    */
    constexpr std::array<uint8_t,24> mezzanineType6InvertedMap{23, 17, 11, 5,
                                                               22, 16, 10, 4,
                                                                21, 15, 9, 3,
                                                                20, 14, 8, 2,
                                                                19, 13, 7, 1,
                                                                18, 12, 6, 0};
    
    constexpr std::array<uint8_t,24> mezzanineType6bInvertedMap{250, 250, 250,250,
                                                                22, 16, 10, 4,
                                                                21, 15, 9, 3,
                                                                20, 14, 8, 2,
                                                                19, 13, 7, 1,
                                                                18, 12, 6, 0};
    constexpr std::array<uint8_t,24> mezzanineType6b2InvertedMap{250, 250, 250, 250,
                                                                 23, 17, 11, 5,
                                                                 22, 16, 10, 4,
                                                                 21, 15, 9, 3,
                                                                 20, 14, 8, 2,
                                                                 19, 13, 7, 1};
    
    constexpr std::array<uint8_t,24> mezzanineType7Map{12, 18, 0, 6,
                                                       13, 19, 1, 7,
                                                       14, 20, 2, 8,
                                                       15, 21, 3, 9,
                                                       16, 22, 4, 10,
                                                       17, 23, 5, 11};
    
    constexpr std::array<uint8_t,24> mezzanineType7bMap{12, 18, 0, 6,
                                                       13, 19, 1, 7,
                                                       14, 20, 2, 8,
                                                       15, 21, 3, 9,
                                                       16, 22, 4, 10,
                                                       250, 250, 250, 250};
    
    constexpr std::array<uint8_t,24> mezzanineType7b2Map{ 13, 19, 1, 7,
                                                          14, 20, 2, 8,
                                                          15, 21, 3, 9,
                                                          16, 22, 4, 10,
                                                          17, 23, 5, 11,
                                                          250, 250, 250, 250};
 
 
  
    // Let's add the four types to the cached cards if not already present so we can reference them later
    std::map<std::string, MdtMezzanineCard> mezzTypes{};
    mezzTypes.emplace(std::make_pair("mezzanineType6", MdtMezzanineCard(mezzanineType6InvertedMap, 4, 0)));
    mezzTypes.emplace(std::make_pair("mezzanineType6b", MdtMezzanineCard(mezzanineType6bInvertedMap, 4, 0)));
    mezzTypes.emplace(std::make_pair("mezzanineType6b2", MdtMezzanineCard(mezzanineType6b2InvertedMap, 4, 0)));
    mezzTypes.emplace(std::make_pair("mezzanineType7", MdtMezzanineCard(mezzanineType7Map, 4, 0)));
    mezzTypes.emplace(std::make_pair("mezzanineType7b", MdtMezzanineCard(mezzanineType7bMap, 4, 0)));
    mezzTypes.emplace(std::make_pair("mezzanineType7b2", MdtMezzanineCard(mezzanineType7b2Map, 4, 0)));
    std::map<std::string, int> mezzTypeIds{};
    for(auto& [type,card] : mezzTypes){
      auto itr = std::ranges::find_if(cached_cards,
        [&card](const MdtMezzanineCard& cachedCard){
          if(cachedCard.numTubeLayers() != card.numTubeLayers()) {
              return false;
          }
          for(size_t ch = 0; ch < card.tdcToTubeMap().size(); ++ch){
            if(cachedCard.tdcToTubeMap()[ch] != card.tdcToTubeMap()[ch]) {
                return false;
            }
          }
          return true;
        });
      if(itr != cached_cards.end()){
        mezzTypeIds[type] = itr->id();
      } else {
        cached_cards.emplace_back(card.tdcToTubeMap(), card.numTubeLayers(), cached_cards.size() + 10);
        mezzTypeIds[type] = cached_cards.back().id();
      } 
    }
    
    for(auto& [type,card] : mezzTypes){
      ATH_MSG_VERBOSE("mezzz type" << type << " has id " << mezzTypeIds[type]);
    }
    /* We have to keep track of the MROD and CSM numbers for each partition.
    Since there are no MRODs in phase two we can assign the new BIS stations to fantasy MROD numbers while we wait for the new FELIX based cabling scheme. Since there are up to 51 MRODs per partition in RUN3 the new chambers will be assigned to MRDOs starting from 55 upwards. Each MROD can read up to 6 CSMs so the CSM counter will be increased accordingly when we assign the cabling for each chamber.
    */
    std::map<uint8_t, uint8_t> mrodMap, csmMap;
    mrodMap[97] = 55; csmMap[97] = 0; // A-side barrel
    mrodMap[98] = 55; csmMap[98] = 0; // C-side barrel
 
    for (auto det_itr = idHelper.detectorElement_begin(); det_itr != idHelper.detectorElement_end(); ++det_itr){
      const MuonGM::MdtReadoutElement* readEle = detectorMgr->getMdtReadoutElement(*det_itr);
      if (!readEle) {
          ATH_MSG_DEBUG("Detector element does not exist. ");
          continue;
      }
      // since BIS 1-6 chambers have identical layouts for multilayer 1 and 2 and for BIS78C the second multilayer has more tubes, so we will loop only over the multilayer 2 readout elements here and take care of the multilayer 1 cabling at the same time 
      if(readEle->getMultilayer()!=2)  {
        continue;
      }
      const Identifier station_id = idHelper.elementID(readEle->identify());
      if(  !m_insertBISCabling || m_idHelperSvc->stationNameString(station_id)!="BIS" || readEle->getStationEta() >= 7) {
        continue;
      }
      //find the multilayer 1 readout element to get the number of tubes per layer
      Identifier id_ml1 = idHelper.multilayerID(readEle->identify(),1);
      const MuonGM::MdtReadoutElement* readEle_ml1 = detectorMgr->getMdtReadoutElement(id_ml1);
      int tubesPerLayer_ml1 = readEle_ml1->getNtubesperlayer(); 
 
      ATH_MSG_INFO("Will insert cabling for station "<<m_idHelperSvc->toString(station_id));
      
      const int nTubes{readEle->getNtubesperlayer()};
      uint8_t tdcId{0};
      uint8_t subdetectorId{0};
      if(readEle->getStationEta() > 0){
           subdetectorId = 97; 
      } else{
            subdetectorId = 98;
      }
      for(int tubeInLayer = 1 ; tubeInLayer <= nTubes; 
              tubeInLayer +=   mezzTypes.at("mezzanineType6").numTubesPerLayer()){
        
        if(tubeInLayer==1 && std::abs(readEle->getStationEta()) !=7){tubeInLayer -=1;} // account for the fact that the first mezzanine card only has 5 tubes per layer for BIS 1-6
        MdtCablingData mezzCablingId_ml1{};
        if(tubeInLayer <= tubesPerLayer_ml1){
          mezzCablingId_ml1.stationIndex = m_idHelperSvc->stationName(station_id);
          mezzCablingId_ml1.eta = readEle->getStationEta();
          mezzCablingId_ml1.phi = readEle->getStationPhi();
          mezzCablingId_ml1.multilayer = 1;
          mezzCablingId_ml1.subdetectorId = subdetectorId;
          mezzCablingId_ml1.csm = csmMap[subdetectorId];
          mezzCablingId_ml1.mrod = mrodMap[subdetectorId];
          mezzCablingId_ml1.tdcId = tdcId;
          mezzCablingId_ml1.tube = tubeInLayer;
        }
 
        
        MdtCablingData mezzCablingId_ml2{};
        mezzCablingId_ml2.stationIndex = m_idHelperSvc->stationName(station_id);
        mezzCablingId_ml2.eta = readEle->getStationEta();
        mezzCablingId_ml2.phi = readEle->getStationPhi();
        mezzCablingId_ml2.multilayer = 2;
        mezzCablingId_ml2.subdetectorId = subdetectorId;
        mezzCablingId_ml2.csm = csmMap[subdetectorId];
        mezzCablingId_ml2.mrod = mrodMap[subdetectorId];
        mezzCablingId_ml2.tdcId = tdcId +1;
        mezzCablingId_ml2.tube = tubeInLayer;
 
 
        // A-side type chamber
        std::string mezzEndString = "";
        //The first and the last mezzanine cards of BIS1-6 chambers have a 4x5 configuration
        if(std::abs(readEle->getStationEta()) !=7 &&  tubeInLayer + 6 > nTubes) { mezzEndString = "b"; }
        if(std::abs(readEle->getStationEta()) !=7 &&  tubeInLayer==0  ) { mezzEndString = "b2"; }
        ATH_MSG_VERBOSE("For station "<<m_idHelperSvc->toString(station_id) << " tube " << tubeInLayer<< " assigning mezzanine types for tubeInLayer "<<tubeInLayer <<" mezzEndString "<<mezzEndString);
        if(readEle->getStationEta() > 0 ){
          // swapped chamber in sectors 12 and 16
          if(readEle->getStationPhi() == 6 || readEle->getStationPhi() == 8){
            if(tubeInLayer <= tubesPerLayer_ml1) {
                mezzCablingId_ml1.mezzanine_type = mezzTypeIds["mezzanineType7"+mezzEndString]; // type 7 is used for type C-side chambers
            }
            mezzCablingId_ml2.mezzanine_type = mezzTypeIds["mezzanineType7"+mezzEndString]; // type 7 is used for type C-side chambers
          } else {
            if(tubeInLayer <= tubesPerLayer_ml1) {
               mezzCablingId_ml1.mezzanine_type = mezzTypeIds["mezzanineType6"+mezzEndString];
            }
            mezzCablingId_ml2.mezzanine_type = mezzTypeIds["mezzanineType6"+mezzEndString];
          }
        } else {
          // swapped chamber in sectors 12 and 16
          if(readEle->getStationPhi() == 6 || readEle->getStationPhi() == 8){
            if(tubeInLayer <= tubesPerLayer_ml1){
              mezzCablingId_ml1.mezzanine_type = mezzTypeIds["mezzanineType6"+mezzEndString]; // type 6 is used for type A-side chambers
            }
            mezzCablingId_ml2.mezzanine_type = mezzTypeIds["mezzanineType6"+mezzEndString]; // type 6 is used for type A-side chambers
          } else {
            if(tubeInLayer <= tubesPerLayer_ml1){
               mezzCablingId_ml1.mezzanine_type = mezzTypeIds["mezzanineType7"+mezzEndString];
            }
            mezzCablingId_ml2.mezzanine_type = mezzTypeIds["mezzanineType7"+mezzEndString];
          }
        }
        if(tubeInLayer <= tubesPerLayer_ml1) {
          cached_chnls.insert(std::move(mezzCablingId_ml1));
        }
        cached_chnls.insert(std::move(mezzCablingId_ml2));
        // Update identifier
        tdcId+=2;
        // since BIS1 and BIS78C chambers are read out by two CSMs we need to update the tdcId and CSM numbers accordingly
        if((std::abs(readEle->getStationEta()) == 1  && tdcId >=12) ||
                   (readEle->getStationEta() == -7  && tdcId >=17) ){  
         
          tdcId = 0;
          csmMap[subdetectorId] +=1;
          if(csmMap[subdetectorId]>5){
            csmMap[subdetectorId]=0;
            mrodMap[subdetectorId]+=1;
          }
        }
      }
      // Update CSM and MROD for next chamber
      csmMap[subdetectorId] +=1;
      if(csmMap[subdetectorId]>5){
        csmMap[subdetectorId]=0;
        mrodMap[subdetectorId]+=1;
      }
    }
  }    
    {
        std::ofstream summary{m_summaryTxt};
        if (!summary.good()) {
          ATH_MSG_ERROR("Failed to write "<<m_summaryTxt);
          return StatusCode::FAILURE;
        }
        summary<<"Extracted "<<cached_cards.size()<<" mezzanine card layouts and "
              <<cached_chnls.size()<<" chamber channels. \n\n\n";
        for (const MdtMezzanineCard& card : cached_cards) {
          summary<<card;
          MdtCablingOffData chamb{};
          for (const MdtCablingData& cabling : cached_chnls) {
              if (cabling.mezzanine_type != card.id()) continue;
              if (chamb != cabling) {
                chamb = cabling;
                summary<<std::endl<<" *** "<<idHelper.stationNameString(chamb.stationIndex);             
                summary<<static_cast<int>(std::abs(chamb.eta));
                summary<<(chamb.eta > 0 ? "A" : "C");
                summary<<static_cast<int>(chamb.phi);
                summary<<"M"<<static_cast<int>(chamb.multilayer);
                summary<<" --- tdcs: ";
              }
              summary<<static_cast<int>(cabling.tdcId)<<", ";
          }
          summary<<"\n\n"
                  <<"##############################################################\n";
        }
    }
    {
      std::ofstream mezz_json{m_mezzJSON};
      if (!mezz_json.good()) {
          ATH_MSG_ERROR("Failed to write "<<m_summaryTxt);
          return StatusCode::FAILURE;
      }
      mezz_json<<"["<<std::endl;
      for (size_t i = 0; i < cached_cards.size() ; ++i) {
         const MdtMezzanineCard& card  = cached_cards[i];
         mezz_json<<"     {\n";
         mezz_json<<"       \"mezzId\": "<<static_cast<int>(card.id())<<",\n";
         mezz_json<<"       \"nTubeLayer\": "<<static_cast<int>(card.numTubeLayers())<<",\n";
         mezz_json<<"       \"tdcToTubeMap\": [";
         for (size_t ch = 0 ; ch < card.tdcToTubeMap().size(); ++ch) {
           mezz_json<<static_cast<int>(card.tdcToTubeMap()[ch]);
           if (ch + 1 != card.tdcToTubeMap().size())mezz_json<<",";
         }
         mezz_json<<"]\n";
         mezz_json<<"     }";
         if (i +1 != cached_cards.size()) mezz_json<<",";
         mezz_json<<"\n";
      }
      mezz_json<<"]";
    }
    {
      std::ofstream chamb_json{m_cablingJSON};
      if (!chamb_json.good()) {
         ATH_MSG_FATAL("Failed to write "<<m_cablingJSON);
         return StatusCode::FAILURE;
      }
      chamb_json<<"[\n";
      size_t i =0;
      for (const MdtCablingData& chamb : cached_chnls){ 
        chamb_json<<"    {\n";
        chamb_json<<"     \"station\": \""<<idHelper.stationNameString(chamb.stationIndex)<<"\",\n";
        chamb_json<<"     \"eta\": "<<static_cast<int>(chamb.eta)<<",\n";
        chamb_json<<"     \"phi\": "<<static_cast<int>(chamb.phi)<<",\n";
        chamb_json<<"     \"ml\": "<<static_cast<int>(chamb.multilayer)<<",\n";
        chamb_json<<"     \"subDet\": "<<static_cast<int>(chamb.subdetectorId)<<",\n";
        chamb_json<<"     \"csm\": "<<static_cast<int>(chamb.csm)<<",\n";
        chamb_json<<"     \"mrod\": "<<static_cast<int>(chamb.mrod)<<",\n";
        chamb_json<<"     \"tdcId\": "<<static_cast<int>(chamb.tdcId)<<",\n";
        chamb_json<<"     \"mezzId\": "<<static_cast<int>(chamb.mezzanine_type)<<",\n";
        chamb_json<<"     \"tubeZero\": "<<static_cast<int>(chamb.tube)<<"\n";
        chamb_json<<"    }";
        if (i +1 != cached_chnls.size()) chamb_json<<",";
        chamb_json<<"\n";
        ++i;
      }
      chamb_json<<"]"<<std::endl;
    }
 
    
  return StatusCode::SUCCESS;
} 

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return Algorithm::extraOutputDeps();
}

filterPassed() [1/2]

bool AthAlgorithm::filterPassed ( ) const

inherited

Definition at line 94 of file AthAlgorithm.cxx.

                                      {
  return filterPassed( Gaudi::Hive::currentContext() );
}

filterPassed() [2/2]

bool AthAlgorithm::filterPassed ( const EventContext & ctx ) const

inherited

Definition at line 98 of file AthAlgorithm.cxx.

                                                             {
  return execState( ctx ).filterPassed();
}

getContext()

const EventContext & AthAlgorithm::getContext ( ) const

inherited

Deprecated methods (use the ones with EventContext).

Definition at line 90 of file AthAlgorithm.cxx.

                                                   {
  return Gaudi::Hive::currentContext();
}

initialize()

StatusCode MdtCablingJsonDumpAlg::initialize ( )

overridevirtual

Definition at line 37 of file MdtCablingJsonDumpAlg.cxx.

                                            {
  ATH_CHECK(m_idHelperSvc.retrieve());
  ATH_CHECK(m_DetectorManagerKey.initialize());
  ATH_CHECK(m_cablingKey.initialize());
  return StatusCode::SUCCESS;
} 

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

isReEntrant()

virtual bool AthAlgorithm::isReEntrant ( ) const

inlinefinaloverrideprotectedvirtualinherited

Legacy algorithms are not thread-safe.

Definition at line 111 of file AthAlgorithm.h.

{ return false; }

msg()

MsgStream & AthCommonMsg< Gaudi::Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

setFilterPassed() [1/2]

void AthAlgorithm::setFilterPassed ( bool state ) const

inherited

Definition at line 102 of file AthAlgorithm.cxx.

                                                     {
  setFilterPassed( state, Gaudi::Hive::currentContext() );
}

setFilterPassed() [2/2]

void AthAlgorithm::setFilterPassed ( bool state, const EventContext & ctx ) const

inherited

Definition at line 106 of file AthAlgorithm.cxx.

                                                                              {
  execState( ctx ).setFilterPassed(state);
}

sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc.

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, AthHistogramAlgorithm, and PyAthena::Alg.

Definition at line 66 of file AthAlgorithm.cxx.

                                       {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<Algorithm>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;
}

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

updateVHKA()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_cablingJSON

Gaudi::Property<std::string> MdtCablingJsonDumpAlg::m_cablingJSON {this, "OutCablingJSON", "MdtCabling.json", "Cabling JSON"}

private

Definition at line 42 of file MdtCablingJsonDumpAlg.h.

{this, "OutCablingJSON", "MdtCabling.json", "Cabling JSON"};

m_cablingKey

SG::ReadCondHandleKey<MuonMDT_CablingMap> MdtCablingJsonDumpAlg::m_cablingKey {this, "ReadKey", "MuonMDT_CablingMap", "Key of input MDT cabling map"}

private

Definition at line 38 of file MdtCablingJsonDumpAlg.h.

{this, "ReadKey", "MuonMDT_CablingMap", "Key of input MDT cabling map"};

m_DetectorManagerKey

SG::ReadCondHandleKey<MuonGM::MuonDetectorManager> MdtCablingJsonDumpAlg::m_DetectorManagerKey

private

Initial value:

{this, "DetectorManagerKey", "MuonDetectorManager",
                                                                                "Key of input MuonDetectorManager condition data"}

Definition at line 35 of file MdtCablingJsonDumpAlg.h.

                                                                       {this, "DetectorManagerKey", "MuonDetectorManager",
                                                                                "Key of input MuonDetectorManager condition data"};

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects

privateinherited

Definition at line 114 of file AthAlgorithm.h.

m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> MdtCablingJsonDumpAlg::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

private

Definition at line 33 of file MdtCablingJsonDumpAlg.h.

{this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};

m_insertBISCabling

Gaudi::Property<bool> MdtCablingJsonDumpAlg::m_insertBISCabling {this, "insertBISCabling", false, "override the BIS cabling in the JSON file"}

private

Definition at line 43 of file MdtCablingJsonDumpAlg.h.

{this, "insertBISCabling", false, "override the BIS cabling in the JSON file"};

m_mezzJSON

Gaudi::Property<std::string> MdtCablingJsonDumpAlg::m_mezzJSON {this, "OutMezzanineJSON", "MezzMapping.json", "Mezzanine JSON"}

private

Definition at line 41 of file MdtCablingJsonDumpAlg.h.

{this, "OutMezzanineJSON", "MezzMapping.json", "Mezzanine JSON"};

m_summaryTxt

Gaudi::Property<std::string> MdtCablingJsonDumpAlg::m_summaryTxt {this, "SummaryFile", "SummaryFile.txt", "Summary of the extracted mapping"}

private

Definition at line 40 of file MdtCablingJsonDumpAlg.h.

{this, "SummaryFile", "SummaryFile.txt", "Summary of the extracted mapping"};

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

The documentation for this class was generated from the following files:

• MdtCablingJsonDumpAlg.h
• MdtCablingJsonDumpAlg.cxx