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)
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...
LocalError	localError (float strip, float stripErr2) const override
LocalError	localError (const MeasurementPoint &, const MeasurementError &) const override
float	localPitch (const LocalPoint &) const override
virtual LocalPoint	localPosition (float strip, const Topology::LocalTrackPred &) const
	conversion taking also the predicted track state More...
virtual LocalPoint	localPosition (float strip) const=0
LocalPoint	localPosition (float strip) const override
LocalPoint	localPosition (const MeasurementPoint &) const override
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 (float strip, float stripErr2, const Topology::LocalTrackPred &) const
	conversion taking also the angle from the predicted track state More...
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 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.

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

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

GEMStripTopology() [2/2]

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

Definition at line 23 of file GEMStripTopology.cc.

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

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

~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.

References angularWidth_.

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

{ return angularWidth_; }

centreToIntersection()

float GEMStripTopology::centreToIntersection ( ) const

inline

Definition at line 53 of file GEMStripTopology.h.

References centreToIntersection_.

Referenced by localError().

{ return centreToIntersection_; }

channel()

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

overridevirtual

Implements Topology.

Definition at line 107 of file GEMStripTopology.cc.

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

Referenced by GEMSignalModel::simulateClustering().

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

detHeight()

float GEMStripTopology::detHeight ( ) const

inline

Definition at line 51 of file GEMStripTopology.h.

References detHeight_.

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

{ return detHeight_; }

localError() [1/4]

virtual LocalError StripTopology::localError

localError() [2/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);
  }

localError() [3/4]

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

overridevirtual

Implements StripTopology.

Definition at line 45 of file GEMStripTopology.cc.

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.

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

localError() [4/4]

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

overridevirtual

Implements Topology.

Definition at line 55 of file GEMStripTopology.cc.

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

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

localPitch()

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

overridevirtual

Implements StripTopology.

Definition at line 111 of file GEMStripTopology.cc.

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

Referenced by pitch().

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

localPosition() [1/4]

virtual LocalPoint StripTopology::localPosition

localPosition() [2/4]

virtual LocalPoint StripTopology::localPosition

inline

conversion taking also the predicted track state

Definition at line 20 of file StripTopology.h.

                                                                                       {
    return localPosition(strip);
  }

localPosition() [3/4]

LocalPoint GEMStripTopology::localPosition ( float  strip ) const

overridevirtual

Implements StripTopology.

Definition at line 34 of file GEMStripTopology.cc.

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

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

localPosition() [4/4]

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

overridevirtual

Implements Topology.

Definition at line 38 of file GEMStripTopology.cc.

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

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

localStripLength()

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

overridevirtual

Implements StripTopology.

Definition at line 123 of file GEMStripTopology.cc.

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

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

measurementError()

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

overridevirtual

Implements Topology.

Definition at line 91 of file GEMStripTopology.cc.

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

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

measurementPosition()

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

overridevirtual

Implements Topology.

Definition at line 84 of file GEMStripTopology.cc.

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

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

nearestStrip()

int GEMStripTopology::nearestStrip

(

const LocalPoint &

lp

)

const

Definition at line 80 of file GEMStripTopology.cc.

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

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

nstrips()

int GEMStripTopology::nstrips ( ) const

inlineoverridevirtual

Implements StripTopology.

Definition at line 44 of file GEMStripTopology.h.

References numberOfStrips_.

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

{ return numberOfStrips_; }

originToIntersection()

float GEMStripTopology::originToIntersection ( ) const

inline

Definition at line 55 of file GEMStripTopology.h.

References centreToIntersection_, and yCentre_.

Referenced by localPosition(), and yDistanceToIntersection().

{ return (centreToIntersection_ - yCentre_); }

phiOfOneEdge()

float GEMStripTopology::phiOfOneEdge ( ) const

inline

Definition at line 60 of file GEMStripTopology.h.

References phiOfOneEdge_.

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

{ return phiOfOneEdge_; }

phiPitch()

float GEMStripTopology::phiPitch ( void  ) const

inline

Definition at line 36 of file GEMStripTopology.h.

References angularWidth().

{ return angularWidth(); }

pitch()

float GEMStripTopology::pitch ( ) const

overridevirtual

Implements StripTopology.

Definition at line 109 of file GEMStripTopology.cc.

References localPitch().

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

radius()

float GEMStripTopology::radius ( void  ) const

inline

Definition at line 54 of file GEMStripTopology.h.

References centreToIntersection_.

{ return centreToIntersection_; }

strip()

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

overridevirtual

Implements StripTopology.

Definition at line 73 of file GEMStripTopology.cc.

References angularWidth(), SiStripPI::max, SiStripPI::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().

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

stripAngle()

float GEMStripTopology::stripAngle ( float  strip ) const

overridevirtual

Implements StripTopology.

Definition at line 119 of file GEMStripTopology.cc.

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

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

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

stripLength()

float GEMStripTopology::stripLength ( ) const

inlineoverridevirtual

Implements StripTopology.

Definition at line 46 of file GEMStripTopology.h.

References detHeight_.

Referenced by GEMDQMBase::readGeometryRadiusInfoChamber().

{ return detHeight_; }

xOfStrip()

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

Definition at line 132 of file GEMStripTopology.cc.

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

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

yAxisOrientation()

float GEMStripTopology::yAxisOrientation ( ) const

inline

Definition at line 59 of file GEMStripTopology.h.

References yAxisOrientation_.

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

{ return yAxisOrientation_; }

yCentreOfStripPlane()

float GEMStripTopology::yCentreOfStripPlane ( ) const

inline

Definition at line 61 of file GEMStripTopology.h.

References yCentre_.

Referenced by localPosition(), and measurementPosition().

{ return yCentre_; }

yDistanceToIntersection()

float GEMStripTopology::yDistanceToIntersection

(

float

y

)

const

Definition at line 128 of file GEMStripTopology.cc.

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

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

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

yExtentOfStripPlane()

float GEMStripTopology::yExtentOfStripPlane ( ) const

inline

Definition at line 52 of file GEMStripTopology.h.

References detHeight_.

{ return 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().