<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">And thanks to Andreas as well (Luc missed his mail ;-)<div class=""> <div><blockquote type="cite" class=""><div class="">Le 2 févr. 2017 à 17:23, Luc SIMON <<a href="mailto:uhqd75@gmail.com" class="">uhqd75@gmail.com</a>> a écrit :</div> <div class=""><div dir="ltr" class=""><div class="">Dear all and in particular: <a href="mailto:david.sarrut@creatis.insa-lyon.fr" class="">david.sarrut@creatis.insa-lyon.fr</a> <a href="mailto:gergely.patay@gmail.com" class="">gergely.patay@gmail.com</a> <a href="mailto:alexvergaragil@gmail.com" class="">alexvergaragil@gmail.com</a> <a href="mailto:thomas@deschler.fr" class="">thomas@deschler.fr</a> </div>and all the authors of <a href="http://doseactor.cc" class="">DoseActor.cc</a> ... <div class=""> </div><div class="">Following the mail of Tony Younes we investigated this bug. </div><div class="">We found that the NaN values seem to come from <a href="http://dosetoactor.cc" class="">DoseToActor.cc</a> at this level (l 360) // Other particles should be taken into account (Helium etc), but bug ? FIXME if (PartName== "proton" || PartName== "e-" || PartName== "e+" || PartName== "deuteron") { DEDX = emcalc->ComputeTotalDEDX(Energy, PartName, material, cut); DEDX_Water = emcalc->ComputeTotalDEDX(Energy, PartName, "G4_WATER", cut); doseToWater=edep/density/Volume/gray*(DEDX_Water/1.)/(DEDX/(density*e_SI)); } else { DEDX = emcalc->ComputeTotalDEDX(100, "proton", material, cut); DEDX_Water = emcalc->ComputeTotalDEDX(100, "proton", "G4_WATER", cut); doseToWater=edep/density/Volume/gray*(DEDX_Water/1.)/(DEDX/(density*e_SI)); } </div><div class="">In this "if" section sometimes DEDX and 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...) </div><div class="">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... ) </div><div class="">Moreover, photons should not contribute to "dose". So why do we need this "else" part ???? </div><div class="">Thank you for your help. </div><div class="">Luc & Tony </div><div class=""> </div><div class=""> </div><div class=""> </div><div class=""> </div><div class=""> </div></div><div class="gmail_extra"> <div class="gmail_quote">2017-01-31 16:22 GMT+01:00 Grevillot Loic <<a href="mailto:loic.grevillot@medaustron.at" target="_blank" class="">loic.grevillot@medaustron.at</a>>: <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"> <div bgcolor="white" link="blue" vlink="purple" lang="EN-US" class=""> <div class="m_-6703093031588237583WordSection1">Hi,Could you send macros and output files?Best,Loïc <div class=""> <div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0in 0in 0in" class="">From: tony youness [mailto:<a href="mailto:tony_youness@live.com" target="_blank" class="">tony_youness@live.com</a>] Sent: Tuesday, January 31, 2017 3:14 PM To: Grevillot Loic; <a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank" class="">gate-users@lists.opengatecollaboration.org</a> Subject: Re: Gate dose to water problem </div> </div> Hi,I am simulating high energy photon beams (6 and 10 MeV) for external radiotherapy applications. 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.Regards, Tony <div class="">On 31/01/17 8:37 am, Grevillot Loic wrote: </div> <blockquote style="margin-top:5.0pt;margin-bottom:5.0pt" class="">Hello, Which beam are you simulating? (photons, electrons, light ions, etc.)Which water did you define in your simulation? “Water” is different from “GEANT4_WATER” -> could you try using GEANT4_WATER? Best,Loïc <div class=""> <div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0in 0in 0in" class="">From: Gate-users [<a href="mailto:gate-users-bounces@lists.opengatecollaboration.org" target="_blank" class="">mailto:gate-users-bounces@lists.opengatecollaboration.org</a>] On Behalf Of tony youness Sent: Monday, January 30, 2017 1:04 PM To: <a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank" class="">gate-users@lists.opengatecollaboration.org</a> Subject: [Gate-users] Gate dose to water problem </div> </div> Dear All,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 when in the water region. Has anyone already faced this problem? Any help would be much appreciated.Kind Regards,Tony Younes Disclaimer: Please notice our E-Mail Disclaimer <a href="http://www.medaustron.at/email-disclaimer/" target="_blank" class="">http://www.medaustron.at/email-disclaimer/</a> </blockquote> </div> <div class="m_-6703093031588237583WordSection1">Disclaimer: Please notice our E-Mail Disclaimer <a href="http://www.medaustron.at/email-disclaimer/" target="_blank" class="">http://www.medaustron.at/email-disclaimer/</a> </div> </div> _______________________________________________ Gate-users mailing list <a href="mailto:Gate-users@lists.opengatecollaboration.org" class="">Gate-users@lists.opengatecollaboration.org</a> <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" rel="noreferrer" target="_blank" class="">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a> </blockquote></div> -- <div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class="">___________________ SIMON Luc, Physicien Médical, PhD Département d'Ingénierie et de Physique Médicale - IUCT-Oncopole - Toulouse, FRANCE</div><div class="">INSERM UMR 1037 (CRCT) - Equipe 15 - Toulouse, FRANCE</div></div></div></div></div> </div> _______________________________________________ Gate-users mailing list <a href="mailto:Gate-users@lists.opengatecollaboration.org" class="">Gate-users@lists.opengatecollaboration.org</a> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</div></blockquote></div> </div></body></html>