MJGeneratorPNNLDecayChain.hh

Go to the documentation of this file.

// This code implementation is the intellectual property of
// the RD44 GEANT4 collaboration.
//
// By copying, distributing or modifying the Program (or any work
// based on the Program) you indicate your acceptance of this statement,
// and all its terms.
//
// $Id: MJGeneratorPNNLDecayChain.hh,v 1.2 2004/11/09 13:42:39 xliu Exp $
// GEANT4 tag $Name:  $
//
// File:        MJGeneratorPNNLDecayChain.hh
// Description: Determines relative freqency of radioisotope decays and
//              stores radioisotope objects in PNNL radiactive decay model
// Version:     1.0
// Created:     2004-05-08
// Author:      David Jordan
// mail:        david.jordan@pnl.gov
// phone:       509-376-0115

#ifndef MJGeneratorPNNLDecayChain_h
#define MJGeneratorPNNLDecayChain_h 1

#include "generators/MJGeneratorPNNLRadioisotope.hh"
#include "generators/MJGeneratorPNNLLookup.hh"
#include "globals.hh"

class MJGeneratorPNNLDecayChain 
{
  public:
    MJGeneratorPNNLDecayChain();
    MJGeneratorPNNLDecayChain(G4String decay_chain_filename, 
                              G4double source_age);
    MJGeneratorPNNLDecayChain(const MJGeneratorPNNLDecayChain& DecayChain); 
    ~MJGeneratorPNNLDecayChain();
    void operator=(const MJGeneratorPNNLDecayChain& DecayChain); 

    void SetDecayChainFilename(G4String decay_chain_filename) {
                        fDecayChain_filename = decay_chain_filename; };
    void SetSourceAge(G4double source_age) {fSource_age = source_age;};

    G4String GetDecayChainFilename() const {return fDecayChain_filename; };
    G4double GetSourceAge()          const {return fSource_age/kYear; };
    G4int    GetIndexOfLastIsotopeSampled() const 
                                     {return fLastIsotopeSampled;};

    MJGeneratorPNNLCascadeEvent DoCascadeEvent();
    void DoInitializeGenerator();
    void DoPrintRadioisotopeList() const;

  private:
    G4String   fDecayChain_filename;
    G4int      fNum_RadioIsotopes;  
    G4double*  fActivity_list;
    G4double*  f_pT_half_list;
    char*      f_pT_half_units_list;
    G4double*  f_pLamda_list;
    G4String*  fRadioIsotope_label_list;
    G4double   fSource_age;
    G4int      fLastIsotopeSampled;
  
    MJGeneratorPNNLRadioisotope** fRadioIsotope_list;
    MJGeneratorPNNLLookup* f_pRadioIsotopeSelector;

    void CalculateActivityList();

    const static G4int kSecond;
    const static G4int kMinute;
    const static G4int kHour;
    const static G4int kDay;
    const static G4int kYear;
};

#endif

---