Point in space, time and velocity. More...

#include <gparticle.h>

Public Member Functions
	stvpoint ()=default
	Default constructor.

	stvpoint (const point &fpt, const vec &fdir, vfloat fspeed, manip_absvol_treeid &ftid, vfloat fprange, vfloat ftime, int fsb, int fs_ent, manip_absvol *faeid)
	Constructor.

	stvpoint (const stvpoint &pstv, const trajestep &ts, vfloat mrange, int fsb, int fs_ent, manip_absvol *faeid)

	stvpoint (const stvpoint &pstv, const trajestep &ts, int fsb, int fs_ent, manip_absvol *faeid)

	stvpoint (const stvpoint &fp)=default
	Copy constructor.

void	print (std::ostream &file, int l) const

absvol *	volume ()

Public Attributes
point	pt
	Coordinates in the first system in the tree.

vec	dir
	Unit vector, in the first system in the tree.

point	ptloc
	Coordinates in the local system (last system in the tree).

vec	dirloc
	Unit vector, in the local system (last system in the tree).

vfloat	speed = 0.
	Longitudinal velocity.

manip_absvol_treeid	tid

int	sb = 0

int	s_ent = 0

manip_absvol *	next_eid = nullptr

vfloat	prange = 0.
	Range from previous point.

vfloat	time = 0.

Detailed Description

Point in space, time and velocity.

Definition at line 20 of file gparticle.h.

Constructor & Destructor Documentation

◆ stvpoint() [1/5]

Heed::stvpoint::stvpoint ( )

default

Default constructor.

◆ stvpoint() [2/5]

Heed::stvpoint::stvpoint	(	const point &	fpt,
		const vec &	fdir,
		vfloat	fspeed,
		manip_absvol_treeid &	ftid,
		vfloat	fprange,
		vfloat	ftime,
		int	fsb,
		int	fs_ent,
		manip_absvol *	faeid
	)

inline

Constructor.

Definition at line 52 of file gparticle.h.

      : pt(fpt),
        dir(unit_vec(fdir)),
        speed(fspeed),
        tid(ftid),
        sb(fsb),
        s_ent(fs_ent),
        next_eid(faeid),
        prange(fprange),
        time(ftime) {
    ptloc = pt;
    tid.up_absref(&ptloc);
    dirloc = dir;
    tid.up_absref(&dirloc);
  }

◆ stvpoint() [3/5]

Heed::stvpoint::stvpoint	(	const stvpoint &	pstv,
		const trajestep &	ts,
		vfloat	mrange,
		int	fsb,
		int	fs_ent,
		manip_absvol *	faeid
	)

inline

Constructor.

Parameters

pstv	previous point
ts	trajectory step (in the local system)
mrange	step length (may be less than the one in ts)
fsb	position flag
fs_ent	"entering" flag
faeid	next volume

Definition at line 77 of file gparticle.h.

      : pt(),
        dir(),
        ptloc(),
        dirloc(),
        speed(pstv.speed),
        tid(pstv.tid),
        sb(fsb),
        s_ent(fs_ent),
        next_eid(faeid),
        prange(mrange) {
    if (pstv.speed == 0) {
      time = pstv.time;  // just to put anything
    } else {
      // If speed is changed, this time is to be corrected in derivative class
      time = pstv.time + mrange / pstv.speed;
    }
    ts.Gnextpoint(mrange, ptloc, dirloc);
    pt = ptloc;
    tid.down_absref(&pt);
    dir = dirloc;
    tid.down_absref(&dir);
  }

◆ stvpoint() [4/5]

Heed::stvpoint::stvpoint	(	const stvpoint &	pstv,
		const trajestep &	ts,
		int	fsb,
		int	fs_ent,
		manip_absvol *	faeid
	)

inline

Constructor.

Parameters

pstv	previous point
ts	trajectory step (in the local system)
fsb	position flag
fs_ent	"entering" flag
faeid	next volume

Definition at line 108 of file gparticle.h.

      : pt(),
        dir(),
        ptloc(),
        dirloc(),
        speed(pstv.speed),
        tid(pstv.tid),
        sb(fsb),
        s_ent(fs_ent),
        next_eid(faeid),
        prange(ts.mrange) {
    if (pstv.speed == 0) {
      time = pstv.time;  // just to put anything
    } else {
      // If speed is changed, this time is to be corrected in derivative class
      time = pstv.time + ts.mrange / pstv.speed;
    }
    ptloc = ts.mpoint;
    point temp;
    ts.Gnextpoint(ts.mrange, temp, dirloc);
    pt = ptloc;
    tid.down_absref(&pt);
    dir = dirloc;
    tid.down_absref(&dir);
  }

◆ stvpoint() [5/5]

Heed::stvpoint::stvpoint ( const stvpoint & fp )

default

Copy constructor.

Member Function Documentation

◆ print()

void Heed::stvpoint::print	(	std::ostream &	file,
		int	l
	)		const

Definition at line 16 of file gparticle.cpp.

                                                  {
  if (l < 0) return;
  Ifile << "stvpoint: sb=" << sb << " s_ent=" << s_ent << " prange=" << prange
        << " time=" << time << '\n';
  indn.n += 2;
  Ifile << "position:\n" << pt << ptloc;
  Ifile << "direction:\n" << dir << dirloc;
  Ifile << "speed=" << speed << '\n';
  if (tid.eid.empty()) {
    Ifile << "point is outside universe\n";
    file.flush();
    indn.n -= 2;
    return;
  }
  tid.print(file, 1);
  char s[100];
  if (sb == 2) {
    next_eid->m_chname(s);
    Ifile << "next volume name " << s << '\n';
  }
  indn.n -= 2;
  file.flush();
}