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