HelixArbitraryPlaneCrossing.h

Go to the documentation of this file.

#ifndef HELIXARBITRARYPLANECROSSING_H_
#define HELIXARBITRARYPLANECROSSING_H_
#include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"
#include "TrackingTools/GeomPropagators/interface/HelixPlaneCrossing.h"
#include "TrackingTools/GeomPropagators/interface/HelixArbitraryPlaneCrossing2Order.h"

class HelixArbitraryPlaneCrossing : public HelixPlaneCrossing {
public:
  HelixArbitraryPlaneCrossing(const PositionType& point,
                    const DirectionType& direction,
                    const float curvature,
                        const PropagationDirection propDir = alongMomentum);
  // destructor
  virtual ~HelixArbitraryPlaneCrossing() {}

  virtual std::pair<bool,double> pathLength(const Plane& plane);

  virtual PositionType position(double s) const;

  virtual DirectionType direction(double s) const;
  //
  // double precision vectors for internal use
  //
  typedef Basic3DVector<double>  PositionTypeDouble;
  typedef Basic3DVector<double>  DirectionTypeDouble;

  PositionTypeDouble positionInDouble(double s) const;

  DirectionTypeDouble directionInDouble(double s) const;

private:
  inline bool notAtSurface (const Plane&,
                const PositionTypeDouble&,
                const float) const;

private:
  HelixArbitraryPlaneCrossing2Order theQuadraticCrossingFromStart;


  const double theX0,theY0,theZ0;
  double theCosPhi0,theSinPhi0;
  double theCosTheta,theSinTheta;
  const double theRho;

  const PropagationDirection thePropDir;

  mutable double theCachedS;
  mutable double theCachedDPhi;
  mutable double theCachedSDPhi;
  mutable double theCachedCDPhi;

  static const float theNumericalPrecision;
  static const float theMaxDistToPlane;

};

#endif