<div dir="ltr">Maybe this figure is better to illustrate my problem. I expected to get an annulus source after simultaneously moving phantom and colde_area to a new location, and attaching the source to the phantom, however, it looks after the movement, the source didn't dump the portion in colde_area. <div><img src="cid:ii_k8hy06ke0" alt="movement_of_annulus_source1.png" width="354" height="337"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Apr 1, 2020 at 3:23 PM Zhengzhi Liu <<a href="mailto:zliu36@stanford.edu">zliu36@stanford.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear Gate users,<br><div><br></div><div>I am trying to move NEMA IQ phantom with the source along the axial direction back and forth to simulate a real experiment that we have done. I started with a much simpler case where there is only an annulus phantom with an annulus source in the simulated PET system. Phantom as a volume is very easy to be translated. According to the user manual, to move a source, the source had better be attached to a volume. For common source without cold_area, it works very well. However, I found it did not work in the case where the source is determined by a volume and a cold_area (which is also a volume), such as the annulus source. Could anyone teach me how to do so if it is possible?</div><div><br></div><div>Following is my code, I moved the phantom and cold_area simultaneously, whether or not I attached the source to phantom, after the translation, GATE did not return an annulus source rather a cylindrical source same as the phantom as shown in the attached screenshot. <br></div><div><br></div><div>phantom.mac<br></div><div>#=====================================================<br># WATER PHANTOM<br>#=====================================================<br>/gate/world/daughters/name Normal_AP<br>/gate/world/daughters/insert cylinder<br>/gate/Normal_AP/geometry/setRmax 30.051 cm<br>/gate/Normal_AP/geometry/setRmin 30.0 cm<br>/gate/Normal_AP/geometry/setHeight 6 cm<br>/gate/Normal_AP/placement/setTranslation 0.0 0.0 0.0 cm<br>/gate/Normal_AP/moves/insert genericMove<br>/gate/Normal_AP/genericMove/setPlacementsFilename ./macros/phantomMovement0.placements<br>/gate/Normal_AP/setMaterial Water<br>/gate/Normal_AP/vis/forceSolid<br>/gate/Normal_AP/vis/setColor white<br>/gate/Normal_AP/vis/setVisible 1<br><br>/gate/world/daughters/name cold_area<br>/gate/world/daughters/insert cylinder<br>/gate/cold_area/geometry/setRmax 30.00 cm<br>/gate/cold_area/geometry/setHeight 6 cm<br>/gate/cold_area/placement/setTranslation 0.0 0.0 0.0 cm<br>/gate/cold_area/moves/insert genericMove<br>/gate/cold_area/genericMove/setPlacementsFilename ./macros/phantomMovement0.placements<br>/gate/cold_area/setMaterial Air<br>/gate/cold_area/vis/forceWireframe<br>/gate/cold_area/vis/setColor white</div><div><br></div><div><br></div><div>source.mac<br></div><div>#=====================================================<br># FASTER PRODUCTION MODULE FOR GE-68 ISOTOPES<br>#=====================================================<br>/gate/source/addSource Ge68<br>/gate/source/Ge68/setActivity 37000 Bq #1. uci<br>/gate/source/Ge68/gps/particle gamma<br>/gate/source/Ge68/setType backtoback<br>/gate/source/Ge68/gps/type Volume<br>/gate/source/Ge68/gps/shape Cylinder<br>/gate/source/Ge68/gps/radius 30.051 cm<br>/gate/source/Ge68/gps/halfz 3 cm<br>/gate/source/Ge68/gps/centre 0. 0. 0. cm<br>#/gate/source/Ge68/attachTo Normal_AP<br>/gate/source/Ge68/gps/energytype Mono<br>/gate/source/Ge68/gps/monoenergy 0.511 MeV<br>/gate/source/Ge68/gps/angtype iso<br>/gate/source/Ge68/setForcedUnstableFlag false<br> #/gate/source/Ge68/setForcedHalfLife 6586 s<br>/gate/source/Ge68/gps/Forbid cold_area<br>/gate/source/Ge68/dump 1<br><br>/gate/source/list</div><div><br></div><div><br></div><div>Thanks a lot for your reading and any help.</div><div>Sincerely,</div><div>Zhengzhi </div><div><br></div></div>
</blockquote></div>