TkRadialStripTopology.h

Go to the documentation of this file.

#ifndef _TkRADIAL_STRIP_TOPOLOGY_H_
#define _TkRADIAL_STRIP_TOPOLOGY_H_
 
#include "Geometry/CommonTopologies/interface/RadialStripTopology.h"
 
class TkRadialStripTopology final : public RadialStripTopology {
public:
  TkRadialStripTopology(int ns, float aw, float dh, float r, int yAx = 1, float yMid = 0.);
 
  ~TkRadialStripTopology() override {}
 
  // =========================================================
  // StripTopology interface - implement pure methods
  // =========================================================
 
  LocalPoint localPosition(float strip) const override;
 
  LocalPoint localPosition(const MeasurementPoint&) const override;
 
  LocalError localError(float strip, float stripErr2) const override;
 
  LocalError localError(const MeasurementPoint&, const MeasurementError&) const override;
 
  float strip(const LocalPoint&) const override;
 
  // the number of strip span by the segment between the two points..
  float coveredStrips(const LocalPoint& lp1, const LocalPoint& lp2) const override;
 
  float localPitch(const LocalPoint&) const override;
 
  float stripAngle(float strip) const override {
    return yAxisOrientation() * (phiOfOneEdge() + strip * angularWidth());
  }
 
  int nstrips() const override { return theNumberOfStrips; }
 
  float stripLength() const override { return theDetHeight; }
 
  float localStripLength(const LocalPoint&) const override;
 
  // =========================================================
  // Topology interface (not already implemented for
  // StripTopology interface)
  // =========================================================
 
  MeasurementPoint measurementPosition(const LocalPoint&) const override;
 
  MeasurementError measurementError(const LocalPoint&, const LocalError&) const override;
 
  int channel(const LocalPoint&) const override;
 
  // =========================================================
  // RadialStripTopology interface itself
  // =========================================================
 
  float angularWidth() const override { return theAngularWidth; }
 
  float phiPitch(void) const override { return angularWidth(); }
 
  float detHeight() const override { return theDetHeight; }
 
  float yExtentOfStripPlane() const override { return theDetHeight; }  // same as detHeight()
 
  float centreToIntersection() const override { return theCentreToIntersection; }
 
  float originToIntersection() const override { return (theCentreToIntersection - yCentre); }
 
  float phiOfOneEdge() const override { return thePhiOfOneEdge; }
 
  float xOfStrip(int strip, float y) const override;
 
  int nearestStrip(const LocalPoint&) const override;
 
  float yAxisOrientation() const override { return theYAxisOrientation; }
 
  float yCentreOfStripPlane() const override { return yCentre; }
 
  float yDistanceToIntersection(float y) const override;
 
private:
  int theNumberOfStrips;          // total no. of strips in plane of strips
  float theAngularWidth;          // angle subtended by each strip = phi pitch
  float theAWidthInverse;         // inverse of above
  float theTanAW;                 // its tangent
  float theDetHeight;             // length of long symmetry axis = twice the apothem of the enclosing trapezoid
  float theCentreToIntersection;  // distance centre of detector face to intersection of edge strips (projected)
  float thePhiOfOneEdge;          // local 'phi' of one edge of plane of strips (I choose it negative!)
  float theTanOfOneEdge;          // the positive tangent of the above...
  float theYAxisOrientation;      // 1 means y axis going from smaller to larger side, -1 means opposite direction
  float yCentre;  // Non-zero if offset in local y between midpoint of detector (strip plane) extent and local origin.
  double theRadialSigma;  // radial sigma^2( uniform prob density along strip)
};
 
#endif