Muon::MuonIntersectGeoData Class Reference

#include <MuonIntersectGeoData.h>

Collaboration diagram for Muon::MuonIntersectGeoData:

Public Member Functions
	~MuonIntersectGeoData ()
	MuonIntersectGeoData ()
	MuonIntersectGeoData (MsgStream &msg, const MuonGM::MuonDetectorManager *detMgr, const IMuonIdHelperSvc *idHelperSvc, const MdtCondDbData *cond_data)
Muon::MuonStationIntersect	tubesCrossedByTrack (const MuonGM::MuonDetectorManager *detMgr, const Identifier &id, const Amg::Vector3D &pos, const Amg::Vector3D &dir) const
std::vector< std::shared_ptr< const MdtIntersectGeometry > >	getStationGeometry (const Identifier &id) const
	get geometry description of the given chamber + neighbouring chambers More...
std::shared_ptr< const MdtIntersectGeometry >	getChamber (const Identifier &id) const
	get a pointer to the cached chamber More...

Private Member Functions
std::vector< Identifier >	binPlusneighbours (const Identifier &id) const

Private Attributes
std::vector< std::shared_ptr< MdtIntersectGeometry > >	m_geometry {}
const IMuonIdHelperSvc *	m_idHelperSvc {nullptr}
const MdtCondDbData *	m_dbData {nullptr}

Detailed Description

Definition at line 22 of file MuonIntersectGeoData.h.

Constructor & Destructor Documentation

~MuonIntersectGeoData()

Muon::MuonIntersectGeoData::~MuonIntersectGeoData ( )

default

MuonIntersectGeoData() [1/2]

Muon::MuonIntersectGeoData::MuonIntersectGeoData ( )

default

MuonIntersectGeoData() [2/2]

Muon::MuonIntersectGeoData::MuonIntersectGeoData	(	MsgStream &	msg,
		const MuonGM::MuonDetectorManager *	detMgr,
		const IMuonIdHelperSvc *	idHelperSvc,
		const MdtCondDbData *	cond_data
	)

Definition at line 18 of file MuonIntersectGeoData.cxx.

                                                                                                                 :
        m_idHelperSvc{idHelperSvc}, m_dbData{dbData} {
        m_geometry.resize(m_idHelperSvc->mdtIdHelper().detectorElement_hash_max());
        for (unsigned int n = 0; n < m_geometry.size(); ++n) {
            IdentifierHash id_hash{n};
            const MuonGM::MdtReadoutElement* mdt_ele = detMgr->getMdtReadoutElement(id_hash);
            // The Mdt intersectionGeometry relies on the first multi layer. We can skip the other layers
            if (!mdt_ele || mdt_ele->getMultilayer() != 1) continue;
            m_geometry[id_hash] = std::make_unique<Muon::MdtIntersectGeometry>(log, mdt_ele->identify(), idHelperSvc, detMgr, dbData);
        }    
    } 

Member Function Documentation

binPlusneighbours()

std::vector< Identifier > Muon::MuonIntersectGeoData::binPlusneighbours ( const Identifier &  id ) const

private

Definition at line 68 of file MuonIntersectGeoData.cxx.

                                                                                            {
        std::vector<Identifier> chIds;
        int stName = m_idHelperSvc->mdtIdHelper().stationName(id);
        int stPhi = m_idHelperSvc->mdtIdHelper().stationPhi(id);
        int stEta = m_idHelperSvc->mdtIdHelper().stationEta(id);
        bool isBarrel = m_idHelperSvc->mdtIdHelper().isBarrel(id);
        auto [stEtaMin, stEtaMax] = m_idHelperSvc->mdtIdHelper().stationEtaMinMax(id);
 
        int chEtaLeft = stEta - 1;
        int chEtaRight = stEta + 1;
 
        
 
        // chamber with smallest eta
        if (chEtaLeft < stEtaMin) chEtaLeft = -999;
 
        // chamber with largest eta
        if (chEtaRight > stEtaMax) chEtaRight = -999;
 
        using namespace MuonStationIndex;
        ChIndex chIndex = m_idHelperSvc->chamberIndex(id);
 
        // special treatment of EOS chambers
        if (chIndex == ChIndex::EOS) {
            chEtaRight = -999;
            chEtaLeft = -999;
        }
 
        if (isBarrel) {
            // eta = 0 doesn't exist take next
            if (chEtaLeft == 0) chEtaLeft -= 1;
            if (chEtaRight == 0) chEtaRight += 1;
        } else {
            // do not combined positive and negative endcaps
            if (chEtaLeft == 0) chEtaLeft = -999;
            if (chEtaRight == 0) chEtaRight = -999;
        }
 
        // no neighbours for BIS8
        if (chIndex == ChIndex::BIS && std::abs(stEta) == 8) {
            chEtaLeft = -999;
            chEtaRight = -999;
        }
 
        // BIS 8 never neighbour of a chamber
        if (chIndex == ChIndex::BIS) {
            if (std::abs(chEtaLeft) == 8) chEtaLeft = -999;
            if (std::abs(chEtaRight) == 8) chEtaRight = -999;
        }
        
        if ((chIndex == ChIndex::EIS || chIndex == ChIndex::EIL) && !m_idHelperSvc->hasCSC()) {
            //Chambers can only be -5, -4, 4, 5
            if( std::abs(chEtaLeft) != 5 && std::abs(chEtaLeft) != 4 ) chEtaLeft = -999;
            if( std::abs(chEtaRight) != 5 && std::abs(chEtaRight) != 4 ) chEtaRight = -999;
        }
 
        if (chEtaLeft != -999 &&
            m_idHelperSvc->mdtIdHelper().validElement(m_idHelperSvc->mdtIdHelper().elementID(stName, chEtaLeft, stPhi)))
            chIds.push_back(m_idHelperSvc->mdtIdHelper().elementID(stName, chEtaLeft, stPhi));
        chIds.push_back(m_idHelperSvc->mdtIdHelper().elementID(id));
        if (chEtaRight != -999 &&
            m_idHelperSvc->mdtIdHelper().validElement(m_idHelperSvc->mdtIdHelper().elementID(stName, chEtaRight, stPhi)))
            chIds.push_back(m_idHelperSvc->mdtIdHelper().elementID(stName, chEtaRight, stPhi));
 
        // if (msgLvl(MSG::VERBOSE)) {
        //     ATH_MSG_VERBOSE(" returning chambers for id " << m_idHelperSvc->toString(id) << " ids " << chIds.size() << "   eta " <<
        //     chEtaLeft
        //                                                   << "     " << stEta << "     " << chEtaRight << " chambers: ");
        //     for (unsigned int i = 0; i < chIds.size(); ++i) ATH_MSG_VERBOSE(m_idHelperSvc->toString(chIds[i]));
        // }
 
        return chIds;
    }

getChamber()

std::shared_ptr< const MdtIntersectGeometry > Muon::MuonIntersectGeoData::getChamber

(

const Identifier &

id

)

const

get a pointer to the cached chamber

Definition at line 46 of file MuonIntersectGeoData.cxx.

                                                                                                           {
         IdentifierHash hash{0};
         const MdtIdHelper& idHelper{m_idHelperSvc->mdtIdHelper()};
         if (idHelper.get_detectorElement_hash(idHelper.multilayerID(chId), hash)) return nullptr;
         return hash < m_geometry.size() ? m_geometry[hash] : nullptr;
    }

getStationGeometry()

std::vector< std::shared_ptr< const Muon::MdtIntersectGeometry > > Muon::MuonIntersectGeoData::getStationGeometry

(

const Identifier &

id

)

const

get geometry description of the given chamber + neighbouring chambers

Definition at line 31 of file MuonIntersectGeoData.cxx.

                                                                                                                                {
        // get ids of geometry associated with identifier
        std::vector<Identifier> chambers = binPlusneighbours(id);
 
        // vector to store result
        std::vector<std::shared_ptr<const Muon::MdtIntersectGeometry>> stations;
        stations.reserve(chambers.size());
        // loop over bins, retrieve geometry
        for (const Identifier& chId : chambers) {
            if (m_dbData && !m_dbData->isGood(chId)) { continue; }
            std::shared_ptr<const MdtIntersectGeometry> chamb = getChamber(chId);
            if (chamb) stations.push_back(std::move(chamb));
        }
        return stations;
    }

tubesCrossedByTrack()

Muon::MuonStationIntersect Muon::MuonIntersectGeoData::tubesCrossedByTrack	(	const MuonGM::MuonDetectorManager *	detMgr,
		const Identifier &	id,
		const Amg::Vector3D &	pos,
		const Amg::Vector3D &	dir
	)		const

Definition at line 53 of file MuonIntersectGeoData.cxx.

                                                                                                       {
        std::vector<std::shared_ptr<const Muon::MdtIntersectGeometry>> stations = getStationGeometry(id);
 
        Muon::MuonStationIntersect::TubeIntersects tubeIntersects;
        for (std::shared_ptr<const Muon::MdtIntersectGeometry>& it : stations) {
            Muon::MuonStationIntersect intersect = it->intersection(detMgr, pos, dir);
            tubeIntersects.insert(tubeIntersects.end(), intersect.tubeIntersects().begin(), intersect.tubeIntersects().end());
        }
 
        Muon::MuonStationIntersect intersection{std::move(tubeIntersects)};
        return intersection;
    }

Member Data Documentation

m_dbData

const MdtCondDbData* Muon::MuonIntersectGeoData::m_dbData {nullptr}

private

Definition at line 42 of file MuonIntersectGeoData.h.

m_geometry

std::vector<std::shared_ptr<MdtIntersectGeometry> > Muon::MuonIntersectGeoData::m_geometry {}

private

Definition at line 40 of file MuonIntersectGeoData.h.

m_idHelperSvc

const IMuonIdHelperSvc* Muon::MuonIntersectGeoData::m_idHelperSvc {nullptr}

private

Definition at line 41 of file MuonIntersectGeoData.h.

---

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

• MuonIntersectGeoData.h
• MuonIntersectGeoData.cxx