<div dir="ltr">Thanks, Bryan for the clarification. <div><br></div><div>I definitely misunderstood the command "confine" earlier. Now it seems clear to me, but it comes with a couple of concerns. It looks like you created a very large Lu176 source to encompass/embody your crystals. But by confining the source to crystal_P, does it mean that the source is only activated within crystal_P and outside crystal_P the source is not activated?  If so, how can the source strength inside crystal_P be correctly calculated and implemented? </div><div><br></div><div>My second question is, why does "confine" recognize crystal_P as the assembly of all crystals created with repeaters, but "attachTo" command does not recognize crystal_P at all?  Meanwhile, crystal only represents the very first one crystal created with a repeater. Is this an internal bug?  If "attachTo" recognized crystal_P, then none of these would have become a problem. I really want to bring this up to the attention of the developers. </div><div><br></div><div>Thanks a lot for your help.</div><div>Albert </div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jun 12, 2020 at 5:31 AM Bryan McIntosh <<a href="mailto:mcintoshster@gmail.com">mcintoshster@gmail.com</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"><div dir="ltr"><div>Hi Albert,</div><div><br></div><div>Ashok's code looks very similar to what I had sent, and it only uses the confine command, NOT attach. The whole idea behind the way that I implemented this method for simulating a Lu-176 ion source was to make the source large enough to fill the entire scanner volume, then confine the activity to the crystal only; again, I have not tested it under a newer version of GATE to verify that it still works, and it may not.</div><div><br></div><div>Again, be advised that this does increase the simulation time by a significant amount since you're adding an ion source that has a lot of beta and gamma interactions inside the crystal volume; if you can simulate the noise due to background Lu-176 using the "noise" modules as Emilie and Maxime suggested earlier, I highly recommend doing that instead. I know that computers have significantly faster single-thread performance than in 2010 and 2011, but the fewer particles that need to be simulated the better!</div><div><br></div><div>Best regards,</div><div><br></div><div>-Bryan McIntosh<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jun 11, 2020 at 6:20 PM Albert Grace Lieu <<a href="mailto:albertnew2018@gmail.com" target="_blank">albertnew2018@gmail.com</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">Hi Bryan, <div><br></div><div>Thank you very much for your paper and help. I tried your code for a couple of times, it does not confine/attach Lu176 to the crystal. Instead,  there is only one source showed up at the center of the world. Also, your code does not look like confining a source to crystals, because of the shape and size of the source. It is much larger than any possible detector crystal. My problem now is that I can only attach Lu176 to only one crystal, not all of them. Could you please give me any further help, please? </div><div><br></div><div>Thanks a lot. </div><div>Zhengzhi </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jun 11, 2020 at 3:02 PM Bryan McIntosh <<a href="mailto:mcintoshster@gmail.com" target="_blank">mcintoshster@gmail.com</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"><div>Hello Zhengzhi,</div><div><br></div><div>I wrote the second paper that you linked, and it looks like your initial idea was in the right direction; I made a Lu-176 source that encompassed the entire detector. The difference is that I used the "confine" command for the source, not "attach;" the complete code was as follows:</div><div><br></div><div>#Single intrinsic source; confined to crystal<br>/gate/source/addSource intrinsic<br>/gate/source/intrinsic/gps/particle ion<br>/gate/source/intrinsic/gps/ion 71 176 0 0<br>/gate/source/intrinsic/gps/monoenergy 0. keV #Super important!<br>/gate/source/intrinsic/gps/angtype iso<br>/gate/source/intrinsic/gps/type Volume<br>/gate/source/intrinsic/gps/shape Cylinder<br>/gate/source/intrinsic/gps/radius 14.5 cm<br>/gate/source/intrinsic/gps/halfz 70.1 mm<br>/gate/source/intrinsic/setActivity 20192.56 Bq<br></div><div>/gate/source/intrinsic/gps/confine crystal_P<br></div><div><br></div><div>This worked for a system where the crystal volume was named called "crystal," and was duplicated using a cubic array. Be advised that the activity and dimensions of the source were set up to match the Siemens Inveon dPET system; adjust the size and activity accordingly for your particular system. I used an activity of 276 Bq/cm^3 of LSO crystal, calculated from a 2.6% concentration of Lu-176 in the Lu atoms making up the crystal. This is different than the 250 Bq/cm^3 that is often used, since that number was from a measurement that was 250 counts per second/cm^3, not necessarily Bq.</div><div><br></div><div>This code worked when I ran it back in 2011; I haven't tested it on a GATE version newer than 6.2. I THINK it will work with newer versions of GATE, but your results may vary!</div><div><br></div><div>Good luck!</div><div><br></div><div>-Bryan McIntosh</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jun 11, 2020 at 1:37 PM <<a href="mailto:gate-users-request@lists.opengatecollaboration.org" target="_blank">gate-users-request@lists.opengatecollaboration.org</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">------------------------------<br>
<br>
Message: 2<br>
Date: Wed, 10 Jun 2020 22:09:13 -0700<br>
From: Albert Grace Lieu <<a href="mailto:albertnew2018@gmail.com" target="_blank">albertnew2018@gmail.com</a>><br>
To: gate-users <<a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank">gate-users@lists.opengatecollaboration.org</a>><br>
Cc: Peter Olcott <<a href="mailto:polcott@reflexion.com" target="_blank">polcott@reflexion.com</a>><br>
Subject: [Gate-users] How to simulate Lu-176 in LSO crystal<br>
Message-ID:<br>
        <CAGLxnT-Qdh-NZ8KoBOQuH=<a href="mailto:zszHP%2Bzw5f_aeD5ADLNDP9zYzG6Q@mail.gmail.com" target="_blank">zszHP+zw5f_aeD5ADLNDP9zYzG6Q@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear Gate users and developers,<br>
<br>
There are a couple of papers [1<br>
<<a href="https://aapm.onlinelibrary.wiley.com/doi/abs/10.1118/1.4824694" rel="noreferrer" target="_blank">https://aapm.onlinelibrary.wiley.com/doi/abs/10.1118/1.4824694</a>>,2<br>
<<a href="https://ieeexplore.ieee.org/document/5751205" rel="noreferrer" target="_blank">https://ieeexplore.ieee.org/document/5751205</a>>] addressing the importance<br>
of adding the simulation of Lu-176 in Lutetium based crystal in GATE<br>
simulation. I try to do this to my simulated PET system that has thousands<br>
of pieces of crystals (created with linear and cubic array repeater). My<br>
general idea of implementing this is to create a Lu176 ion source and then<br>
attach it to my crystals. However,  when I do:<br>
<br>
/gate/source/Lu176/attachTo Crystal<br>
<br>
the Lu176 is only attached to a single piece of crystal (which can be<br>
visualized with /gate/source/Lu176/visualize 1000 red 3).  Has anyone<br>
implemented this before? Could you please share your thoughts?<br>
<br>
Thanks a lot. Greatly appreciated.<br>
<br>
Zhengzhi<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200610/8e812275/attachment-0001.html" rel="noreferrer" target="_blank">http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200610/8e812275/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
Gate-users mailing list<br>
<a href="mailto:Gate-users@lists.opengatecollaboration.org" target="_blank">Gate-users@lists.opengatecollaboration.org</a><br>
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" rel="noreferrer" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br>
<br>
------------------------------<br>
<br>
End of Gate-users Digest, Vol 169, Issue 16<br>
*******************************************<br>
</blockquote></div></div>
</blockquote></div>
</blockquote></div></div>
</blockquote></div>