<div dir="ltr"><div><div><div><div><div>Hi Ottavia,<br><br>Keep in mind that your phantom is the daughter of the world and your cold area is the daughter of your phantom.<br>Thus, you have to attach phantomSD to both cold area and cylinder phantom .<br> <br>Here are some codes tested successfully and used by my previous project, hope it will help you.<br><br>In this code, there are two cold areas inside a water phantom.<br><br>#=====================================================<br> # P H A N T O M<br>#=====================================================<br><br>/gate/world/daughters/name phantom<br>/gate/world/daughters/insert cylinder<br>/gate/phantom/setMaterial Water<br>/gate/phantom/vis/setColor grey<br> /gate/phantom/geometry/setRmin 0. mm<br>/gate/phantom/geometry/setRmax 6. mm<br>/gate/phantom/geometry/setHeight 0.3 mm<br>/gate/phantom/attachPhantomSD<br><br>/gate/world/daughters/name cold_area <br>/gate/world/daughters/insert cylinder <br> /gate/cold_area/vis/forceWireframe <br>/gate/cold_area/vis/setColor green <br>/gate/cold_area/geometry/setRmax 1.0 mm<br>/gate/cold_area/geometry/setRmin 0.0 mm<br>/gate/cold_area/geometry/setHeight 0.3 mm<br>/gate/cold_area/placement/setTranslation -1.0 0.0 0.0 mm<br> /gate/cold_area/setMaterial Water<br>/gate/cold_area/attachPhantomSD<br><br>/gate/world/daughters/name cold_area_a <br>/gate/world/daughters/insert cylinder <br>/gate/cold_area_a/vis/forceWireframe <br>/gate/cold_area_a/vis/setColor green <br> /gate/cold_area_a/geometry/setRmax 0.5 mm<br>/gate/cold_area_a/geometry/setRmin 0.0 mm<br>/gate/cold_area_a/geometry/setHeight 0.3 mm<br>/gate/cold_area_a/placement/setTranslation 1.4 0.0 0.0 mm<br>/gate/cold_area_a/setMaterial Water<br> /gate/cold_area_a/attachPhantomSD<br><br><br>#=====================================================<br># P A R T I C L E S O U R C E<br>#=====================================================<br>/gate/source/addSource number1 <br> /gate/source/number1/setActivity 37000000. becquerel <br>/gate/source/number1/gps/particle e+ <br>/gate/source/number1/gps/energytype Fluor18 <br>/gate/source/number1/setForcedUnstableFlag true<br> /gate/source/number1/setForcedHalfLife 6586 s<br>/gate/source/number1/gps/type Volume <br>/gate/source/number1/gps/shape Cylinder <br>/gate/source/number1/gps/radius 2. mm <br>/gate/source/number1/gps/halfz 0.15 mm <br> /gate/source/number1/gps/centre 0. 0. 0. mm <br>/gate/source/number1/visualize 2000 red 3 <br>/gate/source/number1/gps/angtype iso <br>/gate/source/number1/gps/Forbid cold_area <br>/gate/source/number1/gps/Forbid cold_area_a<br> /gate/source/number1/dump 2 <br><br>/gate/source/list<br><br><br><br><br></div>Also, VERY important, you need to run /gate/run/initialize before the PARTICLE SOURCE part.<br><br><br></div>Cheers,<br>Kaiyang Li<br><br></div> Centre for Medical Radiation Physics<br></div>University of Wollongong<br></div>Australia<br><div><div><div><div><div><br><br><br></div></div></div></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote"> 2013/9/11 Ottavia Bertolli <span dir="ltr"><<a href="mailto:ottavia.bertolli@gmail.com" target="_blank">ottavia.bertolli@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"> <div dir="ltr"><div><div><div>Dear Gate users,<br><br></div>I would like to define a Volume source that has the shape of an annular cylinder.<br></div>From the Users Guide it seems impossible to do this just defining the shape, since "annulus" is a Plane source, thus I thought about exploiting the Forbid command, to define a cold (i.e. with no activity) volume in a phantom, that in my case should be a smaller cylinder (filled with water) to be located within the bigger cylinder source.<br> <br></div>This is my code:<br><div><br>===============================================================<br><br></div><div># COLD AREA<br></div><div><br>/gate/world/daughters/name cold_area<br>/gate/world/daughters/insert cylinder<br> /gate/cold_area/placement/setTranslation 0.0 0.0 0.0 mm<br>/gate/cold_area/geometry/setRmax 0.5 cm<br>/gate/cold_area/geometry/setRmin 0.0 cm<br>/gate/cold_area/geometry/setHeight 4.0 cm<br> /gate/cold_area/setMaterial Water<br>#/gate/cold_area/vis/forceWireframe<br>/gate/cold_area/vis/forceSolid<br>/gate/cold_area/vis/setColor green<br><br><br># S O U R C E<br> /gate/source/addSource oxy15<br> /gate/source/oxy15/setActivity 3700000 becquerel <br>/gate/source/oxy15/gps/particle ion<br>/gate/source/oxy15/gps/ion 8 15 0 0 <br>/gate/source/oxy15/gps/energytype Mono<br>/gate/source/oxy15/gps/monoenergy 0. MeV<br>/gate/source/oxy15/gps/angtype iso<br> /gate/source/oxy15/gps/number 1<br>/gate/source/oxy15/gps/centre 0. 0. 0. cm<br>/gate/source/oxy15/gps/type Volume<br>/gate/source/oxy15/gps/shape Cylinder<br>/gate/source/oxy15/gps/radius 1 cm<br>/gate/source/oxy15/gps/halfz 2 cm<br> /gate/source/oxy15/gps/Forbid cold_area<br>/gate/source/oxy15/dump 1 <br><br>/gate/source/list<br><br>==================================================================<br><br></div><div>When I enable the "/gate/source/oxy15/gps/Forbid cold_area" line command, I get this error:<br> <br>*** Break *** segmentation violation<br>Attaching to program: /proc/2223/exe, process 2223<br>[Thread debugging using libthread_db enabled]<br>0x0054d422 in __kernel_vsyscall ()<br>error detected on stdin<br>A debugging session is active.<br> <br> Inferior 1 [process 2223] will be detached.<br><br></div><div>What is it that I do wrong? Does anybody know the right way to define a cold region inside a bigger source?<br><br>Thanks in advance!<span class="HOEnZb"><font color="#888888"><br> <br></font></span></div><span class="HOEnZb"><font color="#888888"><div>Ottavia Bertolli<br> </div></font></span></div> <br>_______________________________________________<br> Gate-users mailing list<br> <a href="mailto:Gate-users@lists.opengatecollaboration.org">Gate-users@lists.opengatecollaboration.org</a><br> <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br></blockquote></div><br></div>