59 std::vector< G4QMDParticipant* >::iterator
it;
61 p += (*it)->Get4Momentum();
72 std::vector< G4QMDParticipant* >::iterator
it;
92 std::vector< G4QMDParticipant* >::iterator
it;
125 - Asym * ( N -
Z )* ( N - Z ) /
A;
151 for (
G4int i= 0; i <
n ; i++ )
155 G4double trans = gamma / ( gamma + 1.0 ) * p_i * beta;
156 pcm[i] = p_i - trans*beta;
161 pcm0 = pcm0 / double ( n );
165 for (
G4int i= 0; i <
n ; i++ )
177 for (
G4int i= 0; i <
n ; i++ )
180 G4double trans = gamma / ( gamma + 1.0 ) * ri * beta;
184 rcm[i] = ri + trans*beta;
186 rcm0 +=
rcm[i]*
es[i];
193 for (
G4int i= 0; i <
n ; i++ )
202 for (
G4int i= 0; i <
n ; i++ )
204 rl +=
rcm[i].cross (
pcm[i] );
207 jj =
int ( std::sqrt ( rl*rl /
hbc ) + 0.5 );
213 for (
G4int i= 0; i <
n ; i++ )