80 if (anA >= 65 && ExEnergy > 0.0) {
108 if (U <= pcorr) {
return Fragment1; }
129 G4double FragmentsExcitationEnergy = 0.0;
130 G4double FragmentsKineticEnergy = 0.0;
143 if (A2 < 1 || Z2 < 0 || Z2 > A2) {
144 FragmentsExcitationEnergy = -1.0;
149 G4double Tmax = M + U - M1 - M2 - pcorr;
153 FragmentsExcitationEnergy = -1.0;
168 FragmentsExcitationEnergy =
170 Tmax - FragmentsKineticEnergy + pcorr;
173 }
while (FragmentsExcitationEnergy < 0.0
176 if (FragmentsExcitationEnergy <= 0.0) {
178 "G4CompetitiveFission::BreakItUp: Excitation energy for fragments < 0.0!");
182 M1 += FragmentsExcitationEnergy * A1/
static_cast<G4double>(
A);
184 M2 += FragmentsExcitationEnergy * A2/
static_cast<G4double>(
A);
188 G4double etot1 = ((M - M2)*(M + M2) + M1*M1)/(2*M);
195 Fragment1 =
new G4Fragment( A1, Z1, FourMomentum1);
196 theNucleusMomentum -= FourMomentum1;
219 if (w > 1000.0 ) { C2 = C2S; }
220 else if (w < 0.001) { C2 = C2A; }
240 if (Mass2 > MassMax) { MassMax = Mass2; }
241 if (Mass3 > MassMax) { MassMax = Mass3; }
242 if (Mass4 > MassMax) { MassMax = Mass4; }
243 if (Mass5 > MassMax) { MassMax = Mass5; }
271 + 0.5*(
G4Exp(-0.5*z1*z1) +
G4Exp(-0.5*z2*z2));
275 if (w > 1000) { res = Xsym; }
276 else if (w < 0.001) { res = Xasym; }
277 else { res = w*Xsym+Xasym; }
286 if (Af >= 134.0) { DeltaZ = -0.45; }
287 else if (Af <= (A-134.0)) { DeltaZ = 0.45; }
288 else { DeltaZ = -0.45*(Af-A*0.5)/(134.0-A*0.5); }
296 }
while (theZ < 1.0 || theZ > (Z-1.0) || theZ > Af);
316 Pas = 0.5*
G4Exp(-0.5*x1*x1) +
G4Exp(-0.5*x2*x2);
349 TaverageAfMax = (Eaverage + 12.5 * Xsy) * (PPas/ScaleFactor) *
355 TaverageAfMax = (Eaverage - 12.5*
CLHEP::MeV*Xas)
366 if (++i > 100)
return Eaverage;
368 }
while (KineticEnergy < Eaverage-3.72*ESigma ||
369 KineticEnergy > Eaverage+3.72*ESigma ||
370 KineticEnergy > Tmax);
372 return KineticEnergy;
G4double AsymmetricRatio(G4int A, G4double A11)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
T max(const T t1, const T t2)
brief Return the largest of the two arguments
const G4LorentzVector & GetMomentum() const
G4double GetFissionPairingCorrection(G4int A, G4int Z) const
G4double Z13(G4int Z) const
G4double GetSigma1(void) const
Float_t y1[n_points_granero]
G4double LevelDensityParameter
G4double GetSigma2(void) const
virtual G4double LevelDensityParameter(G4int A, G4int Z, G4double U) const =0
G4bool MyOwnFissionProbability
Float_t x1[n_points_granero]
void SetZandA_asInt(G4int Znew, G4int Anew)
static G4NuclearLevelData * GetInstance()
G4int FissionCharge(G4int A, G4int Z, G4double Af)
G4double MaximalKineticEnergy
G4ThreeVector G4RandomDirection()
G4double GetGroundStateMass() const
G4FissionParameters theParam
Float_t y2[n_points_geant4]
G4double GetSigmaS(void) const
static G4Pow * GetInstance()
G4VEmissionProbability * theFissionProbabilityPtr
G4double FissionKineticEnergy(G4int A, G4int Z, G4int Af1, G4int Zf1, G4int Af2, G4int Zf2, G4double U, G4double Tmax)
static constexpr double MeV
virtual ~G4CompetitiveFission()
double A(double temperature)
G4PairingCorrection * pairingCorrection
G4VLevelDensityParameter * theLevelDensityPtr
void SetMomentum(const G4LorentzVector &value)
static G4double GetNuclearMass(const G4double A, const G4double Z)
G4double GetAs(void) const
void DefineParameters(G4int A, G4int Z, G4double ExEnergy, G4double FissionBarrier)
G4VFissionBarrier * theFissionBarrierPtr
ThreeVector shoot(const G4int Ap, const G4int Af)
G4double FissionProbability
virtual G4double FissionBarrier(G4int A, G4int Z, G4double U) const =0
virtual G4double GetEmissionProbability(G4Fragment *theNucleus)
Hep3Vector boostVector() const
virtual G4Fragment * EmittedFragment(G4Fragment *theNucleus)
G4int FissionAtomicNumber(G4int A)
G4double MassDistribution(G4double x, G4int A)
G4PairingCorrection * GetPairingCorrection()
Float_t x2[n_points_geant4]
G4double GetExcitationEnergy() const
G4double GetW(void) const
G4double SymmetricRatio(G4int A, G4double A11)
HepLorentzVector & boost(double, double, double)
virtual G4double EmissionProbability(const G4Fragment &fragment, G4double anEnergy)=0
G4bool MyOwnFissionBarrier