In This Package:

Line.h

Go to the documentation of this file.

// $Id: Line.h,v 1.8 2008/06/24 12:24:27 jpalac Exp $
// ============================================================================
#ifndef LHCBMATH_LINE_H 
#define LHCBMATH_LINE_H 1
// ============================================================================
// Include files
// ============================================================================
#include <iostream>
// ============================================================================
namespace Gaudi
{
  namespace Math
  {
    
    // ==========================================================================
    template<typename aPoint, typename aVector>
    class Line 
    {
    public:
      typedef aPoint  Point;
      typedef aVector Vector;   
    public:
      Line() {}
      Line ( const aPoint& p0 , const aVector& v0 ) : m_p0 ( p0 ) , m_v0 ( v0      ) {}
      Line ( const aPoint& p0 , const aPoint&  p1 ) : m_p0 ( p0 ) , m_v0 ( p1 - p0 ) {}
      const aPoint&  beginPoint() const { return m_p0 ; }
      const aVector& direction()  const { return m_v0 ; }
      aPoint position ( const double mu ) const 
      { return beginPoint() + direction() * mu ; }
      aPoint operator() ( const double mu ) const 
      { return beginPoint() + direction() * mu ; }
    public:
      inline std::ostream& fillStream ( std::ostream& os ) const 
      {
        os << "\np0 (" 
           << m_p0.x() << " " << m_p0.y() << " " << m_p0.z()
           << ") direction (" 
           << m_v0.x()<< " " << m_v0.y() << " " << m_v0.z() << ")\n"
           << std::endl;
        return os;
      }
    private:
      aPoint  m_p0; // the start point on the line 
      aVector m_v0; // the direction vector of the line
    };

    // ========================================================================
  } // end of namespace Math

  // ==========================================================================
} // end of namespace Gaudi

// ============================================================================
template<typename aPoint, typename aVector>
inline std::ostream& operator<<
  (std::ostream&                      os  , 
   const Gaudi::Math::Line<aPoint,aVector>& rhs ) 
{ 
  return rhs.fillStream(os); 
}
// ==========================================================================
// The END 
// ==========================================================================
#endif // LHCBMATH_LINE_H
// ==========================================================================

---