#include <eiparticle.h>

Inheritance diagram for Heed::eiparticle:

Public Member Functions
	eiparticle (void)
	Constructors.

	eiparticle (manip_absvol primvol, const point &pt, const vec &vel, vfloat time, particle_def fpardef, int fs_add_loss=0)

	AnyType_copy (eiparticle, gparticle)

virtual	~eiparticle ()
	Destructor.

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

double	Bethe_Bloch_en_loss (void)
	Energy loss per unit length.

virtual void	physics_after_new_speed (void)

Public Member Functions inherited from Heed::eparticle
	eparticle (void)
	Constructors.

	eparticle (manip_absvol primvol, const point &pt, const vec &vel, vfloat time, particle_def fpardef)

	eparticle (manip_absvol primvol, const point &pt, const vec &vel, vfloat time, particle_def fpardef, double gamma_1)

	AnyType_copy (eparticle, gparticle)

virtual	~eparticle ()
	Destructor.

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

virtual int	force (const point &pt, vec &f, vec &f_perp, vfloat &mrange)

Public Member Functions inherited from Heed::mparticle
void	check_consistency (void) const

virtual void	step (void)

virtual void	curvature (int &fs_cf, vec &frelcen, vfloat &fmrange, vfloat prec)

virtual void	physics_after_new_speed (void)

virtual void	physics (void)

virtual int	force (const point &pt, vec &f, vec &f_perp, vfloat &mrange)

void	new_speed (void)

	mparticle (void)

	mparticle (gparticle const &gp, double fmass)

	mparticle (gparticle const &gp, double fmass, double gamma_1)

	mparticle (manip_absvol *primvol, const point &pt, const vec &vel, vfloat time, double fmass, double gamma_1)

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

	macro_copy_total (gparticle)

virtual	~mparticle ()

Public Member Functions inherited from gparticle
	gparticle (void)

	gparticle (const stvpoint &sp)

	gparticle (manip_absvol *primvol, const point &pt, const vec &vel, vfloat time)

virtual void	step (void)

virtual void	change_vol (void)

virtual void	curvature (int &fs_cf, vec &frelcen, vfloat &fmrange, vfloat prec)

virtual void	physics_after_new_speed (void)

virtual void	physics (void)

virtual void	physics_mrange (double &fmrange)

virtual stvpoint	calc_step_to_bord ()

stvpoint	switch_new_vol (void)

virtual void	fly (void)

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

	macro_copy_total (gparticle)

virtual	~gparticle ()

Public Member Functions inherited from Heed::particle_type
	particle_type (void)

	particle_type (particle_def *f)

	particle_type (const char *name, int s=0)

int	operator== (const particle_type &f)

int	operator!= (const particle_type &f)

void	print_notation (std::ostream &file) const

Public Attributes
int	s_add_loss

double	total_loss

Public Attributes inherited from Heed::mparticle
double	mass
	Mass (not mass * speed_of_light^2)

double	orig_kin_energy

double	orig_gamma_1

double	prev_kin_energy

double	prev_gamma_1

double	curr_kin_energy

double	curr_gamma_1

Public Attributes inherited from gparticle
int	s_life

long	nstep

double	total_range_from_origin

long	n_zero_step

stvpoint	origin

stvpoint	prevpos

stvpoint	currpos

stvpoint	nextpos

vec	curr_relcen

Public Attributes inherited from Heed::particle_type
PassivePtr< particle_def >	pardef

Additional Inherited Members
Static Public Attributes inherited from gparticle
static long	max_q_zero_step = 100

Detailed Description

Definition at line 16 of file eiparticle.h.

Constructor & Destructor Documentation

◆ eiparticle() [1/2]

Heed::eiparticle::eiparticle ( void )

inline

Constructors.

Definition at line 19 of file eiparticle.h.

19: eparticle() { ; }

Heed::eparticle::eparticle

eparticle(void)

Constructors.

Definition: eparticle.h:19

◆ eiparticle() [2/2]

Heed::eiparticle::eiparticle	(	manip_absvol *	primvol,
		const point &	pt,
		const vec &	vel,
		vfloat	time,
		particle_def *	fpardef,
		int	fs_add_loss = `0`
	)

inline

Definition at line 20 of file eiparticle.h.

      : eparticle(primvol, pt, vel, time, fpardef),
        s_add_loss(fs_add_loss),
        total_loss(0.0) {}

◆ ~eiparticle()

virtual Heed::eiparticle::~eiparticle ( )

inlinevirtual

Destructor.

Definition at line 27 of file eiparticle.h.

27{}

Member Function Documentation

◆ AnyType_copy()

Heed::eiparticle::AnyType_copy	(	eiparticle	,
		gparticle
	)

◆ Bethe_Bloch_en_loss()

double Heed::eiparticle::Bethe_Bloch_en_loss ( void )

Energy loss per unit length.

Definition at line 14 of file eiparticle.cpp.

                                       {
 
  mfunname("double eiparticle::Bethe_Bloch_energy_loss()");
  // Get least address of volume
  const absvol* av = currpos.G_lavol();
  const MatterType* amt = dynamic_cast<const MatterType*>(av);
  if (amt == NULL) return 0.;
  MatterDef* amd = amt->matdef.get();
  const double beta = lorbeta(curr_gamma_1);
  const double loss =
      Bethe_Bloch_energy_loss(amd->Z_mean() / amd->A_mean(), amd->I_eff(), beta,
                              pardef->charge / eplus);
  return loss * amd->density();
}

Referenced by physics_after_new_speed().

◆ physics_after_new_speed()

void Heed::eiparticle::physics_after_new_speed ( void )

virtual

Reimplemented from Heed::mparticle.

Definition at line 29 of file eiparticle.cpp.

                                         {
  mfunname("void eiparticle::physics_after_new_speed(void)");
  double loss = Bethe_Bloch_en_loss();
  loss *= currpos.prange;
  total_loss += loss;
  if (s_add_loss == 0) {
    curr_kin_energy -= loss;
    if (curr_kin_energy <= 0.0) {
      curr_kin_energy = 0.0;
      curr_gamma_1 = 0.0;
      currpos.speed = 0.0;
      s_life = 0;
    } else {
      const double resten = mass * c_squared;
      curr_gamma_1 = curr_kin_energy / resten;
      currpos.speed = c_light * lorbeta(curr_gamma_1);
    }
  } else {
    curr_kin_energy += loss;
    const double resten = mass * c_squared;
    curr_gamma_1 = curr_kin_energy / resten;
    currpos.speed = c_light * lorbeta(curr_gamma_1);
  }
}

◆ print()

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

virtual

Reimplemented from Heed::eparticle.

Definition at line 54 of file eiparticle.cpp.

                                                    {
  if (l >= 0) {
    Ifile << "eiparticle: s_add_loss=" << s_add_loss
          << " total_loss/keV=" << total_loss / keV << '\n';
    eparticle::print(file, l);
  }
}

Member Data Documentation

◆ s_add_loss

int Heed::eiparticle::s_add_loss

Flag for handling energy loss 0 - loss is subtracted from current particle energy 1 - loss is added (simulation of particle passage to opposite

Definition at line 39 of file eiparticle.h.

Referenced by physics_after_new_speed(), and print().

◆ total_loss

double Heed::eiparticle::total_loss

Definition at line 40 of file eiparticle.h.