NrpcDigitToNrpcRDO.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
 
#include "NrpcDigitToNrpcRDO.h"
 
#include <algorithm>
#include "GaudiKernel/PhysicalConstants.h"
#include "MuonDigitContainer/RpcDigit.h"
#include "MuonDigitContainer/RpcDigitCollection.h"
#include "MuonReadoutGeometry/RpcReadoutElement.h"
#include "StoreGate/StoreGateSvc.h"
#include "TrigT1RPClogic/ShowData.h"
#include "xAODMuonRDO/NRPCRDOAuxContainer.h"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
namespace Muon{
StatusCode NrpcDigitToNrpcRDO::initialize() {
    ATH_MSG_DEBUG(" in initialize()");
    
    ATH_CHECK(m_idHelperSvc.retrieve());
    ATH_CHECK( m_NrpcContainerKey.initialize() );
 
    ATH_MSG_DEBUG("Tag info filled successfully");
 
    ATH_CHECK(m_cablingKey.initialize());
    ATH_MSG_DEBUG("Cabling info initialized");
 
    ATH_CHECK(m_digitContainerKey.initialize());
    ATH_MSG_VERBOSE("Initialized ReadHandleKey: " << m_digitContainerKey);
    
    for (const std::string& statName : m_convStat){
        m_selectedStations.insert(m_idHelperSvc->rpcIdHelper().stationNameIndex(statName));
    }    
    return StatusCode::SUCCESS;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
StatusCode NrpcDigitToNrpcRDO::execute(const EventContext& ctx) const {
        ATH_MSG_DEBUG("in execute()");
 
        SG::ReadHandle container(m_digitContainerKey, ctx);
        if (!container.isValid()) {
            ATH_MSG_FATAL("Could not find RpcDigitContainer called " << container.name() << " in store " << container.store());
            return StatusCode::FAILURE;
        }
        SG::ReadCondHandle cabling(m_cablingKey, ctx);
        if (!cabling.isValid()) {
            ATH_MSG_FATAL("Could not find MuonNRPC_CablingMap " );
            return StatusCode::FAILURE;
        }
        ATH_MSG_DEBUG("Found MuonNRPC_CablingMap ");

        SG::WriteHandle nrpcRdoData(m_NrpcContainerKey, ctx);
        ATH_CHECK(nrpcRdoData.record(std::make_unique<xAOD::NRPCRDOContainer>(), 
                                     std::make_unique<xAOD::NRPCRDOAuxContainer>()));
 
 
         const IdContext rpcContext = m_idHelperSvc->rpcIdHelper().module_context();
    
        // loop over digit collections
        for (const RpcDigitCollection* rpcCollection : *container) {
            ATH_MSG_DEBUG("RPC Digit -> Pad loop :: digitCollection at " << rpcCollection);
 
            IdentifierHash moduleHash = rpcCollection->identifierHash();
            Identifier moduleId{0};
 
            if (m_idHelperSvc->rpcIdHelper().get_id(moduleHash, moduleId, &rpcContext)) {
                ATH_MSG_WARNING("Failed to translate the "<<moduleHash<<" to a valid identifier");
                continue;
            }
            if (m_selectedStations.size() &&
                !m_selectedStations.count(m_idHelperSvc->stationName(moduleId))) {
                ATH_MSG_DEBUG("Detector element "<<m_idHelperSvc->toString(moduleId)
                             <<" is not considered to be a small gap RPC");
                continue;
            }
 
            // loop over digit 
            for (const RpcDigit* rpcDigit : *rpcCollection) {
                const Identifier channelId = rpcDigit->identify();
 
                ATH_MSG_DEBUG("Convert RPC digit "<<m_idHelperSvc->toString(channelId));
                      
                RpcCablingData translateCache{};
                if (!cabling->convert(channelId, translateCache, rpcDigit->stripSide())) {
                    ATH_MSG_FATAL("Found a non NRPC identifier " << m_idHelperSvc->toString(channelId));
                    return StatusCode::FAILURE;
                }
            
                if (!cabling->getOnlineId(translateCache, msgStream())) {
                    ATH_MSG_ERROR("Offline to Online Id conversion for NRPC chamber." 
                                  << m_idHelperSvc->toString(channelId) << " failed.");
                    return StatusCode::FAILURE;
                }
                const float rdo_time = rpcDigit->time();
 
                const float the_timeoverthr = rpcDigit->ToT();
                uint32_t the_bcid= rdo_time > 0 ? rdo_time / 25. : 0;
 
                xAOD::NRPCRDO* NrpcRdo = nrpcRdoData->push_back(std::make_unique<xAOD::NRPCRDO>());
                NrpcRdo->setBcid(the_bcid);
                NrpcRdo->setTime(rdo_time);
                NrpcRdo->setSubdetector(translateCache.subDetector);
                NrpcRdo->setBoardsector(translateCache.boardSector);
                NrpcRdo->setBoard(translateCache.board);
                NrpcRdo->setChannel(translateCache.channelId);
                NrpcRdo->setTimeoverthr(the_timeoverthr);
           }
        }
        return StatusCode::SUCCESS;
    }
}