#include <trajestep.h>

Inheritance diagram for Heed::trajestep:

Public Member Functions
void	Gnextpoint (vfloat frange, point &fpos, vec &fdir) const

void	Gnextpoint1 (vfloat frange, point &fpos, vec &fdir, vec &frelcen) const

	trajestep (trajestep_limit *ftl, const point &fcurrpos, const vec &fdir, int fs_cf, const vec &frelcen, vfloat fmrange, vfloat prec)

	trajestep (const trajestep &fts, vfloat fmrange)

	trajestep ()
	Default constructor.

virtual	~trajestep ()
	Destructor.

Public Member Functions inherited from Heed::absref
virtual	~absref ()
	Destructor.

virtual void	down (const abssyscoor *fasc)
	Convert numbering representation of object to basical system of fasc.

virtual void	up (const abssyscoor *fasc)
	Convert numbering representation of objects to new system.

virtual void	turn (const vec &dir, vfloat angle)
	Turn around axis doing via center of coordinate system along dir.

virtual void	shift (const vec &dir)

Public Attributes
PassivePtr< trajestep_limit >	tl

point	currpos

vec	dir
	Unit vector.

int	s_cf

vec	relcen

int	s_range_cf

int	s_prec

vfloat	mrange
	Maximal possible range.

point	mpoint

Protected Member Functions
virtual void	get_components (ActivePtr< absref_transmit > &aref_tran)

Static Protected Attributes
static absrefabsref::*[4]	aref

Detailed Description

Trajectory step of any object (particle, photon, ...). Here we interested in geometrical parameters only. The time, speed, acceleration, mass, and forces are not interesting here. The object can move by straight line or by curved line. The real trajectory is approximated by little steps. Each step can be straight or with constant curvature (circumference). If the path is always straight (for example for light or neutral particles), the step length is limited by extra conditions, typically edge of a volume. If the path is curved, the step length is limited first by length at which the curvature is changed, then by precision of approximation of the real shape of curvature.

Definition at line 65 of file trajestep.h.

Constructor & Destructor Documentation

◆ trajestep() [1/3]

Heed::trajestep::trajestep	(	trajestep_limit *	ftl,
		const point &	fcurrpos,
		const vec &	fdir,
		int	fs_cf,
		const vec &	frelcen,
		vfloat	fmrange,
		vfloat	prec
	)

Constructor. Here prec is used to check if frelcen is perp. to dir. If it is not perpendicular with this precision, the function terminates the program. To reduce range fmrange may be used.

Definition at line 38 of file trajestep.cpp.

    : tl(ftl),
      currpos(fcurrpos),
      dir(),
      s_cf(fs_cf),
      relcen(frelcen),
      s_prec(1),
      mrange(fmrange) {
  pvecerror("trajestep::trajestep(...)");
  if (fdir == dv0) {
    dir = dv0;
    mrange = 0;
  } else {
    dir = unit_vec(fdir);
    if (s_cf == 1) {
      check_econd11a(check_perp(dir, relcen, prec), != 1,
                     "dir=" << dir << "relcen=" << relcen
                            << "fcurrpos=" << fcurrpos << "fdir=" << fdir,
                     mcerr);
    }
    tl->range(s_cf, relcen.length(), s_range_cf, mrange);
  }
}

◆ trajestep() [2/3]

Heed::trajestep::trajestep	(	const trajestep &	fts,
		vfloat	fmrange
	)

Constructor to continue propagation from the end point of another step.

Parameters

fts	old step to continue
fmrange	new range to travel

Definition at line 64 of file trajestep.cpp.

                                                         {
  mfunname("trajestep::trajestep(const trajestep& fts, vfloat fmrange)");
  // Continue propagation from the end point of the old step.
  point fpos;
  vec fdir;
  vec frelcen;
  fts.Gnextpoint1(fts.mrange, fpos, fdir, frelcen);
  vfloat prec = 0.1;  // not important here
  *this =
      trajestep(fts.tl.getver(), fpos, fdir, fts.s_cf, frelcen, fmrange, prec);
}

◆ trajestep() [3/3]

Heed::trajestep::trajestep ( )

inline

Default constructor.

Definition at line 113 of file trajestep.h.

              : 
        tl(),
        currpos(),
        dir(),
        s_cf(0),
        relcen(),
        s_range_cf(0),
        s_prec(0),
        mrange(0),
        mpoint() {}

Referenced by trajestep().

◆ ~trajestep()

virtual Heed::trajestep::~trajestep ( )

inlinevirtual

Destructor.

Definition at line 124 of file trajestep.h.

124{}

Member Function Documentation

◆ get_components()

void Heed::trajestep::get_components ( ActivePtr< absref_transmit > & aref_tran )

protectedvirtual

Reimplemented from Heed::absref.

Definition at line 34 of file trajestep.cpp.

                                                                    {
  aref_tran.pass(new absref_transmit(4, aref));
}

◆ Gnextpoint()

void Heed::trajestep::Gnextpoint	(	vfloat	frange,
		point &	fpos,
		vec &	fdir
	)		const

Definition at line 76 of file trajestep.cpp.

                                                                      {
  pvecerror("int trajestep::Gnextpoint(vfloat frange, point& fpos, vec& fdir)");
  check_econd12(frange, >, mrange, mcerr);
  if (s_range_cf == 0) {
    // interpolation by straight line
    fpos = currpos + frange * dir;
    if (s_cf == 0) {
      // no curvature
      fdir = dir;
      return;
    } else {
      vfloat ang = frange / relcen.length();
      fdir = dir;
      fdir.turn(dir || relcen, ang);
      return;
    }
  } else {
    vfloat ang = frange / relcen.length();  // angle to turn
    fdir = dir;
    fdir.turn(dir || relcen, ang);  // direction at the end
    vec frelcen = relcen;
    frelcen.turn(dir || relcen, ang);
    fpos = currpos + relcen - frelcen;
    return;
  }
}

Referenced by Heed::stvpoint::stvpoint().

◆ Gnextpoint1()

void Heed::trajestep::Gnextpoint1	(	vfloat	frange,
		point &	fpos,
		vec &	fdir,
		vec &	frelcen
	)		const

Definition at line 103 of file trajestep.cpp.

                                                {
  pvecerror(
      "int trajestep::Gnextpoint(vfloat frange, point& fpos, vec& fdir, "
      "vec& frelcen)");
  check_econd12(frange, >, mrange, mcerr);
  if (s_range_cf == 0) {
    // interpolation by straight line
    fpos = currpos + frange * dir;
    if (s_cf == 0) {
      // no curvature
      fdir = dir;
      frelcen = relcen;  // whatever it is
      return;
    } else {
      vfloat ang = frange / relcen.length();
      fdir = dir;
      fdir.turn(dir || relcen, ang);
      frelcen = relcen;
      frelcen.turn(dir || relcen, ang);
      return;
    }
  } else {
    vfloat ang = frange / relcen.length();  // angle to turn
    fdir = dir;
    fdir.turn(dir || relcen, ang);  // direction at the end
    frelcen = relcen;
    frelcen.turn(dir || relcen, ang);
    fpos = currpos + relcen - frelcen;
    return;
  }
}

Referenced by trajestep().

Member Data Documentation

◆ aref

absref absref::* Heed::trajestep::aref

staticprotected

Initial value:

= {(absref absref::*)&trajestep::currpos,
                                        (absref absref::*)&trajestep::dir,
                                        (absref absref::*)&trajestep::relcen,
                                        (absref absref::*)&trajestep::mpoint}

Definition at line 128 of file trajestep.h.

Referenced by get_components().