How do I replace multiple Coulomb scattering by single scattering?
To activate single scattering everywhere, use the MULSOPT command with the following values:
WHAT(1): 0.0 WHAT(2): 0.0 WHAT(3): 0.0 WHAT(4): 1.0 WHAT(5): 1.0 WHAT(6): > 1000.0 SDUM: GLOBAL (GLOBEMF applies only to electrons and positrons, GLOBHAD only tho hadrons and muons)
Note that this choice is generally extremely demanding in terms of CPU time, except for particles of very low energy (a few keV), which have a very short history anyway. In such cases, the single scattering option is even recommended.