---

MyHeSource Class Reference

MyHeSource generates delayed neutrons from radioactive He decays. More...

#include <MyHeSource.hh>

List of all members.

Public Methods
	MyHeSource (int newIsotope=ISOTOPE, double newRmaxgen=RMAXGEN)
	MyHeSource constructor. More...
	~MyHeSource ()
	MyHeSource destructor.
	MyHeSource (const MyHeSource &HeSource)
	MyHeSource copy constructor.
MyHeSource &	operator= (const MyHeSource &rhs)
	MyHeSource overloaded = operator.
double	GetHeNeutronEnergy () const
	Returns the energy of the neutrons produced by each He decay.
double	GetHeBetaEnergy () const
	Returns the energy of the betas produced by each He decay.
double	GetHeGammaEnergy () const
	Returns the energy of the photons produced by each He decay.
double	GetHeWeight () const
	Returns the branching ratio weight for each He decay.
double	GetHeDecayTime () const
	Returns the decay time for each He decay.
Static Public Methods
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

MyHeSource generates delayed neutrons from radioactive He decays.

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

Definition at line 18 of file MyHeSource.hh.

---

Constructor & Destructor Documentation

MyHeSource::MyHeSource (  int    newIsotope = ISOTOPE, double    newRmaxgen = RMAXGEN )

MyHeSource constructor. Called with the desired He isotope: 8 and the maximum allowed radius for generating particles in the simulation, Rmaxgen (default = 200 cm). Definition at line 17 of file MyHeSource.cpp. { Isotope = newIsotope; Rmaxgen = newRmaxgen; //cout << "MyHeSource::MyHeSource(" << Isotope << "," << Rmaxgen // << ")" << endl; // 8He if(Isotope == 8){ // 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[4] = {0.080,0.160,0.169,1.0000}; // Q-value for beta decay float static betaQ[4] = {7.442,5.252,1.852,0.000}; // Q-value for the total decay chain float static totQ; totQ = 8.619; // neutron energy float static neuE[4] = {0.697,2.887,6.727,0.000}; // prepare the beta decay functions bool static first=1; float static fspace[3][200]; if(first){ float xlo = 0.; funlxp_(BetaHighQ,fspace[0],xlo,betaQ[0]); funlxp_(BetaMiddleQ,fspace[1],xlo,betaQ[1]); funlxp_(BetaLowQ,fspace[2],xlo,betaQ[2]); first = 0; } // // let the 8He decay // int ntrial = 0; bool passed = false; while(!passed){ const int len=1; float rv[len]; ranlux_(rv,len); // for(int i=0;i<4; i++){ if(rv[0] < br[i] && !passed){ // pick an energy for each neutron neutronE = neuE[i]; if(neutronE){ // get the energy of the beta int len2 = 1; funlux_(fspace[i],betaE,len2); passed = true; } } } ntrial++; } Weight = 1./float(ntrial); //cout << " 8He neutron energy = " << neutronE << " MeV" << endl; //cout << " 8He beta energy = " << betaE << " MeV" << endl; //cout << " 8He weight = " << Weight << endl; // pick a decay time in ns double HalfLife = 1.2e8; const int len=1; float rv[len]; ranlux_(rv,len); DecayTime = HalfLife*log(1/rv[0]); } // done with 8He }

---

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

• MyHeSource.hh
• MyHeSource.cpp

---