46 theFS.
Init(A, Z, M, dirName, aFSType, projectile);
47 theFC.
Init(A, Z, M, dirName, aFSType, projectile);
48 theSC.
Init(A, Z, M, dirName, aFSType, projectile);
49 theTC.
Init(A, Z, M, dirName, aFSType, projectile);
50 theLC.
Init(A, Z, M, dirName, aFSType, projectile);
52 theFF.
Init(A, Z, M, dirName, aFSType, projectile);
55 G4cout <<
"Fission fragment production is now activated in HP package for "
57 <<
", A = " << (
G4int)A
60 G4cout <<
"As currently modeled this option precludes production of delayed neutrons from fission fragments." <<
G4endl;
111 theNeutron.
Lorentz(theNeutron, -1*theTarget);
138 if(random<xSec[i]/xSec[3])
break;
146 G4int Prompt=0, delayed=0, all=0;
152 if(Prompt==0&&delayed==0) Prompt=all;
158 if(Prompt==0&&delayed==0) Prompt=all;
163 if(Prompt==0&&delayed==0) Prompt=all;
168 if(Prompt==0&&delayed==0) Prompt=all;
183 if( theNeutrons != 0)
185 theDecayConstants =
new G4double[delayed];
190 for(i=0; i<theNeutrons->size(); i++)
199 for(i=0; i<theDelayed->size(); i++)
212 theDecayConstants =
new G4double[delayed];
213 if(Prompt==0&&delayed==0) Prompt=all;
219 for(i0=0; i0<Prompt; i0++)
225 for(i0=Prompt; i0<Prompt+delayed; i0++)
229 if ( theDecayConstants[i0-Prompt] > 1.0
e-30 ) {
233 ed <<
" theDecayConstants[i0-Prompt]=" << theDecayConstants[i0-Prompt]
234 <<
" -> cannot sample the time : set it to 0.0 !" <<
G4endl;
244 delete [] theDecayConstants;
246 unsigned int nPhotons = 0;
249 nPhotons = thePhotons->size();
250 for(i=0; i<nPhotons; i++)
299 G4ThreeVector direction(sinth*std::cos(phi), sinth*std::sin(phi), costheta);