69 outFile <<
"G4LFission is one of the Low Energy Parameterized\n"
70 <<
"(LEP) models used to implement neutron-induced fission of\n"
71 <<
"nuclei. It is a re-engineered version of the GHEISHA code\n"
72 <<
"of H. Fesefeldt which emits neutrons and gammas but no\n"
73 <<
"nuclear fragments. The model is applicable to all incident\n"
74 <<
"neutron energies.\n";
83 for (i = 2; i <= 10; i++) {
85 xxx = std::sqrt(2.29*xx);
88 for (i = 1; i <= 10; i++) {
114 G4cout <<
"G4LFission:ApplyYourself: incident particle:" <<
G4endl;
125 G4cout <<
"G4LFission:ApplyYourself: material:" <<
G4endl;
128 G4cout <<
"atomic mass " <<
134 if (E02 < 0) E0 = -E0;
147 if (e1 < 1.) e1 = 1.;
153 if (!photofission) avern = 2.569 + 0.900*
G4Log(e1);
163 G4int ng =
static_cast<G4int>(averg + ran*3. + 0.5);
172 for (i = 1; i <=
nn; i++) {
175 for (j = 1; j <= 10; j++) {
176 if (ran <
spneut[j-1])
goto label12;
191 for (i = 1; i <= ng; i++) {
205 for (i = 1; i <= nn + ng; i++) {
215 G4double px = pp*sint*std::sin(phi);
216 G4double py = pp*sint*std::cos(phi);
224 a = (a/(E + E0) - e)/E0;
229 G4double p2 = px*px + py*py + pz*pz;
231 e = std::sqrt(e0*e0 + p2);
254 if (ia < 1)
return 0;
256 if (iz < 0)
return 0;
257 if (iz > ia)
return 0;
260 if (iz == 0)
return rmn;
261 if (iz == 1)
return rmp + rmel;
263 else if (ia == 2 && iz == 1) {
266 else if (ia == 4 && iz == 2) {
273 + 93.15*(A/2. -
Z)*(A/2. - Z)/A
274 + 0.6984523*Z*Z/Pow->
A13(A);
275 G4int ipp = (ia - iz)%2;
277 if (ipp == izz) mass = mass + (ipp + izz -1)*12.*Pow->
powA(A, -0.5);
285 return std::pair<G4double, G4double>(5*
perCent,250*
GeV);