Essentially in my simulation, it appears that my 3MeV gamma ray under goes Compton scattering almost exclusively.
For every scattering position, GEANT4 does not however always produce an electron it appears. I’ve been thinking this is more down to my understanding of the physics, but if a photon undergoes compton scattering, shouldn’t an electron always be produced (assuming it is a loosely bound electron)?
I understand that in the energy region 2-4 MeV this corresponds to a region where the absorption and scattering probabilities for compton scattering are close together as in the picture below:
How do you have your production cuts set? If the secondary electron (from “ionization”, maybe not from Compton scattering) is too low in energy, it may have an estimated range in your target below your production cut. In that case, it will not be created as a “tracked particule”, but it’s energy will be assigned to the energy deposit of the photon’s step.
Thanks @mkelsey, so my electron has a production cut of 1m, so it shouldn’t be this I assume. I show an example event to illustrate the problem (gamma = green, electron=red, opticalphotons=grey)
Please, notice that, by default, production cut is applied only to electrons from ionization and photons from Bremstrahlung. Compton scattering itself in independent of production cut.
Below, a macro for example TestEm14, its printout, and the plot of e- spectrum.
It shows that Compton scattering works as expected, even with a cut of 1 meter.
There is an option (/process/em/applyCuts) to generalize production threshold to all electromagnetic processes. This option is not by default, and you should have a strong reason to use it.