Line2D.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2015  Belle II Collaboration                              *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Kirill Chilikin                                          *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
/* KLM headers. */
#include <klm/eklm/geometry/Circle2D.h>
#include <klm/eklm/geometry/Arc2D.h>
 
/* CLHEP headers. */
#include <CLHEP/Geometry/Point3D.h>
#include <CLHEP/Geometry/Vector3D.h>
 
namespace Belle2 {
  namespace EKLM {
 
    class Line2D {
 
    public:
 
      Line2D(double x, double y, double vecx, double vecy);
 
      ~Line2D();
 
      const HepGeom::Point3D<double>& getInitialPoint() const
      {
        return m_Point;
      }
 
      const HepGeom::Vector3D<double>& getVector() const
      {
        return m_Vector;
      }
 
      int findIntersection(const Line2D& line,
                           HepGeom::Point3D<double>* intersection) const;
 
      int findIntersection(const Circle2D& circle,
                           HepGeom::Point3D<double> intersections[2]) const;
 
      int findIntersection(const Arc2D& arc,
                           HepGeom::Point3D<double> intersections[2]) const;
 
      int findIntersection(const Line2D& line,
                           HepGeom::Point3D<double>* intersection,
                           double t[2]) const;
 
    protected:
 
      int findIntersection(const Circle2D& circle,
                           HepGeom::Point3D<double> intersections[2],
                           double t[2], double angles[2]) const;
 
      int selectIntersections(HepGeom::Point3D<double>* intersections,
                              bool* condition, int n) const;
 
      HepGeom::Point3D<double> m_Point;
 
      HepGeom::Vector3D<double> m_Vector;
 
    };
 
  }
 
}