[Gate-users] Z position
venialgo_e
venialgo_e at secyt.frba.utn.edu.ar
Thu Feb 14 17:18:53 CET 2008
I changed Z lenght by:
/gate/cristalIoduro/geometry/setZLength 2 cm
/gate/cristalIoduro/placement/setTranslation 0 0 -1 cm
and I have exactly the same histogram.
P.D.: I attached the new macro file and histogram.
Esteban Venialgo
"Chen, Yu" <Yu.Chen at umassmed.edu> ha escrito:
[Ocultar texto citado]
The peak in bottom is due to your geometry setting: the sensitive
detector is the same
as
world size at -z, then no photons can escape your SD. You need to set
your world size
in z
slightly bigger.
None of your three plots makes sense to me. Maybe next time you can
look only those with
small globalPosX, Y events.
Yu Chen, Ph.D.
University of Massachusetts Medical School
Division of Nuclear Medicine
55 Lake Avenue North
Worcester, MA 01655-0243
Phone: (508) 856-6123
Fax: (508) 856-4572
-----Original Message-----
From: gate-users-bounces at lists.healthgrid.org on behalf of venialgo_e
Sent: Thu 2/14/2008 8:16 AM
To: gate-users at lists.healthgrid.org
Subject: [Gate-users] Z position
Hello Gate users,
I'm studying the effect of Depth Of Interaction in a crystal of NaI. I
simulated a detector of 30cm x 40 cm x 2.5 cm using a source with
several emission angles.
the detector is:
/gate/scanner/daughters/name cristalIoduro
/gate/scanner/daughters/insert box
/gate/cristalIoduro/geometry/setXLength 5 cm
/gate/cristalIoduro/geometry/setYLength 5 cm
/gate/cristalIoduro/geometry/setZLength 2.5 cm
/gate/cristalIoduro/placement/setTranslation 0 0 -1.25 cm
/gate/cristalIoduro/setMaterial NaI
/gate/cristalIoduro/vis/setColor red
/gate/cristalIoduro/attachCrystalSD
the source is:
/gate/source/addSource cesio
/gate/source/cesio/setActivity 20000 becquerel
/gate/source/cesio/gps/centre 0. 0. 1 cm
/gate/source/cesio/gps/particle gamma
/gate/source/cesio/gps/energytype Mono
/gate/source/cesio/gps/monoenergy 0.622 MeV
/gate/source/cesio/gps/type Volume
/gate/source/cesio/gps/shape Sphere
/gate/source/cesio/gps/radius 0 cm
/gate/source/cesio/gps/angtype iso
1)In the first case the emission angle is:
/gate/source/cesio/gps/mintheta 0. deg
/gate/source/cesio/gps/maxtheta 0. deg
/gate/source/cesio/gps/minphi 0. deg
/gate/source/cesio/gps/maxphi 360. deg
2)In the second, it is:
/gate/source/cesio/gps/mintheta 0. deg
/gate/source/cesio/gps/maxtheta 3. deg
/gate/source/cesio/gps/minphi 0. deg
/gate/source/cesio/gps/maxphi 360. deg
3)In the last, it is:
/gate/source/cesio/gps/mintheta 0. deg
/gate/source/cesio/gps/maxtheta 20. deg
/gate/source/cesio/gps/minphi 0. deg
/gate/source/cesio/gps/maxphi 360. deg
And The questions about Z positions of the singles are:
1)In the first case , why are there a lot of counts only on the
crystal's top face?.
2)In the second case, why are there more counts in one half of the
crystal than the other?
3)In the third case. why is there a peak of counts, in the crystal's
bottom face ?.
P.D.: I attached the macro files for each case, and gif files with the
Z position's histograms of the singles, also.
Thanks in advance,
Esteban Venialgo
----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.
----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.
----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.
-------------- next part --------------
# ------------------------------------------------------------------------------
# - -
# - Archivo: case3.mac -
# - -
# - -
# - Comentario: Macro file del gate que simula una cristal detector -
# - con una fuente gamma de 622 keV -
# - -
# - Autor: Esteban (Tetu) Venialgo -
# - -
# - Plataforma/Gate: ROOT 5.14 / CLHEP 2.0.3.2 / GEANT 4.9.0 (PATCH) -
# / GATE 3.1.2 -
# - O.S.: Debian 4.0r1 AMD64 Kernel 2.6.18-4 -
# - -
# - Version: 001 -
# - -
# - Ultima modificacion: 11/02/08 -
# - -
# ------------------------------------------------------------------------------
# ------------------------------------------------------------------------------
# - COMANDOS DE VISUALIZACION -
# ------------------------------------------------------------------------------
# ------------------------------------------------------------------------------
# - COMANDOS DE VISUALIZACION -
# ------------------------------------------------------------------------------
/vis/open OGLSX
/vis/viewer/reset
/vis/viewer/zoom 1.5
/vis/viewer/set/viewpointThetaPhi 60 60
/vis/viewer/set/style surface
/vis/drawVolume
/tracking/storeTrajectory 1
/vis/scene/endOfEventAction accumulate
/gate/geometry/enableAutoUpdate
/vis/viewer/set/globalLineWidthScale 4
# /vis/disable
# ------------------------------------------------------------------------------
# - DEFINICIONES DE GEOMETRIA -
# ------------------------------------------------------------------------------
# Defino el world volume (es el volumen global que crea el gate por defecto, que
# es el padre de todas las definciones geometricas, y todos los volumenes,
# y particulas deben estar ahi adentro, sino se anulan).
# WORLD
/gate/world/vis/setColor blue
/gate/world/geometry/setXLength 5 cm
/gate/world/geometry/setYLength 5 cm
/gate/world/geometry/setZLength 5 cm
/gate/world/describe
# SYSTEM
/gate/world/daughters/name scanner
/gate/world/daughters/insert box
/gate/scanner/geometry/setXLength 5 cm
/gate/scanner/geometry/setYLength 5 cm
/gate/scanner/geometry/setZLength 5 cm
/gate/scanner/vis/forceWireframe
/gate/scanner/vis/setColor red
/gate/scanner/setMaterial Air
# EJES
/gate/world/daughters/insert 3axes
# Defino el cristal de NaI(Tl)
/gate/scanner/daughters/name cristalIoduro
/gate/scanner/daughters/insert box
/gate/cristalIoduro/geometry/setXLength 5 cm
/gate/cristalIoduro/geometry/setYLength 5 cm
/gate/cristalIoduro/geometry/setZLength 2 cm
/gate/cristalIoduro/placement/setTranslation 0 0 -1 cm
/gate/cristalIoduro/setMaterial NaI
/gate/cristalIoduro/vis/setColor red
/gate/cristalIoduro/attachCrystalSD
# ------------------------------------------------------------------------------
# - SETTING UP SYSTEM
# ------------------------------------------------------------------------------
/gate/systems/scanner/level1/attach cristalIoduro
/gate/cristalIoduro/attachCrystalSD
# ------------------------------------------------------------------------------
# - SETTING UP PHYSICS
# ------------------------------------------------------------------------------
/gate/physics/gamma/selectRayleigh inactive
/gate/physics/gamma/selectPhotoelectric standard
/gate/physics/gamma/selectCompton standard
/gate/physics/gamma/selectGammaConversion inactive
# INACTIVE SECONDARY ELECTRONS
# /gate/physics/setElectronCut 100. cm
# INACTIVE X - RAYS
# /gate/physics/setXRayCut 1. GeV
# /gate/physics/setDeltaRayCut 1. GeV
# /gate/physics/verbose 2
# Despues de esto el proceso fisico, y las geometrias no se pueden tocar
/run/initialize
# ------------------------------------------------------------------------------
# - DEFINICIONES DE LA FUENTE DE CESIO 137
# ------------------------------------------------------------------------------
# Defino una fuente de Cesio de 622 keV puntual
/gate/source/addSource cesio
/gate/source/cesio/setActivity 20000 becquerel
/gate/source/cesio/gps/centre 0. 0. 1 cm
/gate/source/cesio/gps/particle gamma
/gate/source/cesio/gps/energytype Mono
/gate/source/cesio/gps/monoenergy 0.622 MeV
/gate/source/cesio/gps/type Volume
/gate/source/cesio/gps/shape Sphere
/gate/source/cesio/gps/radius 0 cm
/gate/source/cesio/gps/angtype iso
/gate/source/cesio/gps/mintheta 0. deg
/gate/source/cesio/gps/maxtheta 20. deg
/gate/source/cesio/gps/minphi 0. deg
/gate/source/cesio/gps/maxphi 360. deg
/gate/source/list
# ------------------------------------------------------------------------------
# - DEFINICIONES DE LOS ARCHIVOS DE SALIDA
# ------------------------------------------------------------------------------
# Solamente voy a usar la salida ASCII
/gate/output/ascii/enable
/gate/output/ascii/setOutFileHitsFlag 0
/gate/output/ascii/setOutFileSinglesFlag 1
/gate/output/ascii/setSingleMask 0 0 0 0 0 0 0 1 1 1 1 1
# /gate/output/root/disable
/gate/output/root/enable
/gate/output/root/setFileName tetu
# ------------------------------------------------------------------------------
# - ARRANCO LA SIMULACION
# ------------------------------------------------------------------------------
/gate/application/setTimeSlice 3. s
/gate/application/setTimeStart 0. s
/gate/application/setTimeStop 3. s
/gate/application/startDAQ
-------------- next part --------------
A non-text attachment was scrubbed...
Name: case3moidf.gif
Type: image/gif
Size: 8883 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20080214/1ce46cf2/attachment.gif>
More information about the Gate-users
mailing list