InDetDD::HGTD_DetectorElement Class Reference

Class to hold geometrical description of an HGTD detector element. More...

#include <HGTD_DetectorElement.h>

Inheritance diagram for InDetDD::HGTD_DetectorElement:

Collaboration diagram for InDetDD::HGTD_DetectorElement:

Public Member Functions
	HGTD_DetectorElement (const Identifier &id, const HGTD_ModuleDesign *design, const GeoVFullPhysVol *geophysvol, const SiCommonItems *commonItems, const GeoAlignmentStore *geoAlignStore=nullptr)
	Constructor:
virtual	~HGTD_DetectorElement ()
	Destructor:
std::string	detectorTypeString () const
	Returns a string of the Detector element type.
Identification
Methods to identify the element and identifier manipulation.
Identifier	identifierFromCellId (const SiCellId &cellId) const override final
	Identifier from SiCellId (ie pixel eta_index,phi_index)
SiCellId	cellIdFromIdentifier (const Identifier &identifier) const override final
	SiCellId from Identifier.
Transformation/Orientation
virtual const std::vector< const Trk::Surface * > &	surfaces () const
	Returns the full list of surfaces associated to this detector element.
Element Extent
Methods to get extent of element in r,phi and z.
double	get_rz () const override final
Design methods
const HGTD_ModuleDesign &	design () const override final
	access to the local description:
InDetDD::CarrierType	carrierType () const
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)
virtual Trk::DetectorElemType	detectorType () const override
	TrkDetElementBase interface detectorType.
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)
double	phiPitch () const
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
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
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.
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.
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 an HGTD detector element.

A detector element in this case is a single module in the HGTD.

Definition at line 40 of file HGTD_DetectorElement.h.

Constructor & Destructor Documentation

HGTD_DetectorElement()

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

Constructor:

Definition at line 17 of file HGTD_DetectorElement.cxx.

                                                                                   :
  SolidStateDetectorElementBase(id, design, geophysvol, commonItems, geoAlignStore)
{
    const HGTD_ID* hgtdId = dynamic_cast<const HGTD_ID *>(getIdHelper());
    m_idHash = hgtdId->wafer_hash(m_id);
 
    if (!m_idHash.is_valid()) throw std::runtime_error("HGTD_DetectorElement: Unable to set IdentifierHash");
}

~HGTD_DetectorElement()

InDetDD::HGTD_DetectorElement::~HGTD_DetectorElement ( )

virtual

Destructor:

Definition at line 31 of file HGTD_DetectorElement.cxx.

{
}

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::HGTD_DetectorElement::carrierType ( ) const

inline

Definition at line 126 of file HGTD_DetectorElement.h.

{
  return design().carrierType();
}

cellIdFromIdentifier()

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

finaloverridevirtual

SiCellId from Identifier.

Implements InDetDD::SolidStateDetectorElementBase.

Definition at line 51 of file HGTD_DetectorElement.cxx.

{
    SiCellId cellId; // Initialized in invalid state.
    // If something fails it returns the cellId in an invalid state.
    if (identifier.is_valid()) {
        const HGTD_ID * hgtdIdHelper = static_cast<const HGTD_ID *>(getIdHelper());
        if (hgtdIdHelper) {
            cellId = SiCellId(hgtdIdHelper->phi_index(identifier), hgtdIdHelper->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::SolidStateDetectorElementBase::commonConstructor ( )

privateinherited

Common code for constructors.

Definition at line 257 of file SolidStateDetectorElementBase.cxx.

  {
    //Do we need this at all in base class, or just in derived?
    if (!m_id.is_valid()) throw std::runtime_error("SolidStateDetectorElementBase: Invalid identifier");
 
    // 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);
  }

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()

const HGTD_ModuleDesign & InDetDD::HGTD_DetectorElement::design ( ) const

inlinefinaloverridevirtual

access to the local description:

Reimplemented from InDetDD::SolidStateDetectorElementBase.

Definition at line 121 of file HGTD_DetectorElement.h.

{
  return static_cast<const HGTD_ModuleDesign&> (SolidStateDetectorElementBase::design());
}

detectorType()

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

overridevirtualinherited

TrkDetElementBase interface detectorType.

Implements Trk::TrkDetElementBase.

Reimplemented in InDetDD::SiDetectorElement, and SiDetectorElement.

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

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.

get_rz()

double InDetDD::HGTD_DetectorElement::get_rz ( ) const

finaloverridevirtual

Implements InDetDD::SolidStateDetectorElementBase.

Definition at line 80 of file HGTD_DetectorElement.cxx.

{
    // Calculate z (endcap)
    // Use center of sensor ((0,0,0) in local coordinates) for determining this.
    return center().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::HGTD_DetectorElement::identifierFromCellId ( const SiCellId & cellId ) const

finaloverridevirtual

Identifier from SiCellId (ie pixel eta_index,phi_index)

Implements InDetDD::SolidStateDetectorElementBase.

Definition at line 36 of file HGTD_DetectorElement.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 HGTD_ID * hgtdIdHelper = static_cast<const HGTD_ID *>(getIdHelper());
        if (hgtdIdHelper) {
            id = hgtdIdHelper->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)

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)

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.

msg()

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

privateinherited

Declaring the Message method for further use (inline)

msgLvl()

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

privateinherited

Declaring the Method providing Verbosity Level (inline)

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

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()

double InDetDD::SolidStateDetectorElementBase::phiPitch ( ) const

inherited

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

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::HGTD_DetectorElement::surfaces ( ) const

virtual

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

Definition at line 65 of file HGTD_DetectorElement.cxx.

{
    // This method is needed to satisfy inheritance from TrkDetElementBase
    // so just return the one surface
    if (!m_surfaces.isValid()) {
      std::vector<const Trk::Surface*> s;
      // get this surface
      s.push_back(&surface());
      m_surfaces.set (std::move (s));
    }
 
    // return the surfaces
    return *m_surfaces.ptr();
}

swapEtaReadoutDirection()

bool InDetDD::SolidStateDetectorElementBase::swapEtaReadoutDirection ( ) const

inherited

For eta_index (only relevant for pixel) (inline)

swapPhiReadoutDirection()

bool InDetDD::SolidStateDetectorElementBase::swapPhiReadoutDirection ( ) const

inherited

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_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::HGTD_DetectorElement::m_surfaces

protected

Definition at line 113 of file HGTD_DetectorElement.h.

---

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

• HGTD_DetectorElement.h
• HGTD_DetectorElement.cxx