CSCWireTopology.h

Go to the documentation of this file.

#ifndef CSC_WIRE_TOPOLOGY_H
#define CSC_WIRE_TOPOLOGY_H

#include "Geometry/CSCGeometry/interface/WireTopology.h"
#include "Geometry/CSCGeometry/interface/CSCWireGeometry.h"
#include "Geometry/CSCGeometry/src/CSCWireGrouping.h"
#include "Geometry/CSCGeometry/src/CSCWireGroupPackage.h"
#include <utility> // for std::pair

class CSCWireTopology : public WireTopology {

 public:

  virtual ~CSCWireTopology();
  
  CSCWireTopology( const CSCWireGroupPackage& wg,
             double yOfFirstWire,
                     float wireAngleInDegrees );
  CSCWireTopology( const CSCWireTopology& );

  CSCWireTopology& operator=( const CSCWireTopology& );

  virtual LocalPoint localPosition( const MeasurementPoint& ) const;
  virtual LocalError localError( const MeasurementPoint&, const MeasurementError& ) const;
  virtual MeasurementPoint measurementPosition( const LocalPoint&) const;
  virtual MeasurementError measurementError( const LocalPoint&, const LocalError& ) const;

  virtual int channel( const LocalPoint& p) const;

  double wireSpacing() const {
    return theWireGeometry->wireSpacing(); }

  float wirePitch() const {
    return static_cast<float>(wireSpacing()); }

  float wireAngle() const { 
    return theWireGeometry->wireAngle(); }

  int nearestWire(const LocalPoint& lp) const {
    return theWireGeometry->nearestWire( lp ); }

  float yOfWire(float wire, float x=0.) const {
    return theWireGeometry->yOfWire(wire,x); }

  double narrowWidthOfPlane() const {
    return theWireGeometry->narrowWidthOfPlane(); }

  double wideWidthOfPlane() const {
    return theWireGeometry->wideWidthOfPlane(); }

  double lengthOfPlane() const {
    return theWireGeometry->lengthOfPlane(); }


  int numberOfWires() const {
    return theWireGrouping->numberOfWires(); }

  int numberOfWireGroups() const {
    return theWireGrouping->numberOfWireGroups(); }

  int numberOfWiresPerGroup( int wireGroup ) const {
    return theWireGrouping->numberOfWiresPerGroup( wireGroup ); }

  int wireGroup(int wire) const {
    return theWireGrouping->wireGroup( wire ); }

  float middleWireOfGroup( int wireGroup ) const {
    return theWireGrouping->middleWireOfGroup( wireGroup ); }

  float yOfWireGroup(int wireGroup, float x=0.) const;
 
  float yResolution( int wireGroup=1 ) const;

  double extentOfWirePlane() const {
    return wireSpacing() * (numberOfWires() - 1) ; }

  std::pair< LocalPoint, LocalPoint > wireEnds( float wire ) const {
    return theWireGeometry->wireEnds( wire ); }

  std::vector<float> wireValues( float wire ) const {
    return theWireGeometry->wireValues( wire ); }

  std::pair<float, float> equationOfWire( float wire ) const;
  
  float restrictToYOfWirePlane( float y ) const;
  
  bool insideYOfWirePlane( float y ) const;
  
 private:

  CSCWireGrouping* theWireGrouping; // handles grouping of wires for read out
  CSCWireGeometry* theWireGeometry; // handles non-zero angle w.r.t x axis

  double theAlignmentPinToFirstWire; //@@ Not sure this is actually required!

};

#endif