p_raytr_hits.H

/* p_raytr_hits.H */
/* $Id: p_raytr_hits.H,v 2.7 2003/10/01 17:40:05 fern Exp $ */
#ifndef __p_raytr_hits_H
#define __p_raytr_hits_H

#include <mix/wenv.h>
#include <mix/ObjectsAr.H>
#include <mix/geometry.H>
#include <mth/simplex.H>

class ray;
class ray_parameters;
// ºÛÐááë ÔÛï ãÝØäØÚÐæØØ ßàÞæÕÔãàë ÝÐáâàÕÛÐ
// ºÞÛØçÕáâÒÞ ØâÕàÐæØÙ Ò àÐÙÞÝÕ _ÝãÛÕÒÞÓÞ_ ãÓÛÐ
//#define MN_CHINA_D 5
#define MN_CHINA_D 0      // fern
// _ÝãÛÕÒÞÙ_ ãÓÞÛ
#define MN_PHI_D 0.03

// ºÞÛØçÕáâÒÞ ÝÕãÔÐçÝëå ßÞßëâÞÚ ßÞÔàïÔ ßàØÑÛØ×Øâìáï Ú ßàØÕÜÝØÚã
#define MILNUM2_D 10

#define PM_ALPHA_STOP_D 0.0001
#define PM_MAX_ITERATIONS_COUNT_D 50

#define PM_ALPHA_LESS_D 0.7
#define PM_ALPHA_MORE_D 1.3

// ¾ÑéÐï çÐáâì ÒëáâàÕÛÐ ÔÛï ÒáÕå ÝÐáâàÕÛÞÒ
class general_hit_info{
protected:
  double dLess,dMore;
  int mn_china;
  double mn_phi;
  int milnum2;
  double alpha_stop;
  int max_iterations;
  
public:
  general_hit_info(double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
           int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
           double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);

  virtual ray*  shoot(ray_parameters&,point&,Logic&,double);
  virtual ray*  shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);

  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const crossing &cClosest,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);

  // ÄãÝÚæØØ ÔÛï Ø×ÞÑàÐÖÕÝØï ×ÐßØåÐÝÝëå ÚÞÝáâÐÝâ
  double lesser_alpha(){return dLess;};
  double increase_alpha(){return dMore;};
  int n_of_iterations(){return max_iterations;};
  double alpha_min(){return alpha_stop;};
  int mn2(){return milnum2;};
  double phi(){return mn_phi;};
  double iphi(){return mn_china;};
};

class hit_point_info:public general_hit_info{
protected:
  double dfl,dal;
public:
  hit_point_info(double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
         int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
         double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D)
    :general_hit_info(dl,dm,i1,d1,i2,d2,i3){dfl=0.15;dal=0.1;};
  
  virtual void
  get_ray_parameters(const point3d &pWhere,const crossing &cClosest,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
};

class straight_hit_info:public hit_point_info{
protected:
  rays_shoot* pr;
  int iBody;
  double dAcc;
  double dALast;
  IntAr iaTemplate;

  Logic verify_ray(ray*)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  straight_hit_info(rays_shoot*,int,double,
            ray*,int,// Template
            double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
            int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
            double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
};


class reflected_hit_info:public hit_point_info{
protected:
  IntAr iaTemplate;
  rays_shoot* pr;
  int iBody;
  double dAcc;
  double dALast;
  int iBorder;
  Logic verify_ray(ray*)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  reflected_hit_info(rays_shoot*,int,int,double,
             ray*,int,// Template
             double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
             int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
             double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
};

// ¿ÞßÐÔÐÝØÕ Ò âÞçÚã áØÜßÛÕÚá-ÜÕâÞÔÞÜ
class ray_hit_simplex:public simplex2D{
protected:
  general_hit_info& gs;
  point pWhere;
  IntAr iaTemplate;
  Logic verify_ray(ray*)const;
public:
  virtual double value(const point2D&);
  ray_hit_simplex(general_hit_info& ghiShot,const point&pW,
          const point2D&p1,const point2D&p2,const double&d1,
          const double&d2,
          ray*,// Template
          double le=DEFAULT_LESSER_VALUE,double ee=DEFAULT_ENLARGE_VALUE,
          double _a=DEFAULT_A,double _b=DEFAULT_B,double _g=DEFAULT_G);
};

#endif /* p_raytr_hits.H */

/*  Melnikov George Yu. 24-Dec-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 
/*          ...          */ 


// ºÞÛØçÕáâÒÞ ØâÕàÐæØÙ Ò àÐÙÞÝÕ _ÝãÛÕÒÞÓÞ_ ãÓÛÐ//#define MN_CHINA_D 5
#define MN_CHINA_D 0      // fern
// _ÝãÛÕÒÞÙ_ ãÓÞÛ
#define MN_PHI_D 0.03

// ºÞÛØçÕáâÒÞ ÝÕãÔÐçÝëå ßÞßëâÞÚ ßÞÔàïÔ ßàØÑÛØ×Øâìáï Ú ßàØÕÜÝØÚã
#define MILNUM2_D 10

#define PM_ALPHA_STOP_D 0.0001
#define PM_MAX_ITERATIONS_COUNT_D 50

#define PM_ALPHA_LESS_D 0.7
#define PM_ALPHA_MORE_D 1.3

// ¾ÑéÐï çÐáâì ÒëáâàÕÛÐ ÔÛï ÒáÕå ÝÐáâàÕÛÞÒ
class general_hit_info{
protected:
  double dLess,dMore;
  int mn_china;
  double mn_phi;
  int milnum2;
  double alpha_stop;
  int max_iterations;
  
public:
  general_hit_info(double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
           int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
           double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);

  virtual ray*  shoot(ray_parameters&,point&,Logic&,double);
  virtual ray*  shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);

  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const crossing &cClosest,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);

  // ÄãÝÚæØØ ÔÛï Ø×ÞÑàÐÖÕÝØï ×ÐßØåÐÝÝëå ÚÞÝáâÐÝâ
  double lesser_alpha(){return dLess;};
  double increase_alpha(){return dMore;};
  int n_of_iterations(){return max_iterations;};
  double alpha_min(){return alpha_stop;};
  int mn2(){return milnum2;};
  double phi(){return mn_phi;};
  double iphi(){return mn_china;};
};

class hit_point_info:public general_hit_info{
protected:
  double dfl,dal;
public:
  hit_point_info(double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
         int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
         double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D)
    :general_hit_info(dl,dm,i1,d1,i2,d2,i3){dfl=0.15;dal=0.1;};
  
  virtual void
  get_ray_parameters(const point3d &pWhere,const crossing &cClosest,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
};

class straight_hit_info:public hit_point_info{
protected:
  rays_shoot* pr;
  int iBody;
  double dAcc;
  double dALast;
  IntAr iaTemplate;

  Logic verify_ray(ray*)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  straight_hit_info(rays_shoot*,int,double,
            ray*,int,// Template
            double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
            int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
            double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
};


class reflected_hit_info:public hit_point_info{
protected:
  IntAr iaTemplate;
  rays_shoot* pr;
  int iBody;
  double dAcc;
  double dALast;
  int iBorder;
  Logic verify_ray(ray*)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  reflected_hit_info(rays_shoot*,int,int,double,
             ray*,int,// Template
             double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
             int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
             double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
};

// ¿ÞßÐÔÐÝØÕ Ò âÞçÚã áØÜßÛÕÚá-ÜÕâÞÔÞÜ
class ray_hit_simplex:public simplex2D{
protected:
  general_hit_info& gs;
  point pWhere;
  IntAr iaTemplate;
  Logic verify_ray(ray*)const;
public:
  virtual double value(const point2D&);
  ray_hit_simplex(general_hit_info& ghiShot,const point&pW,
          const point2D&p1,const point2D&p2,const double&d1,
          const double&d2,
          ray*,// Template
          double le=DEFAULT_LESSER_VALUE,double ee=DEFAULT_ENLARGE_VALUE,
          double _a=DEFAULT_A,double _b=DEFAULT_B,double _g=DEFAULT_G);
};

#endif /* p_raytr_hits.H */

/*  Melnikov George Yu. 24-Dec-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 
/*          ...          */ 

//#define MN_CHINA_D 5
#define MN_CHINA_D 0      // fern
// _ÝãÛÕÒÞÙ_ ãÓÞÛ#define MN_PHI_D 0.03

// ºÞÛØçÕáâÒÞ ÝÕãÔÐçÝëå ßÞßëâÞÚ ßÞÔàïÔ ßàØÑÛØ×Øâìáï Ú ßàØÕÜÝØÚã
#define MILNUM2_D 10

#define PM_ALPHA_STOP_D 0.0001
#define PM_MAX_ITERATIONS_COUNT_D 50

#define PM_ALPHA_LESS_D 0.7
#define PM_ALPHA_MORE_D 1.3

// ¾ÑéÐï çÐáâì ÒëáâàÕÛÐ ÔÛï ÒáÕå ÝÐáâàÕÛÞÒ
class general_hit_info{
protected:
  double dLess,dMore;
  int mn_china;
  double mn_phi;
  int milnum2;
  double alpha_stop;
  int max_iterations;
  
public:
  general_hit_info(double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
           int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
           double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);

  virtual ray*  shoot(ray_parameters&,point&,Logic&,double);
  virtual ray*  shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);

  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const crossing &cClosest,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);

  // ÄãÝÚæØØ ÔÛï Ø×ÞÑàÐÖÕÝØï ×ÐßØåÐÝÝëå ÚÞÝáâÐÝâ
  double lesser_alpha(){return dLess;};
  double increase_alpha(){return dMore;};
  int n_of_iterations(){return max_iterations;};
  double alpha_min(){return alpha_stop;};
  int mn2(){return milnum2;};
  double phi(){return mn_phi;};
  double iphi(){return mn_china;};
};

class hit_point_info:public general_hit_info{
protected:
  double dfl,dal;
public:
  hit_point_info(double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
         int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
         double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D)
    :general_hit_info(dl,dm,i1,d1,i2,d2,i3){dfl=0.15;dal=0.1;};
  
  virtual void
  get_ray_parameters(const point3d &pWhere,const crossing &cClosest,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
};

class straight_hit_info:public hit_point_info{
protected:
  rays_shoot* pr;
  int iBody;
  double dAcc;
  double dALast;
  IntAr iaTemplate;

  Logic verify_ray(ray*)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  straight_hit_info(rays_shoot*,int,double,
            ray*,int,// Template
            double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
            int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
            double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
};


class reflected_hit_info:public hit_point_info{
protected:
  IntAr iaTemplate;
  rays_shoot* pr;
  int iBody;
  double dAcc;
  double dALast;
  int iBorder;
  Logic verify_ray(ray*)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  reflected_hit_info(rays_shoot*,int,int,double,
             ray*,int,// Template
             double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
             int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
             double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
};

// ¿ÞßÐÔÐÝØÕ Ò âÞçÚã áØÜßÛÕÚá-ÜÕâÞÔÞÜ
class ray_hit_simplex:public simplex2D{
protected:
  general_hit_info& gs;
  point pWhere;
  IntAr iaTemplate;
  Logic verify_ray(ray*)const;
public:
  virtual double value(const point2D&);
  ray_hit_simplex(general_hit_info& ghiShot,const point&pW,
          const point2D&p1,const point2D&p2,const double&d1,
          const double&d2,
          ray*,// Template
          double le=DEFAULT_LESSER_VALUE,double ee=DEFAULT_ENLARGE_VALUE,
          double _a=DEFAULT_A,double _b=DEFAULT_B,double _g=DEFAULT_G);
};

#endif /* p_raytr_hits.H */

/*  Melnikov George Yu. 24-Dec-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 
/*          ...          */ 

#define MN_PHI_D 0.03

// ºÞÛØçÕáâÒÞ ÝÕãÔÐçÝëå ßÞßëâÞÚ ßÞÔàïÔ ßàØÑÛØ×Øâìáï Ú ßàØÕÜÝØÚã#define MILNUM2_D 10

#define PM_ALPHA_STOP_D 0.0001
#define PM_MAX_ITERATIONS_COUNT_D 50

#define PM_ALPHA_LESS_D 0.7
#define PM_ALPHA_MORE_D 1.3

// ¾ÑéÐï çÐáâì ÒëáâàÕÛÐ ÔÛï ÒáÕå ÝÐáâàÕÛÞÒ
class general_hit_info{
protected:
  double dLess,dMore;
  int mn_china;
  double mn_phi;
  int milnum2;
  double alpha_stop;
  int max_iterations;
  
public:
  general_hit_info(double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
           int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
           double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);

  virtual ray*  shoot(ray_parameters&,point&,Logic&,double);
  virtual ray*  shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);

  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const crossing &cClosest,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);

  // ÄãÝÚæØØ ÔÛï Ø×ÞÑàÐÖÕÝØï ×ÐßØåÐÝÝëå ÚÞÝáâÐÝâ
  double lesser_alpha(){return dLess;};
  double increase_alpha(){return dMore;};
  int n_of_iterations(){return max_iterations;};
  double alpha_min(){return alpha_stop;};
  int mn2(){return milnum2;};
  double phi(){return mn_phi;};
  double iphi(){return mn_china;};
};

class hit_point_info:public general_hit_info{
protected:
  double dfl,dal;
public:
  hit_point_info(double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
         int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
         double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D)
    :general_hit_info(dl,dm,i1,d1,i2,d2,i3){dfl=0.15;dal=0.1;};
  
  virtual void
  get_ray_parameters(const point3d &pWhere,const crossing &cClosest,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
};

class straight_hit_info:public hit_point_info{
protected:
  rays_shoot* pr;
  int iBody;
  double dAcc;
  double dALast;
  IntAr iaTemplate;

  Logic verify_ray(ray*)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  straight_hit_info(rays_shoot*,int,double,
            ray*,int,// Template
            double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
            int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
            double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
};


class reflected_hit_info:public hit_point_info{
protected:
  IntAr iaTemplate;
  rays_shoot* pr;
  int iBody;
  double dAcc;
  double dALast;
  int iBorder;
  Logic verify_ray(ray*)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  reflected_hit_info(rays_shoot*,int,int,double,
             ray*,int,// Template
             double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
             int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
             double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
};

// ¿ÞßÐÔÐÝØÕ Ò âÞçÚã áØÜßÛÕÚá-ÜÕâÞÔÞÜ
class ray_hit_simplex:public simplex2D{
protected:
  general_hit_info& gs;
  point pWhere;
  IntAr iaTemplate;
  Logic verify_ray(ray*)const;
public:
  virtual double value(const point2D&);
  ray_hit_simplex(general_hit_info& ghiShot,const point&pW,
          const point2D&p1,const point2D&p2,const double&d1,
          const double&d2,
          ray*,// Template
          double le=DEFAULT_LESSER_VALUE,double ee=DEFAULT_ENLARGE_VALUE,
          double _a=DEFAULT_A,double _b=DEFAULT_B,double _g=DEFAULT_G);
};

#endif /* p_raytr_hits.H */

/*  Melnikov George Yu. 24-Dec-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 
/*          ...          */ 

#define MILNUM2_D 10

#define PM_ALPHA_STOP_D 0.0001
#define PM_MAX_ITERATIONS_COUNT_D 50

#define PM_ALPHA_LESS_D 0.7
#define PM_ALPHA_MORE_D 1.3

// ¾ÑéÐï çÐáâì ÒëáâàÕÛÐ ÔÛï ÒáÕå ÝÐáâàÕÛÞÒ
class general_hit_info{
protected:
  double dLess,dMore;
  int mn_china;
  double mn_phi;
  int milnum2;
  double alpha_stop;
  int max_iterations;
  
public:
  general_hit_info(double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
           int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
           double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);

  virtual ray*  shoot(ray_parameters&,point&,Logic&,double);
  virtual ray*  shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);

  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const crossing &cClosest,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);

  // ÄãÝÚæØØ ÔÛï Ø×ÞÑàÐÖÕÝØï ×ÐßØåÐÝÝëå ÚÞÝáâÐÝâ
  double lesser_alpha(){return dLess;};
  double increase_alpha(){return dMore;};
  int n_of_iterations(){return max_iterations;};
  double alpha_min(){return alpha_stop;};
  int mn2(){return milnum2;};
  double phi(){return mn_phi;};
  double iphi(){return mn_china;};
};

class hit_point_info:public general_hit_info{
protected:
  double dfl,dal;
public:
  hit_point_info(double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
         int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
         double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D)
    :general_hit_info(dl,dm,i1,d1,i2,d2,i3){dfl=0.15;dal=0.1;};
  
  virtual void
  get_ray_parameters(const point3d &pWhere,const crossing &cClosest,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
  virtual void
  get_ray_parameters(const point3d &pWhere,const ray* rLast,double &dFi,
             double &dAz,double &dDFi,double& dDAz,double &dRo);
};

class straight_hit_info:public hit_point_info{
protected:
  rays_shoot* pr;
  int iBody;
  double dAcc;
  double dALast;
  IntAr iaTemplate;

  Logic verify_ray(ray*)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  straight_hit_info(rays_shoot*,int,double,
            ray*,int,// Template
            double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
            int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
            double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
};


class reflected_hit_info:public hit_point_info{
protected:
  IntAr iaTemplate;
  rays_shoot* pr;
  int iBody;
  double dAcc;
  double dALast;
  int iBorder;
  Logic verify_ray(ray*)const;
public:
  virtual ray* shoot(ray_parameters&,point&,Logic&,double);
  virtual ray* shoot(double,double,point&,Logic&,double);
  virtual Logic good_enough(double,ray*);
  reflected_hit_info(rays_shoot*,int,int,double,
             ray*,int,// Template
             double dl=PM_ALPHA_LESS_D,double dm=PM_ALPHA_MORE_D,
             int i1=MN_CHINA_D,double d1=MN_PHI_D,int i2=MILNUM2_D,
             double d2=PM_ALPHA_STOP_D,int i3=PM_MAX_ITERATIONS_COUNT_D);
};

// ¿ÞßÐÔÐÝØÕ Ò âÞçÚã áØÜßÛÕÚá-ÜÕâÞÔÞÜ
class ray_hit_simplex:public simplex2D{
protected:
  general_hit_info& gs;
  point pWhere;
  IntAr iaTemplate;
  Logic verify_ray(ray*)const;
public:
  virtual double value(const point2D&);
  ray_hit_simplex(general_hit_info& ghiShot,const point&pW,
          const point2D&p1,const point2D&p2,const double&d1,
          const double&d2,
          ray*,// Template
          double le=DEFAULT_LESSER_VALUE,double ee=DEFAULT_ENLARGE_VALUE,
          double _a=DEFAULT_A,double _b=DEFAULT_B,double _g=DEFAULT_G);
};

#endif /* p_raytr_hits.H */

/*  Melnikov George Yu. 24-Dec-1999  */ 
/*          ...          */ 
/*  Melnikov George Yu.  1-Feb-2000  */ 
/*          ...          */