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

[Matthew Strugari]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20211213/9505ab7a/attachment-0001.html>