HelixArbitraryPlaneCrossing2Order.h

Go to the documentation of this file.

#ifndef HELIXARBITRARYPLANECROSSING2ORDER_H_
#define HELIXARBITRARYPLANECROSSING2ORDER_H_
#include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"
#include "TrackingTools/GeomPropagators/interface/HelixPlaneCrossing.h"
#include "FWCore/Utilities/interface/Visibility.h"

class HelixArbitraryPlaneCrossing2Order final : public HelixPlaneCrossing {
public:
  HelixArbitraryPlaneCrossing2Order(const PositionType& point,
                                    const DirectionType& direction,
                                    const float curvature,
                                    const PropagationDirection propDir = alongMomentum);
  HelixArbitraryPlaneCrossing2Order(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)
      : theX0(x0),
        theY0(y0),
        theZ0(z0),
        theCosPhi0(cosPhi0),
        theSinPhi0(sinPhi0),
        theCosTheta(cosTheta),
        theSinThetaI(1. / sinTheta),
        theRho(rho),
        thePropDir(propDir) {}

  // destructor
  ~HelixArbitraryPlaneCrossing2Order() override {}

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

  PositionType position(double s) const override;

  DirectionType direction(double s) const override;
  //
  // double precision vectors
  //
  typedef Basic3DVector<double> PositionTypeDouble;
  typedef Basic3DVector<double> DirectionTypeDouble;

  PositionTypeDouble positionInDouble(double s) const;

  DirectionTypeDouble directionInDouble(double s) const;

  inline double smallestPathLength(const double firstPathLength, const double secondPathLength) const {
    return fabs(firstPathLength) < fabs(secondPathLength) ? firstPathLength : secondPathLength;
  }

private:
  std::pair<bool, double> solutionByDirection(const double dS1, const double dS2) const dso_internal;

private:
  const double theX0, theY0, theZ0;
  double theCosPhi0, theSinPhi0;
  double theCosTheta, theSinThetaI;
  const double theRho;
  const PropagationDirection thePropDir;
};

#endif