Geant4中如何设置各向同性的源

	// ge xiang tong xing
	G4double Position[3];
	G4double R = 0.5 * micrometer;
	
	G4double th = acos(1.0-2*G4UniformRand());
	G4double ph = G4UniformRand()*2.0*CLHEP::pi;

	Position[0] = sin(th)*cos(ph);
	Position[1] = sin(th)*sin(ph);
	Position[2] = cos(th);

	G4ThreeVector positionPri = G4ThreeVector(Position[0],Position[1],Position[2]);
	
	//define particle position
	fpParticleGun->SetParticlePosition(R*positionPri);
	
	//get particle momentum direction
	G4double Direction[3];
	
	G4double theta = acos(1.0-G4UniformRand());
	G4double phi = G4UniformRand()*2.0*CLHEP::pi;

	Direction[0] = sin(theta)*cos(phi);
	Direction[1] = sin(theta)*sin(phi);
	Direction[2] = cos(theta);

	G4ThreeVector directPri = G4ThreeVector(Direction[0],Direction[1],Direction[2]);
	directPri.rotateUz(-positionPri);

	//define particle position
	fpParticleGun->SetParticleMomentumDirection(directPri);
	fpParticleGun->GeneratePrimaryVertex(anEvent);