[Gate-users] FastY90 spectrum and Ion Source spectrum don't match

[Nicola Rizzi]

Hello Gaters,

I'm working on GATE for my Bachelor final project and I have a problem on
choosing the right source of Y90 for a SPECT simulation. The problem is
that I obtained two different spectra for simulations running,
respectively, an Ion source and a fast Y90 source. I found an old
submission on this topic but I was not able to solve my problem anyway.
Moreover, I cannot find a complete documentation about fastY90 source,
except the few lines in the User Guide.

Hope you can help me understand the issue.

Here's the code:

ION SOURCE:
#=====================================================
# GEOMETRY
#=====================================================

/gate/geometry/setMaterialDatabase data/GateMaterials.db

#World

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


#Waterbox

/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        30 cm
/gate/waterbox/setMaterial                         Water

/gate/waterbox/vis/setVisible             1
/gate/waterbox/vis/setColor             blue

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

/gate/physics/addPhysicsList             LBE
/gate/physics/addPhysicsList             emlivermore
/gate/physics/addProcess                  RadioactiveDecay

/gate/physics/Gamma/SetCutInRegion      waterbox  1 km
/gate/physics/Electron/SetCutInRegion      waterbox  1 km
/gate/physics/Positron/SetCutInRegion      waterbox  1 km

/gate/physics/displayCuts

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

/gate/actor/addActor
EnergySpectrumActor    ESA
/gate/actor/ESA/attachTo
waterbox
/gate/actor/ESA/save
output/energy-spectrum1.root
/gate/actor/ESA/energySpectrum/setEmin                     1 keV

/gate/actor/ESA/energySpectrum/setEmax                    2.5 MeV
/gate/actor/ESA/energySpectrum/setNumberOfBins     100

/gate/actor/addActor                                          DoseActor
dosedepth
/gate/actor/dosedepth/attachTo                        waterbox
/gate/actor/dosedepth/stepHitType                    random
/gate/actor/dosedepth/setPosition                     0 0 0 cm
/gate/actor/dosedepth/setResolution                 1 1 100
/gate/actor/dosedepth/setSize                           30 30 5 cm
/gate/actor/dosedepth/saveEveryNSeconds      20
/gate/actor/dosedepth/enableEdep                    false
/gate/actor/dosedepth/enableUncertaintyEdep false
/gate/actor/dosedepth/enableDose                    true
/gate/actor/dosedepth/enableUncertaintyDose true
/gate/actor/dosedepth/enableNumberOfHits     false
/gate/actor/dosedepth/save
output/dosedepth1.txt

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

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

/gate/run/initialize

#=====================================================
# SOURCE
#=====================================================

 /gate/source/addSource                                     Yttrium
 /gate/source/Yttrium/setActivity                        1.3 MBq
 /gate/source/Yttrium/gps/particle                       ion
 /gate/source/Yttrium/gps/ion                              39 90 0 0
 /gate/source/Yttrium/setForcedUnstableFlag     true
 /gate/source/Yttrium/useDefaultHalfLife
 /gate/source/Yttrium/gps/type                            Volume
 /gate/source/Yttrium/gps/shape                         Sphere
 /gate/source/Yttrium/gps/radius                         1.5 mm
 /gate/source/Yttrium/gps/centre                         0. 0. 0. mm
 /gate/source/Yttrium/gps/energy                        0 MeV

 /gate/source/Yttrium/visualize 1000 red 1

#=====================================================
# START BEAM
#=====================================================

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

 /gate/application/noGlobalOutput

 /gate/application/setTimeStart 0 s
 /gate/application/setTimeSlice 1 s
 /gate/application/setTimeStop  1 s

 /gate/application/start
#======================================================

fastY90 SOURCE (IDENTICAL EXCEPT FOR THE SOURCE PART):

 /gate/source/addSource                              Yttrium fastY90
 /gate/source/Yttrium/setActivity                1.3 MBq
 /gate/source/Yttrium/setMinBremEnergy   50 keV

 /gate/source/Yttrium/visualize 100 red 1

Best wishes,

Nicola Rizzi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20170718/abd382df/attachment.html>