Hello Kaiyang Li,<br><br>I'm still struggling with "forbid".<br><br>I tried the simple geometry, too, but my simulation stops with "***Break*** segmentation violation" and I cannot proceed to even see the reconstructed image... this didn't happen before...<br>
<br>Anyway, it seems the command is working somehow. <br>When I tried to forbid the whole source by surrounding with the "forbidding" volume and start simulation, GATE gave me the message of "primaries were always generated in forbidden volumes!" and the simulation was shut down.<br>
<br>Also, you can visualize the source by<br>/gate/source/FictiveSource/visualize "count" "color" "size"<br>This command should reflect the "forbid" as well, so we are supposed to be able to see the "forbidden" source... but...<br>
<br>You can check how the "forbid" command works in GATE in GateSPSPosDistribution.cc line:150<br><br>After all, I cannot give you any solution to this problem.<br>I'll still work on with this, so if you find anything new, please let me know.<br>
<br>Taku<br><br><div class="gmail_quote">2012/12/10 Kaiyang Li <span dir="ltr"><<a href="mailto:kl780@uowmail.edu.au" target="_blank">kl780@uowmail.edu.au</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">Hi Taku,<br>
<br>
My name is Kaiyang Li, a GATE user.<br>
I have the problem of using "forbid" in my simulation. It seems that the forbid area is simply<br>
not forbidden at all.<br>
<br>
I have read all other questions, including yours. Have you figure it out yet?<br>
I have used GATE for half a year and run some other complex simulations successfully but for this forbid one, it seems something is wrong.<br>
<br>
Here is my simple code followed by the GATE instruction:<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 2.0 mm <br>
/gate/cold_area/geometry/setHeight 0.3 mm<br>
/gate/cold_area/setMaterial Water<br>
<br>
/gate/run/initialize<br>
<br>
/gate/source/addSource number1 <br>
/gate/source/number1/setActivity 37000000. becquerel <br>
/gate/source/number1/gps/particle gamma <br>
/gate/source/number1/setType backtoback <br>
/gate/source/number1/gps/monoenergy 511. keV <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 4. mm <br>
/gate/source/number1/gps/halfz 0.15 mm <br>
/gate/source/number1/gps/centre 0. 0. 0. mm <br>
/gate/source/number1/gps/angtype iso <br>
/gate/source/number1/gps/Forbid cold_area <br>
/gate/source/number1/dump 1 <br>
/gate/source/list<br>
<br>
The reconstructed image shows that there is no forbidden area.<br>
Probably the forbidden command is not working, what do you think?<br>
<br>
<br>
Kind Regards,<br>
Kaiyang Li <br>
</div>
</div>
</blockquote></div><br><br clear="all"><br>-- <br>Hokkaido University, Graduate school of medicine, medical physics and engineering D1<br>北海道大学 大学院 医学研瘢雹究科 医学物理工学分野 D1<br>Taku Yamanaka<br>山中 琢<br><br><a href="mailto:takuyama@med.hokudai.ac.jp" target="_blank">takuyama@med.hokudai.ac.jp</a><br>
--<br>