[Gate-users] Fwd: All Gate users, d

[Subin Kim]

Hi, all Gate users,

I simulated mouse whole body dosimetry with
68Ga-labeled radiotracer for 10^8 number of events. and I got too many NA
values(or value=0) in the 3D output ( (i.e., Dose and Edep).

I discussed various causes affecting these results with my co-worker.
I guessed there is the wrong code in "geometry step", because the black
dots  (NA values or value=0) are dispersed on the air materials.
But i couldn't find the exact ways to solve these problems.

 The examples of the black dots are below.
[image: image.png]


My questions are obvious:

1. Why did i get these NA values in every 3D output?
2. How can i solve this problem?


#=====================================================
# VEROBOSITY
#=====================================================

/control/execute  mac/verbose.mac

#=====================================================
# VISUALISATION
#=====================================================

/control/execute mac/visu.mac

#=====================================================
# GEOMETRY
#=====================================================

/gate/geometry/setMaterialDatabase data/GateMaterials.db

# World
/gate/world/geometry/setXLength 10 cm
/gate/world/geometry/setYLength 10 cm
/gate/world/geometry/setZLength 10 cm
/gate/world/setMaterial Air
/gate/world/vis/setVisible 1
/gate/world/vis/forceWireframe

# Voxel phantom
/gate/world/daughters/name                      ct
/gate/world/daughters/insert
 ImageRegularParametrisedVolume
/gate/ct/geometry/setRangeToMaterialFile     data/range_geo.dat
/gate/ct/geometry/setImage                   data/ct-gate.mhd
/gate/ct/setSkipEqualMaterials 1
/gate/ct/vis/setVisible 1
# /gate/ct/vis/setColor blue
/gate/ct/vis/forceWireframe

# /gate/geometry/rebuild

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

/gate/physics/addPhysicsList emstandard_opt3
/gate/physics/addProcess RadioactiveDecay #to generate decay process of F-18

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

/gate/actor/addActor                                DoseActor
 doseDistribution
/gate/actor/doseDistribution/attachTo       ct
/gate/actor/doseDistribution/stepHitType            random
/gate/actor/doseDistribution/setSize                30 28.8 99 mm #Size in
mm
/gate/actor/doseDistribution/setPosition            0 0 0 cm # Origin
/gate/actor/doseDistribution/setResolution          58 48 165 #Size (pixel)
/gate/actor/doseDistribution/enableEdep             true
/gate/actor/doseDistribution/enableUncertaintyEdep  true
/gate/actor/doseDistribution/enableSquaredEdep      false
/gate/actor/doseDistribution/enableDose             true #from true
/gate/actor/doseDistribution/setDoseAlgorithm VolumeWeighting
# /gate/actor/doseDistribution/enableDoseToWater      true #new set
/gate/actor/doseDistribution/enableUncertaintyDose  true #from true
/gate/actor/doseDistribution/enableSquaredDose      false
/gate/actor/doseDistribution/enableNumberOfHits     true
/gate/actor/doseDistribution/save                   output/output.mhd

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

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

/gate/run/initialize

#=====================================================
# SOURCE :
#=====================================================

#Description

#Ion
/gate/source/addSource ion voxel
/gate/source/ion/reader/insert image
/gate/source/ion/imageReader/translator/insert linear
/gate/source/ion/imageReader/linearTranslator/setScale 1. Bq
/gate/source/ion/imageReader/readFile data/pt-gate.mhd
/gate/source/ion/imageReader/verbose 1
/gate/source/ion/setPosition -15.0 -14.4 -49.5 mm

/gate/source/ion/gps/particle ion
/gate/source/ion/gps/ion 9 18 0 0 #F-18
/gate/source/ion/gps/energytype Mono
/gate/source/ion/gps/monoenergy 0. keV
/gate/source/ion/gps/angtype iso
/gate/source/ion/setForcedUnstableFlag true
/gate/source/ion/useDefaultHalfLife 6586 s


/gate/source/ion/dump                                   1

# /gate/source/list

#=====================================================
# RANDOM
#=====================================================
/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto


#=====================================================
# START
#=====================================================
/gate/application/noGlobalOutput
/gate/application/setTotalNumberOfPrimaries  100000000
/gate/application/setTimeStart 0s
/gate/application/setTimeSlice 1s
/gate/application/setTimeStop 60s
/gate/application/start




Thank you for your help

Best regards.

Subin Kim
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200904/8c78c053/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 148837 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200904/8c78c053/attachment-0001.png>