/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/RecoTracker/NuclearSeedGenerator/interface/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);

     bool isValid() const { return isValid(); } 

 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