[Gate-users] Radioactive decay with general ion source and new physics list

[Grevillot Loic]

Hi Albert,

A physics-list includes a list of physics processes (hadronics and electromagnetics) and parameters (cut, step lengths, etc.) adapted to a type of application: medical, micro-dosimetry, LHC, etc.
The list of physics-processes included in a Geant4 physics-list (builder) shall be presented on Geant4 website (and probably some documentation also on Gate website).
Apparently, the process you need in that case is not part of this physics-list, but as you did you can add it manually in Gate.
Otherwise, you might look for a Geant4 physics-list including directly this process, or build your own...

Sincerely,

Loïc



-----Original Message-----
From: gate-users-bounces at lists.opengatecollaboration.org [mailto:gate-users-bounces at lists.opengatecollaboration.org] On Behalf Of Albert Hirtl
Sent: Dienstag, 2. September 2014 13:46
To: gate-users
Subject: [Gate-users] Radioactive decay with general ion source and new physics list

Dear Colleagues,

I have a question regarding the new physics list mechanism in GATE v7.0.
If I understood it correctly I do not need to specify the various processes explicitly if I use, e.g.,

/gate/physics/addPhysicsList emstandard_opt3

However, if I run a simulation with the general ion source

/gate/source/addSource fluor18
/gate/source/fluor18/setActivity 500 becquerel /gate/source/fluor18/gps/particle ion /gate/source/fluor18/gps/ion 9 18 0 0 /gate/source/fluor18/gps/energytype Mono /gate/source/fluor18/gps/monoenergy 0. MeV /gate/source/fluor18/setForcedUnstableFlag true /gate/source/fluor18/setForcedHalfLife 6586 s

and if I start the simulation with

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

I do not get any decaying particles. In order to get a signal I have to add

/gate/physics/addProcess RadioactiveDecay

So, isn't the advantage of the new mechanism, that I do not have to define the processes explicitly? I guess if I ran the simulation with

/gate/application/setTotalNumberOfPrimaries

everything would be fine, right?

Additionally, if I use

/gate/physics/addPhysicsList  QGSP_BERT_EMY

I do not get any decay products no matter what I do...

Any comment would be very much appreciated.

Cheers,
Albert

--
Medizinische Universität Wien
Universitätsklinik für Radiologie und Nuklearmedizin Währinger Gürtel 18-20
1090 Wien
Austria

Tel.: +43 (0) 1 / 40400 55450

Web:
http://www.meduniwien.ac.at/user/albert.hirtl
http://www.meduniwien.ac.at/moccamed/

Disclaimer:
Please notice our E-Mail Disclaimer http://www.medaustron.at/email-disclaimer/