---

MyTlSource Class Reference

MyTlSource generates photons from radioactive Thallium decays. More...

#include <MyTlSource.hh>

List of all members.

Public Methods
	MyTlSource (int newIsotope=k.Isotope)
	MyTlSource constructor. More...
	~MyTlSource ()
	MyTlSource destructor.
	MyTlSource (const MyTlSource &TlSource)
	MyTlSource copy constructor.
MyTlSource &	operator= (const MyTlSource &rhs)
	MyTlSource overloaded = operator.
int	GetNumTlGammas () const
	Returns the number of photons produced by the Tl decay.
float	GetTlGammaEnergy (int &n) const
	Returns the energy of each photon produced by the Tl decay. More...

---

Detailed Description

MyTlSource generates photons from radioactive Thallium decays.

The supported isotopes are Tl208 and Tl210. The decay table is read in from tl_bi_decays.txt and parsed by MyTlSource to create an array of the decay branch probabilites and the number and energies of the photons generated by each branch.

Definition at line 18 of file MyTlSource.hh.

---

Constructor & Destructor Documentation

MyTlSource::MyTlSource (  int    newIsotope = k.Isotope )

MyTlSource constructor. Called with the desired Thallium isotope: 208 or 210. Default is the isotope given in ReactorConstants, 252, which does nothing. Definition at line 17 of file MyTlSource.cpp. { Isotope = newIsotope; //cout << "MyTlSource::MyTlSource(" << Isotope << ")" << endl; // input the source number and position fstream fin("data/tl_bi_decays.txt", ios::in); //cout << "Tl input file: data/tl_bi_decays.txt" << endl; char paramLine[256]; string dummyString; int dummyInt; float dummyFloat; bool done = false; bool Set = false; int Z, A; int Nbranch; float branchProb[100]; int Ngamma[100]; float gammaE[100][6]; // Tl208 if(Isotope == 208){ while(!done){ // '#' specifies a comment line -- ignore it if(fin.peek() == '#'){ fin.getline(paramLine, 256); continue; } // 'E'ND or 'e'nd specifies end of file -- set done flag to true if(fin.peek() == 'e' || fin.peek() == 'E'){ done = true; continue; } // 'N'UCLEUS = or 'n'ucleus = specifies Z and A of daughter nucleus if(fin.peek() == 'n' || fin.peek() == 'N'){ if(Set){ fin.getline( paramLine, 256 ); continue; } fin >> dummyString; while(fin.peek() == ' ' || fin.peek() == '=') fin.ignore(1); fin >> Z >> A; //cout << Z << " " << A << endl; if(A == 208){ // read in Tl208 decay table, 40 branches Nbranch = 40; for(int n=0; n<Nbranch; n++){ fin >> branchProb[n] >> dummyInt >> dummyFloat >> Ngamma[n]; for(int i=0; i<6; i++) fin >> gammaE[n][i]; fin >> dummyFloat >> dummyFloat >> dummyFloat; /* cout << " branchProb[" << n << "] " << branchProb[n] << " and Ngamma[" << n << "] " << Ngamma[n] << endl; for(int i=0; i<Ngamma[n]; i++){ cout << " gammaE[" << n << "][" << i << "] " << gammaE[n][i]; } cout << endl; */ } done = true; continue; Set = true; } continue; } // ignore lines starting with numbers if they are not in our decay table if(fin.peek() == '0' || fin.peek() == '1'){ fin.getline(paramLine, 256); continue; } // ignore extra whitespace if(fin.peek() == ' ' || fin.peek() == '\n'){ fin.ignore(1); continue; } } // pick a decay branch and get the gamma energies int len = 1; float r[len]; ranlux_(r,len); bool passed = false; for(int n=0;n<Nbranch; n++){ if(branchProb[n] > r[0] && !passed){ // save the number and energies of the gammas numGamma = Ngamma[n]; for(int i=0;i<numGamma; i++) gammaEnergy[i] = gammaE[n][i]; passed = true; } } } // done with Tl208 // Tl210 if(Isotope == 210){ while(!done){ // '#' specifies a comment line -- ignore it if(fin.peek() == '#'){ fin.getline(paramLine, 256); continue; } // 'E'ND or 'e'nd specifies end of file -- set done flag to true if(fin.peek() == 'e' || fin.peek() == 'E'){ done = true; continue; } // 'N'UCLEUS = or 'n'ucleus = specifies Z and A of daughter nucleus if(fin.peek() == 'n' || fin.peek() == 'N'){ if(Set){ fin.getline( paramLine, 256 ); continue; } fin >> dummyString; while(fin.peek() == ' ' || fin.peek() == '=') fin.ignore(1); fin >> Z >> A; //cout << Z << " " << A << endl; if(A == 210){ // read in Tl210 decay table, 12 branches Nbranch = 12; for(int n=0; n<Nbranch; n++){ fin >> branchProb[n] >> dummyInt >> dummyFloat >> Ngamma[n]; for(int i=0; i<6; i++) fin >> gammaE[n][i]; fin >> dummyFloat >> dummyFloat >> dummyFloat; /* cout << " branchProb[" << n << "] " << branchProb[n] << " and Ngamma[" << n << "] " << Ngamma[n] << endl; for(int i=0; i<Ngamma[n]; i++){ cout << " gammaE[" << n << "][" << i << "] " << gammaE[n][i]; } cout << endl; */ } done = true; continue; Set = true; } continue; } // ignore lines starting with numbers if they are not in our decay table if(fin.peek() == '0' || fin.peek() == '1'){ fin.getline(paramLine, 256); continue; } // ignore extra whitespace if(fin.peek() == ' ' || fin.peek() == '\n'){ fin.ignore(1); continue; } } } // done with Tl210 }

---

Member Function Documentation

float MyTlSource::GetTlGammaEnergy (  int &    n )  const [inline]

Returns the energy of each photon produced by the Tl decay. Called with the integer index for each photon in the gammaEnergy array, from 0 to the total number of photons. Definition at line 38 of file MyTlSource.hh. {return gammaEnergy[n];}

---

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

• MyTlSource.hh
• MyTlSource.cpp

---