Geant 4 射线源定义（初学笔记）

一、面源

源代码可参考basic/B1，表面 z=z0（x0/y0为一定范围内的随机数）

注：如果未指定方向，将沿（1., 0., 0.）即 X 轴

  // default particle kinematic
  // gamma particle
  G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
  G4String particleName;
  G4ParticleDefinition* particle
    = particleTable->FindParticle(particleName="gamma");
  fParticleGun->SetParticleDefinition(particle);

  // energy
  fParticleGun->SetParticleEnergy(6.*MeV);

  // position
  G4double x0 = G4UniformRand()-0.5;
  G4double y0 = G4UniformRand()-0.5;
  G4double z0 = -0.5;
  
  fParticleGun->SetParticlePosition(G4ThreeVector(x0,y0,z0));

  // direction
  fParticleGun->SetParticleMomentumDirection(G4ThreeVector(0.,0.,1.));

  fParticleGun->GeneratePrimaryVertex(anEvent);

二、点源

1. 各向同性：thetaMax = 180;

2. Co-60源：1.17/1.33 MeV；锥束的圆锥半角：21度。

  // default particle kinematic
  // gamma partical
  G4ParticleDefinition* particleDefinition 
    = G4ParticleTable::GetParticleTable()->FindParticle("gamma");
  fParticleGun->SetParticleDefinition(particleDefinition);

  // energy: Co-60 source
  G4double Flag = G4UniformRand();
  if(Flag<0.5){
      fParticleGun->SetParticleEnergy(1173.2*keV);}
  else{
      fParticleGun->SetParticleEnergy(1332.5*keV);
  }

  // position
  fParticleGun->SetParticlePosition(G4ThreeVector(0,0,-14*cm));

  // direction
  G4double phi = G4UniformRand()*360.*deg;                                                       
  G4double thetaMax = G4UniformRand()*21*deg;
  G4double mu = std::cos(thetaMax);
  
  G4double v1=std::sqrt(1-mu*mu)*std::cos(phi) ;
  G4double v2=std::sqrt(1-mu*mu)*std::sin(phi) ;
  G4double v3=mu;

  fParticleGun->SetParticleMomentumDirection(G4ThreeVector(v1,v2,v3));

  fParticleGun->GeneratePrimaryVertex(anEvent);

方向定义部分等同于：

  // direction
  G4double phi = G4UniformRand()*360.*deg;                                                       
  G4double thetaMax=21*deg;
  G4double muMin=std::cos(thetaMax);
  G4double mu=muMin+G4UniformRand()*(1.0-muMin);