Distribution uniforme en 4pi

Bonjour,


J'ai une question concernant les distributions uniformes. Voila je suis entrain de programmer en c++ et geant4 (pour faire des simulations d'interaction de particules au travers de la matière). Et je voudrais avoir une distribution uniforme en 4pi
voici ma distribution en 2pi :

G4double phi = 2*pi*G4UniformRand();

G4double cosTheta = G4UniformRand();
G4double sinTheta = std::sqrt(1. - cosTheta*cosTheta);
G4double ux = sinTheta*std::cos(phi),
uy = sinTheta*std::sin(phi),
uz = -cosTheta;

comment faire en 4pi? si vous pouvez m'aider c'est gentil.