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"
#include "FWCore/Utilities/interface/Visibility.h"

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

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

  PositionType position(double s) const override;

  DirectionType direction(double s) const override;
  //
  // 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 dso_internal;

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