---

MyTlSource.cpp

Go to the documentation of this file.

#include <math.h>
#include <iostream.h>
#include <fstream.h>   // file I/O
#include <iomanip.h>   // format manipulation
#include <string>
#include "ReactorDetector.hh"
#include "ReactorConstants.hh"
#include "ReactorFortran.hh"
#include "MyTlSource.hh"

MyTlSource::MyTlSource(int newIsotope){

  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

}

MyTlSource::~MyTlSource(){
  ;
}

MyTlSource::MyTlSource(const MyTlSource& TlSource){

  Isotope = TlSource.Isotope;

  numGamma = TlSource.numGamma;

  for(int n=0;n<6;n++){
    gammaEnergy[n] = TlSource.gammaEnergy[n];
  }
}    

MyTlSource& MyTlSource::operator=(const MyTlSource& rhs){

  if(this == &rhs)return *this;

  Isotope = rhs.Isotope;  

  numGamma = rhs.numGamma;

  for(int n=0;n<6;n++){
    gammaEnergy[n] = rhs.gammaEnergy[n];
  }

  return *this;

}    

---