<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:rgb(32,18,77)">Hello David,<br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:rgb(32,18,77)">Thank you ! I had a feeling actors were playing a role in this but I suppose I read the chapter a bit too fast.<br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:rgb(32,18,77)">One thing I hadn't tried, because I had a odd fear that someting was going to crash, is to edit an output file which is being updated during the simulation (stat.txt for instance) : for those who had the same feeling, you can do it it's safe ! :) You just have to reload it when you want updates. (geany automatically asks if you want to reload it )<br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:rgb(32,18,77)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:rgb(32,18,77)">And thanks for the tip about reaching a desired uncertainty in a chosen region, I'm sure I'll be using this soon.<br><br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:rgb(32,18,77)">Another trick was sent to me in private message : using the timeSlice macos <br><br>As I'm working on static dosimetry applications, I thought this was dedicated to dynamic-imaging, but you can make a good use of it to know where you are in the simulation process.<br><br><div>
<div>/gate/application/setTimeSlice 1 s</div>
<div>/gate/application/setTimeStart 0 s</div>
<div>/gate/application/setTimeStop 5 s<br></div><div><br>As it is said in the User Guide : <br><br>"<i>User can fix the total number of events. In this case, the number of particles is splitted between slices in function of the <span>time</span> of each slice</i>."<br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:rgb(32,18,77)">For instance if you send 50 000 particles, with a setTimeStop of 5 sec and a setTimeSlice of 1 sec, you'll have a console output every 10 000 particles.<br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:rgb(32,18,77)">Thanks again for the help<br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:rgb(32,18,77)">Matthieu<br></div><br><br><br><div class="gmail_quote">2015-04-02 10:56 GMT+02:00 David Sarrut <span dir="ltr"><<a href="mailto:David.Sarrut@creatis.insa-lyon.fr" target="_blank">David.Sarrut@creatis.insa-lyon.fr</a>></span>:<br><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">Hello Matthieu, <div><br></div><div>- use the option saveEveryNSeconds of some actor to allow dump results at regular interval, like : '/gate/actor/mydoseactor/saveEveryNSeconds 60'</div><div><br></div><div>- no way to stop at a given uncertainty. A good strategy : perform 3 or 4 simulations with low and increasing nb of particles. For each, compute uncertainty in the region you are interested in. Plot and fit a 1/sqrt(n) curve. Look the n you need to reach your desired uncertainty. Relauch simulation with n. </div><div><br></div><div>- for activity mode, you can add a simple SimulationStatisticActor to get some useful information at regular time interval</div><div><div>/gate/actor/addActor SimulationStatisticActor stat</div><div>/gate/actor/stat/save output/stat.txt</div><div>/gate/actor/stat/saveEveryNSeconds 60</div></div><div><br></div><div>David</div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Wed, Apr 1, 2015 at 3:32 PM, Matthieu Dondey <span dir="ltr"><<a href="mailto:m.dondey@gmail.com" target="_blank">m.dondey@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div><div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;color:rgb(32,18,77)"><span style="font-family:"Trebuchet MS","sans-serif";color:rgb(32,18,77)" lang="EN-US">Dear Gaters,<br><br></span>
<p class="MsoNormal" style="margin-bottom:12pt;text-align:left;line-height:normal" align="left"><span style="font-family:"Trebuchet MS","sans-serif";color:rgb(32,18,77)" lang="EN-US">Maybe it's a silly question, or it has already been
answered but I'm running long simulations and I can't find any clues about the
progression of an ongoing simulation.</span></p>
<p class="MsoNormal" style="margin-bottom:12pt;text-align:left;line-height:normal" align="left"><span style="font-family:"Trebuchet MS","sans-serif";color:rgb(32,18,77)" lang="EN-US">Does someone have a trick to get this information ?
maybe in % relatively to the TotalNumberOfPrimaries ?</span></p>
<p class="MsoNormal" style="margin-bottom:12pt;text-align:left;line-height:normal" align="left"><span style="font-family:"Trebuchet MS","sans-serif";color:rgb(32,18,77)" lang="EN-US">I've tried using the verbose commands, and the
/gate/verbose Beam command gives output of every generated particle ... but
flooding the console with outputs doesn't help with the simulation time :/</span></p>
<p class="MsoNormal" style="margin-bottom:12pt;text-align:left;line-height:normal" align="left"><span style="font-family:"Trebuchet MS","sans-serif";color:rgb(32,18,77)" lang="EN-US">I've read in <i>"A review of the use and
potential of the GATE Monte Carlo simulation code for radiation therapy and
dosimetry applications" by David Sarrut</i> a table with time and number
of particles needed to get a fixed uncertainty (2% in this table).</span></p>
<p class="MsoNormal" style="margin-bottom:12pt;text-align:left;line-height:normal" align="left"><span style="font-family:"Trebuchet MS","sans-serif";color:rgb(32,18,77)" lang="EN-US">For cases where the source isn't defined with
TotalNumberOfPrimaries, but in activity mode, I'd also be happy to know how to
get this type of information (time and number of particles needed to get a
fixed uncertainty) if anyone has a simple way to get it :)</span></p><p class="MsoNormal" style="margin-bottom:12pt;text-align:left;line-height:normal" align="left"><span style="font-family:"Trebuchet MS","sans-serif";color:rgb(32,18,77)" lang="EN-US"></span></p><p class="MsoNormal" style="margin-bottom:12pt;text-align:left;line-height:normal" align="left"><span style="font-family:"Trebuchet MS","sans-serif";color:rgb(32,18,77)" lang="EN-US">Thanks in advance<span><font color="#888888"><br></font></span></span></p><span><font color="#888888"><p class="MsoNormal" style="margin-bottom:12pt;text-align:left;line-height:normal" align="left"><span style="font-family:"Trebuchet MS","sans-serif";color:rgb(32,18,77)" lang="EN-US"></span></p><p class="MsoNormal" style="margin-bottom:12pt;text-align:left;line-height:normal" align="left"><span style="font-family:"Trebuchet MS","sans-serif";color:rgb(32,18,77)" lang="EN-US">Matthieu<br></span></p>
</font></span></div></div>
<br></div></div><span>_______________________________________________<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" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br></span></blockquote></div><span><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div><div dir="ltr">David Sarrut, Phd<br>Directeur de recherche CNRS<br>CREATIS, UMR CNRS 5220, Inserm U 1044<div>Centre de lutte contre le cancer Léon Bérard<br>28 rue Laënnec, 69373 Lyon cedex 08<br>Tel : 04 78 78 51 51 / 06 74 72 05 42<br><a href="http://www.creatis.insa-lyon.fr/~dsarrut" target="_blank">http://www.creatis.insa-lyon.fr/~dsarrut</a><br>_________________________________</div><div> "2 + 2 = 5, for extremely large values of 2"<br>_________________________________</div></div></div>
</font></span></div>
</blockquote></div><br></div></div>