HelixExtrapolatorToLine2Order.h

Go to the documentation of this file.

#ifndef HelixExtrapolatorToLine2Order_h_
#define HelixExtrapolatorToLine2Order_h_

#include "TrackingTools/GeomPropagators/interface/HelixLineExtrapolation.h"
#include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"
#include "FWCore/Utilities/interface/Visibility.h"

class HelixExtrapolatorToLine2Order final : public HelixLineExtrapolation {
public:
  HelixExtrapolatorToLine2Order(const PositionType& point,
                                const DirectionType& direction,
                                const float curvature,
                                const PropagationDirection propDir = alongMomentum);

  HelixExtrapolatorToLine2Order(const double& x0,
                                const double& y0,
                                const double& z0,
                                const double& cosPhi0,
                                const double& sinPhi0,
                                const double& cosTheta,
                                const double& sinTheta,
                                const double& rho,
                                const PropagationDirection propDir = alongMomentum)
      : thePosition(x0, y0, z0),
        theDirection(cosPhi0, sinPhi0, cosTheta / sinTheta),
        theSinTheta(sinTheta),
        theRho(rho),
        thePropDir(propDir) {}

  // destructor
  ~HelixExtrapolatorToLine2Order() override {}

  std::pair<bool, double> pathLength(const GlobalPoint& point) const override;

  std::pair<bool, double> pathLength(const Line& line) const override;

  PositionType position(double s) const override;

  DirectionType direction(double s) const override;

  PositionTypeDouble positionInDouble(double s) const;

  DirectionTypeDouble directionInDouble(double s) const;

private:
  virtual std::pair<bool, double> pathLengthFromCoefficients(const double ceq[4]) const dso_internal;
  int solve3rdOrder(const double ceq[], double sol[]) const dso_internal;
  int solve2ndOrder(const double ceq[], double sol[]) const dso_internal;

private:
  const PositionTypeDouble thePosition;
  DirectionTypeDouble theDirection;
  double theSinTheta;
  const double theRho;
  const PropagationDirection thePropDir;
};

#endif