nucltransKL.h

Go to the documentation of this file.

#ifndef GENBB_DECAY0_NUCLTRANSKL_H_
#define GENBB_DECAY0_NUCLTRANSKL_H_ 1

// This project:
#include <genbb_help/decay0/common.h>

namespace mygsl{
  class rng;
}

namespace genbb {

  class primary_event;

  namespace decay0 {

    // Subroutine nucltransKL choise one of the three concurent processes
    // by which the transition from one nuclear state to another is
    // occured: gamma-ray emission, internal conversion and internal
    // pair creation. Conversion electrons are emitted with two fixed energies
    // (Egamma-E(K)_binding_energy and Egamma-E(L)_binding_energy).
    // Call  : common/genevent/tevst,npfull,npgeant(100),pmoment(3,100),ptime(100)
    //         call nucltransKL(Egamma,EbindeK,conveK,EbindeL,conveL,convp,
    //                          tclev,thlev,tdlev)
    // Input:  Egamma  - gamma-ray energy (MeV) [=difference in state energies];
    //         EbindeK - binding energy of electron (MeV) on K-shell;
    //         conveK  - internal conversion coefficient [=Nelectron/Ngamma] from
    //                   K-shell;
    //         EbindeL - binding energy of electron (MeV) on L-shell;
    //         conveL  - internal conversion coefficient [=Nelectron/Ngamma] from
    //                   L-shell;
    //         convp   - pair conversion coefficient [=Npair/Ngamma];
    //         tclev   - time of creation of level from which particle will be
    //                   emitted (sec);
    //         thlev   - level halflife (sec).
    // Output: tdlev               - time of decay of level (sec);
    // VIT, 5.07.1995.
    void decay0_nucltransKL(mygsl::rng & prng, primary_event & event,
          double Egamma,double EbindeK,double conveK,
          double EbindeL,double conveL,double convp,
          double tclev,double thlev,
          double & tdlev);

  } // end of namespace decay0
} // end of namespace genbb

#endif // GENBB_DECAY0_NUCLTRANSKL_H_