Heed::HeedParticle Class Reference

#include <HeedParticle.h>

Inheritance diagram for Heed::HeedParticle:

Public Member Functions
	HeedParticle (void)
	Constructors.
	HeedParticle (manip_absvol *primvol, const point &pt, const vec &vel, vfloat time, particle_def *fpardef, int fs_loss_only=0, int fs_print_listing=0)
	macro_copy_total (HeedParticle)
virtual	~HeedParticle ()
	Destructor.
virtual void	physics (void)
virtual void	print (std::ostream &file, int l) const
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_print_listing
long	particle_number
double	transferred_energy_in_step
long	qtransfer
int	s_loss_only
BlkArr< double >	transferred_energy
BlkArr< long >	natom
BlkArr< long >	nshell
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 30 of file HeedParticle.h.

Constructor & Destructor Documentation

HeedParticle() [1/2]

Heed::HeedParticle::HeedParticle ( void  )

inline

Constructors.

Definition at line 33 of file HeedParticle.h.

: eparticle() {}

HeedParticle() [2/2]

Heed::HeedParticle::HeedParticle	(	manip_absvol *	primvol,
		const point &	pt,
		const vec &	vel,
		vfloat	time,
		particle_def *	fpardef,
		int	fs_loss_only = 0,
		int	fs_print_listing = 0
	)

Definition at line 19 of file HeedParticle.cpp.

    : eparticle(primvol, pt, vel, time, fpardef),
      s_print_listing(fs_print_listing),
      transferred_energy_in_step(0.0),
      qtransfer(0),
      s_loss_only(fs_loss_only) {
 
  mfunname("HeedParticle::HeedParticle(...)");
  particle_number = last_particle_number;
  last_particle_number++;
  transferred_energy.allocate_block(100);
  natom.allocate_block(100);
  nshell.allocate_block(100);
 
}

~HeedParticle()

virtual Heed::HeedParticle::~HeedParticle ( )

inlinevirtual

Destructor.

Definition at line 42 of file HeedParticle.h.

{}

Member Function Documentation

macro_copy_total()

Heed::HeedParticle::macro_copy_total

(

HeedParticle

)

physics()

void Heed::HeedParticle::physics ( void  )

virtual

Reimplemented from Heed::mparticle.

Definition at line 37 of file HeedParticle.cpp.

                               {
  mfunname("void HeedParticle::physics(void)");
  if (s_print_listing == 1) {
    mcout << "HeedParticle::physics is started\n";
    Iprintn(mcout, currpos.prange);
  }
  transferred_energy_in_step = 0.0;
  qtransfer = 0;
  transferred_energy.allocate_block(100);
  natom.allocate_block(100);
  nshell.allocate_block(100);
  if (currpos.prange <= 0.0) return;
  // Get least address of volume
  const absvol* av = currpos.G_lavol();
  const EnTransfCSType* etcst = dynamic_cast<const EnTransfCSType*>(av);
  // Check if dynamic cast was successful.
  if (etcst == NULL) return;
 
  EnTransfCS* aetcs = etcst->etcs.getver();
  HeedMatterDef* ahmd = aetcs->hmd.getver();
  MatterDef* amatter = ahmd->matter.getver();
  EnergyMesh* a_energy_mesh = ahmd->energy_mesh.getver();
  const double* aetemp = ahmd->energy_mesh->get_ae();
  PointCoorMesh<double, const double*> pcm_e(a_energy_mesh->get_q() + 1,
                                             &(aetemp));
  long qa = amatter->qatom();
  if (s_print_listing == 1) Iprintn(mcout, qa);
  basis tempbas(currpos.dir, "tempbas");
  for (long na = 0; na < qa; ++na) {
    if (s_print_listing == 1) Iprintn(mcout, na);
    long qs = ahmd->apacs[na]->get_qshell();
    for (long ns = 0; ns < qs; ++ns) {
      if (s_print_listing == 1) Iprintn(mcout, ns);
      long qt = 0;
#ifdef SINGLE_TRANSFER
      if (aetcs == aetcs_single_transf && na == na_single_transf &&
          ns == ns_single_transf) {
        qt = 1;
      }
#else
      if (aetcs->quan[na][ns] > 0.0) {
        int ierror;
        qt = pois(aetcs->quan[na][ns] * currpos.prange / cm, ierror);
        check_econd11a(
            ierror, == 1,
            " aetcs->quan[na][ns]=" << aetcs->quan[na][ns]
                                    << " currpos.prange/cm=" << currpos.prange /
                                                                    cm << '\n',
            mcerr);
      }
#endif
      if (s_print_listing == 1) Iprintn(mcout, qt);
      if (qt > 0) {
        point curpt = prevpos.pt;
        vec dir = unit_vec(currpos.pt - prevpos.pt);
        double range = length(currpos.pt - prevpos.pt);
        if (s_print_listing == 1) {
          Iprint(mcout, curpt);
          Iprint(mcout, dir);
          Iprintn(mcout, range);
        }
        for (long nt = 0; nt < qt; ++nt) {
#ifdef SINGLE_TRANSFER
          transferred_energy.append(ener_single_transf);
#else
          double rn = SRANLUX();
          if (s_print_listing == 1) {
            Iprintn(mcout, rn);
            Iprintn(mcout, aetcs);
            Iprintn(mcout, aetcs->fadda[na][ns][1]);
          }
          double r = t_hisran_step_ar<double, DynLinArr<double>,
                                      PointCoorMesh<double, const double*> >(
              pcm_e, aetcs->fadda[na][ns], rn);
 
          // Convert to internal units.
          transferred_energy.append(r * MeV);
#endif
          if (s_print_listing == 1) {
            Iprint2n(mcout, nt, transferred_energy[qtransfer]);
          }
          transferred_energy_in_step += transferred_energy[qtransfer];
          natom.append(na);
          nshell.append(ns);
#ifdef SINGLE_TRANSFER
          double arange = 0.5 * range;
#else
          double arange = SRANLUX() * range;
#endif
          point pt = curpt + dir * arange;
          point ptloc = pt;
          prevpos.tid.up_absref(&ptloc);
          qtransfer++;
          if (s_loss_only == 0) {
            if (s_print_listing == 1) {
              mcout << "generating new cluster\n";
            }
            cluster_bank.append(HeedCluster(transferred_energy[qtransfer - 1],
                                            0, pt, ptloc, prevpos.tid, na, ns));
 
            double Ep0 = mass * c_squared + curr_kin_energy;
            double Ep1 = Ep0 - transferred_energy[qtransfer - 1];
            double Mp = mass;
            double Mt = electron_def.mass;
            double theta_p, theta_t;
            theta_two_part(Ep0, Ep1, Mp, Mt, theta_p, theta_t);
            vec vel;
            vel.random_conic_vec(fabs(theta_t));
            vel.down(&tempbas);  // direction is OK
            vel *= c_light;
            // HS
            double speed = length(vel);
            double time = arange / speed;
            if (s_print_listing == 1) {
              mcout << "generating new virtual photon\n";
            }
            HeedPhoton hp(currpos.tid.eid[0].amvol.getver(), pt, vel, time,
                          particle_number, transferred_energy[qtransfer - 1],
                          0);
            hp.s_photon_absorbed = 1;
            hp.s_delta_generated = 0;
            hp.na_absorbing = na;
            hp.ns_absorbing = ns;
            ActivePtr<gparticle> ac;
            ac.put(&hp);
            particle_bank.insert_after(particle_bank.get_last_node(), ac);
          }
        }
      }
    }
  }
  if (s_print_listing == 1) {
    Iprint2n(mcout, qtransfer, transferred_energy_in_step);
    mcout << "Exiting HeedParticle::physics\n";
  }
}

print()

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

virtual

Reimplemented from Heed::eparticle.

Definition at line 174 of file HeedParticle.cpp.

                                                      {
  if (l >= 0) {
    Ifile << "HeedParticle (l=" << l << "): particle_number=" << particle_number
          << " type=";
    print_notation(file);
    file << std::endl;
    if (l == 1) return;
    mparticle::print(file, l - 1);
    Iprintn(mcout, transferred_energy_in_step);
    Iprintn(mcout, qtransfer);
    if (l >= 5) {
      Ifile << "   nt  natom nshell transferred_energy\n";
      for (long nt = 0; nt < qtransfer; nt++) {
        Ifile << std::setw(3) << nt << ' ' << std::setw(3) << natom[nt] << ' '
              << std::setw(3) << nshell[nt] << ' ' << std::setw(12)
              << transferred_energy[nt] << '\n';
      }
    }
  }
 
}

Member Data Documentation

natom

BlkArr<long> Heed::HeedParticle::natom

Definition at line 54 of file HeedParticle.h.

Referenced by HeedParticle(), physics(), and print().

nshell

BlkArr<long> Heed::HeedParticle::nshell

Definition at line 55 of file HeedParticle.h.

Referenced by HeedParticle(), physics(), and print().

particle_number

long Heed::HeedParticle::particle_number

Definition at line 48 of file HeedParticle.h.

Referenced by HeedParticle(), physics(), and print().

qtransfer

long Heed::HeedParticle::qtransfer

Definition at line 51 of file HeedParticle.h.

Referenced by physics(), and print().

s_loss_only

int Heed::HeedParticle::s_loss_only

Definition at line 52 of file HeedParticle.h.

Referenced by physics().

s_print_listing

int Heed::HeedParticle::s_print_listing

Definition at line 47 of file HeedParticle.h.

Referenced by physics().

transferred_energy

BlkArr<double> Heed::HeedParticle::transferred_energy

Definition at line 53 of file HeedParticle.h.

Referenced by HeedParticle(), physics(), and print().

transferred_energy_in_step

double Heed::HeedParticle::transferred_energy_in_step

Definition at line 50 of file HeedParticle.h.

Referenced by physics(), and print().

---

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

• HeedParticle.h
• HeedParticle.cpp