TangentCircle.h

Go to the documentation of this file.

#ifndef _TangentCircle_H_
#define _TangentCircle_H_
 
#include "RecoTracker/TkSeedGenerator/interface/FastCircle.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
 
class TangentCircle {
  //   TODO : for speed-up :
  //   typedef Point2DBase< float, GlobalTag>     Global2DPoint;
  //   typedef Vector2DBase< float, GlobalTag>    Global2DVector;
 
public:
  TangentCircle()
      : theInnerPoint(),
        theOuterPoint(),
        theVertexPoint(),
        theDirectionAtVertex(),
        theX0(0),
        theY0(0),
        theRho(0),
        theVertexError(0),
        valid(false),
        theCharge(0) {}
 
  TangentCircle(const GlobalVector& direction, const GlobalPoint& innerPoint, const GlobalPoint& outerPoint);
 
  TangentCircle(const GlobalPoint& outerPoint, const GlobalPoint& innerPoint, const GlobalPoint& vertexPoint);
 
  TangentCircle(const TangentCircle& primCircle, const GlobalPoint& outerPoint, const GlobalPoint& innerPoint);
 
  GlobalVector directionAtVertex();
 
  double x0() const { return theX0; }
 
  double y0() const { return theY0; }
 
  double rho() const { return theRho; }
 
  GlobalPoint outerPoint() const { return theOuterPoint; }
 
  GlobalPoint innerPoint() const { return theInnerPoint; }
 
  GlobalPoint vertexPoint() const { return theVertexPoint; }
 
  double vertexError() const { return theVertexError; }
 
  double curvatureError();
 
  int charge(float magz);
 
private:
  GlobalPoint theInnerPoint;
  GlobalPoint theOuterPoint;
  GlobalPoint theVertexPoint;
 
  GlobalVector theDirectionAtVertex;
 
  double theX0;  
  double theY0;  
  double theRho; 
  double theVertexError; 
  bool valid;
  int theCharge;
 
  double isTangent(const TangentCircle& primCircle, const TangentCircle& secCircle) const;
  GlobalPoint getPosition(const TangentCircle& circle,
                          const GlobalPoint& initalPosition,
                          double theta,
                          int direction) const;
  int chargeLocally(float magz, GlobalVector v) const;
  GlobalVector direction(const GlobalPoint& point) const;
};
 
#endif