Heed::AtomicSecondaryProducts Class Reference

#include <PhotoAbsCS.h>

Inheritance diagram for Heed::AtomicSecondaryProducts:

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

Protected Attributes
DynLinArr< double >	channel_prob_dens
DynLinArr< DynLinArr< double > >	electron_energy
DynLinArr< DynLinArr< double > >	photon_energy

Detailed Description

Definition at line 290 of file PhotoAbsCS.h.

Constructor & Destructor Documentation

AtomicSecondaryProducts()

Heed::AtomicSecondaryProducts::AtomicSecondaryProducts ( )

inline

Definition at line 296 of file PhotoAbsCS.h.

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

~AtomicSecondaryProducts()

virtual Heed::AtomicSecondaryProducts::~AtomicSecondaryProducts ( )

inlinevirtual

Definition at line 300 of file PhotoAbsCS.h.

{}

Member Function Documentation

add_channel()

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

Definition at line 875 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.get_qel();
  long q_new = q_old + 1;
  channel_prob_dens.put_qel(q_new);
  electron_energy.put_qel(q_new);
  photon_energy.put_qel(q_new);
  //Iprintn(mcout, s);
  //check_econd11( s , > 1.0 , mcerr);
  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);
  }
}

Referenced by Heed::generate_Argon_PACS_mod_esc().

get_channel()

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

Definition at line 914 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.get_qel());
#endif
  int ir = 0;
  if (channel_prob_dens.get_qel() > 0) {
    double rn = SRANLUX();
#ifdef DEBUG_PRINT_get_escape_particles
    Iprintn(mcout, rn);
#endif
    if (channel_prob_dens.get_qel() == 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.get_qel();
      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

Definition at line 958 of file PhotoAbsCS.cpp.

                                                                 {
  if (l > 0) {
    Ifile << "AtomicSecondaryProducts(l=" << l << "):\n";
    long q = channel_prob_dens.get_qel();
    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].get_qel();
      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].get_qel();
      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

DynLinArr<double> Heed::AtomicSecondaryProducts::channel_prob_dens

protected

Definition at line 316 of file PhotoAbsCS.h.

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

electron_energy

DynLinArr<DynLinArr<double> > Heed::AtomicSecondaryProducts::electron_energy

protected

Definition at line 318 of file PhotoAbsCS.h.

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

photon_energy

DynLinArr<DynLinArr<double> > Heed::AtomicSecondaryProducts::photon_energy

protected

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