// 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);
Geant4中如何设置各向同性的源
猜你喜欢
转载自blog.csdn.net/baidu_29950065/article/details/79098045
今日推荐
周排行