---

RadialStripTopology Class Reference

A StripTopology in which the component strips subtend a constant angular width, and, if projected, intersect at a point. More...

#include <Geometry/CommonTopologies/interface/RadialStripTopology.h>

Inheritance diagram for RadialStripTopology:

List of all members.

Public Member Functions
float	angularWidth () const
	Angular width of a each strip.
float	centreToIntersection () const
	Distance from the intersection of the projections of the extreme edges of the two extreme strips to the symmetry centre of the plane of strips.
virtual int	channel (const LocalPoint &) const
	Channel number corresponding to a given LocalPoint.
float	detHeight () const
	Length of long symmetry axis of plane of strips.
virtual LocalError	localError (const MeasurementPoint &, const MeasurementError &) const
	LocalError for a given MeasurementPoint with known MeasurementError.
virtual LocalError	localError (float strip, float stripErr2) const
	LocalError for a pure strip measurement, where 'strip' is the (float) position (a 'phi' angle wrt y axis) and stripErr2 is the sigma-squared.
virtual float	localPitch (const LocalPoint &) const
	Pitch (strip width) at a given LocalPoint.
virtual LocalPoint	localPosition (const MeasurementPoint &) const
	LocalPoint for a given MeasurementPoint What's a MeasurementPoint? In analogy with that used with TrapezoidalStripTopology objects, a MeasurementPoint is a 2-dim object.
virtual LocalPoint	localPosition (float strip) const
	LocalPoint on x axis for given 'strip' 'strip' is a float in units of the strip (angular) width.
virtual float	localStripLength (const LocalPoint &) const
	Length of a strip passing through a given LocalPpoint.
virtual MeasurementError	measurementError (const LocalPoint &, const LocalError &) const
virtual MeasurementPoint	measurementPosition (const LocalPoint &) const
virtual int	nearestStrip (const LocalPoint &) const
	Nearest strip to given LocalPoint.
virtual int	nstrips () const
	Total number of strips.
float	originToIntersection () const
	(y) distance from intersection of the projections of the strips to the local coordinate origin.
float	phiOfOneEdge () const
	Convenience function to access azimuthal angle of extreme edge of first strip measured relative to long symmetry axis of the plane of strips.
virtual float	phiPitch (void) const
	Phi pitch of each strip (= angular width!).
virtual float	pitch () const
	BEWARE: calling pitch() throws an exception.
	RadialStripTopology (int ns, float aw, float dh, float r, int yAx=1, float yMid=0.)
	Constructor from:.
virtual float	strip (const LocalPoint &) const
	Strip in which a given LocalPoint lies.
virtual float	stripAngle (float strip) const
	Angle between strip and symmetry axis (=local y axis) for given strip.
virtual float	stripLength () const
	Height of detector (= length of long symmetry axis of the plane of strips).
float	xOfStrip (int strip, float y) const
	Local x where centre of strip intersects input local y 'strip' should be in range 1 to nstrips() .
int	yAxisOrientation () const
	y axis orientation, 1 means detector width increases with local y
float	yCentreOfStripPlane () const
	Offset in local y between midpoint of detector (strip plane) extent and local origin.
float	yDistanceToIntersection (float y) const
	Distance in local y from a hit to the point of intersection of projected strips.
float	yExtentOfStripPlane () const
	y extent of strip plane
virtual	~RadialStripTopology ()
	Destructor.
Private Attributes
float	theAngularWidth
float	theCentreToIntersection
float	theDetHeight
int	theNumberOfStrips
float	thePhiOfOneEdge
int	theYAxisOrientation
float	yCentre
Friends
std::ostream &	operator<< (std::ostream &, const RadialStripTopology &)

---

Detailed Description

A StripTopology in which the component strips subtend a constant angular width, and, if projected, intersect at a point.

Author:

Tim Cox

WARNING! Wherever 'float strip' is used the units of 'strip' are angular widths of each strip. The range is from 0.0 at the extreme edge of the 'first' strip at one edge of the detector, to nstrip*angular width at the other edge.
The centre of the first strip is at strip = 0.5
The centre of the last strip is at strip = 0.5 + (nstrip-1)
This is for consistency with CommonDet usage of 'float strip' (but where units are strip pitch rather than strip angular width.)

WARNING! If the mid-point along local y of the plane of strips does not correspond to the local coordinate origin, set the final ctor argument appropriately.

Definition at line 26 of file RadialStripTopology.h.

---

Constructor & Destructor Documentation

RadialStripTopology::RadialStripTopology	(	int	ns,
		float	aw,
		float	dh,
		float	r,
		int	yAx = 1,
		float	yMid = 0.
	)

Constructor from:.

Parameters:

	ns	number of strips
	aw	angular width of a strip
	dh	detector height (usually 2 x apothem of TrapezoidalPlaneBounds)
	r	radial distance from symmetry centre of detector to the point at which the outer edges of the two extreme strips (projected) intersect.
	yAx	orientation of local y axis: 1 means pointing from the smaller side of the module to the larger side (along apothem), and -1 means in the opposite direction, i.e. from the larger side along the apothem to the smaller side. Default value is 1.
	yMid	local y offset if mid-point of detector (strip plane) does not coincide with local origin. This decouples the extent of strip plane from the boundary of the detector in which the RST is embedded.

Definition at line 9 of file RadialStripTopology.cc.

References LogTrace, theAngularWidth, theNumberOfStrips, thePhiOfOneEdge, theYAxisOrientation, and yCentre.

                                                                                                 :
  theNumberOfStrips(ns), theAngularWidth(aw), 
  theDetHeight(dh), theCentreToIntersection(r),
  theYAxisOrientation(yAx), yCentre( yMid) {   
  // Angular offset of extreme edge of detector, so that angle is
  // zero for a strip lying along local y axis = long symmetry axis of plane of strips
  thePhiOfOneEdge = -(theNumberOfStrips/2.) * theAngularWidth * yAx;
  
  LogTrace("RadialStripTopology") << "RadialStripTopology: constructed with"
        << " strips = " << ns
        << " width = " << aw << " rad "
        << " det_height = " << dh
        << " ctoi = " << r 
        << " phi_edge = " << thePhiOfOneEdge << " rad "
        << " y_ax_ori = " << theYAxisOrientation
        << " y_det_centre = " << yCentre 
        << "\n";
}    

virtual RadialStripTopology::~RadialStripTopology

(

)

[inline, virtual]

Destructor.

Definition at line 48 of file RadialStripTopology.h.

{}

---

Member Function Documentation

float RadialStripTopology::angularWidth

(

)

const [inline]

Angular width of a each strip.

Definition at line 169 of file RadialStripTopology.h.

References theAngularWidth.

Referenced by localError(), localPitch(), OffsetRadialStripTopology::localPosition(), measurementError(), measurementPosition(), OffsetRadialStripTopology::OffsetRadialStripTopology(), phiPitch(), OffsetRadialStripTopology::strip(), strip(), OffsetRadialStripTopology::stripAngle(), and stripAngle().

{ return theAngularWidth;}

float RadialStripTopology::centreToIntersection

(

)

const [inline]

Distance from the intersection of the projections of the extreme edges of the two extreme strips to the symmetry centre of the plane of strips.

Definition at line 191 of file RadialStripTopology.h.

References theCentreToIntersection.

Referenced by localError(), and operator<<().

{ return theCentreToIntersection; }

int RadialStripTopology::channel

(

const LocalPoint &

lp

)

const [virtual]

Channel number corresponding to a given LocalPoint.

This is effectively an integer version of strip(), with range 0 to nstrips-1.
LocalPoints outside the detector strip plane will be considered as contributing to the edge channels 0 or nstrips-1.

Implements Topology.

Reimplemented in OffsetRadialStripTopology, CSCGangedStripTopology, and CSCUngangedStripTopology.

Definition at line 182 of file RadialStripTopology.cc.

References min, strip(), and theNumberOfStrips.

Referenced by CSCUngangedStripTopology::channel().

                                                       {
  return std::min( int( strip(lp) ), theNumberOfStrips-1 );
}

float RadialStripTopology::detHeight

(

)

const [inline]

Length of long symmetry axis of plane of strips.

Definition at line 179 of file RadialStripTopology.h.

References theDetHeight.

Referenced by localError(), OffsetRadialStripTopology::localPosition(), localPosition(), localStripLength(), measurementError(), measurementPosition(), and operator<<().

{ return theDetHeight;}

LocalError RadialStripTopology::localError	(	const MeasurementPoint &	mp,
		const MeasurementError &	merr
	)			const [virtual]

LocalError for a given MeasurementPoint with known MeasurementError.

This may be used in Kalman filtering and hence must allow possible correlations between the components.

Implements StripTopology.

Definition at line 83 of file RadialStripTopology.cc.

References funct::A, angularWidth(), c1, c2, centreToIntersection(), funct::D, detHeight(), L, phi, s1, s2, funct::sin(), funct::sqrt(), stripAngle(), MeasurementError::uu(), MeasurementError::uv(), MeasurementError::vv(), PV2DBase< T, PVType, FrameType >::x(), PV2DBase< T, PVType, FrameType >::y(), and yAxisOrientation().

                                                                    {
  // Here we need to allow the possibility of correlated errors, since
  // that may happen during Kalman filtering

  //  float phi = phiOfOneEdge() + yAxisOrientation() * mp.x() * angularWidth();
  float phi = stripAngle( mp.x() );

  //  float t = tan( phi );                        // tan(angle between strip and y)
  //float c2 = 1./(1. + t*t);                    // cos(angle)**2
  // float cs = t*c2;                             // sin(angle)*cos(angle); tan carries sign of sin!
  // float s2 = t*t * c2;                         // sin(angle)**2
  

  float s1 = sin(phi);
  float s2 = s1*s1;
  float c2 = 1.-s2;
  float c1 = std::sqrt(c2);
  float cs = c1*s1;

  float A  = angularWidth();
  float A2 = A * A;

  // D is distance from intersection of edges to hit on strip
  float D = (centreToIntersection() + yAxisOrientation() * mp.y() * detHeight()) /c1;
  float D2 = D * D;

  // L is length of strip across face of chamber
  float L = detHeight()/ c1;  
  float L2  = L*L; 

  // MeasurementError elements are already squared
  // but they're normalized to products of A and L 
  // (N.B. uses L=length of strip, and not D=distance to intersection)
  // Remember to ensure measurement error components are indeed normalized like this!

  float SA2 = merr.uu() * A2;
  float SD2 = merr.vv() * L2; // Note this norm uses stripLength**2
  float RHOSASR = merr.uv() * A * L;

  float sx2 = SA2 * D2 * c2  +  2. * RHOSASR * D * cs  +  SD2 * s2;
  float sy2 = SA2 * D2 * s2  -  2. * RHOSASR * D * cs  +  SD2 * c2;
  float rhosxsy = cs * ( SD2 - D2 * SA2 )  +  RHOSASR * D * ( c2 - s2 );

  return LocalError(sx2, rhosxsy, sy2);
}

LocalError RadialStripTopology::localError	(	float	strip,
		float	stripErr2
	)			const [virtual]

LocalError for a pure strip measurement, where 'strip' is the (float) position (a 'phi' angle wrt y axis) and stripErr2 is the sigma-squared.

Both quantities are expressed in units of theAngularWidth of a strip.

Implements StripTopology.

Definition at line 55 of file RadialStripTopology.cc.

References angularWidth(), c2, centreToIntersection(), detHeight(), s2, stripAngle(), t, and funct::tan().

                                                                  {
  // Consider measurement as strip (phi) position and mid-point of strip
  // Since 'strip' is in units of angular strip-widths, stripErr2 is
  // required to be in corresponding units.

  float t = std::tan( stripAngle( strip ) );        // tan(angle between strip and y)
  float c2 = 1./(1. + t*t);                    // cos(angle)**2
  float cs = t*c2;                             // sin(angle)*cos(angle); tan carries sign of sin!
  float s2 = t*t * c2;                         // sin(angle)**2

  float D2 = centreToIntersection()*centreToIntersection() / c2; // (mid pt of strip to intersection)**2
  float L2 = detHeight()*detHeight() / c2;   // length**2 of strip across detector
  float A2 = angularWidth()*angularWidth();

  // The error**2 we're assigning is L2/12 wherever we measure the position along the strip
  // from... we just know the measurement is somewhere ON the strip. 

  float SD2 = L2 / 12.;       // SD = Sigma-Distance-along-strip
  float SA2 = A2 * stripErr2; // SA = Sigma-Angle-of-strip

  float sx2 = c2 * D2 * SA2 + s2 * SD2;
  float sy2 = s2 * D2 * SA2 + c2 * SD2;
  float rhosxsy = cs * ( SD2 - D2 * SA2 );

  return LocalError(sx2, rhosxsy, sy2);
}

float RadialStripTopology::localPitch

(

const LocalPoint &

lp

)

const [virtual]

Pitch (strip width) at a given LocalPoint.

BEWARE: are you sure you really want to call this for a RadialStripTopology?

Implements StripTopology.

Definition at line 197 of file RadialStripTopology.cc.

References angularWidth(), funct::cos(), nstrips(), funct::sin(), strip(), stripAngle(), PV3DBase< T, PVType, FrameType >::y(), and yDistanceToIntersection().

Referenced by TrackerValidationVariables::fillHitQuantities(), and CSCLayerGeometry::stripPitch().

                                                          {
  // The local pitch is the local x width of the strip at the local (x,y)

  // Calculating it is a nightmare...deriving the expression below is left
  // as a exercise for the reader. 

  float fstrip = strip(lp); // position in strip units
  int istrip = static_cast<int>(fstrip) + 1; // which strip number
  if (istrip>nstrips() ) istrip = nstrips(); // enforce maximum
  float fangle = stripAngle(static_cast<float>(istrip) - 0.5); // angle of strip centre
  float localp = yDistanceToIntersection( lp.y() ) * sin(angularWidth()) /
    ( cos(fangle-0.5*angularWidth())*cos(fangle+0.5*angularWidth()) );
  return localp;
}

LocalPoint RadialStripTopology::localPosition

(

const MeasurementPoint &

mp

)

const [virtual]

LocalPoint for a given MeasurementPoint
What's a MeasurementPoint?
In analogy with that used with TrapezoidalStripTopology objects, a MeasurementPoint is a 2-dim object.

The first dimension measures the angular position wrt central line of symmetry of detector, in units of strip (angular) widths (range 0 to total angle subtended by a detector).
The second dimension measures the fractional position along the strip (range -0.5 to +0.5).
BEWARE! The components are not Cartesian.

Implements StripTopology.

Reimplemented in OffsetRadialStripTopology.

Definition at line 41 of file RadialStripTopology.cc.

References detHeight(), phi, stripAngle(), funct::tan(), PV2DBase< T, PVType, FrameType >::x(), x, PV2DBase< T, PVType, FrameType >::y(), y, yAxisOrientation(), yCentreOfStripPlane(), and yDistanceToIntersection().

                                                                   {
  float phi = stripAngle( mp.x() );
  // 2nd component of MP is fractional position along strip, range -0.5 to +0.5 in direction of y
  // so distance along strip (measured from mid-point of length of strip) is
  //     mp.y() * (length of strip). 
  // Thus distance along y direction, from mid-point of strip, is 
  //    mp.y() * (length of strip) * cos(phi)
  // But (length of strip) = theDetHeight/cos(phi), so
  float y =  mp.y() * detHeight() + yCentreOfStripPlane();
  float x = yAxisOrientation() * yDistanceToIntersection( y ) * tan ( phi );
  return LocalPoint( x, y );
}

LocalPoint RadialStripTopology::localPosition

(

float

strip

)

const [virtual]

LocalPoint on x axis for given 'strip' 'strip' is a float in units of the strip (angular) width.

Implements StripTopology.

Reimplemented in OffsetRadialStripTopology.

Definition at line 36 of file RadialStripTopology.cc.

References originToIntersection(), stripAngle(), funct::tan(), and yAxisOrientation().

Referenced by OffsetRadialStripTopology::localPosition().

                                                    {
  return LocalPoint( yAxisOrientation() * originToIntersection() * tan( stripAngle(strip) ), 0.0 );
}

float RadialStripTopology::localStripLength

(

const LocalPoint &

lp

)

const [virtual]

Length of a strip passing through a given LocalPpoint.

Implements StripTopology.

Definition at line 217 of file RadialStripTopology.cc.

References c2, detHeight(), funct::sqrt(), t, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and yDistanceToIntersection().

                                                                      {
  float yHitToInter = yDistanceToIntersection( lp.y() );
  // since we're dealing with magnitudes, sign is unimportant
  float t  = lp.x() / yHitToInter;    // tan(angle between strip and y)
  float c2 = 1./(1. + t*t);           // cos(angle)**2
  return detHeight() / sqrt(c2);
}

MeasurementError RadialStripTopology::measurementError	(	const LocalPoint &	lp,
		const LocalError &	lerr
	)			const [virtual]

Implements Topology.

Definition at line 152 of file RadialStripTopology.cc.

References funct::A, angularWidth(), c2, funct::D, detHeight(), L, s2, funct::sqrt(), t, PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), LocalError::xy(), PV3DBase< T, PVType, FrameType >::y(), yAxisOrientation(), yDistanceToIntersection(), and LocalError::yy().

Referenced by TrackerValidationVariables::fillHitQuantities().

                                {

  float yHitToInter = yDistanceToIntersection( lp.y() );
  float t  = yAxisOrientation() * lp.x() / yHitToInter; // tan(angle between strip and y) 
  float c2 = 1./(1. + t*t);  // cos(angle)**2
  float cs = t*c2;           // sin(angle)*cos(angle); tan carries sign of sin!
  float s2 = 1. - c2;        // sin(angle)**2

  float A  = angularWidth();

  // L is length of strip across face of chamber
  float L2 = detHeight()*detHeight() / c2;  
  float L  = std::sqrt(L2); 

  // D is distance from intersection of edges to hit on strip
  float D2 = lp.x()*lp.x() + yHitToInter*yHitToInter;
  float D = std::sqrt(D2);

  float LP = D * A;
  float LP2 = LP * LP;

  float SA2 = ( c2 * lerr.xx() - 2. * cs * lerr.xy() + s2 * lerr.yy() ) / LP2;
  float SD2 = ( s2 * lerr.xx() + 2. * cs * lerr.xy() + c2 * lerr.yy() ) / L2;
  float RHOSASR = ( cs * ( lerr.xx() - lerr.yy() ) + ( c2 - s2 ) * lerr.xy() ) / (LP*L);

  return MeasurementError(SA2, RHOSASR, SD2);
}

MeasurementPoint RadialStripTopology::measurementPosition

(

const LocalPoint &

lp

)

const [virtual]

Implements Topology.

Reimplemented in OffsetRadialStripTopology.

Definition at line 143 of file RadialStripTopology.cc.

References angularWidth(), detHeight(), phi, phiOfOneEdge(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), yAxisOrientation(), yCentreOfStripPlane(), and yDistanceToIntersection().

Referenced by TrackerValidationVariables::fillHitQuantities(), and OffsetRadialStripTopology::measurementPosition().

                                                                   {
  // Note that this phi is (pi/2 - conventional local phi)
  // This means use atan2(x,y) rather than more usual atan2(y,x)
  float phi = yAxisOrientation() * atan2( lp.x(), yDistanceToIntersection( lp.y() ) );
  return MeasurementPoint( yAxisOrientation()*( phi-phiOfOneEdge() )/angularWidth(),
                          (lp.y() - yCentreOfStripPlane())/detHeight() );
}

int RadialStripTopology::nearestStrip

(

const LocalPoint &

lp

)

const [virtual]

Nearest strip to given LocalPoint.

Definition at line 225 of file RadialStripTopology.cc.

References nstrips(), and strip().

Referenced by CSCLayerGeometry::nearestStrip().

{
  // xxxStripTopology::strip() is expected to have range 0. to 
  // float(no of strips), but be extra careful and enforce that

  float fstrip = this->strip(lp);
  int n_strips = this->nstrips();
  fstrip = ( fstrip>=0. ? fstrip :  0. ); // enforce minimum 0.
  int near = static_cast<int>( fstrip ) + 1; // first strip is 1
  near = ( near<=n_strips ? near : n_strips ); // enforce maximum at right edge
  return near;
}

virtual int RadialStripTopology::nstrips

(

)

const [inline, virtual]

Total number of strips.

Implements StripTopology.

Definition at line 130 of file RadialStripTopology.h.

References theNumberOfStrips.

Referenced by localPitch(), nearestStrip(), CSCLayerGeometry::numberOfStrips(), operator<<(), and OffsetRadialStripTopology::strip().

{ return theNumberOfStrips; }

float RadialStripTopology::originToIntersection

(

)

const [inline]

(y) distance from intersection of the projections of the strips to the local coordinate origin.

Same as centreToIntersection() if symmetry centre of strip plane coincides with local origin.

Definition at line 198 of file RadialStripTopology.h.

References theCentreToIntersection, and yCentre.

Referenced by CSCStripTopology::equationOfStrip(), TrackerValidationVariables::fillHitQuantities(), OffsetRadialStripTopology::localPosition(), localPosition(), OffsetRadialStripTopology::strip(), OffsetRadialStripTopology::toLocal(), OffsetRadialStripTopology::toPrime(), and yDistanceToIntersection().

{ return (theCentreToIntersection - yCentre); }

float RadialStripTopology::phiOfOneEdge

(

)

const [inline]

Convenience function to access azimuthal angle of extreme edge of first strip measured relative to long symmetry axis of the plane of strips.

WARNING! This angle is measured clockwise from the local y axis which means it is in the conventional azimuthal phi plane, but azimuth is of course measured from local x axis not y. The range of this angle is -(full angle)/2 to +(full angle)/2.
where (full angle) = nstrips() * angularWidth().

Definition at line 212 of file RadialStripTopology.h.

References thePhiOfOneEdge.

Referenced by OffsetRadialStripTopology::localPosition(), measurementPosition(), operator<<(), OffsetRadialStripTopology::strip(), strip(), OffsetRadialStripTopology::stripAngle(), and stripAngle().

{ return thePhiOfOneEdge; }

virtual float RadialStripTopology::phiPitch

(

void

)

const [inline, virtual]

Phi pitch of each strip (= angular width!).

Definition at line 174 of file RadialStripTopology.h.

References angularWidth().

Referenced by operator<<(), and CSCLayerGeometry::stripPhiPitch().

{ return angularWidth(); }

float RadialStripTopology::pitch

(

)

const [virtual]

BEWARE: calling pitch() throws an exception.

Pitch is conventional name for width of something, but this is not sensible for a RadialStripTopology since strip widths vary with local y. Use localPitch(.) instead.

Implements StripTopology.

Definition at line 187 of file RadialStripTopology.cc.

                                 { 
  // BEWARE: this originally returned the pitch at the local origin 
  // but Tracker prefers throwing an exception since the strip is 
  // not constant along local x axis for a RadialStripTopology.

  throw Genexception("RadialStripTopology::pitch() called - makes no sense, use localPitch(.) instead.");
  return 0.;
}

float RadialStripTopology::strip

(

const LocalPoint &

lp

)

const [virtual]

Strip in which a given LocalPoint lies.

This is a float which represents the fractional strip position within the detector.
Returns zero if the LocalPoint falls at the extreme low edge of the detector or BELOW, and float(nstrips) if it falls at the extreme high edge or ABOVE.

Implements StripTopology.

Reimplemented in OffsetRadialStripTopology.

Definition at line 131 of file RadialStripTopology.cc.

References angularWidth(), phi, phiOfOneEdge(), theNumberOfStrips, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), yAxisOrientation(), and yDistanceToIntersection().

Referenced by CSCGangedStripTopology::channel(), channel(), localPitch(), and nearestStrip().

                                                     {
  // Note that this phi is measured from y axis, so sign of angle <-> sign of x * yAxisOrientation
  // Use atan2(x,y) rather than more usual atan2(y,x)

  float phi = atan2( lp.x(), yDistanceToIntersection( lp.y() ) );
  float aStrip = ( phi - yAxisOrientation() * phiOfOneEdge() )/angularWidth();
  if (aStrip < 0. ) aStrip = 0.;
  else if (aStrip > theNumberOfStrips)  aStrip = theNumberOfStrips;
  return aStrip;
}

float RadialStripTopology::stripAngle

(

float

strip

)

const [virtual]

Angle between strip and symmetry axis (=local y axis) for given strip.

This is like a phi angle but measured clockwise from y axis rather than counter clockwise from x axis. Note that 'strip' is a float with a continuous range from 0 to float(nstrips) to cover the whole detector, and the centres of strips correspond to half-integer values 0.5, 1.5, ..., nstrips-0.5 whereas values 1, 2, ... nstrips correspond to the upper phi edges of the strips.

Implements StripTopology.

Reimplemented in OffsetRadialStripTopology.

Definition at line 213 of file RadialStripTopology.cc.

References angularWidth(), phiOfOneEdge(), and yAxisOrientation().

Referenced by localError(), localPitch(), localPosition(), and xOfStrip().

                                                 {
  return ( phiOfOneEdge() + yAxisOrientation() * strip * angularWidth() );
}

virtual float RadialStripTopology::stripLength

(

)

const [inline, virtual]

Height of detector (= length of long symmetry axis of the plane of strips).

Implements StripTopology.

Definition at line 135 of file RadialStripTopology.h.

References theDetHeight.

{ return theDetHeight; }

float RadialStripTopology::xOfStrip	(	int	strip,
		float	y
	)			const

Local x where centre of strip intersects input local y
'strip' should be in range 1 to nstrips()
.

Definition at line 29 of file RadialStripTopology.cc.

References stripAngle(), funct::tan(), yAxisOrientation(), and yDistanceToIntersection().

Referenced by CSCLayerGeometry::xOfStrip().

                                                      {
  // Expect input 'strip' to be in range 1 to nstrips()
  float tanPhi = std::tan( stripAngle(static_cast<float>(strip) - 0.5 ) );
  return yAxisOrientation()* yDistanceToIntersection( y ) * tanPhi;
}

int RadialStripTopology::yAxisOrientation

(

)

const [inline]

y axis orientation, 1 means detector width increases with local y

Definition at line 228 of file RadialStripTopology.h.

References theYAxisOrientation.

Referenced by localError(), localPosition(), measurementError(), measurementPosition(), operator<<(), strip(), stripAngle(), xOfStrip(), and yDistanceToIntersection().

{ return theYAxisOrientation; }

float RadialStripTopology::yCentreOfStripPlane

(

)

const [inline]

Offset in local y between midpoint of detector (strip plane) extent and local origin.

Definition at line 233 of file RadialStripTopology.h.

References yCentre.

Referenced by OffsetRadialStripTopology::localPosition(), localPosition(), measurementPosition(), operator<<(), and CSCStripTopology::yLimitsOfStripPlane().

{ return yCentre; }

float RadialStripTopology::yDistanceToIntersection

(

float

y

)

const

Distance in local y from a hit to the point of intersection of projected strips.

Definition at line 238 of file RadialStripTopology.cc.

References originToIntersection(), and yAxisOrientation().

Referenced by TrackerValidationVariables::fillHitQuantities(), localPitch(), localPosition(), localStripLength(), measurementError(), measurementPosition(), strip(), and xOfStrip().

                                                                  {
  return yAxisOrientation() * y + originToIntersection();
}

float RadialStripTopology::yExtentOfStripPlane

(

)

const [inline]

y extent of strip plane

Definition at line 184 of file RadialStripTopology.h.

References theDetHeight.

Referenced by CSCStripTopology::yLimitsOfStripPlane().

{ return theDetHeight; } // same as detHeight()

---

Friends And Related Function Documentation

std::ostream& operator<<	(	std::ostream &	os,
		const RadialStripTopology &	rst
	)			[friend]

Definition at line 242 of file RadialStripTopology.cc.

{
  os  << "RadialStripTopology " << std::endl
      << " " << std::endl
      << "number of strips          " << rst.nstrips() << std::endl
      << "centre to whereStripsMeet " << rst.centreToIntersection() << std::endl
      << "detector height in y      " << rst.detHeight() << std::endl
      << "angular width of strips   " << rst.phiPitch() << std::endl
      << "phi of one edge           " << rst.phiOfOneEdge() << std::endl
      << "y axis orientation        " << rst.yAxisOrientation() << std::endl
      << "y of centre of strip plane " << rst.yCentreOfStripPlane() << std::endl;
  return os;
}

---

Member Data Documentation

float RadialStripTopology::theAngularWidth [private]

Definition at line 245 of file RadialStripTopology.h.

Referenced by angularWidth(), and RadialStripTopology().

float RadialStripTopology::theCentreToIntersection [private]

Definition at line 247 of file RadialStripTopology.h.

Referenced by centreToIntersection(), and originToIntersection().

float RadialStripTopology::theDetHeight [private]

Definition at line 246 of file RadialStripTopology.h.

Referenced by detHeight(), stripLength(), and yExtentOfStripPlane().

int RadialStripTopology::theNumberOfStrips [private]

Definition at line 244 of file RadialStripTopology.h.

Referenced by channel(), nstrips(), RadialStripTopology(), and strip().

float RadialStripTopology::thePhiOfOneEdge [private]

Definition at line 248 of file RadialStripTopology.h.

Referenced by phiOfOneEdge(), and RadialStripTopology().

int RadialStripTopology::theYAxisOrientation [private]

Definition at line 249 of file RadialStripTopology.h.

Referenced by RadialStripTopology(), and yAxisOrientation().

float RadialStripTopology::yCentre [private]

Definition at line 250 of file RadialStripTopology.h.

Referenced by originToIntersection(), RadialStripTopology(), and yCentreOfStripPlane().

---

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

• Geometry/CommonTopologies/interface/RadialStripTopology.h
• Geometry/CommonTopologies/src/RadialStripTopology.cc

---