---

MyLiSource Class Reference

MyLiSource generates delayed neutrons from radioactive Li decays. More...

#include <MyLiSource.hh>

List of all members.

Public Methods
	MyLiSource (int newIsotope=k.Isotope, int newRmaxgen=k.Rmaxgen)
	MyLiSource constructor. More...
	~MyLiSource ()
	MyLiSource destructor.
	MyLiSource (const MyLiSource &LiSource)
	MyLiSource copy constructor.
MyLiSource &	operator= (const MyLiSource &rhs)
	MyLiSource overloaded = operator.
double	GetLiNeutronEnergy () const
	Returns the energy of the neutrons produced by each Li decay.
double	GetLiBetaEnergy () const
	Returns the energy of the betas produced by each Li decay.
double	GetLiGammaEnergy () const
	Returns the energy of the photons produced by each Li decay.
double	GetLiWeight () const
	Returns the branching ratio weight for each Li decay.
double	GetLiDecayTime () const
	Returns the decay time for each Li decay.
Static Public Methods
float	BetaHighestQ (const float &x)
float	BetaHighQ (const float &x)
float	BetaMiddleQ (const float &x)
float	BetaLowQ (const float &x)
double	Fermi (const int &Z, const int &A, const float &x)

---

Detailed Description

MyLiSource generates delayed neutrons from radioactive Li decays.

The supported isotope is 9. The neutrons are generated by the Li9NeutronSpectrum function, which the probability density of the produced neutrons as a function of neutron energy.

Definition at line 18 of file MyLiSource.hh.

---

Constructor & Destructor Documentation

MyLiSource::MyLiSource (  int    newIsotope = k.Isotope, int    newRmaxgen = k.Rmaxgen )

MyLiSource constructor. Called with the desired Li isotope: 9 and the maximum allowed radius for generating particles in the simulation, Rmaxgen (default = 200 cm). Definition at line 18 of file MyLiSource.cpp. { Isotope = newIsotope; Rmaxgen = newRmaxgen; //cout << "MyLiSource::MyLiSource(" << Isotope << "," << Rmaxgen // << ")" << endl; // 9Li if(Isotope == 9){ // setup the decay parameters (from Jon Link) // Nneutron is always 1 and Ngamma always 0 Nneutron = 1; Ngamma = 0; // branching fraction (the 1.0000 bin contains all decays // which do nothing: neutron E = 0, beta Q = 0) float static br[8] = {0.0238,0.3400,0.4400,0.4475, 0.4550,0.4894,0.4950,1.0000}; // Q-value for beta decay float static betaQ[8] = {11.176,11.176,10.826,5.666, 5.666,2.323,2.323,0.000}; // Q-value for the total decay chain float static totQ = 12.033; // neutron energy float static neuE[8] = {0.764,0.856,1.115,3.235, 5.473,8.816,6.578,0.000}; // prepare the beta decay functions bool static first=1; float static fspace[7][200]; if(first){ float xlo = 0.; funlxp_(BetaHighestQ,fspace[0],xlo,betaQ[0]); funlxp_(BetaHighestQ,fspace[1],xlo,betaQ[1]); funlxp_(BetaHighQ,fspace[2],xlo,betaQ[2]); funlxp_(BetaMiddleQ,fspace[3],xlo,betaQ[3]); funlxp_(BetaMiddleQ,fspace[4],xlo,betaQ[4]); funlxp_(BetaLowQ,fspace[5],xlo,betaQ[5]); funlxp_(BetaLowQ,fspace[6],xlo,betaQ[6]); first = 0; } // // let the 9Li decay // int ntrial = 0; bool passed = false; while(!passed){ int len=1; float rv[len]; ranlux_(rv,len); // for(int i=0;i<8; i++){ if(rv[0] < br[i] && !passed){ // pick an energy for each neutron neutronE = neuE[i]; if(neutronE){ // get the energy of the beta len = 1; funlux_(fspace[i],betaE,len); passed = true; } } } ntrial++; } Weight = 1./float(ntrial); //cout << " 9Li neutron energy = " << neutronE << " MeV" << endl; //cout << " 9Li beta energy = " << betaE << " MeV" << endl; //cout << " 9Li weight = " << Weight << endl; // pick a decay time in ns double HalfLife = 1.8e8; int len=1; float rv[len]; ranlux_(rv,len); DecayTime = HalfLife*log(1/rv[0]); } // done with 9Li }

---

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

• MyLiSource.hh
• MyLiSource.cpp

---