Line.hh

#ifndef Line_hh
#define Line_hh

#include <utility>
using std::pair;

#include "StThreeVectorD.hh"

class Line {
public:
  Line() {}
  Line(const StThreeVectorD& o, const StThreeVectorD& d);

  StThreeVectorD origin() const;

  StThreeVectorD direction() const;

  StThreeVectorD at(double pathlength) const;

  StThreeVectorD perigee(const StThreeVectorD& point) const;

  StThreeVectorD dca(const Line& line) const;

  double pathlength(const StThreeVectorD& point) const;

  pair<double, double> pathlengths(const Line& line) const;

  void setOrigin(const StThreeVectorD& o);

  void setDirection(const StThreeVectorD& d);

private:
  StThreeVectorD mOrigin;
  StThreeVectorD mDirection;
};

inline StThreeVectorD Line::origin() const { return mOrigin; }
inline StThreeVectorD Line::direction() const { return mDirection; }
inline void Line::setOrigin(const StThreeVectorD& o) { mOrigin = o; }
inline void Line::setDirection(const StThreeVectorD& d) { mDirection = d.unit(); }

#endif