contain application tools for (muon) tracking geometry More...

#include <MuonTGMeasurementTool.h>

Inheritance diagram for Muon::MuonTGMeasurementTool:

Collaboration diagram for Muon::MuonTGMeasurementTool:

Public Types
typedef Trk::Track	Track
typedef std::pair< const Trk::Layer , std::vector< Identifier > >	PairOfLayerID
typedef std::pair< const Trk::Layer , std::vector< const Trk::PrepRawData > * >	PairOfLayerPrd

Public Member Functions
	MuonTGMeasurementTool (const std::string &type, const std::string &name, const IInterface *)
	Constructor with AlgTool parameters.
virtual	~MuonTGMeasurementTool ()=default
virtual StatusCode	initialize () override
const Trk::TrackParameters *	layerToDetEl (const Trk::Layer , const Trk::TrackParameters , Identifier) const override
const Trk::TrackParameters *	detElToLayer (const Trk::Layer , const Trk::TrackParameters , Identifier) const override
const Trk::RIO_OnTrack *	measToLayer (const Trk::Layer , const Trk::TrackParameters , const Trk::RIO_OnTrack *) const override
double	residual (const Trk::Layer , const Trk::TrackParameters , const Trk::RIO_OnTrack *) const override
double	residual (const Trk::Layer , const Trk::TrackParameters , Identifier) const override
double	residual (const Trk::TrackParameters , const Trk::RIO_OnTrack ) const override
double	residual (const Trk::TrackParameters *, Identifier &) const override
const Identifier	nearestDetEl (const Trk::Layer , const Trk::TrackParameters , bool measPhi, double &pitch) const override
const Trk::Layer *	associatedLayer (Identifier id, Amg::Vector3D &gp) const override
const Trk::Layer *	associatedLayer (Identifier id, const Trk::TrackingVolume *vol) const override
const Trk::Layer *	match (Identifier id, const Trk::Layer *lay) const override
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	sysInitialize () override
	Perform system initialization for an algorithm.
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

Static Public Member Functions
static const InterfaceID &	interfaceID ()
	AlgTool interface method.

Protected Member Functions
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
const Trk::TrackingGeometry *	getGeometry () const
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
ServiceHandle< Trk::ITrackingGeometrySvc >	m_trackingGeometrySvc
SG::ReadCondHandleKey< Trk::TrackingGeometry >	m_trackingGeometryReadKey
const MuonGM::MuonDetectorManager *	m_muonDetMgr
std::unique_ptr< AmgMatrix(5, 5)>	m_tgcProjEta
std::unique_ptr< AmgMatrix(5, 5)>	m_tgcProjPhi
std::unique_ptr< AmgMatrix(5, 5)>	m_rpcProjEta
std::unique_ptr< AmgMatrix(5, 5)>	m_rpcProjPhi
Gaudi::Property< bool >	m_alignedMode {this, "AlignedMode", true}
Gaudi::Property< bool >	m_useDSManager {this, "UseDSManager", false}
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

contain application tools for (muon) tracking geometry

Author: Sarka.nosp@m..Tod.nosp@m.orova.nosp@m.@cer.nosp@m.n.ch

Definition at line 49 of file MuonTGMeasurementTool.h.

Member Typedef Documentation

◆ PairOfLayerID

typedef std::pair<const Trk::Layer*, std::vector<Identifier>*> Muon::MuonTGMeasurementTool::PairOfLayerID

Definition at line 52 of file MuonTGMeasurementTool.h.

◆ PairOfLayerPrd

typedef std::pair<const Trk::Layer*, std::vector<const Trk::PrepRawData*>*> Muon::MuonTGMeasurementTool::PairOfLayerPrd

Definition at line 53 of file MuonTGMeasurementTool.h.

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

◆ Track

typedef Trk::Track Muon::MuonTGMeasurementTool::Track

Definition at line 51 of file MuonTGMeasurementTool.h.

Constructor & Destructor Documentation

◆ MuonTGMeasurementTool()

Muon::MuonTGMeasurementTool::MuonTGMeasurementTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent )

Constructor with AlgTool parameters.

Definition at line 34 of file MuonTGMeasurementTool.cxx.

                                                                                                                       :
    AthAlgTool(type, name, parent), m_muonDetMgr(nullptr) {
    declareInterface<Muon::IMuonTGMeasTool>(this);
}

◆ ~MuonTGMeasurementTool()

virtual Muon::MuonTGMeasurementTool::~MuonTGMeasurementTool ( )

virtualdefault

Member Function Documentation

◆ associatedLayer() [1/2]

const Trk::Layer * Muon::MuonTGMeasurementTool::associatedLayer	(	Identifier	id,
		Amg::Vector3D &	gp ) const

overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1112 of file MuonTGMeasurementTool.cxx.

                                                                                               {
    // Get the messaging service, print where you are
    ATH_MSG_DEBUG("MuonTGMeasurementTool::associatedLayer");
    const Trk::Layer* lay = nullptr;
    // check input
    if (!id.get_identifier32().get_compact()) return lay;
 
    // rely on having misalignment uncertainty covered by span safety marge ( don't loose station from static volume
    //  when misaligned
    const Trk::TrackingVolume* staticVol = getGeometry()->lowestStaticTrackingVolume(gp);
    const Trk::DetachedTrackingVolume* station = nullptr;
    if (staticVol && !staticVol->confinedDetachedVolumes().empty()) {
        Trk::ArraySpan<const Trk::DetachedTrackingVolume* const>  detTV = staticVol->confinedDetachedVolumes();
        for (const auto *i : detTV) {
            if (i->layerRepresentation() && i->layerRepresentation()->layerType() > 0) {
                Identifier stId(i->layerRepresentation()->layerType());
                if (m_idHelperSvc->mdtIdHelper().stationName(stId) == m_idHelperSvc->mdtIdHelper().stationName(id) &&
                    m_idHelperSvc->mdtIdHelper().stationEta(stId) == m_idHelperSvc->mdtIdHelper().stationEta(id) &&
                    m_idHelperSvc->mdtIdHelper().stationPhi(stId) == m_idHelperSvc->mdtIdHelper().stationPhi(id)) {
                    station = i;
                    break;
                }
            }
        }
    }
 
    if (station) lay = associatedLayer(id, station->trackingVolume());
 
    return lay;
}

◆ associatedLayer() [2/2]

const Trk::Layer * Muon::MuonTGMeasurementTool::associatedLayer	(	Identifier	id,
		const Trk::TrackingVolume *	vol ) const

overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1143 of file MuonTGMeasurementTool.cxx.

                                                                                                            {
    const Trk::Layer* lay = nullptr;
    // check input
    if (!vol) return lay;
 
    if (vol->confinedVolumes()) {
        std::span<Trk::TrackingVolume const * const > subVols = vol->confinedVolumes()->arrayObjects();
        std::span<Trk::TrackingVolume const * const >::iterator iter = subVols.begin();
        while (!lay && iter != subVols.end()) {
            lay = associatedLayer(id, *iter);
            if (lay) break;
            ++iter;
        }
        if (lay) return lay;
    }
 
    if (vol->confinedLayers()) {
        std::span<Trk::Layer const * const > ordLay = vol->confinedLayers()->arrayObjects();
        std::span<Trk::Layer const * const >::iterator iter = ordLay.begin();
        while (!lay && iter != ordLay.end()) {
            lay = match(id, *iter);
            if (lay) break;
            ++iter;
        }
        if (lay) return lay;
    }
 
    if (!vol->confinedArbitraryLayers().empty()) {
        Trk::ArraySpan<const Trk::Layer* const> unOrdLay = vol->confinedArbitraryLayers();
        Trk::ArraySpan<const Trk::Layer* const>::iterator iter = unOrdLay.begin();
        while (!lay && iter != unOrdLay.end()) {
            lay = match(id, *iter);
            if (lay) break;
            ++iter;
        }
        if (lay) return lay;
    }
 
    return lay;
}

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::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< AlgTool > >::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());
  }

◆ detElToLayer()

const Trk::TrackParameters * Muon::MuonTGMeasurementTool::detElToLayer	(	const Trk::Layer *	lay,
		const Trk::TrackParameters *	parm,
		Identifier	id ) const

overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 366 of file MuonTGMeasurementTool.cxx.

                                                                                           {
    // Get the messaging service, print where you are
    ATH_MSG_DEBUG("MuonTGMeasurementTool::detElToLayer");
    const Trk::TrackParameters* projPar = nullptr;
 
    // check input
    if (!lay || !parm || !(id.get_identifier32().get_compact() > 0)) return projPar;
 
    // check compatibility of layer info and required id ? this was already done when associating !
    if (!lay->layerType()) return projPar;
    Identifier layId(lay->layerType());
 
    ATH_MSG_DEBUG("MuonTGMeasurementTool::input ok");
 
    unsigned int hitType = 0;
    if (m_idHelperSvc->isMdt(id)) hitType = 1;
    if (m_idHelperSvc->isRpc(id)) hitType = 2;
    if (m_idHelperSvc->isCsc(id)) hitType = 3;
    if (m_idHelperSvc->isTgc(id)) hitType = 4;
 
    unsigned int layType = 0;
    if (m_idHelperSvc->isMdt(layId)) layType = 1;
    if (m_idHelperSvc->isRpc(layId)) layType = 2;
    if (m_idHelperSvc->isCsc(layId)) layType = 3;
    if (m_idHelperSvc->isTgc(layId)) layType = 4;
 
    if (layType != hitType) return projPar;
 
    if (hitType == 0) {
        ATH_MSG_DEBUG("unknown hit technology");
        return projPar;
    }
 
    if (hitType == 1) {
        // local position of the tube
        int tube = m_idHelperSvc->mdtIdHelper().tube(id);
        Amg::Vector2D locWire(0., lay->getRef() + (tube - 1) * 30.035);
        if (fabs(lay->getRef()) > 10e6) {
            double sign = (lay->getRef() > 0) ? 1. : -1.;
            int dec = int(lay->getRef() / 1.e5);
            double ref0 = dec / 1.e5;
            double ref1 = lay->getRef() - dec * 1e5 - 0.5 * (sign + 1) * 1e5;
            locWire[0] = ref0;
            locWire[1] = ref1;
            int tube = m_idHelperSvc->mdtIdHelper().tube(id);
            int tubeMax = m_idHelperSvc->mdtIdHelper().tubeMax(id);
            if (tube > 6 && tubeMax - tube > 5) locWire[0] = 0.;
        }
        if (sqrt(locWire[0] * locWire[0] + locWire[1] * locWire[1]) > 2000.) {
            ATH_MSG_WARNING(name() << " wire shift out bounds for MDT tube :" << m_idHelperSvc->mdtIdHelper().stationName(id) << ","
                                   << m_idHelperSvc->mdtIdHelper().stationEta(id) << "," << m_idHelperSvc->mdtIdHelper().stationPhi(id)
                                   << ": abandon projection");
            return projPar;
        }
        // local position (tube)
        // const Amg::Vector2D locPos = parm->localPosition();
        // direction at the layer
        Amg::Vector3D dir = parm->momentum().unit();
        Amg::Vector3D locDir = (lay->surfaceRepresentation().transform().inverse()) * dir;
        Amg::Vector3D normal = lay->surfaceRepresentation().normal();
        Amg::Vector3D wireDir(1., 0., 0.);
        const Trk::PlaneSurface* laySurf = dynamic_cast<const Trk::PlaneSurface*>(&(lay->surfaceRepresentation()));
        if (!laySurf) return projPar;
        //
        double ND = dir.dot(normal);
        double DL = locDir.dot(wireDir);
        double A = sqrt(1. - DL * DL);
        AmgMatrix(5, 5) mdtProj;
        mdtProj.setIdentity();
        //
        double ori = (laySurf->transform().inverse() * parm->associatedSurface().transform()).rotation()(1, 1) > 0. ? -1. : 1.;
        if (A > 10e-6) {
            mdtProj(0, 1) = ND / A;
            mdtProj(1, 0) = ori;
            mdtProj(1, 1) = DL / A * sqrt(A * A - ND * ND) / A;
            mdtProj(0, 0) = 0.;
        } else {
            mdtProj(0, 1) = 1.;
            mdtProj(1, 0) = ori;
            mdtProj(1, 1) = 0.;
            mdtProj(0, 0) = 0.;
        }
        Amg::VectorX locPar = parm->parameters();
        AmgMatrix(5, 5) mdtProjInv = mdtProj.inverse();
        Amg::VectorX pPar = mdtProjInv * locPar;
        pPar[0] += locWire[0];
        pPar[1] += locWire[1];
        ATH_MSG_DEBUG("back projected parameters(MDT->layer):" << pPar);
 
        std::optional<AmgMatrix(5, 5)> projEM = parm->covariance()->similarity(mdtProjInv);
 
        ATH_MSG_DEBUG("back projected covariance(MDT->layer):" << (*projEM));
        projPar = new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *laySurf, projEM);
        // delete projEM;
    }
 
    if (hitType == 2) {
        // local position at detEl
        // const Amg::Vector2D locPos = parm->localPosition();
        const Trk::PlaneSurface* laySurf = dynamic_cast<const Trk::PlaneSurface*>(&(lay->surfaceRepresentation()));
        if (!laySurf) return projPar;
        //
        Amg::VectorX locPar = parm->parameters();
        AmgMatrix(5, 5)* pMx = nullptr;
        if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
            pMx = m_rpcProjPhi.get();
        else
            pMx = m_rpcProjEta.get();
 
        double eta = 1.;
        double sign = (m_idHelperSvc->rpcIdHelper().measuresPhi(id) && m_idHelperSvc->rpcIdHelper().doubletPhi(id) == 2) ? -1. : 1.;
 
        double ref = (lay->getRef() > 10000.) ? lay->getRef() - 20000. : lay->getRef();
        double zswap = (lay->getRef() > 10000.) ? -1. : 1.;
 
        if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
            locPar[0] *= eta;
        else
            locPar[1] *= -eta;
        if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
            locPar[1] *= zswap;
        else
            locPar[0] *= zswap;
 
        Amg::VectorX pPar = (*pMx) * locPar;
        pPar[0] += sign * ref;
 
        //
        ATH_MSG_DEBUG("back projected parameters(RPC->layer):" << m_idHelperSvc->rpcIdHelper().measuresPhi(id) << "," << pPar);
 
        std::optional<AmgMatrix(5, 5)> projEM = parm->covariance()->similarityT(*pMx);
 
        ATH_MSG_DEBUG("back projected covariance(RPC->layer):" << (*projEM));
        projPar = new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *laySurf, projEM);
        // delete projEM;
 
        if ((parm->position() - projPar->position()).mag() > 0.001) {
            ATH_MSG_DEBUG("geometrical RPC projection (detElToLayer) for hit : "
                          << m_idHelperSvc->rpcIdHelper().measuresPhi(id) << "," << m_idHelperSvc->rpcIdHelper().stationName(id) << ","
                          << m_idHelperSvc->rpcIdHelper().stationEta(id) << "," << m_idHelperSvc->rpcIdHelper().stationPhi(id) << ","
                          << m_idHelperSvc->rpcIdHelper().doubletPhi(id) << "," << m_idHelperSvc->rpcIdHelper().doubletR(id) << ","
                          << m_idHelperSvc->rpcIdHelper().doubletZ(id));
            Amg::Vector2D locPos;
            const Amg::Vector3D& globPos = parm->position();
            bool onSurface = laySurf->globalToLocal(globPos, globPos, locPos);
            if (onSurface) {
                pPar[0] = locPos[0];
                pPar[1] = locPos[1];
                delete projPar;
                projPar = new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *laySurf, projEM);
                // delete projEM;
            } else {
                delete projPar;
                return nullptr;
            }
        }
    }
 
    if (hitType == 3) {
        // local position of detEl
        // const Amg::Vector2D locPos = parm->localPosition();
        const Trk::PlaneSurface* laySurf = dynamic_cast<const Trk::PlaneSurface*>(&(lay->surfaceRepresentation()));
        if (!laySurf) return projPar;
        // dealing with parallel displaced planes
        double diff = m_idHelperSvc->cscIdHelper().measuresPhi(id) ? 1.55 : -1.55;
        //
        Amg::Vector3D layNormal = lay->surfaceRepresentation().normal();
        double DN = parm->momentum().dot(layNormal);
        double t = diff / DN;
        // displacement of planes
        const Amg::Vector2D csc_shift(0., lay->getRef());
        // projection : take into account possible misalignment ;
        AmgMatrix(5, 5)* pMx = nullptr;
        if (m_idHelperSvc->cscIdHelper().measuresPhi(id))
            pMx = m_tgcProjPhi.get();
        else
            pMx = m_tgcProjEta.get();
        AmgMatrix(5, 5) pMxInv = pMx->inverse();
        Amg::VectorX parProj = pMxInv * parm->parameters();
        Amg::Vector3D corrLocPos = lay->surfaceRepresentation().center() - t * parm->momentum() + t * DN * layNormal;
        Amg::Vector2D locCorrLay;
        bool onSurface = lay->surfaceRepresentation().globalToLocal(corrLocPos, corrLocPos, locCorrLay);
        if (onSurface && locCorrLay.size() > 0) {
            parProj[0] += locCorrLay[Trk::locX] + csc_shift[Trk::locX];
            parProj[1] += locCorrLay[Trk::locY] + csc_shift[Trk::locY];
            ATH_MSG_DEBUG("back projected parameters(CSC->layer):" << m_idHelperSvc->cscIdHelper().measuresPhi(id) << "," << parProj);
        }
        std::optional<AmgMatrix(5, 5)> projEM = parm->covariance()->similarity(pMxInv);
 
        ATH_MSG_DEBUG("back projected covariance(CSC->layer):" << (*projEM));
        projPar = new Trk::AtaPlane(parProj[0], parProj[1], parProj[2], parProj[3], parProj[4], *laySurf, projEM);
 
        ATH_MSG_DEBUG("test CSC projection:detElToLayer:" << parm->position() << "," << projPar->position());
    }
 
    if (hitType == 4) {
        const Trk::PlaneSurface* laySurf = dynamic_cast<const Trk::PlaneSurface*>(&(lay->surfaceRepresentation()));
        if (!laySurf) { return projPar; }
 
        AmgMatrix(5, 5)* pMx = nullptr;
        if (m_idHelperSvc->tgcIdHelper().isStrip(id))
            pMx = m_tgcProjPhi.get();
        else
            pMx = m_tgcProjEta.get();
        AmgMatrix(5, 5) pMxInv = pMx->inverse();
        Amg::VectorX locPar = pMxInv * parm->parameters();
        ATH_MSG_DEBUG("back projected parameters(TGC->layer):" << m_idHelperSvc->tgcIdHelper().isStrip(id) << "," << locPar);
 
        std::optional<AmgMatrix(5, 5)> projEM = parm->covariance()->similarity(pMxInv);
 
        ATH_MSG_DEBUG("back projected covariance(TGC->layer):" << (*projEM));
        projPar = new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *laySurf, projEM);
 
        // verify
        if (m_alignedMode && (parm->position() - projPar->position()).mag() > 0.001) {
            ATH_MSG_DEBUG("geometrical TGC projection ( detEl2Layer ):"
                          << m_idHelperSvc->tgcIdHelper().stationName(id) << "," << m_idHelperSvc->tgcIdHelper().stationEta(id) << ","
                          << m_idHelperSvc->tgcIdHelper().stationPhi(id) << "," << m_idHelperSvc->tgcIdHelper().isStrip(id));
            Amg::Vector2D locPos;
            const Amg::Vector3D& globPos = parm->position();
            bool onSurface = laySurf->globalToLocal(globPos, globPos, locPos);
            if (onSurface) {
                locPar[0] = locPos[0];
                locPar[1] = locPos[1];
                delete projPar;
                projPar = nullptr;
                projPar = new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *laySurf, projEM);
            } else {
                delete projPar;
                return nullptr;
            }
        }
    }
    return projPar;
}

◆ detStore()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

95{ return m_detStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

85{ return m_evtStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::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

◆ getGeometry()

const Trk::TrackingGeometry * Muon::MuonTGMeasurementTool::getGeometry ( ) const

inlineprivate

Good old way of retrieving the volume via the geometry service

Definition at line 98 of file MuonTGMeasurementTool.h.

                                                              {
            if (m_trackingGeometryReadKey.empty()) { return m_trackingGeometrySvc->trackingGeometry(); }
            SG::ReadCondHandle<Trk::TrackingGeometry> handle(m_trackingGeometryReadKey, Gaudi::Hive::currentContext());
            if (!handle.isValid()) {
                ATH_MSG_WARNING("Could not retrieve a valid tracking geometry");
                return nullptr;
            }
            return handle.cptr();
        }

◆ initialize()

StatusCode Muon::MuonTGMeasurementTool::initialize ( )

overridevirtual

Definition at line 40 of file MuonTGMeasurementTool.cxx.

                                                 {
    // Get the messaging service, print where you are
    ATH_MSG_INFO("MuonTGMeasurementTool::initialize()");
 
    ATH_CHECK(m_idHelperSvc.retrieve());
    ATH_CHECK(m_DetectorManagerKey.initialize(!m_useDSManager));
    if (m_useDSManager) { ATH_CHECK(detStore()->retrieve(m_muonDetMgr)); }
 
    // define projection matrices
    m_tgcProjEta = std::make_unique<AmgMatrix(5, 5)>();
    m_tgcProjEta->setIdentity();
    (*m_tgcProjEta)(0, 0) = 0.;
    (*m_tgcProjEta)(1, 1) = 0.;
    (*m_tgcProjEta)(0, 1) = 1.;
    (*m_tgcProjEta)(1, 0) = -1.;
    m_tgcProjPhi = std::make_unique<AmgMatrix(5, 5)>();
    m_tgcProjPhi->setIdentity();
 
    m_rpcProjEta = std::make_unique<AmgMatrix(5, 5)>();
    m_rpcProjEta->setIdentity();
    (*m_rpcProjEta)(0, 0) = 0.;
    (*m_rpcProjEta)(1, 1) = 0.;
    (*m_rpcProjEta)(0, 1) = 1.;
    (*m_rpcProjEta)(1, 0) = 1.;
    m_rpcProjPhi = std::make_unique<AmgMatrix(5, 5)>();
    m_rpcProjPhi->setIdentity();
 
    if (!m_trackingGeometryReadKey.empty()) {
        ATH_CHECK(m_trackingGeometryReadKey.initialize());
    } else {
        ATH_CHECK(m_trackingGeometrySvc.retrieve());
    }
    return StatusCode::SUCCESS;
}

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::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.

◆ interfaceID()

const InterfaceID & Muon::IMuonTGMeasTool::interfaceID ( )

inlinestaticinherited

AlgTool interface method.

Definition at line 29 of file IMuonTGMeasTool.h.

                                                {
            static const InterfaceID IID_IMuonTGMeasTool("IMuonTGMeasTool", 1, 0);
            return IID_IMuonTGMeasTool;
        };

◆ layerToDetEl()

const Trk::TrackParameters * Muon::MuonTGMeasurementTool::layerToDetEl	(	const Trk::Layer *	lay,
		const Trk::TrackParameters *	parm,
		Identifier	id ) const

overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 74 of file MuonTGMeasurementTool.cxx.

                                                                                           {
    const MuonGM::MuonDetectorManager* MuonDetMgr = m_muonDetMgr;
    if (!m_useDSManager) {
        SG::ReadCondHandle<MuonGM::MuonDetectorManager> DetectorManagerHandle{m_DetectorManagerKey};
        MuonDetMgr = DetectorManagerHandle.cptr();
        if (MuonDetMgr == nullptr) {
            ATH_MSG_ERROR("Null pointer to the read MuonDetectorManager conditions object");
            // return StatusCode::FAILURE;
        }
    }
 
    // Get the messaging service, print where you are
    ATH_MSG_DEBUG("MuonTGMeasurementTool::layerToDetEl");
    const Trk::TrackParameters* projPar = nullptr;
    // const Amg::Vector2D* locPos = 0;
    // check input
    if (!lay || !parm || !id.get_identifier32().get_compact()) return projPar;
 
    // get tracking geometry
 
    // check compatibility of layer info and required id ? this was already done when associating !
    if (!lay->layerType()) return projPar;
    Identifier layId(lay->layerType());
 
    unsigned int hitType = 0;
    if (m_idHelperSvc->isMdt(id)) hitType = 1;
    if (m_idHelperSvc->isRpc(id)) hitType = 2;
    if (m_idHelperSvc->isCsc(id)) hitType = 3;
    if (m_idHelperSvc->isTgc(id)) hitType = 4;
 
    unsigned int layType = 0;
    if (m_idHelperSvc->isMdt(layId)) layType = 1;
    if (m_idHelperSvc->isRpc(layId)) layType = 2;
    if (m_idHelperSvc->isCsc(layId)) layType = 3;
    if (m_idHelperSvc->isTgc(layId)) layType = 4;
 
    if (layType != hitType) return projPar;
 
    if (hitType == 0) {
        ATH_MSG_ERROR("unknown hit technology");
        return projPar;
    } else {
        ATH_MSG_DEBUG("hit technology:" << hitType);
    }
 
    ATH_MSG_DEBUG("extrapolated covariance:" << parm->covariance());
 
    if (hitType == 1) {
        const MuonGM::MdtReadoutElement* mdtROE = MuonDetMgr->getMdtReadoutElement(id);
        if (!mdtROE) {
            ATH_MSG_WARNING(name() << "MDT readout element not found");
            return projPar;
        }
        // local position of tube
        int tube = m_idHelperSvc->mdtIdHelper().tube(id);
        Amg::Vector2D locWire(0., lay->getRef() + (tube - 1) * 30.035);
        if (fabs(lay->getRef()) > 10e6) {
            double sign = (lay->getRef() > 0) ? 1. : -1.;
            int dec = int(lay->getRef() / 1.e5);
            double ref0 = dec / 1.e5;
            double ref1 = lay->getRef() - dec * 1e5 - 0.5 * (sign + 1) * 1e5;
            locWire[0] = ref0;
            locWire[1] = ref1;
            int tube = m_idHelperSvc->mdtIdHelper().tube(id);
            int tubeMax = m_idHelperSvc->mdtIdHelper().tubeMax(id);
            if (tube > 6 && tubeMax - tube > 5) locWire[0] = 0.;
        }
        if (sqrt(locWire[0] * locWire[0] + locWire[1] * locWire[1]) > 2000.) {
            ATH_MSG_WARNING(name() << " wire shift out bounds for MDT tube :" << m_idHelperSvc->mdtIdHelper().stationName(id) << ","
                                   << m_idHelperSvc->mdtIdHelper().stationEta(id) << "," << m_idHelperSvc->mdtIdHelper().stationPhi(id)
                                   << ": abandon projection");
            return projPar;
        }
        // direction at the layer
        Amg::Vector3D dir = parm->momentum().unit();
        Amg::Vector3D locDir = (lay->surfaceRepresentation().transform().inverse()) * dir;
        Amg::Vector3D wireDir(1., 0., 0.);
        //
        double ND = dir.dot(lay->surfaceRepresentation().normal());
        double DL = locDir.dot(wireDir);
        double A = sqrt(1. - DL * DL);
        AmgMatrix(5, 5) mdtProj;
        mdtProj.setIdentity();
        //
        const Trk::StraightLineSurface* tubeSurf = dynamic_cast<const Trk::StraightLineSurface*>(&(mdtROE->surface(id)));
        double ori = (lay->surfaceRepresentation().transform().inverse() * tubeSurf->transform()).rotation()(1, 1) > 0. ? -1. : 1.;
        if (A > 10e-6) {
            mdtProj(0, 1) = ND / A;
            mdtProj(1, 0) = ori;
            mdtProj(1, 1) = DL / A * sqrt(A * A - ND * ND) / A;
            mdtProj(0, 0) = 0.;
        } else {
            mdtProj(0, 1) = 1.;
            mdtProj(1, 0) = ori;
            mdtProj(1, 1) = 0.;
            mdtProj(0, 0) = 0.;
        }
        //
        Amg::VectorX locPar = parm->parameters();
        locPar[0] -= locWire[0];
        locPar[1] -= locWire[1];
        Amg::VectorX pPar = mdtProj * locPar;
        ATH_MSG_DEBUG("projected parameters(layer->MDT):" << pPar);
 
        if (parm->covariance()) {
            std::optional<AmgMatrix(5, 5)> projEM = parm->covariance()->similarity(mdtProj);
            ATH_MSG_DEBUG("projected covariance(layer->MDT):" << (*projEM));
            projPar = new Trk::AtaStraightLine(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *tubeSurf, projEM);
 
        } else {
            projPar = new Trk::AtaStraightLine(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *tubeSurf);
        }
    }
 
    if (hitType == 2) {
        //
        const MuonGM::RpcReadoutElement* rpcROE = MuonDetMgr->getRpcReadoutElement(id);
        if (!rpcROE) return projPar;
        const Trk::PlaneSurface* stripSurf = dynamic_cast<const Trk::PlaneSurface*>(&(rpcROE->surface(id)));
        if (!stripSurf) return projPar;
        // decode ref position
        Amg::VectorX locPar = parm->parameters();
        // projection matrix
        AmgMatrix(5, 5)* pMx = nullptr;
        if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
            pMx = m_rpcProjPhi.get();
        else
            pMx = m_rpcProjEta.get();
        // projected parameters
        double eta = 1.;
        double sign = (m_idHelperSvc->rpcIdHelper().measuresPhi(id) && m_idHelperSvc->rpcIdHelper().doubletPhi(id) == 2) ? -1. : 1.;
        double zswap = (lay->getRef() > 10000.) ? -1. : 1.;
        double ref = (zswap < 0.) ? lay->getRef() - 20000. : lay->getRef();
        locPar[0] -= sign * ref;
        Amg::VectorX pPar = (*pMx) * locPar;
        if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
            pPar[0] *= eta;
        else
            pPar[1] *= -eta;
        if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
            pPar[1] *= zswap;
        else
            pPar[0] *= zswap;
 
        std::optional<AmgMatrix(5, 5)> projEM;
        projEM.emplace();
        if (parm->covariance()) {
            projEM = parm->covariance()->similarity(*pMx);
            projPar = new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *stripSurf, projEM);
 
        } else {
            projPar = new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *stripSurf);
            projEM = std::nullopt;
        }
 
        if (m_alignedMode && (parm->position() - projPar->position()).mag() > 0.001) {
            ATH_MSG_DEBUG("geometrical RPC projection (layerToDetEl) for hit : "
                          << m_idHelperSvc->rpcIdHelper().measuresPhi(id) << "," << m_idHelperSvc->rpcIdHelper().stationName(id) << ","
                          << m_idHelperSvc->rpcIdHelper().stationEta(id) << "," << m_idHelperSvc->rpcIdHelper().stationPhi(id) << ","
                          << m_idHelperSvc->rpcIdHelper().doubletPhi(id) << "," << m_idHelperSvc->rpcIdHelper().doubletR(id) << ","
                          << m_idHelperSvc->rpcIdHelper().doubletZ(id));
            Amg::Vector2D locPos;
            const Amg::Vector3D& globPos = parm->position();
            bool onSurface = stripSurf->globalToLocal(globPos, globPos, locPos);
            if (onSurface) {
                pPar[0] = locPos[0];
                pPar[1] = locPos[1];
                delete projPar;
                projPar = nullptr;
                projPar = projEM ? new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *stripSurf, projEM)
                                 : new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *stripSurf);
            } else {
                delete projPar;
                return nullptr;
            }
        }
    }
 
    if (hitType == 3) {
        // local position of detEl
        const MuonGM::CscReadoutElement* cscROE = MuonDetMgr->getCscReadoutElement(id);
        if (!cscROE) {
            ATH_MSG_WARNING(name() << "CSC readout element not found");
            return projPar;
        }
 
        const Trk::PlaneSurface* stripSurf = dynamic_cast<const Trk::PlaneSurface*>(&(cscROE->surface(id)));
        if (!stripSurf) return projPar;
        // dealing with displaced planes, possibly sligthly displaced
        double diff = m_idHelperSvc->cscIdHelper().measuresPhi(id) ? 1.55 : -1.55;
        // distance between planes (assuming parallel planes)
        Amg::Vector3D layNormal = lay->surfaceRepresentation().normal();
        double DN = parm->momentum().dot(layNormal);
        double t = diff / DN;
        // displacement ( as projection of strip surface center on the layer )
        const Amg::Vector2D csc_shift(0., lay->getRef());
        // correct local position due to plane distance : use method independent on misalignment
        Amg::Vector3D corrLocPos = parm->position() - t * parm->momentum() + t * DN * layNormal;
        Amg::Vector2D locCorrLay;
        bool onSurface = lay->surfaceRepresentation().globalToLocal(corrLocPos, corrLocPos, locCorrLay);
 
        if (!onSurface) {
            ATH_MSG_WARNING(name() << ": misplaced CSC " << id);
            return projPar;
        }
 
        Amg::VectorX parProj(5);
        parProj[0] = locCorrLay[Trk::locX] - csc_shift[Trk::locX];
        parProj[1] = locCorrLay[Trk::locY] - csc_shift[Trk::locY];
        parProj[2] = parm->parameters()[Trk::phi];
        parProj[3] = parm->parameters()[Trk::theta];
        parProj[4] = parm->parameters()[Trk::qOverP];
        //
        AmgMatrix(5, 5)* pMx = nullptr;
        if (m_idHelperSvc->cscIdHelper().measuresPhi(id))
            pMx = m_tgcProjPhi.get();
        else
            pMx = m_tgcProjEta.get();
        Amg::VectorX locPar = (*pMx) * parProj;
        ATH_MSG_DEBUG("projected parameters (layer->CSC):" << m_idHelperSvc->cscIdHelper().measuresPhi(id) << "," << locPar);
 
        if (parm->covariance()) {
            std::optional<AmgMatrix(5, 5)> projEM = parm->covariance()->similarity(*pMx);
 
            ATH_MSG_DEBUG("projected covariance (layer->CSC):" << *projEM);
 
            projPar = new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *stripSurf, projEM);
 
        } else
            projPar = new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *stripSurf);
 
        ATH_MSG_DEBUG("test CSC projection:layerToDetEl:" << parm->position() << "," << projPar->position());
    }
 
    if (hitType == 4) {
        // local position at layer
        const MuonGM::TgcReadoutElement* tgcROE = MuonDetMgr->getTgcReadoutElement(id);
        if (!tgcROE) {
            ATH_MSG_WARNING(name() << "TGC readout element not found");
            return projPar;
        }
        const Trk::PlaneSurface* stripSurf = dynamic_cast<const Trk::PlaneSurface*>(&(tgcROE->surface(id)));
        if (!stripSurf) return projPar;
        //
        AmgMatrix(5, 5)* pMx = nullptr;
        if (m_idHelperSvc->tgcIdHelper().isStrip(id))
            pMx = m_tgcProjPhi.get();
        else
            pMx = m_tgcProjEta.get();
        Amg::VectorX locPar = (*pMx) * parm->parameters();
        ATH_MSG_DEBUG("projected parameters (layer->TGC):" << m_idHelperSvc->tgcIdHelper().isStrip(id) << "," << locPar << ","
                                                           << stripSurf);
 
        std::optional<AmgMatrix(5, 5)> projEM = std::nullopt;
        bool bcov = false;
 
        if (parm->covariance()) {
            projEM = parm->covariance()->similarity(*pMx);
            bcov = true;
            ATH_MSG_DEBUG("projected covariance (layer->TGC):" << (*projEM));
            projPar = new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *stripSurf, projEM);
 
        } else {
            projPar = new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *stripSurf);
        }
 
        // verify
        if (m_alignedMode && (parm->position() - projPar->position()).mag() > 0.001) {
            ATH_MSG_DEBUG("geometrical TGC projection ( layer2detEl ):"
                          << m_idHelperSvc->tgcIdHelper().stationName(id) << "," << m_idHelperSvc->tgcIdHelper().stationEta(id) << ","
                          << m_idHelperSvc->tgcIdHelper().stationPhi(id) << "," << m_idHelperSvc->tgcIdHelper().isStrip(id));
            Amg::Vector2D locPos;
            const Amg::Vector3D& globPos = parm->position();
            bool onSurface = stripSurf->globalToLocal(globPos, globPos, locPos);
            if (onSurface) {
                locPar[0] = locPos[0];
                locPar[1] = locPos[1];
                delete projPar;
                projPar = nullptr;
                projPar = bcov ? new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *stripSurf, projEM)
                               : new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *stripSurf);
            } else {
                delete projPar;
                return nullptr;
            }
        }
    }
 
    return projPar;
}

◆ match()

const Trk::Layer * Muon::MuonTGMeasurementTool::match	(	Identifier	id,
		const Trk::Layer *	lay ) const

overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1184 of file MuonTGMeasurementTool.cxx.

                                                                                         {
    const Trk::Layer* mLay = nullptr;
    if (!id.get_identifier32().get_compact() || !lay || !lay->layerType()) return mLay;
 
    Identifier layId(lay->layerType());
 
    if (m_idHelperSvc->isMdt(id) && m_idHelperSvc->isMdt(layId)) {
        if (m_idHelperSvc->mdtIdHelper().multilayer(layId) == m_idHelperSvc->mdtIdHelper().multilayer(id) &&
            m_idHelperSvc->mdtIdHelper().tubeLayer(layId) == m_idHelperSvc->mdtIdHelper().tubeLayer(id))
            return lay;
    }
 
    if (m_idHelperSvc->isRpc(id) && m_idHelperSvc->isRpc(layId)) {
        if (m_idHelperSvc->rpcIdHelper().doubletR(layId) == m_idHelperSvc->rpcIdHelper().doubletR(id) &&
            m_idHelperSvc->rpcIdHelper().doubletZ(layId) == m_idHelperSvc->rpcIdHelper().doubletZ(id) &&
            m_idHelperSvc->rpcIdHelper().gasGap(layId) == m_idHelperSvc->rpcIdHelper().gasGap(id))
            return lay;
    }
 
    if (m_idHelperSvc->isTgc(id) && m_idHelperSvc->isTgc(layId)) {
        if (m_idHelperSvc->tgcIdHelper().gasGap(layId) == m_idHelperSvc->tgcIdHelper().gasGap(id)) return lay;
    }
 
    if (m_idHelperSvc->isCsc(id) && m_idHelperSvc->isCsc(layId)) {
        if (m_idHelperSvc->cscIdHelper().chamberLayer(layId) == m_idHelperSvc->cscIdHelper().chamberLayer(id) &&
            m_idHelperSvc->cscIdHelper().wireLayer(layId) == m_idHelperSvc->cscIdHelper().wireLayer(id))
            return lay;
    }
 
    return mLay;
}

◆ measToLayer()

const Trk::RIO_OnTrack * Muon::MuonTGMeasurementTool::measToLayer	(	const Trk::Layer *	lay,
		const Trk::TrackParameters *	parm,
		const Trk::RIO_OnTrack *	rio ) const

overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 603 of file MuonTGMeasurementTool.cxx.

                                                                                                  {
    const MuonGM::MuonDetectorManager* MuonDetMgr = m_muonDetMgr;
    if (!m_useDSManager) {
        SG::ReadCondHandle<MuonGM::MuonDetectorManager> DetectorManagerHandle{m_DetectorManagerKey};
        MuonDetMgr = DetectorManagerHandle.cptr();
        if (MuonDetMgr == nullptr) {
            ATH_MSG_ERROR("Null pointer to the read MuonDetectorManager conditions object");
            // return StatusCode::FAILURE;
        }
    }
    // Get the messaging service, print where you are
    ATH_MSG_DEBUG("MuonTGMeasurementTool::measToLayer");
    const Trk::RIO_OnTrack* projRIO = nullptr;
 
    // check input
    if (!lay || !parm || !rio) return projRIO;
 
    // check compatibility of layer info and required id ? this was already done when associating !
    Identifier id = rio->identify();
 
    ATH_MSG_DEBUG("MuonTGMeasurementTool::input ok");
 
    unsigned int hitType = 0;
    if (m_idHelperSvc->isMdt(id)) hitType = 1;
    if (m_idHelperSvc->isRpc(id)) hitType = 2;
    if (m_idHelperSvc->isCsc(id)) hitType = 3;
    if (m_idHelperSvc->isTgc(id)) hitType = 4;
 
    if (hitType == 1) {
        // local position of the tube
        int tube = m_idHelperSvc->mdtIdHelper().tube(id);
        Amg::Vector2D locWire(0., lay->getRef() + (tube - 1) * 30.035);
        if (fabs(lay->getRef()) > 10e6) {
            double sign = (lay->getRef() > 0) ? 1. : -1.;
            int dec = int(lay->getRef() / 1.e5);
            double ref0 = dec / 1.e5;
            double ref1 = lay->getRef() - dec * 1e5 - 0.5 * (sign + 1) * 1e5;
            locWire[0] = ref0;
            locWire[1] = ref1;
            int tube = m_idHelperSvc->mdtIdHelper().tube(id);
            int tubeMax = m_idHelperSvc->mdtIdHelper().tubeMax(id);
            if (tube > 6 && tubeMax - tube > 5) locWire[0] = 0.;
        }
        // direction at the layer
        Amg::Vector3D dir = parm->momentum().unit();
        Amg::Vector3D locDir = (lay->surfaceRepresentation().transform().inverse()) * dir;
        Amg::Vector3D normal = lay->surfaceRepresentation().normal();
        Amg::Vector3D wireDir(1., 0., 0.);
        //
        double ND = dir.dot(normal);
        double DL = locDir.dot(wireDir);
        double A = sqrt(1. - DL * DL);
        double A_ND = A / ND;
        //
        double locLay = A_ND * rio->localParameters()[Trk::locR] + locWire[1];
        // create (fake!) rio ( rio image on TG layer )
        IdentifierHash idHash(0);
        const MuonGM::MdtReadoutElement* mdtROE = MuonDetMgr->getMdtReadoutElement(id);
        auto cov = Amg::MatrixX();
        cov = A_ND * A_ND * rio->localCovariance();
        Muon::MdtDriftCircleStatus status = Muon::MdtStatusDriftTime;
        const Muon::MdtPrepData* mdtPrd =
          new Muon::MdtPrepData(id,
                                Amg::Vector2D(locLay, 0.),
                                cov,
                                mdtROE,
                                0,
                                0,
                                status);
        Muon::MuonDriftCircleErrorStrategyInput stratbits(
          MuonDriftCircleErrorStrategy::UnknownStrategy);
        Muon::MuonDriftCircleErrorStrategy strat(stratbits);
        const Muon::MdtDriftCircleOnTrack* mdtRio =
            new Muon::MdtDriftCircleOnTrack(mdtPrd, Trk::LocalParameters(Trk::DefinedParameter(locLay, Trk::locY)),
                                            Amg::MatrixX(mdtPrd->localCovariance()), 0., Trk::DECIDED, 0.0, strat);
        return mdtRio;
    }
 
    else if (hitType == 2) {
        //
        double eta = (m_idHelperSvc->rpcIdHelper().stationEta(id) < 0) ? -1. : 1.;
        double sign = (m_idHelperSvc->rpcIdHelper().measuresPhi(id) && m_idHelperSvc->rpcIdHelper().doubletPhi(id) == 2) ? -1. : 1.;
 
        double locPos = rio->localParameters()[Trk::locX];
 
        if (m_idHelperSvc->rpcIdHelper().measuresPhi(id)) locPos *= eta;
 
        double ref = (lay->getRef() > 10000.) ? lay->getRef() - 20000. : lay->getRef();
        double zswap = (lay->getRef() > 10000.) ? -1. : 1.;
 
        if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
            locPos += sign * ref;
        else
            locPos *= zswap;
        //
        const MuonGM::RpcReadoutElement* rpcROE = MuonDetMgr->getRpcReadoutElement(id);
        const IdentifierHash idHash(0);
        std::vector<Identifier> rdoList;
        rdoList.push_back(id);
        const Muon::RpcPrepData* rpcPrd = new Muon::RpcPrepData(id, idHash, Amg::Vector2D(locPos, 0.), rdoList,
                                                                Amg::MatrixX(rio->localCovariance()), rpcROE, float(0.), 0, 0);
        const Muon::RpcClusterOnTrack* rpcRio = nullptr;
        if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
            rpcRio = new Muon::RpcClusterOnTrack(rpcPrd, Trk::LocalParameters(Trk::DefinedParameter(locPos, Trk::locX)),
                                                 Amg::MatrixX(rpcPrd->localCovariance()), parm->localPosition()[Trk::locY]);
        else
            rpcRio = new Muon::RpcClusterOnTrack(rpcPrd, Trk::LocalParameters(Trk::DefinedParameter(locPos, Trk::locY)),
                                                 Amg::MatrixX(rpcPrd->localCovariance()), parm->localPosition()[Trk::locX]);
        return rpcRio;
    }
 
    else if (hitType == 3) {
        // dealing with parallel displaced planes
        double diff = m_idHelperSvc->cscIdHelper().measuresPhi(id) ? 1.55 : -1.55;
        //
        Amg::Vector3D layNormal = lay->surfaceRepresentation().normal();
        double DN = parm->momentum().dot(layNormal);
        double t = diff / DN;
        // displacement of planes
        const Amg::Vector2D csc_shift(0., lay->getRef());
        Amg::Vector3D corrLocPos = lay->surfaceRepresentation().center() - t * parm->momentum() + t * DN * layNormal;
        Amg::Vector2D locCorrLay;
        bool onSurface = lay->surfaceRepresentation().globalToLocal(corrLocPos, corrLocPos, locCorrLay);
        //
        // if( !locCorrLay.size()==2) { ??
        if (!onSurface) { return projRIO; }
        double locPos;
        if (m_idHelperSvc->cscIdHelper().measuresPhi(id)) {
            locPos = rio->localParameters()[Trk::locX] + locCorrLay[Trk::locX] + csc_shift[Trk::locX];
        } else {
            locPos = rio->localParameters()[Trk::locX] + locCorrLay[Trk::locY] + csc_shift[Trk::locY];
        }
        //
        const MuonGM::CscReadoutElement* cscROE = MuonDetMgr->getCscReadoutElement(id);
        IdentifierHash idHash(0);
        std::vector<Identifier> rdoList;
        rdoList.push_back(id);
        Muon::CscClusterStatus status = Muon::CscStatusSimple;
        const Muon::CscPrepData* cscPrd =
          new Muon::CscPrepData(id,
                                idHash,
                                Amg::Vector2D(locPos, 0.),
                                rdoList,
                                Amg::MatrixX(rio->localCovariance()),
                                cscROE,
                                0,
                                0.,
                                status);
        const Muon::CscClusterOnTrack* cscRio = nullptr;
        if (m_idHelperSvc->cscIdHelper().measuresPhi(id))
            cscRio = new Muon::CscClusterOnTrack(cscPrd, Trk::LocalParameters(Trk::DefinedParameter(locPos, Trk::locX)),
                                                 Amg::MatrixX(cscPrd->localCovariance()), parm->localPosition()[Trk::locY], cscPrd->status());
        else
            cscRio = new Muon::CscClusterOnTrack(cscPrd, Trk::LocalParameters(Trk::DefinedParameter(locPos, Trk::locY)),
                                                 Amg::MatrixX(cscPrd->localCovariance()), parm->localPosition()[Trk::locX], cscPrd->status());
        return cscRio;
    }
 
    else if (hitType == 4) {
        //
        double locPos = rio->localParameters()[Trk::locX];
        //
        const MuonGM::TgcReadoutElement* tgcROE = MuonDetMgr->getTgcReadoutElement(id);
        IdentifierHash idHash(0);
        std::vector<Identifier> rdoList;
        rdoList.push_back(id);
        const Muon::TgcPrepData* tgcPrd =
          new Muon::TgcPrepData(id,
                                idHash,
                                Amg::Vector2D(locPos, 0.),
                                rdoList,
                                Amg::MatrixX(rio->localCovariance()),
                                tgcROE);
        const Muon::TgcClusterOnTrack* tgcRio = nullptr;
        if (m_idHelperSvc->tgcIdHelper().isStrip(id)) {
            Amg::Vector2D loc(locPos, parm->localPosition()[Trk::locY]);
            tgcRio = new Muon::TgcClusterOnTrack(tgcPrd, Trk::LocalParameters(Trk::DefinedParameter(locPos, Trk::locX)),
                                                 Amg::MatrixX(tgcPrd->localCovariance()), parm->localPosition()[Trk::locY]);
        } else {
            tgcRio = new Muon::TgcClusterOnTrack(tgcPrd, Trk::LocalParameters(Trk::DefinedParameter(locPos, Trk::locY)),
                                                 Amg::MatrixX(tgcPrd->localCovariance()), parm->localPosition()[Trk::locX]);
        }
        return tgcRio;
    }
 
    ATH_MSG_DEBUG("unknown hit technology");
    return projRIO;
}

◆ msg()

MsgStream & AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

◆ msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

◆ nearestDetEl()

const Identifier Muon::MuonTGMeasurementTool::nearestDetEl	(	const Trk::Layer *	lay,
		const Trk::TrackParameters *	parm,
		bool	measPhi,
		double &	pitch ) const

overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 793 of file MuonTGMeasurementTool.cxx.

                                                                                {
    const MuonGM::MuonDetectorManager* MuonDetMgr = m_muonDetMgr;
    if (!m_useDSManager) {
        SG::ReadCondHandle<MuonGM::MuonDetectorManager> DetectorManagerHandle{m_DetectorManagerKey};
        MuonDetMgr = DetectorManagerHandle.cptr();
        if (MuonDetMgr == nullptr) { ATH_MSG_ERROR("Null pointer to the read MuonDetectorManager conditions object"); }
    }
    // Get the messaging service, print where you are
    ATH_MSG_DEBUG("MuonTGMeasurementTool::nearestDetEl");
    Identifier nid(0);
    // check input
    if (!lay || !parm || !lay->layerType()) return nid;
 
    // check compatibility of layer info and required id ? this was already done when associating !
    Identifier layId(lay->layerType());
 
    unsigned int hitType = 0;
    if (m_idHelperSvc->isMdt(layId)) hitType = 1;
    if (m_idHelperSvc->isRpc(layId)) hitType = 2;
    if (m_idHelperSvc->isCsc(layId)) hitType = 3;
    if (m_idHelperSvc->isTgc(layId)) hitType = 4;
 
    if (hitType == 0) {
        ATH_MSG_DEBUG("unknown hit technology");
        return nid;
    }
 
    if (hitType == 1) {
        const MuonGM::MdtReadoutElement* mdtROE = MuonDetMgr->getMdtReadoutElement(layId);
        if (!mdtROE) return nid;
        int tMax = mdtROE->getNtubesperlayer();
        // local position at layer
        Amg::Vector2D locLay;
        const Amg::Vector3D& globPos = parm->position();
        bool onSurface = lay->surfaceRepresentation().globalToLocal(globPos, globPos, locLay);
        if (!onSurface || measPhi) { return nid; }
        // nearest tube index
        double refL = lay->getRef();
        if (fabs(refL) > 10e6) {
            double sign = (lay->getRef() > 0) ? 1. : -1.;
            int dec = int(lay->getRef() / 1.e5);
            refL = lay->getRef() - dec * 1e5 - 0.5 * (sign + 1) * 1e5;
        }
        double dloc = locLay[Trk::locY] - refL + 15.0175;
        int itube = int(dloc / 30.035) + 1;
        if (itube < 1 || itube > tMax || dloc < 0.) { return nid; }
        // tube id
        Identifier nearId = m_idHelperSvc->mdtIdHelper().channelID(
            m_idHelperSvc->mdtIdHelper().stationName(layId), m_idHelperSvc->mdtIdHelper().stationEta(layId),
            m_idHelperSvc->mdtIdHelper().stationPhi(layId), m_idHelperSvc->mdtIdHelper().multilayer(layId),
            m_idHelperSvc->mdtIdHelper().tubeLayer(layId), itube);
        // check if position within active volume
        if (fabs(locLay[Trk::locX]) > 0.5 * mdtROE->getActiveTubeLength(m_idHelperSvc->mdtIdHelper().tubeLayer(layId), itube)) {
            return nid;
        }
        //
        if (m_idHelperSvc->mdtIdHelper().valid(nearId))
            return nearId;
        else
            return nid;
    }
 
    if (hitType == 2) {
        // identify phi doublet first
        Identifier phiId(0);
        // both phi and eta needs to be within bounds !!!
        int doubletPhi = 1;
        bool foundDoubletPhi = false;
        const MuonGM::RpcReadoutElement* rpcROE = nullptr;
        while (!foundDoubletPhi && doubletPhi < 3) {
            Identifier refPhi1 = m_idHelperSvc->rpcIdHelper().channelID(
                m_idHelperSvc->rpcIdHelper().stationName(layId), m_idHelperSvc->rpcIdHelper().stationEta(layId),
                m_idHelperSvc->rpcIdHelper().stationPhi(layId), m_idHelperSvc->rpcIdHelper().doubletR(layId),
                m_idHelperSvc->rpcIdHelper().doubletZ(layId), doubletPhi, m_idHelperSvc->rpcIdHelper().gasGap(layId), 1, 1);
            rpcROE = MuonDetMgr->getRpcReadoutElement(refPhi1);
            if (!rpcROE) return nid;
            if (!m_idHelperSvc->rpcIdHelper().valid(refPhi1)) return nid;
            int nStripPhi = rpcROE->Nstrips(1);
 
            Identifier refPhiN = m_idHelperSvc->rpcIdHelper().channelID(
                m_idHelperSvc->rpcIdHelper().stationName(layId), m_idHelperSvc->rpcIdHelper().stationEta(layId),
                m_idHelperSvc->rpcIdHelper().stationPhi(layId), m_idHelperSvc->rpcIdHelper().doubletR(layId),
                m_idHelperSvc->rpcIdHelper().doubletZ(layId), doubletPhi, m_idHelperSvc->rpcIdHelper().gasGap(layId), 1, nStripPhi);
            //
            Amg::Vector2D loc1;
            const Amg::Vector3D globPos1 = rpcROE->stripPos(refPhi1);
            bool onSurface1 = rpcROE->surface(refPhi1).globalToLocal(globPos1, globPos1, loc1);
            Amg::Vector2D locN;
            const Amg::Vector3D globPosN = rpcROE->stripPos(refPhiN);
            bool onSurfaceN = rpcROE->surface(refPhiN).globalToLocal(globPosN, globPosN, locN);
            Amg::Vector2D loc;
            const Amg::Vector3D& glob = parm->position();
            bool onSurface = rpcROE->surface(refPhiN).globalToLocal(glob, glob, loc);
 
            int strip = 0;
            if (onSurface1 && onSurfaceN && onSurface) {
                pitch = (locN[Trk::locX] - loc1[Trk::locX]) / fmax(1, nStripPhi - 1);
                double dstrip = (loc[Trk::locX] - loc1[Trk::locX]) / pitch + 0.5;
                strip = dstrip >= 0. ? int(dstrip) + 1 : 0;
 
                if (strip > 0 && strip <= nStripPhi) {
                    // correct doublet
                    phiId = m_idHelperSvc->rpcIdHelper().channelID(
                        m_idHelperSvc->rpcIdHelper().stationName(layId), m_idHelperSvc->rpcIdHelper().stationEta(layId),
                        m_idHelperSvc->rpcIdHelper().stationPhi(layId), m_idHelperSvc->rpcIdHelper().doubletR(layId),
                        m_idHelperSvc->rpcIdHelper().doubletZ(layId), doubletPhi, m_idHelperSvc->rpcIdHelper().gasGap(layId), 1, strip);
                    foundDoubletPhi = true;
 
                } else {
                    doubletPhi += 1;
                }
            } else
                doubletPhi += 1;
        }
        // no valid phi strip - skip search
        if (!foundDoubletPhi) return nid;
 
        // eta strip
        Identifier etaId(0);
        int doubletZ = 1;
        while (doubletZ < 4) {
            Identifier refEta1 = m_idHelperSvc->rpcIdHelper().channelID(
                m_idHelperSvc->rpcIdHelper().stationName(layId), m_idHelperSvc->rpcIdHelper().stationEta(layId),
                m_idHelperSvc->rpcIdHelper().stationPhi(layId), m_idHelperSvc->rpcIdHelper().doubletR(layId), doubletZ, doubletPhi,
                m_idHelperSvc->rpcIdHelper().gasGap(layId), 0, 1);
            rpcROE = MuonDetMgr->getRpcReadoutElement(refEta1);
            if (!rpcROE) return nid;
            if (!m_idHelperSvc->rpcIdHelper().valid(refEta1)) return nid;
            int nStrips = rpcROE->Nstrips(0);
 
            Identifier refEtaN = m_idHelperSvc->rpcIdHelper().channelID(
                m_idHelperSvc->rpcIdHelper().stationName(layId), m_idHelperSvc->rpcIdHelper().stationEta(layId),
                m_idHelperSvc->rpcIdHelper().stationPhi(layId), m_idHelperSvc->rpcIdHelper().doubletR(layId), doubletZ, doubletPhi,
                m_idHelperSvc->rpcIdHelper().gasGap(layId), 0, nStrips);
            //
 
            Amg::Vector2D loc1;
            const Amg::Vector3D globPos1 = rpcROE->stripPos(refEta1);
            bool onSurface1 = rpcROE->surface(refEta1).globalToLocal(globPos1, globPos1, loc1);
            Amg::Vector2D locN;
            const Amg::Vector3D globPosN = rpcROE->stripPos(refEtaN);
            bool onSurfaceN = rpcROE->surface(refEtaN).globalToLocal(globPosN, globPosN, locN);
            Amg::Vector2D loc;
            const Amg::Vector3D& glob = parm->position();
            bool onSurface = rpcROE->surface(refEtaN).globalToLocal(glob, glob, loc);
 
            int strip = 0;
            if (onSurface && onSurface1 && onSurfaceN) {
                pitch = (locN[Trk::locX] - loc1[Trk::locX]) / fmax(1, nStrips - 1);
                // strip = int( (refPar->localPosition()[Trk::locX]-loc1[Trk::locX])/pitch+0.5 )+1;
                double dstrip = (loc[Trk::locX] - loc1[Trk::locX]) / pitch + 0.5;
                strip = dstrip >= 0. ? int(dstrip) + 1 : 0;
 
                if (strip > 0 && strip <= nStrips) {
                    etaId = m_idHelperSvc->rpcIdHelper().channelID(
                        m_idHelperSvc->rpcIdHelper().stationName(layId), m_idHelperSvc->rpcIdHelper().stationEta(layId),
                        m_idHelperSvc->rpcIdHelper().stationPhi(layId), m_idHelperSvc->rpcIdHelper().doubletR(layId), doubletZ, doubletPhi,
                        m_idHelperSvc->rpcIdHelper().gasGap(layId), 0, strip);
                    if (measPhi)
                        return phiId;
                    else
                        return etaId;
                }
            }
            doubletZ += 1;
        }
 
        return nid;
    }
 
    if (hitType == 3) {
        // ref id
        Identifier refId = m_idHelperSvc->cscIdHelper().channelID(
            m_idHelperSvc->cscIdHelper().stationName(layId), m_idHelperSvc->cscIdHelper().stationEta(layId),
            m_idHelperSvc->cscIdHelper().stationPhi(layId), m_idHelperSvc->cscIdHelper().chamberLayer(layId),
            m_idHelperSvc->cscIdHelper().wireLayer(layId), measPhi, 1);
        if (!m_idHelperSvc->cscIdHelper().valid(refId)) return nid;
        // residual in ref frame
        const Trk::TrackParameters* refPar = layerToDetEl(lay, parm, refId);
        if (!refPar) return nid;
        //
        const MuonGM::CscReadoutElement* cscROE = MuonDetMgr->getCscReadoutElement(refId);
        if (!cscROE) {
            delete refPar;
            return nid;
        }
        pitch = cscROE->StripPitch(measPhi);
        int nStrips = m_idHelperSvc->cscIdHelper().stripMax(refId);
        if (nStrips < 1) {
            delete refPar;
            return nid;
        }
 
        Identifier refIdN = m_idHelperSvc->cscIdHelper().channelID(
            m_idHelperSvc->cscIdHelper().stationName(refId), m_idHelperSvc->cscIdHelper().stationEta(refId),
            m_idHelperSvc->cscIdHelper().stationPhi(refId), m_idHelperSvc->cscIdHelper().chamberLayer(refId),
            m_idHelperSvc->cscIdHelper().wireLayer(refId), measPhi, nStrips);
        //
        Amg::Vector3D loc1 = cscROE->surface(refId).transform().inverse() * cscROE->stripPos(refId);
        Amg::Vector3D locN = cscROE->surface(refIdN).transform().inverse() * cscROE->stripPos(refIdN);
        int strip = 0;
 
        pitch = (locN[0] - loc1[0]) / fmax(1, nStrips - 1);
        strip = int((refPar->localPosition()[Trk::locX] - loc1[0]) / pitch + 0.5) + 1;
 
        delete refPar;
        refPar = nullptr;
        if (strip > 0 && strip <= nStrips) {
            // strip id
            Identifier nearId = m_idHelperSvc->cscIdHelper().channelID(
                m_idHelperSvc->cscIdHelper().stationName(layId), m_idHelperSvc->cscIdHelper().stationEta(layId),
                m_idHelperSvc->cscIdHelper().stationPhi(layId), m_idHelperSvc->cscIdHelper().chamberLayer(layId),
                m_idHelperSvc->cscIdHelper().wireLayer(layId), measPhi, strip);
            if (fabs(residual(parm, nearId)) > 0.5 * pitch)
                ATH_MSG_DEBUG("nearest CSC channel residual too large: " << residual(parm, nearId));
            return nearId;
        }
    }
 
    if (hitType == 4) {  // ref id
        if (measPhi && m_idHelperSvc->tgcIdHelper().gasGap(layId) == 2 && m_idHelperSvc->tgcIdHelper().gasGapMax(layId) == 3)
            return nid;  // no phi strips here
                         // ref id
        Identifier refId = layId;
        if (measPhi)
            refId = m_idHelperSvc->tgcIdHelper().channelID(
                m_idHelperSvc->tgcIdHelper().stationName(layId), m_idHelperSvc->tgcIdHelper().stationEta(layId),
                m_idHelperSvc->tgcIdHelper().stationPhi(layId), m_idHelperSvc->tgcIdHelper().gasGap(layId), measPhi, 1);
        if (!m_idHelperSvc->tgcIdHelper().valid(refId)) return nid;
        // residual in ref frame
        const Trk::TrackParameters* refPar = layerToDetEl(lay, parm, refId);
        if (!refPar) return nid;
        //
        const MuonGM::TgcReadoutElement* tgcROE = MuonDetMgr->getTgcReadoutElement(layId);
        if (!tgcROE) {
            delete refPar;
            return nid;
        }
        int nStrips = m_idHelperSvc->tgcIdHelper().channelMax(refId);
 
        if (nStrips < 1) {
            delete refPar;
            return nid;
        }
 
        Identifier refIdN = m_idHelperSvc->tgcIdHelper().channelID(
            m_idHelperSvc->tgcIdHelper().stationName(layId), m_idHelperSvc->tgcIdHelper().stationEta(layId),
            m_idHelperSvc->tgcIdHelper().stationPhi(layId), m_idHelperSvc->tgcIdHelper().gasGap(layId), measPhi, nStrips);
        //
        Amg::Vector2D loc1;
        const Amg::Vector3D glob1 = tgcROE->channelPos(refId);
        bool onSurface1 = tgcROE->surface(refId).globalToLocal(glob1, glob1, loc1);
        Amg::Vector2D locN;
        const Amg::Vector3D globN = tgcROE->channelPos(refIdN);
        bool onSurfaceN = tgcROE->surface(refIdN).globalToLocal(globN, globN, locN);
 
        int strip = 0;
        if (onSurface1 && onSurfaceN) {
            pitch = (locN[Trk::locX] - loc1[Trk::locX]) / fmax(1, nStrips - 1);
            strip = int((refPar->localPosition()[Trk::locX] - loc1[Trk::locX]) / pitch + 0.5) + 1;
        } else {
            ATH_MSG_DEBUG("local position of boundary elements not retrieved, return 0 ");
            delete refPar;
            return nid;
        }
 
        if (strip > 0 && strip <= nStrips) {
            // check second coordinate for active volume
            if (!measPhi && fabs(refPar->localPosition()[Trk::locY]) >
                                tgcROE->gangCentralWidth(m_idHelperSvc->tgcIdHelper().gasGap(layId), strip)) {
                delete refPar;
                return nid;
            }
            Identifier nearId = m_idHelperSvc->tgcIdHelper().channelID(
                m_idHelperSvc->tgcIdHelper().stationName(layId), m_idHelperSvc->tgcIdHelper().stationEta(layId),
                m_idHelperSvc->tgcIdHelper().stationPhi(layId), m_idHelperSvc->tgcIdHelper().gasGap(layId), measPhi, strip);
            Amg::Vector3D stripposition = tgcROE->surface(nearId).transform().inverse() * tgcROE->channelPos(nearId);
            Amg::Vector3D localhit = tgcROE->surface(nearId).transform().inverse() * parm->position();
 
            int plane = m_idHelperSvc->tgcIdHelper().gasGap(nearId);
            if (m_idHelperSvc->tgcIdHelper().isStrip(nearId))
                pitch = tgcROE->stripPitch(plane, m_idHelperSvc->tgcIdHelper().channel(nearId), localhit[1]);
            int last = 0;
            while (fabs(stripposition[0] - localhit[0]) > 0.5 * pitch) {
                if (stripposition[0] < localhit[0]) {
                    if (last != -1) {
                        strip += 1;
                        last = 1;
                    } else
                        break;
                } else {
                    if (last != 1) {
                        strip -= 1;
                        last = -1;
                    } else
                        break;
                }
                if (strip < 1 || strip > nStrips) break;
                nearId = m_idHelperSvc->tgcIdHelper().channelID(
                    m_idHelperSvc->tgcIdHelper().stationName(layId), m_idHelperSvc->tgcIdHelper().stationEta(layId),
                    m_idHelperSvc->tgcIdHelper().stationPhi(layId), m_idHelperSvc->tgcIdHelper().gasGap(layId), measPhi, strip);
                stripposition = tgcROE->surface(nearId).transform().inverse() * tgcROE->channelPos(nearId);
                localhit = tgcROE->surface(nearId).transform().inverse() * parm->position();
                if (m_idHelperSvc->tgcIdHelper().isStrip(nearId))
                    pitch = tgcROE->stripPitch(plane, m_idHelperSvc->tgcIdHelper().channel(nearId), localhit[1]);
            }
            delete refPar;
            if (strip < 1 || strip > nStrips) return nid;
            if (m_idHelperSvc->tgcIdHelper().valid(nearId))
                return nearId;
            else
                return nid;
        }
        delete refPar;
    }
    return nid;
}

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::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< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

◆ residual() [1/4]

double Muon::MuonTGMeasurementTool::residual	(	const Trk::Layer *	layer,
		const Trk::TrackParameters *	layPar,
		const Trk::RIO_OnTrack *	rio ) const

overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1238 of file MuonTGMeasurementTool.cxx.

                                                                              {
    double res = 10000.;
    if (!layer || !layPar || !rio) return res;
 
    const Trk::TrackParameters* detElPar = layerToDetEl(layer, layPar, rio->identify());
    if (!detElPar) return res;
    if (m_idHelperSvc->isMdt(rio->identify())) {
        res = fabs(detElPar->localPosition()[Trk::locR] - rio->localParameters()[Trk::locR]);
    } else {
        res = fabs(detElPar->localPosition()[Trk::locX] - rio->localParameters()[Trk::locX]);
    }
    delete detElPar;
    return res;
}

◆ residual() [2/4]

double Muon::MuonTGMeasurementTool::residual	(	const Trk::Layer *	layer,
		const Trk::TrackParameters *	layPar,
		Identifier	id ) const

overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1254 of file MuonTGMeasurementTool.cxx.

                                                                                                                       {
    const MuonGM::MuonDetectorManager* MuonDetMgr = m_muonDetMgr;
    if (!m_useDSManager) {
        SG::ReadCondHandle<MuonGM::MuonDetectorManager> DetectorManagerHandle{m_DetectorManagerKey};
        MuonDetMgr = DetectorManagerHandle.cptr();
        if (MuonDetMgr == nullptr) {
            ATH_MSG_ERROR("Null pointer to the read MuonDetectorManager conditions object");
            // return StatusCode::FAILURE;
        }
    }
 
    double res = 10000.;
    if (!layer || !layPar || !id.get_identifier32().get_compact()) return res;
 
    const Trk::TrackParameters* detElPar = layerToDetEl(layer, layPar, id);
    if (!detElPar) return res;
    if (m_idHelperSvc->isMdt(id)) {
        res = detElPar->localPosition()[Trk::locR];
    } else if (m_idHelperSvc->isRpc(id)) {
        const MuonGM::RpcReadoutElement* rpcROE = MuonDetMgr->getRpcReadoutElement(id);
        if (rpcROE)
            res = detElPar->localPosition()[Trk::locX] -
                  (detElPar->associatedSurface().transform().inverse() * (rpcROE->stripPos(id)))[Trk::locX];
    } else if (m_idHelperSvc->isCsc(id)) {
        const MuonGM::CscReadoutElement* cscROE = MuonDetMgr->getCscReadoutElement(id);
        if (cscROE)
            res = detElPar->localPosition()[Trk::locX] - (detElPar->associatedSurface().transform().inverse() * (cscROE->stripPos(id)))[0];
    } else if (m_idHelperSvc->isTgc(id)) {
        if (m_idHelperSvc->tgcIdHelper().isStrip(id) && m_idHelperSvc->tgcIdHelper().gasGap(id) == 2 &&
            m_idHelperSvc->tgcIdHelper().gasGapMax(id) == 3) {
            delete detElPar;
            return res;  // no phi strips here
        }
        const MuonGM::TgcReadoutElement* tgcROE = MuonDetMgr->getTgcReadoutElement(id);
        if (tgcROE) {
            Amg::Vector2D locPos;
            const Amg::Vector3D globPos = tgcROE->channelPos(id);
            bool onSurface = detElPar->associatedSurface().globalToLocal(globPos, globPos, locPos);
            if (onSurface) res = detElPar->localPosition()[Trk::locX] - locPos[Trk::locX];
        }
    }
    delete detElPar;
    return res;
}

◆ residual() [3/4]

double Muon::MuonTGMeasurementTool::residual	(	const Trk::TrackParameters *	layPar,
		const Trk::RIO_OnTrack *	rio ) const

overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1216 of file MuonTGMeasurementTool.cxx.

                                                                                                            {
    double res = 10000.;
    if (!layPar || !rio) return res;
 
    Amg::Vector3D gp = layPar->position();
    const Trk::Layer* layer = associatedLayer(rio->identify(), gp);
    if (!layer) return res;
 
    return residual(layer, layPar, rio);
}

◆ residual() [4/4]

double Muon::MuonTGMeasurementTool::residual	(	const Trk::TrackParameters *	layPar,
		Identifier &	id ) const

overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1227 of file MuonTGMeasurementTool.cxx.

                                                                                                 {
    double res = 10000.;
    if (!layPar || !id.get_identifier32().get_compact()) return res;
 
    Amg::Vector3D gp = layPar->position();
    const Trk::Layer* layer = associatedLayer(id, gp);
    if (!layer) return res;
 
    return residual(layer, layPar, id);
}

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::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< AlgTool > >::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_alignedMode

Gaudi::Property<bool> Muon::MuonTGMeasurementTool::m_alignedMode {this, "AlignedMode", true}

private

Definition at line 95 of file MuonTGMeasurementTool.h.

95{this, "AlignedMode", true};

◆ m_DetectorManagerKey

SG::ReadCondHandleKey<MuonGM::MuonDetectorManager> Muon::MuonTGMeasurementTool::m_DetectorManagerKey

private

Initial value:

{this, "DetectorManagerKey", "MuonDetectorManager",

"Key of input MuonDetectorManager condition data"}

Definition at line 77 of file MuonTGMeasurementTool.h.

77 {this, "DetectorManagerKey", "MuonDetectorManager",

78 "Key of input MuonDetectorManager condition data"};

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_idHelperSvc

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

private

Definition at line 75 of file MuonTGMeasurementTool.h.

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

◆ m_muonDetMgr

const MuonGM::MuonDetectorManager* Muon::MuonTGMeasurementTool::m_muonDetMgr

private

Definition at line 86 of file MuonTGMeasurementTool.h.

◆ m_rpcProjEta

std::unique_ptr<AmgMatrix(5, 5)> Muon::MuonTGMeasurementTool::m_rpcProjEta

private

Definition at line 91 of file MuonTGMeasurementTool.h.

◆ m_rpcProjPhi

std::unique_ptr<AmgMatrix(5, 5)> Muon::MuonTGMeasurementTool::m_rpcProjPhi

private

Definition at line 92 of file MuonTGMeasurementTool.h.

◆ m_tgcProjEta

std::unique_ptr<AmgMatrix(5, 5)> Muon::MuonTGMeasurementTool::m_tgcProjEta

private

Definition at line 89 of file MuonTGMeasurementTool.h.

◆ m_tgcProjPhi

std::unique_ptr<AmgMatrix(5, 5)> Muon::MuonTGMeasurementTool::m_tgcProjPhi

private

Definition at line 90 of file MuonTGMeasurementTool.h.

◆ m_trackingGeometryReadKey

SG::ReadCondHandleKey<Trk::TrackingGeometry> Muon::MuonTGMeasurementTool::m_trackingGeometryReadKey

private

Initial value:

{this, "TrackingGeometryReadKey", "",

"Key of input TrackingGeometry"}

Definition at line 83 of file MuonTGMeasurementTool.h.

83 {this, "TrackingGeometryReadKey", "",

84 "Key of input TrackingGeometry"};

◆ m_trackingGeometrySvc

ServiceHandle<Trk::ITrackingGeometrySvc> Muon::MuonTGMeasurementTool::m_trackingGeometrySvc

private

Initial value:

{this, "TrackingGeometrySvc",

"TrackingGeometrySvc/AtlasTrackingGeometrySvc"}

Definition at line 80 of file MuonTGMeasurementTool.h.

80 {this, "TrackingGeometrySvc",

81 "TrackingGeometrySvc/AtlasTrackingGeometrySvc"};

◆ m_useDSManager

Gaudi::Property<bool> Muon::MuonTGMeasurementTool::m_useDSManager {this, "UseDSManager", false}

private

Definition at line 96 of file MuonTGMeasurementTool.h.

96{this, "UseDSManager", false};

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ PairOfLayerID

◆ PairOfLayerPrd

◆ StoreGateSvc_t

◆ Track

Constructor & Destructor Documentation

◆ MuonTGMeasurementTool()

◆ ~MuonTGMeasurementTool()

Member Function Documentation

◆ associatedLayer() [1/2]

◆ associatedLayer() [2/2]

◆ declareGaudiProperty()

◆ declareProperty()

◆ detElToLayer()

◆ detStore()

◆ evtStore()

◆ extraDeps_update_handler()

◆ getGeometry()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ layerToDetEl()

◆ match()

◆ measToLayer()

◆ msg()

◆ msgLvl()

◆ nearestDetEl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ residual() [1/4]

◆ residual() [2/4]

◆ residual() [3/4]

◆ residual() [4/4]

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_alignedMode

◆ m_DetectorManagerKey

◆ m_detStore

◆ m_evtStore

◆ m_idHelperSvc

◆ m_muonDetMgr

◆ m_rpcProjEta

◆ m_rpcProjPhi

◆ m_tgcProjEta

◆ m_tgcProjPhi

◆ m_trackingGeometryReadKey

◆ m_trackingGeometrySvc

◆ m_useDSManager

◆ m_varHandleArraysDeclared

◆ m_vhka