[Gate-users] SIMULATION SPECT

[LORENA VÁZQUEZ CANELAS]

Hi Gaters,

I'm doing a simulation of SPECT with size of crystal of (59.1x44.5) cm and
thickness 0.95 cm.
I want to measure the sensitivity of SPECT, but it doesn't measure well in
ROOT and I can't see any archives.

My big problem is the design of the collimator. I'm going to use one of
type LEHR with the following characteristics:
- Number of holes: 148000 and hole shape hexagonal.
- Hole Length: 2.405 cm = 24.05 mm
- Hole diameter: 0.111 cm = 1.11 mm

The code that i wrote for the collimator is:
# Create a full volume defining the shape of the collimator
/gate/SPECThead/daughters/name                 collimator
/gate/SPECThead/daughters/insert                  box
/gate/collimator/geometry/setXLength            2.405 cm
/gate/collimator/geometry/setYLength            59.1 cm
/gate/collimator/geometry/setZLength            44.5 cm
/gate/collimator/placement/setTranslation      -2. 0. 0. cm
/gate/collimator/setMaterial                              Lead
/gate/collimator/vis/setColor                            black
#/gate/collimator/vis/forceWireframe

# Insert the first hole of air in the collimator
/gate/collimator/daughters/name                     hole
/gate/collimator/daughters/insert                     hexagone
*/gate/hole/geometry/setHeight              2.405 cm*
*/gate/hole/geometry/setRadius              0.0555 cm*
/gate/hole/placement/setRotationAxis             0 1 0
/gate/hole/placement/setRotationAngle          90 deg
/gate/hole/setMaterial                                      Air

# Repeat the hole in an array
/gate/hole/repeaters/insert                                 cubicArray
/gate/hole/cubicArray/setRepeatNumberX         1
/gate/hole/cubicArray/setRepeatNumberY         164
/gate/hole/cubicArray/setRepeatNumberZ         71
/gate/hole/cubicArray/setRepeatVector             0. 0.36  0.624 cm

# Repeat these holes in a linear
/gate/hole/repeaters/insert                       linear
/gate/hole/linear/setRepeatNumber         2
/gate/hole/linear/setRepeatVector            0. 0.18 0.312 cm

If i change the value 2.405 cm by 3cm and 0.0555cm by 0.15 cm, i can see
all archives in ROOT, but i I have no such measures!!
You know if RADIUS is half the diameter of the hole?

Thanks in advance

-- 

LORENA VÁZQUEZ CANELAS
PhD student, Research staff
LRI - Datación (Despacho Técnicos)
Calle Espejo, s/n, 37008 - Salamanca (Spain)
TEL: +34 923 294930 ext:4932 email: lorena.vazquez at usal.es
web: http://campus.usal.es/~lri/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20150929/67634e5f/attachment.html>