<html><head></head><body><div class="ydp5c8b0e00yahoo-style-wrap" style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:16px;"><div></div>
<div dir="ltr" data-setdir="false">Hi Maria, </div><div dir="ltr" data-setdir="false"><br></div><div dir="ltr" data-setdir="false">I would check the source spectrum. I think the endpoint energy of Y90 beta is 2.28 MeV and is not listed in your source spectrum. Lanconelli et al used the decay data from the Brookhaven National Laboratory database (www.doseinfo-radar.com, RADAR), so check whether you use the same source spectrum for a fair comparison. Have you only used a beta spectrum in simulation, not photon emission x-rays, gammas, and Auger electrons? And a similar value of range or energy cuts?</div><div><br></div><div dir="ltr" data-setdir="false">Regards, </div><div dir="ltr" data-setdir="false">Ashok</div>
</div><div id="yahoo_quoted_3516025315" class="yahoo_quoted">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div>
On Wednesday, 10 February 2021, 10:30:50 am GMT-6, Maria Inês Ribeiro <mid.ribeiro@campus.fct.unl.pt> wrote:
</div>
<div><br></div>
<div><br></div>
<div><div id="yiv4927210009"><div><div dir="ltr"><div>Hello <span class="yiv4927210009gmail-gD">陈戈飞 and Ashok, <br clear="none"></span></div><div><span class="yiv4927210009gmail-gD"><br clear="none"></span></div><div><span class="yiv4927210009gmail-gD">I tried my best to implement your suggestions (thank you sooo much). I found reasonable s-values comparing to the one indicates at <a shape="rect" href="http://medphys.it" target="_blank" rel="noreferrer noopener">medphys.it</a></span></div><div><span class="yiv4927210009gmail-gD">For example, my center voxel s-value is of 4.5 mGy/(MBq s) and theirs is 3 mGy/(MBq s). <br clear="none"></span></div><div><span class="yiv4927210009gmail-gD">However, 2 things seem a bit strange:<br clear="none"></span></div><div><span class="yiv4927210009gmail-gD">1) when I load the dose map at MATLAB, I get a lot of zeros in the border voxels of my kernel. I am using a 11x11x11 kernel, each with 2,21 mm. This is not supposed to be happening.<br clear="none"></span></div><div><span class="yiv4927210009gmail-gD">2) The values are a bit different from the reference values. Do you know if the source should have any additional spectrums that I should worry about? anyway if I add more spectrums, the dose will only increase, which is not good.... <br clear="none"></span></div><div><span class="yiv4927210009gmail-gD">A simulates my source as follows: <br clear="none"></span></div><div><span class="yiv4927210009gmail-gD"><br clear="none"></span></div><div><span class="yiv4927210009gmail-gD">/gate/source/addSource Y90Source<br clear="none">/gate/source/Y90Source/setActivity 1 MBq<br clear="none">/gate/source/Y90Source/gps/particle e-<br clear="none">/gate/source/Y90Source/gps/energytype Arb<br clear="none">/gate/source/Y90Source/gps/hist/type arb<br clear="none">/gate/source/Y90Source/gps/hist/point 0.0571 0.0426<br clear="none">/gate/source/Y90Source/gps/hist/point 0.1713 0.0518<br clear="none">/gate/source/Y90Source/gps/hist/point 0.2855 0.0594<br clear="none">/gate/source/Y90Source/gps/hist/point 0.3997 0.0649<br clear="none">/gate/source/Y90Source/gps/hist/point 0.5139 0.0686<br clear="none">/gate/source/Y90Source/gps/hist/point 0.6281 0.0708<br clear="none">/gate/source/Y90Source/gps/hist/point 0.7423 0.0717<br clear="none">/gate/source/Y90Source/gps/hist/point 0.8565 0.0715<br clear="none">/gate/source/Y90Source/gps/hist/point 0.9707 0.0704<br clear="none">/gate/source/Y90Source/gps/hist/point 1.0849 0.0685<br clear="none">/gate/source/Y90Source/gps/hist/point 1.1991 0.0657<br clear="none">/gate/source/Y90Source/gps/hist/point 1.3133 0.0619<br clear="none">/gate/source/Y90Source/gps/hist/point 1.4275 0.0569<br clear="none">/gate/source/Y90Source/gps/hist/point 1.5417 0.0507<br clear="none">/gate/source/Y90Source/gps/hist/point 1.6559 0.0430<br clear="none">/gate/source/Y90Source/gps/hist/point 1.7701 0.0342<br clear="none">/gate/source/Y90Source/gps/hist/point 1.8843 0.0246<br clear="none">/gate/source/Y90Source/gps/hist/point 1.9985 0.0150<br clear="none">/gate/source/Y90Source/gps/hist/point 2.1127 0.00643<br clear="none">/gate/source/Y90Source/gps/hist/point 2.2269 0.00113<br clear="none">/gate/source/Y90Source/gps/hist/inter Lin<br clear="none">/gate/source/Y90Source/gps/pos/type Point<br clear="none">/gate/source/Y90Source/gps/angtype iso<br clear="none">/gate/source/Y90Source/setForcedUnstableFlag true<br clear="none">/gate/source/Y90Source/setForcedHalfLife 230760. s</span></div><div><span class="yiv4927210009gmail-gD"><br clear="none"></span></div><div><span class="yiv4927210009gmail-gD">Thank you so much for your help again!</span></div><div><span class="yiv4927210009gmail-gD"><br clear="none"></span></div><div><span class="yiv4927210009gmail-gD">Maria<br clear="none"></span></div></div><br clear="none"><div class="yiv4927210009gmail_quote"><div class="yiv4927210009gmail_attr" dir="ltr">Ashok Tiwari <<a shape="rect" ymailto="mailto:tiwarias@yahoo.com" href="mailto:tiwarias@yahoo.com" target="_blank" rel="noreferrer noopener">tiwarias@yahoo.com</a>> escreveu no dia terça, 9/02/2021 à(s) 16:49:<br clear="none"></div><div class="yiv4927210009yqt1653335508" id="yiv4927210009yqt21932"><blockquote class="yiv4927210009gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;"><div><div style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:16px;"><div></div>
<div dir="ltr">Hi Maria, </div><div dir="ltr"><br clear="none"></div><div dir="ltr">You could start a simulation by setting a total number of primaries instead of time. Based on my experience with GATEv8.1, the transport of 20 million particles was found to provide sufficient statistical precision. See details in: <a shape="rect" href="https://meridian.allenpress.com/radiation-research/article-abstract/193/6/531/436368/The-Impact-of-Tissue-Type-and-Density-on-Dose?redirectedFrom=fulltext" target="_blank" rel="noreferrer noopener">The Impact of Tissue Type and Density on Dose Point Kernels for Patient-Specific Voxel-Wise Dosimetry: A Monte Carlo Investigation | Radiation Research | Allen Press</a></div><div><br clear="none"></div><div id="yiv4927210009gmail-m_-2011264584781472718ydp3b9ead41enhancr_card_3691519457" style="max-width:400px;font-family:YahooSans, Helvetica Neue, Segoe UI, Helvetica, Arial, sans-serif;"><a shape="rect" href="https://meridian.allenpress.com/radiation-research/article-abstract/193/6/531/436368/The-Impact-of-Tissue-Type-and-Density-on-Dose?redirectedFrom=fulltext" style="text-decoration:none;color:rgb(0,0,0);" target="_blank" rel="noreferrer noopener"><table style="max-width:400px;" cellspacing="0" cellpadding="0" border="0"><tbody><tr><td colspan="1" rowspan="1" width="400"><table style="max-width:400px;border-width:1px;border-style:solid;border-color:rgb(224,228,233);border-radius:2px;" width="100%" cellspacing="0" cellpadding="0" border="0"><tbody><tr><td colspan="1" rowspan="1"><table style="background:rgb(255,255,255) none repeat scroll 0% 0%;width:100%;max-width:400px;border-radius:0px 0px 2px 2px;border-top:1px solid rgb(224,228,233);" cellspacing="0" cellpadding="0" border="0"><tbody><tr><td colspan="1" rowspan="1" style="background-color:rgb(255,255,255);padding:16px 0px 16px 12px;vertical-align:top;border-radius:0px 0px 0px 2px;"></td><td colspan="1" rowspan="1" style="vertical-align:middle;padding:12px 24px 16px 12px;width:99%;font-family:YahooSans, Helvetica Neue, Segoe UI, Helvetica, Arial, sans-serif;border-radius:0px 0px 2px;"><h2 style="font-size:14px;line-height:19px;margin:0px 0px 6px;font-family:YahooSans, Helvetica Neue, Segoe UI, Helvetica, Arial, sans-serif;color:rgb(38,40,42);max-width:314px;">The Impact of Tissue Type and Density on Dose Point Kernels for Patient-...</h2><p style="font-size:12px;line-height:16px;margin:0px;color:rgb(151,155,167);"></p></td></tr></tbody></table></td></tr></tbody></table></td></tr></tbody></table></a></div><div><br clear="none"></div><div dir="ltr">Also, I recommend using ~1 keV energy cuts in dose point kernel simulations. Cuts can be set by using:</div><div dir="ltr"><br clear="none"></div><div dir="ltr"><div><div>/gate/physics/Gamma/SetCutInRegion world ... mm</div><div>/gate/physics/Electron/SetCutInRegion world ... mm</div><div><br clear="none"></div><div>/gate/physics/Gamma/SetCutInRegion my_phantom ... mm </div><div>/gate/physics/Electron/SetCutInRegion my_phantom .. mm</div><div><br clear="none"></div><div dir="ltr">Another important thing is an steplimiter, you could set it by using:</div><div dir="ltr"><br clear="none"></div><div>/gate/physics/SetMaxStepSizeInRegion my_phantom ... mm</div><div>/gate/physics/ActivateStepLimiter e-</div><div><br clear="none"></div></div><div dir="ltr">All other infos are described in the paper.</div><div dir="ltr"><br clear="none"></div><div dir="ltr">Regards,</div><div dir="ltr">Ashok</div><br clear="none"></div>
</div><div id="yiv4927210009gmail-m_-2011264584781472718yahoo_quoted_3771814349">
<div style="font-family:Helvetica, Arial, sans-serif;font-size:13px;color:rgb(38,40,42);">
<div>
On Tuesday, 9 February 2021, 09:19:23 am GMT-6, Maria Inês Ribeiro <<a shape="rect" ymailto="mailto:mid.ribeiro@campus.fct.unl.pt" href="mailto:mid.ribeiro@campus.fct.unl.pt" target="_blank" rel="noreferrer noopener">mid.ribeiro@campus.fct.unl.pt</a>> wrote:
</div>
<div><br clear="none"></div>
<div><br clear="none"></div>
<div><div id="yiv4927210009gmail-m_-2011264584781472718yiv8564264232"><div dir="ltr">Thank you so much. <div><br clear="none"></div><div>It was supposed to be 100000 s. I just killed the simulation because it was taking too long. </div><div>I swapped to the normal ion source definition since fastY90 can't be used for dosimetry. </div><div><br clear="none"></div><div>Do you recommend decreasing cut energy for photons and electrons?</div><div><br clear="none"></div><div>Kind regards, </div><div>Maria </div></div><br clear="none"><div><div dir="ltr">陈戈飞 <<a shape="rect" ymailto="mailto:xiaofeizaina@126.com" href="mailto:xiaofeizaina@126.com" target="_blank" rel="noreferrer noopener">xiaofeizaina@126.com</a>> escreveu no dia terça, 9/02/2021 à(s) 15:13:<br clear="none"></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;"><div><div id="yiv4927210009gmail-m_-2011264584781472718yiv8564264232gmail-m_-973855092368961664contentDescription" style="line-height:1.5;text-align:justify;">
<div id="yiv4927210009gmail-m_-2011264584781472718yiv8564264232gmail-m_-973855092368961664imail_extraContent">Hi maria
<br clear="none">The total time of your simulations is 10000s，right? If so, you only need to multiply by 1000(Gy->mGy) and divide 10000(time to 1s).
<br clear="none">And I recommend setting cut energy for beta and photon. Default is 10 mm.
<br clear="none">You'd better recheck the definition of fastY90. It can't be used for dose purposes.
<br clear="none">Thank you
<br clear="none"></div>
<div><br clear="none"></div>
<div><br clear="none"></div>
<div>
<div id="yiv4927210009gmail-m_-2011264584781472718yiv8564264232gmail-m_-973855092368961664imail_signature">
<div style="margin:0px;padding:0px;">
<a shape="rect" href="https://maas.mail.163.com/dashi-web-extend/html/proSignature.html?iconUrl=https%3A%2F%2Fmail-online.nosdn.127.net%2Fqiyelogo%2FdefaultAvatar.png&name=%E9%99%88%E6%88%88%E9%A3%9E&uid=example%40163.com&ftlId=3&items=%5B%22%E9%82%AE%E7%AE%B1%EF%BC%9Axiaofeizaina%40126.com%22%5D" style="display:block;max-width:400px;background:rgb(255,255,255) none repeat scroll 0% 0%;padding:15px 0px 10px;text-decoration:none;outline:currentcolor none medium;" target="_blank" rel="noreferrer noopener">
<table style="width:100%;max-width:100%;table-layout:fixed;border-collapse:collapse;border-spacing:0px;line-height:1.3;color:rgb(155,158,161);font-size:14px;" cellpadding="0"><tbody style="font-family:Sans Micro verdana;"><tr><td colspan="1" rowspan="1" style="padding:0px 7px 0px 0px;width:45px;" width="45">
<img style="width:38px;min-height:38px;border-radius:50%;" src="https://mail-online.nosdn.127.net/qiyelogo/defaultAvatar.png" width="38" height="38">
</td><td colspan="1" rowspan="1" style="padding:0px 0px 0px 7px;">
<div style="max-width:380px;">
<div style="padding-right:35px;font-size:16px;margin-bottom:5px;color:rgb(49,53,59);font-weight:bold;width:100%;white-space:nowrap;overflow:hidden;">陈戈飞</div>
<div style="font-size:0px;line-height:0;">
</div>
<div>
邮箱：xiaofeizaina@126.com
</div>
</div>
</td></tr></tbody></table>
</a>
</div>
<div><p style="border-top:1px solid rgb(229,229,229);padding-top:8px;font-size:12px;color:rgb(182,184,187);line-height:1.833;">签名由 <a shape="rect" href="https://mail.163.com/dashi/dlpro.html?from=mail88" style="color:rgb(106,168,246);text-decoration:none;" target="_blank" rel="noreferrer noopener">网易邮箱大师</a> 定制</p></div></div>
</div>
<div style="background-color:rgb(242,242,242);color:black;padding-top:6px;padding-bottom:6px;border-radius:3px;margin-top:45px;margin-bottom:20px;"><div style="font-size:14px;line-height:1.5;margin-left:10px;margin-right:10px;">On <span>02/09/2021 22:38</span>, <a shape="rect" style="text-decoration:none;color:rgb(42,151,255);" ymailto="mailto:mid.ribeiro@campus.fct.unl.pt" href="mailto:mid.ribeiro@campus.fct.unl.pt" target="_blank" rel="noreferrer noopener">Maria Inês Ribeiro</a> wrote: </div></div><blockquote id="yiv4927210009gmail-m_-2011264584781472718yiv8564264232gmail-m_-973855092368961664ntes-iosmail-quote" style="margin:0px;"><div dir="ltr">Hi everyone, <div><br clear="none"></div><div>I am attempting to compute a dose point kernel in soft tissue for a point source of Y90. I am using the fastY90 with an activity of 1 MBq. </div><div>I then ran 40.000 slices of 1 second and stored the dose. </div><div><br clear="none"></div><div>My question is how do I convert my dose map (Gy) into the s values (mGy/(MBq.s)) for each voxel? </div><div><br clear="none"></div><div>Here is my could in case it could help to support my question:</div><div><br clear="none"></div><div>#===============================================================================<br clear="none"># GEOMETRY<br clear="none">#===============================================================================<br clear="none"><br clear="none">/gate/geometry/setMaterialDatabase data/GateMaterials.db<br clear="none"><br clear="none">#WORLD ======================================<br clear="none">/gate/world/setMaterial Air<br clear="none">/gate/world/geometry/setXLength 50 mm<br clear="none">/gate/world/geometry/setYLength 50 mm<br clear="none">/gate/world/geometry/setZLength 50 mm<br clear="none"><br clear="none">#LIVER ======================================<br clear="none">/gate/world/daughters/name liverPhantom<br clear="none">/gate/world/daughters/insert box<br clear="none">/gate/liverPhantom/setMaterial Liver<br clear="none">/gate/liverPhantom/geometry/setXLength 24.31 mm<br clear="none">/gate/liverPhantom/geometry/setYLength 24.31 mm<br clear="none">/gate/liverPhantom/geometry/setZLength 24.31 mm<br clear="none"><br clear="none">#===============================================================================<br clear="none"># PHYSICS<br clear="none">#===============================================================================<br clear="none"><br clear="none">/gate/physics/addPhysicsList emstandard_opt3<br clear="none"><br clear="none">#===============================================================================<br clear="none"># DETECTORS<br clear="none">#===============================================================================<br clear="none"><br clear="none">#DOSE ACTOR ==============================<br clear="none">#builds 3D images of:<br clear="none"># - the energy deposited (edep in MeV)<br clear="none"># - dose deposited (Gy)<br clear="none"># - the number of hits in a volume<br clear="none"><br clear="none">/gate/actor/addActor DoseActor VSV<br clear="none">/gate/actor/VSV/attachTo liverPhantom<br clear="none">/gate/actor/VSV/stepHitType random<br clear="none">/gate/actor/VSV/setSize 24.31 24.31 24.31 mm<br clear="none">/gate/actor/VSV/setPosition 0 0 0 cm<br clear="none">/gate/actor/VSV/setResolution 11 11 11<br clear="none">/gate/actor/VSV/enableEdep true<br clear="none">/gate/actor/VSV/enableUncertaintyEdep true<br clear="none">/gate/actor/VSV/enableSquaredEdep true<br clear="none">/gate/actor/VSV/enableDose true<br clear="none">/gate/actor/VSV/enableUncertaintyDose true<br clear="none">/gate/actor/VSV/enableSquaredDose true<br clear="none">/gate/actor/VSV/enableNumberOfHits true<br clear="none">/gate/actor/VSV/save output/VSV.mhd<br clear="none"><br clear="none">#SIMULATION STATISTIC =====================<br clear="none">#this actos counts the number of steps, tracks, events, runs in the simulation<br clear="none"><br clear="none">/gate/actor/addActor SimulationStatisticActor stat<br clear="none">/gate/actor/stat/save output/stat-VSV.txt<br clear="none">/gate/actor/stat/saveEveryNSeconds 60<br clear="none"><br clear="none"><br clear="none">#===============================================================================<br clear="none"># INITIALISATION<br clear="none">#===============================================================================<br clear="none"><br clear="none">/gate/run/initialize<br clear="none"><br clear="none">/gate/random/setEngineName MersenneTwister<br clear="none">/gate/random/setEngineSeed auto<br clear="none"><br clear="none">#===========================================================<br clear="none"># SOURCE <br clear="none">#===========================================================<br clear="none"><br clear="none"># fastY90<br clear="none">/gate/source/addSource y90 fastY90<br clear="none">/gate/source/y90/setActivity 1. MBq<br clear="none"><br clear="none">#=====================================================<br clear="none"># START<br clear="none">#=====================================================<br clear="none"><br clear="none">/gate/application/setTimeSlice 1. s<br clear="none">/gate/application/setTimeStart 0. s<br clear="none">/gate/application/setTimeStop 1000000. s <br clear="none">/gate/application/start<br clear="none"></div><div><br clear="none"></div><div><br clear="none"></div><div>Thanks a lot for your attention, </div><div><br clear="none"></div><div>Maria Inês</div><div><br clear="none"></div><div><br clear="none"></div><div><br clear="none"></div><div><br clear="none"></div><div><br clear="none"></div></div>
</blockquote>
</div></div></blockquote></div>
</div>_______________________________________________<br clear="none">Gate-users mailing list<br clear="none"><a shape="rect" ymailto="mailto:Gate-users@lists.opengatecollaboration.org" href="mailto:Gate-users@lists.opengatecollaboration.org" target="_blank" rel="noreferrer noopener">Gate-users@lists.opengatecollaboration.org</a><br clear="none"><a shape="rect" href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank" rel="noreferrer noopener">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a></div>
</div>
</div></div></blockquote></div></div>
</div></div></div>
</div>
</div></body></html>