<div dir="ltr"><div dir="ltr">Hi Loïc,<div><br></div><div>Thank you very much for actively answering DoseActor related questions in the past ten years. I saw your answers in the GATE archive quite often. I did read through almost all example codes in GATE-RTiON, but I did not find the example where DoseActor is used to record energy deposited by ions for dosimetry. </div><div><br></div><div>Here is the problem I have found with Edep in DoseActor,<a href="http://lists.opengatecollaboration.org/pipermail/gate-users/2009-July/003406.html" target="_blank"> a similar question</a> was also asked many years ago. No matter what source I use, like I-131, Lu-177,

Edep in DoseActor always returns <b>1 (Mev) per disintegration</b>. Dose in DoseActor also utilizes the value returned Edep, which leads to a wrong result. I would rather believe that I made a mistake somewhere than think GATE's DoseActor failed in my case. </div><div><br></div><div>Could you or any other GATE user have a look at my code? It is self-containt code. You will finish the entire verification in less than 5 minutes. Here is what the GATE code does: </div><div>
<div>
<ul type="disc">
<li class="MsoNormal">
Phantom: 101*101*101 voxelized tissue ( water equivalent) phantom, with a side of 2 mm for each voxel <u></u><u></u></li><li class="MsoNormal">
Source: 10 mm radius sphere voxelized Lu177 source with a side of 2 mm for each voxel, and being placed in the center of the phantom<u></u><u></u></li><li class="MsoNormal">
DoseActor: record dose, Edep, and export a mass image of the phantom (all 3D array)<u></u><u></u></li><li class="MsoNormal">
ProductionAndStoppingActor: record the number of Lu177 generated in each voxel, indicated by Lu177-Stop (3D array)
<u></u><u></u></li></ul>
<div>
<p class="MsoNormal">        According to the user manual, "The unit of edep is MeV and the unit of dose is Gy". </p><p class="MsoNormal"><br></p><p class="MsoNormal">You will also find a python script (result-analysis.ipynb) to do the analysis:</p><p class="MsoNormal"></p><ul><li>reads dose, edep, mass, and particle-prod</li><li>calculate the average energy deposited by a single disintegration.</li></ul><div>Thank you very much for any help and verification. </div><div><br></div><div>Albert</div><p></p></div></div></div><div><br></div><div><br></div></div><div><br></div><div><br></div><div><br></div><div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Mar 5, 2021 at 12:48 AM Grevillot Loic <<a href="mailto:loic.grevillot@medaustron.at" target="_blank">loic.grevillot@medaustron.at</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 lang="DE">
<div>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">FYI,<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">DoseActor validation tests were provided for the
<a href="https://gate.uca.fr/download/gate-rtion" target="_blank">GATE-RTiON</a> release, in the context of light ion beam therapy under the
<a href="https://github.com/OpenGATE/GateContrib/tree/master/GATE-RTion" target="_blank">gitcontrib</a> repository.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">The tests can be applied in principle to any GATE release and may help you as a starting point.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">It would be nice to increase the collection of tests for other clinical applications such as  yours!<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="FR" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Best,<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="FR" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Loïc<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif">From:</span></b><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif"> Gate-users [mailto:<a href="mailto:gate-users-bounces@lists.opengatecollaboration.org" target="_blank">gate-users-bounces@lists.opengatecollaboration.org</a>]
<b>On Behalf Of </b>Albert Grace Lieu<br>
<b>Sent:</b> Freitag, 5. März 2021 05:59<br>
<b>To:</b> gate-users <<a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank">gate-users@lists.opengatecollaboration.org</a>><br>
<b>Subject:</b> [Gate-users] Edep returned by DoseActor is about 7 times bigger then hand calculated value<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">Hi Gate users, <u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I am validating the validity of DoseActor in GATE for radionuclide dosimetry application. Surprisingly, I found the Edep returned by DoseActor is about 7 times bigger than my hand calculated value, so does the dose returned by DoseActor.
 I could be wrong, but I double-checked everything I did and had been as prudent as I can before asking for suggestions here. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Here is what I did in GATE:<u></u><u></u></p>
</div>
<div>
<ul type="disc">
<li class="MsoNormal">
Phantom: 101*101*101 voxelized tissue ( water equivalent) phantom, with a side of 2 mm for each voxel <u></u><u></u></li><li class="MsoNormal">
Source: 10 mm radius sphere voxelized Lu177 source with a side of 2 mm for each voxel, and being placed in the center of the phantom<u></u><u></u></li><li class="MsoNormal">
DoseActor: record dose, Edep, and export a mass image of the phantom (all 3D array)<u></u><u></u></li><li class="MsoNormal">
ProductionAndStoppingActor: record the number of Lu177 generated in each voxel, indicated by Lu177-Stop (3D array)
<u></u><u></u></li></ul>
<div>
<p class="MsoNormal">        According to the user manual, "The unit of edep is MeV and the unit of dose is Gy". <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Here is how I calculate the average of the energy of Lu177 in water, Lu177 is a beta emitter, the range of electron in water is way smaller than 2mm, thus it is stopped where is created. Given the branching ratio of Lu177, on average a
 single Lu177 disintegration would deposit 134.16 Kev (electrons only) or 162.95 Kev (electrons and gammas).<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal"><img border="0" width="279" height="181" id="m_7937145216447734549gmail-m_745001018199457084_x0000_i1025" src="cid:17804b8b9eb4cff311" alt="RTEmagicC_Lu-177_decays.png.png"><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">However, when I divide Edep by Lu177-Stop, the energy per Lu177 disintegration returned by DoseActor is
<span style="color:red">~1.0 Mev</span>, which <span style="color:red">7</span> times bigger than my hand calculated value. Is this an internal issue? I attach my Gate code for your review.  <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Thanks a lot for your help. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Sincere wishes, <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Stay safe and healthy<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Albert <u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
<div>
<p class="MsoNormal"><small><span><u></u><u></u></span><b><span style="font-family:Verdana,sans-serif;color:rgb(17,35,89)" lang="EN-GB">Disclaimer:</span></b><br>
<span style="font-family:Verdana,sans-serif;color:rgb(17,35,89)" lang="EN-GB">Please notice our E-Mail Disclaimer <a href="http://www.medaustron.at/email-disclaimer/" target="_blank"><span style="color:rgb(17,35,89)">http://www.medaustron.at/email-disclaimer/</span></a></span></small><span><u></u><u></u></span><u></u> <u></u></p>
</div>
<br>
</div>

</blockquote></div></div>