[Gate-users] PBS Source-Surface-Distance (SSD)
Awad Almarhaby
mph-awad at hotmail.com
Tue Oct 11 16:10:43 CEST 2016
Dears Gate users,
I tried to run simulations for one proton PBS at different source to surface distances (nozzle-patient-distances) to investigate the effect of SSD variation on the dose distribution in water phantom. SSD were set to be (100,75,50,25 and 5 cm) at every run I changed the source position in the z-direction. Dose was normalised to the maximum dose. But, I got the same distribution for the first four values which I believe is not correct. Also, I did a similar simulation for proton beam with a fix energy and different diameters (changing sigma_x and sigma_y) and I got the same dose distribution with a slight change. Until now I am unable to detect the error in my set-up.
I attached a copy of the main macro file and I would be grateful if you could help me to find out the error source.
Thanks in advance,
Awad
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
# Example of electron beam in phantom (water)
#-------------------oooooOOOOO00000OOOOOooooo---------------------#
#=====================================================
# VERBOSE and VISUALISATION
#=====================================================
/control/execute mac/verbose.mac
#/control/execute mac/visu.mac
#=====================================================
# GEOMETRY
#=====================================================
/gate/geometry/setMaterialDatabase data/GateMaterials.db
/gate/geometry/setMaterialDatabase data/patient-HUmaterials.db
# WORLD
/gate/world/setMaterial Air
/gate/world/geometry/setXLength 5.0 m
/gate/world/geometry/setYLength 5.0 m
/gate/world/geometry/setZLength 5.0 m
#/gate/world/placement/setTranslation 0 0 0 mm
# WATERBOX
/gate/world/daughters/name water_box
/gate/world/daughters/insert box
/gate/water_box/setMaterial Water
/gate/water_box/geometry/setXLength 30 cm
/gate/water_box/geometry/setYLength 30 cm
/gate/water_box/geometry/setZLength 50 cm
/gate/water_box/placement/setTranslation 0.0 0.0 0.0 cm
/gate/water_box/vis/setColor blue
#=====================================================
# PHYSICS
#=====================================================
/gate/physics/addPhysicsList emstandard_opt3
/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 water_box 0.1 mm
/gate/physics/Electron/SetCutInRegion water_box 0.1 mm
/gate/physics/Positron/SetCutInRegion water_box 0.1 mm
/gate/physics/SetMaxStepSizeInRegion world 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 proton
/gate/physics/ActivateStepLimiter e-
#=====================================================
# Statistics actor
#=====================================================
/gate/actor/addActor SimulationStatisticActor stat
/gate/actor/stat/saveEveryNSeconds 60
/gate/actor/stat/save output/statistics.txt
#=====================================================
# Dose Actor
#=====================================================
/gate/actor/addActor DoseActor dose
/gate/actor/dose/save output/120.txt
/gate/actor/dose/attachTo water_box
/gate/actor/dose/stepHitType random
/gate/actor/dose/setPosition 0 0 0 cm
/gate/actor/dose/setResolution 50 50 500
/gate/actor/dose/setSize 5 5 50 cm
/gate/actor/dose/saveEveryNSeconds 60
/gate/actor/dose/enableEdep true
/gate/actor/dose/enableUncertaintyEdep true
/gate/actor/dose/enableDose true
/gate/actor/dose/enableUncertaintyDose true
/gate/actor/dose/enableNumberOfHits true
#=====================================
# SOURCE
#=====================================
/gate/source/addSource PBS PencilBeam
/gate/source/PBS/setPosition 0 0 -100 cm
/gate/source/PBS/setEnergy 120 MeV
/gate/source/PBS/setSigmaEnergy 0.1 MeV
/gate/source/PBS/setSigmaX 3 mm
/gate/source/PBS/setSigmaY 3 mm
/gate/source/PBS/setSigmaTheta 1 mrad
/gate/source/PBS/setSigmaPhi 1 mrad
/gate/source/PBS/setEllipseXThetaEmittance 3 mm*mrad
/gate/source/PBS/setEllipseXThetaRotationNorm positive
/gate/source/PBS/setEllipseYPhiEmittance 3 mm*mrad
/gate/source/PBS/setEllipseYPhiRotationNorm positive
/gate/source/PBS/setParticleType proton
#/gate/source/PBS/setRotationAxis 1 0 0
#/gate/source/PBS/setRotationAngle 90 deg
#=====================================================
# INITIALISATION
#=====================================================
/gate/run/initialize
#=====================================================
# START BEAMS
#=====================================================
/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto
# When only using the EmCalculator actor for various configurations
# (without caring of the dose deposition)
#/gate/application/setTotalNumberOfPrimaries 1
# Otherwise
/gate/application/setTotalNumberOfPrimaries 1000000
/gate/application/start
Sent from Outlook<http://aka.ms/weboutlook>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20161011/e1053262/attachment.html>
More information about the Gate-users
mailing list