/data/refman/pasoursint/CMSSW_4_1_8_patch13/src/TrackingTools/GeomPropagators/interface/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"

class HelixArbitraryPlaneCrossing2Order : 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
  virtual ~HelixArbitraryPlaneCrossing2Order() {}

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

  virtual PositionType position(double s) const;

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

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

};

#endif