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 39 of file G4WendtFissionFragmentGenerator.hh.

Constructor & Destructor Documentation

G4WendtFissionFragmentGenerator::~G4WendtFissionFragmentGenerator

(

)

Definition at line 183 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

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

Definition at line 54 of file G4WendtFissionFragmentGenerator.cc.

References G4FFG_FUNCTIONENTER__, G4FFG_FUNCTIONLEAVE__, G4FissionFragmentGenerator::G4MakeIsotopeCode(), and stopAndKill.

Referenced by G4NeutronHPChannel::ApplyYourself().

{
G4FFG_FUNCTIONENTER__

    G4HadFinalState* finalState = NULL;
    G4DynamicParticleVector* finalParticles = NULL;
    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)
            {
                finalParticles = fissionGenerator->second->G4GenerateFission(projectile);
            }

            break;
        }
    }

    if(finalParticles)
    {
        finalState = new G4HadFinalState();

        for(unsigned int i = 0; i < finalParticles->size(); ++i)
        {
            finalState->AddSecondary((*finalParticles)[i]);
        }
    }

    //TK modified 131108 add next line 
    finalState->SetStatusChange(stopAndKill);
G4FFG_FUNCTIONLEAVE__
    return finalState;
}

G4WendtFissionFragmentGenerator * G4WendtFissionFragmentGenerator::GetInstance ( void  )

static

Definition at line 46 of file G4WendtFissionFragmentGenerator.cc.

{
    static G4WendtFissionFragmentGenerator newMe;

    return &newMe;
}

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

Definition at line 98 of file G4WendtFissionFragmentGenerator.cc.

References G4FissionFragmentGenerator::G4MakeIsotopeCode(), G4FissionFragmentGenerator::G4SetCause(), G4FissionFragmentGenerator::G4SetIncidentEnergy(), G4FissionFragmentGenerator::G4SetIsotope(), G4FissionFragmentGenerator::G4SetMetaState(), G4FissionFragmentGenerator::G4SetSamplingScheme(), G4FissionFragmentGenerator::G4SetYieldType(), G4NeutronHPDataUsed::GetA(), G4NeutronHPManager::GetDataStream(), G4NeutronHPManager::GetInstance(), G4NeutronHPNames::GetName(), G4NeutronHPDataUsed::GetName(), G4NeutronHPDataUsed::GetZ(), G4FFGEnumerations::GROUND_STATE, G4FFGEnumerations::INDEPENDENT, G4FissionFragmentGenerator::InitializeFissionProductYieldClass(), G4FFGEnumerations::META_1, G4FFGEnumerations::META_2, G4FFGEnumerations::NEUTRON_INDUCED, and G4FFGEnumerations::NORMAL.

Referenced by G4NeutronHPChannel::UpdateData().

{
//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*)NULL));

    if(newIsotope.second || newIsotope.first->second == NULL)
    {
        // Get the data file
        G4bool flag;
        G4NeutronHPDataUsed 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);
        G4NeutronHPManager::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;
        }

        G4FissionFragmentGenerator* 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__
}

---

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

• G4WendtFissionFragmentGenerator.hh
• G4WendtFissionFragmentGenerator.cc