Heed::AtomicSecondaryProducts Class Reference

#include <PhotoAbsCS.h>

Public Member Functions
	AtomicSecondaryProducts ()
	Constructor.
	~AtomicSecondaryProducts ()
	Destructor.
int	get_channel (std::vector< double > &felectron_energy, std::vector< double > &fphoton_energy) const
void	add_channel (double fchannel_prob_dens, const std::vector< double > &felectron_energy, const std::vector< double > &fphoton_energy, int s_all_rest=0)
void	print (std::ostream &file, int l) const

Protected Attributes
std::vector< double >	channel_prob_dens
std::vector< std::vector< double > >	electron_energy
std::vector< std::vector< double > >	photon_energy

Detailed Description

Sampling of secondary particles. The initial photo-electron is not included. The photo-electron is completely independent on the secondary channel and it always has an energy equal to transferred energy minus shell energy. This class proposes particles which can be emitted at filling this shell. This class can be interpreted as an additional channels for the default one. The sum of channel_prob_dens may be less than unity. The concrete channel for the particular event is sampled by get_channel(). If the random sampling does not point to one of the channels registering in this class, get_channel returns 0, which must be interpreted as the use of standard channel.

Definition at line 250 of file PhotoAbsCS.h.

Constructor & Destructor Documentation

AtomicSecondaryProducts()

Heed::AtomicSecondaryProducts::AtomicSecondaryProducts ( )

inline

Constructor.

Definition at line 253 of file PhotoAbsCS.h.

      : channel_prob_dens(), electron_energy(), photon_energy() {}

~AtomicSecondaryProducts()

Heed::AtomicSecondaryProducts::~AtomicSecondaryProducts ( )

inline

Destructor.

Definition at line 256 of file PhotoAbsCS.h.

{}

Member Function Documentation

add_channel()

void Heed::AtomicSecondaryProducts::add_channel	(	double	fchannel_prob_dens,
		const std::vector< double > &	felectron_energy,
		const std::vector< double > &	fphoton_energy,
		int	s_all_rest = 0
	)

Add new decay channel. Should be used at initialization.

Parameters

fchannel_prob_dens	probability for this channel.
felectron_energy	electron energies [MeV]
fphoton_energy	photon energies [MeV]
s_all_rest	if 1, the probability of this channel is assigned to that what is left to 1. fchannel_prob_dens is then ignored, it can be just 0.

Definition at line 469 of file PhotoAbsCS.cpp.

                                                             {
  mfunnamep("void AtomicSecondaryProducts::add_channel(...)");
  check_econd21(fchannel_prob_dens, < 0.0 ||, > 1.0, mcerr);
  long q_old = channel_prob_dens.size();
  long q_new = q_old + 1;
  channel_prob_dens.resize(q_new);
  electron_energy.resize(q_new);
  photon_energy.resize(q_new);
  if (s_all_rest == 1) {
    double s = 0.0;
    for (long n = 0; n < q_old; ++n) {
      s += channel_prob_dens[n];
    }
    check_econd21(s, < 0.0 ||, > 1.0, mcerr);
    fchannel_prob_dens = 1.0 - s;
  }
  channel_prob_dens[q_old] = fchannel_prob_dens;
  electron_energy[q_old] = felectron_energy;
  photon_energy[q_old] = fphoton_energy;
  double s = 0.0;
  for (long n = 0; n < q_new; ++n) {
    s += channel_prob_dens[n];
  }
  if (s > 1.0) {
    funnw.ehdr(mcerr);
    mcerr << "s > 1.0, s=" << s << '\n';
    Iprintn(mcerr, q_new);
    for (long n = 0; n < q_new; ++n) {
      mcerr << "n=" << n << " channel_prob_dens[n]=" << channel_prob_dens[n]
            << '\n';
    }
    spexit(mcerr);
  }
}

get_channel()

int Heed::AtomicSecondaryProducts::get_channel	(	std::vector< double > &	felectron_energy,
		std::vector< double > &	fphoton_energy
	)		const

Sample a channel (photon and electron energies in MeV). Return 1 if the channel is generated and 0 if not.

Definition at line 506 of file PhotoAbsCS.cpp.

          {
  mfunname("int AtomicSecondaryProducts::get_channel(...)");
#ifdef DEBUG_PRINT_get_escape_particles
  mcout << "AtomicSecondaryProducts::get_channel is started\n";
  Iprintn(mcout, channel_prob_dens.size());
#endif
  if (channel_prob_dens.empty()) return 0;
  int ir = 0;
  double rn = SRANLUX();
#ifdef DEBUG_PRINT_get_escape_particles
  Iprintn(mcout, rn);
#endif
  if (channel_prob_dens.size() == 1) {
    if (rn < channel_prob_dens[0]) {
      felectron_energy = electron_energy[0];
      fphoton_energy = photon_energy[0];
      ir = 1;
    }
  } else {
    long q = channel_prob_dens.size();
    double s = 0.0;
    for (long n = 0; n < q; ++n) {
      s += channel_prob_dens[n];
      if (rn <= s) {
        felectron_energy = electron_energy[n];
        fphoton_energy = photon_energy[n];
        ir = 1;
#ifdef DEBUG_PRINT_get_escape_particles
        Iprint2n(mcout, n, s);
#endif
        break;
      }
    }
  }
#ifdef DEBUG_PRINT_get_escape_particles
  mcout << "AtomicSecondaryProducts::get_channel is finishing\n";
  Iprintn(mcout, ir);
#endif
  return ir;
}

print()

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

Definition at line 549 of file PhotoAbsCS.cpp.

                                                                 {
  if (l <= 0) return;
  Ifile << "AtomicSecondaryProducts(l=" << l << "):\n";
  const long q = channel_prob_dens.size();
  Ifile << "number of channels=" << q << '\n';
  indn.n += 2;
  for (long n = 0; n < q; ++n) {
    Ifile << "n_channel=" << n << " probability=" << channel_prob_dens[n]
          << '\n';
    indn.n += 2;
    long qel = electron_energy[n].size();
    Ifile << "number of electrons=" << qel << '\n';
    indn.n += 2;
    for (long nel = 0; nel < qel; ++nel) {
      Ifile << "nel=" << nel << " electron_energy=" << electron_energy[n][nel]
            << '\n';
    }
    indn.n -= 2;
    long qph = photon_energy[n].size();
    Ifile << "number of photons=" << qph << '\n';
    indn.n += 2;
    for (long nph = 0; nph < qph; ++nph) {
      Ifile << "nph=" << nph << " photon_energy=" << photon_energy[n][nph]
            << '\n';
    }
    indn.n -= 2;
    indn.n -= 2;
  }
  indn.n -= 2;
}

Member Data Documentation

channel_prob_dens

std::vector<double> Heed::AtomicSecondaryProducts::channel_prob_dens

protected

Definition at line 279 of file PhotoAbsCS.h.

Referenced by add_channel(), get_channel(), and print().

electron_energy

std::vector<std::vector<double> > Heed::AtomicSecondaryProducts::electron_energy

protected

Definition at line 281 of file PhotoAbsCS.h.

Referenced by add_channel(), get_channel(), and print().

photon_energy

std::vector<std::vector<double> > Heed::AtomicSecondaryProducts::photon_energy

protected

Definition at line 282 of file PhotoAbsCS.h.

Referenced by add_channel(), get_channel(), and print().

---

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

• PhotoAbsCS.h
• PhotoAbsCS.cpp