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(double const x,
             double const y,
             double const z,
             double const dirx,
             double const diry,
             double const dirz);

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

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

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

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

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

  public:
    //
    // Template
    //

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

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

    template <class T>
    void Dir(T const& 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