Get a Gaussian Energy deposited Spectrum

So is there any other way to get a Gaussian kind of spectrum while decaying a source? Like tracking the energy deposited by secondaries something like that…

For example, how did they get such kind of spectrum in the following?