Antihydrogen beam simulation

Particles, Track, Event, Run and Biasing
1

Geant4 Version: 11.1.2
Operating System: Linux Server
Compiler/Version: gcc-11.3.0
CMake Version:

Hi!
I am trying to simulate a beam of antihydrogen, but as far as I can see neither hydrogen or antihydrogen is available in the particle gun list. How would I have to define it?

Thanks in advance,
Viktoria

2

Well, the particle list includes proton, but that’s just the nucleus of hydrogen (and it will always have charge +1). For a neutral hydrogen beam, you would use

/gun/particle ion
/gun/ion 1 0 0      # Z A Q

while for an H- beam (like from a tandem), you would use

/gun/particle ion
/gun/ion 1 0 -1     # Z A Q

I don’t believe that Geant4 supports antihydrogen as a G4Ions particle, just anti_protons. It does support beams of low-mass anti-nuclei (use /gun/list to see the full list of particles), like anti_deuteron, anti_He3.

3

Thank you @mkelsey !
How would one go about defining the particle and then adding it to the G4ParticleTable? Or is that even possible?

4

As part of your application, you would write a G4ParticleDefinition subclass for the particle of interest. You would probably need to write processes for that particle’s unique interactions. Then in your PhysicsList you would need to add those processes to the particle’s ProcessManager.

Take a look at examples/extended/exoticphysics/monopole/.