MuonGM::MdtReadoutElement Class Reference

#include <MdtReadoutElement.h>

Inheritance diagram for MuonGM::MdtReadoutElement:

Classes
struct	GeoInfo

Public Member Functions
	MdtReadoutElement (GeoVFullPhysVol *pv, const std::string &stName, MuonDetectorManager *mgr)
	~MdtReadoutElement ()=default
unsigned int	nMDTinStation () const
	How many MDT chambers are in the station.
void	setNMdtInStation (unsigned int numMdt)
bool	barrel () const
	Returns whether the chamber is in the barrel (Assement on first later in stationName)
bool	endcap () const
	Returns whether the chamber is in the endcap.
int	getMultilayer () const
	Returns the multilayer represented by the readout element.
int	getNLayers () const
	Returns the number of tube layers inside the multilayer.
int	getNtubesperlayer () const
	Returns the number of tubes in each tube layer.
int	getNtubesinastep () const
	Returns the number of tubes in the endcap trapezoid sharing the same length.
void	setMultilayer (const int ml)
	Sets the multilayer number.
void	setNLayers (const int nl)
	Sets the number of layers.
bool	getWireFirstLocalCoordAlongZ (int tubeLayer, double &coord) const
bool	getWireFirstLocalCoordAlongR (int tubeLayer, double &coord) const
bool	containsId (const Identifier &id) const override
double	tubeLength (const int tubeLayer, const int tube) const
double	getActiveTubeLength (const int tubeLayer, const int tube) const
double	getWireLength (const int tubeLayer, const int tube) const
double	tubeLength (const Identifier &id) const
double	innerTubeRadius () const
	Returns the inner tube radius excluding the aluminium walls.
double	outerTubeRadius () const
	Returns the tube radius taking the thickness of the tubes into account.
double	tubePitch () const
	Returns the distance between 2 tubes in a tube layer.
const Amg::Transform3D &	localToGlobalTransf (const Identifier &id) const
const Amg::Transform3D &	localToGlobalTransf (const int tubeLayer, const int tube) const
Amg::Transform3D	nodeform_localToGlobalTransf (const Identifier &id) const
Amg::Transform3D	nodeform_localToGlobalTransf (const int tubeLayer, const int tube) const
Amg::Transform3D	globalToLocalTransf (const int tubeLayer, const int tube) const
Amg::Transform3D	globalToLocalTransf (const Identifier &id) const
Amg::Transform3D	nodeform_globalToLocalTransf (const Identifier &id) const
Amg::Transform3D	nodeform_globalToLocalTransf (const int tubeLayer, const int tube) const
Amg::Vector3D	localTubePos (const Identifier &id) const
Amg::Vector3D	localTubePos (const int tubelayer, const int tube) const
Amg::Vector3D	nodeform_localTubePos (const Identifier &id) const
Amg::Vector3D	nodeform_localTubePos (const int tubelayer, const int tube) const
Amg::Vector3D	tubePos (const Identifier &id) const
	Returns the global position of the given tube.
Amg::Vector3D	tubePos (const int tubelayer, const int tube) const
Amg::Vector3D	nodeform_tubePos (const Identifier &id) const
	Returns the global position of the tube excluding the B-line & As-built corrections.
Amg::Vector3D	nodeform_tubePos (const int tubelayer, const int tube) const
double	signedRODistanceFromTubeCentre (const Identifier &id) const
double	signedRODistanceFromTubeCentre (const int tubeLayer, const int tube) const
double	RODistanceFromTubeCentre (const Identifier &id) const
double	RODistanceFromTubeCentre (const int tubeLayer, const int tube) const
double	distanceFromRO (const Amg::Vector3D &GlobalHitPosition, const Identifier &id) const
double	distanceFromRO (const Amg::Vector3D &GlobalHitPosition, const int tubelayer, const int tube) const
int	isAtReadoutSide (const Amg::Vector3D &GlobalHitPosition, const Identifier &id) const
int	isAtReadoutSide (const Amg::Vector3D &GlobalHitPosition, const int tubelayer, const int tube) const
Amg::Vector3D	localROPos (const Identifier &id) const
Amg::Vector3D	localROPos (const int tubelayer, const int tube) const
Amg::Vector3D	ROPos (const int tubelayer, const int tube) const
Amg::Vector3D	ROPos (const Identifier &id) const
Amg::Vector3D	tubeFrame_localROPos (const int tubelayer, const int tube) const
Amg::Vector3D	tubeFrame_localROPos (const Identifier &id) const
void	setBLinePar (const BLinePar *bLine)
void	clearBLinePar ()
const BLinePar *	getBLinePar () const
void	clearCache () override final
void	fillCache () override final
virtual const Trk::Surface &	surface () const override final
	Return surface associated with this detector element.
virtual const Trk::StraightLineSurface &	surface (const Identifier &id) const override final
	Return surface associated with this identifier, which should come from the PrepRawData object (i.e.
virtual const Trk::StraightLineSurface &	surface (const int tubeLayer, const int tube) const
virtual const Trk::SurfaceBounds &	bounds () const override final
	Return the boundaries of the element.
virtual const Trk::CylinderBounds &	bounds (const Identifier &id) const override final
	Return the boundaries of the surface associated with this identifier.
virtual const Trk::CylinderBounds &	bounds (const int tubeLayer, const int tube) const
virtual const Amg::Transform3D &	transform (const Identifier &id) const override final
	Return local to global transform associated with this identifier.
virtual const Amg::Transform3D &	transform () const override final
	Return local to global transform.
const Amg::Transform3D &	transform (const int tubeLayer, const int tube) const
virtual const Amg::Vector3D &	center (const Identifier &) const override final
	Return the center of the surface associated with this identifier In the case of silicon it returns the same as center()
virtual const Amg::Vector3D &	center () const override final
	Return the center of the element.
const Amg::Vector3D &	center (const int tubeLayer, const int tube) const
virtual const Amg::Vector3D &	normal (const Identifier &) const override final
	Return the normal of the surface associated with this identifier In the case of silicon it returns the same as normal()
virtual const Amg::Vector3D &	normal () const override final
	Return the normal of the element.
std::vector< const Trk::Surface * >	surfaces () const
	returns all the surfaces contained in this detector element
Amg::Transform3D	tubeToMultilayerTransf (const Identifier &id) const
Amg::Transform3D	tubeToMultilayerTransf (const int tubeLayer, const int tube) const
Identifier	identify () const override final
	Returns the ATLAS Identifier of the MuonReadOutElement.
IdentifierHash	identifyHash () const override final
	Returns the IdentifierHash of the MuonStation, i.e.
IdentifierHash	detectorElementHash () const
	Returns the IdentifierHash of the detector element.
void	setIdentifier (const Identifier &id)
	Sets the Identifier, hashes & station names.
int	getStationIndex () const
int	getStationEta () const
int	getStationPhi () const
const Amg::Vector3D	globalPosition () const
const Amg::Transform3D &	absTransform () const
const Amg::Transform3D &	defTransform () const
virtual Amg::Vector3D	AmdbLRSToGlobalCoords (const Amg::Vector3D &x) const
virtual Amg::Transform3D	AmdbLRSToGlobalTransform () const
virtual Amg::Vector3D	GlobalToAmdbLRSCoords (const Amg::Vector3D &x) const
virtual Amg::Transform3D	GlobalToAmdbLRSTransform () const
std::string	getTechnologyType () const
const std::string &	getTechnologyName () const
void	setTechnologyName (const std::string &str)
std::string	getStationType () const
const std::string &	getStationName () const
void	setStationName (const std::string &)
double	getSsize () const
double	getRsize () const
double	getZsize () const
double	getLongSsize () const
double	getLongRsize () const
double	getLongZsize () const
double	getStationS () const
	Seems to be exclusively used by the MDTs --> Move it to MdtReadoutElement.
void	setStationS (double)
void	setLongSsize (double)
void	setLongRsize (double)
void	setLongZsize (double)
void	setSsize (double)
void	setRsize (double)
void	setZsize (double)
bool	largeSector () const
bool	smallSector () const
bool	sideA () const
bool	sideC () const
const MuonStation *	parentMuonStation () const
void	setParentMuonStation (const MuonStation *)
Amg::Transform3D	toParentStation () const
Amg::Vector3D	parentMuonStationPos () const
bool	hasCutouts () const
void	setHasCutouts (bool flag)
void	refreshCache ()
const Muon::IMuonIdHelperSvc *	idHelperSvc () const
Trk::DetectorElemType	detectorType () const override final
	Return the Detector element type.
std::string	detectorTypeString () const
	Returns a string of the Detector element type.
bool	msgLvl (const MSG::Level lvl) const
	Test the output level.
MsgStream &	msg () const
	The standard message stream.
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream.
void	setLevel (MSG::Level lvl)
	Change the current logging level.

Protected Member Functions
const MuonDetectorManager *	manager () const

Protected Attributes
double	m_Ssize {-9999.}
double	m_Rsize {-9999.}
double	m_Zsize {-9999.}
double	m_LongSsize {-9999.}
double	m_LongRsize {-9999.}
double	m_LongZsize {-9999.}
	size in the specified direction
std::string	m_techname {"TTT0"}
	MDT or RPC or TGC or CSC plus a two digits subtype; example RPC17.
std::string	m_statname {"XXX0"}
	examples are BMS5, CSS1, EML1
int	m_caching {-1}
	0 if we want to avoid caching geometry info for tracking interface
bool	m_hasCutouts {false}
	true is there are cutouts in the readdout-element

Private Member Functions
void	geoInitDone ()
double	getTubeLengthForCaching (const int tubeLayer, const int tube) const
double	getNominalTubeLengthWoCutouts (const int tubeLayer, const int tube) const
Amg::Vector3D	localNominalTubePosWoCutouts (const int tubelayer, const int tube) const
const Amg::Transform3D &	fromIdealToDeformed (const int tubelayer, const int tube) const
Amg::Vector3D	posOnDefChamWire (const Amg::Vector3D &locAMDBPos, const double width_narrow, const double width_wide, const double height, const double thickness, const Amg::Vector3D &fixedPoint) const
void	wireEndpointsAsBuilt (Amg::Vector3D &locAMDBWireEndP, Amg::Vector3D &locAMDBWireEndN, const int tubelayer, const int tube) const
Amg::Transform3D	globalTransform (const Amg::Vector3D &tubePos, const Amg::Transform3D &toDeform) const
std::unique_ptr< GeoInfo >	makeGeoInfo (const int tubelayer, const int tube) const
const GeoInfo &	geoInfo (const int tubeLayer, const int tube) const
Amg::Transform3D	deformedTransform (const int tubelayer, const int tube) const
void	initMessaging () const
	Initialize our message level and MessageSvc.

Static Private Member Functions
static Amg::Transform3D	tubeToMultilayerTransf (const Amg::Vector3D &tubePos, const Amg::Transform3D &toDeform)

Private Attributes
const MdtIdHelper &	m_idHelper {idHelperSvc()->mdtIdHelper()}
const int	m_stIdx_BIS {m_idHelper.stationNameIndex("BIS")}
const int	m_stIdx_BOL {m_idHelper.stationNameIndex("BOL")}
unsigned int	m_nMDTinStation {0}
int	m_multilayer {0}
int	m_nlayers {-1}
double	m_tubepitch {-9999.}
double	m_tubelayerpitch {-9999.}
int	m_ntubesperlayer {-1}
int	m_nsteps {-1}
int	m_ntubesinastep {-1}
double	m_tubelenStepSize {-9999.}
double	m_cutoutShift {-9999.}
double	m_endpluglength {-9999.}
double	m_deadlength {-9999.}
bool	m_inBarrel {false}
std::array< double, maxnsteps >	m_tubelength {make_array<double, maxnsteps>(-9999.)}
std::array< double, maxnlayers >	m_firstwire_x {make_array<double, maxnlayers>(-9999.)}
std::array< double, maxnlayers >	m_firstwire_y {make_array<double, maxnlayers>(-9999.)}
double	m_innerRadius {-9999.}
double	m_tubeWallThickness {-9999.}
CxxUtils::CachedValue< int >	m_zsignRO_tubeFrame {}
std::vector< CxxUtils::CachedUniquePtr< GeoInfo > >	m_tubeGeo {}
std::vector< CxxUtils::CachedUniquePtr< GeoInfo > >	m_backupTubeGeo {}
std::vector< CxxUtils::CachedUniquePtr< Amg::Transform3D > >	m_deformTransf {}
const BLinePar *	m_BLinePar {nullptr}
CxxUtils::CachedValue< Amg::Vector3D >	m_elemNormal {}
std::vector< CxxUtils::CachedUniquePtr< Trk::StraightLineSurface > >	m_tubeSurfaces {}
std::vector< CxxUtils::CachedUniquePtr< Trk::CylinderBounds > >	m_tubeBounds {}
std::atomic< bool >	m_haveTubeSurfaces {false}
	Flag whether any elements have been inserted into the corresponding vectors.
std::atomic< bool >	m_haveTubeGeo {false}
std::atomic< bool >	m_haveTubeBounds {false}
std::atomic< bool >	m_haveDeformTransf {false}
CxxUtils::CachedUniquePtr< Trk::Surface >	m_associatedSurface {}
CxxUtils::CachedUniquePtr< Trk::SurfaceBounds >	m_associatedBounds {}
ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc {"Muon::MuonIdHelperSvc/MuonIdHelperSvc", "MuonDetectorManager"}
Trk::DetectorElemType	m_type {Trk::DetectorElemType::SolidState}
Identifier	m_id {}
	extended data-collection identifier
IdentifierHash	m_idhash {0}
	data-collection hash identifier
IdentifierHash	m_detectorElIdhash {0}
	detector element hash identifier
double	m_stationS {0.}
int	m_stIdx {-1}
	Identifier field of the station index.
int	m_eta {-1}
	Identifier field of the station eta.
int	m_phi {-1}
	Identifier field of the station phi.
const MuonStation *	m_parentMuonStation {nullptr}
MuonDetectorManager *	m_muon_mgr {nullptr}
std::string	m_nm
	Message source name.
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels)
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer.
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level.
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging)

Friends
class	MuonGMR4::ReadoutGeomCnvAlg
class	MuonChamber
class	MuonChamberLite

Detailed Description

Definition at line 51 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

Constructor & Destructor Documentation

MdtReadoutElement()

MuonGM::MdtReadoutElement::MdtReadoutElement	(	GeoVFullPhysVol *	pv,
		const std::string &	stName,
		MuonDetectorManager *	mgr )

Definition at line 60 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                               :
        MuonReadoutElement(pv, mgr, Trk::DetectorElemType::Mdt) {
        // get the setting of the caching flag from the manager
        m_inBarrel = stName[0]== 'B';
 
        setStationName(stName);
    }

~MdtReadoutElement()

MuonGM::MdtReadoutElement::~MdtReadoutElement ( )

default

Member Function Documentation

absTransform()

const Amg::Transform3D & MuonGM::MuonReadoutElement::absTransform ( ) const

inlineinherited

Definition at line 199 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return getMaterialGeom()->getAbsoluteTransform(); }

AmdbLRSToGlobalCoords()

Amg::Vector3D MuonGM::MuonReadoutElement::AmdbLRSToGlobalCoords ( const Amg::Vector3D & x ) const

virtualinherited

Reimplemented in MuonGM::MMReadoutElement, and MuonGM::sTgcReadoutElement.

Definition at line 86 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

                                                                                    {       
        return AmdbLRSToGlobalTransform() * x;
    }

AmdbLRSToGlobalTransform()

Amg::Transform3D MuonGM::MuonReadoutElement::AmdbLRSToGlobalTransform ( ) const

virtualinherited

Reimplemented in MuonGM::MMReadoutElement, and MuonGM::sTgcReadoutElement.

Definition at line 90 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

                                                                      {
        return parentMuonStation()->getTransform() * parentMuonStation()->getNativeToAmdbLRS().inverse();
    }

barrel()

bool MuonGM::MdtReadoutElement::barrel

(

)

const

Returns whether the chamber is in the barrel (Assement on first later in stationName)

bounds() [1/3]

const Trk::SurfaceBounds & MuonGM::MdtReadoutElement::bounds ( ) const

finaloverridevirtual

Return the boundaries of the element.

Implements Trk::TrkDetElementBase.

Definition at line 910 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                            {
        if (!m_associatedBounds) {
            if (barrel()) {
                m_associatedBounds.set(
                    std::make_unique<Trk::RectangleBounds>(getSsize() / 2., getZsize() / 2.));
            } else {
                m_associatedBounds.set(std::make_unique<Trk::TrapezoidBounds>(
                    getSsize() / 2., getLongSsize() / 2., getRsize() / 2.));
            }
        }
        return *m_associatedBounds;
    }

bounds() [2/3]

virtual const Trk::CylinderBounds & MuonGM::MdtReadoutElement::bounds ( const Identifier & id ) const

finaloverridevirtual

Return the boundaries of the surface associated with this identifier.

For TRT this will be a straw, while for silicon it returns the same as bounds()

Implements Trk::TrkDetElementBase.

bounds() [3/3]

const Trk::CylinderBounds & MuonGM::MdtReadoutElement::bounds ( const int tubeLayer, const int tube ) const

virtual

Definition at line 850 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                {
        int istep = 0;
        int ntot_steps = m_nsteps;
 
        if (hasCutouts() && manager()->MinimalGeoFlag() == 0) {
            ntot_steps = m_nlayers * m_ntubesperlayer;
            istep = (tubeLayer - 1) * m_ntubesperlayer + tube - 1;
        } else {
            if (endcap()) istep = int((tube - 1) / m_ntubesinastep);
 
            if (istep < 0 || istep >= ntot_steps) {
                ATH_MSG_WARNING( "bounds for Element named "<< " with tech. " << getTechnologyType()
                    << " DEid = " << idHelperSvc()->toStringDetEl(identify()) << " called with: tubeL, tube " << tubeLayer << " " << tube
                    << "; step " << istep << " out of range 0-" << m_nsteps - 1 << " m_ntubesinastep " << m_ntubesinastep );
                ATH_MSG_WARNING( "Please run in DEBUG mode to get extra diagnostic; setting istep = 0" );
            }
        }
        if ((unsigned int)istep >= m_tubeBounds.size()) {
            THROW_EXCEPTION(__func__<<"("<<tubeLayer<<","<<tube<<") but m_tubeBounds.size()="<<m_tubeBounds.size()<<" for "<<
                             idHelperSvc()->toStringDetEl(identify()));
        }
        const CxxUtils::CachedUniquePtr<Trk::CylinderBounds>& ptr = m_tubeBounds.at(istep);
        if (!ptr) {
            double tubelength = getTubeLengthForCaching(tubeLayer, tube);
            ptr.set(std::make_unique<Trk::CylinderBounds>(innerTubeRadius(), 0.5 * tubelength - m_deadlength));
            if (!m_haveTubeBounds) m_haveTubeBounds = true;
        }
        return *ptr;
    }

center() [1/3]

const Amg::Vector3D & MuonGM::MdtReadoutElement::center ( ) const

finaloverridevirtual

Return the center of the element.

Implements Trk::TrkDetElementBase.

Definition at line 908 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

{ return surface().center(); }

center() [2/3]

virtual const Amg::Vector3D & MuonGM::MdtReadoutElement::center ( const Identifier & id ) const

finaloverridevirtual

Return the center of the surface associated with this identifier In the case of silicon it returns the same as center()

Implements Trk::TrkDetElementBase.

center() [3/3]

const Amg::Vector3D & MuonGM::MdtReadoutElement::center	(	const int	tubeLayer,
		const int	tube ) const

Definition at line 879 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

{ return geoInfo(tubeLayer, tube).m_center; }

clearBLinePar()

void MuonGM::MdtReadoutElement::clearBLinePar

(

)

Definition at line 69 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

{ m_BLinePar = nullptr; }

clearCache()

void MuonGM::MdtReadoutElement::clearCache ( )

finaloverridevirtual

Implements MuonGM::MuonReadoutElement.

Definition at line 922 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                       {
        ATH_MSG_DEBUG( "Clearing cache for ReadoutElement " << idHelperSvc()->toStringDetEl(identify()) );
        if (m_associatedSurface) {
            m_associatedSurface.release();
        }
        else  ATH_MSG_VERBOSE( "no associated surface to be deleted" );
 
        if (m_associatedBounds) {
            m_associatedBounds.release();
        } ATH_MSG_VERBOSE( "no associated bounds to be deleted" );
        m_elemNormal.reset();
        if (m_haveTubeSurfaces) {
            m_haveTubeSurfaces = false;
            for (auto& s : m_tubeSurfaces) { s.release(); }
        }
        if (m_haveTubeGeo) {
            m_haveTubeGeo = false;
            for (auto& g : m_tubeGeo) { g.release(); }
        }
        if (m_haveTubeBounds) {
            m_haveTubeBounds = false;
            for (auto& b : m_tubeBounds) { b.release(); }
        }
        // reset here the deform-related transforms        
        if (m_haveDeformTransf) {
            m_haveDeformTransf = false;
            for (auto& d : m_deformTransf) { d.release(); }
        }
    }

containsId()

bool MuonGM::MdtReadoutElement::containsId ( const Identifier & id ) const

overridevirtual

Implements MuonGM::MuonReadoutElement.

Definition at line 1023 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                 {
        if (idHelperSvc()->detElId(id) != identify()) return false;
        int layer = m_idHelper.tubeLayer(id);
        if (layer < 1 || layer > getNLayers()) return false;
        int tube = m_idHelper.tube(id);
        return tube >= 1 && tube <= getNtubesperlayer();
    }

deformedTransform()

Amg::Transform3D MuonGM::MdtReadoutElement::deformedTransform ( const int tubelayer, const int tube ) const

private

Definition at line 461 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                      {
 
        const MuonStation* ms = parentMuonStation();
        if (!ms->hasBLines() && !ms->hasMdtAsBuiltParams()) {
            return Amg::Transform3D::Identity();
        }
        const Amg::Vector3D fixedPoint = ms->getBlineFixedPointInAmdbLRS();
 
        // Chamber parameters
        double moduleWidthS = m_Ssize;
        double moduleWidthL = m_LongSsize;
        double height = barrel() ? ms->ZsizeMdtStation() : ms->RsizeMdtStation();
        double thickness = barrel() ? ms->RsizeMdtStation() : ms->ZsizeMdtStation();
       
        ATH_MSG_VERBOSE("Calculate deformed transform "<<idHelperSvc()->toStringDetEl(identify())
                      <<", layer: "<<tubeLayer<<", tube: "<<tube<<", fixedPoint: "<<Amg::toString(fixedPoint)<<" / "
                      <<Amg::toString(ms->getGeoTransform()->getDefTransform() *ms->getNativeToAmdbLRS().inverse()* fixedPoint)
                      <<", height: "<<height<<", thickness: "<<thickness
                      <<", ideal tube: "<<Amg::toString(localNominalTubePosWoCutouts(tubeLayer,tube)));
 
#ifndef NDEBUG
        double heightML = barrel() ? m_Zsize : m_Rsize;
        double thicknessML = barrel() ? m_Rsize : m_Zsize;
        if (std::abs(height - heightML) > 10.) {
            ATH_MSG_DEBUG( "RE " << idHelperSvc()->toStringDetEl(identify())
                << "Different ML and MDTStation length in the precision coord.direction  ---  MultiLayerHeight, MDTstationHeigh "
                << heightML << " " << height << " MultiLayerThickness, MDTstationThickness " << thicknessML << " " << thickness);
        }
#endif
        // Chamber dimension in Z is defined with extraneous glue width. Correct for it
        double glue = (tubePitch() - 2. * outerTubeRadius());
        height -= glue;
 
        // Calculate transformation from native to AMDB.
        const Amg::Transform3D toAMDB = ms->getNativeToAmdbLRS() * toParentStation();
        const Amg::Transform3D fromAMDB = toAMDB.inverse();
 
        // Get positions of the wire center and end without deformations
        Amg::Vector3D pt_center = localNominalTubePosWoCutouts(tubeLayer, tube);
        double halftubelen = 0.5 * getNominalTubeLengthWoCutouts(tubeLayer, tube);
 
        // Compute tube ends in AMDB coordinates
        Amg::Vector3D pt_end1 = toAMDB * pt_center + halftubelen * Amg::Vector3D::UnitX();  // s>0 side
        Amg::Vector3D pt_end2 = toAMDB * pt_center - halftubelen * Amg::Vector3D::UnitX();  // s>0 side
 
        Amg::Vector3D pt_end1_new = pt_end1;
        Amg::Vector3D pt_end2_new = pt_end2;
 
        // if there are as built parameters ... apply them here
        if (ms->hasMdtAsBuiltParams()) {
            wireEndpointsAsBuilt(pt_end1_new, pt_end2_new, tubeLayer, tube);
        }
 
        // if there are deformation parameters ... apply them here
        if (ms->hasBLines()) {
 
            // Get positions after the deformations applied
            // first wire end point
            pt_end1_new = posOnDefChamWire(pt_end1_new, moduleWidthS, moduleWidthL, height, thickness, fixedPoint);
 
            // second wire end point
            pt_end2_new = posOnDefChamWire(pt_end2_new, moduleWidthS, moduleWidthL, height, thickness, fixedPoint);
        }
 
        // Switch tube ends back to MGM coordinates
        pt_end1 = fromAMDB * pt_end1;
        pt_end2 = fromAMDB * pt_end2;
        pt_end1_new = fromAMDB * pt_end1_new;
        pt_end2_new = fromAMDB * pt_end2_new;
 
        // Calculate deformation. Make sure that the wire length stays the same.
        // Code in positionOnDeformedChamber does not provide this by default.
        // Break transformation into translation of the wire center and the rotation of the wire
        // Move to the coordinate system originated at the wire center, then rotate the wire, then
        // move wire center to the new position
        const Amg::Vector3D pt_center_new = 0.5 * (pt_end1_new + pt_end2_new);
        const Amg::Transform3D to_center{Amg::getTranslate3D(-pt_center)};
        const Amg::Transform3D from_center{Amg::getTranslate3D(pt_center_new)};
        const Amg::Vector3D old_direction = (pt_end2 - pt_end1).unit();
        const Amg::Vector3D new_direction = (pt_end2_new - pt_end1_new).unit();
        const Amg::Vector3D rotation_vector = old_direction.cross(new_direction);
        
 
        Amg::Transform3D deformedTransform{Amg::Transform3D::Identity()};
        if (rotation_vector.mag() > 10. * std::numeric_limits<double>::epsilon()) {
            const Amg::AngleAxis3D wire_rotation(std::asin(rotation_vector.mag()), rotation_vector.unit());
            deformedTransform = from_center * wire_rotation * to_center;
        } else {
            deformedTransform = from_center * to_center;
        }
        ATH_MSG_VERBOSE("To center "<<GeoTrf::toString(to_center)<<" from: "<<GeoTrf::toString(from_center)<<
                  " -- direction: "<<GeoTrf::toString(old_direction)<<" vs. "<<GeoTrf::toString(new_direction)
                  <<" --> rot: "<<GeoTrf::toString(rotation_vector)<<" ==> "<<GeoTrf::toString(deformedTransform,true));
        return deformedTransform;
    }

defTransform()

const Amg::Transform3D & MuonGM::MuonReadoutElement::defTransform ( ) const

inlineinherited

Definition at line 201 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return getMaterialGeom()->getDefAbsoluteTransform(); }

detectorElementHash()

IdentifierHash MuonGM::MuonReadoutElement::detectorElementHash ( ) const

inlineinherited

Returns the IdentifierHash of the detector element.

Unlike the identifyHash, this hash is unique for each readout element

Definition at line 175 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_detectorElIdhash; }

detectorType()

Trk::DetectorElemType MuonGM::MuonReadoutElement::detectorType ( ) const

inlinefinaloverridevirtualinherited

Return the Detector element type.

Implements Trk::TrkDetElementBase.

Definition at line 132 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_type; }

detectorTypeString()

std::string Trk::TrkDetElementBase::detectorTypeString ( ) const

inherited

Returns a string of the Detector element type.

Definition at line 10 of file TrkDetElementBase.cxx.

                                                         {
  auto type = detectorType();
  if (type == Trk::DetectorElemType::SolidState)
    return "SolidState";
  else if (type == Trk::DetectorElemType::Silicon)
    return "Silicon";
  else if (type == Trk::DetectorElemType::TRT)
    return "TRT";
  else if (type == Trk::DetectorElemType::Csc)
    return "Csc";
  else if (type == Trk::DetectorElemType::Mdt)
    return "Mdt";
  else if (type == Trk::DetectorElemType::Rpc)
    return "Rpc";
  else if (type == Trk::DetectorElemType::Tgc)
    return "Tgc";
  else if (type == Trk::DetectorElemType::sTgc)
    return "sTgc";
  else if (type == Trk::DetectorElemType::MM)
    return "Mm";
  return "Unknown";
}

distanceFromRO() [1/2]

double MuonGM::MdtReadoutElement::distanceFromRO	(	const Amg::Vector3D &	GlobalHitPosition,
		const Identifier &	id ) const

distanceFromRO() [2/2]

double MuonGM::MdtReadoutElement::distanceFromRO	(	const Amg::Vector3D &	GlobalHitPosition,
		const int	tubelayer,
		const int	tube ) const

Definition at line 189 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                {
        // x is given in the global reference frame
        const Amg::Vector3D cPos = center(tubeLayer, tube);
        const Amg::Vector3D roPos = ROPos(tubeLayer, tube);
        const Amg::Vector3D c_ro = cPos - roPos;
        const Amg::Vector3D x_ro = x - roPos;
 
        double scalprod = c_ro.dot(x_ro);
        double wlen = getWireLength(tubeLayer, tube);
        if (wlen > 10. * CLHEP::mm)
            scalprod = std::abs(2. * scalprod / getWireLength(tubeLayer, tube));
        else {
            ATH_MSG_WARNING( " Distance of Point " <<Amg::toString(x) << " from RO side cannot be calculated (=0) since wirelength = " << wlen );
            scalprod = 0.;
        }
        return scalprod;
    }

endcap()

bool MuonGM::MdtReadoutElement::endcap

(

)

const

Returns whether the chamber is in the endcap.

fillCache()

void MuonGM::MdtReadoutElement::fillCache ( )

finaloverridevirtual

Implements MuonGM::MuonReadoutElement.

Definition at line 961 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                      {
        ATH_MSG_DEBUG( "Filling cache for ReadoutElement " << idHelperSvc()->toStringDetEl(identify()));
 
        const Trk::PlaneSurface* tmpSurface = dynamic_cast<const Trk::PlaneSurface*>(&surface());  //<! filling m_associatedSurface
        const Trk::SurfaceBounds* tmpBounds = nullptr;                         //<! filling m_associatedBounds
        if (barrel())
            tmpBounds = dynamic_cast<const Trk::RectangleBounds*>(&bounds());
        else
            tmpBounds = dynamic_cast<const Trk::TrapezoidBounds*>(&bounds());
        ATH_MSG_VERBOSE( "global Surface / Bounds pointers " << tmpSurface << " " << tmpBounds );
        ATH_MSG_VERBOSE( "global Normal " << normal() );
 
        const Trk::CylinderBounds* tmpCil = nullptr;
        const Trk::StraightLineSurface* tmpSaggL = nullptr;
        Amg::Vector3D myPoint{Amg::Vector3D::Zero()};
        Amg::Transform3D myTransform{Amg::Transform3D::Identity()};
        for (int tl = 1; tl <= getNLayers(); ++tl) {
            for (int tube = 1; tube <= getNtubesperlayer(); ++tube) {
                // in case of BMG chambers, do not check the 'dead' tubes
                // (the tubes are numbered from 1-54 for each however there are cutouts for the
                // alignment system where no tubes are built-in, meaning, those tubes do not exist/are 'dead')
                if (manager()->mdtIdHelper()->isBMG(identify())) {
                    PVConstLink cv = getMaterialGeom();
                    // usually the tube number corresponds to the child number, however for
                    // BMG chambers full tubes are skipped during the building process
                    // therefore the matching needs to be done via the volume ID
                    int packed_id = tube + maxNTubesPerLayer * tl;
                    bool found = false;
                    geoGetIds(
                        [&](int id) {
                            if (!found && id == packed_id) {
                                myTransform = transform(tl, tube);                                           //<! filling m_tubeTransf
                                myPoint = center(tl, tube);                                                  //<! filling m_tubeCenter
                                tmpCil = dynamic_cast<const Trk::CylinderBounds*>(&bounds(tl, tube));        //<! filling m_tubeBounds
                                tmpSaggL = dynamic_cast<const Trk::StraightLineSurface*>(&surface(tl, tube));  //<! filling m_tubeSurfaces
                                found = true;
                            }
                        },
                        &*cv);
                    if (found) {
                        ATH_MSG_VERBOSE( "tubeLayer/tube " << tl << " " << tube << " transform at origin  "
                            << Amg::toString(myTransform.linear()) );
                        ATH_MSG_VERBOSE( "tubeLayer/tube " << tl << " " << tube << " tube center          " << myPoint );
                        ATH_MSG_VERBOSE( "tubeLayer/tube " << tl << " " << tube << " tube bounds pointer  " << tmpCil );
                        ATH_MSG_VERBOSE( "tubeLayer/tube " << tl << " " << tube << " tube surface pointer " << tmpSaggL );
                    }
                } else {
                    // print in order to compute !!!
                    myTransform = transform(tl, tube);                                           //<! filling m_tubeTransf
                    myPoint = center(tl, tube);                                                  //<! filling m_tubeCenter
                    tmpCil = dynamic_cast<const Trk::CylinderBounds*>(&bounds(tl, tube));        //<! filling m_tubeBounds
                    tmpSaggL = dynamic_cast<const Trk::StraightLineSurface*>(&surface(tl, tube));  //<! filling m_tubeSurfaces
                    ATH_MSG_VERBOSE( "tubeLayer/tube " << tl << " " << tube << " transform at origin  "
                                    << Amg::toString(myTransform.translation()) );
                    ATH_MSG_VERBOSE( "tubeLayer/tube " << tl << " " << tube << " tube center          " << myPoint );
                    ATH_MSG_VERBOSE( "tubeLayer/tube " << tl << " " << tube << " tube bounds pointer  " << tmpCil );
                    ATH_MSG_VERBOSE( "tubeLayer/tube " << tl << " " << tube << " tube surface pointer " << tmpSaggL );
                }
            }
        }
    }

fromIdealToDeformed()

const Amg::Transform3D & MuonGM::MdtReadoutElement::fromIdealToDeformed ( const int tubelayer, const int tube ) const

private

Definition at line 433 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                          {
        size_t itube = (tubeLayer - 1) * m_ntubesperlayer + tube - 1;
        if (itube >= m_deformTransf.size()) {
            ATH_MSG_WARNING(__func__<<"() :"<<__LINE__<< " called with tubeLayer or tube out of range in chamber "
                << idHelperSvc()->toStringDetEl(identify()) << " : layer " << tubeLayer << " max " << m_nlayers << " tube " << tube
                << " max " << m_ntubesperlayer << " will compute deformation for first tube in this chamber" );
            ATH_MSG_WARNING( "Please run in DEBUG mode to get extra diagnostic" );
            itube = 0;
        }
 
        const CxxUtils::CachedUniquePtr<Amg::Transform3D>& ptr = m_deformTransf.at(itube);
        if (!ptr) {
            Amg::Transform3D trans = deformedTransform(tubeLayer, tube);
            ptr.set(std::make_unique<Amg::Transform3D>(std::move(trans)));
            if (!m_haveDeformTransf) m_haveDeformTransf = true;
        }
        return *ptr;
    }

geoInfo()

const MdtReadoutElement::GeoInfo & MuonGM::MdtReadoutElement::geoInfo ( const int tubeLayer, const int tube ) const

private

Definition at line 807 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                        {
        size_t itube = (tubeLayer - 1) * m_ntubesperlayer + tube - 1;
        if (itube >= m_tubeGeo.size()) {
           ATH_MSG_WARNING(__func__<<"() :"<<__LINE__<<"  called with tubeLayer or tube out of range in chamber "
                << idHelperSvc()->toStringDetEl(identify()) << " : layer " << tubeLayer << " max " << m_nlayers << " tube " << tube
                << " max " << m_ntubesperlayer << " will compute transform for first tube in this chamber" );
            ATH_MSG_WARNING( "Please run in DEBUG mode to get extra diagnostic" );
            itube = 0;
        }
 
        const CxxUtils::CachedUniquePtr<GeoInfo>& ptr = m_tubeGeo.at(itube);
        if (!ptr) {
            ptr.set(makeGeoInfo(tubeLayer, tube));
            if (!m_haveTubeGeo) m_haveTubeGeo = true;
        }
        return *ptr;
    }

geoInitDone()

void MuonGM::MdtReadoutElement::geoInitDone ( )

private

Definition at line 93 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                        {
        m_tubeGeo.resize(m_nlayers * m_ntubesperlayer);
        m_deformTransf.resize(m_nlayers * m_ntubesperlayer);
        m_tubeSurfaces.resize(m_nlayers * m_ntubesperlayer);
 
        int ntot_steps = m_nsteps;
        if (hasCutouts() && manager()->MinimalGeoFlag() == 0) { ntot_steps = m_nlayers * m_ntubesperlayer; }
        m_tubeBounds.resize(ntot_steps);
    }

getActiveTubeLength()

double MuonGM::MdtReadoutElement::getActiveTubeLength	(	const int	tubeLayer,
		const int	tube ) const

getBLinePar()

const BLinePar * MuonGM::MdtReadoutElement::getBLinePar ( ) const

inline

Definition at line 156 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{ return m_BLinePar; }

getLongRsize()

double MuonGM::MuonReadoutElement::getLongRsize ( ) const

inlineinherited

Definition at line 189 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_LongRsize; }

getLongSsize()

double MuonGM::MuonReadoutElement::getLongSsize ( ) const

inlineinherited

Definition at line 188 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_LongSsize; }

getLongZsize()

double MuonGM::MuonReadoutElement::getLongZsize ( ) const

inlineinherited

Definition at line 190 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_LongZsize; }

getMultilayer()

int MuonGM::MdtReadoutElement::getMultilayer

(

)

const

Returns the multilayer represented by the readout element.

getNLayers()

int MuonGM::MdtReadoutElement::getNLayers

(

)

const

Returns the number of tube layers inside the multilayer.

getNominalTubeLengthWoCutouts()

double MuonGM::MdtReadoutElement::getNominalTubeLengthWoCutouts ( const int tubeLayer, const int tube ) const

private

Definition at line 407 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                     {
        if (barrel())
            return m_Ssize;
        else {
            int istep = int((tube - 1) / m_ntubesinastep);
            if (istep < 0 || istep >= m_nsteps) {
                ATH_MSG_WARNING( "getNominalTubeLengthWoCutouts for Element named " << idHelperSvc()->toStringDetEl(identify())
                         << " called with: tubeL, tube " << tubeLayer
                    << " " << tube << "; step " << istep << " out of range 0-" << m_nsteps - 1 << " m_ntubesinastep " << m_ntubesinastep);
                istep = 0;
            }
            return m_tubelength[istep];
        }
    }

getNtubesinastep()

int MuonGM::MdtReadoutElement::getNtubesinastep

(

)

const

Returns the number of tubes in the endcap trapezoid sharing the same length.

getNtubesperlayer()

int MuonGM::MdtReadoutElement::getNtubesperlayer

(

)

const

Returns the number of tubes in each tube layer.

getRsize()

double MuonGM::MuonReadoutElement::getRsize ( ) const

inlineinherited

Definition at line 186 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_Rsize; }

getSsize()

double MuonGM::MuonReadoutElement::getSsize ( ) const

inlineinherited

Definition at line 185 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_Ssize; }

getStationEta()

int MuonGM::MuonReadoutElement::getStationEta ( ) const

inlineinherited

Definition at line 182 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_eta; }

getStationIndex()

int MuonGM::MuonReadoutElement::getStationIndex ( ) const

inlineinherited

Definition at line 181 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_stIdx; }

getStationName()

const std::string & MuonGM::MuonReadoutElement::getStationName ( ) const

inlineinherited

Definition at line 179 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_statname; }

getStationPhi()

int MuonGM::MuonReadoutElement::getStationPhi ( ) const

inlineinherited

Definition at line 183 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_phi; }

getStationS()

double MuonGM::MuonReadoutElement::getStationS ( ) const

inlineinherited

Seems to be exclusively used by the MDTs --> Move it to MdtReadoutElement.

Definition at line 191 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_stationS; }

getStationType()

std::string MuonGM::MuonReadoutElement::getStationType ( ) const

inlineinherited

Definition at line 178 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_statname.substr(0,3); }

getTechnologyName()

const std::string & MuonGM::MuonReadoutElement::getTechnologyName ( ) const

inlineinherited

Definition at line 177 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_techname; }

getTechnologyType()

std::string MuonGM::MuonReadoutElement::getTechnologyType ( ) const

inlineinherited

Definition at line 176 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_techname.substr(0, 3); }

getTubeLengthForCaching()

double MuonGM::MdtReadoutElement::getTubeLengthForCaching ( const int tubeLayer, const int tube ) const

private

Definition at line 103 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                               {
        double nominalTubeLength = 0.;
        if (barrel())
            nominalTubeLength = m_Ssize;
        else {
            int istep = int((tube - 1) / m_ntubesinastep);
            if (istep < 0 || istep >= m_nsteps) {
                ATH_MSG_WARNING( "getTubeLenght for Element with tech. " << getTechnologyType()
                    << " DEid = " << idHelperSvc()->toStringDetEl(identify()) << " called with: tubeL, tube " << tubeLayer << " " << tube
                    << "; step " << istep << " out of range 0-" << m_nsteps - 1 << " m_ntubesinastep " << m_ntubesinastep );
                istep = 0;
            }
            nominalTubeLength = m_tubelength[istep];
        }
 
        double tlength = nominalTubeLength;
 
        if (hasCutouts()) {
            if (manager()->MinimalGeoFlag() == 0) {
                ATH_MSG_VERBOSE( " MdtReadoutElement " <<idHelperSvc()->toStringDetEl(identify())
                                << " has cutouts, check for real tube length for tubeLayer, tube " << tubeLayer << " " << tube );
                PVConstLink cv = getMaterialGeom();  // it is "Multilayer"
                int nGrandchildren = cv->getNChildVols();
                if (nGrandchildren <= 0) return tlength;
                // child vol 0 is foam; 1 to (nGrandchildren-1) should be tubes
                int ii = (tubeLayer - 1) * m_ntubesperlayer + tube;
                // BIS78 only (the BIS7 of Run1/2 has no cutouts, thus, this block won't be reached)
                if ((getStationIndex() == m_stIdx_BIS && std::abs(getStationEta()) == 7)) --ii;
                if (m_idHelper.isBMG(identify())) {
                    // usually the tube number corresponds to the child number, however for
                    // BMG chambers full tubes are skipped during the building process
                    // therefore the matching needs to be done via the volume ID
                    int packed_id = tube + maxNTubesPerLayer * tubeLayer;
                    int kk = 0;
                    bool found = false;
                    geoGetIds(
                        [&](int id) {
                            if (!found && id == packed_id) {
                                ii = kk;
                                found = true;
                            }
                            ++kk;
                        }, cv);
                    if (found) {
                        ATH_MSG_DEBUG( " MdtReadoutElement tube match found for BMG - input : tube(" << tube << "), layer("
                            << tubeLayer << ") - output match : tube(" << ii % maxNTubesPerLayer << "), layer(" << ii / maxNTubesPerLayer
                            << ")" );
                    }
                }
                if (ii >= nGrandchildren) {
                    ATH_MSG_WARNING( " MdtReadoutElement " << idHelperSvc()->toStringDetEl(identify()) << " has cutouts, nChild = " << nGrandchildren
                        << " --- getTubeLength is looking for child with index ii=" << ii << " for tubeL and tube = " << tubeLayer << " "
                        << tube );
                    ATH_MSG_WARNING( "returning nominalTubeLength " );
                    return tlength;
                }
                if (ii < 0) {
                    ATH_MSG_WARNING( " MdtReadoutElement " << idHelperSvc()->toStringDetEl(identify()) << " has cutouts, nChild = " << nGrandchildren
                        << " --- getTubeLength is looking for child with index ii=" << ii << " for tubeL and tube = " << tubeLayer << " "
                        << tube );
                    ATH_MSG_WARNING( "returning nominalTubeLength " );
                    return tlength;
                }
                PVConstLink physChild = cv->getChildVol(ii);
                const GeoShape* shape = physChild->getLogVol()->getShape();
                if (shape == nullptr) return tlength;
                const GeoTube* theTube = dynamic_cast<const GeoTube*>(shape);
                if (theTube != nullptr)
                    tlength = 2. * theTube->getZHalfLength();
                else
                    ATH_MSG_WARNING( "PhysChild with index " << ii
                        << " out of (tubeLayer-1)*m_ntubesperlayer+tube with tl=" << tubeLayer << " tubes/lay=" << m_ntubesperlayer
                        << " t=" << tube << " for  MdtReadoutElement " << idHelperSvc()->toStringDetEl(identify()) );
            }
            if (std::abs(tlength - nominalTubeLength) > 0.1) {
                ATH_MSG_VERBOSE( "Tube " << tube << " in tubeLayer = " << tubeLayer
                    << " is affected by a cutout: eff. length =  " << tlength << " while nominal = " << nominalTubeLength
                    << " in station " << idHelperSvc()->toStringDetEl(identify()));
            } else {
                ATH_MSG_VERBOSE( "Tube " << tube << " in tubeLayer = " << tubeLayer
                        << " is NOT affected by the cutout: eff. length =  " << tlength << " while nominal = " << nominalTubeLength);
            }
        }
        return tlength;
    }

getWireFirstLocalCoordAlongR()

bool MuonGM::MdtReadoutElement::getWireFirstLocalCoordAlongR	(	int	tubeLayer,
		double &	coord ) const

Definition at line 86 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                           {
        coord = -9999.;
        if (tubeLayer > getNLayers() || tubeLayer < 1) return false;
        coord = m_firstwire_y[tubeLayer - 1];
        return true;
    }

getWireFirstLocalCoordAlongZ()

bool MuonGM::MdtReadoutElement::getWireFirstLocalCoordAlongZ	(	int	tubeLayer,
		double &	coord ) const

Definition at line 80 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                           {
        coord = -9999.;
        if (tubeLayer > getNLayers() || tubeLayer < 1) return false;
        coord = m_firstwire_x[tubeLayer - 1];
        return true;
    }

getWireLength()

double MuonGM::MdtReadoutElement::getWireLength	(	const int	tubeLayer,
		const int	tube ) const

getZsize()

double MuonGM::MuonReadoutElement::getZsize ( ) const

inlineinherited

Definition at line 187 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_Zsize; }

globalPosition()

const Amg::Vector3D MuonGM::MuonReadoutElement::globalPosition ( ) const

inherited

Definition at line 47 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

{ return absTransform().translation(); }

GlobalToAmdbLRSCoords()

Amg::Vector3D MuonGM::MuonReadoutElement::GlobalToAmdbLRSCoords ( const Amg::Vector3D & x ) const

virtualinherited

Reimplemented in MuonGM::MMReadoutElement, and MuonGM::sTgcReadoutElement.

Definition at line 94 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

                                                                                    {       
        return GlobalToAmdbLRSTransform() * x;
    }

GlobalToAmdbLRSTransform()

Amg::Transform3D MuonGM::MuonReadoutElement::GlobalToAmdbLRSTransform ( ) const

virtualinherited

Reimplemented in MuonGM::MMReadoutElement, and MuonGM::sTgcReadoutElement.

Definition at line 98 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

                                                                      {
        return parentMuonStation()->getNativeToAmdbLRS() * parentMuonStation()->getTransform().inverse();
    }

globalToLocalTransf() [1/2]

Amg::Transform3D MuonGM::MdtReadoutElement::globalToLocalTransf

(

const Identifier &

id

)

const

globalToLocalTransf() [2/2]

Amg::Transform3D MuonGM::MdtReadoutElement::globalToLocalTransf	(	const int	tubeLayer,
		const int	tube ) const

Definition at line 400 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                   {
        return transform(tubeLayer, tube).inverse();
    }

globalTransform()

Amg::Transform3D MuonGM::MdtReadoutElement::globalTransform ( const Amg::Vector3D & tubePos, const Amg::Transform3D & toDeform ) const

private

Definition at line 391 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                                      {
        return transform() * tubeToMultilayerTransf(tubePos, toDeform);
    }

hasCutouts()

bool MuonGM::MuonReadoutElement::hasCutouts ( ) const

inlineinherited

Definition at line 120 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_hasCutouts; }

identify()

Identifier MuonGM::MuonReadoutElement::identify ( ) const

inlinefinaloverridevirtualinherited

Returns the ATLAS Identifier of the MuonReadOutElement.

Usually the Identifier corresponds to the first channel covered by the Element

Implements Trk::TrkDetElementBase.

Definition at line 173 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_id; }

identifyHash()

IdentifierHash MuonGM::MuonReadoutElement::identifyHash ( ) const

inlinefinaloverridevirtualinherited

Returns the IdentifierHash of the MuonStation, i.e.

2 readoutelements with same sationIndex, stationEta, stationPhi have the same identifiyHash

Implements Trk::TrkDetElementBase.

Definition at line 174 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_idhash; }

idHelperSvc()

const Muon::IMuonIdHelperSvc * MuonGM::MuonReadoutElement::idHelperSvc ( ) const

inlineinherited

Definition at line 130 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_idHelperSvc.get(); }

initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

{
  m_imsg = Athena::getMessageSvc();
  // If user did not set an explicit level, set a default
  if (m_lvl == MSG::NIL) {
    m_lvl = m_imsg ?
      static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
      MSG::INFO;
  }
}

innerTubeRadius()

double MuonGM::MdtReadoutElement::innerTubeRadius

(

)

const

Returns the inner tube radius excluding the aluminium walls.

isAtReadoutSide() [1/2]

int MuonGM::MdtReadoutElement::isAtReadoutSide	(	const Amg::Vector3D &	GlobalHitPosition,
		const Identifier &	id ) const

isAtReadoutSide() [2/2]

int MuonGM::MdtReadoutElement::isAtReadoutSide	(	const Amg::Vector3D &	GlobalHitPosition,
		const int	tubelayer,
		const int	tube ) const

Definition at line 207 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                              {
        double distance = distanceFromRO(GlobalHitPosition, tubeLayer, tube);
        if (distance < 0) {
            ATH_MSG_WARNING( "isAtReadoutSide() - GlobalHitPosition appears to be outside the tube volume " << distance );
            return 1;
        } else if (distance <= getWireLength(tubeLayer, tube) / 2.)
            return 1;
        else if (distance < getWireLength(tubeLayer, tube))
            return -1;
        else {
            ATH_MSG_WARNING( "isAtReadoutSide() - GlobalHitPosition appears to be outside the tube volume " << distance );
            return -1;
        }
    }

largeSector()

bool MuonGM::MuonReadoutElement::largeSector ( ) const

inherited

Definition at line 49 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

                                               {
        // this doesn't apply to TGC
        if(m_statname.size() >= 3){
            char c = m_statname[2];
            if (c == 'L')
                return true;
            else if (c == 'S')
                return false;
            else {
                if (c == 'E' || c == 'F' || c == 'G') return false;
                if (c == 'M' || c == 'R') return true;
            }
        }
        ATH_MSG_FATAL(__FILE__<<":"<<__LINE__<<" largeSector() - is this station a larger sector answer is no for readout element "<<m_idHelperSvc->toStringDetEl(identify()));
        throw std::runtime_error("Unknown sector");
        return false;
    }

localNominalTubePosWoCutouts()

Amg::Vector3D MuonGM::MdtReadoutElement::localNominalTubePosWoCutouts ( const int tubelayer, const int tube ) const

private

Definition at line 421 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                         {
        double xtube{0.}, ztube{0.};
        if (barrel()) {
            xtube = -m_Rsize / 2. + m_firstwire_y[tubeLayer - 1];
            ztube = -m_Zsize / 2. + m_firstwire_x[tubeLayer - 1] + (tube - 1) * m_tubepitch;
        } else {
            xtube = -m_Zsize / 2. + m_firstwire_y[tubeLayer - 1];
            ztube = -m_Rsize / 2. + m_firstwire_x[tubeLayer - 1] + (tube - 1) * m_tubepitch;
        }
        return Amg::Vector3D{xtube, 0., ztube};
    }

localROPos() [1/2]

Amg::Vector3D MuonGM::MdtReadoutElement::localROPos

(

const Identifier &

id

)

const

localROPos() [2/2]

Amg::Vector3D MuonGM::MdtReadoutElement::localROPos	(	const int	tubelayer,
		const int	tube ) const

Definition at line 285 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                       {
        return tubeToMultilayerTransf(tubeLayer, tube) * tubeFrame_localROPos(tubeLayer, tube);
    }

localToGlobalTransf() [1/2]

const Amg::Transform3D & MuonGM::MdtReadoutElement::localToGlobalTransf

(

const Identifier &

id

)

const

localToGlobalTransf() [2/2]

const Amg::Transform3D & MuonGM::MdtReadoutElement::localToGlobalTransf	(	const int	tubeLayer,
		const int	tube ) const

localTubePos() [1/2]

Amg::Vector3D MuonGM::MdtReadoutElement::localTubePos

(

const Identifier &

id

)

const

localTubePos() [2/2]

Amg::Vector3D MuonGM::MdtReadoutElement::localTubePos	(	const int	tubelayer,
		const int	tube ) const

Definition at line 291 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                         {
        return fromIdealToDeformed(tubeLayer, tube) * nodeform_localTubePos(tubeLayer, tube);
    }

makeGeoInfo()

std::unique_ptr< MdtReadoutElement::GeoInfo > MuonGM::MdtReadoutElement::makeGeoInfo ( const int tubelayer, const int tube ) const

private

Definition at line 802 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                                    {
        Amg::Transform3D transform = globalTransform(nodeform_localTubePos(tubeLayer, tube), fromIdealToDeformed(tubeLayer, tube));
        return std::make_unique<GeoInfo>(std::move(transform));
    }

manager()

const MuonDetectorManager * MuonGM::MuonReadoutElement::manager ( ) const

inlineprotectedinherited

Definition at line 193 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return m_muon_mgr; }

msg() [1/2]

MsgStream & AthMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 163 of file AthMessaging.h.

{
  MsgStream* ms = m_msg_tls.get();
  if (!ms) {
    if (!m_initialized.test_and_set()) initMessaging();
    ms = new MsgStream(m_imsg,m_nm);
    m_msg_tls.reset( ms );
  }
 
  ms->setLevel (m_lvl);
  return *ms;
}

msg() [2/2]

MsgStream & AthMessaging::msg ( const MSG::Level lvl ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 178 of file AthMessaging.h.

{ return msg() << lvl; }

msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns

boolean Indicating if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

{
  if (m_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

nMDTinStation()

unsigned int MuonGM::MdtReadoutElement::nMDTinStation ( ) const

inline

How many MDT chambers are in the station.

Definition at line 62 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

                                           { 
             return m_nMDTinStation;
        }

nodeform_globalToLocalTransf() [1/2]

Amg::Transform3D MuonGM::MdtReadoutElement::nodeform_globalToLocalTransf

(

const Identifier &

id

)

const

nodeform_globalToLocalTransf() [2/2]

Amg::Transform3D MuonGM::MdtReadoutElement::nodeform_globalToLocalTransf	(	const int	tubeLayer,
		const int	tube ) const

Definition at line 404 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                            {
        return nodeform_localToGlobalTransf(tubeLayer, tube).inverse();
    }

nodeform_localToGlobalTransf() [1/2]

Amg::Transform3D MuonGM::MdtReadoutElement::nodeform_localToGlobalTransf

(

const Identifier &

id

)

const

nodeform_localToGlobalTransf() [2/2]

Amg::Transform3D MuonGM::MdtReadoutElement::nodeform_localToGlobalTransf	(	const int	tubeLayer,
		const int	tube ) const

Definition at line 397 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                            {
        return globalTransform(nodeform_localTubePos(tubeLayer, tube), Amg::Transform3D::Identity());
    }

nodeform_localTubePos() [1/2]

Amg::Vector3D MuonGM::MdtReadoutElement::nodeform_localTubePos

(

const Identifier &

id

)

const

nodeform_localTubePos() [2/2]

Amg::Vector3D MuonGM::MdtReadoutElement::nodeform_localTubePos	(	const int	tubelayer,
		const int	tube ) const

Definition at line 295 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                  {
        ATH_MSG_VERBOSE( " Computing LocalTubePos for "<< idHelperSvc()->toStringDetEl(identify())
                                                       << "/" << tubeLayer << "/" << tube );
        double xtube{0.}, ytube{0.}, ztube{0.};
        if (barrel()) {
            xtube = -m_Rsize / 2. + m_firstwire_y[tubeLayer - 1];
            ztube = -m_Zsize / 2. + m_firstwire_x[tubeLayer - 1] + (tube - 1) * m_tubepitch;
        } else {
            xtube = -m_Zsize / 2. + m_firstwire_y[tubeLayer - 1];
            ztube = -m_Rsize / 2. + m_firstwire_x[tubeLayer - 1] + (tube - 1) * m_tubepitch;
        }
        Amg::Vector3D tubePos{xtube, ytube, ztube};
        if (hasCutouts()) {
            if (manager()->MinimalGeoFlag() == 0) {
                ATH_MSG_DEBUG( " MdtReadoutElement " << idHelperSvc()->toStringDetEl(identify()) << " has cutouts, check for real position of tubes " );
                PVConstLink cv = getMaterialGeom();  // it is "Multilayer"
                int nGrandchildren = cv->getNChildVols();
                // child vol 0 is foam; 1 to (nGrandchildren-1) should be tubes
                int ii = (tubeLayer - 1) * m_ntubesperlayer + tube;
 
                // BIS78 only (the BIS7 of Run1/2 has no cutouts, thus, this block won't be reached)
                if ((getStationIndex() == m_stIdx_BIS && std::abs(getStationEta()) == 7)) --ii;
                if (m_idHelper.isBMG(identify())) {
                    // usually the tube number corresponds to the child number, however for
                    // BMG chambers full tubes are skipped during the building process
                    // therefore the matching needs to be done via the volume ID
                    int packed_id = tube + maxNTubesPerLayer * tubeLayer;
                    int kk = 0;
                    bool found = false;
                    geoGetIds(
                        [&](int id) {
                            if (!found && id == packed_id) {
                                ii = kk;
                                found = true;
                            }
                            ++kk;
                        },
                        &*cv);
                    if (found) {
                        ATH_MSG_DEBUG( " MdtReadoutElement tube match found for BMG - input : tube(" << tube << "), layer("
                            << tubeLayer << ") - output match : tube(" << ii % maxNTubesPerLayer << "), layer(" << ii / maxNTubesPerLayer
                            << ")" );
                    }
                }
                GeoTrf::Transform3D tubeTrans = cv->getXToChildVol(ii);
                PVConstLink tv = cv->getChildVol(ii);
                constexpr double maxtol = 0.0000001;
 
                if (std::abs(xtube - tubeTrans(0, 3)) > maxtol || std::abs(ztube - tubeTrans(2, 3)) > maxtol) {
                    THROW_EXCEPTION(__func__<<"("<<tubeLayer<<","<<tube<<")"<<" - mismatch between local from tube-id/pitch/cutout position"<<
                           Amg::toString(tubePos)<<" and GeoModel "<< Amg::toString(tubeTrans.linear().col(3))<<" for detector element "<<
                            idHelperSvc()->toStringDetEl(identify()) <<"There are "<<nGrandchildren<<" child volumes and "<<
                                       (m_ntubesperlayer * m_nlayers)<<" are expected. There should be "<<m_nlayers<<" and "<<
                                           m_ntubesperlayer<<" tubes per layer");
                }
                if (tubeTrans(1, 3) > maxtol) {
 
                    ATH_MSG_DEBUG( "This a tube with cutout stName/Eta/Phi/ml/tl/t = " << idHelperSvc()->toStringDetEl(identify())
                                         << "/" << tubeLayer << "/" << tube);
                    // check only for tubes actually shifted
                    if (std::abs(m_cutoutShift - tubeTrans(1, 3)) > maxtol) {
                        THROW_EXCEPTION(__func__<<"("<<tubeLayer<<","<<tube<<")"<<" - mismatch between local from tube-id/pitch/cutout position"<<
                           Amg::toString(tubePos)<<" and GeoModel "<< Amg::toString(tubeTrans.linear().col(3))<<" for detector element "<<
                            idHelperSvc()->toStringDetEl(identify()) <<"There are "<<nGrandchildren<<" child volumes and "<<
                                       (m_ntubesperlayer * m_nlayers)<<" are expected. There should be "<<m_nlayers<<" and "<<
                                           m_ntubesperlayer<<" tubes per layer");
                    }
                }
 
                Amg::Vector3D x = tubeTrans.translation();
                if (tube > m_ntubesperlayer || tubeLayer > m_nlayers) { x = tubePos; }
                return x;
            }
        }
        return tubePos;
    }

nodeform_tubePos() [1/2]

Amg::Vector3D MuonGM::MdtReadoutElement::nodeform_tubePos

(

const Identifier &

id

)

const

Returns the global position of the tube excluding the B-line & As-built corrections.

nodeform_tubePos() [2/2]

Amg::Vector3D MuonGM::MdtReadoutElement::nodeform_tubePos	(	const int	tubelayer,
		const int	tube ) const

Definition at line 372 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                 {
        return transform() * nodeform_localTubePos(tubeLayer, tube);
    }

normal() [1/2]

const Amg::Vector3D & MuonGM::MdtReadoutElement::normal ( ) const

finaloverridevirtual

Return the normal of the element.

Implements Trk::TrkDetElementBase.

Definition at line 880 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                       {
        if (!m_elemNormal.isValid()) { m_elemNormal.set(transform().linear() * Amg::Vector3D::UnitX()); }
        return *m_elemNormal.ptr();
    }

normal() [2/2]

virtual const Amg::Vector3D & MuonGM::MdtReadoutElement::normal ( const Identifier & id ) const

finaloverridevirtual

Return the normal of the surface associated with this identifier In the case of silicon it returns the same as normal()

Implements Trk::TrkDetElementBase.

outerTubeRadius()

double MuonGM::MdtReadoutElement::outerTubeRadius

(

)

const

Returns the tube radius taking the thickness of the tubes into account.

parentMuonStation()

const MuonStation * MuonGM::MuonReadoutElement::parentMuonStation ( ) const

inherited

Definition at line 80 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

{ return m_parentMuonStation; }

parentMuonStationPos()

Amg::Vector3D MuonGM::MuonReadoutElement::parentMuonStationPos ( ) const

inherited

Definition at line 82 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

                                                               {
        return parentMuonStation()->getTransform().translation();       
    }

posOnDefChamWire()

Amg::Vector3D MuonGM::MdtReadoutElement::posOnDefChamWire ( const Amg::Vector3D & locAMDBPos, const double width_narrow, const double width_wide, const double height, const double thickness, const Amg::Vector3D & fixedPoint ) const

private

Definition at line 578 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                           {
 
        using Parameter = BLinePar::Parameter;
        const double sp = m_BLinePar->getParameter(Parameter::sp);
        const double sn = m_BLinePar->getParameter(Parameter::sn);
        const double tw = m_BLinePar->getParameter(Parameter::tw);
        const double eg = m_BLinePar->getParameter(Parameter::eg);
        double ep = m_BLinePar->getParameter(Parameter::ep);
        double en = m_BLinePar->getParameter(Parameter::en);
        // S.Spagnolo Feb.6, 2011, modified by P.F Giraud, 2015-01-17
        // This version does not implement deformations modifying only the second
        // coordinate, or leaving the end-plug positions unchanged.
        // This version should be called only with the coordinates of the end-plugs
        // as argument, as is done in fromIdealToDeformed
 
        // MDT deformations like tube bow: bz,bp,bn bend the tube while the wire endpoints are not affected
        //                                 => the wire keeps it's nominal straight line trajectory but it is not concentric to the tube
        //                                 ==> in this function bz, bp and bn are ignored or set to 0
        // MDT deformations that extend/contract the wire longitudinally (while keeping it straight):
        //                                 delta_s from eg and tr are irrelevant for the tube geometry
        //                                 and for the wire trajectory => set to 0 here
        //                                 (should be applied as a correction to the
        //                                 tube lenght => tube surface bounds
        //                                 =++>>>> effect in tracking just through the gravitational sagging TOTALLY NEGLIGIBLE=> ignore)
        // pg is irrelevant for tracking purposes and (at least for the endcaps) is applies to the internal bars only, not to the tubes !!!
        //                                 =++>>>> IGNORE IT
        // ep,en: bend the tube by moving (differently) the endplugs ===> the wire straight trajectory is moved w.r.t. the nominal one
        //                                 in addition the tubes keep their nominal position at the center => the wire is not concentric to
        //                                 the tube delta_s from ep,en must also be considered for the implementation of the realistic tube
        //                                 trajectory induced by ep,en
        // tw,sp,sn,pg (for deltaT and deltaZ) are geometrical effects, that impact on tracking and keep the wire straight.
 
 
 
        Amg::Vector3D deformPos(locAMDBPos);
 
        // NOTE s0,z0,t0 are the coord. in the amdb frame of this point: the origin of the frame can be different than the fixed point for
        // deformations s0mdt,z0mdt,t0mdt
        //    (always equal to the point at lowest t,z and s=0 of the MDT stack)
        ATH_MSG_VERBOSE( "** In "<<__func__<<" - moduleWidthS " << moduleWidthS<<", moduleWidthL: "
                       <<moduleWidthL <<", height: "<<height<<", thickness: " <<thickness << "." );
        ATH_MSG_VERBOSE( "** In "<<__func__<<" - going to correct for B-line the position of Point at " <<Amg::toString(locAMDBPos)
                       << " in the amdb-szt frame" );
 
        double s0mdt = locAMDBPos.x();  // always I think !
        if (std::abs(fixedPoint.x()) > 0.01) {
            s0mdt = locAMDBPos.x() - fixedPoint.x();
        }
        double z0mdt = locAMDBPos.y();  
        // unless in the D section of this station there's a dy diff. from 0 for the innermost MDT multilayer (sometimes in the barrel)
        if (std::abs(fixedPoint.y()) > 0.01) {
            z0mdt = locAMDBPos.y() - fixedPoint.y();
        }
        double t0mdt = locAMDBPos.z();
              // unless in the D section of this station there's a dz diff. from 0 for the innermost MDT multilayer (often in barrel)
        if (std::abs(fixedPoint.z()) > 0.01) t0mdt = locAMDBPos.z() - fixedPoint.z();
        if (z0mdt < 0 || t0mdt < 0) {
            ATH_MSG_WARNING(""<<__func__<<": correcting the local position of a point outside the mdt station (2 multilayers) volume -- RE "
                << idHelperSvc()->toStringDetEl(identify()) << " local point: szt=" << Amg::toString(locAMDBPos)
                << " fixedPoint " <<Amg::toString(fixedPoint) );
        }
        ATH_MSG_VERBOSE( "** In "<<__func__<<" - correct for offset of B-line fixed point " << s0mdt << " " << z0mdt << " " << t0mdt);
 
        double ds{0.},dz{0.},dt{0.};
        double width_actual = moduleWidthS + (moduleWidthL - moduleWidthS) * (z0mdt / height);
        double s_rel = s0mdt / (width_actual / 2.);
        double z_rel = (z0mdt - height / 2.) / (height / 2.);
        double t_rel = (t0mdt - thickness / 2.) / (thickness / 2.);
 
        ATH_MSG_VERBOSE( "** In "<<__func__<<" - width_actual: "<<width_actual<<", s_rel: "<<s_rel<<", z_rel: "
                    <<z_rel<<", t_rel: "<<t_rel);
 
        // sp, sn - cross plate sag out of plane
        if ((sp != 0) || (sn != 0)) {
            double ztmp = z_rel * z_rel - 1;
            dt += 0.5 * (sp + sn) * ztmp + 0.5 * (sp - sn) * ztmp * s_rel;
        }
 
        // tw     - twist
        if (tw != 0) {
            dt -= tw * s_rel * z_rel;
            dz += tw * s_rel * t_rel * thickness / height;
            ATH_MSG_VERBOSE( "** In "<<__func__<<": tw=" << tw << " dt, dz " << dt << " " << dz );
        }
 
        // eg     - global expansion
        if (eg != 0) {
            double egppm = eg / 1000.;
            ds += 0.;
            dz += z0mdt * egppm;
            dt += t0mdt * egppm;
        }
 
        // ep, en - local expansion
        //
        // Imporant note: the chamber height and length, as they denoted in Christoph's talk,
        // correspond to thickness and height parameters of this function;
        //
 
        if ((ep != 0) || (en != 0)) {
            ep = ep / 1000.;
            en = en / 1000.;
            double phi = 0.5 * (ep + en) * s_rel * s_rel + 0.5 * (ep - en) * s_rel;
            double localDt = phi * (t0mdt - thickness / 2.);
            double localDz = phi * (z0mdt - height / 2.);
            dt += localDt;
            dz += localDz;
        }
 
        ATH_MSG_VERBOSE( "posOnDefChamStraighWire: ds,z,t = " << ds << " " << dz << " " << dt );
        deformPos[0] = locAMDBPos[0] + ds;
        deformPos[1] = locAMDBPos[1] + dz;
        deformPos[2] = locAMDBPos[2] + dt;
 
        return deformPos;
    }

refreshCache()

void MuonGM::MuonReadoutElement::refreshCache ( )

inherited

Definition at line 122 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

                                          {
        clearCache();
        fillCache();
    }

RODistanceFromTubeCentre() [1/2]

double MuonGM::MdtReadoutElement::RODistanceFromTubeCentre

(

const Identifier &

id

)

const

RODistanceFromTubeCentre() [2/2]

double MuonGM::MdtReadoutElement::RODistanceFromTubeCentre	(	const int	tubeLayer,
		const int	tube ) const

Definition at line 221 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                {
        return getWireLength(tubeLayer, tube) / 2.;
    }

ROPos() [1/2]

Amg::Vector3D MuonGM::MdtReadoutElement::ROPos

(

const Identifier &

id

)

const

ROPos() [2/2]

Amg::Vector3D MuonGM::MdtReadoutElement::ROPos	(	const int	tubelayer,
		const int	tube ) const

Definition at line 288 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                  {
        return transform(tubeLayer, tube) * tubeFrame_localROPos(tubeLayer, tube);
    }

setBLinePar()

void MuonGM::MdtReadoutElement::setBLinePar

(

const BLinePar *

bLine

)

Definition at line 952 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                             {
        ATH_MSG_DEBUG( "Setting B-line for " << idHelperSvc()->toStringDetEl(identify()) );
        if (m_BLinePar == bLine) {
            return;
        }
        m_BLinePar = bLine;
        refreshCache();
    }

setHasCutouts()

void MuonGM::MuonReadoutElement::setHasCutouts ( bool flag )

inlineinherited

Definition at line 121 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ m_hasCutouts = flag; }

setIdentifier()

void MuonGM::MuonReadoutElement::setIdentifier ( const Identifier & id )

inherited

Sets the Identifier, hashes & station names.

Definition at line 101 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

                                                               {
        m_id = id;
        if (!m_idHelperSvc->isMuon(id)) {
            ATH_MSG_FATAL("The Identifier "<<m_idHelperSvc->toString(id)<<" is not a muon one.");
            throw std::runtime_error("Invalid Identifier set");
        }
        m_stIdx = m_idHelperSvc->stationName(id);
        m_eta = m_idHelperSvc->stationEta(id);
        m_phi = m_idHelperSvc->stationPhi(id);
        m_idhash = m_idHelperSvc->moduleHash(id);
        m_detectorElIdhash = m_idHelperSvc->detElementHash(id);
    }

setLevel()

void AthMessaging::setLevel ( MSG::Level lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

{
  m_lvl = lvl;
}

setLongRsize()

void MuonGM::MuonReadoutElement::setLongRsize ( double v )

inherited

Definition at line 117 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

{ m_LongRsize = v; }

setLongSsize()

void MuonGM::MuonReadoutElement::setLongSsize ( double v )

inherited

Definition at line 116 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

{ m_LongSsize = v; }

setLongZsize()

void MuonGM::MuonReadoutElement::setLongZsize ( double v )

inherited

Definition at line 118 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

{ m_LongZsize = v; }

setMultilayer()

void MuonGM::MdtReadoutElement::setMultilayer

(

const int

ml

)

Sets the multilayer number.

Definition at line 67 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

{ m_multilayer = ml; }

setNLayers()

void MuonGM::MdtReadoutElement::setNLayers

(

const int

nl

)

Sets the number of layers.

Definition at line 68 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

{ m_nlayers = nl; }

setNMdtInStation()

void MuonGM::MdtReadoutElement::setNMdtInStation ( unsigned int numMdt )

inline

Definition at line 65 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{ m_nMDTinStation = numMdt;}

setParentMuonStation()

void MuonGM::MuonReadoutElement::setParentMuonStation ( const MuonStation * mstat )

inherited

Definition at line 76 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

                                                                          { 
        m_parentMuonStation = mstat; 
    }

setRsize()

void MuonGM::MuonReadoutElement::setRsize ( double v )

inherited

Definition at line 120 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

{ m_Rsize = v; }

setSsize()

void MuonGM::MuonReadoutElement::setSsize ( double v )

inherited

Definition at line 119 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

{ m_Ssize = v; }

setStationName()

void MuonGM::MuonReadoutElement::setStationName ( const std::string & str )

inherited

Definition at line 114 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

{ m_statname = str; }

setStationS()

void MuonGM::MuonReadoutElement::setStationS ( double v )

inherited

Definition at line 115 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

{ m_stationS = v; }

setTechnologyName()

void MuonGM::MuonReadoutElement::setTechnologyName ( const std::string & str )

inherited

Definition at line 113 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

{ m_techname = str; }

setZsize()

void MuonGM::MuonReadoutElement::setZsize ( double v )

inherited

Definition at line 121 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

{ m_Zsize = v; }

sideA()

bool MuonGM::MuonReadoutElement::sideA ( ) const

inlineinherited

Definition at line 196 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return (getStationEta() > 0); }

sideC()

bool MuonGM::MuonReadoutElement::sideC ( ) const

inlineinherited

Definition at line 197 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{ return (getStationEta() < 0); }

signedRODistanceFromTubeCentre() [1/2]

double MuonGM::MdtReadoutElement::signedRODistanceFromTubeCentre

(

const Identifier &

id

)

const

signedRODistanceFromTubeCentre() [2/2]

double MuonGM::MdtReadoutElement::signedRODistanceFromTubeCentre	(	const int	tubeLayer,
		const int	tube ) const

Definition at line 224 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                      {
        // it is a signed quantity:
        // the sign corresponds to the sign of the z coordinate of the RO endplug in the tube
        // reference frame
        int amdb_plus_minus1 = 1;
        if (!m_zsignRO_tubeFrame.isValid()) {
            const MuonStation* ms = parentMuonStation();
            if (std::abs(ms->xAmdbCRO()) > 10.) {
                Amg::Vector3D tem = ms->xAmdbCRO()* Amg::Vector3D::UnitX();
                Amg::Transform3D amdbToGlobal{ms->getAmdbLRSToGlobal()};
                Amg::Vector3D temGlo = amdbToGlobal * tem;
                Amg::Vector3D ROtubeFrame = nodeform_globalToLocalTransf(tubeLayer, tube) * temGlo;
                if (ROtubeFrame.z() < 0)
                    m_zsignRO_tubeFrame.set(-1);
                else
                    m_zsignRO_tubeFrame.set(1);
            }
        }
        // if no CRO in a chamber in AMDB (BIS in layout R), use the standard convention for RO-HV side
        if (!m_zsignRO_tubeFrame.isValid()) {  
            int sign = 0;
            if (barrel()) {
                if (sideA()) {
                    if (largeSector())
                        sign = -1;
                    else
                        sign = 1;
                } else {
                    if (largeSector())
                        sign = 1;
                    else
                        sign = -1;
                }
                // a special case is BIS in sector 12
                if (getStationName().substr(0, 3) == "BIS" && getStationPhi() == 6) sign = -sign;
            } else {
                if (sideA()) {
                    if (largeSector())
                        sign = 1;
                    else
                        sign = -1;
                } else {
                    if (largeSector())
                        sign = -1;
                    else
                        sign = 1;
                }
            }
            m_zsignRO_tubeFrame.set(sign);
        }
        amdb_plus_minus1 = *m_zsignRO_tubeFrame.ptr();
        if (amdb_plus_minus1 == 0) {
            ATH_MSG_WARNING( "Unable to get the sign of RO side; signedRODistancefromTubeCenter returns 0" );
        }
 
        return amdb_plus_minus1 * getWireLength(tubeLayer, tube) / 2.;
    }

smallSector()

bool MuonGM::MuonReadoutElement::smallSector ( ) const

inherited

Definition at line 67 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

                                               {
        // this doesn't apply to TGC
        return (!largeSector());
    }

surface() [1/3]

const Trk::Surface & MuonGM::MdtReadoutElement::surface ( ) const

finaloverridevirtual

Return surface associated with this detector element.

Implements Trk::TrkDetElementBase.

Definition at line 885 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                       {
        if (!m_associatedSurface) {
            Amg::RotationMatrix3D muonTRotation(transform().rotation());
            Amg::RotationMatrix3D surfaceTRotation;
            surfaceTRotation.col(0) = muonTRotation.col(1);
            surfaceTRotation.col(1) = muonTRotation.col(2);
            surfaceTRotation.col(2) = muonTRotation.col(0);
 
            Amg::Transform3D trans3D(surfaceTRotation);
            trans3D.pretranslate(transform().translation());
 
            if (barrel()) {
                m_associatedSurface.set(std::make_unique<Trk::PlaneSurface>(Amg::Transform3D(trans3D), getSsize() * 0.5,
                                                                            getZsize() * 0.5));
            } else {
                m_associatedSurface.set(std::make_unique<Trk::PlaneSurface>(Amg::Transform3D(trans3D), getSsize() * 0.5,
                                                                            getLongSsize() * 0.5,
                                                                            getRsize() * 0.5));
            }
        }
        return *m_associatedSurface;
    }

surface() [2/3]

virtual const Trk::StraightLineSurface & MuonGM::MdtReadoutElement::surface ( const Identifier & id ) const

finaloverridevirtual

Return surface associated with this identifier, which should come from the PrepRawData object (i.e.

Surface s = PRD.detElement().surface( PRD.identify() ). This is only really relevant for the TRT (where there are several surfaces per detector element). For other detector elements it will just return surface() ... the Identifier is ignored.

Implements Trk::TrkDetElementBase.

surface() [3/3]

const Trk::StraightLineSurface & MuonGM::MdtReadoutElement::surface ( const int tubeLayer, const int tube ) const

virtual

Definition at line 829 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                      {
 
        int ntot_tubes = m_nlayers * m_ntubesperlayer;
        int itube = (tubeLayer - 1) * m_ntubesperlayer + tube - 1;
        // consistency checks
        if (itube >= ntot_tubes) {
            ATH_MSG_WARNING( "surface called with tubeLayer or tube out of range in chamber "
                << idHelperSvc()->toStringDetEl(identify()) << " : layer " << tubeLayer << " max " << m_nlayers << " tube " << tube
                << " max " << m_ntubesperlayer << " will compute surface for first tube in this chamber" );
            ATH_MSG_WARNING( "Please run in DEBUG mode to get extra diagnostic" );
            itube = 0;
        }
 
        const CxxUtils::CachedUniquePtr<Trk::StraightLineSurface>& ptr = m_tubeSurfaces.at(itube);
        if (!ptr) {
            Identifier id = m_idHelper.channelID(identify(), getMultilayer(), tubeLayer, tube);
            ptr.set(std::make_unique<Trk::StraightLineSurface>(*this, id));
            if (!m_haveTubeSurfaces) m_haveTubeSurfaces = true;
        }
        return *ptr;
    }

surfaces()

std::vector< const Trk::Surface * > MuonGM::MdtReadoutElement::surfaces

(

)

const

returns all the surfaces contained in this detector element

Definition at line 70 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                   {
        std::vector<const Trk::Surface*> elementSurfaces;
        elementSurfaces.reserve(m_tubeSurfaces.size() + 1);
        if (m_associatedSurface) { elementSurfaces.push_back(m_associatedSurface.get()); }
        for (const auto& s : m_tubeSurfaces) {
            if (s) elementSurfaces.push_back(s.get());
        }
        return elementSurfaces;
    }

toParentStation()

Amg::Transform3D MuonGM::MuonReadoutElement::toParentStation ( ) const

inherited

Definition at line 72 of file MuonDetDescr/MuonReadoutGeometry/src/MuonReadoutElement.cxx.

                                                             {
        return getMaterialGeom()->getX();
    }

transform() [1/3]

virtual const Amg::Transform3D & MuonGM::MdtReadoutElement::transform ( ) const

finaloverridevirtual

Return local to global transform.

Implements Trk::TrkDetElementBase.

transform() [2/3]

virtual const Amg::Transform3D & MuonGM::MdtReadoutElement::transform ( const Identifier & id ) const

finaloverridevirtual

Return local to global transform associated with this identifier.

Implements Trk::TrkDetElementBase.

transform() [3/3]

const Amg::Transform3D & MuonGM::MdtReadoutElement::transform	(	const int	tubeLayer,
		const int	tube ) const

Definition at line 825 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                {
        return geoInfo(tubeLayer, tube).m_transform;
    }

tubeFrame_localROPos() [1/2]

Amg::Vector3D MuonGM::MdtReadoutElement::tubeFrame_localROPos

(

const Identifier &

id

)

const

tubeFrame_localROPos() [2/2]

Amg::Vector3D MuonGM::MdtReadoutElement::tubeFrame_localROPos	(	const int	tubelayer,
		const int	tube ) const

Definition at line 282 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                 {
        return signedRODistanceFromTubeCentre(tubeLayer, tube) * Amg::Vector3D::UnitZ();
    }

tubeLength() [1/2]

double MuonGM::MdtReadoutElement::tubeLength

(

const Identifier &

id

)

const

tubeLength() [2/2]

double MuonGM::MdtReadoutElement::tubeLength	(	const int	tubeLayer,
		const int	tube ) const

tubePitch()

double MuonGM::MdtReadoutElement::tubePitch

(

)

const

Returns the distance between 2 tubes in a tube layer.

tubePos() [1/2]

Amg::Vector3D MuonGM::MdtReadoutElement::tubePos

(

const Identifier &

id

)

const

Returns the global position of the given tube.

tubePos() [2/2]

Amg::Vector3D MuonGM::MdtReadoutElement::tubePos	(	const int	tubelayer,
		const int	tube ) const

Definition at line 375 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                        {
      return transform() * localTubePos(tubeLayer, tube);
    }

tubeToMultilayerTransf() [1/3]

Amg::Transform3D MuonGM::MdtReadoutElement::tubeToMultilayerTransf ( const Amg::Vector3D & tubePos, const Amg::Transform3D & toDeform )

staticprivate

Definition at line 394 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                                       {
        return toDeform * Amg::Translation3D{tubePos} * Amg::getRotateX3D(90. * CLHEP::deg);
    }

tubeToMultilayerTransf() [2/3]

Amg::Transform3D MuonGM::MdtReadoutElement::tubeToMultilayerTransf

(

const Identifier &

id

)

const

tubeToMultilayerTransf() [3/3]

Amg::Transform3D MuonGM::MdtReadoutElement::tubeToMultilayerTransf	(	const int	tubeLayer,
		const int	tube ) const

Definition at line 378 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                                      {
        return tubeToMultilayerTransf(nodeform_localTubePos(tubeLayer, tube),
                                      fromIdealToDeformed(tubeLayer, tube));
    }

wireEndpointsAsBuilt()

void MuonGM::MdtReadoutElement::wireEndpointsAsBuilt ( Amg::Vector3D & locAMDBWireEndP, Amg::Vector3D & locAMDBWireEndN, const int tubelayer, const int tube ) const

private

Definition at line 712 of file MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx.

                                                                                       {
        const MdtAsBuiltPar* params = parentMuonStation()->getMdtAsBuiltParams();
        if (!params) {
            ATH_MSG_WARNING( "Cannot find Mdt AsBuilt parameters for chamber " << idHelperSvc()->toStringDetEl(identify()) );
            return;
        }
 
        ATH_MSG_VERBOSE( "Applying as-built parameters for chamber " << idHelperSvc()->toStringDetEl(identify())
                        << " tubeLayer=" << tubeLayer << " tube=" << tube);
 
        constexpr int nsid = 2;
        using tubeSide_t = MdtAsBuiltPar::tubeSide_t;
        using multilayer_t = MdtAsBuiltPar::multilayer_t;
        std::array<MdtAsBuiltPar::tubeSide_t, nsid> tubeSide{tubeSide_t::POS, tubeSide_t::NEG};
        std::array<Amg::Vector3D, nsid> wireEnd{locAMDBWireEndP, locAMDBWireEndN};
        multilayer_t ml = (getMultilayer() == 1) ? multilayer_t::ML1 : multilayer_t::ML2;
 
        const double xmin = *std::min_element(m_firstwire_x.begin(), m_firstwire_x.begin() + m_nlayers) - outerTubeRadius();
        const int ref_layer = (ml == multilayer_t::ML1) ? m_nlayers : 1;
        const double y_offset = (ml == multilayer_t::ML1) ? outerTubeRadius() : -outerTubeRadius();
 
        for (int isid = 0; isid < nsid; ++isid) {  // first s>0 then s<0
            // Compute the reference for the as-built parameters
            double xref{0.}, yref{0.}, zref{0.};
            if (barrel()) {
                xref = -m_Rsize / 2. + m_firstwire_y[ref_layer - 1] + y_offset;
                zref = -m_Zsize / 2. + xmin;
            } else {
                xref = -m_Zsize / 2. + m_firstwire_y[ref_layer - 1] + y_offset;
                zref = -m_Rsize / 2. + xmin;
            }
            const Amg::Transform3D toAMDB = parentMuonStation()->getNativeToAmdbLRS() * toParentStation();
            Amg::Vector3D reference_point{xref, yref, zref};
 
            reference_point = toAMDB * reference_point;
 
            if (isid == 0)
                reference_point = reference_point + 0.5 * getNominalTubeLengthWoCutouts(ref_layer, 1) * Amg::Vector3D::UnitX();
            else
                reference_point = reference_point -0.5 * getNominalTubeLengthWoCutouts(ref_layer, 1) * Amg::Vector3D::UnitX();
 
            ATH_MSG_VERBOSE("AMDB transform "<<idHelperSvc()->toStringDetEl(identify())<<
                            " "<<GeoTrf::toString(toAMDB, true)<<", reference point: "<<Amg::toString(reference_point));
 
            int layer_delta = tubeLayer;
            if (ml == multilayer_t::ML1) layer_delta = m_nlayers + 1 - tubeLayer;
 
            // Get the As-Built parameters for this ML and side of the chamber
            double zpitch = params->zpitch(ml, tubeSide[isid]);
            double ypitch = params->ypitch(ml, tubeSide[isid]);
            double stagg = (double)params->stagg(ml, tubeSide[isid]);
            double alpha = params->alpha(ml, tubeSide[isid]);
            double y0 = params->y0(ml, tubeSide[isid]);
            double z0 = params->z0(ml, tubeSide[isid]);
 
            // Find the vector from the reference_point to the endplug
            int do_stagg = (layer_delta - 1) % 2;  // 0 for layer 1 and 3, 1 for layer 2 and 4
            double offset_stagg = ((double)do_stagg) * 0.5 * zpitch * stagg;
            Amg::Vector3D end_plug(0., (tube - 1.0) * zpitch + offset_stagg, (layer_delta - 1) * ypitch);
 
            const Amg::Transform3D amgTransf{Amg::getRotateX3D(-alpha)};
            end_plug = amgTransf * end_plug;
 
            // Calculate x position, which varies for endcap chambers
            double xshift = 0.;
            if (isid == 0)
                xshift = 0.5 * getNominalTubeLengthWoCutouts(tubeLayer, tube) - 0.5 * getNominalTubeLengthWoCutouts(ref_layer, 1);
            else
                xshift = -0.5 * getNominalTubeLengthWoCutouts(tubeLayer, tube) + 0.5 * getNominalTubeLengthWoCutouts(ref_layer, 1);
 
            Amg::Vector3D ret(reference_point.x() + xshift, reference_point.y() + z0 + end_plug.y(),
                              reference_point.z() + y0 + end_plug.z());
 
            // Warn if result of calculation is too far off
            // BIL1A13 has as-built parameters up to 3mm off, giving the size of the cut
            if ((ret - wireEnd[isid]).mag() > 3. * CLHEP::mm) {
                ATH_MSG_WARNING( "Large as-built correction for chamber " << idHelperSvc()->toStringDetEl(identify()) << ", side "
                    << isid << ", Delta " << Amg::toString(ret - wireEnd[isid]) );
            }
            ATH_MSG_VERBOSE(( tubeSide[isid] == tubeSide_t::POS ? "positive" : "negative")<<" wire end has moved from "
                        <<Amg::toString(wireEnd[isid])<<" to "<<Amg::toString(ret));
 
            // Save the result
            if (tubeSide[isid] == tubeSide_t::POS)
                locAMDBWireEndP = ret;
            else
                locAMDBWireEndN = ret;
        }
    }

Friends And Related Symbol Documentation

MuonChamber

friend class MuonChamber

friend

Definition at line 53 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

MuonChamberLite

friend class MuonChamberLite

friend

Definition at line 54 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

MuonGMR4::ReadoutGeomCnvAlg

friend class MuonGMR4::ReadoutGeomCnvAlg

friend

Definition at line 52 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

Member Data Documentation

ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging)

Definition at line 141 of file AthMessaging.h.

m_associatedBounds

CxxUtils::CachedUniquePtr<Trk::SurfaceBounds> MuonGM::MdtReadoutElement::m_associatedBounds {}

private

Definition at line 272 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{};

m_associatedSurface

CxxUtils::CachedUniquePtr<Trk::Surface> MuonGM::MdtReadoutElement::m_associatedSurface {}

private

Definition at line 271 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{};

m_backupTubeGeo

std::vector<CxxUtils::CachedUniquePtr<GeoInfo> > MuonGM::MdtReadoutElement::m_backupTubeGeo {}

private

Definition at line 253 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{};                // one per tube

m_BLinePar

const BLinePar* MuonGM::MdtReadoutElement::m_BLinePar {nullptr}

private

Definition at line 256 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{nullptr};

m_caching

int MuonGM::MuonReadoutElement::m_caching {-1}

protectedinherited

0 if we want to avoid caching geometry info for tracking interface

Definition at line 150 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{-1};

m_cutoutShift

double MuonGM::MdtReadoutElement::m_cutoutShift {-9999.}

private

Definition at line 224 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{-9999.};

m_deadlength

double MuonGM::MdtReadoutElement::m_deadlength {-9999.}

private

Definition at line 226 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{-9999.};

m_deformTransf

std::vector<CxxUtils::CachedUniquePtr<Amg::Transform3D> > MuonGM::MdtReadoutElement::m_deformTransf {}

private

Definition at line 254 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{};        // one per tube

m_detectorElIdhash

IdentifierHash MuonGM::MuonReadoutElement::m_detectorElIdhash {0}

privateinherited

detector element hash identifier

Definition at line 159 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{0};  

m_elemNormal

CxxUtils::CachedValue<Amg::Vector3D> MuonGM::MdtReadoutElement::m_elemNormal {}

private

Definition at line 257 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{};                               // one

m_endpluglength

double MuonGM::MdtReadoutElement::m_endpluglength {-9999.}

private

Definition at line 225 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{-9999.};

m_eta

int MuonGM::MuonReadoutElement::m_eta {-1}

privateinherited

Identifier field of the station eta.

Definition at line 165 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{-1};

m_firstwire_x

std::array<double, maxnlayers> MuonGM::MdtReadoutElement::m_firstwire_x {make_array<double, maxnlayers>(-9999.)}

private

Definition at line 229 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{make_array<double, maxnlayers>(-9999.)};

m_firstwire_y

std::array<double, maxnlayers> MuonGM::MdtReadoutElement::m_firstwire_y {make_array<double, maxnlayers>(-9999.)}

private

Definition at line 230 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{make_array<double, maxnlayers>(-9999.)};

m_hasCutouts

bool MuonGM::MuonReadoutElement::m_hasCutouts {false}

protectedinherited

true is there are cutouts in the readdout-element

Definition at line 152 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{false};            

m_haveDeformTransf

std::atomic<bool> MuonGM::MdtReadoutElement::m_haveDeformTransf {false}

mutableprivate

Definition at line 268 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{false};

m_haveTubeBounds

std::atomic<bool> MuonGM::MdtReadoutElement::m_haveTubeBounds {false}

mutableprivate

Definition at line 267 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{false};

m_haveTubeGeo

std::atomic<bool> MuonGM::MdtReadoutElement::m_haveTubeGeo {false}

mutableprivate

Definition at line 266 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{false};

m_haveTubeSurfaces

std::atomic<bool> MuonGM::MdtReadoutElement::m_haveTubeSurfaces {false}

mutableprivate

Flag whether any elements have been inserted into the corresponding vectors.

Used to speed up the clear-cache operations for the case where the vectors are empty.

Definition at line 265 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{false};

m_id

Identifier MuonGM::MuonReadoutElement::m_id {}

privateinherited

extended data-collection identifier

Definition at line 157 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{};                    

m_idhash

IdentifierHash MuonGM::MuonReadoutElement::m_idhash {0}

privateinherited

data-collection hash identifier

Definition at line 158 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{0};            

m_idHelper

const MdtIdHelper& MuonGM::MdtReadoutElement::m_idHelper {idHelperSvc()->mdtIdHelper()}

private

Definition at line 211 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{idHelperSvc()->mdtIdHelper()};

m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> MuonGM::MuonReadoutElement::m_idHelperSvc {"Muon::MuonIdHelperSvc/MuonIdHelperSvc", "MuonDetectorManager"}

privateinherited

Definition at line 155 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

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

m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

m_inBarrel

bool MuonGM::MdtReadoutElement::m_inBarrel {false}

private

Definition at line 227 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{false};

m_innerRadius

double MuonGM::MdtReadoutElement::m_innerRadius {-9999.}

private

Definition at line 231 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{-9999.};

m_LongRsize

double MuonGM::MuonReadoutElement::m_LongRsize {-9999.}

protectedinherited

Definition at line 141 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{-9999.};

m_LongSsize

double MuonGM::MuonReadoutElement::m_LongSsize {-9999.}

protectedinherited

Definition at line 140 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{-9999.};

m_LongZsize

double MuonGM::MuonReadoutElement::m_LongZsize {-9999.}

protectedinherited

size in the specified direction

Definition at line 142 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{-9999.};

m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

{ MSG::NIL };

m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

m_multilayer

int MuonGM::MdtReadoutElement::m_multilayer {0}

private

Definition at line 216 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{0};

m_muon_mgr

MuonDetectorManager* MuonGM::MuonReadoutElement::m_muon_mgr {nullptr}

privateinherited

Definition at line 170 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{nullptr};

m_nlayers

int MuonGM::MdtReadoutElement::m_nlayers {-1}

private

Definition at line 217 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{-1};

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_nMDTinStation

unsigned int MuonGM::MdtReadoutElement::m_nMDTinStation {0}

private

Definition at line 215 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{0};

m_nsteps

int MuonGM::MdtReadoutElement::m_nsteps {-1}

private

Definition at line 221 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{-1};

m_ntubesinastep

int MuonGM::MdtReadoutElement::m_ntubesinastep {-1}

private

Definition at line 222 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{-1};

m_ntubesperlayer

int MuonGM::MdtReadoutElement::m_ntubesperlayer {-1}

private

Definition at line 220 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{-1};

m_parentMuonStation

const MuonStation* MuonGM::MuonReadoutElement::m_parentMuonStation {nullptr}

privateinherited

Definition at line 169 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{nullptr};

m_phi

int MuonGM::MuonReadoutElement::m_phi {-1}

privateinherited

Identifier field of the station phi.

Definition at line 167 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{-1};

m_Rsize

double MuonGM::MuonReadoutElement::m_Rsize {-9999.}

protectedinherited

Definition at line 138 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{-9999.};

m_Ssize

double MuonGM::MuonReadoutElement::m_Ssize {-9999.}

protectedinherited

Definition at line 137 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{-9999.};

m_stationS

double MuonGM::MuonReadoutElement::m_stationS {0.}

privateinherited

Definition at line 161 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{0.};

m_statname

std::string MuonGM::MuonReadoutElement::m_statname {"XXX0"}

protectedinherited

examples are BMS5, CSS1, EML1

Definition at line 148 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{"XXX0"};     

m_stIdx

int MuonGM::MuonReadoutElement::m_stIdx {-1}

privateinherited

Identifier field of the station index.

Definition at line 163 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{-1};

m_stIdx_BIS

const int MuonGM::MdtReadoutElement::m_stIdx_BIS {m_idHelper.stationNameIndex("BIS")}

private

Definition at line 212 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{m_idHelper.stationNameIndex("BIS")};

m_stIdx_BOL

const int MuonGM::MdtReadoutElement::m_stIdx_BOL {m_idHelper.stationNameIndex("BOL")}

private

Definition at line 213 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{m_idHelper.stationNameIndex("BOL")};

m_techname

std::string MuonGM::MuonReadoutElement::m_techname {"TTT0"}

protectedinherited

MDT or RPC or TGC or CSC plus a two digits subtype; example RPC17.

Definition at line 146 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{"TTT0"};

m_tubeBounds

std::vector<CxxUtils::CachedUniquePtr<Trk::CylinderBounds> > MuonGM::MdtReadoutElement::m_tubeBounds {}

private

Definition at line 259 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{};       // one per step in tube-length

m_tubeGeo

std::vector<CxxUtils::CachedUniquePtr<GeoInfo> > MuonGM::MdtReadoutElement::m_tubeGeo {}

private

Definition at line 252 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{};                      // one per tube

m_tubelayerpitch

double MuonGM::MdtReadoutElement::m_tubelayerpitch {-9999.}

private

Definition at line 219 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{-9999.};

m_tubelength

std::array<double, maxnsteps> MuonGM::MdtReadoutElement::m_tubelength {make_array<double, maxnsteps>(-9999.)}

private

Definition at line 228 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{make_array<double, maxnsteps>(-9999.)};

m_tubelenStepSize

double MuonGM::MdtReadoutElement::m_tubelenStepSize {-9999.}

private

Definition at line 223 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{-9999.};

m_tubepitch

double MuonGM::MdtReadoutElement::m_tubepitch {-9999.}

private

Definition at line 218 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{-9999.};

m_tubeSurfaces

std::vector<CxxUtils::CachedUniquePtr<Trk::StraightLineSurface> > MuonGM::MdtReadoutElement::m_tubeSurfaces {}

private

Definition at line 258 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{};  // one per tube

m_tubeWallThickness

double MuonGM::MdtReadoutElement::m_tubeWallThickness {-9999.}

private

Definition at line 232 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{-9999.};

m_type

Trk::DetectorElemType MuonGM::MuonReadoutElement::m_type {Trk::DetectorElemType::SolidState}

privateinherited

Definition at line 156 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{Trk::DetectorElemType::SolidState};

m_zsignRO_tubeFrame

CxxUtils::CachedValue<int> MuonGM::MdtReadoutElement::m_zsignRO_tubeFrame {}

private

Definition at line 233 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h.

{};  // comes from AMDB CRO location in the station

m_Zsize

double MuonGM::MuonReadoutElement::m_Zsize {-9999.}

protectedinherited

Definition at line 139 of file MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h.

{-9999.};

---

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

• MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h
• MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx