InDetDD::SolidStateDetectorElementBase Class Reference

Class to hold geometrical description of a solid state detector element. More...

#include <SolidStateDetectorElementBase.h>

Inheritance diagram for InDetDD::SolidStateDetectorElementBase:

Collaboration diagram for InDetDD::SolidStateDetectorElementBase:

Classes
struct	AxisDir
struct	CachedVals

Public Member Functions
	SolidStateDetectorElementBase (const Identifier &id, const DetectorDesign *design, const GeoVFullPhysVol *geophysvol, const SiCommonItems *commonItems, const GeoAlignmentStore *geoAlignStore=nullptr)
	Constructor with parameters.
virtual	~SolidStateDetectorElementBase ()
	Destructor.
	SolidStateDetectorElementBase ()=delete
	Don't allow no-argument constructor.
	SolidStateDetectorElementBase (const SolidStateDetectorElementBase &)=delete
	Don't allow copy constructor.
SolidStateDetectorElementBase &	operator= (const SolidStateDetectorElementBase &)=delete
	Don't allow assignment operator.
	SolidStateDetectorElementBase (SolidStateDetectorElementBase &&)=delete
	Don't allow move constructor.
SolidStateDetectorElementBase &	operator= (SolidStateDetectorElementBase &&)=delete
	Don't allow move assignment operator.
std::string	detectorTypeString () const
	Returns a string of the Detector element type.
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)
virtual Identifier	identifierFromCellId (const SiCellId &cellId) const =0
	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 =0
	SiCellId from Identifier.
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.
virtual double	get_rz () const =0
Design methods
virtual const DetectorDesign &	design () const
	access to the local description (inline):
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
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
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
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 a solid state detector element.

This base class is shared between the InnerDetetector/ITk/HGTD since there is a lot of commonality.

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
• Transformation/Orientation
• Module Frame
• Element Extent
• Design methods
• Intersection Tests
• Lorentz Correction
• Readout cell id
• Miscellaneous
• Cache handling.

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

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 132 of file SolidStateDetectorElementBase.h.

Constructor & Destructor Documentation

SolidStateDetectorElementBase() [1/4]

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

Constructor with parameters.

Definition at line 20 of file SolidStateDetectorElementBase.cxx.

                                                                               :
    TrkDetElementBase(geophysvol),
    m_id(id),
    m_design(design),
    m_commonItems(commonItems),
    m_surface(),
    m_geoAlignStore(geoAlignStore)
  {
    m_hitEta = m_design->etaAxis();
    m_hitPhi = m_design->phiAxis();
    m_hitDepth = m_design->depthAxis();
 
    commonConstructor();
  }

~SolidStateDetectorElementBase()

InDetDD::SolidStateDetectorElementBase::~SolidStateDetectorElementBase ( )

virtualdefault

Destructor.

SolidStateDetectorElementBase() [2/4]

InDetDD::SolidStateDetectorElementBase::SolidStateDetectorElementBase ( )

delete

Don't allow no-argument constructor.

SolidStateDetectorElementBase() [3/4]

InDetDD::SolidStateDetectorElementBase::SolidStateDetectorElementBase ( const SolidStateDetectorElementBase & )

delete

Don't allow copy constructor.

SolidStateDetectorElementBase() [4/4]

InDetDD::SolidStateDetectorElementBase::SolidStateDetectorElementBase ( SolidStateDetectorElementBase && )

delete

Don't allow move constructor.

Member Function Documentation

bounds() [1/2]

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

finaloverridevirtual

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

finaloverridevirtual

TrkDetElementBase interface (inline)

Implements Trk::TrkDetElementBase.

carrierType()

InDetDD::CarrierType InDetDD::SolidStateDetectorElementBase::carrierType

(

)

const

carrier type for readout.

ie holes for SCT and electrons for pixels. inline

cellIdFromIdentifier()

virtual SiCellId InDetDD::SolidStateDetectorElementBase::cellIdFromIdentifier ( const Identifier & identifier ) const

pure virtual

SiCellId from Identifier.

Implemented in HGTD_DetectorElement, InDetDD::HGTD_DetectorElement, InDetDD::SiDetectorElement, and SiDetectorElement.

cellIdOfPosition()

SiCellId InDetDD::SolidStateDetectorElementBase::cellIdOfPosition

(

const Amg::Vector2D &

localPos

)

const

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

finaloverridevirtual

Center in global coordinates.

Implements Trk::TrkDetElementBase.

center() [2/2]

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

finaloverridevirtual

TrkDetElementBase interface (inline)

Implements Trk::TrkDetElementBase.

commonConstructor()

void InDetDD::SolidStateDetectorElementBase::commonConstructor ( )

private

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

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

Definition at line 60 of file SolidStateDetectorElementBase.cxx.

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

defTransformCLHEP()

const HepGeom::Transform3D InDetDD::SolidStateDetectorElementBase::defTransformCLHEP

(

)

const

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

depthDirection()

bool InDetDD::SolidStateDetectorElementBase::depthDirection

(

)

const

Return information on orientation.

design()

virtual const DetectorDesign & InDetDD::SolidStateDetectorElementBase::design ( ) const

virtual

access to the local description (inline):

Reimplemented in HGTD_DetectorElement, InDetDD::HGTD_DetectorElement, InDetDD::SiDetectorElement, and SiDetectorElement.

detectorType()

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

overridevirtual

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

etaAxis()

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

(

)

const

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

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

etaPitch()

double InDetDD::SolidStateDetectorElementBase::etaPitch

(

)

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.

get_rz()

virtual double InDetDD::SolidStateDetectorElementBase::get_rz ( ) const

pure virtual

Implemented in HGTD_DetectorElement, InDetDD::HGTD_DetectorElement, InDetDD::SiDetectorElement, and SiDetectorElement.

getCommonItems()

const SiCommonItems * InDetDD::SolidStateDetectorElementBase::getCommonItems

(

)

const

getCorners()

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

private

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

private

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

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

private

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

Returns the id helper (inline)

globalPosition() [1/3]

Amg::Vector3D InDetDD::SolidStateDetectorElementBase::globalPosition

(

const Amg::Vector2D &

localPos

)

const

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

globalPosition() [2/3]

Amg::Vector3D InDetDD::SolidStateDetectorElementBase::globalPosition

(

const Amg::Vector3D &

localPos

)

const

globalPosition() [3/3]

HepGeom::Point3D< double > InDetDD::SolidStateDetectorElementBase::globalPosition

(

const HepGeom::Point3D< double > &

localPos

)

const

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

globalPositionHit() [1/2]

Amg::Vector3D InDetDD::SolidStateDetectorElementBase::globalPositionHit

(

const Amg::Vector3D &

simulationLocalPos

)

const

globalPositionHit() [2/2]

HepGeom::Point3D< double > InDetDD::SolidStateDetectorElementBase::globalPositionHit

(

const HepGeom::Point3D< double > &

simulationLocalPos

)

const

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

hitDepthDirection()

double InDetDD::SolidStateDetectorElementBase::hitDepthDirection

(

)

const

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

Returns +/-1. inline

hitEtaDirection()

double InDetDD::SolidStateDetectorElementBase::hitEtaDirection

(

)

const

See previous method.

inline

hitLocalToLocal()

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

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

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

See previous method.

inline

identifierFromCellId()

virtual Identifier InDetDD::SolidStateDetectorElementBase::identifierFromCellId ( const SiCellId & cellId ) const

pure virtual

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

Implemented in InDetDD::HGTD_DetectorElement, InDetDD::SiDetectorElement, and SiDetectorElement.

identifierOfPosition()

Identifier InDetDD::SolidStateDetectorElementBase::identifierOfPosition

(

const Amg::Vector2D &

localPos

)

const

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

finaloverridevirtual

identifier of this detector element (inline)

Implements Trk::TrkDetElementBase.

identifyHash()

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

finaloverridevirtual

identifier hash (inline)

Implements Trk::TrkDetElementBase.

inDetector() [1/2]

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

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

Definition at line 211 of file SolidStateDetectorElementBase.cxx.

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

invalidate()

void InDetDD::SolidStateDetectorElementBase::invalidate

(

)

Signal that cached values are no longer valid.

Invalidate general cache (inline)

length()

double InDetDD::SolidStateDetectorElementBase::length

(

)

const

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

localPosition() [1/2]

Amg::Vector2D InDetDD::SolidStateDetectorElementBase::localPosition

(

const Amg::Vector3D &

globalPosition

)

const

localPosition() [2/2]

Amg::Vector2D InDetDD::SolidStateDetectorElementBase::localPosition

(

const HepGeom::Point3D< double > &

globalPosition

)

const

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

maxWidth()

double InDetDD::SolidStateDetectorElementBase::maxWidth

(

)

const

Max width.

Needed for the SCT endcap.

minWidth()

double InDetDD::SolidStateDetectorElementBase::minWidth

(

)

const

Min width.

Needed for the SCT endcap.

msg()

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

private

Declaring the Message method for further use (inline)

msgLvl()

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

private

Declaring the Method providing Verbosity Level (inline)

normal() [1/2]

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

finaloverridevirtual

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

finaloverridevirtual

TrkDetElementBase interface (inline)

Implements Trk::TrkDetElementBase.

numberOfConnectedCells()

int InDetDD::SolidStateDetectorElementBase::numberOfConnectedCells

(

const SiCellId

cellId

)

const

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]

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

delete

Don't allow assignment operator.

operator=() [2/2]

SolidStateDetectorElementBase & InDetDD::SolidStateDetectorElementBase::operator= ( SolidStateDetectorElementBase && )

delete

Don't allow move assignment operator.

phiAngle()

double InDetDD::SolidStateDetectorElementBase::phiAngle

(

)

const

phiAxis()

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

(

)

const

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

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

phiMax()

double InDetDD::SolidStateDetectorElementBase::phiMax

(

)

const

phiMin()

double InDetDD::SolidStateDetectorElementBase::phiMin

(

)

const

phiPitch()

double InDetDD::SolidStateDetectorElementBase::phiPitch

(

)

const

rawLocalPositionOfCell() [1/2]

Amg::Vector2D InDetDD::SolidStateDetectorElementBase::rawLocalPositionOfCell

(

const Identifier &

id

)

const

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

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

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

rMin()

double InDetDD::SolidStateDetectorElementBase::rMin

(

)

const

surface() [1/3]

Trk::Surface & InDetDD::SolidStateDetectorElementBase::surface

(

)

Element Surface.

surface() [2/3]

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

finaloverridevirtual

Return surface associated with this detector element.

Implements Trk::TrkDetElementBase.

surface() [3/3]

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

finaloverridevirtual

TrkDetElementBase interface (inline)

Implements Trk::TrkDetElementBase.

swapEtaReadoutDirection()

bool InDetDD::SolidStateDetectorElementBase::swapEtaReadoutDirection

(

)

const

For eta_index (only relevant for pixel) (inline)

swapPhiReadoutDirection()

bool InDetDD::SolidStateDetectorElementBase::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

transform() [1/2]

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

finaloverridevirtual

Return local to global transform.

Implements Trk::TrkDetElementBase.

transform() [2/2]

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

finaloverridevirtual

TrkDetElementBase interface (inline)

Implements Trk::TrkDetElementBase.

transformCLHEP()

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

(

)

const

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

Local (simulation/hit frame) to global transform.

updateCache()

void InDetDD::SolidStateDetectorElementBase::updateCache ( ) const

virtual

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

Methods from design (inline)

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

zMax()

double InDetDD::SolidStateDetectorElementBase::zMax

(

)

const

zMin()

double InDetDD::SolidStateDetectorElementBase::zMin

(

)

const

Member Data Documentation

m_axisDir

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

protected

Definition at line 675 of file SolidStateDetectorElementBase.h.

m_cache

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

protected

Definition at line 702 of file SolidStateDetectorElementBase.h.

m_commonItems

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

protected

Definition at line 624 of file SolidStateDetectorElementBase.h.

{nullptr};

m_design

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

protected

local description of this detector element

Definition at line 622 of file SolidStateDetectorElementBase.h.

{nullptr};

m_geoAlignStore

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

protected

Definition at line 626 of file SolidStateDetectorElementBase.h.

{};

m_hitDepth

DetectorDesign::Axis InDetDD::SolidStateDetectorElementBase::m_hitDepth

protected

Definition at line 633 of file SolidStateDetectorElementBase.h.

m_hitEta

DetectorDesign::Axis InDetDD::SolidStateDetectorElementBase::m_hitEta

protected

Axes.

Definition at line 631 of file SolidStateDetectorElementBase.h.

m_hitPhi

DetectorDesign::Axis InDetDD::SolidStateDetectorElementBase::m_hitPhi

protected

Definition at line 632 of file SolidStateDetectorElementBase.h.

m_id

Identifier InDetDD::SolidStateDetectorElementBase::m_id {}

protected

identifier of this detector element

Definition at line 617 of file SolidStateDetectorElementBase.h.

{};

m_idHash

IdentifierHash InDetDD::SolidStateDetectorElementBase::m_idHash {}

protected

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

protected

Definition at line 625 of file SolidStateDetectorElementBase.h.

---

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

• SolidStateDetectorElementBase.h
• SolidStateDetectorElementBase.cxx