[Gate-users] how to get proportionate number of photons detection

Ashok Tiwari tiwarias at yahoo.com
Fri Apr 8 07:03:13 CEST 2022


Dear GATE users,
 I tried to simulate the I-123 photons spectrum (a) using only 159 keV photon (83.3%) in one simulation and (b) all other high energy photons greater than 159 keV (~3%) in another simulation. 
In both of my simulations, I simulated a fixed number of primaries (decays) so that I could estimate the contribution of 159 keV and high energy photons separately. However, by defining the sources (see below), a similar number of photons have been detected in both simulations! This is undesirable as the intensity of their emissions are hugely different, as defined in the lines below. It looks like GATE normalized the intensity value (probability of emission); as a result of that, I may end up getting similar results for both simulations. Now, my question is, how can I define the sources in (b) to get a proportionate number of photons detection? If my assumption is right about normalization, is it possible to force GATE not to normalize the intensity values in (b)?

For (a) I define the following lines of code in my simulation settings:
/gate/source/addSource  src1            gps
/gate/source/src1/gps/particle            gamma/gate/source/src1/gps/energytype      Mono/gate/source/src1/gps/ene/mono        158.97 keV
/gate/source/src1/setIntensity             0.833/gate/source/src1/gps/angtype           iso/gate/source/src1/gps/type                 Point/gate/source/src1/gps/centre              0. 0. 0. mm/gate/source/src1/attachTo                 src_1/gate/source/list

For (b) I selected 5 major gamma emissions having energy greater than 159 keV:/gate/source/addSource  src1          gps/gate/source/src1/gps/particle          gamma/gate/source/src1/gps/energytype    Mono/gate/source/src1/gps/ene/mono      346.35 keV
/gate/source/src1/setIntensity           0.00125783/gate/source/src1/gps/angtype          iso/gate/source/src1/gps/type                Point/gate/source/src1/gps/centre             0. 0. 0. mm/gate/source/src1/attachTo                src_1
/gate/source/addSource  src2          gps/gate/source/src2/gps/particle          gamma/gate/source/src2/gps/energytype    Mono/gate/source/src2/gps/ene/mono      440.02 keV
/gate/source/src2/setIntensity          0.00428162/gate/source/src2/gps/angtype         iso/gate/source/src2/gps/type               Point/gate/source/src2/gps/centre            0. 0. 0. mm/gate/source/src2/attachTo               src_1
/gate/source/addSource  src3          gps/gate/source/src3/gps/particle          gamma/gate/source/src3/gps/energytype    Mono/gate/source/src3/gps/ene/mono      505.33 keV
/gate/source/src3/setIntensity          0.00315707/gate/source/src3/gps/angtype         iso/gate/source/src3/gps/type               Point/gate/source/src3/gps/centre            0. 0. 0. mm/gate/source/src3/attachTo               src_1
/gate/source/addSource  src4          gps/gate/source/src4/gps/particle          gamma/gate/source/src4/gps/energytype    Mono/gate/source/src4/gps/ene/mono      528.96 keV
/gate/source/src4/setIntensity           0.013911/gate/source/src4/gps/angtype         iso/gate/source/src4/gps/type               Point/gate/source/src4/gps/centre            0. 0. 0. mm/gate/source/src4/attachTo               src_1
/gate/source/addSource  src5          gps/gate/source/src5/gps/particle          gamma/gate/source/src5/gps/energytype    Mono/gate/source/src5/gps/ene/mono      538.54 keV
/gate/source/src5/setIntensity           0.00381514/gate/source/src5/gps/angtype         iso/gate/source/src5/gps/type               Point/gate/source/src5/gps/centre            0. 0. 0. mm/gate/source/src5/attachTo               src_1
Any help is highly appreciated. Thanks in advance.
Regards,Ashok 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20220408/267776b1/attachment.html>


More information about the Gate-users mailing list