Heed::AtomicSecondaryProducts Class Reference

#include <PhotoAbsCS.h>

Inheritance diagram for Heed::AtomicSecondaryProducts:

Public Member Functions
	AtomicSecondaryProducts ()
	Constructor.
virtual	~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)
virtual void	print (std::ostream &file, int l) const
Public Member Functions inherited from Heed::RegPassivePtr
	RegPassivePtr (void)
	RegPassivePtr (char fs_ban_del, char fs_ban_sub, char fs_ban_cop=0)
	RegPassivePtr (const RegPassivePtr &f)
RegPassivePtr &	operator= (const RegPassivePtr &f)
CountPP_ns::CountPassivePtr *	book (void) const
void	clear_pointers (void) const
virtual RegPassivePtr *	copy () const
virtual	~RegPassivePtr ()
virtual void	print (std::ostream &file, int l=1) const
void	set_s_ban_del (char fs_ban_del)
char	get_s_ban_del (void) const
void	set_s_ban_sub (char fs_ban_sub)
char	get_s_ban_sub (void) const
void	set_s_ban_cop (char fs_ban_cop)
char	get_s_ban_cop (void) const
void	set_s_allow_del_at_zero_count (char fs_allow_del_at_zero_count)
char	get_s_allow_del_at_zero_count (void) const
long	get_total_number_of_references (void) const

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

Additional Inherited Members
Static Public Member Functions inherited from Heed::RegPassivePtr
static void	set_s_ban_del_ignore (char fs_ban_del_ignore)
static char	get_s_ban_del_ignore (void)
static void	set_s_print_adr_cpp (char fs_print_adr_cpp)
static char	get_s_print_adr_cpp (void)

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 249 of file PhotoAbsCS.h.

Constructor & Destructor Documentation

AtomicSecondaryProducts()

Heed::AtomicSecondaryProducts::AtomicSecondaryProducts ( )

inline

Constructor.

Definition at line 252 of file PhotoAbsCS.h.

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

~AtomicSecondaryProducts()

virtual Heed::AtomicSecondaryProducts::~AtomicSecondaryProducts ( )

inlinevirtual

Destructor.

Definition at line 255 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

virtual

Reimplemented from Heed::RegPassivePtr.

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 278 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 280 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 281 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