[Gate-users] Asking for help on GATE program used for Linear Accelerator radiation therapy shielding study
Heidie Jayo
jayoheidie at gmail.com
Tue May 7 11:52:07 CEST 2019
Good Day Everyone,
I would like to ask for help regarding the code attached below. I am using
it to do a GATE Monte Carlo Simulation on the bunker shielding of our
Linear Accelerator with 10 MV Photon energy. I would like to get the
results of the dose outside the primary and secondary barrier walls but I
am having difficulty with getting the results.
I would like to ask for your advice if the settings used for the physics,
cut, output-detector and beams are correct.
Hoping for your help. Thank you very much.
Sincerely,
Heidie Jayo
Below is the code:
#=====================================================
# VERBOSE and VISUALISATION
#=====================================================
/control/execute verbose.mac
/control/execute visu.mac
#=====================================================
# GEOMETRY
#=====================================================
/gate/geometry/setMaterialDatabase data/GateMaterials.db
# WORLD
/gate/world/setMaterial Air
/gate/world/geometry/setXLength 35.0 m
/gate/world/geometry/setYLength 30.0 m
/gate/world/geometry/setZLength 17.0 m
# vPHANTOM
/gate/world/daughters/name vphantom
/gate/world/daughters/insert box
/gate/vphantom/setMaterial Air
/gate/vphantom/geometry/setXLength 30.0 m
/gate/vphantom/geometry/setYLength 25.0 m
/gate/vphantom/geometry/setZLength 15.0 m
/gate/vphantom/placement/setTranslation 0 0 0 cm
#/gate/vphantom/vis/forceSolid
/gate/vphantom/vis/setColor blue
# WATERPHANTOM
#/gate/vphantom/daughters/name waterphantom
#/gate/vphantom/daughters/insert box
#/gate/waterphantom/setMaterial Water
#/gate/waterphantom/geometry/setXLength 40.0 cm
#/gate/waterphantom/geometry/setYLength 40.0 cm
#/gate/waterphantom/geometry/setZLength 40.0 cm
#/gate/waterphantom/placement/setTranslation 0 0 0 cm
#/gate/waterphantom/vis/forceSolid
#/gate/waterphantom/vis/setColor yellow
# Volume 1
/gate/vphantom/daughters/name volume1
/gate/vphantom/daughters/insert box
/gate/volume1/setMaterial Concrete
/gate/volume1/geometry/setXLength 6.2 m
/gate/volume1/geometry/setYLength 2.53 m
/gate/volume1/geometry/setZLength 6 m
/gate/volume1/placement/setTranslation 0 4.915 -0.65 m
/gate/volume1/vis/forceSolid
/gate/volume1/vis/setColor white
[Note: I have excluded the other volumes of the bunker shielding walls
since they are mostly of similar geometry]
#=====================================================
# PHYSICS
#=====================================================
/gate/physics/addPhysicsList QGSP_BIC_HP
/gate/physics/displayCuts
/gate/physics/print output/physics.txt
#=====================================================
# CUTS
#=====================================================
#/gate/physics/Gamma/SetCutInRegion world 1 mm
#/gate/physics/Electron/SetCutInRegion world 1 mm
#/gate/physics/Positron/SetCutInRegion world 1 mm
#/gate/physics/Gamma/SetCutInRegion target 0.1 mm
#/gate/physics/Electron/SetCutInRegion target 0.1 mm
#/gate/physics/Positron/SetCutInRegion target 0.1 mm
/gate/physics/SetMaxStepSizeInRegion vphantom 1 mm
#/gate/physics/ActivateStepLimiter proton
#/gate/physics/ActivateStepLimiter deuteron
#/gate/physics/ActivateStepLimiter triton
#/gate/physics/ActivateStepLimiter alpha
#/gate/physics/ActivateStepLimiter GenericIon
/gate/physics/ActivateStepLimiter gamma
/gate/physics/ActivateStepLimiter e-
/gate/physics/ActivateStepLimiter e+
/gate/physics/ActivateStepLimiter neutron
# Verbose physic list (uncomment to display)
/gate/physics/displayCuts
/gate/physics/processList Enabled
/gate/physics/processList Initialized
#=====================================================
# OUTPUT - DETECTOR
#=====================================================
/gate/actor/addActor SimulationStatisticActor stat
/gate/actor/stat/save output/stat.txt
#/gate/actor/stat/saveEveryNSeconds 60
# ----------------------------------------------------
# the following actor stores 3D distributions of
# dose/edep/uncertainty/nbHit into files (mhd image file format)
/gate/actor/addActor DoseActor
doseDistribution
/gate/actor/doseDistribution/save output/dose.mhd
/gate/actor/doseDistribution/attachTo vphantom
/gate/actor/doseDistribution/stepHitType random
#/gate/actor/doseDistribution/setSize 20000 20000 20000 mm
/gate/actor/doseDistribution/setVoxelSize 20 20 20 cm
/gate/actor/doseDistribution/setPosition 0 0 0 mm
/gate/actor/doseDistribution/enableEdep false
/gate/actor/doseDistribution/enableUncertaintyEdep false
/gate/actor/doseDistribution/enableDose true
/gate/actor/doseDistribution/enableSquaredDose false
/gate/actor/doseDistribution/enableUncertaintyDose true
/gate/actor/doseDistribution/enableNumberOfHits false
/gate/actor/doseDistribution/saveEveryNSeconds 20
/gate/actor/addActor ProductionAndStoppingActor sGammaMap
/gate/actor/sGammaMap/save output/sGamma.root
/gate/actor/sGammaMap/attachTo vphantom
/gate/actor/sGammaMap/stepHitType random
/gate/actor/sGammaMap/setVoxelSize 18000 18000 18000 mm
/gate/actor/sGammaMap/addFilter particleFilter
/gate/actor/sGammaMap/particleFilter/addParticle gamma
/gate/actor/sGammaMap/saveEveryNSeconds 20
/gate/actor/addActor ProductionAndStoppingActor sElectronMap
/gate/actor/sElectronMap/save output/sElectron.root
/gate/actor/sElectronMap/attachTo vphantom
/gate/actor/sElectronMap/stepHitType random
/gate/actor/sElectronMap/setVoxelSize 200 200 200 mm
/gate/actor/sElectronMap/addFilter particleFilter
/gate/actor/sElectronMap/particleFilter/addParticle e-
/gate/actor/sElectronMap/saveEveryNSeconds 20
/gate/actor/addActor ProductionAndStoppingActor sPositronMap
/gate/actor/sPositronMap/save output/sPositron.root
/gate/actor/sPositronMap/attachTo vphantom
/gate/actor/sPositronMap/stepHitType random
/gate/actor/sPositronMap/setVoxelSize 18000 18000 18000 mm
/gate/actor/sPositronMap/addFilter particleFilter
/gate/actor/sPositronMap/particleFilter/addParticle e+
/gate/actor/sPositronMap/saveEveryNSeconds 20
/gate/actor/addActor ProductionAndStoppingActor sNeutronMap
/gate/actor/sNeutronMap/save output/sNeutron.root
/gate/actor/sNeutronMap/attachTo vphantom
/gate/actor/sNeutronMap/stepHitType random
/gate/actor/sNeutronMap/setVoxelSize 18000 18000 18000 mm
/gate/actor/sNeutronMap/addFilter particleFilter
/gate/actor/sNeutronMap/particleFilter/addParticle neutron
/gate/actor/sNeutronMap/saveEveryNSeconds 20
### PhaseSpaceActor
/gate/actor/addActor PhaseSpaceActor
phantomPhaseSpaceActor
/gate/actor/phantomPhaseSpaceActor/save
output/phantomPhaseSpace.root
/gate/actor/phantomPhaseSpaceActor/attachTo vphantom
/gate/actor/phantomPhaseSpaceActor/storeSecondaries true
/gate/actor/phantomPhaseSpaceActor/enableEkine true
/gate/actor/phantomPhaseSpaceActor/enableXPosition true
/gate/actor/phantomPhaseSpaceActor/enableYPosition true
/gate/actor/phantomPhaseSpaceActor/enableZPosition true
/gate/actor/phantomPhaseSpaceActor/enableXDirection true
/gate/actor/phantomPhaseSpaceActor/enableYDirection true
/gate/actor/phantomPhaseSpaceActor/enableZDirection true
/gate/actor/phantomPhaseSpaceActor/enableProductionVolume true
/gate/actor/phantomPhaseSpaceActor/enableProductionProcess true
/gate/actor/phantomPhaseSpaceActor/enableParticleName true
/gate/actor/phantomPhaseSpaceActor/enableWeight false
/gate/actor/phantomPhaseSpaceActor/enableTime false
#/gate/actor/phantomPhaseSpaceActor/storeAllStep true ##
comment when using storeOutgoingParticles
#=====================================================
# INITIALISATION
#=====================================================
/gate/run/initialize
#=====================================================
# BEAMS
#=====================================================
# Gamma beam
/gate/source/addSource mybeam gps
/gate/source/mybeam/gps/particle gamma
/gate/source/mybeam/gps/ene/type Gauss
/gate/source/mybeam/gps/ene/mono 10 MeV
/gate/source/mybeam/gps/ene/sigma 0.1274 MeV
/gate/source/mybeam/gps/pos/centre 0 0 110 cm
/gate/source/mybeam/gps/direction 0 0 -1
/gate/source/mybeam/visualize 1 green
#=====================================================
# START BEAMS
#=====================================================
/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto
/gate/application/setTotalNumberOfPrimaries 100000000
/gate/application/start
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20190507/2eddd867/attachment-0001.html>
More information about the Gate-users
mailing list