/*! 
  \brief Abstract class for flat surface representation, 

  defined by a sequence
  of point going arround the surface.
*/


#ifndef PNDDRCSURFPOLYFLAT_H
#define PNDDRCSURFPOLYFLAT_H



#include "PndDrcSurfAbs.h"
#include <vector>
using std::vector;
#include "Math/Transform3D.h"
using ROOT::Math::Transform3D;

class PndDrcSurfPolyFlat : public PndDrcSurfAbs
{
 public:

  PndDrcSurfPolyFlat();                    //!< Empty constructor.

  /*! \brief Copy constructor.
    \param s Object to copy.
  */
  PndDrcSurfPolyFlat(const PndDrcSurfPolyFlat& s);

  /*! \brief Assignment operator.
    \param s Object to assign.
  */
  PndDrcSurfPolyFlat& operator=(const PndDrcSurfPolyFlat& s);

  virtual ~PndDrcSurfPolyFlat(){};         //!< Destructor.

  /*! \brief Add a surface boundary point
    \param point The point
  */
  void     AddPoint(XYZPoint point);       

  /*! \brief Test if point within surface.

  \param point Point to test.
  \return true if within surface.
  */
  bool WithinSurface(XYZPoint& point) const;

  // implemetation of pure virtual functions.
  virtual  PndDrcSurfPolyFlat* Clone() const;
  XYZVector Normal(const XYZPoint& point)   const; 
  bool     SurfaceHit(PndDrcPhoton& ph, 
		      XYZPoint&  pos_new, 
		      double&    path_length) const;
  void     Print(fstream& stream)          const; 
  void     Print() const;
  void     AddTransform(const Transform3D& trans);
  virtual bool IsFlat() const {return true;};

  /*! \brief Set radius

  for usage in derived class
  \param r Radius
  */
  void SetRadius(double r){fRadius=r;};

  /*! \brief Set conical constant

  for usage in derived class
  The parameter k determines the shape. For
  \verbatim
  hyperbola k < -1 
  parabola  k = -1 
  ellipse   -1 < k < 0
  sphere    k = 0 
  ellipse   k > 0 
  \endverbatim
  \param c Conical constant
  */
  void SetConicalConstant(double c){fConConst=c;};

  /*! Point limiting the sphere.

  for usage in derived class

  The spheres construction points are not the same as the limiting points 
  of the sphere itself. This function returns the limiting points.
    \param i Index of point starting with 0
    \return The point.
  */
  virtual XYZPoint LimitingPoint(unsigned int i){return XYZPoint(0,0,0);};
 
 protected:

  vector<XYZPoint> fP;                  //!<  Vector of surface boundary points.             
  double fRadius;                       //!< Radius for curved surfaces.
  double fConConst;                     //!< Conial constant for curved surfaces.

 private:

  XYZVector        fNormal;             //!< Normal vector.           

  /*! \brief Auxiliary function for assignment operator and copy constructor..
    \param s The object to copy.
  */
  void Copy(const PndDrcSurfPolyFlat& s);         

  /*! \brief Line intersection

  between two lines each define by two points. The intersection has to be
  between the points.
  \param p1a 1st point of 1st line
  \param p1b 2nd point of 1st line
  \param p2a 1st point of 2nd line
  \param p2b 2nd point of 2nd line
  \return Flag if intersection exists.
  */
  bool LineCross(const XYZPoint& p1a, const XYZPoint& p1b,
		 const XYZPoint& p2a, const XYZPoint& p2b) const;



};
#endif