#include <StripStereoAnnulusDesign.h>

Inheritance diagram for InDetDD::StripStereoAnnulusDesign:

Collaboration diagram for InDetDD::StripStereoAnnulusDesign:

Public Types
enum	Axis { xAxis =0 , yAxis , zAxis }

Public Member Functions
	StripStereoAnnulusDesign (const SiDetectorDesign::Axis &stripDirection, const SiDetectorDesign::Axis &thicknessDirection, const double &thickness, const int &readoutSide, const InDetDD::CarrierType &carrier, const int &nRows, const std::vector< int > &nStrips, const std::vector< double > &pitch, const std::vector< double > &stripStart, const std::vector< double > &stripEnd, const double &stereoAngle, const double &centreR, const bool &usePC, InDetDD::DetectorType detectorType=InDetDD::Undefined)
	StripStereoAnnulusDesign (const SiDetectorDesign::Axis &stripDirection, const SiDetectorDesign::Axis &thicknessDirection, const double &thickness, const int &readoutSide, const InDetDD::CarrierType &carrier, const int &nRows, const std::vector< int > &nStrips, const std::vector< double > &pitch, const std::vector< double > &stripStart, const std::vector< double > &stripEnd, const double &stereoAngle, const double &centreR, const double &waferCentreR, const bool &usePC, InDetDD::DetectorType detectorType=InDetDD::Undefined)
	~StripStereoAnnulusDesign ()=default
SiLocalPosition	beamToStrip (const SiLocalPosition &pos) const
SiLocalPosition	beamToStripPC (const SiLocalPosition &pos) const
SiLocalPosition	beamToStripPCpolar (const SiLocalPosition &pos) const
	Version of StripStereoAnnulusDesign::beamToStripPC transform based exclusively in a polar system.
SiLocalPosition	stripToBeam (const SiLocalPosition &pos) const
SiLocalPosition	stripToBeamPC (const SiLocalPosition &pos) const
SiLocalPosition	stripToBeamPCpolar (const SiLocalPosition &pos) const
	Version of StripStereoAnnulusDesign::stripToBeamPC transform based exclusively in a polar system.
virtual Amg::Vector3D	sensorCenter () const override
	Return the centre of a sensor in the local reference frame.
	StripStereoAnnulusDesign (const StripStereoAnnulusDesign &design)
StripStereoAnnulusDesign &	operator= (const StripStereoAnnulusDesign &design)
std::pair< int, int >	getStripRow (SiCellId cellId) const final
	Get the strip and row number of the cell.
virtual int	strip1Dim (int strip, int row) const override
	only relevant for SCT.
SiLocalPosition	stripPosAtR (int strip, int row, double r) const
virtual int	diodesInRow (const int row) const override
virtual void	distanceToDetectorEdge (const SiLocalPosition &localPosition, double &etaDist, double &phiDist) const override
	Returns distance to nearest detector active edge +ve = inside -ve = outside.
virtual bool	inActiveArea (const SiLocalPosition &chargePos, bool checkBondGap=true) const override
	check if the position is in active area
virtual const Trk::SurfaceBounds &	bounds () const override
	Get a reference to the module bounds object.
virtual std::pair< SiLocalPosition, SiLocalPosition >	endsOfStrip (const SiLocalPosition &position) const override
	Give end points of the strip that covers the given position.
virtual double	stripPitch (const SiLocalPosition &localPosition) const override
	give the strip pitch (dependence on position needed for forward)
double	stripPitch (const SiCellId &cellId) const
virtual double	stripPitch () const override
	give the strip pitch (For Forward returns pitch at center)
virtual double	phiPitch (const SiLocalPosition &localPosition) const override
	Pitch in phi direction.
double	phiPitch (const SiCellId &cellId) const
virtual double	phiPitch () const override
	Pitch in phi direction.
double	phiPitchPhi (const SiLocalPosition &localPosition) const
double	phiPitchPhi (const SiCellId &cellId) const
double	phiPitchPhi () const
virtual double	scaledDistanceToNearestDiode (const SiLocalPosition &chargePos) const override
	give distance to the nearest diode in units of pitch, from 0.0 to 0.5, this method should be fast as it is called for every surface charge in the SCT_SurfaceChargesGenerator
virtual SiDiodesParameters	parameters (const SiCellId &cellId) const override
	Return strip width, centre, length etc. Hard to find if this is used or not.
virtual SiLocalPosition	localPositionOfCell (const SiCellId &cellId) const override
	id -> position
virtual SiLocalPosition	localPositionOfCluster (const SiCellId &cellId, int clusterSize) const override
SiLocalPosition	localPositionOfCellPC (const SiCellId &cellId) const
	This is for debugging only.
SiLocalPosition	localPositionOfClusterPC (const SiCellId &cellId, int clusterSize) const
	This is for debugging only.
virtual SiCellId	cellIdOfPosition (const SiLocalPosition &localPos) const override
	position -> id
SiLocalPosition	positionFromStrip (const SiCellId &cellId) const
virtual SiLocalPosition	positionFromStrip (const int stripNumber) const override
	gives position of strip center ALTERNATIVE/PREFERED way is to use localPositionOfCell(const SiCellId & cellId) or rawLocalPositionOfCell method in SiDetectorElement.
virtual int	row (int stripId1Dim) const override
virtual int	strip (int stripId1Dim) const override
virtual void	neighboursOfCell (const SiCellId &cellId, std::vector< SiCellId > &neighbours) const override
	Get the neighbouring diodes of a given diode: Cell for which the neighbours must be found List of cells which are neighbours of the given one.
virtual SiCellId	cellIdInRange (const SiCellId &) const override
	DEPRECATED: only used in a stupid example (2014) Check if cell is in range.
virtual SiReadoutCellId	readoutIdOfCell (const SiCellId &cellId) const override
	diode id -> readout id
double	waferCentreR () const
double	centreR () const
InDetDD::DetectorType	type () const override final
	Type of element.
virtual HepGeom::Vector3D< double >	phiMeasureSegment (const SiLocalPosition &position) const override
	Helper method for stereo angle computation, DEPRECATED.
virtual double	length () const override
	Method to calculate length of a module.
virtual double	width () const override
	Method to calculate average width of a module.
virtual double	minWidth () const override
	Method to calculate minimum width of a module.
virtual double	maxWidth () const override
	Method to calculate maximum width of a module.
virtual double	etaPitch () const override
virtual bool	swapHitPhiReadoutDirection () const override
	Return true if hit local direction is the same as readout direction.
virtual bool	swapHitEtaReadoutDirection () const override
virtual bool	nearBondGap (const SiLocalPosition &, double) const override
	Test if near bond gap within tolerances, only relevant for SCT.
virtual DetectorShape	shape () const override
	Shape of element.
virtual double	sinStripAngleReco (double phiCoord, double etaCoord) const override
	Give strip angle in the reco frame.
double	pitch (const SiCellId &cellId) const
double	stripLength (const SiCellId &cellId) const
double	minR () const
double	maxR () const
double	phiWidth () const
double	stereo () const
virtual double	deadAreaUpperBoundary () const override
	DEPRECATED for StripStereoAnnulusDesign; no dead area.
virtual double	deadAreaLowerBoundary () const override
	give lower boundary of dead area
virtual double	deadAreaLength () const override
	give length of dead area
virtual int	numberOfConnectedCells (const SiReadoutCellId &readoutId) const override
	readout id -> id of connected diodes Not particularly useful for the SCT but implemented to keep a uniform interface.
virtual SiCellId	connectedCell (const SiReadoutCellId &readoutId, int number) const override
	readout id -> id of connected diodes.
virtual SiCellId	gangedCell (const SiCellId &cellId) const override
	If cell is ganged return the other cell, otherwise return an invalid id.
virtual SiReadoutCellId	readoutIdOfPosition (const SiLocalPosition &localPos) const override
	position -> id
int	crystals () const
	number of crystals within module side:
int	diodes () const
	number of strips within crystal:
int	cells () const
	number of readout stips within module side:
int	shift () const
	number of edge strips before first readout strip.
void	setMother (SCT_ModuleSideDesign *mother)
const SCT_ModuleSideDesign *	getMother () const
const std::map< int, const SCT_ModuleSideDesign * > &	getChildren () const
void	addChildDesign (int index, const SCT_ModuleSideDesign *element)
virtual SiIntersect	inDetector (const SiLocalPosition &localPosition, double phiTol, double etaTol) const override
	Test if point is in the active part of the detector with specified tolerances.
void	setSymmetry (bool phiSymmetric, bool etaSymmetric, bool depthSymmetric)
	Override default symmetries to prevent swapping of axes.
Axis	etaAxis () const
	local axis corresponding to eta direction:
Axis	phiAxis () const
	local axis corresponding to phi direction:
Axis	depthAxis () const
	local axis corresponding to depth direction:
double	thickness () const
	Method which returns thickness of the silicon wafer.
InDetDD::CarrierType	carrierType () const
	Return carrier type (ie electrons or holes)
bool	phiSymmetric () const
bool	etaSymmetric () const
bool	depthSymmetric () const
int	readoutSide () const
	ReadoutSide.
virtual const Amg::Transform3D	moduleShift () const

Protected Attributes
SCT_ReadoutScheme	m_scheme
InDetDD::DetectorType	m_detectorType {Undefined}

Private Attributes
const int	m_nRows
const std::vector< int >	m_nStrips
std::vector< int >	m_firstStrip
const std::vector< double >	m_pitch
const std::vector< double >	m_stripStartRadius
const std::vector< double >	m_stripEndRadius
const double	m_stereo
const double	m_R
const double	m_waferCentreR
const double	m_lengthBF
std::unique_ptr< Trk::SurfaceBounds >	m_bounds
const double	m_sinStereo
const double	m_cosStereo
const double	m_sinNegStereo
const double	m_cosNegStereo
const bool	m_usePC
bool	m_swapStripReadout
const SCT_ModuleSideDesign *	m_motherDesign {nullptr}
std::map< int, const SCT_ModuleSideDesign * >	m_childDesigns
Axis	m_etaAxis
Axis	m_phiAxis
Axis	m_depthAxis
double	m_thickness
InDetDD::CarrierType	m_carrierType
bool	m_phiSymmetric
bool	m_etaSymmetric
bool	m_depthSymmetric
bool	m_readoutSidePosDepth

Detailed Description

Definition at line 50 of file StripStereoAnnulusDesign.h.

Member Enumeration Documentation

◆ Axis

enum InDetDD::DetectorDesign::Axis

inherited

Enumerator
xAxis
yAxis
zAxis

Definition at line 59 of file DetectorDesign.h.

              {
        xAxis=0, yAxis, zAxis
    };

Constructor & Destructor Documentation

◆ StripStereoAnnulusDesign() [1/3]

InDetDD::StripStereoAnnulusDesign::StripStereoAnnulusDesign	(	const SiDetectorDesign::Axis &	stripDirection,
		const SiDetectorDesign::Axis &	thicknessDirection,
		const double &	thickness,
		const int &	readoutSide,
		const InDetDD::CarrierType &	carrier,
		const int &	nRows,
		const std::vector< int > &	nStrips,
		const std::vector< double > &	pitch,
		const std::vector< double > &	stripStart,
		const std::vector< double > &	stripEnd,
		const double &	stereoAngle,
		const double &	centreR,
		const bool &	usePC,
		InDetDD::DetectorType	detectorType = InDetDD::Undefined )

Definition at line 109 of file StripStereoAnnulusDesign.cxx.

                                                                :
    StripStereoAnnulusDesign(stripDirection,thicknessDirection,thickness,readoutSide,carrier,nRows,nStrips,
                             pitch,stripStartRadius,stripEndRadius,stereoAngle,centreR,centreR,usePC,detectorType){
//assuming here that centreR==waferCentreR
}

◆ StripStereoAnnulusDesign() [2/3]

InDetDD::StripStereoAnnulusDesign::StripStereoAnnulusDesign	(	const SiDetectorDesign::Axis &	stripDirection,
		const SiDetectorDesign::Axis &	thicknessDirection,
		const double &	thickness,
		const int &	readoutSide,
		const InDetDD::CarrierType &	carrier,
		const int &	nRows,
		const std::vector< int > &	nStrips,
		const std::vector< double > &	pitch,
		const std::vector< double > &	stripStart,
		const std::vector< double > &	stripEnd,
		const double &	stereoAngle,
		const double &	centreR,
		const double &	waferCentreR,
		const bool &	usePC,
		InDetDD::DetectorType	detectorType = InDetDD::Undefined )

Definition at line 15 of file StripStereoAnnulusDesign.cxx.

                                                                 :
    SCT_ModuleSideDesign(thickness, false, false, true, 1, 0, 0, 0, false, carrier,
                         readoutSide, stripDirection, thicknessDirection),
  m_nRows(nRows),
  m_nStrips(nStrips),
  m_pitch(pitch),
  m_stripStartRadius(stripStartRadius),
  m_stripEndRadius(stripEndRadius),
  m_stereo(stereoAngle),
  m_R(centreR),
  m_waferCentreR(waferCentreR),//if not specified in constructor, wafer centre assumed to simply be element centre 
  m_lengthBF(2. * waferCentreR * std::sin(stereoAngle*0.5)),// Eq. 5 p. 7
  m_sinStereo(std::sin(m_stereo)),
  m_cosStereo(std::cos(m_stereo)),
  m_sinNegStereo(-m_sinStereo),
  m_cosNegStereo(m_cosStereo),
  m_usePC(usePC)
{
    
    m_detectorType = detectorType;
 
    if (nRows < 0) {
        throw std::runtime_error(
                  "ERROR: StripStereoAnnulusDesign called with negative number of rows");
    }
 
    if (pitch.size() != (unsigned) nRows) {
        throw std::runtime_error(
                  "ERROR: StripStereoAnnulusDesign called with insufficiant pitch values for no. of rows");
    }
 
    if (stripStartRadius.size() != (unsigned) nRows) {
        throw std::runtime_error(
                  "ERROR: StripStereoAnnulusDesign called with insufficiant strip start-positions for no. of rows");
    }
 
    if (stripEndRadius.size() != (unsigned) nRows) {
        throw std::runtime_error(
                  "ERROR: StripStereoAnnulusDesign called with insufficiant strip end-positions for no. of rows");
    }
 
    int startStrip = 0;
    for (int r = 0; r < nRows; ++r) {
        m_firstStrip.push_back(startStrip);
        startStrip += m_nStrips[r];
    }
    int totalStrips = startStrip;
    m_firstStrip.push_back(totalStrips); // Add one more as total number of strips, so you can safely use row + 1 in a loop
 
    m_scheme.setCells(totalStrips);
    m_scheme.setDiodes(totalStrips);
 
 
    double phi = m_nStrips[0] * m_pitch[0];
 
    if (m_usePC) {
      // Maths for calculating PC bounds is based off
      // Trk::AnnulusBoundsPC::fromCartesian() and detailed at
      // https://hep.ph.liv.ac.uk/~jsmith/dropbox/AnnulusBoundsPC_Constructor_Maths.pdf.
      // This will later form some kind of INT Note or report or similar
      // internal documentation,
      Amg::Vector2D origin(m_R * (1.0 - std::cos(m_stereo)),
                           m_R * std::sin(-m_stereo));
      m_bounds = std::make_unique<Trk::AnnulusBoundsPC>(
        Trk::AnnulusBoundsPC(m_stripStartRadius[0],
                             m_stripEndRadius.back(),
                             phi * -0.5,
                             phi * 0.5,
                             origin,
                             0));
    } else {
      m_bounds = std::make_unique<Trk::AnnulusBounds>(
        Trk::AnnulusBounds(m_stripStartRadius[0],
                           m_stripEndRadius.back(),
                           m_waferCentreR,
                           phi,
                           m_stereo));
    }
}

◆ ~StripStereoAnnulusDesign()

InDetDD::StripStereoAnnulusDesign::~StripStereoAnnulusDesign ( )

default

◆ StripStereoAnnulusDesign() [3/3]

InDetDD::StripStereoAnnulusDesign::StripStereoAnnulusDesign ( const StripStereoAnnulusDesign & design )

Member Function Documentation

◆ addChildDesign()

void InDetDD::SCT_ModuleSideDesign::addChildDesign	(	int	index,
		const SCT_ModuleSideDesign *	element )

inlineinherited

Definition at line 292 of file SCT_ModuleSideDesign.h.

                                                                                               {
   m_childDesigns.emplace(index,child);
 }

◆ beamToStrip()

SiLocalPosition InDetDD::StripStereoAnnulusDesign::beamToStrip ( const SiLocalPosition & pos ) const

Definition at line 128 of file StripStereoAnnulusDesign.cxx.

                                                                                      {
  if (m_usePC) return beamToStripPC(pos);
 
  const double x_beam = pos.xEta();
  const double y_beam = pos.xPhi();
 
  //    Transform to strip frame SF (eq. 36 in ver G, combined with eq. 2 since we are in beam frame)
  //
  //    x_strip       cos(-m_stereo)  -sin(-m_stereo)  x-R     R
  //   (       )  =  [                               ](   ) + ( )
  //    y_strip       sin(-m_stereo)   cos(-m_stereo)   y      0
  const double x_strip = m_cosNegStereo * (x_beam - m_waferCentreR) - m_sinNegStereo * y_beam + m_waferCentreR;
  const double y_strip = m_sinNegStereo * (x_beam - m_waferCentreR) + m_cosNegStereo * y_beam;
  
  return {x_strip, y_strip};
}

◆ beamToStripPC()

SiLocalPosition InDetDD::StripStereoAnnulusDesign::beamToStripPC ( const SiLocalPosition & pos ) const

Definition at line 145 of file StripStereoAnnulusDesign.cxx.

                                                                                        {
  const double phi_beam = pos.xPhi();
  const double rad_beam = pos.xEta(); 
 
  // Convert to cart and use cartesian transform
  const double x_beam = rad_beam * std::cos(phi_beam);
  const double y_beam = rad_beam * std::sin(phi_beam);
 
  //    Transform to strip frame SF (eq. 36 in ver G, combined with eq. 2 since we are in beam frame)
  //
  //    x_strip       cos(-m_stereo)  -sin(-m_stereo)  x-R     R
  //   (       )  =  [                               ](   ) + ( )
  //    y_strip       sin(-m_stereo)   cos(-m_stereo)   y      0
  const double x_strip = m_cosNegStereo * (x_beam - m_waferCentreR) - m_sinNegStereo * y_beam + m_waferCentreR;
  const double y_strip = m_sinNegStereo * (x_beam - m_waferCentreR) + m_cosNegStereo * y_beam;
 
  const double phi_strip = std::atan2(y_strip, x_strip);
  const double rad_strip = std::hypot(x_strip, y_strip);
 
  return {rad_strip, phi_strip};
}

◆ beamToStripPCpolar()

SiLocalPosition InDetDD::StripStereoAnnulusDesign::beamToStripPCpolar ( const SiLocalPosition & pos ) const

Version of StripStereoAnnulusDesign::beamToStripPC transform based exclusively in a polar system.

Not currently used. Requires debugging and profiling

Parameters

pos

Returns: SiLocalPosition

Definition at line 175 of file StripStereoAnnulusDesign.cxx.

                                                                                             {
  const double phi_beam = pos.xPhi();
  const double rad_beam = pos.xEta(); 
 
  // Exclusively-polar transform:
  // Trig which can be precalculated in the future
  const double stereo_2 = m_stereo*0.5;
  const double cosStereo_2 = std::cos(stereo_2);
  const double sin_plus_stereo_2 = std::sin(stereo_2);
 
  // calculation for strip-frame radius uses cosine law - see ATL-COM-ITK-2021-048
  const double rad_conv = 2.0*m_waferCentreR*std::abs(sin_plus_stereo_2);
  const double rad_strip = std::sqrt( std::pow(rad_beam,2.0) + std::pow(rad_conv,2.0) + 2.0*rad_beam*rad_conv*cosStereo_2);
  
  // calculation for strip-frame angle uses sine law - see ATL-COM-ITK-2021-048
  const double phi_strip = M_PI_2 - stereo_2 - std::asin((rad_beam*std::sin(M_PI + phi_beam - stereo_2))/rad_strip);
 
  return {rad_strip, phi_strip};
}

◆ bounds()

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

overridevirtual

Get a reference to the module bounds object.

Returns: const Trk::SurfaceBounds& The module bounds.

Implements InDetDD::DetectorDesign.

Definition at line 326 of file StripStereoAnnulusDesign.cxx.

                                                               {
    return *(m_bounds); // Equivalent but more explicit than *m_bounds - 
                              // gets the normal pointer from the unique then dereferences it.
}

◆ carrierType()

InDetDD::CarrierType InDetDD::DetectorDesign::carrierType ( ) const

inlineinherited

Return carrier type (ie electrons or holes)

Definition at line 275 of file DetectorDesign.h.

                                                            {
    return m_carrierType;
}

◆ cellIdInRange()

SiCellId InDetDD::StripStereoAnnulusDesign::cellIdInRange ( const SiCellId & cellId ) const

overridevirtual

DEPRECATED: only used in a stupid example (2014) Check if cell is in range.

Returns the original cellId if it is in range, otherwise it

Reimplemented from InDetDD::SCT_ModuleSideDesign.

Definition at line 532 of file StripStereoAnnulusDesign.cxx.

                                                                             {
 
    if (!cellId.isValid()) {
        return {}; // Invalid
    }
    int row = cellId.etaIndex();
    int strip = cellId.phiIndex();
    if (strip < 0 || row < 0 || row >= m_nRows || strip >= m_nStrips[row]) {
        return {}; // Invalid
    }
    return cellId;
}

◆ cellIdOfPosition()

SiCellId InDetDD::StripStereoAnnulusDesign::cellIdOfPosition ( const SiLocalPosition & localPos ) const

overridevirtual

position -> id

Implements InDetDD::DetectorDesign.

Definition at line 331 of file StripStereoAnnulusDesign.cxx.

                                                                                    {
//
//    Find the row
//
    double rad_beam = (m_usePC) ? pos.xEta() : pos.r();
    if (rad_beam < m_stripStartRadius[0] || rad_beam >= m_stripEndRadius.back()) { 
        return {}; // return an invalid id
    }
 
    std::vector<double>::const_iterator endPtr = upper_bound(m_stripStartRadius.begin(), m_stripStartRadius.end(), rad_beam);
    int row = distance(m_stripStartRadius.begin(), endPtr) - 1;
    // Following should never happen, check is done on r above
    if (row < 0 || row >= m_nRows) {
      return {}; // return an invalid id
    }
    
    //    Transform to strip frame SF (eq. 36 in ver G, combined with eq. 2 since we are in beam frame)
    SiLocalPosition pos_strip = beamToStrip(pos);
    double phi_strip = (m_usePC) ? pos_strip.xPhi() : std::atan2(pos_strip.xPhi(), pos_strip.xEta()); 
    int strip = std::floor(phi_strip / m_pitch[row]) + m_nStrips[row] *0.5;
    if (strip < 0) { // Outside
      return {}; // return an invalid id
    }
    if (strip >= m_nStrips[row]) { // Outside
      return {}; // return an invalid id
    }
    
    int strip1D = strip1Dim(strip, row);
    return {strip1D, 0};
}

◆ cells()

int InDetDD::SCT_ModuleSideDesign::cells ( ) const

inlineinherited

number of readout stips within module side:

Definition at line 228 of file SCT_ModuleSideDesign.h.

                                             {
    return m_scheme.cells();
}

◆ centreR()

double InDetDD::StripStereoAnnulusDesign::centreR ( ) const

inline

Definition at line 366 of file StripStereoAnnulusDesign.h.

                                                      {
    return m_R;
}

◆ connectedCell()

SiCellId InDetDD::SCT_ModuleSideDesign::connectedCell	(	const SiReadoutCellId &	readoutId,
		int	number ) const

inlineoverridevirtualinherited

readout id -> id of connected diodes.

parameter number = 0 will refer to the primary diode and in general the cell number will be the same. NB. SiCellId cellId = connectedCell(readoutId, 0) will in general be equivalent to SiCellId cellId = readoutId

Implements InDetDD::DetectorDesign.

Definition at line 241 of file SCT_ModuleSideDesign.h.

                                                                      {
    return m_scheme.connectedCell(readoutId, number);
}

◆ crystals()

int InDetDD::SCT_ModuleSideDesign::crystals ( ) const

inlineinherited

number of crystals within module side:

Definition at line 216 of file SCT_ModuleSideDesign.h.

                                                {
    return m_scheme.crystals();
}

◆ deadAreaLength()

double InDetDD::StripStereoAnnulusDesign::deadAreaLength ( ) const

inlineoverridevirtual

give length of dead area

Implements InDetDD::SCT_ModuleSideDesign.

Definition at line 383 of file StripStereoAnnulusDesign.h.

                                                             {
    return 0.;
}

◆ deadAreaLowerBoundary()

double InDetDD::StripStereoAnnulusDesign::deadAreaLowerBoundary ( ) const

inlineoverridevirtual

give lower boundary of dead area

Implements InDetDD::SCT_ModuleSideDesign.

Definition at line 379 of file StripStereoAnnulusDesign.h.

                                                                    {
    return 0.;
}

◆ deadAreaUpperBoundary()

double InDetDD::StripStereoAnnulusDesign::deadAreaUpperBoundary ( ) const

inlineoverridevirtual

DEPRECATED for StripStereoAnnulusDesign; no dead area.

Implements InDetDD::SCT_ModuleSideDesign.

Definition at line 375 of file StripStereoAnnulusDesign.h.

                                                                    {
    return 0.;
}

◆ depthAxis()

DetectorDesign::Axis InDetDD::DetectorDesign::depthAxis ( ) const

inlineinherited

local axis corresponding to depth direction:

Definition at line 267 of file DetectorDesign.h.

                                                          {
    return m_depthAxis;
}

◆ depthSymmetric()

bool InDetDD::DetectorDesign::depthSymmetric ( ) const

inlineinherited

Definition at line 287 of file DetectorDesign.h.

                                                 {
    return m_depthSymmetric;
}

◆ diodes()

int InDetDD::SCT_ModuleSideDesign::diodes ( ) const

inlineinherited

number of strips within crystal:

Definition at line 220 of file SCT_ModuleSideDesign.h.

                                              {
    return m_scheme.diodes();
}

◆ diodesInRow()

int InDetDD::StripStereoAnnulusDesign::diodesInRow ( const int row ) const

inlineoverridevirtual

Reimplemented from InDetDD::SCT_ModuleSideDesign.

Definition at line 251 of file StripStereoAnnulusDesign.h.

                                                                     {
    return m_nStrips[row];
}

◆ distanceToDetectorEdge()

void InDetDD::StripStereoAnnulusDesign::distanceToDetectorEdge	(	const SiLocalPosition &	localPosition,
		double &	etaDist,
		double &	phiDist ) const

overridevirtual

Returns distance to nearest detector active edge +ve = inside -ve = outside.

Implements InDetDD::DetectorDesign.

Definition at line 579 of file StripStereoAnnulusDesign.cxx.

                                                                                                                           {
  // For eta, we use the Strip frame. This is centred at the beamline, x along eta, y along phi, z along depth
  // Happens to coincide with SiLocalPosition; no transform needed.
  double rInner = m_stripStartRadius[0];
  double rOuter = m_stripEndRadius[m_nRows - 1];
 
  double rad_beam = (m_usePC) ? pos.xEta() : std::hypot(pos.xEta(), pos.xPhi());
 
  if (rad_beam < rInner)
    etaDist = rad_beam - rInner;
  else if (rad_beam > rOuter)
    etaDist = rOuter - rad_beam;
  else
    etaDist = std::min(rOuter - rad_beam, rad_beam - rInner);
 
  // For phi, we use the Strip frame. Transform to Strip-frame:
  SiLocalPosition pos_strip = beamToStrip(pos);
 
  // Put these into polar coordinates
  double rad_strip = (m_usePC) ? pos_strip.xEta() : std::hypot(pos_strip.xEta(), pos_strip.xPhi());
  double phi_strip = (m_usePC) ? pos_strip.xPhi() : std::atan2(pos_strip.xPhi(), pos_strip.xEta());
 
  double phiAngleMax = m_pitch[0] * m_nStrips[0]*0.5;
  double phiAngleMin = -phiAngleMax;
 
  if (phi_strip < phiAngleMin)
    phiDist = rad_strip * (phi_strip - phiAngleMin);
  else if (phi_strip > phiAngleMax)
    phiDist = rad_strip * (phiAngleMax - phi_strip);
  else
    phiDist = rad_strip * std::min(phiAngleMax - phi_strip, phi_strip - phiAngleMin);
}

◆ endsOfStrip()

std::pair< SiLocalPosition, SiLocalPosition > InDetDD::StripStereoAnnulusDesign::endsOfStrip ( const SiLocalPosition & position ) const

overridevirtual

Give end points of the strip that covers the given position.

Implements InDetDD::SCT_ModuleSideDesign.

Definition at line 474 of file StripStereoAnnulusDesign.cxx.

                                                                                                                {
 
    SiCellId cellId = cellIdOfPosition(pos);
 
    auto [strip, row] = getStripRow(cellId);
 
    SiLocalPosition innerEnd = stripPosAtR(strip, row, m_stripStartRadius[row]);
 
    SiLocalPosition outerEnd = stripPosAtR(strip, row, m_stripEndRadius[row]);
 
    return std::pair<SiLocalPosition, SiLocalPosition>(innerEnd, outerEnd);
}

◆ etaAxis()

DetectorDesign::Axis InDetDD::DetectorDesign::etaAxis ( ) const

inlineinherited

local axis corresponding to eta direction:

Definition at line 259 of file DetectorDesign.h.

                                                        {
    return m_etaAxis;
}

◆ etaPitch()

double InDetDD::StripStereoAnnulusDesign::etaPitch ( ) const

overridevirtual

Implements InDetDD::DetectorDesign.

Definition at line 569 of file StripStereoAnnulusDesign.cxx.

                                                {
// Return average strip length
    return length() / m_stripStartRadius.size();
}

◆ etaSymmetric()

bool InDetDD::DetectorDesign::etaSymmetric ( ) const

inlineinherited

Definition at line 283 of file DetectorDesign.h.

                                               {
    return m_etaSymmetric;
}

◆ gangedCell()

SiCellId InDetDD::SCT_ModuleSideDesign::gangedCell ( const SiCellId & cellId ) const

inlineoverridevirtualinherited

If cell is ganged return the other cell, otherwise return an invalid id.

Not particularly useful for the SCT but implemented to keep a uniform interface. For SCT always returns an invalid id.

Implements InDetDD::SiDetectorDesign.

Definition at line 246 of file SCT_ModuleSideDesign.h.

                                                                       {
    // No ganged cells in SCT so always return invalid id.
    return SiCellId();
}

◆ getChildren()

const std::map< int, const SCT_ModuleSideDesign * > & InDetDD::SCT_ModuleSideDesign::getChildren ( ) const

inlineinherited

Definition at line 287 of file SCT_ModuleSideDesign.h.

                                                                                                 {
   return m_childDesigns;
 }

◆ getMother()

const SCT_ModuleSideDesign * InDetDD::SCT_ModuleSideDesign::getMother ( ) const

inlineinherited

Definition at line 283 of file SCT_ModuleSideDesign.h.

                                                                           {
   return m_motherDesign;
 }

◆ getStripRow()

std::pair< int, int > InDetDD::StripStereoAnnulusDesign::getStripRow ( SiCellId cellId ) const

finalvirtual

Get the strip and row number of the cell.

Can be used as auto [strip, row] = getStripRow(cellId);

Parameters

cellId The SiCellId

Returns: std::pai<int,int> A pair of ints representing the strip ID and row ID

Reimplemented from InDetDD::SCT_ModuleSideDesign.

Definition at line 287 of file StripStereoAnnulusDesign.cxx.

                                                                            {
    int strip1D = cellId.phiIndex();
    int rowNum = row(strip1D);
    int strip2D = strip1D - m_firstStrip[rowNum];
    return {strip2D,rowNum};
}

◆ inActiveArea()

bool InDetDD::StripStereoAnnulusDesign::inActiveArea	(	const SiLocalPosition &	chargePos,
		bool	checkBondGap = true ) const

overridevirtual

check if the position is in active area

Implements InDetDD::SCT_ModuleSideDesign.

Definition at line 487 of file StripStereoAnnulusDesign.cxx.

                                                                                    {
 
    SiCellId id = cellIdOfPosition(pos);
    bool inside = id.isValid();
 
    return inside;
}

◆ inDetector()

SiIntersect InDetDD::SCT_ModuleSideDesign::inDetector	(	const SiLocalPosition &	localPosition,
		double	phiTol,
		double	etaTol ) const

overridevirtualinherited

Test if point is in the active part of the detector with specified tolerances.

Reimplemented from InDetDD::DetectorDesign.

Definition at line 131 of file SCT_ModuleSideDesign.cxx.

                                                                                                                       {
 
    double etaDist = 0;
    double phiDist = 0;
 
    distanceToDetectorEdge(localPosition, etaDist, phiDist);
 
    SiIntersect state;
 
    if (phiDist < -phiTol || etaDist < -etaTol) {
        state.setOut();
        return state;
    }
 
    if (phiDist > phiTol && etaDist > etaTol) {
        state.setIn();
        return state;
    }
 
    // Near boundary.
    state.setNearBoundary();
    return state;
}

◆ length()

double InDetDD::StripStereoAnnulusDesign::length ( ) const

overridevirtual

Method to calculate length of a module.

Implements InDetDD::DetectorDesign.

Definition at line 545 of file StripStereoAnnulusDesign.cxx.

                                              {
// Return the total length of all strips, i.e. the active area length.
    return m_stripEndRadius.back() - m_stripStartRadius[0];
}

◆ localPositionOfCell()

SiLocalPosition InDetDD::StripStereoAnnulusDesign::localPositionOfCell ( const SiCellId & cellId ) const

overridevirtual

id -> position

Implements InDetDD::SCT_ModuleSideDesign.

Definition at line 362 of file StripStereoAnnulusDesign.cxx.

                                                                                          {
    auto [strip, row] = getStripRow(cellId);
    double r = (m_stripEndRadius[row] + m_stripStartRadius[row])*0.5;
 
    // Modified StripPosAtR? 
    // get phi of strip in the strip system
    double phi_strip = (strip - m_nStrips[row] / 2. + 0.5) * m_pitch[row];
 
    double b = -2. * m_lengthBF * std::sin(m_stereo/2. + phi_strip);
    double c = m_lengthBF * m_lengthBF - r * r;
    // this is the radius in the strip system
    double rad_strip = (-b + std::sqrt(b * b - 4. * c))/2.;
    
    if (m_usePC) return stripToBeam({rad_strip,phi_strip}); 
    else {
      //else use cart
      double x_strip = rad_strip * std::cos(phi_strip);
      double y_strip = rad_strip * std::sin(phi_strip);
      //Return strip pos not beam pos?
      return stripToBeam({x_strip, y_strip});
    }
}

◆ localPositionOfCellPC()

SiLocalPosition InDetDD::StripStereoAnnulusDesign::localPositionOfCellPC ( const SiCellId & cellId ) const

This is for debugging only.

Call to explicitly use PC cells.

Deprecated: No longer required as localPositionOfCell() will automatically use the appropriate transform.

Parameters

cellId

Returns: SiLocalPosition

Definition at line 428 of file StripStereoAnnulusDesign.cxx.

                                                                                            {
  
    auto [strip, row] = getStripRow(cellId);
    // this is the radius in the module / radial system
    double r = (m_stripEndRadius[row] + m_stripStartRadius[row])*0.5;
 
    // get phi of strip in the strip system
    double phiPrime = (strip - m_nStrips[row]*0.5 + 0.5) * m_pitch[row];
 
    double b = -2. * m_lengthBF * std::sin(m_stereo*0.5 + phiPrime);
    double c = m_lengthBF * m_lengthBF - r * r;
    // this is the radius in the strip systems
    double rPrime = (-b + std::sqrt(b * b - 4. * c))*0.5;
 
    // flip this, since coordinate system is defined the other way round
    double phi = -1*phiPrime;
 
    // xEta => r, xPhi = phi
    return SiLocalPosition(rPrime, phi); 
}

◆ localPositionOfCluster()

SiLocalPosition InDetDD::StripStereoAnnulusDesign::localPositionOfCluster	(	const SiCellId &	cellId,
		int	clusterSize ) const

overridevirtual

Implements InDetDD::SCT_ModuleSideDesign.

Definition at line 401 of file StripStereoAnnulusDesign.cxx.

                                                                                                              {
//
//    Return the average centre-position of the first and last strips.
//
 
    SiLocalPosition startPos = localPositionOfCell(cellId); // Should automatically detect PC Usage
 
    if (clusterSize <= 1) {
        return startPos;
    }
 
    auto [strip, row] = getStripRow(cellId);
    int stripEnd = strip + clusterSize - 1;
    SiCellId endId = strip1Dim(stripEnd, row);
    SiLocalPosition endPos = localPositionOfCell(endId); // Should automatically detect PC Usage
 
    return (startPos + endPos)*0.5;
}

◆ localPositionOfClusterPC()

SiLocalPosition InDetDD::StripStereoAnnulusDesign::localPositionOfClusterPC	(	const SiCellId &	cellId,
		int	clusterSize ) const

This is for debugging only.

Call to explicitly use PC clusters.

Deprecated: No longer required as localPositionOfCluster() will automatically use the appropriate transform.

Parameters

cellId
clusterSize

Returns: SiLocalPosition

Definition at line 458 of file StripStereoAnnulusDesign.cxx.

                                                                                                                {
    SiLocalPosition startPos = localPositionOfCellPC(cellId);
 
    if (clusterSize <= 1) {
        return startPos;
    }
 
    auto [strip, row] = getStripRow(cellId);
    int stripEnd = strip + clusterSize - 1;
    SiCellId endId = strip1Dim(stripEnd, row);
    SiLocalPosition endPos = localPositionOfCellPC(endId);
 
    return (startPos + endPos)*0.5;
}

◆ maxR()

double InDetDD::StripStereoAnnulusDesign::maxR ( ) const

inline

Definition at line 326 of file StripStereoAnnulusDesign.h.

                                                   {
    return m_stripEndRadius.back();
}

◆ maxWidth()

double InDetDD::StripStereoAnnulusDesign::maxWidth ( ) const

overridevirtual

Method to calculate maximum width of a module.

Implements InDetDD::DetectorDesign.

Definition at line 565 of file StripStereoAnnulusDesign.cxx.

                                                {
    return 2. * tan((m_pitch.back() * m_nStrips.back()) * 0.5) * m_stripEndRadius.back();
}

◆ minR()

double InDetDD::StripStereoAnnulusDesign::minR ( ) const

inline

Definition at line 322 of file StripStereoAnnulusDesign.h.

                                                   {
    return m_stripStartRadius[0];
}

◆ minWidth()

double InDetDD::StripStereoAnnulusDesign::minWidth ( ) const

overridevirtual

Method to calculate minimum width of a module.

Implements InDetDD::DetectorDesign.

Definition at line 561 of file StripStereoAnnulusDesign.cxx.

                                                {
    return 2. * tan((m_pitch[0] * m_nStrips[0]) * 0.5 ) * m_stripStartRadius[0];
}

◆ moduleShift()

const Amg::Transform3D InDetDD::DetectorDesign::moduleShift ( ) const

virtualinherited

Reimplemented in InDetDD::StripBoxDesign.

Definition at line 130 of file DetectorDesign.cxx.

                                                         {
    return Amg::Transform3D::Identity();
  }

◆ nearBondGap()

bool InDetDD::StripStereoAnnulusDesign::nearBondGap	(	const SiLocalPosition &	localPosition,
		double	etaTol ) const

inlineoverridevirtual

Test if near bond gap within tolerances, only relevant for SCT.

Implements InDetDD::SiDetectorDesign.

Definition at line 309 of file StripStereoAnnulusDesign.h.

                                                                                       {
// No bond gap in strip modules
    return false;
}

◆ neighboursOfCell()

void InDetDD::StripStereoAnnulusDesign::neighboursOfCell	(	const SiCellId &	cellId,
		std::vector< SiCellId > &	neighbours ) const

overridevirtual

Get the neighbouring diodes of a given diode: Cell for which the neighbours must be found List of cells which are neighbours of the given one.

Reimplemented from InDetDD::SCT_ModuleSideDesign.

Definition at line 299 of file StripStereoAnnulusDesign.cxx.

                                                                                                             {
 
 
    neighbours.clear();
 
    if (!cellId.isValid()) {
        return;
    }
 
    auto [strip, row] = getStripRow(cellId);
    int stripM = strip - 1;
    int stripP = strip + 1;
 
    if (stripM >= m_firstStrip[row]) {
        neighbours.emplace_back(stripM);
    }
    if (stripP < m_firstStrip[row] + m_nStrips[row]) {
        neighbours.emplace_back(stripP);
    }
 
}

◆ numberOfConnectedCells()

int InDetDD::SCT_ModuleSideDesign::numberOfConnectedCells ( const SiReadoutCellId & readoutId ) const

inlineoverridevirtualinherited

readout id -> id of connected diodes Not particularly useful for the SCT but implemented to keep a uniform interface.

Implements InDetDD::DetectorDesign.

Definition at line 236 of file SCT_ModuleSideDesign.h.

      {
    return m_scheme.numberOfConnectedCells(readoutId);
}

◆ operator=()

StripStereoAnnulusDesign & InDetDD::StripStereoAnnulusDesign::operator= ( const StripStereoAnnulusDesign & design )

◆ parameters()

SiDiodesParameters InDetDD::StripStereoAnnulusDesign::parameters ( const SiCellId & cellId ) const

overridevirtual

Return strip width, centre, length etc. Hard to find if this is used or not.

Implements InDetDD::DetectorDesign.

Definition at line 514 of file StripStereoAnnulusDesign.cxx.

                                                                                {
    throw std::runtime_error("Call to StripStereoAnnulusDesign::parameters; not yet implemented");
}

◆ phiAxis()

DetectorDesign::Axis InDetDD::DetectorDesign::phiAxis ( ) const

inlineinherited

local axis corresponding to phi direction:

Definition at line 263 of file DetectorDesign.h.

                                                        {
    return m_phiAxis;
}

◆ phiMeasureSegment()

HepGeom::Vector3D< double > InDetDD::StripStereoAnnulusDesign::phiMeasureSegment ( const SiLocalPosition & position ) const

overridevirtual

Helper method for stereo angle computation, DEPRECATED.

Implements InDetDD::SiDetectorDesign.

Definition at line 574 of file StripStereoAnnulusDesign.cxx.

      {
    throw std::runtime_error("Call to phiMeasureSegment, DEPRECATED, not implemented.");
}

◆ phiPitch() [1/3]

double InDetDD::StripStereoAnnulusDesign::phiPitch ( ) const

inlineoverridevirtual

Pitch in phi direction.

Implements InDetDD::DetectorDesign.

Definition at line 284 of file StripStereoAnnulusDesign.h.

                                                       {
// Return pitch in mm for the row just below or including the centre
    const int middleRow = m_nRows*0.5;
    const int middleStrip = m_nStrips[middleRow]*0.5; 
    return phiPitch(SiCellId(middleStrip, middleRow));
}

◆ phiPitch() [2/3]

double InDetDD::StripStereoAnnulusDesign::phiPitch ( const SiCellId & cellId ) const

inline

Definition at line 278 of file StripStereoAnnulusDesign.h.

                                                                             {
// Return pitch in mm for centre of this strip.
    const int row = cellId.etaIndex();
    return m_pitch[row] * (m_stripStartRadius[row] + m_stripEndRadius[row])*0.5;
}

◆ phiPitch() [3/3]

double InDetDD::StripStereoAnnulusDesign::phiPitch ( const SiLocalPosition & localPosition ) const

inlineoverridevirtual

Pitch in phi direction.

Implements InDetDD::DetectorDesign.

Definition at line 270 of file StripStereoAnnulusDesign.h.

                                                                                 {
// Return pitch in mm for the strip at this position, at this point's distance along the strip.
    const SiCellId cellId = cellIdOfPosition(pos);
    const int row = cellId.etaIndex();
    const double radius = (m_usePC) ? pos.xEta() : std::hypot(pos.xEta(), pos.xPhi());
    return m_pitch[row] * radius;
}

◆ phiPitchPhi() [1/3]

double InDetDD::StripStereoAnnulusDesign::phiPitchPhi ( ) const

inline

Definition at line 303 of file StripStereoAnnulusDesign.h.

                                                          {
  const int middleRow = m_nRows *0.5;
  const int middleStrip = m_nStrips[middleRow] *0.5; 
  return phiPitchPhi(SiCellId(middleStrip, middleRow));
}

◆ phiPitchPhi() [2/3]

double InDetDD::StripStereoAnnulusDesign::phiPitchPhi ( const SiCellId & cellId ) const

inline

Definition at line 298 of file StripStereoAnnulusDesign.h.

                                                                                {
  const int row = cellId.etaIndex();
  return m_pitch[row];
}

◆ phiPitchPhi() [3/3]

double InDetDD::StripStereoAnnulusDesign::phiPitchPhi ( const SiLocalPosition & localPosition ) const

inline

Definition at line 291 of file StripStereoAnnulusDesign.h.

                                                                                    {
  // even though SiLocalPosition might not be in PC, the etaIndex should be correct
  const SiCellId cellId = cellIdOfPosition(pos);
  const int row = cellId.etaIndex();
  return m_pitch[row];
}

◆ phiSymmetric()

bool InDetDD::DetectorDesign::phiSymmetric ( ) const

inlineinherited

Definition at line 279 of file DetectorDesign.h.

                                               {
    return m_phiSymmetric;
}

◆ phiWidth()

double InDetDD::StripStereoAnnulusDesign::phiWidth ( ) const

inline

Definition at line 331 of file StripStereoAnnulusDesign.h.

                                                       {
    return m_nStrips[0] * m_pitch[0];
}

◆ pitch()

double InDetDD::StripStereoAnnulusDesign::pitch ( const SiCellId & cellId ) const

◆ positionFromStrip() [1/2]

SiLocalPosition InDetDD::StripStereoAnnulusDesign::positionFromStrip ( const int stripNumber ) const

overridevirtual

gives position of strip center ALTERNATIVE/PREFERED way is to use localPositionOfCell(const SiCellId & cellId) or rawLocalPositionOfCell method in SiDetectorElement.

DEPRECATED

Reimplemented from InDetDD::SCT_ModuleSideDesign.

Definition at line 525 of file StripStereoAnnulusDesign.cxx.

                                                                                       {
    return localPositionOfCell(SiCellId(stripNumber, 0));
}

◆ positionFromStrip() [2/2]

SiLocalPosition InDetDD::StripStereoAnnulusDesign::positionFromStrip ( const SiCellId & cellId ) const

Definition at line 519 of file StripStereoAnnulusDesign.cxx.

                                                                                        {
//    throw std::runtime_error("Deprecated positionFromStrip called.");
    return localPositionOfCell(cellId);
}

◆ readoutIdOfCell()

SiReadoutCellId InDetDD::StripStereoAnnulusDesign::readoutIdOfCell ( const SiCellId & cellId ) const

inlineoverridevirtual

diode id -> readout id

Reimplemented from InDetDD::SCT_ModuleSideDesign.

Definition at line 314 of file StripStereoAnnulusDesign.h.

                                                                                             {
    int strip = cellId.phiIndex();
    int row = cellId.etaIndex();
 
    return SiReadoutCellId(strip, row);
}

◆ readoutIdOfPosition()

SiReadoutCellId InDetDD::SCT_ModuleSideDesign::readoutIdOfPosition ( const SiLocalPosition & localPos ) const

inlineoverridevirtualinherited

position -> id

Implements InDetDD::DetectorDesign.

Definition at line 256 of file SCT_ModuleSideDesign.h.

                                           {
    // Get the cellId then convert it to a readoutId.
    // NEED to think about active area and bond gap check
    return readoutIdOfCell(cellIdOfPosition(localPos));
}

◆ readoutSide()

int InDetDD::DetectorDesign::readoutSide ( ) const

inlineinherited

ReadoutSide.

+1 = postive depth side, -1 = negative depth side.

Definition at line 291 of file DetectorDesign.h.

                                             {
    return (m_readoutSidePosDepth) ? +1 : -1;
}

◆ row()

int InDetDD::StripStereoAnnulusDesign::row ( int stripId1Dim ) const

inlineoverridevirtual

Reimplemented from InDetDD::SCT_ModuleSideDesign.

Definition at line 339 of file StripStereoAnnulusDesign.h.

                                                              {
 
  //This method is scheduled for deletion
    std::vector<int>::const_iterator endPtr = std::upper_bound(m_firstStrip.begin(), m_firstStrip.end(), stripId1Dim);
    int rowNum = std::distance(m_firstStrip.begin(), endPtr) - 1;
    if (rowNum < 0 || rowNum >= m_nRows) {
      const std::string errMsg=std::string("StripId1Dim index out of acceptable range ") + __FILE__+std::string(": ")+std::to_string(__LINE__);
        throw std::runtime_error(errMsg);
    }
    return rowNum;
}

◆ scaledDistanceToNearestDiode()

double InDetDD::StripStereoAnnulusDesign::scaledDistanceToNearestDiode ( const SiLocalPosition & chargePos ) const

overridevirtual

give distance to the nearest diode in units of pitch, from 0.0 to 0.5, this method should be fast as it is called for every surface charge in the SCT_SurfaceChargesGenerator

Implements InDetDD::SCT_ModuleSideDesign.

Definition at line 496 of file StripStereoAnnulusDesign.cxx.

                                                                                              {
 
  
  //    Get stripframe phi of pos
  SiLocalPosition pos_stripframe = beamToStrip(pos);
  double pos_phi_stripframe = (m_usePC) ? pos_stripframe.xPhi() : std::atan2(pos_stripframe.xPhi(),pos_stripframe.xEta());
 
  //    Get stripframe phi of strip
  SiCellId cellId = cellIdOfPosition(pos);
  SiLocalPosition stripPos_beamframe  = localPositionOfCell(cellId);
  SiLocalPosition stripPos_stripframe = beamToStrip(stripPos_beamframe);
  double strip_phi_stripframe = (m_usePC) ? stripPos_stripframe.xPhi() : std::atan2(stripPos_stripframe.xPhi(),stripPos_stripframe.xEta());
  
  auto [strip, row] = getStripRow(cellId);
  return std::abs(pos_phi_stripframe - strip_phi_stripframe) / m_pitch[row];
}

◆ sensorCenter()

Amg::Vector3D InDetDD::StripStereoAnnulusDesign::sensorCenter ( ) const

overridevirtual

Return the centre of a sensor in the local reference frame.

For most designs, this is the origin of the local reference frame. For StripStereoAnnulusDesign, this is not the case (coordinate origin is "on the beamline")

Reimplemented from InDetDD::DetectorDesign.

Definition at line 262 of file StripStereoAnnulusDesign.cxx.

                                                         {
  return Amg::Vector3D(m_R, 0., 0.);
}

◆ setMother()

void InDetDD::SCT_ModuleSideDesign::setMother ( SCT_ModuleSideDesign * mother )

inherited

Definition at line 109 of file SCT_ModuleSideDesign.cxx.

                                                                   {
    if(m_motherDesign){
      const std::string errMsg=std::string("SCT_ModuleSideDesign already has a mother set!");
    throw std::runtime_error(errMsg);
    } 
    m_motherDesign = mother;
 
  }

◆ setSymmetry()

void InDetDD::DetectorDesign::setSymmetry	(	bool	phiSymmetric,
		bool	etaSymmetric,
		bool	depthSymmetric )

inherited

Override default symmetries to prevent swapping of axes.

NB. Flags can be changed from true to false but not false to true.

Definition at line 65 of file DetectorDesign.cxx.

                                                      {
    // Flags can be changed from true to false but not false to true.
    if (m_phiSymmetric) {
        m_phiSymmetric = phiSymmetric;
    }
    else if (phiSymmetric) {
        std::cout <<
            "SiDetectorDesign: WARNING! Attempt to allow swapping of xPhi axis direction ignored."
                  << std::endl;
    }
 
    if (m_etaSymmetric) {
        m_etaSymmetric = etaSymmetric;
    }
    else if (etaSymmetric) {
        std::cout <<
            "SiDetectorDesign: WARNING! Attempt to allow swapping of xEta axis direction ignored."
                  << std::endl;
    }
 
    if (m_depthSymmetric) {
        m_depthSymmetric = depthSymmetric;
    }
    else if (depthSymmetric) {
        std::cout <<
            "SiDetectorDesign: WARNING! Attempt to allow swapping of xDepth axis direction ignored."
                  << std::endl;
    }
}

◆ shape()

DetectorShape InDetDD::StripStereoAnnulusDesign::shape ( ) const

overridevirtual

Shape of element.

Reimplemented from InDetDD::DetectorDesign.

Definition at line 612 of file StripStereoAnnulusDesign.cxx.

 {
   if (m_usePC) return InDetDD::PolarAnnulus;
   else         return InDetDD::Annulus;
 }

◆ shift()

int InDetDD::SCT_ModuleSideDesign::shift ( ) const

inlineinherited

number of edge strips before first readout strip.

Definition at line 232 of file SCT_ModuleSideDesign.h.

                                             {
    return m_scheme.shift();
}

◆ sinStripAngleReco()

double InDetDD::StripStereoAnnulusDesign::sinStripAngleReco	(	double	phiCoord,
		double	etaCoord ) const

overridevirtual

Give strip angle in the reco frame.

Reimplemented from InDetDD::SiDetectorDesign.

Definition at line 266 of file StripStereoAnnulusDesign.cxx.

                                                                                         {
//
//    Transform to strip frame SF (eq. 36 in ver G, combined with eq. 2 since we are in beam frame)
//
    SiLocalPosition pos_strip = beamToStrip({etaCoord, phiCoord});
    
    double phi_strip = (m_usePC) ? pos_strip.xPhi() : std::atan2(pos_strip.xPhi(), pos_strip.xEta());
 
    // The minus sign below is because this routine is called by tracking software, which swaps x and y, then measures angles from y 
    // to x
    return -std::sin(phi_strip + m_stereo);
}

◆ stereo()

double InDetDD::StripStereoAnnulusDesign::stereo ( ) const

inline

Definition at line 335 of file StripStereoAnnulusDesign.h.

                                                     {
    return m_stereo;
}

◆ strip()

int InDetDD::StripStereoAnnulusDesign::strip ( int stripId1Dim ) const

inlineoverridevirtual

Reimplemented from InDetDD::SCT_ModuleSideDesign.

Definition at line 351 of file StripStereoAnnulusDesign.h.

                                                                {
    int rowNum = row(stripId1Dim);
 
    int strip2D = stripId1Dim - m_firstStrip[rowNum];
    if (strip2D < 0 || strip2D >= m_firstStrip[rowNum + 1]) {
        std::cout << "str1D " << stripId1Dim << " gives strip " << strip2D << " which is outside range 0 - " << 
                     m_firstStrip[rowNum + 1] << "\n";
    }
    return strip2D;
}

◆ strip1Dim()

int InDetDD::StripStereoAnnulusDesign::strip1Dim	(	int	strip,
		int	row ) const

overridevirtual

only relevant for SCT.

Return strip1Dim(int strip, int row) if SCT; otherwise -1

Reimplemented from InDetDD::SCT_ModuleSideDesign.

Definition at line 294 of file StripStereoAnnulusDesign.cxx.

                                                                 {
 
    return m_firstStrip[row] + strip;
}

◆ stripLength()

double InDetDD::StripStereoAnnulusDesign::stripLength ( const SiCellId & cellId ) const

Definition at line 618 of file StripStereoAnnulusDesign.cxx.

{
  SiLocalPosition lpoc = localPositionOfCell(cellId);
  std::pair<SiLocalPosition, SiLocalPosition> end = endsOfStrip(lpoc);
  double dx = end.second.xEta() - end.first.xEta();
  double dy = end.second.xPhi() - end.first.xPhi();
  return std::sqrt(dx * dx + dy * dy);
}

◆ stripPitch() [1/3]

double InDetDD::StripStereoAnnulusDesign::stripPitch ( ) const

inlineoverridevirtual

give the strip pitch (For Forward returns pitch at center)

Implements InDetDD::SCT_ModuleSideDesign.

Definition at line 266 of file StripStereoAnnulusDesign.h.

                                                         { // Don't use; which strip?
    return phiPitch();
}

◆ stripPitch() [2/3]

double InDetDD::StripStereoAnnulusDesign::stripPitch ( const SiCellId & cellId ) const

inline

Definition at line 262 of file StripStereoAnnulusDesign.h.

                                                                               {
    return phiPitch(cellId);
}

◆ stripPitch() [3/3]

double InDetDD::StripStereoAnnulusDesign::stripPitch ( const SiLocalPosition & chargePos ) const

inlineoverridevirtual

give the strip pitch (dependence on position needed for forward)

Implements InDetDD::SCT_ModuleSideDesign.

Definition at line 258 of file StripStereoAnnulusDesign.h.

                                                                                   {
    return phiPitch(pos);
}

◆ stripPosAtR()

SiLocalPosition InDetDD::StripStereoAnnulusDesign::stripPosAtR	(	int	strip,
		int	row,
		double	r ) const

Definition at line 385 of file StripStereoAnnulusDesign.cxx.

                                                                                        {
 
    double phi_strip = (strip - m_nStrips[row]*0.5 + 0.5) * m_pitch[row];
 
    double b = -2. * m_lengthBF * std::sin(m_stereo*0.5 + phi_strip);
    double c = m_lengthBF * m_lengthBF - r * r;
    double rad_strip = (-b + std::sqrt(b * b - 4. * c))*0.5;
 
    if (m_usePC) return stripToBeam({rad_strip,phi_strip});
 
    double x_strip = rad_strip * std::cos(phi_strip);
    double y_strip = rad_strip * std::sin(phi_strip);
   
    return stripToBeam({x_strip, y_strip});
}

◆ stripToBeam()

SiLocalPosition InDetDD::StripStereoAnnulusDesign::stripToBeam ( const SiLocalPosition & pos ) const

Definition at line 195 of file StripStereoAnnulusDesign.cxx.

                                                                                      {
  if (m_usePC) return stripToBeamPC(pos);
 
  const double x_strip = pos.xEta();
  const double y_strip = pos.xPhi();
 
  //    Transform to beam frame (eq. 36 in ver G, combined with eq. 2 since we are in strip frame)
  //
  //    x_beam       cos(m_stereo)  -sin(m_stereo)  x-R     R
  //   (       )  =  [                               ](   ) + ( )
  //    y_beam       sin(m_stereo)   cos(m_stereo)   y      0
  const double x_beam = m_cosStereo * (x_strip - m_waferCentreR) - m_sinStereo * y_strip + m_waferCentreR;
  const double y_beam = m_sinStereo * (x_strip - m_waferCentreR) + m_cosStereo * y_strip;
  
  return {x_beam, y_beam};
}

◆ stripToBeamPC()

SiLocalPosition InDetDD::StripStereoAnnulusDesign::stripToBeamPC ( const SiLocalPosition & pos ) const

Definition at line 212 of file StripStereoAnnulusDesign.cxx.

                                                                                        {
  const double phi_strip = pos.xPhi();
  const double rad_strip = pos.xEta(); 
 
  // Convert to cart and use previous transform
  const double x_strip = rad_strip * std::cos(phi_strip);
  const double y_strip = rad_strip * std::sin(phi_strip);
 
  //    Transform to beam frame (eq. 36 in ver G, combined with eq. 2 since we are in strip frame)
  //
  //    x_beam        cos(m_stereo)  -sin(m_stereo)  x-R     R
  //   (       )  =  [                             ](   ) + ( )
  //    y_beam        sin(m_stereo)   cos(m_stereo)   y      0
  const double x_beam = m_cosStereo * (x_strip - m_waferCentreR) - m_sinStereo * y_strip + m_waferCentreR;
  const double y_beam = m_sinStereo * (x_strip - m_waferCentreR) + m_cosStereo * y_strip;
 
  const double phi_beam = std::atan2(y_beam, x_beam);
  const double rad_beam = std::hypot(x_beam, y_beam);
 
  return {rad_beam, phi_beam};
}

◆ stripToBeamPCpolar()

SiLocalPosition InDetDD::StripStereoAnnulusDesign::stripToBeamPCpolar ( const SiLocalPosition & pos ) const

Version of StripStereoAnnulusDesign::stripToBeamPC transform based exclusively in a polar system.

Not currently used. Requires debugging and profiling

Parameters

pos

Returns: SiLocalPosition

Definition at line 242 of file StripStereoAnnulusDesign.cxx.

                                                                                             {
  const double phi_strip = pos.xPhi();
  const double rad_strip = pos.xEta(); 
 
  // Exclusively-polar transform:
  // Trig which can be precalculated in the future
  const double stereo_2 = m_stereo*0.5;
  const double sin_plus_stereo_2 = std::sin(stereo_2);
 
  // calculation for beam-frame radius uses cosine law - see ATL-COM-ITK-2021-048
  const double rad_conv = 2.0*m_waferCentreR*std::abs(sin_plus_stereo_2);
  const double rad_beam = std::sqrt( std::pow(rad_strip,2.0) + std::pow(rad_conv,2.0) + 2.0*rad_strip*rad_conv*std::cos(M_PI_2 - stereo_2 - phi_strip));
  
  // calculation for beam-frame angle uses sine law - see ATL-COM-ITK-2021-048
  const double phi_beam = stereo_2 - M_PI_2 + std::asin((rad_strip*std::sin(M_PI_2 - stereo_2 - phi_strip))/rad_beam);
 
  return {rad_beam, phi_beam};
}

◆ swapHitEtaReadoutDirection()

bool InDetDD::StripStereoAnnulusDesign::swapHitEtaReadoutDirection ( ) const

inlineoverridevirtual

Reimplemented from InDetDD::SCT_ModuleSideDesign.

Definition at line 391 of file StripStereoAnnulusDesign.h.

                                                                       {
    return false;
 }

◆ swapHitPhiReadoutDirection()

bool InDetDD::StripStereoAnnulusDesign::swapHitPhiReadoutDirection ( ) const

inlineoverridevirtual

Return true if hit local direction is the same as readout direction.

Reimplemented from InDetDD::SCT_ModuleSideDesign.

Definition at line 387 of file StripStereoAnnulusDesign.h.

                                                                       {
    return false;
}

◆ thickness()

double InDetDD::DetectorDesign::thickness ( ) const

inlineinherited

Method which returns thickness of the silicon wafer.

Definition at line 271 of file DetectorDesign.h.

                                              {
    return m_thickness;
}

◆ type()

InDetDD::DetectorType InDetDD::StripStereoAnnulusDesign::type ( ) const

inlinefinaloverridevirtual

Type of element.

Reimplemented from InDetDD::DetectorDesign.

Definition at line 370 of file StripStereoAnnulusDesign.h.

                                                               {
    return m_detectorType;
}

◆ waferCentreR()

double InDetDD::StripStereoAnnulusDesign::waferCentreR ( ) const

inline

Definition at line 362 of file StripStereoAnnulusDesign.h.

                                                           {
    return m_waferCentreR;
}

◆ width()

double InDetDD::StripStereoAnnulusDesign::width ( ) const

overridevirtual

Method to calculate average width of a module.

Implements InDetDD::DetectorDesign.

Definition at line 550 of file StripStereoAnnulusDesign.cxx.

                                             {
// Return approximate width between the two central rows
float middleRow = m_stripStartRadius.size() * 0.5 - 1;
if (middleRow < 0) {
  //single-row version
  return 2. * tan((m_pitch[0] * m_nStrips[0])*0.5) * ((m_stripStartRadius[0] + m_stripEndRadius[0])*0.5);
  }
else return 2. * tan((m_pitch[middleRow] * m_nStrips[middleRow]) * 0.5) * m_stripEndRadius[middleRow];
 
}

Definition at line 197 of file SCT_ModuleSideDesign.h.

197{Undefined};

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

Public Types

Public Member Functions

Protected Attributes

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ Axis

Constructor & Destructor Documentation

◆ StripStereoAnnulusDesign() [1/3]

◆ StripStereoAnnulusDesign() [2/3]

◆ ~StripStereoAnnulusDesign()

◆ StripStereoAnnulusDesign() [3/3]

Member Function Documentation

◆ addChildDesign()

◆ beamToStrip()

◆ beamToStripPC()

◆ beamToStripPCpolar()

◆ bounds()

◆ carrierType()

◆ cellIdInRange()

◆ cellIdOfPosition()

◆ cells()

◆ centreR()

◆ connectedCell()

◆ crystals()

◆ deadAreaLength()

◆ deadAreaLowerBoundary()

◆ deadAreaUpperBoundary()

◆ depthAxis()

◆ depthSymmetric()

◆ diodes()

◆ diodesInRow()

◆ distanceToDetectorEdge()

◆ endsOfStrip()

◆ etaAxis()

◆ etaPitch()

◆ etaSymmetric()

◆ gangedCell()

◆ getChildren()

◆ getMother()

◆ getStripRow()

◆ inActiveArea()

◆ inDetector()

◆ length()

◆ localPositionOfCell()

◆ localPositionOfCellPC()

◆ localPositionOfCluster()

◆ localPositionOfClusterPC()

◆ maxR()

◆ maxWidth()

◆ minR()

◆ minWidth()

◆ moduleShift()

◆ nearBondGap()

◆ neighboursOfCell()

◆ numberOfConnectedCells()

◆ operator=()

◆ parameters()

◆ phiAxis()

◆ phiMeasureSegment()

◆ phiPitch() [1/3]

◆ phiPitch() [2/3]

◆ phiPitch() [3/3]

◆ phiPitchPhi() [1/3]

◆ phiPitchPhi() [2/3]

◆ phiPitchPhi() [3/3]

◆ phiSymmetric()

◆ phiWidth()

◆ pitch()

◆ positionFromStrip() [1/2]

◆ positionFromStrip() [2/2]

◆ readoutIdOfCell()

◆ readoutIdOfPosition()

◆ readoutSide()

◆ row()

◆ scaledDistanceToNearestDiode()

◆ sensorCenter()

◆ setMother()

◆ setSymmetry()

◆ shape()