[Gate-users] initializing Radioactive Decay in Gate v7.0, Geant4 v.9.6.4

[APSARA TENNAKOON]

Dear GATE users,

I am trying to look at secondary particles produced when a (e.g. brain)
target in a water phantom is irradiated with a proton beam.

I am using Gate v7.0, Geant4 v.9.6.4 on my university's Linux server.

I have chosen the physics list QGSP_BIC_HP (Arce, P., et al., 2020).

A secondary production actor attached to the target shows β+ emitters
produced (O15, N13, C11, etc.). However, I am not getting any "e+ Decay" or
"gamma Decay", which tells me the decay of these isotopes must not have
been simulated.

I have located the dataset 'RadioactiveDecay3.6'
in geant4/geant4.9.6.p04-install/share/Geant4-9.6.4/data/ and have verified
that the environment variable G4RADIOACTIVEDATA points to this.

The list of enabled/ initialized processes for the simulation shows that
RadioactiveDecay has been "enabled", but not "initialized".

I would appreciate any help understanding what might be causing this.

#=====================================================
# PHYSICS
#=====================================================
/gate/physics/addPhysicsList QGSP_BIC_HP

/gate/physics/addProcess    Decay
/gate/physics/addProcess    RadioactiveDecay

/gate/physics/Gamma/SetCutInRegion      world 10 mm
/gate/physics/Electron/SetCutInRegion      world 10 mm
/gate/physics/Positron/SetCutInRegion      world 10 mm
/gate/physics/Proton/SetCutInRegion         world 10 mm

/gate/physics/Gamma/SetCutInRegion    phantom 1 mm
/gate/physics/Electron/SetCutInRegion   phantom 1 mm
/gate/physics/Positron/SetCutInRegion   phantom 1 mm
/gate/physics/Proton/SetCutInRegion      phantom 1 mm

/gate/physics/Gamma/SetCutInRegion   target 0.05 mm
/gate/physics/Electron/SetCutInRegion   target 0.05 mm
/gate/physics/Positron/SetCutInRegion   target 0.05 mm
/gate/physics/Proton/SetCutInRegion      target 0.05 mm

/gate/physics/displayCuts
/gate/physics/processList Enabled
/gate/physics/processList Initialized
/gate/physics/print ../output/physics.txt

#=====================================================
# SecondaryProductionActor
#=====================================================
/gate/actor/addActor SecondaryProductionActor   SecnP
/gate/actor/SecnP/save                  ../output/SecnP.root
/gate/actor/SecnP/attachTo              target

Thank you!

Best regards,
Apsara
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210303/bfc79053/attachment-0001.html>