GeoHalfLine.h

Go to the documentation of this file.

#ifndef BASICTOOL_GEOHALFLINE_H
#define BASICTOOL_GEOHALFLINE_H

#include "larcorealg/GeoAlgo/GeoAlgoException.h"
#include "larcorealg/GeoAlgo/GeoVector.h"

#include "TVector3.h"

namespace geoalgo {
  class HalfLine {

  public:
    HalfLine();

    virtual ~HalfLine(){};

    HalfLine(const double x,
             const double y,
             const double z,
             const double dirx,
             const double diry,
             const double dirz);

    HalfLine(const Point_t& start, const Vector_t& dir);

    const Point_t& Start() const; 
    const Vector_t& Dir() const;  

    void Start(const double x, const double y, const double z); 
    void Dir(const double x, const double y, const double z);   

    void Start(const TVector3& pt); 
    void Dir(const TVector3& dir);  

  protected:
    void Normalize(); 
    Point_t _start;   
    Vector_t _dir;    

  public:
    //
    // Template
    //

    template <class T, class U>
    HalfLine(const T& start, const U& dir) : HalfLine(Point_t(start), Vector_t(dir))
    {}

    template <class T>
    void Start(const T& pos)
    {
      _start = Point_t(pos);
      if (_start.size() != 3)
        throw GeoAlgoException("<<Start>> Only 3 dimensional start point allowed!");
    }

    template <class T>
    void Dir(const T& dir)
    {
      _dir = Vector_t(dir);
      if (_dir.size() != 3)
        throw GeoAlgoException("<<Start>> Only 3 dimensional start point allowed!");
      Normalize();
    }
  };

  typedef HalfLine HalfLine_t;
}
#endif
 // end of doxygen group