ミューオンしたい

/gun/particle mu-
/gun/particle mu+

ミューオンは正ミューオン(mu+)、負ミューオン(mu-)が用意されています。

1G4ParticleTable *particle_table = G4ParticleTable::GetParticleTable();
2G4ParticleDefinition *muon_m = particle_table->FindParticle("mu-");
3G4ParticleDefinition *muon_p = particle_table->FindParticle("mu+");

宇宙線ミューオンしたい

 1PrimaryGeneratorAction::PrimaryGeneratorAction()
 2{
 3    G4int n_particle = 1;
 4    fParticleGun = new G4ParticleGun(n_particle);
 5
 6    G4ParticleTable *particle_table = new G4ParticleTable::GetParticleTable();
 7    G4ParticleDefinition *particle = particle_table->FindParticle("mu-");
 8    fParticleGun->SetParticleDefinition(particle);
 9    fParticleGun->SetParticlePosition(G4ThreeVector(0., 0., 0.));
10    fParticleGun->SetParticleEnergy(500*MeV);
11    fParticleGun->SetParticleMomentumDirection(G4ThreeVector(0., 1., 0.));
12}

宇宙線ミューオンなので、鉛直下向きに打ち込んでいます。