[Gate-users] Gat exercices - photon beam, dose in water - only 0 values as output

Matthew Strugari matthew.strugari at dal.ca
Mon Dec 13 14:40:13 CET 2021


Did you increase your number of primaries from 1?

Matthew

On Dec 13, 2021, at 9:00 AM, Sarah Blind <sarah.blind at univ-lorraine.fr> wrote:


CAUTION: The Sender of this email is not from within Dalhousie.
I still have the same problem even without   /gate/application/noGlobalOutput ...

----- Le 13 Déc 21, à 13:43, Matthew Strugari <matthew.strugari at dal.ca> a écrit :
Hi Sarah,

You have the output disabled. Remove /gate/application/noGlobalOutput from your macro and try again.

Best,
Matthew

On Dec 13, 2021, at 6:19 AM, Sarah Blind <sarah.blind at univ-lorraine.fr> wrote:

CAUTION: The Sender of this email is not from within Dalhousie.
Dear Gate users,

I am trying to follow the exercices (1 and 2) of https://davidsarrut.pages.in2p3.fr/gate-exercices-site/docs/welcome/.
When I open the output files (with ImageJ) I obtain only 0 values.
Does anyone have any idea why those dose distribution values are all zero?


# Visualization
/vis/open OGLIQt
/vis/drawVolume
/vis/viewer/flush
/tracking/storeTrajectory             1
/vis/scene/add/trajectories
/vis/scene/endOfEventAction           accumulate
/vis/modeling/trajectories/create/drawByParticleID

/gate/geometry/setMaterialDatabase ../data/GateMaterials.db

# World
/gate/world/geometry/setXLength 1 m
/gate/world/geometry/setYLength 1 m
/gate/world/geometry/setZLength 1 m
/gate/world/setMaterial Air

# Water Box
/gate/world/daughters/name              waterbox
/gate/world/daughters/insert            box
/gate/waterbox/geometry/setXLength      40 cm
/gate/waterbox/geometry/setYLength      40 cm
/gate/waterbox/geometry/setZLength      40 cm
/gate/waterbox/placement/setTranslation 0 0 25 cm
/gate/waterbox/setMaterial Water
/gate/waterbox/vis/setVisible 1
/gate/waterbox/vis/setColor blue

# aluminium box
/gate/world/daughters/name     alubox
/gate/world/daughters/insert    box
/gate/alubox/geometry/setXLength      5 cm
/gate/alubox/geometry/setYLength      5 cm
/gate/alubox/geometry/setZLength      5 cm
/gate/alubox/placement/setTranslation 0 0 15 cm
/gate/alubox/setMaterial Aluminium
/gate/alubox/vis/setVisible 1
/gate/alubox/vis/setColor yellow

#=========
# PHYSICS
#=========

/gate/physics/addPhysicsList emstandard_opt3

/gate/physics/Gamma/SetCutInRegion      world 10 mm
/gate/physics/Electron/SetCutInRegion   world 10 mm
/gate/physics/Positron/SetCutInRegion   world 10 mm

/gate/physics/Gamma/SetCutInRegion      waterbox 0.01 mm
/gate/physics/Electron/SetCutInRegion   waterbox 0.01 mm
/gate/physics/Positron/SetCutInRegion   waterbox 0.01 mm

/gate/physics/displayCuts


#===========
# DETECTORS
#===========

/gate/actor/addActor                        DoseActor  depthdose
/gate/actor/depthdose/save                  output/gamma-depth.txt
/gate/actor/depthdose/attachTo            waterbox
/gate/actor/depthdose/stepHitType           random
/gate/actor/depthdose/setPosition           0 0 0 cm
/gate/actor/depthdose/setResolution         1 1 100
/gate/actor/depthdose/saveEveryNSeconds     60
/gate/actor/depthdose/enableEdep            true
/gate/actor/depthdose/enableUncertaintyEdep true
/gate/actor/depthdose/enableDose            true
/gate/actor/depthdose/enableUncertaintyDose true
/gate/actor/depthdose/enableNumberOfHits    true

/gate/actor/addActor                          DoseActor  doseprofile
/gate/actor/doseprofile/save                  output/gamma-profile.txt
/gate/actor/doseprofile/attachTo              waterbox
/gate/actor/doseprofile/stepHitType           random
/gate/actor/doseprofile/setPosition           0 0 4 cm
/gate/actor/doseprofile/setResolution         100 1 1
/gate/actor/doseprofile/setSize               5 40 40 cm
/gate/actor/doseprofile/saveEveryNSeconds     60
/gate/actor/doseprofile/enableEdep            true
/gate/actor/doseprofile/enableUncertaintyEdep true
/gate/actor/doseprofile/enableDose            true
/gate/actor/doseprofile/enableUncertaintyDose true
/gate/actor/doseprofile/enableNumberOfHits    true

/gate/actor/addActor               SimulationStatisticActor stat
/gate/actor/stat/save              output/gamma-stat.txt
/gate/actor/stat/saveEveryNSeconds 60

/gate/actor/addActor                                              DoseActor  doseDistribution3D_1
/gate/actor/doseDistribution3D_1/save                               output/doseDistribution3D_1.hdr
/gate/actor/doseDistribution3D_1/attachTo                           waterbox
/gate/actor/doseDistribution3D_1/stepHitType                        random
/gate/actor/doseDistribution3D_1/setPosition                        0.0 0.0 0.0 um
/gate/actor/doseDistribution3D_1/setSize                    40.0 40.0 40.0 um
/gate/actor/doseDistribution3D_1/setResolution                      200 200 200
/gate/actor/doseDistribution3D_1/enableDose                         true
/gate/actor/doseDistribution3D_1/saveEveryNSeconds 60


#=============
# INITIALISATION
#==============

/gate/run/initialize

#================
# BEAMS / SOURCE
#================

/gate/source/addSource mybeam gps
/gate/source/mybeam/gps/particle gamma
/gate/source/mybeam/gps/pos/type Beam
/gate/source/mybeam/gps/pos/rot1 1 1 1
#/gate/source/mybeam/gps/pos/rot2 0 0 0
/gate/source/mybeam/gps/pos/shape Circle
/gate/source/mybeam/gps/pos/centre 0 0 0 mm
/gate/source/mybeam/gps/pos/sigma_x 3 mm
/gate/source/mybeam/gps/pos/sigma_y 3 mm
/gate/source/mybeam/gps/ene/mono 10 MeV
/gate/source/mybeam/gps/ene/type Gauss
/gate/source/mybeam/gps/ene/sigma 1.0 MeV
/gate/source/mybeam/gps/direction 0 0 1


#=============
# START BEAMS
#=============

/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto

/gate/application/noGlobalOutput
/gate/application/setTotalNumberOfPrimaries 1
/gate/application/start



Many thanks for your help,

Sarah Blind
_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users


--
Sarah Blind
Ingénieure de recherche - CNRS
Centre de Recherche en Automatique de Nancy (CRAN)
Département Biologie, Signaux et Systèmes en Cancérologie et Neurosciences (BioSiS)
Site de la faculté de médecine
9, Avenue de la Forêt de Haye
54505 Vandoeuvre les Nancy
Tél : 03.72.74.61.13
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20211213/566906ef/attachment-0001.html>


More information about the Gate-users mailing list