Heed::HeedPhoton Class Reference

#include <HeedPhoton.h>

Inheritance diagram for Heed::HeedPhoton:

Public Member Functions
	HeedPhoton (manip_absvol *primvol, const point &pt, const vec &vel, vfloat time, long fparent_particle_number, double fenergy, int fs_print_listing=0)
void	physics_after_new_speed (void)
virtual void	physics (void)
	HeedPhoton (void)
virtual void	print (std::ostream &file, int l) const
	macro_copy_total (HeedPhoton)
virtual	~HeedPhoton ()
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 Attributes
long	particle_number
long	parent_particle_number
int	s_print_listing
double	energy
int	s_photon_absorbed
long	na_absorbing
long	ns_absorbing
int	s_delta_generated
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

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

Detailed Description

Definition at line 17 of file HeedPhoton.h.

Constructor & Destructor Documentation

HeedPhoton() [1/2]

Heed::HeedPhoton::HeedPhoton	(	manip_absvol *	primvol,
		const point &	pt,
		const vec &	vel,
		vfloat	time,
		long	fparent_particle_number,
		double	fenergy,
		int	fs_print_listing = 0
	)

Definition at line 18 of file HeedPhoton.cpp.

    : gparticle(primvol, pt, vel, time),
      parent_particle_number(fparent_particle_number),
      s_print_listing(fs_print_listing),
      energy(fenergy),
      s_photon_absorbed(0),
#ifdef SFER_PHOTOEL
      s_sfer_photoel(0),
#endif
      s_delta_generated(0) {
  mfunname("HeedPhoton::HeedPhoton(...)");
  particle_number = last_particle_number;
  //set_count_references();
  last_particle_number++;
  double length_vel = length(vel);
  check_econd11(fabs(length_vel - c_light) / (length_vel + c_light), > 1.0e-10,
                mcerr);
}

HeedPhoton() [2/2]

Heed::HeedPhoton::HeedPhoton ( void  )

inline

Definition at line 39 of file HeedPhoton.h.

: gparticle() { ; }

Referenced by physics_after_new_speed().

~HeedPhoton()

virtual Heed::HeedPhoton::~HeedPhoton ( )

inlinevirtual

Definition at line 42 of file HeedPhoton.h.

{ ; }

Member Function Documentation

macro_copy_total()

Heed::HeedPhoton::macro_copy_total

(

HeedPhoton

)

physics()

void Heed::HeedPhoton::physics ( void  )

virtual

Reimplemented from gparticle.

Definition at line 39 of file HeedPhoton.cpp.

                             {
  mfunname("void HeedPhoton::physics(void)");
 
  if (s_print_listing == 1) {
    mcout << "void HeedPhoton::physics(void) is starting\n";
  }
  if (s_photon_absorbed == 0)  // still not absorbed
      {
    if (nextpos.prange > 0.0) {
      const absvol* av = currpos.G_lavol();  // get least address of volume
      HeedMatterDef* hmd = NULL;
      const EnTransfCSType* etcst = dynamic_cast<const EnTransfCSType*>(av);
      if (etcst != NULL) {
        hmd = etcst->etcs->hmd.getver();
      } else {
        const HeedDeltaElectronCSType* hmecst =
            dynamic_cast<const HeedDeltaElectronCSType*>(av);
        if (hmecst != NULL) {
          hmd = hmecst->hdecs->hmd.get();
        }
      }
      if (hmd != NULL) {
        // calculate cross section
        // First to count shells
        long qst = 0;
        long qa = hmd->matter->qatom();
        long na;
        for (na = 0; na < qa; na++) {
          long qs = hmd->apacs[na]->get_qshell();
          qst += qs;
        }
        DynLinArr<double> cs(qst);
        DynLinArr<long> nat(qst);
        DynLinArr<long> nsh(qst);
        double s = 0.0;
        long nst = 0;
        for (na = 0; na < qa; na++) {
          long qs = hmd->apacs[na]->get_qshell();
          long ns;
          double at_weight_quan = hmd->matter->weight_quan(na);
          for (ns = 0; ns < qs; ns++) {
            cs[nst] = hmd->apacs[na]->get_ICS(ns, energy) * at_weight_quan;
            // threshold is taken into account in apacs[na]->get_ACS(ns,..
            nat[nst] = na;
            nsh[nst] = ns;
            s += cs[nst];
            //Imcout<<"na="<<na<<" ns="<<ns<<" cs[nst]="<<cs[nst]<<'\n';
            nst++;
          }
        }
        if (s_print_listing == 1) {
          Iprintn(mcout, s);
        }
        //s=s * hmd->eldens / hmd->matter->Z_mean() * C1_MEV_CM;
        s = s * 1.0e-18 * AVOGADRO / (hmd->matter->A_mean() / (gram / mole)) *
            hmd->matter->density() / (gram / cm3);
        if (s_print_listing == 1) {
          Iprintn(mcout, s);
        }
        double path_length = 1.0 / s;  // cm
        if (s_print_listing == 1) {
          Iprint2n(mcout, energy, path_length);
        }
        double xleng = -path_length * log(1.0 - SRANLUX());
        // computing randon length
        if (s_print_listing == 1) {
          Iprintn(mcout, xleng);
          Iprintn(mcout, nextpos.prange / cm);
        }
        if (xleng * cm < nextpos.prange) {
          s_photon_absorbed = 1;
#ifdef SFER_PHOTOEL
          s_sfer_photoel = 1;  // assumes that virtual photons are already
                               // absorbed and s_sfer_photoel is  0 for them
#endif
 
          //s_life = 0;
          chispre(cs);
          double r = chisran(SRANLUX(), cs);
          //Iprintn(mcout, r);
          long n = long(r);
          if (n < 0) n = 0;
          if (n > nst - 1) n = nst - 1;
          if (s_print_listing == 1) {
            Iprintn(mcout, n);
          }
          na_absorbing = nat[n];
          ns_absorbing = nsh[n];
          nextpos.prange = xleng * cm;
          nextpos.pt = currpos.pt + nextpos.prange * currpos.dir;
          nextpos.ptloc = nextpos.pt;
          nextpos.tid.up_absref(&nextpos.ptloc);
        }
      }
    }
  }
}

physics_after_new_speed()

void Heed::HeedPhoton::physics_after_new_speed ( void  )

virtual

Reimplemented from gparticle.

Definition at line 137 of file HeedPhoton.cpp.

                                             {
  mfunname("void HeedPhoton::physics_after_new_speed(void)");
  if (s_print_listing == 1) {
    mcout << "HeedPhoton::physics_after_new_speed is started\n";
  }
  if (s_photon_absorbed == 1 && s_delta_generated == 0) {
    const absvol* av = currpos.G_lavol();  // get least address of volume
 
    HeedMatterDef* hmd = NULL;
    const EnTransfCSType* etcst = dynamic_cast<const EnTransfCSType*>(av);
    if (etcst != NULL) {
      hmd = etcst->etcs->hmd.getver();
    } else {
      const HeedDeltaElectronCSType* hmecst =
          dynamic_cast<const HeedDeltaElectronCSType*>(av);
      if (hmecst != NULL) {
        hmd = hmecst->hdecs->hmd.get();
      }
    }
    check_econd11(hmd, == NULL, mcerr);
    // generate delta-electrons
    DynLinArr<double> el_energy;
    DynLinArr<double> ph_energy;
    hmd->apacs[na_absorbing]
        ->get_escape_particles(ns_absorbing, energy, el_energy, ph_energy);
    if (s_print_listing == 1) {
      mcout << "The condition:\n";
      Iprint2n(mcout, na_absorbing, ns_absorbing);
      mcout << "The decay products:\n";
      Iprintn(mcout, el_energy);
      Iprintn(mcout, ph_energy);
    }
    long qel = el_energy.get_qel();
    long nel;
    for (nel = 0; nel < qel; nel++) {
      vec vel;
      if (nel == 0) {  // assumed it is photoelectron
#ifdef SFER_PHOTOEL
        if (s_sfer_photoel == 1) {
          vel.random_sfer_vec();
        } else {
          vel = currpos.dir;
        }
#else
        vel = currpos.dir;  // direction is OK
#endif
      } else {
        vel.random_sfer_vec();
      }
      double gam_1 = el_energy[nel] / ELMAS;
      double betta = sqrt(1.0 - pow(1.0 / (gam_1 + 1.0), 2.0));
      double mod_v = betta * c_light;
      vel = vel * mod_v;
      if (s_print_listing == 1) {
        mcout << "Initializing delta electron\n";
        Iprint4n(mcout, el_energy[nel], gam_1, betta, mod_v);
      }
      ActivePtr<gparticle> ac;
      ac.pass(
          new HeedDeltaElectron(currpos.tid.eid[0].amvol.getver(), currpos.pt,
                                vel, currpos.time, particle_number));
      particle_bank.insert_after(particle_bank.get_last_node(), ac);
    }
    long qph = ph_energy.get_qel();
    long nph;
    for (nph = 0; nph < qph; nph++) {
      vec vel;
      vel.random_sfer_vec();
      vel *= c_light;
      if (s_print_listing == 1) {
        mcout << "Initializing photon\n";
        Iprint2n(mcout, el_energy[nph], vel);
      }
      ActivePtr<gparticle> ac;
      ac.pass(new HeedPhoton(currpos.tid.eid[0].amvol.getver(), currpos.pt, vel,
                             currpos.time, particle_number, ph_energy[nph]));
      particle_bank.insert_after(particle_bank.get_last_node(), ac);
    }
    s_delta_generated = 1;
    s_life = 0;
  }
  if (s_print_listing == 1) {
    mcout << "HeedPhoton::physics_after_new_speed is exitted\n";
  }
}

print()

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

virtual

Reimplemented from gparticle.

Definition at line 223 of file HeedPhoton.cpp.

                                                    {
  if (l >= 0) {
    Ifile << "HeedPhoton (l=" << l << "): particle_number=" << particle_number
          << " energy=" << energy << "MeV\n";
    if (l == 1) return;
    indn.n += 2;
    Ifile << "s_photon_absorbed=" << s_photon_absorbed
          << " na_absorbing=" << na_absorbing
          << " ns_absorbing=" << ns_absorbing
          << " s_delta_generated=" << s_delta_generated
#ifdef SFER_PHOTOEL
          << " s_sfer_photoel=" << s_sfer_photoel
#endif
          << " parent_particle_number=" << parent_particle_number
          << " s_print_listing=" << s_print_listing << '\n';
    gparticle::print(file, l - 1);
    indn.n -= 2;
  }
}

Member Data Documentation

energy

double Heed::HeedPhoton::energy

Definition at line 24 of file HeedPhoton.h.

Referenced by Garfield::TrackHeed::GetCluster(), physics(), physics_after_new_speed(), and print().

na_absorbing

long Heed::HeedPhoton::na_absorbing

Definition at line 26 of file HeedPhoton.h.

Referenced by Heed::HeedParticle::physics(), Heed::HeedParticle_BGM::physics(), physics(), physics_after_new_speed(), and print().

ns_absorbing

long Heed::HeedPhoton::ns_absorbing

Definition at line 27 of file HeedPhoton.h.

Referenced by Heed::HeedParticle::physics(), Heed::HeedParticle_BGM::physics(), physics(), physics_after_new_speed(), and print().

parent_particle_number

long Heed::HeedPhoton::parent_particle_number

Definition at line 20 of file HeedPhoton.h.

Referenced by Garfield::TrackHeed::GetCluster(), print(), and Garfield::TrackHeed::TransportPhoton().

particle_number

long Heed::HeedPhoton::particle_number

Definition at line 19 of file HeedPhoton.h.

Referenced by Garfield::TrackHeed::GetCluster(), HeedPhoton(), physics_after_new_speed(), print(), and Garfield::TrackHeed::TransportPhoton().

s_delta_generated

int Heed::HeedPhoton::s_delta_generated

Definition at line 32 of file HeedPhoton.h.

Referenced by Heed::HeedParticle::physics(), Heed::HeedParticle_BGM::physics(), physics_after_new_speed(), and print().

s_photon_absorbed

int Heed::HeedPhoton::s_photon_absorbed

Definition at line 25 of file HeedPhoton.h.

Referenced by Heed::HeedParticle::physics(), Heed::HeedParticle_BGM::physics(), physics(), physics_after_new_speed(), and print().

s_print_listing

int Heed::HeedPhoton::s_print_listing

Definition at line 21 of file HeedPhoton.h.

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

---

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

• HeedPhoton.h
• HeedPhoton.cpp