[Gate-users] Implementation of the NEMA Image Quality source

Ane Etxebeste fplico817 at gmail.com
Wed Feb 14 11:30:18 CET 2018

Dear Pawe Kowalski,

I have not looked carefully into your macro.
But I have seen that you have applied Forbid command several times to
the same volume without specifying it.

For example, you have the following lines in your macro.
 /gate/source/IBox_hot/gps/Forbid              sphere22
 /gate/source/IBox_hot/gps/Forbid              sphere22in
 /gate/source/IBox_hot/gps/Forbid              fill22
 /gate/source/IBox_hot/gps/Forbid              fill22in
 /gate/source/IBox_hot/gps/Forbid              sphere28
 /gate/source/IBox_hot/gps/Forbid              sphere28in
 /gate/source/IBox_hot/gps/Forbid              fill28
 /gate/source/IBox_hot/gps/Forbid              fill28in

In this case you have applied Forbid command 8 times to IBox_hot volume.
Try adding the following command to specify the number of cold volumes
you want to create inside IBox_hot volume
 /gate/source/IBox_hot/dump 8

I hope this helps you.
Best regards,
Ane Etxebeste

Quoting "Pawe  Kowalski" <pawkowalski at yahoo.com>:

Dear GATE users,
I am currently working on simulations of the NEMA characteristics of
a PET scanner.
While it is quite easy to implement phantoms and sources for spatial
resolution, scatter fraction and sensitivity characteristics,  it is
challenging to implement the NEMA Image Quality phantom and the source.
Happily, the macro for the phantom may be found at Github repository
in path:
the phantom seems to be useless without the dedicated
I have searched for many archive messeges from gate-users mailing list and
I have found no working macro with implementation of this source. It is
very strange for me because it seems that the Image Quality is a quite
important characteristic and many people try to perform such simulations.
Has noone really managed to implement the source for the NEMA IQ phantom?
If someone managed to do it, please make it public.
Using some examples that may be found in the network, I have been
trying to implement this source by myself. Unfortunatelly I have
some problems that I cannot solve (and it seems that they were not
solved by the people that tried to implement the phantom in previous
In the attachments there are my macros for the phantom, the source,
and the test source which is a part of the whole source.
Because of the fact that it is not possible to use gps source with volume
which is an angular part of the cylinder (for the phantom there are
commands: setPhiStart and setDeltaPhi) I have tried to implement parts of
the IQ source using cylindrical sources and cold regions. Unortunatelly, it
seems that the Forbid command does not work as it should.There are some
generations in regions, in which theoretically they should not appear.
The macro source_test.mac is a part of the source. The generations
should be only between angles 270 and 360 but they are also in the
rest of the region. I attach 2d histogram of sourcePosX and
sourcePosY as a proof (for simulation with source_test.mac).
What am I doing wrong? Please take a look at my macros. Maybe
someone will know how to correct them in order to make them proper.
I would be very grateful for your help and for any suggestions that
would lead to the solution.
Best regards,Pawe? Kowalski,National Centre for Nuclear Research
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20180214/3cf9803f/attachment.html>

More information about the Gate-users mailing list