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

[Sarah Blind]

Dear Gate users, 

I am trying to follow the exercices (1 and 2) of [ https://davidsarrut.pages.in2p3.fr/gate-exercices-site/docs/welcome/ | 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 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20211213/98a7e566/attachment.html>