<html><head></head><body><div style="color:#000; background-color:#fff; font-family:Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:13px"><div id="yiv4949722788"><div id="yui_3_16_0_ym19_1_1518539686278_2411"><div style="color:#000;background-color:#fff;font-family:Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:13px;" id="yui_3_16_0_ym19_1_1518539686278_2410"><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627">Dear GATE users,</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627">I am currently working on simulations of the NEMA characteristics of a PET scanner.</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627">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.</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">Happily, the macro for the phantom may be found at Github repository in path: GateContrib-master/imaging/PhantomSource/Analytical_Phantom/NEMA_IEC_2001.mac</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627">Unfortunatelly the phantom seems to be useless without the dedicated source.</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">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. </div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">Has noone really managed to implement the source for the NEMA IQ phantom?</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr"><b id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_24456">If someone managed to do it, please make it public.</b></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">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).</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">In the attachments there are my macros for the phantom, the source, and the test source which is a part of the whole source.</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">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.</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">There are some generations in regions, in which theoretically they should not appear. </div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">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).</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">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.</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">I would be very grateful for your help and for any suggestions that would lead to the solution.</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr"><br></div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">Best regards,</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">PaweÅ‚ Kowalski,</div><div id="yiv4949722788yui_3_16_0_ym19_1_1518533684700_20627" dir="ltr">National Centre for Nuclear Research</div></div></div></div></div></body></html>