nucltransKLM_Pb.h

Go to the documentation of this file.

#ifndef GENBB_DECAY0_NUCLTRANSKLM_PB_H_
#define GENBB_DECAY0_NUCLTRANSKLM_PB_H_ 1

// This project:
#include <genbb_help/decay0/nucltransKLM_Pb.h>
#include <genbb_help/decay0/particle.h>
#include <genbb_help/decay0/gamma.h>
#include <genbb_help/decay0/pair.h>
#include <genbb_help/decay0/electron.h>

namespace mygsl{
  class rng;
}

namespace genbb {

  class primary_event;

  namespace decay0 {

    // The same as nucltransKLM but two X rays are emitted after K conversion
    // in deexcitation of 208-Pb in decay 208Tl->208Pb.
    // VIT, 4.02.2009.
    //
    // Subroutine nucltransKLM choises 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 three fixed energies:
    // Egamma-E(K)_binding_energy, Egamma-E(L)_binding_energy and
    // Egamma-E(M)_binding_energy).
    // Call  : common/genevent/tevst,npfull,npgeant(100),pmoment(3,100),ptime(100)
    //         call nucltransKLM(Egamma,EbindeK,conveK,EbindeL,conveL,
    //                           EbindeM,conveM,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;
    //         EbindeM - binding energy of electron (MeV) on M-shell;
    //         convem  - internal conversion coefficient [=Nelectron/Ngamma] from
    //                   M-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, 4.01.2007.
    void decay0_nucltransKLM_Pb(mygsl::rng & prng, primary_event & event,
                                double Egamma,double EbindeK,double conveK,
        double EbindeL,double conveL,
                                double EbindeM,double conveM,
        double convp,
        double tclev,double thlev,double & tdlev);

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

#endif // GENBB_DECAY0_NUCLTRANSKLM_PB_H_