InDetDD::SiDetectorElement Class Reference

Class to hold geometrical description of a silicon detector element. More...

#include <SiDetectorElement.h>

Inheritance diagram for InDetDD::SiDetectorElement:

Collaboration diagram for InDetDD::SiDetectorElement:

Public Member Functions
	SiDetectorElement (const Identifier &id, const SiDetectorDesign *design, const GeoVFullPhysVol *geophysvol, const SiCommonItems *commonItems, const GeoAlignmentStore *geoAlignStore=nullptr)
	Constructor with parameters.
virtual	~SiDetectorElement ()
	Destructor.
	SiDetectorElement ()=delete
	Don't allow no-argument constructor.
	SiDetectorElement (const SiDetectorElement &)=delete
	Don't allow copy constructor.
SiDetectorElement &	operator= (const SiDetectorElement &)=delete
	Don't allow assignment operator.
	SiDetectorElement (SiDetectorElement &&)=delete
	Don't allow move constructor.
SiDetectorElement &	operator= (SiDetectorElement &&)=delete
	Don't allow move assignment operator.
std::string	detectorTypeString () const
	Returns a string of the Detector element type.
Navigation setters
Non-const methods. These are inline methods and set neighbours.
void	setNextInEta (const SiDetectorElement *element)
void	setPrevInEta (const SiDetectorElement *element)
void	setNextInPhi (const SiDetectorElement *element)
void	setPrevInPhi (const SiDetectorElement *element)
void	setOtherSide (const SiDetectorElement *element)
	For SCT only.
Navigation
Inline methods to access neighbours.
const SiDetectorElement *	nextInEta () const
const SiDetectorElement *	prevInEta () const
const SiDetectorElement *	nextInPhi () const
const SiDetectorElement *	prevInPhi () const
const SiDetectorElement *	otherSide () const
	Useful for SCT only.
Identification
Methods to identify the element and identifier manipulation.
bool	isPixel () const
bool	isSCT () const
bool	isPLR () const
bool	isDBM () const
bool	isBarrel () const
bool	isEndcap () const
bool	isBlayer () const
bool	isInnermostPixelLayer () const
bool	isNextToInnermostPixelLayer () const
virtual Identifier	identifierFromCellId (const SiCellId &cellId) const override final
	Identifier <-> SiCellId (ie strip number or pixel eta_index,phi_index) Identifier from SiCellId (ie strip number or pixel eta_index,phi_index)
virtual SiCellId	cellIdFromIdentifier (const Identifier &identifier) const override final
	SiCellId from Identifier.
Surface
const std::vector< const Trk::Surface * > &	surfaces () const
	Returns the full list of surfaces associated to this detector element.
Module Frame
Methods to help work with the module frame. This is mainly of of use in the SCT where the module frame can in general be different from the element frame. It is actully defined as the frame of one of the sides (currently the axial side), but using these methods one does not need to make any assumptions about what the actual frame is. In the following methods the local element frame is the local reconstruction frame of this element.
const Amg::Transform3D &	moduleTransform () const
	Module to global frame transform.
Amg::Transform3D	defModuleTransform () const
	Default module to global frame transform, ie with no misalignment.
Amg::Transform3D	localToModuleFrame (const Amg::Transform3D &localTransform) const
	Take a transform of the local element frame and return its equivalent in the module frame.
Amg::Transform3D	localToModuleTransform () const
	Transformation from local element to module frame.
bool	isModuleFrame () const
	Check if the element and module frame are the same.
Element Extent
Methods to get extent of element in r,phi and z.
virtual double	get_rz () const override final
Angle
double	sinTilt () const
	Compute sin(tilt angle) at a given position: at center.
double	sinTilt (const Amg::Vector2D &localPos) const
	at given local position
double	sinTilt (const Amg::Vector3D &globalPosition) const
	at given global position
double	sinStereo () const
	Compute sin(stereo angle) at a given position: at center.
double	sinStereo (const Amg::Vector2D &localPos) const
	at given local position
double	sinStereo (const Amg::Vector3D &globalPosition) const
	at given global position
double	sinStereoLocal (const Amg::Vector2D &localPos) const
	Angle of strip in local frame with respect to the etaAxis.
double	sinStereoLocal (const Amg::Vector3D &globalPos) const
	See previous method.
bool	isStereo () const
	Check if it is the stereo side (useful for SCT)
Design methods
virtual const SiDetectorDesign &	design () const override final
	access to the local description (inline):
double	phiPitch () const
	Pitch (inline methods)
double	phiPitch (const Amg::Vector2D &localPosition) const
	Useful for SCT Forward.
InDetDD::CarrierType	carrierType () const
	carrier type for readout.
bool	swapPhiReadoutDirection () const
	Determine if readout direction between online and offline needs swapping.
bool	swapEtaReadoutDirection () const
	For eta_index (only relevant for pixel) (inline)
Intersection Tests
bool	nearBondGap (const Amg::Vector2D &localPosition, double etaTol) const
	Test if near bond gap within tolerances.
bool	nearBondGap (const Amg::Vector3D &globalPosition, double etaTol) const
SiCellId	gangedCell (const SiCellId &cellId) const
	If cell is ganged return the id of the other cell which shares the readout for this cell, otherwise return an invalid id.
Miscellaneous
std::pair< Amg::Vector3D, Amg::Vector3D >	endsOfStrip (const Amg::Vector2D &position) const
	Special method for SCT to retrieve the two ends of a "strip" Returned coordinates are in global frame.
virtual Trk::DetectorElemType	detectorType () const override final
	TrkDetElementBase interface detectorType.
Cache handling
Methods to handle invalidating and updating caches. The cached values include values that are affected by alignment
void	invalidate ()
	Signal that cached values are no longer valid.
Common items
Inline method
const SiCommonItems *	getCommonItems () const
Identification
Methods to identify the element and identifier manipulation.
virtual Identifier	identify () const override final
	identifier of this detector element (inline)
virtual IdentifierHash	identifyHash () const override final
	identifier hash (inline)
const AtlasDetectorID *	getIdHelper () const
	Returns the id helper (inline)
Surface
Trk::Surface &	surface ()
	Element Surface.
virtual const Trk::Surface &	surface () const override final
	Return surface associated with this detector element.
virtual const Trk::Surface &	surface (const Identifier &) const override final
	TrkDetElementBase interface (inline)
Transformation
const GeoTrf::Transform3D &	transformHit () const
	Local (simulation/hit frame) to global transform.
const HepGeom::Transform3D &	transformCLHEP () const
	Local (reconstruction frame) to global transform.
virtual const Amg::Transform3D &	transform () const override final
	Return local to global transform.
virtual const Amg::Transform3D &	transform (const Identifier &) const override final
	TrkDetElementBase interface (inline)
const HepGeom::Transform3D	defTransformCLHEP () const
	Default Local (reconstruction frame) to global transform ie with no misalignment.
const Amg::Transform3D	defTransform () const
const HepGeom::Transform3D	recoToHitTransform () const
	Transform to go from local reconstruction frame to local hit frame.
Orientation
double	hitDepthDirection () const
	Directions of hit depth,phi,eta axes relative to reconstruction local position axes (LocalPosition).
double	hitPhiDirection () const
	See previous method.
double	hitEtaDirection () const
	See previous method.
const HepGeom::Vector3D< double > &	phiAxisCLHEP () const
	To determine if readout direction between online and offline needs swapping, see methods swapPhiReadoutDirection() and swapEtaReadoutDirection() below in "Readout Cell id" section.
const Amg::Vector3D &	phiAxis () const
const HepGeom::Vector3D< double > &	etaAxisCLHEP () const
	Get reconstruction local eta axes in global frame.
const Amg::Vector3D &	etaAxis () const
virtual const Amg::Vector3D &	normal () const override final
	Get reconstruction local normal axes in global frame.
virtual const Amg::Vector3D &	normal (const Identifier &) const override final
	TrkDetElementBase interface (inline)
Position
virtual const Amg::Vector3D &	center () const override final
	Center in global coordinates.
virtual const Amg::Vector3D &	center (const Identifier &) const override final
	TrkDetElementBase interface (inline)
HepGeom::Point3D< double >	globalPositionHit (const HepGeom::Point3D< double > &simulationLocalPos) const
	transform a hit local position into a global position (inline):
Amg::Vector3D	globalPositionHit (const Amg::Vector3D &simulationLocalPos) const
HepGeom::Point3D< double >	globalPosition (const HepGeom::Point3D< double > &localPos) const
	transform a reconstruction local position into a global position (inline):
Amg::Vector3D	globalPosition (const Amg::Vector3D &localPos) const
Amg::Vector3D	globalPosition (const Amg::Vector2D &localPos) const
	as in previous method but for 2D local position (inline)
Amg::Vector2D	hitLocalToLocal (double xEta, double xPhi) const
	Simulation/Hit local frame to reconstruction local frame.
HepGeom::Point3D< double >	hitLocalToLocal3D (const HepGeom::Point3D< double > &hitPosition) const
	Same as previuos method but 3D.
Amg::Vector2D	localPosition (const HepGeom::Point3D< double > &globalPosition) const
	transform a global position into a 2D local position (reconstruction frame) (inline)
Amg::Vector2D	localPosition (const Amg::Vector3D &globalPosition) const
Element Extent
Inline methods to get extent of element in r,phi and z.
double	rMin () const
double	rMax () const
double	zMin () const
double	zMax () const
double	phiMin () const
double	phiMax () const
void	getEtaPhiRegion (double deltaZ, double &etaMin, double &etaMax, double &phiMin, double &phiMax, double &rz) const
	Method for building up region of interest table.
Design methods
virtual const Trk::SurfaceBounds &	bounds () const override final
	Return the boundaries of the element.
virtual const Trk::SurfaceBounds &	bounds (const Identifier &) const override final
	TrkDetElementBase interface (inline)
double	width () const
	Methods from design (inline)
double	minWidth () const
	Min width.
double	maxWidth () const
	Max width.
double	length () const
	Length in eta direction (z - barrel, r - endcap)
double	thickness () const
double	etaPitch () const
	Pitch (inline methods)
Intersection Tests
SiIntersect	inDetector (const Amg::Vector2D &localPosition, double phiTol, double etaTol) const
	Test that it is in the active region.
SiIntersect	inDetector (const Amg::Vector3D &globalPosition, double phiTol, double etaTol) const
Readout cell id
Cell id's are the strip number in SCT and phi_index,eta_index in the pixel as defined in the offline identifier. Their direction runs in the distPhi, distEta direction in the Reconstruction local frame. For methods taking a SiCellId (basically phi,eta index for pixel or strip for SCT) you can do the following fro example: For pixel localPositionOfCell(SiCellId(phi_index,eta_index)); For SCT localPositionOfCell(SiCellId(strip));
Identifier	identifierOfPosition (const Amg::Vector2D &localPos) const
	Full identifier of the cell for a given position: assumes a raw local position (no Lorentz shift)
SiCellId	cellIdOfPosition (const Amg::Vector2D &localPos) const
	As in previous method but returns SiCellId.
Amg::Vector2D	rawLocalPositionOfCell (const SiCellId &cellId) const
	Returns position (center) of cell.
Amg::Vector2D	rawLocalPositionOfCell (const Identifier &id) const
	As above.
int	numberOfConnectedCells (const SiCellId cellId) const
	Test if readout cell has more than one diode associated with it.
SiCellId	connectedCell (const SiCellId cellId, int number) const
	Get the cell ids sharing the readout for this cell.
virtual void	updateCache () const
	Recalculate cached values.
bool	depthDirection () const
	Return information on orientation.
bool	etaDirection () const
bool	phiDirection () const
double	depthAngle () const
double	etaAngle () const
double	phiAngle () const

Protected Attributes
Cache vector of surfaces.
CxxUtils::CachedValue< std::vector< const Trk::Surface * > >	m_surfaces
Variables set by constructor
const SiDetectorDesign *	m_siDesign
Variables set by commonConstructor
bool	m_isPixel {false}
bool	m_isSCT {false}
bool	m_isPLR {false}
bool	m_isDBM {false}
bool	m_isBarrel {false}
Variables set by individual setter methods
const SiDetectorElement *	m_nextInEta {nullptr}
const SiDetectorElement *	m_prevInEta {nullptr}
const SiDetectorElement *	m_nextInPhi {nullptr}
const SiDetectorElement *	m_prevInPhi {nullptr}
const SiDetectorElement *	m_otherSide {nullptr}
Variables set by constructor
Protected data:
Identifier	m_id {}
	identifier of this detector element
const DetectorDesign *	m_design {nullptr}
	local description of this detector element
const SiCommonItems *	m_commonItems {nullptr}
std::unique_ptr< Trk::Surface >	m_surface
const GeoAlignmentStore *	m_geoAlignStore {}
DetectorDesign::Axis	m_hitEta
	Axes.
DetectorDesign::Axis	m_hitPhi
DetectorDesign::Axis	m_hitDepth
Variables set by commonConstructor
IdentifierHash	m_idHash {}
	hash id of this detector element
Variables set by updateCache and not invalidated.
Directions of axes. These are true if the hit/simulation and reconstruction local frames are in the same direction and false if they are opposite.
CxxUtils::CachedValue< AxisDir >	m_axisDir
Variables set by updateCache
CxxUtils::CachedValue< CachedVals >	m_cache

Private Member Functions
Private Methods
Do not use locks.
void	commonConstructor ()
	Common code for constructors.
bool	determineStereo () const
	Find isStereo.
double	sinStereoImpl () const
	Private implementation method with no lock at center.
double	sinStereoImpl (const Amg::Vector3D &globalPosition) const
	Private implementation method with no lock at given global position.
Private Methods
Do not use locks.
void	getExtent (CachedVals &cache) const
	Calculate extent in r,z and phi.
void	getCorners (HepGeom::Point3D< double > *corners) const
	Return the four corners of an element in local coordinates.
void	getEtaPhiPoint (const HepGeom::Point3D< double > &point, double deltaZ, double &etaMin, double &etaMax, double &phi) const
	Get eta and phi coresponding to a point in local coordinates.

Cache result of determineStereo().
CxxUtils::CachedValue< bool >	m_isStereo
MsgStream &	msg (MSG::Level lvl) const
	Declaring the Message method for further use (inline)
bool	msgLvl (MSG::Level lvl) const
	Declaring the Method providing Verbosity Level (inline)

Detailed Description

Class to hold geometrical description of a silicon detector element.

A detector element is a module in the pixels and one side of a module in the SCT. This class is shared between the Pixel and SCT detector since there is a lot of commonality. It inherits frm ReadoutGeometryBase/SolidStateDetectorElement, which is also used for HGTD

Coordinate Frames.

The following coordinate frames are used in these elements.

• Global frame:
  Currently global frame in G4/GeoModel. Probably eventually will be global frame most suitable for reconstruction (eg solenoid axis).
• Local hit frame:
  Local frame for hits. It is the same as local frame in G4 and GeoModel. I also refer to this as the local simulation frame. By convention elements are orientated such that:
• hitDepth = local x
• hitPhi = local y
• hitEta = local z

Directions of these correspond to the physical wafer. Consequently hitDepth and hitPhi axes go in different directions depending on the orientation of the module. The readout side is determined from design()->readoutSide().

• Local reconstruction frame:
  
• distPhi = local x
• distEta = local y
• distDepth = local z

The directions of the axes are defined as

• distPhi in direction of increasing phi
• distEta in direction of increasing z in barrel and increasing r in endcap.
• distDepth (normal) choosen to give right-handed coordinate. => away from intersection point for barrel, decreasing z for endcap

Overview of Methods

Methods are grouped into the the following categories

• Identification
• Navigation
• Module Frame
• Element Extent
• Design methods
• Intersection Tests
• Lorentz Correction
• Readout cell id
• Miscellaneous

Author

Grant Gorfine

• modified & maintained: Nick Styles, Andreas Salzburger
• modified Nigel Hessey: get directions from the design instead of hard-wiring them

Some notes on Thread safety for AthenaMT

The private methods of this class do not have locks.

The method updateCache of is of particular interest as it set all cache values.

The const methods call updateCache() when they need to perform lazy initialization

if (!m_cache.isValid()) updateCache();

So as concurrent const operations are valid

• and do not race with each other.

The non-const methods can set the state of the cache or the cache itself (invalidate methods etc)

Note: Synchronisation of creating SiDetElements for different events and accessing for each events can be done via write/read handles or similar EventContext aware framework machinery.

Definition at line 109 of file SiDetectorElement.h.

Constructor & Destructor Documentation

SiDetectorElement() [1/4]

InDetDD::SiDetectorElement::SiDetectorElement	(	const Identifier &	id,
		const SiDetectorDesign *	design,
		const GeoVFullPhysVol *	geophysvol,
		const SiCommonItems *	commonItems,
		const GeoAlignmentStore *	geoAlignStore = nullptr )

Constructor with parameters.

Definition at line 38 of file SiDetectorElement.cxx.

                                                                               :
    SolidStateDetectorElementBase(id,design,geophysvol,commonItems,geoAlignStore),
    m_siDesign(design)
  {
    commonConstructor();
  }

~SiDetectorElement()

InDetDD::SiDetectorElement::~SiDetectorElement ( )

virtualdefault

Destructor.

SiDetectorElement() [2/4]

InDetDD::SiDetectorElement::SiDetectorElement ( )

delete

Don't allow no-argument constructor.

SiDetectorElement() [3/4]

InDetDD::SiDetectorElement::SiDetectorElement ( const SiDetectorElement & )

delete

Don't allow copy constructor.

SiDetectorElement() [4/4]

InDetDD::SiDetectorElement::SiDetectorElement ( SiDetectorElement && )

delete

Don't allow move constructor.

Member Function Documentation

bounds() [1/2]

const Trk::SurfaceBounds & InDetDD::SolidStateDetectorElementBase::bounds ( ) const

finaloverridevirtualinherited

Return the boundaries of the element.

Implements Trk::TrkDetElementBase.

Definition at line 198 of file SolidStateDetectorElementBase.cxx.

  {
    return m_design->bounds();
  }

bounds() [2/2]

virtual const Trk::SurfaceBounds & InDetDD::SolidStateDetectorElementBase::bounds ( const Identifier & ) const

finaloverridevirtualinherited

TrkDetElementBase interface (inline)

Implements Trk::TrkDetElementBase.

carrierType()

InDetDD::CarrierType InDetDD::SiDetectorElement::carrierType

(

)

const

carrier type for readout.

ie holes for SCT and electrons for pixels. inline

cellIdFromIdentifier()

SiCellId InDetDD::SiDetectorElement::cellIdFromIdentifier ( const Identifier & identifier ) const

finaloverridevirtual

SiCellId from Identifier.

Implements InDetDD::SolidStateDetectorElementBase.

Definition at line 120 of file SiDetectorElement.cxx.

  {
    SiCellId cellId; // Initialized in invalid state.
 
    // If something fails it returns the cellId in an invalid state.
 
    if (identifier.is_valid()) {
      const auto *const pAtlasHelper = getIdHelper();
      if (isPixel() and pAtlasHelper->helper() == AtlasDetectorID::HelperType::Pixel) {
        const auto *const pixelIdHelper = static_cast<const PixelID*>(pAtlasHelper);
        cellId = SiCellId(pixelIdHelper->phi_index(identifier), pixelIdHelper->eta_index(identifier));
      } else if (isSCT() and pAtlasHelper->helper() == AtlasDetectorID::HelperType::SCT) {
        const SCT_ID* sctIdHelper = static_cast<const SCT_ID*>(pAtlasHelper);
        //This adds some extra code for supporting rows,
        //but this method is only used in validation-type code
        //So should not add an overhead in normal running
        //(although we perhaps still try to avoid this...)
        int strip = sctIdHelper->strip(identifier);
        int row = sctIdHelper->row(identifier);
        if(row>0){
          const auto &sctDesign = *static_cast<const SiDetectorDesign *>(m_siDesign);
          int strip1D = sctDesign.strip1Dim(strip, row);
          cellId = SiCellId(strip1D);
        }else {
          cellId =  SiCellId(strip);
        }
 
      } else if (isPLR() and pAtlasHelper->helper() == AtlasDetectorID::HelperType::PLR) {
        const PLR_ID* plrIdHelper = static_cast<const PLR_ID*>(pAtlasHelper);
        cellId = SiCellId(plrIdHelper->phi_index(identifier), plrIdHelper->eta_index(identifier));
      }
    }
 
    return cellId;
  }

cellIdOfPosition()

SiCellId InDetDD::SolidStateDetectorElementBase::cellIdOfPosition ( const Amg::Vector2D & localPos ) const

inherited

As in previous method but returns SiCellId.

Definition at line 224 of file SolidStateDetectorElementBase.cxx.

  {
    return m_design->cellIdOfPosition(localPosition);
  }

center() [1/2]

virtual const Amg::Vector3D & InDetDD::SolidStateDetectorElementBase::center ( ) const

finaloverridevirtualinherited

Center in global coordinates.

Implements Trk::TrkDetElementBase.

center() [2/2]

virtual const Amg::Vector3D & InDetDD::SolidStateDetectorElementBase::center ( const Identifier & ) const

finaloverridevirtualinherited

TrkDetElementBase interface (inline)

Implements Trk::TrkDetElementBase.

commonConstructor()

void InDetDD::SiDetectorElement::commonConstructor ( )

private

Common code for constructors.

Definition at line 349 of file SiDetectorElement.cxx.

  {
    if (!m_id.is_valid()) throw std::runtime_error("SiDetectorElement: Invalid identifier");
 
    // Set booleans for wether we are pixel/sct barrel/endcap
    m_isPixel = getIdHelper()->is_pixel(m_id);
    m_isSCT = getIdHelper()->is_sct(m_id);
    // we use is_lumi here instead of is_plr because is_plr is currently only setup
    // for ExpandedIdentifiers and not Identifiers, which is what is needed here.
    m_isPLR = getIdHelper()->is_lumi(m_id);
    if (!m_isPixel && !m_isSCT && !m_isPLR) {
      ATH_MSG_WARNING("Element id is not for pixel, SCT, or PLR");
    }
 
    // Set m_idHash. Also set m_isBarrel.
    if (isPixel() and getIdHelper()->helper() == AtlasDetectorID::HelperType::Pixel) {
      const PixelID* pixelId = static_cast<const PixelID*>(getIdHelper());
      if (pixelId) {
        m_isBarrel = pixelId->is_barrel(m_id);
        m_idHash = pixelId->wafer_hash(m_id);
 
        if (pixelId->is_dbm(m_id)) {
          m_isBarrel = false;
          m_isDBM = true;
        }
      }
    } else if (isSCT() and getIdHelper()->helper() == AtlasDetectorID::HelperType::SCT) {
      const SCT_ID* sctId = static_cast<const SCT_ID*>(getIdHelper());
      if (sctId) {
        m_isBarrel = sctId->is_barrel(m_id);
        m_idHash = sctId->wafer_hash(m_id);
      }
    } else if (isPLR() and getIdHelper()->helper() == AtlasDetectorID::HelperType::PLR) {
      const PLR_ID* plrId = static_cast<const PLR_ID*>(getIdHelper());
      if (plrId) {
        m_isBarrel = plrId->is_barrel(m_id);
        m_idHash = plrId->wafer_hash(m_id);
      }
    }
 
    if (!m_idHash.is_valid()) throw std::runtime_error("SiDetectorElement: Unable to set IdentifierHash");
 
    // Set surface
    m_surface = std::make_unique<Trk::PlaneSurface>(*this);
 
  }

connectedCell()

SiCellId InDetDD::SolidStateDetectorElementBase::connectedCell ( const SiCellId cellId, int number ) const

inherited

Get the cell ids sharing the readout for this cell.

number = 0 will return the primary readout cell id.

Definition at line 250 of file SolidStateDetectorElementBase.cxx.

  {
    SiReadoutCellId readoutId = m_design->readoutIdOfCell(cellId);
    return m_design->connectedCell(readoutId, number);
  }

defModuleTransform()

Amg::Transform3D InDetDD::SiDetectorElement::defModuleTransform

(

)

const

Default module to global frame transform, ie with no misalignment.

The module frame is defined to be the local reconstruction frame of the axial layer in the SCT. For the pixel it is the same as the element.

Definition at line 179 of file SiDetectorElement.cxx.

  {
    return (isModuleFrame()) ? defTransform() : m_otherSide->defTransform();
  }

defTransform()

const Amg::Transform3D InDetDD::SolidStateDetectorElementBase::defTransform ( ) const

inherited

Definition at line 60 of file SolidStateDetectorElementBase.cxx.

  {
    HepGeom::Transform3D tmpTransform = defTransformCLHEP();
    return Amg::CLHEPTransformToEigen(tmpTransform);
  }

defTransformCLHEP()

const HepGeom::Transform3D InDetDD::SolidStateDetectorElementBase::defTransformCLHEP ( ) const

inherited

Default Local (reconstruction frame) to global transform ie with no misalignment.

Definition at line 50 of file SolidStateDetectorElementBase.cxx.

  {
    if (m_geoAlignStore) {
      const GeoTrf::Transform3D* ptrXf = m_geoAlignStore->getDefAbsPosition(getMaterialGeom());
      if (ptrXf) return Amg::EigenTransformToCLHEP(*ptrXf) * recoToHitTransform();
    }
    return Amg::EigenTransformToCLHEP(getMaterialGeom()->getDefAbsoluteTransform()) * recoToHitTransform();
  }

depthAngle()

double InDetDD::SolidStateDetectorElementBase::depthAngle ( ) const

inherited

depthDirection()

bool InDetDD::SolidStateDetectorElementBase::depthDirection ( ) const

inherited

Return information on orientation.

design()

virtual const SiDetectorDesign & InDetDD::SiDetectorElement::design ( ) const

finaloverridevirtual

access to the local description (inline):

Reimplemented from InDetDD::SolidStateDetectorElementBase.

detectorType()

virtual Trk::DetectorElemType InDetDD::SiDetectorElement::detectorType ( ) const

finaloverridevirtual

TrkDetElementBase interface detectorType.

Reimplemented from InDetDD::SolidStateDetectorElementBase.

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";
}

determineStereo()

bool InDetDD::SiDetectorElement::determineStereo ( ) const

private

Find isStereo.

Definition at line 397 of file SiDetectorElement.cxx.

  {
    if (isSCT() and m_otherSide and getIdHelper()->helper() == AtlasDetectorID::HelperType::SCT) {
      double sinStereoThis = std::abs(sinStereoImpl()); // Call the private impl method
      double sinStereoOther = std::abs(m_otherSide->sinStereo());
      if (std::abs(sinStereoThis - sinStereoOther) < 1e-5) {
        // If they happen to be equal then set side0 as axial and side1 as stereo.
        const SCT_ID* sctId = static_cast<const SCT_ID*>(getIdHelper());
        if (sctId) {
          int side = sctId->side(m_id);
         return (side == 1);
        }
      } else {
        // set the stereo side as the one with largest absolute sinStereo.
        return (sinStereoThis > sinStereoOther);
      }
    }
 
    return false;
  }

endsOfStrip()

std::pair< Amg::Vector3D, Amg::Vector3D > InDetDD::SiDetectorElement::endsOfStrip

(

const Amg::Vector2D &

position

)

const

Special method for SCT to retrieve the two ends of a "strip" Returned coordinates are in global frame.

Definition at line 339 of file SiDetectorElement.cxx.

  {
    //again with the casting...
    const std::pair<Amg::Vector2D,Amg::Vector2D> localEnds=
      static_cast<const SiDetectorDesign *>(m_siDesign)->endsOfStrip(position);
    return std::pair<Amg::Vector3D,Amg::Vector3D >(globalPosition(localEnds.first),
                                                   globalPosition(localEnds.second));
  }

etaAngle()

double InDetDD::SolidStateDetectorElementBase::etaAngle ( ) const

inherited

etaAxis()

const Amg::Vector3D & InDetDD::SolidStateDetectorElementBase::etaAxis ( ) const

inherited

Definition at line 88 of file SolidStateDetectorElementBase.cxx.

  {
    if (!m_cache.isValid()) updateCache();
    return m_cache.ptr()->m_etaAxis;
  }

etaAxisCLHEP()

const HepGeom::Vector3D< double > & InDetDD::SolidStateDetectorElementBase::etaAxisCLHEP ( ) const

inherited

Get reconstruction local eta axes in global frame.

In direction of increasing z in the barrel and increasing r in the endcap.

Definition at line 81 of file SolidStateDetectorElementBase.cxx.

  {
    if (!m_cache.isValid()) updateCache();
    return m_cache.ptr()->m_etaAxisCLHEP;
  }

etaDirection()

bool InDetDD::SolidStateDetectorElementBase::etaDirection ( ) const

inherited

etaPitch()

double InDetDD::SolidStateDetectorElementBase::etaPitch ( ) const

inherited

Pitch (inline methods)

NOTE: phiPitch is ambiguous for the Forward SCT where it varies along the strip.
etaPitch is ambiguous for the pixel which has long pixels between FE chips.
For these cases:

phiPitch: For SCT Forward returns pitch at center.
etaPitch: For pixel returns average pitch. (Active_length/number_of_cells)
All return pitch in distance units.

gangedCell()

SiCellId InDetDD::SiDetectorElement::gangedCell

(

const SiCellId &

cellId

)

const

If cell is ganged return the id of the other cell which shares the readout for this cell, otherwise return an invalid id.

This is a more convenient (and slightly faster) alternative than using the above two methods. Only relevant for pixel. For SCT always returns an invalid ID. (inline)

get_rz()

double InDetDD::SiDetectorElement::get_rz ( ) const

finaloverridevirtual

Implements InDetDD::SolidStateDetectorElementBase.

Definition at line 309 of file SiDetectorElement.cxx.

  {
    // Calculate rz = r (barrel) or z (endcap)
    // Use center of sensor ((0,0,0) in local coordinates) for determining this.
    //  HepGeom::Point3D<double> globalCenter = globalPosition(HepGeom::Point3D<double>(0,0,0));
    if (isBarrel()) {
      return center().perp(); // r
    } else {
      return center().z();  // z
    }
  }

getCommonItems()

const SiCommonItems * InDetDD::SolidStateDetectorElementBase::getCommonItems ( ) const

inherited

getCorners()

void InDetDD::SolidStateDetectorElementBase::getCorners ( HepGeom::Point3D< double > * corners ) const

privateinherited

Return the four corners of an element in local coordinates.

Pass it an array of length 4. This function is used by getEtaPhiRegion()

Definition at line 551 of file SolidStateDetectorElementBase.cxx.

  {
    // This makes the assumption that the forward SCT detectors are orientated such that
    // the positive etaAxis corresponds to the top of the detector where the width is largest.
    // This is currently always the case.
    // For the SCT barrel and pixel detectors minWidth and maxWidth are the same and so should
    // work for all orientations.
 
    double tmpMinWidth = minWidth();
    double tmpMaxWidth = maxWidth();
    double tmpLength   = length();
 
    // Lower left
    corners[0][distPhi] = -0.5 * tmpMinWidth;
    corners[0][distEta] = -0.5 * tmpLength;
    corners[0][distDepth] = 0.;
 
    // Lower right
    corners[1][distPhi] =  0.5 * tmpMinWidth;
    corners[1][distEta] = -0.5 * tmpLength;
    corners[1][distDepth] = 0.;
 
    // Upper Right
    corners[2][distPhi] = 0.5 * tmpMaxWidth;
    corners[2][distEta] = 0.5 * tmpLength;
    corners[2][distDepth] = 0.;
 
    // Upper left
    corners[3][distPhi] = -0.5 * tmpMaxWidth;
    corners[3][distEta] =  0.5 * tmpLength;
    corners[3][distDepth] = 0.;
  }

getEtaPhiPoint()

void InDetDD::SolidStateDetectorElementBase::getEtaPhiPoint ( const HepGeom::Point3D< double > & point, double deltaZ, double & etaMin, double & etaMax, double & phi ) const

privateinherited

Get eta and phi coresponding to a point in local coordinates.

Requires as input the vertex spread. Returns etaMin, etaMax, and phi. This function is used by getEtaPhiRegion()

Definition at line 589 of file SolidStateDetectorElementBase.cxx.

  {
    // Get the point in global coordinates.
    HepGeom::Point3D<double> globalPoint = globalPosition(point);
 
    double r = globalPoint.perp();
    double z = globalPoint.z();
 
    double thetaMin = std::atan2(r,(z + deltaZ));
    etaMax = -std::log(tan(0.5 * thetaMin));
    double thetaMax = std::atan2(r,(z - deltaZ));
    etaMin = -std::log(tan(0.5 * thetaMax));
 
    phi = globalPoint.phi();
  }

getEtaPhiRegion()

void InDetDD::SolidStateDetectorElementBase::getEtaPhiRegion ( double deltaZ, double & etaMin, double & etaMax, double & phiMin, double & phiMax, double & rz ) const

inherited

Method for building up region of interest table.

Get eta/phi extent for the element. Returns min/max eta and phi and r (for barrel) or z (for endcap) Takes as input the vertex spread in z (-deltaZ to +deltaZ)

Definition at line 140 of file SolidStateDetectorElementBase.cxx.

  {
    HepGeom::Point3D<double> corners[4];
    getCorners(corners);
 
    bool first = true;
 
    double phiOffset = 0.;
 
    for (auto & corner : corners) {
      double etaMinPoint = 0.;
      double etaMaxPoint = 0.;
      double phiPoint = 0.;
 
      // Get the eta phi value for this corner.
      getEtaPhiPoint(corner, deltaZ, etaMinPoint, etaMaxPoint, phiPoint);
 
      if (first) { // Use the first point to initialize the min/max values.
 
        // Put phi in a range so that we are not near -180/+180 division.
        // Do this by adding an offset if phi > 90 CLHEP::deg or < -90 CLHEP::deg.
        // This offset is later removed.
        if (phiPoint < -0.5 * M_PI) {
          phiOffset = -0.5 * M_PI;
        } else if (phiPoint > 0.5 * M_PI) {
          phiOffset = 0.5 * M_PI;
        }
        phiMin = phiMax = phiPoint - phiOffset;
        etaMin = etaMinPoint;
        etaMax = etaMaxPoint;
      } else {
        phiPoint -= phiOffset;
        // put phi back in -M_PI < phi < +M_PI range
        if (phiPoint < -M_PI) phiPoint += 2. * M_PI;
        if (phiPoint >  M_PI) phiPoint -= 2. * M_PI;
        phiMin = std::min(phiMin, phiPoint);
        phiMax = std::max(phiMax, phiPoint);
        etaMin = std::min(etaMin, etaMinPoint);
        etaMax = std::max(etaMax, etaMaxPoint);
      }
      first = false;
    }
 
    // put phi back in -M_PI < phi < +M_PI range
    phiMin += phiOffset;
    phiMax += phiOffset;
    if (phiMin < -M_PI) phiMin += 2. * M_PI;
    if (phiMin >  M_PI) phiMin -= 2. * M_PI;
    if (phiMax < -M_PI) phiMax += 2. * M_PI;
    if (phiMax >  M_PI) phiMax -= 2. * M_PI;
    //does it make sense having this if the below is pure virtual for this base class?
    //Should this just be postponed to the derived classes?
    rz = get_rz();
 
  }

getExtent()

void InDetDD::SolidStateDetectorElementBase::getExtent ( CachedVals & cache ) const

privateinherited

Calculate extent in r,z and phi.

The values are cached and there are rMin(), rMax etc methods. It is only used from updateCache

Definition at line 483 of file SolidStateDetectorElementBase.cxx.

  {
    Amg::Vector3D sensorCenter = m_design->sensorCenter();
    double radialShift = sensorCenter[0];//in sensor frame, x is radius
 
    HepGeom::Point3D<double> corners[4];
    getCorners(corners);
 
    bool first = true;
 
    double phiOffset = 0.;
 
 
    const HepGeom::Transform3D rShift = HepGeom::TranslateY3D(radialShift);//in local frame, radius is y=distEta
 
    for (auto & corner : corners) {
 
      corner.transform(rShift);
 
      // m_tranform is already there as  part of the cache construction
      // This method seems to be used only as a helper for updateCache
      HepGeom::Point3D<double> globalPoint = cache.m_transformCLHEP * corner;
 
      double rPoint = globalPoint.perp();
      double zPoint = globalPoint.z();
      double phiPoint = globalPoint.phi();
 
      // Use first point to initializa min/max values.
      if (first) {
 
        // Put phi in a range so that we are not near -180/+180 division.
        // Do this by adding an offset if phi > 90 CLHEP::deg or < -90 CLHEP::deg.
        // This offset is later removed.
        if (phiPoint < -0.5 * M_PI) {
          phiOffset = -0.5 * M_PI;
        } else if (phiPoint > 0.5 * M_PI) {
          phiOffset = 0.5 * M_PI;
        }
        cache.m_minPhi = cache.m_maxPhi = phiPoint - phiOffset;
        cache.m_minR = cache.m_maxR = rPoint;
        cache.m_minZ = cache.m_maxZ = zPoint;
 
      } else {
        phiPoint -= phiOffset;
        // put phi back in -M_PI < phi < +M_PI range
        if (phiPoint < -M_PI) phiPoint += 2. * M_PI;
        if (phiPoint > M_PI)  phiPoint -= 2. * M_PI;
        cache.m_minPhi = std::min(cache.m_minPhi, phiPoint);
        cache.m_maxPhi = std::max(cache.m_maxPhi, phiPoint);
        cache.m_minR = std::min(cache.m_minR, rPoint);
        cache.m_maxR = std::max(cache.m_maxR, rPoint);
        cache.m_minZ = std::min(cache.m_minZ, zPoint);
        cache.m_maxZ = std::max(cache.m_maxZ, zPoint);
      }
      first = false;
    }
 
    // put phi back in -M_PI < phi < +M_PI range
    cache.m_minPhi += phiOffset;
    cache.m_maxPhi += phiOffset;
    if (cache.m_minPhi < -M_PI) cache.m_minPhi += 2. * M_PI;
    if (cache.m_minPhi >  M_PI) cache.m_minPhi -= 2. * M_PI;
    if (cache.m_maxPhi < -M_PI) cache.m_maxPhi += 2. * M_PI;
    if (cache.m_maxPhi >  M_PI) cache.m_maxPhi -= 2. * M_PI;
 
  }

getIdHelper()

const AtlasDetectorID * InDetDD::SolidStateDetectorElementBase::getIdHelper ( ) const

inherited

Returns the id helper (inline)

globalPosition() [1/3]

Amg::Vector3D InDetDD::SolidStateDetectorElementBase::globalPosition ( const Amg::Vector2D & localPos ) const

inherited

as in previous method but for 2D local position (inline)

globalPosition() [2/3]

Amg::Vector3D InDetDD::SolidStateDetectorElementBase::globalPosition ( const Amg::Vector3D & localPos ) const

inherited

globalPosition() [3/3]

HepGeom::Point3D< double > InDetDD::SolidStateDetectorElementBase::globalPosition ( const HepGeom::Point3D< double > & localPos ) const

inherited

transform a reconstruction local position into a global position (inline):

globalPositionHit() [1/2]

Amg::Vector3D InDetDD::SolidStateDetectorElementBase::globalPositionHit ( const Amg::Vector3D & simulationLocalPos ) const

inherited

globalPositionHit() [2/2]

HepGeom::Point3D< double > InDetDD::SolidStateDetectorElementBase::globalPositionHit ( const HepGeom::Point3D< double > & simulationLocalPos ) const

inherited

transform a hit local position into a global position (inline):

hitDepthDirection()

double InDetDD::SolidStateDetectorElementBase::hitDepthDirection ( ) const

inherited

Directions of hit depth,phi,eta axes relative to reconstruction local position axes (LocalPosition).

Returns +/-1. inline

hitEtaDirection()

double InDetDD::SolidStateDetectorElementBase::hitEtaDirection ( ) const

inherited

See previous method.

inline

hitLocalToLocal()

Amg::Vector2D InDetDD::SolidStateDetectorElementBase::hitLocalToLocal ( double xEta, double xPhi ) const

inherited

Simulation/Hit local frame to reconstruction local frame.

2D. TODO: Will change order of parameters at some point.

Definition at line 95 of file SolidStateDetectorElementBase.cxx.

  {                                                                               // due to whether module is centred on
    if (!m_axisDir.isValid()) updateCache();                                      // z or y axes
    const AxisDir& dir = *m_axisDir.ptr();
 
    if (!dir.m_etaDirection) xEta = -xEta;
    if (!dir.m_phiDirection) xPhi = -xPhi;
    auto result = Amg::Vector2D(xPhi, xEta);
 
    if (m_design->shape() == InDetDD::PolarAnnulus) { // Do conversion to polar co-ords as well
      double y = result.x(); // Co-ordinate flip from cartesian needs to be temporarily un-done to
      double x = result.y(); // allow atan2 to work in conversion
 
      double r = std::hypot(x,y);
      double phi = std::atan2(y,x);
      result = Amg::Vector2D(phi,r); // now flip again
    }
 
    return result;
  }

hitLocalToLocal3D()

HepGeom::Point3D< double > InDetDD::SolidStateDetectorElementBase::hitLocalToLocal3D ( const HepGeom::Point3D< double > & hitPosition ) const

inherited

Same as previuos method but 3D.

Definition at line 117 of file SolidStateDetectorElementBase.cxx.

  {
    // Equiv to transform().inverse * transformHit() * hitPosition
    if (!m_axisDir.isValid()) updateCache();
    const AxisDir& dir = *m_axisDir.ptr();
 
    double xDepth = hitPosition[m_hitDepth];
    double xPhi = hitPosition[m_hitPhi];
    double xEta = hitPosition[m_hitEta];
    if (!dir.m_depthDirection) xDepth = -xDepth;
    if (!dir.m_phiDirection) xPhi = -xPhi;
    if (!dir.m_etaDirection) xEta = -xEta;
    return {xPhi, xEta, xDepth};
  }

hitPhiDirection()

double InDetDD::SolidStateDetectorElementBase::hitPhiDirection ( ) const

inherited

See previous method.

inline

identifierFromCellId()

Identifier InDetDD::SiDetectorElement::identifierFromCellId ( const SiCellId & cellId ) const

finaloverridevirtual

Identifier <-> SiCellId (ie strip number or pixel eta_index,phi_index) Identifier from SiCellId (ie strip number or pixel eta_index,phi_index)

Implements InDetDD::SolidStateDetectorElementBase.

Definition at line 89 of file SiDetectorElement.cxx.

  {
    Identifier id; // Will be initialized in an invalid state.
 
    // If something fails it returns the id in an invalid state.
 
    if (cellId.isValid()) {
      const auto *const pAtlasHelper = getIdHelper();
      if (isPixel()) {
        if (pAtlasHelper->helper() ==  AtlasDetectorID::HelperType::Pixel) {
          const PixelID* pixelIdHelper = static_cast<const PixelID*>(pAtlasHelper);
          id = pixelIdHelper->pixel_id(m_id, cellId.phiIndex(), cellId.etaIndex());
        }
      } else if (isSCT()) {
 
        if (pAtlasHelper->helper() ==  AtlasDetectorID::HelperType::SCT) {
          const SCT_ID* sctIdHelper = static_cast<const SCT_ID*>(pAtlasHelper);
          id = sctIdHelper->strip_id(m_id, cellId.strip());
        }
      }  else if (isPLR()) {
        if (pAtlasHelper->helper() ==  AtlasDetectorID::HelperType::PLR) {
          const PLR_ID* plrIdHelper = static_cast<const PLR_ID*>(pAtlasHelper);
          id = plrIdHelper->pixel_id(m_id, cellId.phiIndex(), cellId.etaIndex());
        }
      }
    }
 
    return id;
  }

identifierOfPosition()

Identifier InDetDD::SolidStateDetectorElementBase::identifierOfPosition ( const Amg::Vector2D & localPos ) const

inherited

Full identifier of the cell for a given position: assumes a raw local position (no Lorentz shift)

Definition at line 217 of file SolidStateDetectorElementBase.cxx.

  {
    SiCellId cellId = m_design->cellIdOfPosition(localPosition);
    return identifierFromCellId(cellId);
  }

identify()

virtual Identifier InDetDD::SolidStateDetectorElementBase::identify ( ) const

finaloverridevirtualinherited

identifier of this detector element (inline)

Implements Trk::TrkDetElementBase.

identifyHash()

virtual IdentifierHash InDetDD::SolidStateDetectorElementBase::identifyHash ( ) const

finaloverridevirtualinherited

identifier hash (inline)

Implements Trk::TrkDetElementBase.

inDetector() [1/2]

SiIntersect InDetDD::SolidStateDetectorElementBase::inDetector ( const Amg::Vector2D & localPosition, double phiTol, double etaTol ) const

inherited

Test that it is in the active region.

Intersect has 3 states bool SiIntersect::in() const // definitely in
bool SiIntersect::out() const // definitely out
bool SiIntersect::nearBoundary() const // near a boundary within the tolerances
bool SiIntersect::mayIntersect() const // in() OR nearBoundary()

Definition at line 204 of file SolidStateDetectorElementBase.cxx.

  {
    return m_design->inDetector(localPosition, phiTol, etaTol);
  }

inDetector() [2/2]

SiIntersect InDetDD::SolidStateDetectorElementBase::inDetector ( const Amg::Vector3D & globalPosition, double phiTol, double etaTol ) const

inherited

Definition at line 211 of file SolidStateDetectorElementBase.cxx.

  {
    return m_design->inDetector(localPosition(globalPosition), phiTol, etaTol);
  }

invalidate()

void InDetDD::SolidStateDetectorElementBase::invalidate ( )

inherited

Signal that cached values are no longer valid.

Invalidate general cache (inline)

isBarrel()

bool InDetDD::SiDetectorElement::isBarrel

(

)

const

isBlayer()

bool InDetDD::SiDetectorElement::isBlayer

(

)

const

Definition at line 53 of file SiDetectorElement.cxx.

  {
    if (const auto *const pIdHelper= getIdHelper();isPixel() and isBarrel() and
      pIdHelper->helper() ==  AtlasDetectorID::HelperType::Pixel) {
      const PixelID* p_pixelId = static_cast<const PixelID*>(pIdHelper);
      return (0==p_pixelId->layer_disk(m_id));
    } else {
      return false;
    }
  }

isDBM()

bool InDetDD::SiDetectorElement::isDBM

(

)

const

isEndcap()

bool InDetDD::SiDetectorElement::isEndcap

(

)

const

isInnermostPixelLayer()

bool InDetDD::SiDetectorElement::isInnermostPixelLayer

(

)

const

Definition at line 65 of file SiDetectorElement.cxx.

  {
    if (const auto *const pIdHelper= getIdHelper(); isPixel() and isBarrel()
     and pIdHelper->helper() ==  AtlasDetectorID::HelperType::Pixel) {
      const PixelID* p_pixelId = static_cast<const PixelID*>(pIdHelper);
      return ( 0==p_pixelId->layer_disk(m_id));
    } else {
      return false;
    }
  }

isModuleFrame()

bool InDetDD::SiDetectorElement::isModuleFrame

(

)

const

Check if the element and module frame are the same.

Definition at line 212 of file SiDetectorElement.cxx.

  {
    // The module frame is the axial side.
    // NB isStereo returns false for the pixel and so
    // isModuleFrame is always true for the pixel.
 
    return !isStereo();
  }

isNextToInnermostPixelLayer()

bool InDetDD::SiDetectorElement::isNextToInnermostPixelLayer

(

)

const

Definition at line 77 of file SiDetectorElement.cxx.

  {
    if (const auto *const pIdHelper = getIdHelper(); isPixel() and isBarrel() and
      pIdHelper->helper() ==  AtlasDetectorID::HelperType::Pixel) {
      const PixelID* p_pixelId = static_cast<const PixelID*>(pIdHelper);
      return ( 1==p_pixelId->layer_disk(m_id));
    } else {
      return false;
    }
  }

isPixel()

bool InDetDD::SiDetectorElement::isPixel

(

)

const

isPLR()

bool InDetDD::SiDetectorElement::isPLR

(

)

const

isSCT()

bool InDetDD::SiDetectorElement::isSCT

(

)

const

isStereo()

bool InDetDD::SiDetectorElement::isStereo

(

)

const

Check if it is the stereo side (useful for SCT)

Definition at line 300 of file SiDetectorElement.cxx.

  {
    if (!m_isStereo.isValid()) {
      m_isStereo.set (determineStereo());
    }
 
    return *m_isStereo.ptr();
  }

length()

double InDetDD::SolidStateDetectorElementBase::length ( ) const

inherited

Length in eta direction (z - barrel, r - endcap)

localPosition() [1/2]

Amg::Vector2D InDetDD::SolidStateDetectorElementBase::localPosition ( const Amg::Vector3D & globalPosition ) const

inherited

localPosition() [2/2]

Amg::Vector2D InDetDD::SolidStateDetectorElementBase::localPosition ( const HepGeom::Point3D< double > & globalPosition ) const

inherited

transform a global position into a 2D local position (reconstruction frame) (inline)

localToModuleFrame()

Amg::Transform3D InDetDD::SiDetectorElement::localToModuleFrame

(

const Amg::Transform3D &

localTransform

)

const

Take a transform of the local element frame and return its equivalent in the module frame.

Definition at line 192 of file SiDetectorElement.cxx.

  {
    if (isModuleFrame()) {
      return localTransform;
    } else {
      return m_otherSide->transform().inverse() * transform() * localTransform * transform().inverse() * m_otherSide->transform();
    }
  }

localToModuleTransform()

Amg::Transform3D InDetDD::SiDetectorElement::localToModuleTransform

(

)

const

Transformation from local element to module frame.

This can be used to take a local position in the element frame and transform it to a position in the module frame. If one is already in the module frame it will return the Identity transform.

Definition at line 202 of file SiDetectorElement.cxx.

  {
    if (isModuleFrame()) {
      return {}; // Identity
    } else {
      return m_otherSide->transform().inverse() * transform();
    }
  }

maxWidth()

double InDetDD::SolidStateDetectorElementBase::maxWidth ( ) const

inherited

Max width.

Needed for the SCT endcap.

minWidth()

double InDetDD::SolidStateDetectorElementBase::minWidth ( ) const

inherited

Min width.

Needed for the SCT endcap.

moduleTransform()

const Amg::Transform3D & InDetDD::SiDetectorElement::moduleTransform

(

)

const

Module to global frame transform.

Includes misalignment. The module frame is defined to be the local reconstruction frame of the axial layer in the SCT. For the pixel it is the same as the element.

Definition at line 173 of file SiDetectorElement.cxx.

  {
    return (isModuleFrame()) ? transform() : m_otherSide->transform();
  }

msg()

MsgStream & InDetDD::SiDetectorElement::msg ( MSG::Level lvl ) const

protected

Declaring the Message method for further use (inline)

msgLvl()

bool InDetDD::SiDetectorElement::msgLvl ( MSG::Level lvl ) const

protected

Declaring the Method providing Verbosity Level (inline)

nearBondGap() [1/2]

bool InDetDD::SiDetectorElement::nearBondGap	(	const Amg::Vector2D &	localPosition,
		double	etaTol ) const

Test if near bond gap within tolerances.

Definition at line 323 of file SiDetectorElement.cxx.

  {
    //again, can we avoid casting here?
    return static_cast<const SiDetectorDesign *>(m_siDesign)->nearBondGap(localPosition, etaTol);
  }

nearBondGap() [2/2]

bool InDetDD::SiDetectorElement::nearBondGap	(	const Amg::Vector3D &	globalPosition,
		double	etaTol ) const

Definition at line 330 of file SiDetectorElement.cxx.

  {
    return static_cast<const SiDetectorDesign *>(m_siDesign)->nearBondGap(localPosition(globalPosition), etaTol);
  }

nextInEta()

const SiDetectorElement * InDetDD::SiDetectorElement::nextInEta

(

)

const

nextInPhi()

const SiDetectorElement * InDetDD::SiDetectorElement::nextInPhi

(

)

const

normal() [1/2]

virtual const Amg::Vector3D & InDetDD::SolidStateDetectorElementBase::normal ( ) const

finaloverridevirtualinherited

Get reconstruction local normal axes in global frame.

Choosen to give right-handed coordinate frame.

Implements Trk::TrkDetElementBase.

normal() [2/2]

virtual const Amg::Vector3D & InDetDD::SolidStateDetectorElementBase::normal ( const Identifier & ) const

finaloverridevirtualinherited

TrkDetElementBase interface (inline)

Implements Trk::TrkDetElementBase.

numberOfConnectedCells()

int InDetDD::SolidStateDetectorElementBase::numberOfConnectedCells ( const SiCellId cellId ) const

inherited

Test if readout cell has more than one diode associated with it.

Number of cells sharing the same readout as this cell. ie generally 1 except for ganged pixels which will be 2.

Definition at line 243 of file SolidStateDetectorElementBase.cxx.

  {
    SiReadoutCellId readoutId = m_design->readoutIdOfCell(cellId);
    return m_design->numberOfConnectedCells(readoutId);
  }

operator=() [1/2]

SiDetectorElement & InDetDD::SiDetectorElement::operator= ( const SiDetectorElement & )

delete

Don't allow assignment operator.

operator=() [2/2]

SiDetectorElement & InDetDD::SiDetectorElement::operator= ( SiDetectorElement && )

delete

Don't allow move assignment operator.

otherSide()

const SiDetectorElement * InDetDD::SiDetectorElement::otherSide

(

)

const

Useful for SCT only.

phiAngle()

double InDetDD::SolidStateDetectorElementBase::phiAngle ( ) const

inherited

phiAxis()

const Amg::Vector3D & InDetDD::SolidStateDetectorElementBase::phiAxis ( ) const

inherited

Definition at line 74 of file SolidStateDetectorElementBase.cxx.

  {
    if (!m_cache.isValid()) updateCache();
    return m_cache.ptr()->m_phiAxis;
  }

phiAxisCLHEP()

const HepGeom::Vector3D< double > & InDetDD::SolidStateDetectorElementBase::phiAxisCLHEP ( ) const

inherited

To determine if readout direction between online and offline needs swapping, see methods swapPhiReadoutDirection() and swapEtaReadoutDirection() below in "Readout Cell id" section.

Orientation. Directions.

• phiAxis in same direction as increasing phi and identifier phi_index/strip. NB. This requires some flipping of axes with repsect to the hits.
• etaAxis in direction of increasing z in the barrel and increasing r in the endcap.
• normal choosen to give right-handed coordinate frame (x=normal,y=phiAxis,z=etaAxis) NB. This requires some flipping of axes with repsect to the hits.

Get reconstruction local phi axes in global frame. In same direction as increasing phi and identifier phi_index/strip.

Definition at line 67 of file SolidStateDetectorElementBase.cxx.

  {
    if (!m_cache.isValid()) updateCache();
    return m_cache.ptr()->m_phiAxisCLHEP;
  }

phiDirection()

bool InDetDD::SolidStateDetectorElementBase::phiDirection ( ) const

inherited

phiMax()

double InDetDD::SolidStateDetectorElementBase::phiMax ( ) const

inherited

phiMin()

double InDetDD::SolidStateDetectorElementBase::phiMin ( ) const

inherited

phiPitch() [1/2]

double InDetDD::SiDetectorElement::phiPitch

(

)

const

Pitch (inline methods)

NOTE: phiPitch is ambiguous for the Forward SCT where it varies along the strip.
etaPitch is ambiguous for the pixel which has long pixels between FE chips.
For these cases:

phiPitch: For SCT Forward returns pitch at center.
etaPitch: For pixel returns average pitch. (Active_length/number_of_cells)
All return pitch in distance units.

phiPitch() [2/2]

double InDetDD::SiDetectorElement::phiPitch

(

const Amg::Vector2D &

localPosition

)

const

Useful for SCT Forward.

prevInEta()

const SiDetectorElement * InDetDD::SiDetectorElement::prevInEta

(

)

const

prevInPhi()

const SiDetectorElement * InDetDD::SiDetectorElement::prevInPhi

(

)

const

rawLocalPositionOfCell() [1/2]

Amg::Vector2D InDetDD::SolidStateDetectorElementBase::rawLocalPositionOfCell ( const Identifier & id ) const

inherited

As above.

Definition at line 236 of file SolidStateDetectorElementBase.cxx.

  {
    SiCellId cellId = cellIdFromIdentifier(id);
    return m_design->localPositionOfCell(cellId);
  }

rawLocalPositionOfCell() [2/2]

Amg::Vector2D InDetDD::SolidStateDetectorElementBase::rawLocalPositionOfCell ( const SiCellId & cellId ) const

inherited

Returns position (center) of cell.

These are the raw positions NOT corrected for the Lorentz shift

Definition at line 230 of file SolidStateDetectorElementBase.cxx.

  {
    return m_design->localPositionOfCell(cellId);
  }

recoToHitTransform()

const HepGeom::Transform3D InDetDD::SolidStateDetectorElementBase::recoToHitTransform ( ) const

inherited

Transform to go from local reconstruction frame to local hit frame.

Definition at line 607 of file SolidStateDetectorElementBase.cxx.

  {
    if (!m_axisDir.isValid()) updateCache();
    const AxisDir& dir = *m_axisDir.ptr();
 
    //        = transfromHit * hitLocal
    //        = transformHit * recoToHitTransform * recoLocal
    // recoToHitTransform takes recoLocal to hitLocal
    // x,y,z -> y,z,x
    // equiv to a rotation around Y of 90 deg followed by a rotation around X of 90deg
    //
    // recoToHit is static as it needs to be calculated once only.
    // We use the HepGeom::Transform3D constructor which takes one coordinates system to another where the
    // coordinate system is defined by it center and two axes.
    // distPhi, distEta are the reco local axes and hitPhi and hitEta are the hit local axes.
    // It assume phi, eta, depth makes a right handed system which is the case.
    static const HepGeom::Vector3D<double> localAxes[3] = {
      HepGeom::Vector3D<double>(1., 0., 0.),
      HepGeom::Vector3D<double>(0., 1., 0.),
      HepGeom::Vector3D<double>(0., 0., 1.)
    };
 
    //correct phi and eta as necessary - do not change depth, this will be defined by the transform based on the other two
    int signPhi = dir.m_phiDirection? +1:-1;
    int signEta = dir.m_etaDirection? +1:-1;
 
    const HepGeom::Transform3D recoToHit(HepGeom::Point3D<double>(0., 0., 0.),
                                         signPhi * localAxes[distPhi],
                                         signEta *localAxes[distEta],
                                         HepGeom::Point3D<double>(0., 0., 0.),
                                         localAxes[m_hitPhi],
                                         localAxes[m_hitEta]);
 
    return recoToHit ;
 
  }

rMax()

double InDetDD::SolidStateDetectorElementBase::rMax ( ) const

inherited

rMin()

double InDetDD::SolidStateDetectorElementBase::rMin ( ) const

inherited

setNextInEta()

void InDetDD::SiDetectorElement::setNextInEta

(

const SiDetectorElement *

element

)

setNextInPhi()

void InDetDD::SiDetectorElement::setNextInPhi

(

const SiDetectorElement *

element

)

setOtherSide()

void InDetDD::SiDetectorElement::setOtherSide

(

const SiDetectorElement *

element

)

For SCT only.

setPrevInEta()

void InDetDD::SiDetectorElement::setPrevInEta

(

const SiDetectorElement *

element

)

setPrevInPhi()

void InDetDD::SiDetectorElement::setPrevInPhi

(

const SiDetectorElement *

element

)

sinStereo() [1/3]

double InDetDD::SiDetectorElement::sinStereo

(

)

const

Compute sin(stereo angle) at a given position: at center.

Definition at line 269 of file SiDetectorElement.cxx.

  {
    return sinStereoImpl();
  }

sinStereo() [2/3]

double InDetDD::SiDetectorElement::sinStereo

(

const Amg::Vector2D &

localPos

)

const

at given local position

Definition at line 275 of file SiDetectorElement.cxx.

  {
    Amg::Vector3D point=globalPosition(localPos);
    return sinStereoImpl(point);
  }

sinStereo() [3/3]

double InDetDD::SiDetectorElement::sinStereo

(

const Amg::Vector3D &

globalPosition

)

const

at given global position

Definition at line 282 of file SiDetectorElement.cxx.

  {
    return sinStereoImpl(globalPos);
  }

sinStereoImpl() [1/2]

double InDetDD::SiDetectorElement::sinStereoImpl ( ) const

private

Private implementation method with no lock at center.

Definition at line 419 of file SiDetectorElement.cxx.

  {
    // Stereo is the angle between a refVector and a vector along the strip/pixel in eta direction.
    // I'm not sure how the sign should be defined. I've defined it here
    // with rotation sense respect to normal,
    // where normal is away from IP in barrel and in -ve z direction in endcap
 
    // In Barrel refVector = unit vector along z axis,
    // in endcap refVector = unit vector radial.
    //
    // sinStereo = (refVector cross stripAxis) . normal
    //           = (refVector cross etaAxis) . normal
    //           = refVector . (etaAxis cross normal)
    //           = refVector . phiAxis
    //
    // in Barrel we use
    // sinStereo = refVector . phiAxis
    //           = phiAxis.z()
    //
    // in endcap we use
    // sinStereo = (refVector cross etaAxis) . normal
    //           = -(center cross etaAxis) . zAxis
    //           = (etaAxis cross center). z()
 
    //Since these are barrel, endcap, sensor-type, specific, might be better for these to be calculated in the design()
    //However, not clear how one could do that for the annulus calculation which uses global frame
    double sinStereo = 0.;
    auto designShape = m_siDesign->shape();
    if (isBarrel()) {
      sinStereo = this->phiAxis().z();
    } else { // endcap
      if (designShape == InDetDD::Annulus) { //built-in Stereo angle for Annulus shape sensor
        Amg::Vector3D sensorCenter = m_siDesign->sensorCenter();
        //Below retrieved method will return -sin(m_Stereo), thus sinStereolocal = sin(m_Stereo)
        double sinStereoReco = - (m_siDesign->sinStripAngleReco(sensorCenter[1], sensorCenter[0]));
        double cosStereoReco = sqrt(1-sinStereoReco*sinStereoReco);
        double radialShift = sensorCenter[0];
        //The focus of all strips in the local reco frame
        Amg::Vector2D localfocus(-radialShift*sinStereoReco, radialShift - radialShift*cosStereoReco);
        //The focus of all strips in the global frame
        Amg::Vector3D globalfocus(globalPosition(localfocus));
        //The direction of x-axis of the Strip frame in the global frame
              const Amg::Vector3D& center = this->center();
        Amg::Vector3D globalSFxAxis =(center - globalfocus)/radialShift;
        //Stereo angle is the angle between global radial direction and the x-axis of the Strip frame in the global frame
        sinStereo = (center.y() * globalSFxAxis.x() - center.x() * globalSFxAxis.y()) / center.perp();
      }
      // else if (designShape == InDetDD::PolarAnnulus) {} // Polar specialisation in future
      else { // barrel
        const Amg::Vector3D& etaAxis = this->etaAxis();
        const Amg::Vector3D& center = this->center();
        sinStereo = (center.y() * etaAxis.x() - center.x() * etaAxis.y()) / center.perp();
      }
    }
    return sinStereo;
  }

sinStereoImpl() [2/2]

double InDetDD::SiDetectorElement::sinStereoImpl ( const Amg::Vector3D & globalPosition ) const

private

Private implementation method with no lock at given global position.

Definition at line 477 of file SiDetectorElement.cxx.

  {
    //
    // sinStereo =  (refVector cross stripAxis) . normal
    //
    double sinStereo = 0.;
    if (isBarrel()) {
      if (m_siDesign->shape() != InDetDD::Trapezoid) {
        sinStereo = this->phiAxis().z();
      } else { // trapezoid
        assert (minWidth() != maxWidth());
        double radius = width() * length() / (maxWidth() - minWidth());
        const Amg::Vector3D& etaAxis = this->etaAxis();
        const Amg::Vector3D& center = this->center();
        const Amg::Vector3D& normal = this->normal();
    Amg::Vector3D stripAxis = radius * etaAxis + globalPos - center;
        sinStereo = (stripAxis.x() * normal.y() - stripAxis.y() * normal.x()) / stripAxis.mag();
      }
    } else { // endcap
      if (m_siDesign->shape() != InDetDD::Trapezoid) {
        const Amg::Vector3D& etaAxis = this->etaAxis();
        sinStereo = (globalPos.y() * etaAxis.x() - globalPos.x() * etaAxis.y()) / globalPos.perp();
      } else { // trapezoid
        assert (minWidth() != maxWidth());
        const Amg::Vector3D& etaAxis = this->etaAxis();
        const Amg::Vector3D& center = this->center();
        double radius = width() * length() / (maxWidth() - minWidth());
        // Only need projection in xy plane.
        double stripAxisX  = globalPos.x() - center.x() + etaAxis.x()*radius;
        double stripAxisY  = globalPos.y() - center.y() + etaAxis.y()*radius;
        double norm = 1./(radius*sqrt(stripAxisX*stripAxisX + stripAxisY*stripAxisY));
        sinStereo = norm * (stripAxisX * globalPos.y() - stripAxisY * globalPos.x());
      }
    }
    return sinStereo;
  }

sinStereoLocal() [1/2]

double InDetDD::SiDetectorElement::sinStereoLocal

(

const Amg::Vector2D &

localPos

)

const

Angle of strip in local frame with respect to the etaAxis.

Zero for all elements except trapezoidal detectors (ie SCT forward modules).

Definition at line 288 of file SiDetectorElement.cxx.

  {
    return static_cast<const SiDetectorDesign *>(m_design)->sinStripAngleReco(localPos[0], localPos[1]);
  }

sinStereoLocal() [2/2]

double InDetDD::SiDetectorElement::sinStereoLocal

(

const Amg::Vector3D &

globalPos

)

const

See previous method.

Definition at line 294 of file SiDetectorElement.cxx.

  {
    return sinStereoLocal(localPosition(globalPos));
  }

sinTilt() [1/3]

double InDetDD::SiDetectorElement::sinTilt

(

)

const

Compute sin(tilt angle) at a given position: at center.

Definition at line 223 of file SiDetectorElement.cxx.

  {
    // Tilt is defined as the angle between a refVector and the sensor normal.
    // In barrel refVector = unit vector radial.
    // in endcap it is assumed there is no tilt.
    // sinTilt = (refVector cross normal) . z
 
    // tilt angle is not defined for the endcap
    if (isEndcap()) return 0.;
 
    // Angle between normal and radial vector.
    // HepGeom::Vector3D<double> refVector(m_center.x(), m_center.y(), 0);
    // return (refVector.cross(m_normal)).z()/refVector.mag();
    // or the equivalent
    const Amg::Vector3D& normal = this->normal();
    const Amg::Vector3D& center = this->center();
    return (center.x() * normal.y() - center.y() * normal.x()) / center.perp();
  }

sinTilt() [2/3]

double InDetDD::SiDetectorElement::sinTilt

(

const Amg::Vector2D &

localPos

)

const

at given local position

Definition at line 243 of file SiDetectorElement.cxx.

  {
    // tilt angle is not defined for the endcap
    if (isEndcap()) return 0.;
 
    Amg::Vector3D point = globalPosition(localPos);
    return sinTilt(point);
  }

sinTilt() [3/3]

double InDetDD::SiDetectorElement::sinTilt

(

const Amg::Vector3D &

globalPosition

)

const

at given global position

Definition at line 253 of file SiDetectorElement.cxx.

  {
    // It is assumed that the global position is already in the plane of the element.
 
    // tilt angle is not defined for the endcap
    if (isEndcap()) return 0.;
 
    // Angle between normal and radial vector.
    //HepGeom::Vector3D<double> refVector(globalPos.x(), globalPos.y(), 0);
    //return (refVector.cross(m_normal)).z()/refVector.mag();
    // or the equivalent
    const Amg::Vector3D& normal = this->normal();
    return (globalPos.x() * normal.y() - globalPos.y() * normal.x()) / globalPos.perp();
  }

surface() [1/3]

Trk::Surface & InDetDD::SolidStateDetectorElementBase::surface ( )

inherited

Element Surface.

surface() [2/3]

virtual const Trk::Surface & InDetDD::SolidStateDetectorElementBase::surface ( ) const

finaloverridevirtualinherited

Return surface associated with this detector element.

Implements Trk::TrkDetElementBase.

surface() [3/3]

virtual const Trk::Surface & InDetDD::SolidStateDetectorElementBase::surface ( const Identifier & ) const

finaloverridevirtualinherited

TrkDetElementBase interface (inline)

Implements Trk::TrkDetElementBase.

surfaces()

const std::vector< const Trk::Surface * > & InDetDD::SiDetectorElement::surfaces

(

)

const

Returns the full list of surfaces associated to this detector element.

Definition at line 157 of file SiDetectorElement.cxx.

  {
    if (!m_surfaces.isValid()) {
      std::vector<const Trk::Surface*> s;
      s.push_back(&surface());
      if (otherSide()) {
        s.push_back(&(otherSide()->surface()));
      }
      m_surfaces.set (std::move (s));
    }
 
    // return the surfaces
    return *m_surfaces.ptr();
  }

swapEtaReadoutDirection()

bool InDetDD::SiDetectorElement::swapEtaReadoutDirection

(

)

const

For eta_index (only relevant for pixel) (inline)

swapPhiReadoutDirection()

bool InDetDD::SiDetectorElement::swapPhiReadoutDirection

(

)

const

Determine if readout direction between online and offline needs swapping.

Returns true if online and offline numbers run in opposite directions. For strip/phi_index (inline)

thickness()

double InDetDD::SolidStateDetectorElementBase::thickness ( ) const

inherited

transform() [1/2]

virtual const Amg::Transform3D & InDetDD::SolidStateDetectorElementBase::transform ( ) const

finaloverridevirtualinherited

Return local to global transform.

Implements Trk::TrkDetElementBase.

transform() [2/2]

virtual const Amg::Transform3D & InDetDD::SolidStateDetectorElementBase::transform ( const Identifier & ) const

finaloverridevirtualinherited

TrkDetElementBase interface (inline)

Implements Trk::TrkDetElementBase.

transformCLHEP()

const HepGeom::Transform3D & InDetDD::SolidStateDetectorElementBase::transformCLHEP ( ) const

inherited

Local (reconstruction frame) to global transform.

Definition at line 43 of file SolidStateDetectorElementBase.cxx.

  {
    if (!m_cache.isValid()) updateCache();
    return m_cache.ptr()->m_transformCLHEP;
  }

transformHit()

const GeoTrf::Transform3D & InDetDD::SolidStateDetectorElementBase::transformHit ( ) const

inherited

Local (simulation/hit frame) to global transform.

updateCache()

void InDetDD::SolidStateDetectorElementBase::updateCache ( ) const

virtualinherited

Recalculate cached values.

Definition at line 272 of file SolidStateDetectorElementBase.cxx.

  {
    if (m_cache.isValid()) return;
    CachedVals cache;
 
    // use aligned transform if available
    const GeoTrf::Transform3D* ptrXf;
    GeoTrf::Transform3D geotrf;
 
    if (m_geoAlignStore){
      ptrXf = m_geoAlignStore->getAbsPosition(getMaterialGeom());
      if (ptrXf) {
         cache.m_transformHit = (*ptrXf);
         //the below can have an (optional) shift to acccount for transformations
         //applied when constructing the DetectorElements from the GeoPhysVols
         //(At the moment, only used for ITkStrip barrel when creating a DetElement per row)
         geotrf = (*ptrXf) * m_design->moduleShift();
      }
      else {
        geotrf.setIdentity();
      }
    }
    else{
      cache.m_transformHit  = getMaterialGeom()->getAbsoluteTransform();
      //the below can have an (optional) shift to acccount for transformations
      //applied when constructing the DetectorElements from the GeoPhysVols
      //(At the moment, only used for ITkStrip barrel when creating a DetElement per row)
      geotrf = getMaterialGeom()->getAbsoluteTransform() * m_design->moduleShift();
    }
 
    const GeoTrf::Transform3D& geoTransform = geotrf;
 
    cache.m_center = geoTransform * m_design->sensorCenter();
 
    //Is this needed outside e.g. ReadSiDetElements? Maybe candidate for future removal?
    cache.m_centerCLHEP = HepGeom::Point3D<double>(cache.m_center[0],cache.m_center[1],cache.m_center[2]);
 
    Amg::Vector3D centerGeoModel(0., 0., 0.);
    cache.m_origin = geoTransform * centerGeoModel;
 
    //
    // Determine directions depth, eta and phi axis in reconstruction local frame
    // ie depth away from interaction point
    //    phi in direction of increasing phi
    //    eta in direction of increasing z in barrel, and increasing r in endcap
    //
 
    // depthAxis, phiAxis, and etaAxis are defined to be x,y,z respectively for all detectors for hit local frame.
    // depthAxis, phiAxis, and etaAxis are defined to be z,x,y respectively for all detectors for reco local frame.
    static const Amg::Vector3D localAxes[3] = {
        Amg::Vector3D(1.,0.,0.),
        Amg::Vector3D(0.,1.,0.),
        Amg::Vector3D(0.,0.,1.)
    };
 
    static const Amg::Vector3D & localRecoPhiAxis = localAxes[distPhi];     // Defined to be same as x axis
    static const Amg::Vector3D & localRecoEtaAxis = localAxes[distEta];     // Defined to be same as y axis
    static const Amg::Vector3D & localRecoDepthAxis = localAxes[distDepth]; // Defined to be same as z axis
 
    // We only need to calculate the rough orientation once.
    //For it to change would require extreme unrealistic misalignment changes.
#ifndef NDEBUG
    bool setAxisDir = false;
#endif
    if (!m_axisDir.isValid()) {
        AxisDir dir{};
        // Determine the unit vectors in global frame
 
        const Amg::Vector3D &geoModelPhiAxis = localAxes[m_hitPhi];
        const Amg::Vector3D &geoModelEtaAxis = localAxes[m_hitEta];
        const Amg::Vector3D &geoModelDepthAxis = localAxes[m_hitDepth];
 
        Amg::Vector3D globalDepthAxis(geoTransform.linear() * geoModelDepthAxis);
        Amg::Vector3D globalPhiAxis(geoTransform.linear() * geoModelPhiAxis);
        Amg::Vector3D globalEtaAxis(geoTransform.linear() * geoModelEtaAxis);
 
        // unit radial vector
        Amg::Vector3D unitR(cache.m_center.x(), cache.m_center.y(), 0.);
        unitR.normalize();
 
        Amg::Vector3D nominalEta(0.0,0.0,1.0);
        Amg::Vector3D nominalNormal(0.0,0.0,0.0);
        Amg::Vector3D nominalPhi(-unitR.y(), unitR.x(), 0.0);
 
        // In Barrel like geometry, the etaAxis is along increasing z, and normal is in increasing radial direction.
        // In Endcap like geometry, the etaAxis is along increasing r, and normal is in decreasing z direction,
        // We base whether it is barrel like or endcap like by the orientation of the local z axis of the
        // the element. This allows the use of endcap identifiers in a TB setup.
 
 
        nominalEta(2) = 1.0;
        dir.m_barrelLike = true;
 
        if (std::abs(globalEtaAxis.dot(nominalEta)) < 0.5) { // Check that it is in roughly the right direction. Allowed not to be for ITK inclined/barrel ring modules
            dir.m_barrelLike = false;
        }
 
        if (dir.m_barrelLike) {
          nominalEta(2) = 1.0;
          nominalNormal =  unitR;
        } else { // endcap like
          nominalNormal(2) = -1.0;
          nominalEta = unitR;
        }
 
        // Determine if axes are to have there directions swapped.
 
        //
        // Depth axis.
        //
        dir.m_depthAngle = globalDepthAxis.dot(nominalNormal);
        dir.m_depthDirection = true;
        if (dir.m_depthAngle < 0) {
            if (m_design->depthSymmetric()) {
              dir.m_depthDirection = false;
            } else {
          ATH_MSG_DEBUG( "Unable to swap local depth axis.");
            }
        }
        if (std::abs(dir.m_depthAngle) < 0.5 && (m_design->type())!=InDetDD::PLR) { // Check that it is in roughly the right direction. Ignore this for the PLR as it gives values as low as 0.489)
            ATH_MSG_ERROR( "Orientation of local depth axis does not follow correct convention.");
            dir.m_depthDirection = true; // Don't swap.
        }
 
        // for HGTD modules, the check on phi and eta directions don't make sense
        // as the modules do not respect the conventional position for endcap discs:
        // - the local eta axis is never parallel to the radial direction
        // - the local phi axis is never perpendicular to the radial direction
        // hence, removing errors and allowing swap of the axis when needed
        bool isHGTD = this->getIdHelper()->is_hgtd(m_id);
 
        //
        // Phi axis
        //
        dir.m_phiAngle = globalPhiAxis.dot(nominalPhi);
        dir.m_phiDirection = true;
        if (dir.m_phiAngle < 0) {
            if (m_design->phiSymmetric()) {
                dir.m_phiDirection = false;
            } else {
          ATH_MSG_DEBUG("Unable to swap local xPhi axis.");
            }
        }
        if (not isHGTD and std::abs(dir.m_phiAngle) < 0.5) { // Check that it is in roughly the right direction.
      ATH_MSG_ERROR( "Orientation of local xPhi axis does not follow correct convention.");
          dir.m_phiDirection = true; // Don't swap.
        }
 
        //
        // Eta axis
        //
        dir.m_etaAngle = globalEtaAxis.dot(nominalEta);
        dir.m_etaDirection = true;
        if (dir.m_etaAngle < 0) {
            if (m_design->etaSymmetric()) {
                dir.m_etaDirection = false;
            } else {
          ATH_MSG_DEBUG("Unable to swap local xEta axis.");
            }
        }
        if (not isHGTD and std::abs(dir.m_etaAngle) < 0.5) { // Check that it is in roughly the right direction.
      ATH_MSG_ERROR( "Orientation of local xEta axis does not follow correct convention.");
          dir.m_etaDirection = true; // Don't swap
        }
 
        m_axisDir.set (dir);
#ifndef NDEBUG
        setAxisDir = true;
#endif
    } // end if (!m_axisDir.isValid())
 
    cache.m_transformCLHEP = Amg::EigenTransformToCLHEP(geoTransform) * recoToHitTransform();
    cache.m_transform = Amg::CLHEPTransformToEigen(cache.m_transformCLHEP);
 
    #ifndef NDEBUG
    // Check that local frame is right-handed. (ie transform has no reflection)
    // This can be done by checking that the determinant is >0.
    if (setAxisDir) { // Only need to check this once.
        const AxisDir& dir = *m_axisDir.ptr();
        HepGeom::Transform3D & t = cache.m_transformCLHEP;
        double det = t(0,0) * (t(1,1)*t(2,2) - t(1,2)*t(2,1)) -
                     t(0,1) * (t(1,0)*t(2,2) - t(1,2)*t(2,0)) +
                     t(0,2) * (t(1,0)*t(2,1) - t(1,1)*t(2,0));
        if (det < 0) {
      ATH_MSG_DEBUG( "Local frame is left-handed. (hitEtaDirection, hitPhiDirection, hitDepthDirection) = ("
             << dir.m_etaDirection <<", "
             << dir.m_phiDirection <<", "
             << dir.m_depthDirection <<")");
        }
    }
    #endif
 
    // Initialize various cached members, needs to be done here otherwise the necessary transforms are not yet initialized
    // The unit vectors
    cache.m_normal = cache.m_transform.linear() * localRecoDepthAxis;
 
    cache.m_phiAxis = cache.m_transform.linear() * localRecoPhiAxis;
    cache.m_etaAxis = cache.m_transform.linear() * localRecoEtaAxis;
 
    //Check where these are actually needed - candidates for removal?
    cache.m_phiAxisCLHEP = HepGeom::Vector3D<double>(cache.m_phiAxis[0],cache.m_phiAxis[1],cache.m_phiAxis[2]);
    cache.m_etaAxisCLHEP = HepGeom::Vector3D<double>(cache.m_etaAxis[0],cache.m_etaAxis[1],cache.m_etaAxis[2]);
 
    getExtent(cache);
    m_cache.set (std::move (cache));
  }

width()

double InDetDD::SolidStateDetectorElementBase::width ( ) const

inherited

Methods from design (inline)

Width in phi direction. For the SCT endcap it returns the average width.

zMax()

double InDetDD::SolidStateDetectorElementBase::zMax ( ) const

inherited

zMin()

double InDetDD::SolidStateDetectorElementBase::zMin ( ) const

inherited

Member Data Documentation

m_axisDir

CxxUtils::CachedValue<AxisDir> InDetDD::SolidStateDetectorElementBase::m_axisDir

protectedinherited

Definition at line 675 of file SolidStateDetectorElementBase.h.

m_cache

CxxUtils::CachedValue<CachedVals> InDetDD::SolidStateDetectorElementBase::m_cache

protectedinherited

Definition at line 702 of file SolidStateDetectorElementBase.h.

m_commonItems

const SiCommonItems* InDetDD::SolidStateDetectorElementBase::m_commonItems {nullptr}

protectedinherited

Definition at line 624 of file SolidStateDetectorElementBase.h.

{nullptr};

m_design

const DetectorDesign* InDetDD::SolidStateDetectorElementBase::m_design {nullptr}

protectedinherited

local description of this detector element

Definition at line 622 of file SolidStateDetectorElementBase.h.

{nullptr};

m_geoAlignStore

const GeoAlignmentStore* InDetDD::SolidStateDetectorElementBase::m_geoAlignStore {}

protectedinherited

Definition at line 626 of file SolidStateDetectorElementBase.h.

{};

m_hitDepth

DetectorDesign::Axis InDetDD::SolidStateDetectorElementBase::m_hitDepth

protectedinherited

Definition at line 633 of file SolidStateDetectorElementBase.h.

m_hitEta

DetectorDesign::Axis InDetDD::SolidStateDetectorElementBase::m_hitEta

protectedinherited

Axes.

Definition at line 631 of file SolidStateDetectorElementBase.h.

m_hitPhi

DetectorDesign::Axis InDetDD::SolidStateDetectorElementBase::m_hitPhi

protectedinherited

Definition at line 632 of file SolidStateDetectorElementBase.h.

m_id

Identifier InDetDD::SolidStateDetectorElementBase::m_id {}

protectedinherited

identifier of this detector element

Definition at line 617 of file SolidStateDetectorElementBase.h.

{};

m_idHash

IdentifierHash InDetDD::SolidStateDetectorElementBase::m_idHash {}

protectedinherited

hash id of this detector element

Definition at line 643 of file SolidStateDetectorElementBase.h.

{};

m_isBarrel

bool InDetDD::SiDetectorElement::m_isBarrel {false}

protected

Definition at line 469 of file SiDetectorElement.h.

{false};

m_isDBM

bool InDetDD::SiDetectorElement::m_isDBM {false}

protected

Definition at line 468 of file SiDetectorElement.h.

{false};

m_isPixel

bool InDetDD::SiDetectorElement::m_isPixel {false}

protected

Definition at line 465 of file SiDetectorElement.h.

{false};

m_isPLR

bool InDetDD::SiDetectorElement::m_isPLR {false}

protected

Definition at line 467 of file SiDetectorElement.h.

{false};

m_isSCT

bool InDetDD::SiDetectorElement::m_isSCT {false}

protected

Definition at line 466 of file SiDetectorElement.h.

{false};

m_isStereo

CxxUtils::CachedValue<bool> InDetDD::SiDetectorElement::m_isStereo

protected

Definition at line 487 of file SiDetectorElement.h.

m_nextInEta

const SiDetectorElement* InDetDD::SiDetectorElement::m_nextInEta {nullptr}

protected

Definition at line 476 of file SiDetectorElement.h.

{nullptr}; // set by setNextInEta

m_nextInPhi

const SiDetectorElement* InDetDD::SiDetectorElement::m_nextInPhi {nullptr}

protected

Definition at line 478 of file SiDetectorElement.h.

{nullptr}; // set by setNextInPhi

m_otherSide

const SiDetectorElement* InDetDD::SiDetectorElement::m_otherSide {nullptr}

protected

Definition at line 480 of file SiDetectorElement.h.

{nullptr}; // set by setOtherSide

m_prevInEta

const SiDetectorElement* InDetDD::SiDetectorElement::m_prevInEta {nullptr}

protected

Definition at line 477 of file SiDetectorElement.h.

{nullptr}; // set by setPrevInEta

m_prevInPhi

const SiDetectorElement* InDetDD::SiDetectorElement::m_prevInPhi {nullptr}

protected

Definition at line 479 of file SiDetectorElement.h.

{nullptr}; // set by setPrevInPhi

m_siDesign

const SiDetectorDesign* InDetDD::SiDetectorElement::m_siDesign

protected

Definition at line 458 of file SiDetectorElement.h.

m_surface

std::unique_ptr<Trk::Surface> InDetDD::SolidStateDetectorElementBase::m_surface

protectedinherited

Definition at line 625 of file SolidStateDetectorElementBase.h.

m_surfaces

CxxUtils::CachedValue<std::vector<const Trk::Surface*> > InDetDD::SiDetectorElement::m_surfaces

protected

Definition at line 451 of file SiDetectorElement.h.

---

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

• SiDetectorElement.h
• SiDetectorElement.cxx