RadialStripTopology.h

Go to the documentation of this file.

#ifndef _RADIAL_STRIP_TOPOLOGY_H_
#define _RADIAL_STRIP_TOPOLOGY_H_

#include "Geometry/CommonTopologies/interface/StripTopology.h"

class RadialStripTopology : public StripTopology {
 public:

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

  virtual ~RadialStripTopology(){}

  // =========================================================
  // StripTopology interface - implement pure virtual methods
  // =========================================================

  virtual LocalPoint localPosition(float strip) const;

  virtual LocalPoint localPosition(const MeasurementPoint&) const;

  virtual LocalError localError(float strip, float stripErr2) const;

  virtual LocalError localError(const MeasurementPoint&, const MeasurementError&) const;

  virtual float strip(const LocalPoint&) const;


  virtual float pitch() const;

  virtual float localPitch(const LocalPoint&) const;

  virtual float stripAngle(float strip) const;

  virtual int nstrips() const { return theNumberOfStrips; }

  virtual float stripLength() const { return theDetHeight; }

  virtual float localStripLength(const LocalPoint& ) const;


  // =========================================================
  // Topology interface (not already implemented for 
  // StripTopology interface)
  // =========================================================

  virtual MeasurementPoint measurementPosition( const LocalPoint& ) const;

  virtual MeasurementError measurementError( const LocalPoint&, const LocalError& ) const;

  virtual int channel( const LocalPoint& ) const;


  // =========================================================
  // RadialStripTopology interface itself
  // =========================================================

  float angularWidth() const { return theAngularWidth;}

  virtual float phiPitch(void) const { return angularWidth(); }

  float detHeight() const { return theDetHeight;}

  float yExtentOfStripPlane() const { return theDetHeight; } // same as detHeight()

  float centreToIntersection() const { return theCentreToIntersection; }

  float originToIntersection() const { return (theCentreToIntersection - yCentre); }

  float phiOfOneEdge() const { return thePhiOfOneEdge; }

  float xOfStrip(int strip, float y) const;
 
  virtual int nearestStrip(const LocalPoint&) const;

  int yAxisOrientation() const { return theYAxisOrientation; }

  float yCentreOfStripPlane() const { return yCentre; }

  float yDistanceToIntersection( float y ) const;

  friend std::ostream & operator<<(std::ostream&, const RadialStripTopology& );

 private:

  int   theNumberOfStrips; // total no. of strips in plane of strips
  float theAngularWidth;   // angle subtended by each strip = phi pitch
  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!)
  int   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.
};

#endif