Decay0 Class Reference

Ratpac-two: RAT::Decay0 Class Reference

Ratpac-two

The Decay0 Generator for Initial Kinematics in Alpha, Beta and Double Beta Decays. More...

#include <Decay0.hh>

Public Member Functions
	Decay0 ()
	‍************************************************‍/

	Decay0 (const std::string isotope, const int level, const int mode, const double lE, const double hE)
	‍************************************************‍/

	Decay0 (const std::string isotope)
	‍************************************************‍/

	~Decay0 ()
	‍************************************************‍/

void	GenBBTest ()
	‍************************************************‍/

void	GenBackgTest ()
	‍************************************************‍/

void	GenBBDeex ()
	‍************************************************‍/ More...

void	GenEvent ()
	‍************************************************‍/ More...

double	GetRandom ()
	‍************************************************‍/

double	GetMass (int pdg)
	‍************************************************‍/

void	bb ()
	‍************************************************‍/

void	particle (int np, double E1, double E2, double teta1, double teta2, double phi1, double phi2, double tclev, double thlev)

void	pair (double Epair)
	‍************************************************‍/ More...

double	fe1_mod ()
	probability distribution for energy of e-/e+ in doublebeta decay More...

double	fe2_mod ()
	‍************************************************‍/

double	fermi (const double &Z, const double &E)
	‍************************************************‍/ More...

void	beta (double Qbeta, double tcnuc, double thnuc)
	‍************************************************‍/ More...

void	beta1f (double Qbeta, double tcnuc, double thnuc, double c1, double c2, double c3, double c4)
	‍************************************************‍/

void	beta1fu (double Qbeta, double tcnuc, double thnuc, double c1, double c2, double c3, double c4)
	‍************************************************‍/

void	beta2f (double Qbeta, double tcnuc, double thnuc, int kf, double c1, double c2, double c3, double c4)
	‍************************************************‍/

void	tgold (double a, double b, TF1 &f, double eps, int nmin, double &xextt, double &fextr)
	‍************************************************‍/ More...

void	Ti48low ()
	‍************************************************‍/ More...

void	Fe58low ()

void	Se76low ()

void	Ge74low ()

void	Kr82low ()

void	Mo96low ()

void	Zr92low ()

void	Ru100low ()

void	Pd106low ()

void	Sn116low ()

void	Cd112low ()

void	Te124low ()

void	Xe130low ()

void	Ba136low ()

void	Sm148low ()

void	Sm150low ()

void	As79 ()
	functions for decays (alpha/beta) More...

void	At214 ()

void	Ac228 ()

void	Bi207 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Bi210 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Bi212 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Bi214 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Co60 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Cs136 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Eu147 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Eu152 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Eu154 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Gd146 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	I126 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	I133 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	I134 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	I135 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	K40 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	K42 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Pa234m ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Pb211 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Pb212 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Pb214 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Po214 ()

void	Rn218 ()

void	Ra222 ()

void	Ra228 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Rh106 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Sb125 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Sb126 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Sb133 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Sc48 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Ta182 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Te133 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Te133m ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Te134 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Th234 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Tl208 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Xe133 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Xe135 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Y88 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Zn65 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Nb96 ()
	^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

void	Be11 ()

void	PbAtShell (int KLMenergy)
	‍************************************************‍/ More...

void	nucltransK (double Egamma, double Ebinde, double conve, double convp)
	‍************************************************‍/ More...

void	nucltransKL (double Egamma, double EbindeK, double conveK, double EbindeL, double conveL, double convp)
	‍************************************************‍/

void	nucltransKLM (double Egamma, double EbindeK, double conveK, double EbindeL, double conveL, double EbindeM, double conveM, double convp)
	‍************************************************‍/

void	nucltransKLM_Pb (double Egamma, double EbindeK, double conveK, double EbindeL, double conveL, double EbindeM, double conveM, double convp)
	‍************************************************‍/

Double_t	funbeta (Double_t x, Double_t par)
	‍************************************************‍/

Double_t	funbeta1fu (Double_t x, Double_t par)
	‍************************************************‍/

Double_t	funbeta1f (Double_t x, Double_t par)
	‍************************************************‍/

Double_t	funbeta2f (Double_t x, Double_t par)
	‍************************************************‍/

std::complex< double >	cgamma (std::complex< double > z)
	‍************************************************‍/

double	divdif (double xtab[50], double xval)
	‍************************************************‍/

int	GetNbPart ()
	‍************************************************‍/

double	GetPmoment (int i, int j)

double	GetPtime (int i)

int	GetNpGeant (int i)

void	SetCutoffWindow (double time)

double	GetCutoffWindow ()

void	SetTimeCutoff (bool status)

bool	GetTimeCutoff ()

void	SetAlphaCut (bool status)

bool	GetAlphaCut ()

unsigned int	GetParentIdx (unsigned int i)

double	operator() (double x, double par)

Public Attributes
G4ParticleDefinition *	fPartDef

double	fSl [48]

int	fSlSize

double	fC1

double	fC2

double	fC3

double	fC4

double	fKf

int	fNbPart

double	fPmoment [3][100]

double	fPtime [100]

int	fNpGeant [100]

unsigned int	fCurParentIdx

std::vector< unsigned int >	fPparent

Detailed Description

The Decay0 Generator for Initial Kinematics in Alpha, Beta and Double Beta Decays.

Author: V. Tretyak and the SNO+ Collaboration

History:

The Decay0 generator was originally written by V. Tretyak in FORTRAN; see Reference [1]. This code was ported to C++ by members of the SNO+ Collaboration [2] for internal use. This C++ implementation has been made available in RAT-PAC by agreement of the SNO+ Collaboration.

See the COPYING and COPYING.decay0 files, which should be included with this source code distribution, for license details.

[1] O.A.Ponkratenko, V.I.Tretyak, Yu.G.Zdesenko, Phys. At. Nucl. 63 (2000) 1282, (nucl-ex/0104018). [2] http://snoplus.phy.queensu.ca/

Units:

Energy MeV Momentum MeV/c Time seconds

List of the isotopes included in the generator:

Ac228 Ar39 Ar42 As79 Be11 Bi207 Bi208 Bi210 Bi212 Bi214 C10 C11 C14 Ca48 Cd113 Co60 Cs136 Cs137 Eu147 Eu152 Eu154 Gd146 He6 Hf182 I126 I133 I134 I135 K40 K42 Kr81 Kr85 Mn54 Na22 P32 Pa234m Pb210 Pb211 Pb212 Pb214 Ra228 Rb87 Rh106 Sb125 Sb126 Sb133 Sr90 Ta182 Te133 Te133m Te134 Th234 Tl207 Tl208 Xe129m Xe131m Xe133 Xe135 Y88 Y90 Zn65 Zr96

Usage:

In the macro, define:

/generator/add combo decay0:[POSITION]

where POSITION is a generator generator (point, fill, ...).

Examples:

To simulate the double-beta decay events:

/generator/vtx/set 2beta [ISOTOPE] [LEVEL] [MODE] [LoELim] [HiELim]

To simulate background events:

/generator/vtx/set backg [ISOTOPE]

Revision history: 2017/08/29 - A. Mastbaum: Import from SNO+ RAT

Member Function Documentation

◆ Ac228()

void RAT::Decay0::Ac228 ( )

----—end of functions

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ ^^^^^^^ Section of radioactive decays ^^ ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

◆ As79()

void RAT::Decay0::As79 ( )

functions for decays (alpha/beta)

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

◆ At214()

void RAT::Decay0::At214 ( )

◆ Ba136low()

void RAT::Decay0::Ba136low ( )

◆ beta()

void RAT::Decay0::beta	(	double	Qbeta,
		double	tcnuc,
		double	thnuc
	)

‍************************************************‍/

simulation of the angles and energy of beta particles emitted in beta decay of nucleus.

◆ Cd112low()

void RAT::Decay0::Cd112low ( )

◆ fe1_mod()

double RAT::Decay0::fe1_mod ( )

probability distribution for energy of e-/e+ in doublebeta decay

‍************************************************‍/

◆ Fe58low()

void RAT::Decay0::Fe58low ( )

◆ fermi()

double RAT::Decay0::fermi	(	const double &	Z,
		const double &	E
	)

‍************************************************‍/

Function to calculate the traditional function of Fermi in theory of beta decay to take into account the Coulomb correction

◆ Ge74low()

void RAT::Decay0::Ge74low ( )

◆ GenBBDeex()

void RAT::Decay0::GenBBDeex ( )

‍************************************************‍/

Simulation of doublebeta-like radioactive decays

◆ GenEvent()

void RAT::Decay0::GenEvent ( )

‍************************************************‍/

decay of Po212

decay of Po214

◆ Kr82low()

void RAT::Decay0::Kr82low ( )

◆ Mo96low()

void RAT::Decay0::Mo96low ( )

◆ nucltransK()

void RAT::Decay0::nucltransK	(	double	Egamma,
		double	Ebinde,
		double	conve,
		double	convp
	)

‍************************************************‍/

◆ pair()

void RAT::Decay0::pair ( double Epair )

‍************************************************‍/

Generation of e+e- pair in zero-approximation to real subroutine for internal pair creation: 1) energy of e+ is equal to the energy of e-; 2) e+ and e- are emitted in the same direction

◆ particle()

void RAT::Decay0::particle	(	int	np,
		double	E1,
		double	E2,
		double	teta1,
		double	teta2,
		double	phi1,
		double	phi2,
		double	tclev,
		double	thlev
	)

Generation of isotropic emission of particle in the range of energies and angles

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ End of functions ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

◆ PbAtShell()

void RAT::Decay0::PbAtShell ( int KLMenergy )

‍************************************************‍/

describes in some approximation the deexcitation process in atomic shell of Pb after creation of electron vacation in K, L or M shell.

◆ Pd106low()

void RAT::Decay0::Pd106low ( )

◆ Po214()

void RAT::Decay0::Po214 ( )

◆ Ra222()

void RAT::Decay0::Ra222 ( )

◆ Rn218()

void RAT::Decay0::Rn218 ( )

◆ Ru100low()

void RAT::Decay0::Ru100low ( )

◆ Se76low()

void RAT::Decay0::Se76low ( )

◆ Sm148low()

void RAT::Decay0::Sm148low ( )

◆ Sm150low()

void RAT::Decay0::Sm150low ( )

◆ Sn116low()

void RAT::Decay0::Sn116low ( )

◆ Te124low()

void RAT::Decay0::Te124low ( )

◆ tgold()

void RAT::Decay0::tgold	(	double	a,
		double	b,
		TF1 &	f,
		double	eps,
		int	nmin,
		double &	xextt,
		double &	fextr
	)

‍************************************************‍/

determines maximum or minimum of the function f(x) in the interval [a,b] by the gold section method.

◆ Ti48low()

void RAT::Decay0::Ti48low ( )

‍************************************************‍/

functions describe the de-excitation process after 2b-decay to ground and excited 0+ and 2+ levels conditions still hard-coded - work in progress

-—functions describe the deexcitation process after 2b-decay

◆ Xe130low()

void RAT::Decay0::Xe130low ( )

◆ Zr92low()

void RAT::Decay0::Zr92low ( )

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

/home/docs/checkouts/readthedocs.org/user_builds/ratpac/checkouts/latest/src/gen/include/RAT/Decay0.hh
/home/docs/checkouts/readthedocs.org/user_builds/ratpac/checkouts/latest/src/gen/src/Decay0.cc

Decay0 Class Reference

Public Member Functions

Public Attributes

Detailed Description

Member Function Documentation

◆ Ac228()

----—end of functions

◆ As79()

◆ At214()

◆ Ba136low()

◆ beta()

◆ Cd112low()

◆ fe1_mod()

◆ Fe58low()

◆ fermi()

◆ Ge74low()

◆ GenBBDeex()

◆ GenEvent()

◆ Kr82low()

◆ Mo96low()

◆ nucltransK()

◆ pair()

◆ particle()

◆ PbAtShell()

◆ Pd106low()

◆ Po214()

◆ Ra222()

◆ Rn218()

◆ Ru100low()

◆ Se76low()

◆ Sm148low()

◆ Sm150low()

◆ Sn116low()

◆ Te124low()

◆ tgold()

◆ Ti48low()

-—functions describe the deexcitation process after 2b-decay

◆ Xe130low()

◆ Zr92low()