[Gate-users] no scatter with electron source

[Katrin Teichert]

Hello Gaters,

I have the following problem with my SPECT-camera: if I use an
electron source (1 MeV) for bremsstrahl-imaging, I don't record any
Compton scattering, and if I use an ion source with Yttrium 90, I only
record very little scatter, whereas if I use a gamma source (140 keV),
most particles are scattered.
I have also tried with a positron source, and here again, most particles
are scattered.
Can somebody help me with this problem? My source code is below.
Yours, Katrin

#/control/execute /home/teicherk/Gate/visual.mac
/vis/disable
/gate/geometry/enableAutoUpdate


#WORLD

/gate/world/geometry/setXLength 1. m
/gate/world/geometry/setYLength 1. m
/gate/world/geometry/setZLength 1. m


#KAMERA

/gate/world/daughters/name SPECThead
/gate/world/daughters/insert box

/gate/SPECThead/vis/forceWireframe
/gate/SPECThead/vis/setColor blue

/gate/SPECThead/setMaterial Air
/gate/SPECThead/geometry/setXLength 25 cm
/gate/SPECThead/geometry/setYLength 25 cm
/gate/SPECThead/geometry/setZLength 14.6 cm
/gate/SPECThead/placement/setTranslation 0. 0. -20 cm
/gate/SPECThead/attachPhantomSD


# Shielding

/gate/SPECThead/daughters/name shielding
/gate/SPECThead/daughters/insert box
/gate/shielding/geometry/setXLength 25. cm
/gate/shielding/geometry/setYLength 25. cm
/gate/shielding/geometry/setZLength 14.6 cm
/gate/shielding/placement/setTranslation  0. 0. 0. cm
/gate/shielding/setMaterial Lead
/gate/shielding/vis/setColor red
/gate/shielding/vis/forceWireframe
/gate/shielding/attachPhantomSD


#HEGP

/gate/SPECThead/daughters/name collimator
/gate/SPECThead/daughters/insert box
/gate/collimator/geometry/setXLength 23 cm
/gate/collimator/geometry/setYLength 23 cm
/gate/collimator/geometry/setZLength 6 cm
/gate/collimator/setMaterial Lead
/gate/collimator/placement/setTranslation 0. 0. 4.3 cm
/gate/collimator/attachPhantomSD
/gate/collimator/vis/setColor black

/gate/collimator/daughters/name hole
/gate/collimator/daughters/insert hexagone
/gate/hole/placement/setRotationAxis 0 0 1
/gate/hole/placement/setRotationAngle 90 deg
/gate/hole/geometry/setHeight 6. cm
/gate/hole/geometry/setRadius 3.81 mm
/gate/hole/setMaterial Air
/gate/hole/vis/setColor grey

/gate/hole/repeaters/insert cubicArray
/gate/hole/cubicArray/setRepeatNumberX 26
/gate/hole/cubicArray/setRepeatNumberY 15
/gate/hole/cubicArray/setRepeatNumberZ 1
/gate/hole/cubicArray/setRepeatVector 8.236 14.421  0. mm

/gate/hole/repeaters/insert linear
/gate/hole/linear/setRepeatNumber 2
/gate/hole/linear/setRepeatVector 4.118 7.211 0. mm


#CRYSTAL

/gate/SPECThead/daughters/name crystal
/gate/SPECThead/daughters/insert box
/gate/crystal/geometry/setXLength 23 cm
/gate/crystal/geometry/setYLength 23 cm
/gate/crystal/geometry/setZLength 1. cm
/gate/crystal/placement/setTranslation 0. 0. 0.8 cm
/gate/crystal/setMaterial NaI
/gate/crystal/attachCrystalSD
/gate/crystal/vis/setColor yellow
/gate/systems/SPECThead/crystal/attach crystal


#BACK-COMPARTMENT

/gate/SPECThead/daughters/name backcomp
/gate/SPECThead/daughters/insert box
/gate/backcomp/geometry/setXLength 23 cm
/gate/backcomp/geometry/setYLength 23 cm
/gate/backcomp/geometry/setZLength 6.6 cm
/gate/backcomp/placement/setTranslation 0. 0. -3. cm
/gate/backcomp/setMaterial Glass
/gate/backcomp/vis/setColor cyan
/gate/backcomp/attachPhantomSD


#PHANTOM

/gate/world/daughters/name Phantom
/gate/world/daughters/insert cylinder
/gate/Phantom/setMaterial Water
/gate/Phantom/vis/setColor grey
/gate/Phantom/geometry/setRmax 5. cm
/gate/Phantom/geometry/setHeight 20. cm
/gate/Phantom/attachPhantomSD


/vis/viewer/update

#PHYSICS (SONST FEHLER)
/gate/physics/gamma/selectPhotoelectric standard
/gate/physics/gamma/selectCompton standard
/gate/physics/gamma/selectGammaConversion standard
/gate/physics/gamma/selectRayleigh inactive

#CUTS FOR THE PRODUCTION OF SECONDARY PARTICLES- FOR FAST #SIMULATION
#/gate/physics/setElectronCut 1000. km
#/gate/physics/setXRayCut 1. eV
#/gate/physics/setDeltaRayCut 1. GeV


#INITIALIZE
/run/initialize


#SOURCE

/gate/source/addSource Y90beta
/gate/source/Y90beta/setActivity 1000000. becquerel
/gate/source/Y90beta/gps/particle e-
/gate/source/Y90beta/gps/energytype Mono
/gate/source/Y90beta/gps/monoenergy 1.0 MeV
/gate/source/Y90beta/gps/type Volume
/gate/source/Y90beta/gps/shape Cylinder
/gate/source/Y90beta/gps/radius 1. cm
/gate/source/Y90beta/gps/halfz 4. cm
/gate/source/Y90beta/gps/centre 0. 0. 0. cm
/gate/source/Y90beta/gps/confine NULL


#DIGITIZER

/gate/digitizer/Singles/insert adder


#OUTPUT

/gate/output/root/setFileName /home/teicherk/Gate/SPbtest61
/gate/output/root/setRootSinglesAdderFlag 1
/gate/output/root/setRootSinglesFlag 1

/gate/output/ascii/setFileName /home/teicherk/Gate/SPbtest61
/gate/output/ascii/setOutFileHitsFlag 0
/gate/output/ascii/setOutFileSinglesAdderFlag 0
/gate/output/ascii/setOutFileSinglesFlag 0


/gate/output/projection/disable


#START
/gate/application/setTimeSlice 5. s
/gate/application/setTimeStart 0. s
/gate/application/setTimeStop 5. s
/gate/application/startDAQ

-- 
GMX FreeMail: 1 GB Postfach, 5 E-Mail-Adressen, 10 Free SMS.
Alle Infos und kostenlose Anmeldung: http://www.gmx.net/de/go/freemail