Loading...
Searching...
No Matches
G4FPYBiasedLightFragmentDist Class Reference
#include <G4FPYBiasedLightFragmentDist.hh>
Inheritance diagram for G4FPYBiasedLightFragmentDist:
Detailed Description
G4FPYBiasedLightFragmentDist is an inherited class of G4FissionProductYield that only samples the 'light' fission fragments, defined by A < 115
- This inherited class of G4FissionProductYield samples only the lighter fission fragments, defined by A < 115
- This biasing was implemented because of an artifact that is introduced due to random sampling of fission fragments. Typically small fission fragments (neutrons, alphas, gammas) are sampled after the first fragment is sampled. If a heavy fragment (A >= 115) is sampled first then the resulting lighter fission fragment, after all the other particles have been removed from the available mass, will most likely land far off the neutron drip line.
- This implementation reduces the probability that such an improbable nucleus will be created by first sampling the lighter fission fragment and allowing the heavy fission fragment, which has a lot more flexibility for varying neutron populations, to make up the slack.
Definition at line 58 of file G4FPYBiasedLightFragmentDist.hh.
Constructor & Destructor Documentation
◆ G4FPYBiasedLightFragmentDist() [1/2]
| G4FPYBiasedLightFragmentDist::G4FPYBiasedLightFragmentDist | ( | G4int | WhichIsotope, |
| G4FFGEnumerations::MetaState | WhichMetaState, | ||
| G4FFGEnumerations::FissionCause | WhichCause, | ||
| G4FFGEnumerations::YieldType | WhichYieldType, | ||
| std::istringstream & | dataStream | ||
| ) |
Default constructor
- Usage:
WhichIsotope:Isotope number of the element in ZZZAAA formWhichMetaState:GROUND_STATE,META_1, orMETA_2WhichCause:SPONTANEOUSorN_INDUCEDWhichYieldType:INDEPENDENTorCUMULATIVE
- Notes:
Definition at line 42 of file G4FPYBiasedLightFragmentDist.cc.
48: G4FissionProductYieldDist( WhichIsotope,
49 WhichMetaState,
50 WhichCause,
51 WhichYieldType,
52 dataStream)
53{
54 // Initialize the class
55 Initialize();
56}
void Initialize(void)
Definition: G4FPYBiasedLightFragmentDist.cc:77
◆ G4FPYBiasedLightFragmentDist() [2/2]
| G4FPYBiasedLightFragmentDist::G4FPYBiasedLightFragmentDist | ( | G4int | WhichIsotope, |
| G4FFGEnumerations::MetaState | WhichMetaState, | ||
| G4FFGEnumerations::FissionCause | WhichCause, | ||
| G4FFGEnumerations::YieldType | WhichYieldType, | ||
| G4int | Verbosity, | ||
| std::istringstream & | dataStream | ||
| ) |
Overloaded constructor
- Usage:
WhichIsotope:Isotope number of the element in ZZZAAA formWhichMetaState:GROUND_STATE,META_1, orMETA_2WhichCause:SPONTANEOUSorN_INDUCEDWhichYieldType:INDEPENDENTorCUMULATIVEVerbosity:Verbosity level
- Notes:
Definition at line 58 of file G4FPYBiasedLightFragmentDist.cc.
65: G4FissionProductYieldDist( WhichIsotope,
66 WhichMetaState,
67 WhichCause,
68 WhichYieldType,
69 Verbosity,
70 dataStream )
71{
72 // Initialize the class
73 Initialize();
74}
◆ ~G4FPYBiasedLightFragmentDist()
|
virtual |
Default deconstructor. It is a virtual function since G4FPYBiasedLightFragmentDist inherits from G4FissionProductYieldDist
Definition at line 130 of file G4FPYBiasedLightFragmentDist.cc.
131{
133
134 // Empty - all the data elements to be deconstructed are removed by
135 // ~G4FissionProductYieldDist()
137}
Member Function Documentation
◆ GetFissionProduct()
|
protectedvirtual |
Selects a fission product from the probability tree, limited by the number of nucleons available to the system
- Todo:
- Figure out what to do if a suitable particle is not found in 1000 iterations
Implements G4FissionProductYieldDist.
Definition at line 87 of file G4FPYBiasedLightFragmentDist.cc.
89{
91
92 G4Ions* Particle;
93 G4bool IsHeavy;
94 G4bool IsNotFeasable;
95 G4int Counter;
96
97 Counter = 0;
98 do
99 {
100 // If we have sampled too many times then we may need to do something else
101 if(Counter == 1000)
102 {
103 Particle = NULL;
104 break;
105 }
106
107 // Generate a (0, 1] random number, then find the corresponding particle
109
110 // If the daughter is heavier than half of the original nucleus weight
111 // or not physically feasable then set the flags to sample again
112 IsHeavy = (Particle->GetAtomicMass() > HalfWeight_);
113 IsNotFeasable = (Particle->GetAtomicMass() > RemainingA_ + 1
114 || Particle->GetAtomicNumber() > RemainingZ_ + 1);
115
116 Counter++;
117 } while (IsHeavy || IsNotFeasable); // Loop checking, 11.05.2015, T. Koi
118
119 // Something went wrong, so figure out how to fix it
120 if(Particle == NULL)
121 {
122 /// \todo Figure out what to do if a suitable particle is not found in 1000 iterations
123 }
124
125 // Return the G4ParticleDefintion pointer to the random light fragment
127 return Particle;
128}
G4int HalfWeight_
Definition: G4FPYBiasedLightFragmentDist.hh:100
G4FPYSamplingOps * RandomEngine_
Definition: G4FissionProductYieldDist.hh:261
G4Ions * FindParticle(G4double RandomParticle)
Definition: G4FissionProductYieldDist.cc:687
Definition: G4Ions.hh:52
◆ Initialize()
|
protected |
Initialize is a common function called by all constructors.
Definition at line 76 of file G4FPYBiasedLightFragmentDist.cc.
78{
80
81 // Initialize the half-weight for fission product sampling
83
85}
Referenced by G4FPYBiasedLightFragmentDist().
Member Data Documentation
◆ HalfWeight_
|
protected |
Defines the half-weight of the fission isotope
Definition at line 100 of file G4FPYBiasedLightFragmentDist.hh.
Referenced by GetFissionProduct(), and Initialize().
The documentation for this class was generated from the following files:
