GEMStripTopology Class Reference

#include <GEMStripTopology.h>

Inheritance diagram for GEMStripTopology:

Public Member Functions
float	angularWidth () const
float	centreToIntersection () const
int	channel (const LocalPoint &) const override
float	detHeight () const
	GEMStripTopology (int ns, float aw, float dh, float r0)
	GEMStripTopology (int ns, float aw, float dh, float r0, float yAx)
LocalError	localError (const MeasurementPoint &, const MeasurementError &) const override
LocalError	localError (float strip, float stripErr2) const override
virtual LocalError	localError (float strip, float stripErr2) const=0
virtual LocalError	localError (float strip, float stripErr2, const Topology::LocalTrackPred &) const
	conversion taking also the angle from the predicted track state More...
float	localPitch (const LocalPoint &) const override
LocalPoint	localPosition (const MeasurementPoint &) const override
LocalPoint	localPosition (float strip) const override
virtual LocalPoint	localPosition (float strip) const=0
virtual LocalPoint	localPosition (float strip, const Topology::LocalTrackPred &) const
	conversion taking also the predicted track state More...
float	localStripLength (const LocalPoint &) const override
MeasurementError	measurementError (const LocalPoint &, const LocalError &) const override
MeasurementPoint	measurementPosition (const LocalPoint &) const override
int	nearestStrip (const LocalPoint &) const
int	nstrips () const override
float	originToIntersection () const
float	phiOfOneEdge () const
float	phiPitch (void) const
float	pitch () const override
float	radius () const
float	strip (const LocalPoint &) const override
float	stripAngle (float strip) const override
float	stripLength () const override
float	xOfStrip (int strip, float y) const
float	yAxisOrientation () const
float	yCentreOfStripPlane () const
float	yDistanceToIntersection (float y) const
float	yExtentOfStripPlane () const
	~GEMStripTopology () override
Public Member Functions inherited from StripTopology
virtual float	coveredStrips (const LocalPoint &lp1, const LocalPoint &lp2) const
virtual LocalError	localError (const MeasurementPoint &, const MeasurementError &) const=0
virtual LocalError	localError (const MeasurementPoint &mp, const MeasurementError &me, const LocalTrackPred &) const
	conversion taking also the angle from the predicted track state More...
virtual LocalError	localError (float strip, float stripErr2, const Topology::LocalTrackPred &) const
	conversion taking also the angle from the predicted track state More...
virtual float	localPitch (const LocalPoint &lp, const Topology::LocalTrackAngles &) const
	conversion taking also the angle from the track state (LocalTrajectoryParameters) More...
virtual LocalPoint	localPosition (const MeasurementPoint &) const=0
virtual LocalPoint	localPosition (const MeasurementPoint &mp, const LocalTrackPred &) const
	conversion taking also the angle from the predicted track state More...
virtual LocalPoint	localPosition (float strip, const Topology::LocalTrackPred &) const
	conversion taking also the predicted track state More...
virtual float	localStripLength (const LocalPoint &lp, const Topology::LocalTrackAngles &) const
	conversion taking also the angle from the track state (LocalTrajectoryParameters) More...
virtual float	strip (const LocalPoint &lp, const Topology::LocalTrackAngles &) const
	conversion taking also the angle from the track state (LocalTrajectoryParameters) More...
	~StripTopology () override
Public Member Functions inherited from Topology
virtual int	channel (const LocalPoint &lp, const LocalTrackAngles &) const
	conversion taking also the angle from the track state More...
virtual LocalError	localError (const MeasurementPoint &mp, const MeasurementError &me, const LocalTrackPred &) const
	conversion taking also the angle from the predicted track state More...
virtual LocalPoint	localPosition (const MeasurementPoint &mp, const LocalTrackPred &) const
	conversion taking also the angle from the predicted track state More...
virtual MeasurementError	measurementError (const LocalPoint &lp, const LocalError &le, const LocalTrackAngles &) const
	conversion taking also the angle from the track state More...
virtual MeasurementPoint	measurementPosition (const LocalPoint &lp, const LocalTrackAngles &) const
	conversion taking also the angle from the track state More...
virtual	~Topology ()

Private Attributes
float	angularWidth_
float	centreToIntersection_
float	detHeight_
int	numberOfStrips_
float	phiOfOneEdge_
float	yAxisOrientation_
float	yCentre_

Additional Inherited Members
Public Types inherited from Topology
typedef Point2DBase< double, LocalTag >	Local2DPoint
typedef Vector2D::MathVector	MathVector2D
typedef Basic2DVector< double >	Vector2D

Detailed Description

based on CSCRadialStripTopology and TrapezoidalStripTopology

Author

Hyunyong Kim - TAMU

Definition at line 11 of file GEMStripTopology.h.

Constructor & Destructor Documentation

GEMStripTopology() [1/2]

GEMStripTopology::GEMStripTopology	(	int	ns,
		float	aw,
		float	dh,
		float	r0
	)

GEMStripTopology based on CSCRadialStripTopology and TrapezoidalStripTopology

Author

Hyunyong Kim - TAMU

Definition at line 11 of file GEMStripTopology.cc.

    : numberOfStrips_(ns), angularWidth_(aw), detHeight_(dh), centreToIntersection_(r0) {
  assert(angularWidth_ != 0);
  assert(detHeight_ != 0);
  yAxisOrientation_ = 1;
  phiOfOneEdge_ = -(0.5 * numberOfStrips_) * angularWidth_ * yAxisOrientation_;
  yCentre_ = 0;
  LogTrace("GEMStripTopology") << "Constructing GEMStripTopology with"
                               << " nstrips = " << ns << " angular width = " << aw << " det. height = " << dh
                               << " r0 = " << r0 << "\n";
}

References angularWidth_, cms::cuda::assert(), detHeight_, cuy::dh, LogTrace, numberOfStrips_, phiOfOneEdge_, yAxisOrientation_, and yCentre_.

GEMStripTopology() [2/2]

GEMStripTopology::GEMStripTopology	(	int	ns,
		float	aw,
		float	dh,
		float	r0,
		float	yAx
	)

Definition at line 23 of file GEMStripTopology.cc.

    : numberOfStrips_(ns), angularWidth_(aw), detHeight_(dh), centreToIntersection_(r0), yAxisOrientation_(yAx) {
  assert(angularWidth_ != 0);
  assert(detHeight_ != 0);
  phiOfOneEdge_ = -(0.5 * numberOfStrips_) * angularWidth_ * yAxisOrientation_;
  yCentre_ = 0;
  LogTrace("GEMStripTopology") << "Constructing GEMStripTopology with"
                               << " nstrips = " << ns << " angular width = " << aw << " det. height = " << dh
                               << " r0 = " << r0 << " yAxOrientation = " << yAx << "\n";
}

References angularWidth_, cms::cuda::assert(), detHeight_, cuy::dh, LogTrace, numberOfStrips_, phiOfOneEdge_, yAxisOrientation_, and yCentre_.

~GEMStripTopology()

GEMStripTopology::~GEMStripTopology ( )

inlineoverride

Definition at line 15 of file GEMStripTopology.h.

{}

Member Function Documentation

angularWidth()

float GEMStripTopology::angularWidth ( ) const

inline

Definition at line 50 of file GEMStripTopology.h.

{ return angularWidth_; }

References angularWidth_.

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

centreToIntersection()

float GEMStripTopology::centreToIntersection ( ) const

inline

Definition at line 53 of file GEMStripTopology.h.

{ return centreToIntersection_; }

References centreToIntersection_.

Referenced by localError().

channel()

int GEMStripTopology::channel ( const LocalPoint &  lp ) const

overridevirtual

Implements Topology.

Definition at line 107 of file GEMStripTopology.cc.

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

References min(), numberOfStrips_, and strip().

Referenced by GEMSignalModel::simulateClustering().

detHeight()

float GEMStripTopology::detHeight ( ) const

inline

Definition at line 51 of file GEMStripTopology.h.

{ return detHeight_; }

References detHeight_.

Referenced by localError(), localPosition(), localStripLength(), measurementError(), and measurementPosition().

localError() [1/4]

LocalError GEMStripTopology::localError ( const MeasurementPoint &  mp, const MeasurementError &  me  ) const

overridevirtual

Implements Topology.

Definition at line 55 of file GEMStripTopology.cc.

                                                                                                    {
  const double phi(stripAngle(mp.x())), s1(std::sin(phi)), c1(std::cos(phi));
  assert(c1 != 0);
  const double cs(s1 * c1), s2(s1 * s1),
      c2(1 - s2),  // rotation matrix
 
      T(angularWidth() * (centreToIntersection() + yAxisOrientation() * mp.y() * detHeight()) /
        c1),                // tangential measurement unit (local pitch)
      R(detHeight() / c1),  // radial measurement unit (strip length)
      tt(me.uu() * T * T),  // tangential sigma^2
      rr(me.vv() * R * R),  // radial sigma^2
      tr(me.uv() * T * R),
 
      xx(c2 * tt + 2 * cs * tr + s2 * rr), yy(s2 * tt - 2 * cs * tr + c2 * rr), xy(cs * (rr - tt) + tr * (c2 - s2));
 
  return LocalError(xx, xy, yy);
}

References angularWidth(), cms::cuda::assert(), alignmentValidation::c1, centreToIntersection(), funct::cos(), fwrapper::cs, detHeight(), hlt_dqm_clientPB-live_cfg::me, phi, dttmaxenums::R, findQualityFiles::rr, indexGen::s2, funct::sin(), stripAngle(), groupFilesInBlocks::tt, PV2DBase< T, PVType, FrameType >::x(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, PV2DBase< T, PVType, FrameType >::y(), yAxisOrientation(), and geometryCSVtoXML::yy.

localError() [2/4]

LocalError GEMStripTopology::localError ( float  strip, float  stripErr2  ) const

overridevirtual

Implements StripTopology.

Definition at line 45 of file GEMStripTopology.cc.

                                                                          {
  const double phi(stripAngle(strip)), t1(std::tan(phi)), t2(t1 * t1),
      tt(stripErr2 * std::pow(centreToIntersection() * angularWidth(), 2)),  // tangential sigma^2   *c2
      rr(std::pow(detHeight(), 2) * (1.f / 12.f)),  // radial sigma^2( uniform prob density along strip)  *c2
 
      xx(tt + t2 * rr), yy(t2 * tt + rr), xy(t1 * (rr - tt));
 
  return LocalError(xx, xy, yy);
}

References angularWidth(), centreToIntersection(), detHeight(), f, phi, funct::pow(), findQualityFiles::rr, strip(), stripAngle(), RandomServiceHelper::t1, RandomServiceHelper::t2, funct::tan(), groupFilesInBlocks::tt, geometryCSVtoXML::xx, geometryCSVtoXML::xy, and geometryCSVtoXML::yy.

localError() [3/4]

virtual LocalError StripTopology::localError

localError() [4/4]

virtual LocalError StripTopology::localError

inline

conversion taking also the angle from the predicted track state

Definition at line 26 of file StripTopology.h.

                                                                                                     {
    return localError(strip, stripErr2);
  }

localPitch()

float GEMStripTopology::localPitch ( const LocalPoint &  lp ) const

overridevirtual

Implements StripTopology.

Definition at line 111 of file GEMStripTopology.cc.

                                                             {
  const int istrip = std::min(nstrips(), static_cast<int>(strip(lp)) + 1);  // which strip number
  const float fangle = stripAngle(static_cast<float>(istrip) - 0.5);        // angle of strip centre
  assert(std::cos(fangle - 0.5f * angularWidth()) != 0);
  return yDistanceToIntersection(lp.y()) * std::sin(angularWidth()) /
         std::pow(std::cos(fangle - 0.5f * angularWidth()), 2);
}

References angularWidth(), cms::cuda::assert(), funct::cos(), f, min(), nstrips(), funct::pow(), funct::sin(), strip(), stripAngle(), PV3DBase< T, PVType, FrameType >::y(), and yDistanceToIntersection().

Referenced by pitch().

localPosition() [1/4]

LocalPoint GEMStripTopology::localPosition ( const MeasurementPoint &  mp ) const

overridevirtual

Implements Topology.

Definition at line 38 of file GEMStripTopology.cc.

                                                                           {
  const float  // y = (L/cos(phi))*mp.y()*cos(phi)
      y(mp.y() * detHeight() + yCentreOfStripPlane()),
      x(yAxisOrientation() * yDistanceToIntersection(y) * std::tan(stripAngle(mp.x())));
  return LocalPoint(x, y);
}

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

localPosition() [2/4]

LocalPoint GEMStripTopology::localPosition ( float  strip ) const

overridevirtual

Implements StripTopology.

Definition at line 34 of file GEMStripTopology.cc.

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

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

localPosition() [3/4]

virtual LocalPoint StripTopology::localPosition

localPosition() [4/4]

virtual LocalPoint StripTopology::localPosition

inline

conversion taking also the predicted track state

Definition at line 20 of file StripTopology.h.

                                                                                       {
    return localPosition(strip);
  }

localStripLength()

float GEMStripTopology::localStripLength ( const LocalPoint &  lp ) const

overridevirtual

Implements StripTopology.

Definition at line 123 of file GEMStripTopology.cc.

                                                                   {
  assert(yDistanceToIntersection(lp.y()) != 0);
  return detHeight() * std::sqrt(1.f + std::pow(lp.x() / yDistanceToIntersection(lp.y()), 2));
}

References cms::cuda::assert(), detHeight(), f, funct::pow(), mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and yDistanceToIntersection().

measurementError()

MeasurementError GEMStripTopology::measurementError ( const LocalPoint &  p, const LocalError &  e  ) const

overridevirtual

Implements Topology.

Definition at line 91 of file GEMStripTopology.cc.

                                                                                                  {
  const double yHitToInter(yDistanceToIntersection(p.y()));
  assert(yHitToInter != 0);
  const double t(yAxisOrientation() * p.x() / yHitToInter),  // tan(strip angle)
      cs(t / (1 + t * t)), s2(t * cs), c2(1 - s2),           // rotation matrix
 
      T2(1. / (std::pow(angularWidth(), 2) *
               (std::pow(p.x(), 2) + std::pow(yHitToInter, 2)))),  // 1./tangential measurement unit (local pitch) ^2
      R2(c2 / std::pow(detHeight(), 2)),                           // 1./ radial measurement unit (strip length) ^2
 
      uu((c2 * e.xx() - 2 * cs * e.xy() + s2 * e.yy()) * T2), vv((s2 * e.xx() + 2 * cs * e.xy() + c2 * e.yy()) * R2),
      uv((cs * (e.xx() - e.yy()) + e.xy() * (c2 - s2)) * std::sqrt(T2 * R2));
 
  return MeasurementError(uu, uv, vv);
}

References angularWidth(), cms::cuda::assert(), fwrapper::cs, detHeight(), MillePedeFileConverter_cfg::e, AlCaHLTBitMon_ParallelJobs::p, funct::pow(), indexGen::s2, mathSSE::sqrt(), submitPVValidationJobs::t, yAxisOrientation(), and yDistanceToIntersection().

measurementPosition()

MeasurementPoint GEMStripTopology::measurementPosition ( const LocalPoint &  lp ) const

overridevirtual

Implements Topology.

Definition at line 84 of file GEMStripTopology.cc.

                                                                                 {
  const float  // phi is [pi/2 - conventional local phi], use atan2(x,y) rather than atan2(y,x)
      phi(yAxisOrientation() * std::atan2(lp.x(), yDistanceToIntersection(lp.y())));
  return MeasurementPoint(yAxisOrientation() * (phi - phiOfOneEdge()) / angularWidth(),
                          (lp.y() - yCentreOfStripPlane()) / detHeight());
}

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

nearestStrip()

int GEMStripTopology::nearestStrip

(

const LocalPoint &

lp

)

const

Definition at line 80 of file GEMStripTopology.cc.

                                                             {
  return std::min(nstrips(), static_cast<int>(std::max(float(0), strip(lp))) + 1);
}

References SiStripPI::max, min(), nstrips(), and strip().

nstrips()

int GEMStripTopology::nstrips ( ) const

inlineoverridevirtual

Implements StripTopology.

Definition at line 44 of file GEMStripTopology.h.

{ return numberOfStrips_; }

References numberOfStrips_.

Referenced by localPitch(), nearestStrip(), and strip().

originToIntersection()

float GEMStripTopology::originToIntersection ( ) const

inline

Definition at line 55 of file GEMStripTopology.h.

{ return (centreToIntersection_ - yCentre_); }

References centreToIntersection_, and yCentre_.

Referenced by localPosition(), and yDistanceToIntersection().

phiOfOneEdge()

float GEMStripTopology::phiOfOneEdge ( ) const

inline

Definition at line 60 of file GEMStripTopology.h.

{ return phiOfOneEdge_; }

References phiOfOneEdge_.

Referenced by measurementPosition(), strip(), and stripAngle().

phiPitch()

float GEMStripTopology::phiPitch ( void  ) const

inline

Definition at line 36 of file GEMStripTopology.h.

{ return angularWidth(); }

References angularWidth().

pitch()

float GEMStripTopology::pitch ( ) const

overridevirtual

Implements StripTopology.

Definition at line 109 of file GEMStripTopology.cc.

{ return localPitch(LocalPoint(0, 0)); }

References localPitch().

radius()

float GEMStripTopology::radius ( void  ) const

inline

Definition at line 54 of file GEMStripTopology.h.

{ return centreToIntersection_; }

References centreToIntersection_.

strip()

float GEMStripTopology::strip ( const LocalPoint &  lp ) const

overridevirtual

Implements StripTopology.

Definition at line 73 of file GEMStripTopology.cc.

                                                        {
  const float  // phi is measured from y axis --> sign of angle is sign of x * yAxisOrientation --> use atan2(x,y), not atan2(y,x)
      phi(std::atan2(lp.x(), yDistanceToIntersection(lp.y()))),
      aStrip((phi - yAxisOrientation() * phiOfOneEdge()) / angularWidth());
  return std::max(float(0), std::min((float)nstrips(), aStrip));
}

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

Referenced by channel(), localError(), localPitch(), localPosition(), nearestStrip(), stripAngle(), and xOfStrip().

stripAngle()

float GEMStripTopology::stripAngle ( float  strip ) const

overridevirtual

Implements StripTopology.

Definition at line 119 of file GEMStripTopology.cc.

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

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

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

stripLength()

float GEMStripTopology::stripLength ( ) const

inlineoverridevirtual

Implements StripTopology.

Definition at line 46 of file GEMStripTopology.h.

{ return detHeight_; }

References detHeight_.

xOfStrip()

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

Definition at line 132 of file GEMStripTopology.cc.

                                                         {
  return yAxisOrientation() * yDistanceToIntersection(y) * std::tan(stripAngle(static_cast<float>(strip) - 0.5));
}

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

yAxisOrientation()

float GEMStripTopology::yAxisOrientation ( ) const

inline

Definition at line 59 of file GEMStripTopology.h.

{ return yAxisOrientation_; }

References yAxisOrientation_.

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

yCentreOfStripPlane()

float GEMStripTopology::yCentreOfStripPlane ( ) const

inline

Definition at line 61 of file GEMStripTopology.h.

{ return yCentre_; }

References yCentre_.

Referenced by localPosition(), and measurementPosition().

yDistanceToIntersection()

float GEMStripTopology::yDistanceToIntersection

(

float

y

)

const

Definition at line 128 of file GEMStripTopology.cc.

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

References originToIntersection(), y, and yAxisOrientation().

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

yExtentOfStripPlane()

float GEMStripTopology::yExtentOfStripPlane ( ) const

inline

Definition at line 52 of file GEMStripTopology.h.

{ return detHeight_; }

References detHeight_.

Member Data Documentation

angularWidth_

float GEMStripTopology::angularWidth_

private

Definition at line 65 of file GEMStripTopology.h.

Referenced by angularWidth(), and GEMStripTopology().

centreToIntersection_

float GEMStripTopology::centreToIntersection_

private

Definition at line 67 of file GEMStripTopology.h.

Referenced by centreToIntersection(), originToIntersection(), and radius().

detHeight_

float GEMStripTopology::detHeight_

private

Definition at line 66 of file GEMStripTopology.h.

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

numberOfStrips_

int GEMStripTopology::numberOfStrips_

private

Definition at line 64 of file GEMStripTopology.h.

Referenced by channel(), GEMStripTopology(), and nstrips().

phiOfOneEdge_

float GEMStripTopology::phiOfOneEdge_

private

Definition at line 68 of file GEMStripTopology.h.

Referenced by GEMStripTopology(), and phiOfOneEdge().

yAxisOrientation_

float GEMStripTopology::yAxisOrientation_

private

Definition at line 69 of file GEMStripTopology.h.

Referenced by GEMStripTopology(), and yAxisOrientation().

yCentre_

float GEMStripTopology::yCentre_

private

Definition at line 70 of file GEMStripTopology.h.

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