41 #define MAX_SECONDARIES 100
139 G4cout <<
"G4AntiNeutronAnnihilationAtRestProcess::AtRestGetPhysicalInteractionLength ";
176 for (
G4int i1=0; i1 < numberOfElements; i1++ )
178 normalization += theAtomicNumberDensity[i1] ;
183 for (
G4int i2=0; i2 < numberOfElements; i2++ )
185 runningSum += theAtomicNumberDensity[i2];
187 if (random<=runningSum)
193 if (random>runningSum)
200 G4cout <<
"G4AntiNeutronAnnihilationAtRest::AtRestDoIt is invoked " <<
G4endl;
219 G4Track* aNewTrack =
new G4Track( aNewParticle, localtime*
s, position );
280 for (l = 1; l <=
ntot; ++l) {
321 ran1 = xav + ran1 * std::sqrt(xav);
335 for (i = 1; i <= fivex; ++i) {
396 for (i = 2; i <= j; ++i) {
428 G4float ran1, ran2, ekin, tkin;
485 en = ek + (rmnve1 + rmnve2) /
G4float(2.);
486 r__1 = en * en - rmnve1 * rmnve2;
487 pcm = r__1 > 0 ? std::sqrt(r__1) : 0;
527 black = (targ *
G4float(1.25) +
539 for (i = 1; i <= nbl; ++i) {
556 ekin1 = tex - (ekin2 - ekin1);
584 for (i = 1; i <= nt; ++i) {
614 for (i = 1; i <= nbl; ++i) {
631 ekin1 = tex - (ekin2 - ekin1);
658 for (i = 3; i <= nt; ++i) {
672 G4float cfa, gfa, ran1, ran2, ekin1, atno3;
696 cfa =
G4float(.13043478260869565);
725 evapEnergy1 *= ran1 * gfa +
G4float(1.);
726 if (evapEnergy1 < 0.0
f) {
735 while ((ret_val = evapEnergy1 +
evapEnergy3) >= ek1) {