[Gate-users] Fw: O15 production is not enough for SoftTissue.
hajar zarei
hajar_genetic2002 at yahoo.com
Fri Apr 3 05:45:10 CEST 2020
Hello everyone,
I would like to calculate the O15 production from my following input. I'm using Gate8.2I have found the output and tried to compare my output with some documents. I found my results per particle is very poor like 10 times less than from the documents.
Please help me providing the correct actor for production O15, C11 .....
Mymain mac file ....#=====================================================
# GEOMETRY
#=====================================================
/gate/geometry/setMaterialDatabase data/GateMaterials.db
# WORLD(5x5x5m Vacuum box)
/gate/world/setMaterial Air
/gate/world/geometry/setXLength 5 m
/gate/world/geometry/setYLength 5 m
/gate/world/geometry/setZLength 5 m
/gate/world/vis/setVisible 0
# MAIN WATER BOX (30x30x30 cm in blue)
/gate/world/daughters/name waterbox
/gate/world/daughters/insert box
/gate/waterbox/geometry/setXLength 30. cm
/gate/waterbox/geometry/setYLength 30. cm
/gate/waterbox/geometry/setZLength 40. cm
/gate/waterbox/placement/setTranslation 0. 0. 0. cm
/gate/waterbox/setMaterial SoftTissue
/gate/waterbox/vis/forceSolid
/gate/waterbox/vis/setVisible 1
/gate/waterbox/vis/setColor grey
#=====================================================
# PHYSICS
#=====================================================
/gate/physics/addPhysicsList QGSP_BERT_HP_EMY
/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 waterbox 0.1 mm
/gate/physics/Electron/SetCutInRegion waterbox 0.1 mm
/gate/physics/Positron/SetCutInRegion waterbox 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-
# Verbose physic list (uncomment to display)
/gate/physics/displayCuts
# /gate/physics/processList Enabled
# /gate/physics/processList Initialized
#=====================================================
# DOSE ACTOR
#=====================================================
/gate/actor/addActor DoseActor doseDistribution
/gate/actor/doseDistribution/save pencil-output.txt
/gate/actor/doseDistribution/attachTo waterbox
/gate/actor/doseDistribution/stepHitType random
/gate/actor/doseDistribution/setPosition 0 0 0 mm
/gate/actor/doseDistribution/setSize 300 300 400 mm
/gate/actor/doseDistribution/setResolution 1 1 400
/gate/actor/doseDistribution/enableDose true
/gate/actor/doseDistribution/enableSquaredDose false
/gate/actor/doseDistribution/enableUncertaintyDose false
/gate/actor/doseDistribution/enableNumberOfHits false
/gate/actor/doseDistribution/enableEdep false
/gate/actor/doseDistribution/enableUncertaintyEdep false
/gate/actor/doseDistribution/saveEveryNSeconds 300
#=====================================================
# Production ACTOR
#=====================================================
/gate/actor/addActor ProductionActor prodO15
/gate/actor/prodO15/save proton-O15.txt
/gate/actor/prodO15/attachTo waterbox
/gate/actor/prodO15/stepHitType random
/gate/actor/prodO15/setPosition 0 0 0 mm
#/gate/actor/prodO15/setVoxelSize 1 1 1 mm
#/gate/actor/prodO15/setResolution 1 1 400
/gate/actor/prodO15/setSize 1 1 400 mm
/gate/actor/prodO15/setResolution 1 1 400
/gate/actor/prodO15/saveEveryNSeconds 20
/gate/actor/prodO15/stepHitType post
/gate/actor/prodO15/addFilter particleFilter
/gate/actor/prodO15/particleFilter/addParticle O15
#=====================================================
# INITIALISATION
#=====================================================
/gate/run/initialize
#=====================================================
# BEAMS
#=====================================================
/gate/source/addSource PBS PencilBeam
/gate/source/PBS/setParticleType proton
/gate/source/PBS/setEnergy 150.0 MeV
/gate/source/PBS/setSigmaEnergy 1.0 MeV
/gate/source/PBS/setPosition 0 0 -500 mm
/gate/source/PBS/setSigmaX 2 mm
/gate/source/PBS/setSigmaY 4 mm
/gate/source/PBS/setSigmaTheta 3.3 mrad
/gate/source/PBS/setSigmaPhi 3.8 mrad
/gate/source/PBS/setEllipseXThetaEmittance 15 mm*mrad
/gate/source/PBS/setEllipseXThetaRotationNorm negative
/gate/source/PBS/setEllipseYPhiEmittance 20 mm*mrad
/gate/source/PBS/setEllipseYPhiRotationNorm negative
/gate/source/PBS/setRotationAxis 0 0 1
/gate/source/PBS/setRotationAngle 180 deg
#=====================================================
# START BEAMS
#=====================================================
/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto
/gate/application/noGlobalOutput
/gate/application/setTotalNumberOfPrimaries 40000
/gate/application/start
--
With Regards,
M Rafiqul IslamM.Sc.(Physics, R.U ), Advanced MS. (Medical Physics, Italy)MEXT Fellow(Japan), IAEA Fellow. PhD Student in Medical Physics
Tohoku University, JAPAN.#Senior Scientific Officer
#BANGLADESH ATOMIC ENERGY COMMISSIONCell: +81-080-3333-1066.
--
With Regards,
M Rafiqul IslamM.Sc.(Physics, R.U ), Advanced MS. (Medical Physics, Italy)MEXT Fellow(Japan), IAEA Fellow. PhD Student in Medical Physics
Tohoku University, JAPAN.#Senior Scientific Officer
#BANGLADESH ATOMIC ENERGY COMMISSIONCell: +81-080-3333-1066.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200403/a72c6e04/attachment-0001.html>
More information about the Gate-users
mailing list