107 std::vector< G4DynamicParticleVector* > FissionEvent =
G4GenerateFission(1, Projectile);
120 std::vector< G4DynamicParticleVector* > FissionEvent =
G4GenerateFission(1, Projectile);
136 std::vector< G4DynamicParticleVector* > FissionEvents(0);
146 for(
G4long i = 0; i < NumberOfFissions; i++)
154 FissionEvents.push_back(NULL);
158 return FissionEvents;
258 return (A + Z * 1000) * 10 + M;
277 G4cout <<
" -- Alpha production set to " << AlphaProduction_ <<
G4endl;
299 G4cout <<
" -- Ternary fission probability set to " << TernaryProbability_ <<
G4endl;
314 if(!IsSameCause && IsValidCause)
330 CauseString =
"SPONTANEOUS";
333 CauseString =
"NEUTRON_INDUCED";
336 CauseString =
"PROTON_INDUCED";
339 CauseString =
"GAMMA_INDUCED";
352 G4cout <<
" -- Already set to use " << CauseString <<
" as the fission cause. Yield data class will not be reconstructed." <<
G4endl;
355 G4cout <<
" -- Yield data class not yet constructed. " << CauseString <<
" will be applied when it is constructed." <<
G4endl;
369 G4cout <<
" -- Fission cause set to " << CauseString <<
"." <<
G4endl;
392 std::ostringstream EnergyString;
415 G4cout <<
" -- Cannot set a non-zero energy for spontaneous fission" <<
G4endl;
421 G4cout <<
" -- Yield data class not yet constructed. " << EnergyString.str() <<
" will be applied when it is constructed." <<
G4endl;
432 G4cout <<
" -- Incident neutron energy set to " << EnergyString.str() <<
"." <<
G4endl;
461 G4cout <<
" -- Isotope " <<
Isotope_ <<
" already in use. Yield data class will not be reconstructed." <<
G4endl;
467 G4cout <<
" -- Yield data class not yet constructed. The isotope will be set to " <<
Isotope_ <<
" when it is constructed." <<
G4endl;
492 if(!IsSameMetaState && IsValidMetaState)
504 MetaName =
"GROUND_STATE";
521 std::ostringstream Temp;
526 G4cout <<
" -- Already set to use " << MetaName <<
" as the metastable state. Yield data class will not be reconstructed" <<
G4endl;
529 G4cout <<
" -- Yield data class not yet constructed. " << MetaName <<
" will be applied when it is constructed." <<
G4endl;
543 G4cout <<
" -- Metastable state set to " << MetaName <<
"." <<
G4endl;
557 G4bool IsSameScheme = (NewScheme == SamplingScheme_);
559 if(!IsSameScheme && IsValidScheme)
561 SamplingScheme_ = NewScheme;
562 IsReconstructionNeeded_ =
TRUE;
568 switch(SamplingScheme_)
571 SchemeString =
"NORMAL";
575 SchemeString =
"LIGHT_FRAGMENT";
579 SchemeString =
"UNSUPPORTED";
590 if(IsSameScheme && YieldData_ != NULL)
592 G4cout <<
" -- Already set to use " << SchemeString <<
" as the sampling scheme. Yield data class will not be reconstructed." <<
G4endl;
593 }
else if(YieldData_ == NULL)
595 G4cout <<
" -- Yield data class not yet constructed. " << SchemeString <<
" will be applied when it is constructed." <<
G4endl;
609 G4cout <<
" -- Sampling scheme set to " << SchemeString <<
"." <<
G4endl;
625 if(!IsSameYieldType && IsValidYieldType)
637 YieldString =
"INDEPENDENT";
641 YieldString =
"SPONTANEOUS";
645 YieldString =
"UNSUPPORTED";
661 G4cout <<
" -- Yield data class not yet constructed. Yield type " << YieldString <<
" will be applied when it is constructed." <<
G4endl;
675 G4cout <<
" -- Yield type set to " << YieldString <<
G4endl;
711 G4cout <<
" -- Old yield data class deleted." <<
G4endl;
746 G4cout <<
" -- Yield data class constructed with defined values." <<
G4endl;
748 }
catch (std::exception& e)
static const G4double AlphaProduction
G4FFGEnumerations::MetaState MetaState_
G4bool IsReconstructionNeeded_
void G4SetSamplingScheme(G4FFGEnumerations::FissionSamplingScheme NewScheme)
CLHEP::Hep3Vector G4ThreeVector
static const G4FFGEnumerations::FissionSamplingScheme SamplingScheme
static constexpr double keV
void G4SetTernaryProbability(G4double TernaryProbability)
G4FFGEnumerations::MetaState G4GetMetaState(void)
void G4SetYieldType(G4FFGEnumerations::YieldType WhichYieldType)
void G4SetMetaState(G4FFGEnumerations::MetaState WhichMetaState)
G4DynamicParticleVector * G4GenerateFission(void)
G4DynamicParticleVector * G4GetFission(void)
G4FissionProductYieldDist * YieldData_
void G4SetAlphaProduction(G4double WhatAlphaProduction)
static const G4int MetaStateLast
void G4SetVerbosity(G4int WhatVerbosity)
static G4int G4MakeIsotopeCode(G4int Z, G4int A, G4int M)
void G4SetVerbosity(G4int WhatVerbosity)
G4Ions * G4GenerateFissionProduct(void)
bool InitializeFissionProductYieldClass(std::istringstream &dataFile)
static const G4int MetaStateFirst
G4GLOB_DLL std::ostream G4cout
G4double TernaryProbability_
double A(double temperature)
const G4ParticleDefinition * GetDefinition() const
static const G4double TernaryProbability
static constexpr double MeV
static constexpr double eV
std::vector< G4DynamicParticle * > G4DynamicParticleVector
void G4SetIsotope(G4int WhichIsotope)
void G4SetEnergy(G4double WhatIncidentEnergy)
G4FFGEnumerations::YieldType YieldType_
static G4Neutron * Neutron()
G4FFGEnumerations::YieldType G4GetYieldType(void)
static const G4int FissionSamplingSchemeFirst
static const G4int Isotope
static constexpr double eV
void G4SetAlphaProduction(G4double WhatAlphaProduction)
static constexpr double GeV
void G4SetCause(G4FFGEnumerations::FissionCause WhichCause)
G4double AlphaProduction_
void G4SetIncidentEnergy(G4double WhatIncidentEnergy)
G4FFGEnumerations::FissionSamplingScheme SamplingScheme_
static const G4FFGEnumerations::MetaState MetaState
static constexpr double GeV
#define G4FFG_FUNCTIONLEAVE__
G4FFGEnumerations::FissionCause G4GetCause(void)
static constexpr double MeV
static const G4FFGEnumerations::YieldType YieldType
G4double G4GetIncidentEnergy(void)
G4double G4GetTernaryProbability(void)
void G4SetTernaryProbability(G4double WhatTernaryProbability)
static const G4double ThermalNeutronEnergy
G4FFGEnumerations::FissionSamplingScheme G4GetSamplingScheme(void)
static constexpr double keV
G4Ions * G4GetFissionProduct(void)
#define G4FFG_FUNCTIONENTER__
G4double G4GetAlphaProduction(void)
static const G4FFGEnumerations::Verbosity Verbosity
static const G4FFGEnumerations::FissionCause FissionCause
G4FissionFragmentGenerator(void)
G4FFGEnumerations::FissionCause Cause_
~G4FissionFragmentGenerator()
static const G4int FissionSamplingSchemeLast