[Gate-users] Why Y90 has different energy distributions when using ion source and fastY90 source?

[Jinxia Guo]

Hi guys,



Recently I’ve done some simulation with Y90 source in Gate7.2 develop
version. Two settings are with the same camera, digitizer and output but
respectively with normal ion source and fastY90 source.  The energy
distributions of singles in coincidences for each setting are shown below. It
seems that the energy distribution with the normal ion source is very
different with the fastY90 source. I thought they should be the same, but
why the results actually are not? Anything wrong in the simulations
setting? Thanks!


*The energy distributions of each setting.*

(1)   fastY90 source (energy window: [420, inf])                (2)   ion
source (energy window:[420, inf])




       (3) fastY90 source (energy window: [0, inf])
  (4)   ion source (energy window:[0, inf])




（1）Setting1: simulation with fastY90

#-----------------Phantom--------------------------------

/gate/world/daughters/name cylinder_phantom

/gate/world/daughters/insert cylinder

/gate/cylinder_phantom/setMaterial Water

/gate/cylinder_phantom/geometry/setRmax   10  cm

/gate/cylinder_phantom/geometry/setRmin   0  cm

/gate/cylinder_phantom/geometry/setHeight 30.5  cm

/gate/cylinder_phantom/placement/setTranslation 0. 0. 0. mm

/gate/cylinder_phantom/attachPhantomSD



#-----------------Source--------------------------------

/gate/source/addSource ionsrc fastY90

/gate/source/ionsrc/setActivity 0.1 Ci

/gate/source/ionsrc/useDefaultHalfLife

/gate/source/ionsrc/gps/angtype iso

/gate/source/ionsrc/gps/number 1



/gate/source/ionsrc/setMinBremEnergy 50 keV

/gate/source/ionsrc/setPositronProbability 0.00003186



/gate/source/ionsrc/gps/type Volume

/gate/source/ionsrc/gps/shape Cylinder

/gate/source/ionsrc/gps/radius  10 cm

/gate/source/ionsrc/gps/halfz  15.0325 cm

/gate/source/ionsrc/gps/centre 0. 0. 0. Mm

/gate/source/ionsrc/gps/confine cylinder_phantom



#-----------------Physics--------------------------------

/gate/physics/addProcess Decay

/gate/physics/addProcess RadioactiveDecay



/gate/physics/addProcess PhotoElectric

/gate/physics/processes/PhotoElectric/setModel StandardModel



/gate/physics/addProcess Compton

/gate/physics/processes/Compton/setModel StandardModel



/gate/physics/addProcess RayleighScattering

/gate/physics/processes/RayleighScattering/setModel PenelopeModel



/gate/physics/addProcess ElectronIonisation

/gate/physics/processes/ElectronIonisation/setModel StandardModel e-

/gate/physics/processes/ElectronIonisation/setModel StandardModel e+



/gate/physics/addProcess Bremsstrahlung

/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-

/gate/physics/processes/Bremsstrahlung/setModel StandardModel e+



/gate/physics/addProcess PositronAnnihilation

/gate/physics/addProcess MultipleScattering e+

/gate/physics/addProcess MultipleScattering e-



/gate/physics/setEMin 0.1 keV

/gate/physics/setEMax 2.0 GeV

/gate/physics/setDEDXBinning 220

/gate/physics/setLambdaBinning 220



/gate/physics/processList Enabled

/gate/physics/processList Initialized



（2）Setting2: simulation with the ion source

#-----------------Phantom--------------------------------

/gate/world/daughters/name cylinder_phantom

/gate/world/daughters/insert cylinder

/gate/cylinder_phantom/setMaterial Water

/gate/cylinder_phantom/geometry/setRmax   10  cm

/gate/cylinder_phantom/geometry/setRmin   0  cm

/gate/cylinder_phantom/geometry/setHeight 30.065  cm

/gate/cylinder_phantom/placement/setTranslation 0. 0. 0. mm

/gate/cylinder_phantom/attachPhantomSD



#-----------------Source--------------------------------

/gate/source/addSource ionsrc

/gate/source/ionsrc/setActivity 0.10 Ci

/gate/source/ionsrc/gps/particle ion

/gate/source/ionsrc/gps/ion 39 90 0 0

/gate/source/ionsrc/setForcedUnstableFlag true

#/gate/source/ionsrc/setForcedHalfLife 230544 s

/gate/source/ionsrc/useDefaultHalfLife

/gate/source/ionsrc/gps/angtype iso

/gate/source/ionsrc/gps/number 1



/gate/source/ionsrc/gps/type Volume

/gate/source/ionsrc/gps/shape Cylinder

/gate/source/ionsrc/gps/radius  10 cm

/gate/source/ionsrc/gps/halfz  15.0325 cm

/gate/source/ionsrc/gps/centre 0. 0. 0. mm



/gate/source/ionsrc/gps/confine cylinder_phantom



#-----------------Physics--------------------------------

/gate/physics/addProcess Decay

/gate/physics/addProcess RadioactiveDecay



/gate/physics/addProcess PhotoElectric

/gate/physics/processes/PhotoElectric/setModel StandardModel



/gate/physics/addProcess Compton

/gate/physics/processes/Compton/setModel StandardModel



/gate/physics/addProcess RayleighScattering

/gate/physics/processes/RayleighScattering/setModel PenelopeModel



/gate/physics/addProcess GammaConversion

/gate/physics/processes/GammaConversion/setModel StandardModel



/gate/physics/addProcess ElectronIonisation

/gate/physics/processes/ElectronIonisation/setModel StandardModel e-

/gate/physics/processes/ElectronIonisation/setModel StandardModel e+



/gate/physics/addProcess Bremsstrahlung

/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-

/gate/physics/processes/Bremsstrahlung/setModel StandardModel e+



/gate/physics/addProcess PositronAnnihilation

/gate/physics/addProcess MultipleScattering e+

/gate/physics/addProcess MultipleScattering e-



/gate/physics/setEMin 0.1 keV

/gate/physics/setEMax 2.0 GeV

/gate/physics/setDEDXBinning 220

/gate/physics/setLambdaBinning 220



/gate/physics/processList Enabled

/gate/physics/processList Initialized





Regards,

Jinxia
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20170123/64ade27b/attachment-0001.html>