65 if (MaximalKineticEnergy <= 0.0) {
return 0.0; }
72 if(Ucompound < 0.0 || Ufission < 0.0) {
return 0.0; }
79 G4double Cf = 2.0*std::sqrt(afission*MaximalKineticEnergy);
80 G4double Exp1 = (SystemEntropy <= 160.0) ?
G4Exp(-SystemEntropy) : 0.0;
81 G4double Exp2 = (SystemEntropy-Cf <= 160.0) ?
G4Exp(-SystemEntropy+Cf) : 0.0;
85 G4double probability = (Exp1 + (Cf-1.0)*Exp2) / (4.0*
pi*afission);