TrapezoidalStripTopology Class Reference

#include <TrapezoidalStripTopology.h>

Inheritance diagram for TrapezoidalStripTopology:

Public Member Functions
int	channel (const LocalPoint &) const override
LocalError	localError (float strip, float stripErr2) const override
LocalError	localError (const MeasurementPoint &, const MeasurementError &) const override
float	localPitch (const LocalPoint &) const override
LocalPoint	localPosition (float strip) const override
LocalPoint	localPosition (const MeasurementPoint &) const override
float	localStripLength (const LocalPoint &aLP) const override
MeasurementError	measurementError (const LocalPoint &, const LocalError &) const override
MeasurementPoint	measurementPosition (const LocalPoint &) const override
int	nstrips () const override
float	pitch () const override
float	radius () const
float	strip (const LocalPoint &) const override
float	stripAngle (float strip) const override
float	stripLength () const override
	det heigth (strip length in the middle) More...
	TrapezoidalStripTopology (int nstrip, float pitch, float detheight, float r0)
	TrapezoidalStripTopology (int nstrip, float pitch, float detheight, float r0, int yAx)
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 float	localPitch (const LocalPoint &lp, const Topology::LocalTrackAngles &) const
	conversion taking also the angle from the track state (LocalTrajectoryParameters) 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 ()

Protected Member Functions
virtual float	shiftOffset (float pitch_fraction)

Private Attributes
float	theDetHeight
float	theDistToBeam
int	theNumberOfStrips
float	theOffset
float	thePitch
int	theYAxOr

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

Detailed Description

Specialised strip topology for rectangular barrel detectors. The strips are parallel to the local Y axis, so X is the precisely measured coordinate.Specialization of StripTopology for detectors of symmetric trapezoidal shape. The local Y coordinate is parallel to the central strip, and prpendicular to the paralle sides of the trapezoid. The first and last strips are parallel to the corresponding sides of the trapezoid. The pitch is constant at constant y. This topology makes a non-linear transformation: the pitch is not constant along the strips.

Definition at line 21 of file TrapezoidalStripTopology.h.

Constructor & Destructor Documentation

TrapezoidalStripTopology::TrapezoidalStripTopology	(	int	nstrip,
		float	pitch,
		float	detheight,
		float	r0
	)

constructed from: number of strips pitch in the middle det heigth (strip length in the middle) radius of circle passing through the middle of the det with center at the crossing of the two sides.

Definition at line 7 of file TrapezoidalStripTopology.cc.

References gather_cfg::cout, theNumberOfStrips, theOffset, thePitch, and theYAxOr.

    : theNumberOfStrips(ns), thePitch(p), theDistToBeam(r0), theDetHeight(l) {
  theOffset = -theNumberOfStrips / 2. * thePitch;
  theYAxOr = 1;
#ifdef VERBOSE
  cout << "Constructing TrapezoidalStripTopology with"
       << " nstrips = " << ns << " pitch = " << p << " length = " << l << " r0 =" << r0 << endl;
#endif
}

TrapezoidalStripTopology::TrapezoidalStripTopology	(	int	nstrip,
		float	pitch,
		float	detheight,
		float	r0,
		int	yAx
	)

Definition at line 17 of file TrapezoidalStripTopology.cc.

References gather_cfg::cout, theNumberOfStrips, theOffset, and thePitch.

    : theNumberOfStrips(ns), thePitch(p), theDistToBeam(r0), theDetHeight(l), theYAxOr(yAx) {
  theOffset = -theNumberOfStrips / 2. * thePitch;
#ifdef VERBOSE
  cout << "Constructing TrapezoidalStripTopology with"
       << " nstrips = " << ns << " pitch = " << p << " length = " << l << " r0 =" << r0 << " yAxOrientation =" << yAx
       << endl;
#endif
}

Member Function Documentation

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

overridevirtual

Implements Topology.

Definition at line 102 of file TrapezoidalStripTopology.cc.

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

Referenced by GEMSignalModel::simulateClustering().

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

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

overridevirtual

Implements StripTopology.

Definition at line 37 of file TrapezoidalStripTopology.cc.

References theDetHeight, theDistToBeam, theOffset, thePitch, and theYAxOr.

                                                                                  {
  float lt, lc2, ls2, lslc;
  float localL2, localP2;
  float sl2, sp2;
  // angle from strip to local frame (see CMS TN / 95-170)
  lt = -(strip * thePitch + theOffset) * theYAxOr / theDistToBeam;
  lc2 = 1.f / (1. + lt * lt);
  lslc = lt * lc2;
  ls2 = 1.f - lc2;
  localL2 = theDetHeight * theDetHeight / lc2;
  localP2 = thePitch * thePitch * lc2;
  sl2 = localL2 / 12.;
  sp2 = stripErr2 * localP2;
  return LocalError(lc2 * sp2 + ls2 * sl2, lslc * (sp2 - sl2), ls2 * sp2 + lc2 * sl2);
}

LocalError TrapezoidalStripTopology::localError ( const MeasurementPoint &  mp, const MeasurementError &  merr  ) const

overridevirtual

Implements Topology.

Definition at line 53 of file TrapezoidalStripTopology.cc.

References localPitch(), localPosition(), mathSSE::sqrt(), theDetHeight, theDistToBeam, theOffset, thePitch, theYAxOr, MeasurementError::uu(), MeasurementError::uv(), MeasurementError::vv(), and PV2DBase< T, PVType, FrameType >::x().

                                                                                                              {
  float lt, lc2, ls2, lslc;
  float localL, localP;
  float sl2, sp2, spl;
  // angle from strip to local frame (see CMS TN / 95-170)
  lt = -(mp.x() * thePitch + theOffset) * theYAxOr / theDistToBeam;
  lc2 = 1. / (1. + lt * lt);
  lslc = lt * lc2;
  ls2 = 1.f - lc2;
  localL = theDetHeight / std::sqrt(lc2);
  localP = localPitch(localPosition(mp));
  sp2 = merr.uu() * localP * localP;
  sl2 = merr.vv() * localL * localL;
  spl = merr.uv() * localP * localL;
  return LocalError(lc2 * sp2 + ls2 * sl2 - 2 * lslc * spl,
                    lslc * (sp2 - sl2) + (lc2 - ls2) * spl,
                    ls2 * sp2 + lc2 * sl2 + 2 * lslc * spl);
}

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

overridevirtual

Implements StripTopology.

Definition at line 108 of file TrapezoidalStripTopology.cc.

References f, mathSSE::sqrt(), theDistToBeam, thePitch, theYAxOr, x, PV3DBase< T, PVType, FrameType >::x(), y, and PV3DBase< T, PVType, FrameType >::y().

Referenced by localError(), and measurementError().

                                                                     {
  float x = lp.x();
  float y = theYAxOr * lp.y() + theDistToBeam;
  return thePitch * y / (theDistToBeam * std::sqrt(1.f + x * x / (y * y)));
}

LocalPoint TrapezoidalStripTopology::localPosition ( float  strip ) const

overridevirtual

Implements StripTopology.

Definition at line 27 of file TrapezoidalStripTopology.cc.

References theOffset, and thePitch.

Referenced by CSCSegtoRPC::CSCSegtoRPC(), localError(), RPCSimSimple::simulateNoise(), RPCSimParam::simulateNoise(), RPCSimAverage::simulateNoise(), RPCSimAverageNoise::simulateNoise(), RPCSimAverageNoiseEffCls::simulateNoise(), RPCSimModelTiming::simulateNoise(), RPCSimAverageNoiseEff::simulateNoise(), RPCSimAsymmetricCls::simulateNoise(), and TracktoRPC::TracktoRPC().

                                                                    {
  return LocalPoint(strip * thePitch + theOffset, 0.0);
}

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

overridevirtual

Implements Topology.

Definition at line 31 of file TrapezoidalStripTopology.cc.

References theDetHeight, theDistToBeam, theOffset, thePitch, theYAxOr, x, PV2DBase< T, PVType, FrameType >::x(), y, and PV2DBase< T, PVType, FrameType >::y().

                                                                                   {
  float y = mp.y() * theDetHeight;
  float x = (mp.x() * thePitch + theOffset) * (theYAxOr * y + theDistToBeam) / theDistToBeam;
  return LocalPoint(x, y);
}

float TrapezoidalStripTopology::localStripLength ( const LocalPoint &  aLP ) const

overridevirtual

Implements StripTopology.

Definition at line 125 of file TrapezoidalStripTopology.cc.

References f, mathSSE::sqrt(), theDetHeight, theDistToBeam, theYAxOr, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by stripLength().

                                                                           {
  float ltan = -lp.x() / (theYAxOr * lp.y() + theDistToBeam) * theYAxOr;
  float localL = theDetHeight * std::sqrt(1.f + ltan * ltan);
  //  float lcos2 = 1.f/(1.f+ltan*ltan);
  //  float localL = theDetHeight / std::sqrt(lcos2);

  return localL;
}

MeasurementError TrapezoidalStripTopology::measurementError ( const LocalPoint &  lp, const LocalError &  lerr  ) const

overridevirtual

Implements Topology.

Definition at line 86 of file TrapezoidalStripTopology.cc.

References localPitch(), mathSSE::sqrt(), theDetHeight, theDistToBeam, theYAxOr, PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), LocalError::xy(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

                                                                                                              {
  float lt, lc2, ls2, lslc;
  float localL, localP;
  float sl2, sp2, spl;
  lt = -lp.x() / (theYAxOr * lp.y() + theDistToBeam) * theYAxOr;
  lc2 = 1. / (1. + lt * lt);
  lslc = lt * lc2;
  ls2 = 1. - lc2;
  localL = theDetHeight / std::sqrt(lc2);
  localP = localPitch(lp);
  sp2 = lc2 * lerr.xx() + ls2 * lerr.yy() + 2 * lslc * lerr.xy();
  sl2 = ls2 * lerr.xx() + lc2 * lerr.yy() - 2 * lslc * lerr.xy();
  spl = lslc * (lerr.yy() - lerr.xx()) + (lc2 - ls2) * lerr.xy();
  return MeasurementError(sp2 / (localP * localP), spl / (localP * localL), sl2 / (localL * localL));
}

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

overridevirtual

Implements Topology.

Definition at line 81 of file TrapezoidalStripTopology.cc.

References theDetHeight, theDistToBeam, theOffset, thePitch, theYAxOr, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

                                                                                         {
  return MeasurementPoint(((lp.x() * theDistToBeam / (theYAxOr * lp.y() + theDistToBeam)) - theOffset) / thePitch,
                          lp.y() / theDetHeight);
}

int TrapezoidalStripTopology::nstrips ( ) const

overridevirtual

Implements StripTopology.

Definition at line 118 of file TrapezoidalStripTopology.cc.

References theNumberOfStrips.

Referenced by ME0ReDigiProducer::getStripProperties().

{ return theNumberOfStrips; }

float TrapezoidalStripTopology::pitch ( ) const

overridevirtual

Pitch in the middle of the DetUnit

Implements StripTopology.

Definition at line 106 of file TrapezoidalStripTopology.cc.

References thePitch.

Referenced by CSCSegtoRPC::CSCSegtoRPC(), and ME0ReDigiProducer::getStripProperties().

{ return thePitch; }

float TrapezoidalStripTopology::radius ( void  ) const

inline

radius of circle passing through the middle of the det with center at the crossing of the two sides.

Definition at line 69 of file TrapezoidalStripTopology.h.

References shiftOffset(), and theDistToBeam.

Referenced by ME0ReDigiProducer::getStripProperties(), and ME0ReDigiProducer::TemporaryGeometry::TemporaryGeometry().

{ return theDistToBeam; }

float TrapezoidalStripTopology::shiftOffset ( float  pitch_fraction )

protectedvirtual

Definition at line 120 of file TrapezoidalStripTopology.cc.

References theOffset, and thePitch.

Referenced by radius().

                                                                {
  theOffset += thePitch * pitch_fraction;
  return theOffset;
}

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

overridevirtual

Implements StripTopology.

Definition at line 72 of file TrapezoidalStripTopology.cc.

References theDistToBeam, theNumberOfStrips, theOffset, thePitch, theYAxOr, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by channel().

                                                                {
  float aStrip = ((lp.x() * theDistToBeam / (theYAxOr * lp.y() + theDistToBeam)) - theOffset) / thePitch;
  if (aStrip < 0)
    aStrip = 0;
  else if (aStrip > theNumberOfStrips)
    aStrip = theNumberOfStrips;
  return aStrip;
}

float TrapezoidalStripTopology::stripAngle ( float  strip ) const

overridevirtual

angle between strip and symmetry axis

Implements StripTopology.

Definition at line 114 of file TrapezoidalStripTopology.cc.

References theDistToBeam, theOffset, thePitch, and theYAxOr.

                                                            {
  return std::atan(-(strip * thePitch + theOffset) * theYAxOr / theDistToBeam);
}

float TrapezoidalStripTopology::stripLength ( ) const

inlineoverridevirtual

det heigth (strip length in the middle)

Implements StripTopology.

Definition at line 63 of file TrapezoidalStripTopology.h.

References localStripLength(), and theDetHeight.

Referenced by CSCSegtoRPC::CSCSegtoRPC(), RPCSynchronizer::getSimHitBx(), RPCSynchronizer::getSimHitBxAndTimingForIRPC(), ME0ReDigiProducer::getStripProperties(), RPCSimSimple::simulateNoise(), RPCSimParam::simulateNoise(), RPCSimAverage::simulateNoise(), RPCSimAverageNoise::simulateNoise(), RPCSimAverageNoiseEff::simulateNoise(), RPCSimAverageNoiseEffCls::simulateNoise(), RPCSimModelTiming::simulateNoise(), RPCSimAsymmetricCls::simulateNoise(), ME0ReDigiProducer::TemporaryGeometry::TemporaryGeometry(), and TracktoRPC::TracktoRPC().

{ return theDetHeight; }

Member Data Documentation

float TrapezoidalStripTopology::theDetHeight

private

Definition at line 79 of file TrapezoidalStripTopology.h.

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

float TrapezoidalStripTopology::theDistToBeam

private

Definition at line 78 of file TrapezoidalStripTopology.h.

Referenced by localError(), localPitch(), localPosition(), localStripLength(), measurementError(), measurementPosition(), radius(), strip(), and stripAngle().

int TrapezoidalStripTopology::theNumberOfStrips

private

Definition at line 75 of file TrapezoidalStripTopology.h.

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

float TrapezoidalStripTopology::theOffset

private

Definition at line 77 of file TrapezoidalStripTopology.h.

Referenced by localError(), localPosition(), measurementPosition(), shiftOffset(), strip(), stripAngle(), and TrapezoidalStripTopology().

float TrapezoidalStripTopology::thePitch

private

Definition at line 76 of file TrapezoidalStripTopology.h.

Referenced by localError(), localPitch(), localPosition(), measurementPosition(), pitch(), shiftOffset(), strip(), stripAngle(), and TrapezoidalStripTopology().

int TrapezoidalStripTopology::theYAxOr

private

Definition at line 80 of file TrapezoidalStripTopology.h.

Referenced by localError(), localPitch(), localPosition(), localStripLength(), measurementError(), measurementPosition(), strip(), stripAngle(), and TrapezoidalStripTopology().