IterativeHelixExtrapolatorToLine.h

Go to the documentation of this file.

#ifndef IterativeHelixExtrapolatorToLine_h_
#define IterativeHelixExtrapolatorToLine_h_

#include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"
#include "TrackingTools/GeomPropagators/interface/HelixLineExtrapolation.h"
#include "TrackingTools/GeomPropagators/interface/HelixExtrapolatorToLine2Order.h"

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

  virtual std::pair<bool,double> pathLength (const GlobalPoint& point) const {
    return genericPathLength(point);
  }

  virtual std::pair<bool,double> pathLength (const Line& line) const {
    return genericPathLength(line);
  }

  virtual PositionType position (double s) const;

  virtual DirectionType direction (double s) const;

  PositionTypeDouble positionInDouble (double s) const;

  DirectionTypeDouble directionInDouble (double s) const;

private:
  template <class T> 
  std::pair<bool,double> genericPathLength (const T& object) const;

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

  HelixExtrapolatorToLine2Order theQuadraticSolutionFromStart;

  const PropagationDirection thePropDir;

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

#endif