#ifndef CIRCLE_H_SEEN
#define CIRCLE_H_SEEN

#include <math.h>
#include <iostream.h>
#include "vector3d.h"
#include "unitvector.h"
#include "utilgeom.h"

class Circle 
{
  
 public:

  Circle();
  Circle(double theta, double phi, double aperture);
  Circle(double x, double y, double z, double aperture);
  Circle(const Vector3d& v, double aperture=M_PI/2.);
  Circle(const Circle& c);
  virtual ~Circle() {}

  void SetCircle(const UnitVector&, double);
  void SetSpinAxis(double theta, double phi);
  void SetSpinAxis(const Vector3d&);
  void SetSpinAxis(double x, double y, double z);
  void SetApertureAngle(double aperture);
  void SetApertureAngle(const Circle&);
  
  // psi contient les 4 valeurs des angles d'intersection. -1 si les cercles ne se coupent pas
  // voir la numerotation dans le .cc
  bool Intersection(const Circle&, double* psi) const;

  // donne le UnitVector correspondant a une position donnee sur le cercle
  UnitVector ConvCircleSphere(double psi) const;

  // donne le UnitVector correspondant la tangente au cercle en une position donnee sur le cercle
  UnitVector TanOnCircle(double psi) const;

  // donne le vecteur tangent dans le plan (xy) a la sphere en une position donnee sur le cercle
  UnitVector EPhi(double psi) const;

  // donne l'autre vecteur tangent (orthogonal a EPhi)
  UnitVector ETheta(double psi) const;

   // donne l'angle de separation dans [0,2Pi] en une position donnee sur le cercle et EPhi
  double SepAngleTanEPhi02PI(double psi) const;

  double Theta() const {return _spinaxis.Theta();}
  double Phi() const {return _spinaxis.Phi();}
  double ApertureAngle() const {return _angouv;}
  Vector3d Omega() const {return _spinaxis;}
  virtual void Print(ostream&) const;
  Circle& operator=(const Circle&);
  bool operator==(const Circle&) const;
  bool operator!=(const Circle&) const;

 private:

  Vector3d _spinaxis;
  UnitVector _spinunitaxis;
  double _angouv;
  double _theta,_phi;
  double _x,_y,_z;

};

inline ostream& operator<<(ostream& s, const Circle& c) 
{  
  c.Print(s);  
  return s;  
}

#endif