56 return G4double((N-4)*(P-3)*(N-3)*(P-2))*
G4double((N-2)*(P-1)*(N-1)*P)/144.0;
89 if(nCharged >=2 && (nParticles-nCharged) >=2 ) {
91 G4double(nParticles*(nParticles-1)*(nParticles-2)*(nParticles-3));
92 rj = 6.0*nCharged*(nCharged-1)*(nParticles-nCharged)*(nParticles-nCharged-1)