[Gate-users] Progression of an ongoing simulation

[Matthieu Dondey]

Hello David,

Thank you ! I had a feeling actors were playing a role in this but I
suppose I read the chapter a bit too fast.

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 )

And thanks for the tip about reaching a desired uncertainty in a chosen
region, I'm sure I'll be using this soon.


Another trick was sent to me in private message : using the timeSlice
macos

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.

/gate/application/setTimeSlice   1 s
/gate/application/setTimeStart   0   s
/gate/application/setTimeStop    5 s

As it is said in the User Guide :

"*User can fix the total number of events. In this case, the number of
particles is splitted between slices in function of the time of each slice*
."


​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.

Thanks again for the help

Matthieu



2015-04-02 10:56 GMT+02:00 David Sarrut <David.Sarrut at creatis.insa-lyon.fr>:

> Hello Matthieu,
>
> - use the option saveEveryNSeconds of some actor to allow dump results at
> regular interval, like : '/gate/actor/mydoseactor/saveEveryNSeconds     60'
>
> - 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.
>
> - for activity mode, you can add a simple SimulationStatisticActor to get
> some useful information at regular time interval
> /gate/actor/addActor               SimulationStatisticActor stat
> /gate/actor/stat/save              output/stat.txt
> /gate/actor/stat/saveEveryNSeconds 60
>
> David
>
>
>
>
> On Wed, Apr 1, 2015 at 3:32 PM, Matthieu Dondey <m.dondey at gmail.com>
> wrote:
>
>> ​Dear Gaters,
>>
>> 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.
>>
>> Does someone have a trick to get this information ? maybe in % relatively
>> to  the TotalNumberOfPrimaries ?
>>
>> 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 :/
>>
>> I've read in *"A review of the use and potential of the GATE Monte Carlo
>> simulation code for radiation therapy and dosimetry applications" by David
>> Sarrut* a table with time and number of particles needed to get a fixed
>> uncertainty (2% in this table).
>>
>> 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 :)
>>
>> Thanks in advance
>>
>> Matthieu
>>
>> _______________________________________________
>> Gate-users mailing list
>> Gate-users at lists.opengatecollaboration.org
>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>>
>
>
>
> --
> David Sarrut, Phd
> Directeur de recherche CNRS
> CREATIS, UMR CNRS 5220, Inserm U 1044
> Centre de lutte contre le cancer Léon Bérard
> 28 rue Laënnec, 69373 Lyon cedex 08
> Tel : 04 78 78 51 51 / 06 74 72 05 42
> http://www.creatis.insa-lyon.fr/~dsarrut
> _________________________________
>  "2 + 2 = 5,  for extremely large values of 2"
> _________________________________
>
-------------- section suivante --------------
Une pi�ce jointe HTML a �t� nettoy�e...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20150402/e5babbc4/attachment.html>