[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