I want to simulate Cs-137 through the cables lying near Ge-detectors and looking at its energy deposition in Germanium detectors. Basically i’m interested in it’s 662 keV gamma line. I did something like this:
/gun/particle ion
/gun/ion 55 137 0 0
/gun/energy 0 keV
But the problem is that the energy spectrum i got at the end has no 662 keV gamma line at all, although the branching ratio for it quite high (94.7 %).
I got the one at 284 keV but not at 662 keV. I don’t know what is going on. It would be really nice if someone could help me with this.
Here, a macro for example rdecay01, its printout, and the gamma plot.
May I suggest you to run your application for 1 or 2 events, with /tracking/verbose 2, in order to follow what happens to the 662 keV gamma.
Thank you very much for your response. I tried like you said. It worked like this. But i am using MaGe simulation framework for my work. I am trying to incident Cs137 from detector holder to detectors and finding its activity of 662 keV gamma line.
In my case, it seems that it’s not producing daughter nucleus properly so i decided to do something like this. I did one simulation using Cs137 at energy 0 keV and ran other simulation using Excited state of Ba-137 with energy of 662 keV so that i can also get 662 keV gamma line . And then i tried to add energy spectrum of both. But in this case Branching ratio i get for 284 keV(5.30%) and 662 keV(94.7%) are not the way it should be. Figure is attached below.
What is MaGe ?
Are you sure that radioactiveDecay is registered ? it is abnormal that Cs137 do not decay …
What do you plot ? energy of what ?
What is the origin of the 284 keV ray ?