/data/refman/pasoursint/CMSSW_5_2_9/src/Geometry/CommonTopologies/interface/Topology.h

Go to the documentation of this file.

#ifndef Geometry_CommonTopologies_Topology_H
#define Geometry_CommonTopologies_Topology_H

#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/LocalError.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementPoint.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementError.h"
#include "DataFormats/Math/interface/AlgebraicROOTObjects.h"

class GeomDetType;

class Topology {
public:
  typedef Basic2DVector<double> Vector2D;
  typedef Vector2D::MathVector MathVector2D;
  class LocalTrackAngles : public Vector2D {
  public:
    typedef  Basic2DVector<double> Base;
    LocalTrackAngles(const Base & v) : Base(v) {}
    LocalTrackAngles(double dxdz, double dydz) : Base(dxdz,dydz) {}
    double dxdz() const { return x(); }
    double dydz() const { return y(); }
  };
    typedef Point2DBase <double, LocalTag> Local2DPoint;
    class LocalTrackPred {
    public:
    LocalTrackPred(double x, double y, double dxdz, double dydz) : point_(x,y),angles_(dxdz,dydz){}
    LocalTrackPred(const AlgebraicVector5& localTrajPar) :
      point_(localTrajPar[3], localTrajPar[4]), angles_(localTrajPar[1], localTrajPar[2]) {}
      const Local2DPoint& point() const {return point_;}
      const LocalTrackAngles& angles() const {return angles_;}
    private:
      Local2DPoint     point_; 
      LocalTrackAngles angles_;
    };
    
  virtual ~Topology() {}
  
  // Conversion between measurement (strip, pixel, ...) coordinates
  // and local cartesian coordinates

  virtual LocalPoint localPosition( const MeasurementPoint& ) const = 0;

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

  virtual MeasurementPoint measurementPosition( const LocalPoint&) const = 0;

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

  virtual int channel( const LocalPoint& p) const = 0;

  // new sets of methods taking also an angle
  virtual LocalPoint localPosition(const MeasurementPoint &mp, const LocalTrackPred &/*trkPred*/) const {
      return localPosition(mp);
  }

  virtual LocalError 
  localError(const MeasurementPoint &mp, const MeasurementError &me,
             const LocalTrackPred &/*trkPred*/) const {
      return localError(mp,me);
  }

  virtual MeasurementPoint measurementPosition( const LocalPoint &lp, const LocalTrackAngles &/*dir*/) const {
      return measurementPosition(lp);
  }

  virtual MeasurementError 
  measurementError( const LocalPoint &lp, const LocalError &le, const LocalTrackAngles &/*dir*/) const {
      return measurementError(lp, le);
  }

  virtual int channel( const LocalPoint &lp, const LocalTrackAngles &/*dir*/) const {
      return channel(lp);
  }


private:

};

#endif