Heed::trajestep Class Reference

#include <trajestep.h>

Inheritance diagram for Heed::trajestep:

Public Member Functions
	trajestep (const vfloat fmax_range, const vfloat frad_for_straight, const vfloat fmax_straight_arange, const vfloat fmax_circ_arange, const point &fcurrpos, const vec &fdir, const bool fcurved, const vec &frelcen, vfloat fmrange, vfloat prec)
	trajestep (const trajestep &fts, vfloat fmrange)
	trajestep ()=default
	Default constructor.
virtual	~trajestep ()
	Destructor.
void	Gnextpoint (vfloat frange, point &fpos, vec &fdir) const
	Move to the next point.
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
vfloat	max_range = 100. * CLHEP::cm
	Max. step length.
vfloat	rad_for_straight = 1000. * CLHEP::cm
	Radius beyond which to prefer straight lines to reduce calculation time.
vfloat	max_straight_arange = 0.1 * CLHEP::rad
	Angular step for curved lines when using straight-line approximation.
vfloat	max_circ_arange = 0.2 * CLHEP::rad
	Angular step for curved lines.
point	currpos
	Current position.
vec	dir
	Unit vector.
bool	curved = false
	Type of trajectory (curved or straight).
vec	relcen
int	s_range_cf = 0
int	s_prec = 0
vfloat	mrange = 0.
	Maximal possible range.
point	mpoint

Protected Member Functions
virtual absref_transmit	get_components () override

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

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 32 of file trajestep.h.

Constructor & Destructor Documentation

trajestep() [1/3]

Heed::trajestep::trajestep	(	const vfloat	fmax_range,
		const vfloat	frad_for_straight,
		const vfloat	fmax_straight_arange,
		const vfloat	fmax_circ_arange,
		const point &	fcurrpos,
		const vec &	fdir,
		const bool	fcurved,
		const vec &	frelcen,
		vfloat	fmrange,
		vfloat	prec
	)

Constructor.

Parameters

fmax_range	maximum step length
frad_for_straight	radius beyond which to use straight-line steps.
fmax_straight_arange	angular step for straight-line approximation.
fmax_circ_arange	angular step for curved steps.
fcurrpos	initial coordinates.
fdir	initial direction.
fcurved	flag whether the trajectory is curved or straight.
frelcen	centre of rotation (only used for curved lines).
fmrange	can be used for reducing/limiting the step length.
prec	tolerance for checking if frelcen is perpendicular to dir.

Definition at line 25 of file trajestep.cpp.

    : max_range(fmax_range),
      rad_for_straight(frad_for_straight),
      max_straight_arange(fmax_straight_arange),
      max_circ_arange(fmax_circ_arange),
      currpos(fcurrpos),
      dir(),
      curved(fcurved),
      relcen(frelcen),
      s_prec(1),
      mrange(fmrange) {
  pvecerror("trajestep::trajestep(...)");
  if (fdir == dv0) {
    dir = dv0;
    mrange = 0;
  } else {
    dir = unit_vec(fdir);
    if (curved) {
      check_econd11a(check_perp(dir, relcen, prec), != 1,
                     "dir=" << dir << "relcen=" << relcen
                            << "fcurrpos=" << fcurrpos << "fdir=" << fdir,
                     mcerr);
    }
    if (mrange < 0 || mrange > max_range) mrange = max_range;
    s_range_cf = curved;
    if (s_range_cf == 1) {
      const vfloat r = relcen.length();
      if (r >= rad_for_straight) {
        s_range_cf = 0;
        mrange = std::min(mrange, r * max_straight_arange);
      } else {
        mrange = std::min(mrange, r * max_circ_arange);
      }
    }
  }
}

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 66 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.max_range, fts.rad_for_straight,
                fts.max_straight_arange, fts.max_circ_arange,
                fpos, fdir, fts.curved, frelcen, fmrange, prec);
}

trajestep() [3/3]

Heed::trajestep::trajestep ( )

default

Default constructor.

Referenced by trajestep().

~trajestep()

virtual Heed::trajestep::~trajestep ( )

inlinevirtual

Destructor.

Definition at line 58 of file trajestep.h.

{}

Member Function Documentation

get_components()

absref_transmit Heed::trajestep::get_components ( )

overrideprotectedvirtual

Reimplemented from Heed::absref.

Definition at line 21 of file trajestep.cpp.

                                          {
  return absref_transmit(4, aref);
}

Gnextpoint()

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

Move to the next point.

Definition at line 80 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 (!curved) {
      // no curvature
      fdir = dir;
    } else {
      vfloat ang = frange / relcen.length();
      fdir = dir;
      fdir.turn(dir || relcen, ang);
    }
  } 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;
  }
}

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

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 102 of file trajestep.h.

Referenced by get_components().

currpos

point Heed::trajestep::currpos

Current position.

Definition at line 74 of file trajestep.h.

Referenced by Gnextpoint(), Heed::operator<<(), Heed::splane::range(), Heed::box::range_ext(), and Heed::ulsvolume::range_ext().

curved

bool Heed::trajestep::curved = false

Type of trajectory (curved or straight).

Definition at line 78 of file trajestep.h.

Referenced by Gnextpoint(), Heed::operator<<(), and trajestep().

dir

vec Heed::trajestep::dir

Unit vector.

Definition at line 76 of file trajestep.h.

Referenced by Gnextpoint(), Heed::operator<<(), Heed::splane::range(), Heed::ulsvolume::range_ext(), and trajestep().

max_circ_arange

vfloat Heed::trajestep::max_circ_arange = 0.2 * CLHEP::rad

Angular step for curved lines.

Definition at line 71 of file trajestep.h.

Referenced by trajestep().

max_range

vfloat Heed::trajestep::max_range = 100. * CLHEP::cm

Max. step length.

Definition at line 64 of file trajestep.h.

Referenced by trajestep().

max_straight_arange

vfloat Heed::trajestep::max_straight_arange = 0.1 * CLHEP::rad

Angular step for curved lines when using straight-line approximation.

Definition at line 69 of file trajestep.h.

Referenced by trajestep().

mpoint

point Heed::trajestep::mpoint

Definition at line 98 of file trajestep.h.

Referenced by Heed::operator<<(), Heed::ulsvolume::range_ext(), and Heed::stvpoint::stvpoint().

mrange

vfloat Heed::trajestep::mrange = 0.

Maximal possible range.

Definition at line 93 of file trajestep.h.

Referenced by Heed::gparticle::calc_step_to_bord(), Gnextpoint(), Heed::operator<<(), Heed::splane::range(), Heed::box::range_ext(), Heed::ulsvolume::range_ext(), Heed::stvpoint::stvpoint(), and trajestep().

rad_for_straight

vfloat Heed::trajestep::rad_for_straight = 1000. * CLHEP::cm

Radius beyond which to prefer straight lines to reduce calculation time.

Definition at line 67 of file trajestep.h.

Referenced by trajestep().

relcen

vec Heed::trajestep::relcen

Centre of rotation relative to currpos. Used only for curved trajectories. If used, should be perpendicular to dir.

Definition at line 83 of file trajestep.h.

Referenced by Gnextpoint(), Heed::operator<<(), Heed::splane::range(), and trajestep().

s_prec

int Heed::trajestep::s_prec = 0

Definition at line 90 of file trajestep.h.

Referenced by Heed::gparticle::calc_step_to_bord(), Heed::operator<<(), and Heed::ulsvolume::range_ext().

s_range_cf

int Heed::trajestep::s_range_cf = 0

Definition at line 87 of file trajestep.h.

Referenced by Gnextpoint(), Heed::operator<<(), Heed::splane::range(), and trajestep().

---

The documentation for this class was generated from the following files:

• trajestep.h
• trajestep.cpp