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

[Pawe Kowalski]

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: GateContrib-master/imaging/PhantomSource/Analytical_Phantom/NEMA_IEC_2001.macUnfortunatelly the phantom seems to be useless without the dedicated source.
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 years).
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/20180213/f78700d5/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NEMA_IEC_2001_IQ_Phantom.mac
Type: application/octet-stream
Size: 19930 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20180213/f78700d5/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NEMA_IEC_2001_IQ_Source_Gamma.mac
Type: application/octet-stream
Size: 8929 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20180213/f78700d5/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: source_test.mac
Type: application/octet-stream
Size: 727 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20180213/f78700d5/attachment-0005.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sourcePosX_vs_sourcePosY.png
Type: image/png
Size: 39691 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20180213/f78700d5/attachment-0001.png>