64 G4int G4fissionEvent::G4SmpSpNuDistData(
G4int isotope,
G4int Cf252option) {
101 {0.0481677,0.2485215,0.4253044,0.2284094,0.0423438,0.0072533,
102 0.0000000,0.0000000,0.0000000,0.0000000,0.0000000},
104 {0.0631852,0.2319644,0.3333230,0.2528207,0.0986461,0.0180199,
105 0.0020407,0.0000000,0.0000000,0.0000000,0.0000000},
107 {0.0679423,0.2293159,0.3341228,0.2475507,0.0996922,0.0182398,
108 0.0031364,0.0000000,0.0000000,0.0000000,0.0000000},
110 {0.0212550,0.1467407,0.3267531,0.3268277,0.1375090,0.0373815,
111 0.0025912,0.0007551,0.0001867,0.0000000,0.0000000},
113 {0.0150050,0.1161725,0.2998427,0.3331614,0.1837748,0.0429780,
114 0.0087914,0.0002744,0.0000000,0.0000000,0.0000000},
116 {0.0540647,0.2053880,0.3802279,0.2248483,0.1078646,0.0276366,
117 0.0000000,0.0000000,0.0000000,0.0000000,0.0000000},
119 {0.0021100,0.0246700,0.1229000,0.2714400,0.3076300,0.1877000,
120 0.0677000,0.0140600,0.0016700,0.0001000,0.0000000},
122 {0.0020900,0.0262100,0.1262000,0.2752000,0.3018000,0.1846000,
123 0.0668000,0.0150000,0.0021000,0.0000000,0.0000000} };
130 if (isotope == 92238) index = 0;
131 else if (isotope == 94240) index = 1;
132 else if (isotope == 94242) index = 2;
133 else if (isotope == 96242) index = 3;
134 else if (isotope == 96244) index = 4;
135 else if (isotope == 94238) index = 5;
136 else if (isotope == 98252 && Cf252option == 0) index = 6;
137 else if (isotope == 98252 && Cf252option == 1) index = 7;
144 sum = sum + sfnu[
index][i];
145 if (r <= sum || sfnu[index][i+1] == 0.)
return i;
150 G4cout <<
" Random number out of range in SmpSpNuDistData " <<
G4endl;
159 nubar = G4SmpSpNubarData(isotope);
161 return (
G4int) G4SmpTerrell(nubar);
G4GLOB_DLL std::ostream G4cout