[Gate-users] Multiple Fields
Awad Almarhaby
mph-awad at hotmail.com
Wed Sep 14 11:40:40 CEST 2016
Dears Gate users,
I am running an example as exercise for basic dose calculations. It consists of a cylindrical phantom (patient) contains a target treated by ONE proton PBS. Then I summed all the resulted dose in phantom and target separately (phantom 18.81 cGy, target 48.3 cGy). I run the example again with TWO beams and I got the same numbers for both dose actors. The number of primaries was the same in both simulations.
Is the primaries are split on the two beams by default? Am I on the correct way of dose calc.?
Here a copy of the main macro.
Thanks in advance,
Awad
#=====================================================
# 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
# WATERCYLINDER
/gate/world/daughters/name water_cylinder
/gate/world/daughters/insert cylinder
/gate/water_cylinder/setMaterial Water
/gate/water_cylinder/geometry/setRmin 0 cm
/gate/water_cylinder/geometry/setRmax 15 cm
/gate/water_cylinder/geometry/setHeight 40 cm
/gate/water_cylinder/placement/setRotationAxis 1 0 0
/gate/water_cylinder/placement/setRotationAngle 90 deg
/gate/water_cylinder/placement/setTranslation 0.0 0.0 0.0 cm
/gate/water_cylinder/vis/setColor blue
/gate/water_cylinder/vis/forceWireframe
#Tumor
/gate/water_cylinder/daughters/name target
/gate/water_cylinder/daughters/insert sphere
/gate/target/setMaterial Water
/gate/target/geometry/setRmin 0 cm
/gate/target/geometry/setRmax 3 cm
/gate/target/placement/setTranslation 0.0 0.0 0.0 cm
/gate/target/vis/setColor red
#=====================================================
# 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_cylinder 0.1 mm
/gate/physics/Electron/SetCutInRegion water_cylinder 0.1 mm
/gate/physics/Positron/SetCutInRegion water_cylinder 0.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 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/phantom2Fields.txt
/gate/actor/dose/attachTo water_cylinder
/gate/actor/dose/stepHitType random
/gate/actor/dose/setPosition 0 0 0 cm
/gate/actor/dose/setResolution 100 100 300
/gate/actor/dose/setSize 8 8 30 cm
/gate/actor/dose/saveEveryNSeconds 60
/gate/actor/dose/enableEdep false
/gate/actor/dose/enableUncertaintyEdep false
/gate/actor/dose/enableDose true
/gate/actor/dose/enableUncertaintyDose true
/gate/actor/dose/enableNumberOfHits false
/gate/actor/addActor DoseActor dose1
/gate/actor/dose1/save output/target2Fields.txt
/gate/actor/dose1/attachTo target
/gate/actor/dose1/stepHitType random
/gate/actor/dose1/setPosition 0 0 0 cm
/gate/actor/dose1/setResolution 100 100 100
/gate/actor/dose1/setSize 6 6 6 cm
/gate/actor/dose1/saveEveryNSeconds 60
/gate/actor/dose1/enableEdep false
/gate/actor/dose1/enableUncertaintyEdep false
/gate/actor/dose1/enableDose true
/gate/actor/dose1/enableUncertaintyDose true
/gate/actor/dose1/enableNumberOfHits false
#=====================================
# SOURCE
#=====================================
/gate/source/addSource PBS PencilBeam
/gate/source/PBS/setPosition 0 0 -50 cm
/gate/source/PBS/setEnergy 150 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 5 mm*mrad
/gate/source/PBS/setEllipseXThetaRotationNorm positive
/gate/source/PBS/setEllipseYPhiEmittance 5 mm*mrad
/gate/source/PBS/setEllipseYPhiRotationNorm positive
/gate/source/PBS/setParticleType proton
/gate/source/PBS/setRotationAxis 1 0 0
/gate/source/PBS/setRotationAngle 0 deg
/gate/source/addSource PBS2 PencilBeam
/gate/source/PBS2/setPosition 0 0 50 cm
/gate/source/PBS2/setEnergy 150 MeV
/gate/source/PBS2/setSigmaEnergy 0.1 MeV
/gate/source/PBS2/setSigmaX 3 mm
/gate/source/PBS2/setSigmaY 3 mm
/gate/source/PBS2/setSigmaTheta 1 mrad
/gate/source/PBS2/setSigmaPhi 1 mrad
/gate/source/PBS2/setEllipseXThetaEmittance 5 mm*mrad
/gate/source/PBS2/setEllipseXThetaRotationNorm positive
/gate/source/PBS2/setEllipseYPhiEmittance 5 mm*mrad
/gate/source/PBS2/setEllipseYPhiRotationNorm positive
/gate/source/PBS2/setParticleType proton
/gate/source/PBS2/setRotationAxis 1 0 0
/gate/source/PBS2/setRotationAngle 180 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/20160914/ba1a0f66/attachment-0001.html>
More information about the Gate-users
mailing list