<div dir="ltr"><div>Hello Gaters,<br><br></div><div>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.<br><br>Hope you can help me understand the issue.<br><br>Here's the code:<br><br></div><div>ION SOURCE:<br></div><div>#=====================================================<br># GEOMETRY<br>#=====================================================<br><br>/gate/geometry/setMaterialDatabase data/GateMaterials.db<br><br>#World<br><br>/gate/world/setMaterial                          Air<br>/gate/world/geometry/setXLength         1 m<br>/gate/world/geometry/setYLength         1 m<br>/gate/world/geometry/setZLength         1 m<br><br><br>#Waterbox<br><br>/gate/world/daughters/name                       waterbox<br>/gate/world/daughters/insert                       box<br>/gate/waterbox/geometry/setXLength        30 cm<br>/gate/waterbox/geometry/setYLength        30 cm<br>/gate/waterbox/geometry/setZLength        30 cm<br>/gate/waterbox/setMaterial                         Water<br><br>/gate/waterbox/vis/setVisible             1<br>/gate/waterbox/vis/setColor             blue<br><br>#=====================================================<br># PHYSICS<br>#=====================================================<br><br>/gate/physics/addPhysicsList             LBE<br>/gate/physics/addPhysicsList             emlivermore<br>/gate/physics/addProcess                  RadioactiveDecay<br><br>/gate/physics/Gamma/SetCutInRegion      waterbox  1 km<br>/gate/physics/Electron/SetCutInRegion      waterbox  1 km<br>/gate/physics/Positron/SetCutInRegion      waterbox  1 km<br><br>/gate/physics/displayCuts<br><br>#=====================================================<br># DETECTORS<br>#=====================================================<br><br>/gate/actor/addActor                                                       EnergySpectrumActor    ESA<br>/gate/actor/ESA/attachTo                                                waterbox<br>/gate/actor/ESA/save                                                      output/energy-spectrum1.root<br>/gate/actor/ESA/energySpectrum/setEmin                     1 keV             <br>/gate/actor/ESA/energySpectrum/setEmax                    2.5 MeV<br>/gate/actor/ESA/energySpectrum/setNumberOfBins     100<br><br>/gate/actor/addActor                                          DoseActor dosedepth                                             <br>/gate/actor/dosedepth/attachTo                        waterbox<br>/gate/actor/dosedepth/stepHitType                    random<br>/gate/actor/dosedepth/setPosition                     0 0 0 cm<br>/gate/actor/dosedepth/setResolution                 1 1 100<br>/gate/actor/dosedepth/setSize                           30 30 5 cm<br>/gate/actor/dosedepth/saveEveryNSeconds      20<br>/gate/actor/dosedepth/enableEdep                    false<br>/gate/actor/dosedepth/enableUncertaintyEdep false<br>/gate/actor/dosedepth/enableDose                    true<br>/gate/actor/dosedepth/enableUncertaintyDose true<br>/gate/actor/dosedepth/enableNumberOfHits     false<br>/gate/actor/dosedepth/save                               output/dosedepth1.txt<br><br>/gate/actor/addActor               SimulationStatisticActor stat<br>/gate/actor/stat/save              output/gamma-stat1.txt<br>/gate/actor/stat/saveEveryNSeconds 20<br><br>#=====================================================<br># INITIALISATION<br>#=====================================================<br><br>/gate/run/initialize<br><br>#=====================================================<br># SOURCE<br>#=====================================================<br><br> /gate/source/addSource                                     Yttrium <br> /gate/source/Yttrium/setActivity                        1.3 MBq <br> /gate/source/Yttrium/gps/particle                       ion<br> /gate/source/Yttrium/gps/ion                              39 90 0 0<br> /gate/source/Yttrium/setForcedUnstableFlag     true<br> /gate/source/Yttrium/useDefaultHalfLife<br> /gate/source/Yttrium/gps/type                            Volume<br> /gate/source/Yttrium/gps/shape                         Sphere<br> /gate/source/Yttrium/gps/radius                         1.5 mm<br> /gate/source/Yttrium/gps/centre                         0. 0. 0. mm<br> /gate/source/Yttrium/gps/energy                        0 MeV<br><br> /gate/source/Yttrium/visualize 1000 red 1 <br><br>#=====================================================<br># START BEAM<br>#=====================================================<br><br> /gate/random/setEngineName     MersenneTwister<br> /gate/random/setEngineSeed     auto<br><br> /gate/application/noGlobalOutput<br><br> /gate/application/setTimeStart 0 s<br> /gate/application/setTimeSlice 1 s<br> /gate/application/setTimeStop  1 s<br> <br> /gate/application/start<br>#======================================================<br><br></div><div>fastY90 SOURCE (IDENTICAL EXCEPT FOR THE SOURCE PART):<br><br> /gate/source/addSource                              Yttrium fastY90<br> /gate/source/Yttrium/setActivity                1.3 MBq<br> /gate/source/Yttrium/setMinBremEnergy   50 keV<br><br> /gate/source/Yttrium/visualize 100 red 1 <br><br></div><div>Best wishes,<br><br></div><div>Nicola Rizzi<br></div></div>