I have simulated 2 * 10^8, Cs-137 ions from cables near Ge-detectors and looked into the energy depositions by Cs-137 into germanium detectors. The spectrum shown below is what I got. Now spectrum looks familiar and there is also a clear peak at 662 keV energy but I am confused about the broad structure(broad peak-like, marked in red) after 662 keV energy. Can anyone explain to me what this structure represents? Or am I doing something wrong and this shouldn’t be present here?
After reviewing your figure, I want to comment that your spectrum looks good. Perhaps you need to add resolution of the detector in your simulations to make your simulation results valid. This can be done in the EventAction.CC or in the post simulations process.
Red part - > Compare the events with the photo-peak, this is considered as background.
You should first add resolution. As i have no idea about your code therefore, i consider it as background. Also, if you want it on a log scale,
why not scale from 10 to 10^8. Why from 1 or from 0.001.?
Without knowing your geometry, it is difficult to say for certain. However, in the decay of Cs137, there are two betas emitted - a 0.5120 MeV endpoint energy beta with 94.6% probability and a 1.174 MeV beta at 5.4% probability. The shape that you’ve circled looks somewhat like a beta spectrum. Its max energy and energy of max intensity are not inconsistent with the high energy beta. So it could be due to the latter one. Of course, energy deposition of the betas is heavily dependent on the geometry.
One test you could do is to replace the Cs-137 source with a source having only the 0.662 MeV gamma ray and see what the spectrum looks like.
Thank you. I am going to read this. One question though: Giving the mean and sigma would give resolution to only one gamma peak right? I mean i am simulating whole Cs-137 isotope and if i will implement this G4RandGauss::shoot(eDep,eDep_sigma); for 662 keV gama, it will smear the peak at 662 keV only, right?