[Gate-users] Fw: O15 production is not enough for SoftTissue.

[David Boersma]

Hi!


Could you maybe replace the "QGSP_BERT_HP_EMY" physics list by "Shielding_EMY" and check if you get better matching rates for those isotopes?


HTH,

David

________________________________
Von: Gate-users <gate-users-bounces at lists.opengatecollaboration.org> im Auftrag von hajar zarei <hajar_genetic2002 at yahoo.com>
Gesendet: Freitag, 3. April 2020 05:45:10
An: gate-users at lists.opengatecollaboration.org
Betreff: [Gate-users] Fw: O15 production is not enough for SoftTissue.



Hello everyone,

I would like to calculate the O15 production from my following input. I'm using Gate8.2
I 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 Islam
M.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 COMMISSION
Cell: +81-080-3333-1066.


--

With Regards,

M Rafiqul Islam
M.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 COMMISSION
Cell: +81-080-3333-1066.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200405/d05b160a/attachment-0001.html>