TGCRDOVariables.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
 
#include "MuonPRDTest/TGCRDOVariables.h"
#include "MuonSimData/MuonSimDataCollection.h"
#include "StoreGate/ReadCondHandle.h"
#include "MuonReadoutGeometry/TgcReadoutElement.h"
#include "MuonRDO/TgcRdo.h"
 
using namespace Muon;
namespace MuonPRDTest {
    TGCRDOVariables::TGCRDOVariables(MuonTesterTree& tree, 
                                     const std::string& container_name,
                                     const std::string& cabling_key,
                                      MSG::Level msglvl) :
        PrdTesterModule(tree, "RDO_TGC", msglvl), 
        m_key{container_name}, 
        m_tgcCablingKey{cabling_key} {}
    bool TGCRDOVariables::declare_keys() { 
        return declare_dependency(m_key) &&
               declare_dependency(m_tgcCablingKey); 
    }
 
    bool TGCRDOVariables::fill(const EventContext& ctx) {
        ATH_MSG_DEBUG("do fillTGCRDOVariables()");
        const MuonGM::MuonDetectorManager* MuonDetMgr = getDetMgr(ctx);
        if (!MuonDetMgr) { return false; }
        const TgcRdoContainer* tgcrdoContainer{};       
        if (!SG::get(tgcrdoContainer, m_key, ctx).isSuccess()) {
            return false;
        }
        const Muon::TgcCablingMap* cabling{nullptr};
        if (!SG::get(cabling, m_tgcCablingKey, ctx).isSuccess()) {
            return false;
        }
        ATH_MSG_DEBUG("retrieved TGC rdo Container with size " << tgcrdoContainer->size());
 
        if (tgcrdoContainer->size() == 0) ATH_MSG_DEBUG(" TGC rdo Container empty ");
        unsigned int n_rdo{0};
        for (const TgcRdo* coll : *tgcrdoContainer) {
            ATH_MSG_DEBUG("processing collection with size " << coll->size());
            for (const TgcRawData* rdo: *coll) {
 
                bool orFlag = cabling->isOredChannel(rdo->subDetectorId(),
                                                          rdo->rodId(),
                                                          rdo->sswId(),
                                                          rdo->slbId(),
                                                          rdo->bitpos());
 
                Identifier Id;
                bool e_found = cabling->getElementIDfromReadoutID(Id,
                                             rdo->subDetectorId(),
                                             rdo->rodId(),
                                             rdo->sswId(),
                                             rdo->slbId(),
                                             rdo->bitpos(), 
                                             orFlag);
                if (!e_found) {
                    ATH_MSG_ERROR("could not find Identifier from TgcRawData");
                    return false;
                }
 
                const MuonGM::TgcReadoutElement* rdoEl = MuonDetMgr->getTgcReadoutElement(Id);
                if (!rdoEl) {
                   ATH_MSG_ERROR("The TGC hit "<<idHelperSvc()->toString(Id)<<" does not have a detector element attached. That should actually never happen");
                   return false;
                }
 
                m_TGC_rdo_id.push_back(Id);
                Amg::Vector2D localStripPos(0.,0.);
                if ( rdoEl->stripPosition(Id,localStripPos) )  {
                    m_TGC_rdo_localPosX.push_back(localStripPos.x());
                    m_TGC_rdo_localPosY.push_back(localStripPos.y());
                ATH_MSG_DEBUG("TGC RDO: local pos.:  x=" << localStripPos[0] << ",  y=" << localStripPos[1]);
                } else { 
                    ATH_MSG_WARNING("TGC RDO: local Strip position not defined"); 
                }
      
                // asking the detector element to transform this local to the global position
                Amg::Vector3D globalStripPos(0., 0., 0.);
                rdoEl->surface(Id).localToGlobal(localStripPos,Amg::Vector3D(0.,0.,0.),globalStripPos);
                m_TGC_rdo_globalPos.push_back(globalStripPos);
 
                ++n_rdo;
            }
        }
        m_TGC_nRDO = n_rdo;
        ATH_MSG_DEBUG(" finished fillTGCRDOVariables()");
        return true;
    }
}