Hello @anwaets,
Thank you for your remarks on the post Radiative processes included when using the PAIRBREM card. Indeed, the photons I get in my simulation are not originating from proton bremsstrahlung, since I am injecting protons with a power-law energy distribution in the range from 1 MeV to 1 GeV, i.e. with energies below the threshold of 41.9 GeV for protons in H (Radiative processes included when using the PAIRBREM card - #3 by cesc). The photons are also not originating from bremsstrahlung of delta-rays and secondary electrons and positrons produced by other processes, since they are discarded in my simulation (see the attached test.inp file).
As you mentioned, the photons can also be generated by inelastic interactions. Following your suggestion, I will try to run a simulation for power-law protons with energies above the 41.9 GeV threshold and use the mgdraw.f routine to check which processes are generating the photons (any hints are welcomed !)
Cheers,
Sergio
test.flair (3.3 KB)
test.inp (2.8 KB)
source.f (15.5 KB)