<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"> <head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> <meta name="Generator" content="Microsoft Word 14 (filtered medium)"> <style></style> </head> <body lang="EN-US" link="blue" vlink="purple"> <div class="WordSection1"> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Hello,<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Currently the DoseToWater part is incomplete and should contain other particles, e.g. carbon ions, helium, etc.<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">The else part only means at the moment that if the system encounter another type of particle, it will calculate the dose to water assuming this particle is equivalent to a proton with an energy of 100 MeV.<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">If we would remove the if (which could also be an option), then the dose to water would be equal to the dose, as it would not be rescaled by the stopping power ratio.<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">So currently the DoseToWater works for electrons, photons and proton beams as far as I know.<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">I have in my to do list to upgrade it for particles up to carbon ions, but did not have time yet, that being said it should not be very difficult, just require a bit of time and also to find out the correct naming convention used in Geant4 for the different particles that can be generated.<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Hope this help!<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="DE" style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1F497D">Best,<o:p></o:p></span></p> <p class="MsoNormal"><span lang="DE" style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1F497D">Loïc</span><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p> <div> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in"> <p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> Manuel Bardiès [mailto:manuel.bardies@inserm.fr] <br> <b>Sent:</b> Thursday, February 02, 2017 5:27 PM<br> <b>To:</b> Luc SIMON<br> <b>Cc:</b> Grevillot Loic; gate-users@lists.opengatecollaboration.org<br> <b>Subject:</b> Re: [Gate-users] Gate dose to water problem<o:p></o:p></span></p> </div> </div> <p class="MsoNormal"><o:p> </o:p></p> <p class="MsoNormal">And thanks to Andreas as well (Luc missed his mail ;-)<o:p></o:p></p> <div> <p class="MsoNormal"><o:p> </o:p></p> <div> <blockquote style="margin-top:5.0pt;margin-bottom:5.0pt"> <div> <p class="MsoNormal">Le 2 févr. 2017 à 17:23, Luc SIMON <<a href="mailto:uhqd75@gmail.com">uhqd75@gmail.com</a>> a écrit :<o:p></o:p></p> </div> <p class="MsoNormal"><o:p> </o:p></p> <div> <div> <div> <p class="MsoNormal">Dear all and in particular:<br> <b><a href="mailto:david.sarrut@creatis.insa-lyon.fr">david.sarrut@creatis.insa-lyon.fr</a><br> <a href="mailto:gergely.patay@gmail.com">gergely.patay@gmail.com</a><br> <a href="mailto:alexvergaragil@gmail.com">alexvergaragil@gmail.com</a><br> <a href="mailto:thomas@deschler.fr">thomas@deschler.fr</a></b><o:p></o:p></p> </div> <p class="MsoNormal"><b>and all the authors of <a href="http://doseactor.cc">DoseActor.cc</a> ...</b><o:p></o:p></p> <div> <p class="MsoNormal" style="margin-bottom:12.0pt"><o:p> </o:p></p> </div> <div> <p class="MsoNormal"><b>Following the mail of Tony Younes we investigated this bug.</b><o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-bottom:12.0pt"><b>We found that the NaN values seem to come from <a href="http://dosetoactor.cc">DoseToActor.cc</a> at this level (l 360)<br> </b><br> <br> <br> <br> // Other particles should be taken into account (Helium etc), but bug ? FIXME<br> if (PartName== "proton" || PartName== "e-" || PartName== "e+" || PartName== "deuteron") {<br> <br> DEDX = emcalc->ComputeTotalDEDX(Energy, PartName, material, cut);<br> DEDX_Water = emcalc->ComputeTotalDEDX(Energy, PartName, "G4_WATER", cut);<br> doseToWater=edep/density/Volume/gray*(DEDX_Water/1.)/(DEDX/(density*e_SI));<br> }<br> else {<br> DEDX = emcalc->ComputeTotalDEDX(100, "proton", material, cut);<br> DEDX_Water = emcalc->ComputeTotalDEDX(100, "proton", "G4_WATER", cut);<br> doseToWater=edep/density/Volume/gray*(DEDX_Water/1.)/(DEDX/(density*e_SI));<br> }<br> <br> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-bottom:12.0pt"><b>In this "if" section sometimes DEDX <u>and </u>DEDX_Water are 0.0 and the division 0 / 0 is NaN (please note that 3 / 0 is infinite, which is not a NaN value...)</b><o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-bottom:12.0pt">This problem is fixed. But now another question : why in the "else" part DEDX_water is computed for proton ????? (while proton should be considered in the first part... )<br> <br> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-bottom:12.0pt">Moreover, photons should not contribute to "dose". So why do we need this "else" part ????<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-bottom:12.0pt">Thank you for your help.<o:p></o:p></p> </div> <div> <p class="MsoNormal">Luc & Tony<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="margin-bottom:12.0pt"><o:p> </o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal" style="margin-bottom:12.0pt"><br> <br> <o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> <div> <p class="MsoNormal">2017-01-31 16:22 GMT+01:00 Grevillot Loic <<a href="mailto:loic.grevillot@medaustron.at" target="_blank">loic.grevillot@medaustron.at</a>>:<o:p></o:p></p> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Hi,</span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Could you send macros and output files?</span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="DE" style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1F497D">Best,</span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="DE" style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1F497D">Loïc</span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"> </span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"> </span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"> </span><o:p></o:p></p> <div> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in"> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> tony youness [mailto:<a href="mailto:tony_youness@live.com" target="_blank">tony_youness@live.com</a>] <br> <b>Sent:</b> Tuesday, January 31, 2017 3:14 PM<br> <b>To:</b> Grevillot Loic; <a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank"> gate-users@lists.opengatecollaboration.org</a><br> <b>Subject:</b> Re: Gate dose to water problem</span><o:p></o:p></p> </div> </div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Hi,<o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I am simulating high energy photon beams (6 and 10 MeV) for external radiotherapy applications. <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I tried both ("Water" and "G4_WATER") and I am still getting NAN values (only when the absorbed dose is reported as absorbed dose to water). The voxel values between both absorbed dose report (Dose to medium and Dose to water) are also different even though that the medium is WATER.<o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Regards, <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Tony<o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">On 31/01/17 8:37 am, Grevillot Loic wrote:<o:p></o:p></p> </div> <blockquote style="margin-top:5.0pt;margin-bottom:5.0pt"> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Hello,</span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"> </span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Which beam are you simulating? </span><span lang="FR" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">(photons, electrons, light ions, etc.)</span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Which water did you define in your simulation? “Water” is different from “GEANT4_WATER” -> could you try using GEANT4_WATER?</span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"> </span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="DE" style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1F497D">Best,</span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="DE" style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1F497D">Loïc</span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"> </span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"> </span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"> </span><o:p></o:p></p> <div> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in"> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> Gate-users [<a href="mailto:gate-users-bounces@lists.opengatecollaboration.org" target="_blank">mailto:gate-users-bounces@lists.opengatecollaboration.org</a>] <b>On Behalf Of </b>tony youness<br> <b>Sent:</b> Monday, January 30, 2017 1:04 PM<br> <b>To:</b> <a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank"> gate-users@lists.opengatecollaboration.org</a><br> <b>Subject:</b> [Gate-users] Gate dose to water problem</span><o:p></o:p></p> </div> </div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Dear All,<o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I am facing a problem with gate dose actor. When enabling dose to water i get Nan values in some of the voxels in the water region of my phantom. I also get different values in the water voxels between both dose reports (dose to medium and dose to water). It seems like there's a bug when gate tries to converts dose to medium into the dose to water <u>when in the water region</u>. Has anyone already faced this problem? Any help would be much appreciated.<o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Kind Regards,<o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Tony Younes<o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b><span lang="EN-GB" style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#112359">Disclaimer:</span></b><span style="font-size:10.0pt"><br> </span><span lang="EN-GB" style="font-size:10.0pt;font-family:"Verdana","sans-serif"">Please notice our E-Mail Disclaimer <a href="http://www.medaustron.at/email-disclaimer/" target="_blank"><span style="color:#112359">http://www.medaustron.at/email-disclaimer/</span></a></span> <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </blockquote> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b><span lang="EN-GB" style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#112359">Disclaimer:</span></b><span style="font-size:10.0pt"><br> </span><span lang="EN-GB" style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#112359">Please notice our E-Mail Disclaimer <a href="http://www.medaustron.at/email-disclaimer/" target="_blank"><span style="color:#112359">http://www.medaustron.at/email-disclaimer/</span></a></span> <o:p></o:p></p> </div> <p class="MsoNormal"><o:p> </o:p></p> </div> <p class="MsoNormal"><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><o:p></o:p></p> </div> <p class="MsoNormal"><br> <br clear="all"> <br> -- <o:p></o:p></p> <div> <div> <div> <div> <div> <p class="MsoNormal">___________________<br> <br> SIMON Luc, Physicien Médical, PhD<br> Département d'Ingénierie et de Physique Médicale - IUCT-Oncopole - Toulouse, FRANCE<o:p></o:p></p> </div> <div> <p class="MsoNormal">INSERM UMR 1037 (CRCT) - Equipe 15 - Toulouse, FRANCE<o:p></o:p></p> </div> </div> </div> </div> </div> </div> <p class="MsoNormal">_______________________________________________<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">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><o:p></o:p></p> </div> </blockquote> </div> <p class="MsoNormal"><o:p> </o:p></p> </div> </div> <div class="WordSection1"> <p class="MsoNormal"><small><span style=""><o:p></o:p></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/"><span style="color: rgb(17, 35, 89);">http://www.medaustron.at/email-disclaimer/</span></a></span></small><span style=""><o:p></o:p></span><o:p> </o:p></p> </div> <br> </body> </html>