StraightLineCylinderCrossing.h

Go to the documentation of this file.

#ifndef StraightLineCylinderCrossing_H
#define StraightLineCylinderCrossing_H

#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/GeometryVector/interface/LocalVector.h"
#include "DataFormats/GeometryVector/interface/Basic2DVector.h"

#include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"

#include <utility>

class Cylinder;

class StraightLineCylinderCrossing {
public:
  StraightLineCylinderCrossing(const LocalPoint& startingPos,
                               const LocalVector& startingDir,
                               const PropagationDirection propDir = alongMomentum,
                               double tolerance = 0);

  std::pair<bool, double> pathLength(const Cylinder& cyl) const;

  LocalPoint position(const double s) const { return LocalPoint(theX0 + s * theP0); }

private:
  std::pair<bool, double> chooseSolution(const double s1, const double s2) const;

private:
  //
  // single precision is sufficient for intermediate vectors
  //
  typedef LocalPoint PositionType;
  typedef LocalVector DirectionType;
  typedef Basic2DVector<float> PositionType2D;
  typedef Basic2DVector<float> DirectionType2D;

  const PositionType theX0;
  const DirectionType theP0;
  const PropagationDirection thePropDir;
  double theTolerance;
};

#endif