G4WendtFissionFragmentGenerator Class Reference

#include <G4WendtFissionFragmentGenerator.hh>

Public Member Functions
G4HadFinalState *	ApplyYourself (const G4HadProjectile &projectile, G4int Z, G4int A)
void	InitializeANucleus (const G4int A, const G4int Z, const G4int M, const G4String &dataDirectory)
	~G4WendtFissionFragmentGenerator ()

Static Public Member Functions
static G4WendtFissionFragmentGenerator *	GetInstance ()

Detailed Description

Definition at line 41 of file G4WendtFissionFragmentGenerator.hh.

Constructor & Destructor Documentation

~G4WendtFissionFragmentGenerator()

G4WendtFissionFragmentGenerator::~G4WendtFissionFragmentGenerator

(

)

Definition at line 187 of file G4WendtFissionFragmentGenerator.cc.

{
  std::map<const G4int, G4FissionFragmentGenerator*>::iterator fissionGenerator;
 
  for (fissionGenerator = fissionIsotopes.begin(); fissionGenerator != fissionIsotopes.end();
       ++fissionGenerator)
  {
    delete fissionGenerator->second;
  }
}

Member Function Documentation

ApplyYourself()

G4HadFinalState * G4WendtFissionFragmentGenerator::ApplyYourself	(	const G4HadProjectile &	projectile,
		G4int	Z,
		G4int	A )

Definition at line 62 of file G4WendtFissionFragmentGenerator.cc.

{
  G4FFG_FUNCTIONENTER__
 
  G4HadFinalState* finalState = nullptr;
  G4DynamicParticleVector* finalParticles = nullptr;
  G4int isotope;
  std::map<const G4int, G4FissionFragmentGenerator*>::iterator fissionGenerator;
 
  // Look for the first available isomer since no M is provided for ApplyYourself()
  for (unsigned int M = 0; M < 10; ++M) {
    isotope = G4FissionFragmentGenerator::G4MakeIsotopeCode(Z, A, M);
    fissionGenerator = fissionIsotopes.find(isotope);
 
    if (fissionGenerator != fissionIsotopes.end()) {
      // Only generate particles if the generator was constructed
      if (fissionGenerator->second != nullptr) {
        finalParticles = fissionGenerator->second->G4GenerateFission(projectile);
      }
 
      break;
    }
  }
 
  if (finalParticles != nullptr) {
    finalState = new G4HadFinalState();
 
    for (auto& finalParticle : *finalParticles) {
      finalState->AddSecondary(finalParticle, secID);
    }
  }
 
  // TK modified 131108 add next line
  // TK 160112 fix for coverity #53481
  if (finalState != nullptr) finalState->SetStatusChange(stopAndKill);
  G4FFG_FUNCTIONLEAVE__
  return finalState;
}

Referenced by G4ParticleHPChannel::ApplyYourself().

GetInstance()

static G4WendtFissionFragmentGenerator * G4WendtFissionFragmentGenerator::GetInstance ( )

inlinestatic

Definition at line 45 of file G4WendtFissionFragmentGenerator.hh.

    {
      if (instance == nullptr) instance = new G4WendtFissionFragmentGenerator();
      return instance;
    }

Referenced by G4ParticleHPChannel::G4ParticleHPChannel().

InitializeANucleus()

void G4WendtFissionFragmentGenerator::InitializeANucleus	(	const G4int	A,
		const G4int	Z,
		const G4int	M,
		const G4String &	dataDirectory )

Definition at line 102 of file G4WendtFissionFragmentGenerator.cc.

{
  // G4FFG_FUNCTIONENTER__
 
  const G4int isotope = G4FissionFragmentGenerator::G4MakeIsotopeCode(Z, A, M);
  G4FFGEnumerations::MetaState metaState;
  std::pair<std::map<const G4int, G4FissionFragmentGenerator*>::iterator, bool> newIsotope;
 
  // Check to see if the isotope/isomer alread exists in the table
  newIsotope =
    fissionIsotopes.insert(std::make_pair(isotope, (G4FissionFragmentGenerator*)nullptr));
 
  if (newIsotope.second || newIsotope.first->second == NULL) {
    // Get the data file
    G4bool flag;
    G4ParticleHPDataUsed dataFile = fileNames.GetName(A, Z, M, dataDirectory, "FF", flag);
    G4String dataFileName = dataFile.GetName();
 
    // Check if the file exists, and do not create a fission object if it doesn't
    // G4cout << "*** Z = " << Z << "\tA = " << A << "\t\t\t Directory: "<< dataDirectory << " DATA
    // FILE: " << dataFileName << G4endl;
    std::istringstream dataStream(std::ios::in);
    G4ParticleHPManager::GetInstance()->GetDataStream(dataFileName, dataStream);
    if (!dataStream) {
      // G4FFG_FUNCTIONLEAVE__
      //  G4cerr << "*** Stream error" << G4endl;
      return;
    }
 
    // Check the data file parameters
    if (!flag
        || (Z < 2.5
            && ((G4double)abs(dataFile.GetZ() - Z) > 0.001
                || (G4double)abs((G4int)dataFile.GetA() - A) > 0.0001)))
    {
      // G4cerr << "*** Something wrong with the data request.\tFlag :" << flag << G4endl;
      // G4FFG_FUNCTIONLEAVE__
      return;
    }
 
    auto const fissionGenerator = new G4FissionFragmentGenerator();
    newIsotope.first->second = fissionGenerator;
 
    switch (M) {
      case 1:
        metaState = G4FFGEnumerations::META_1;
        break;
 
      case 2:
        metaState = G4FFGEnumerations::META_2;
        break;
 
      default:
        // TODO Display a warning message here indicating that an invalid metastate was passed in
        // Fall through to the ground state by default
      case 0:
        metaState = G4FFGEnumerations::GROUND_STATE;
        break;
    }
 
    fissionGenerator->G4SetIsotope(isotope);
    fissionGenerator->G4SetMetaState(metaState);
    fissionGenerator->G4SetCause(G4FFGEnumerations::NEUTRON_INDUCED);
    // TODO Load all the fission data and use the projectile energy instead
    fissionGenerator->G4SetIncidentEnergy(G4FFGDefaultValues::ThermalNeutronEnergy);
    fissionGenerator->G4SetYieldType(G4FFGEnumerations::INDEPENDENT);
    fissionGenerator->G4SetSamplingScheme(G4FFGEnumerations::NORMAL);
 
    // TODO Remove the need for forcing a load in the initialization phase,
    //      i.e. remove the ability to dynamically change the fission parameters
    //      that cause reload because a G4FissionFragmentGenerator class for
    //      each isotope should be loaded in the initialization phase
    if (!fissionGenerator->InitializeFissionProductYieldClass(dataStream)) {
      // Delete if the initialization fails
      delete fissionGenerator;
 
      fissionIsotopes.erase(newIsotope.first);
    }
  }
 
  // G4FFG_FUNCTIONLEAVE__
}

Referenced by G4ParticleHPChannel::UpdateData().

---

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

• G4WendtFissionFragmentGenerator.hh
• G4WendtFissionFragmentGenerator.cc