<div dir="ltr"><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small">Hi Gaters,<br><br>I'm doing a simulation of SPECT with size of crystal of (59.1x44.5) cm and thickness 0.95 cm. <br>I want to measure the sensitivity of SPECT, but it doesn't measure well in ROOT and I can't see any archives.<br><br>My big problem is the design of the collimator. I'm going to use one of type LEHR with the following characteristics:<br>- Number of holes: 148000 and hole shape hexagonal.<br>- Hole Length: 2.405 cm = 24.05 mm<br>- Hole diameter: 0.111 cm = 1.11 mm<br><br>The code that i wrote for the collimator is:<br># Create a full volume defining the shape of the collimator<br>/gate/SPECThead/daughters/name                 collimator<br>/gate/SPECThead/daughters/insert                  box<br>/gate/collimator/geometry/setXLength            2.405 cm<br>/gate/collimator/geometry/setYLength            59.1 cm<br>/gate/collimator/geometry/setZLength            44.5 cm<br>/gate/collimator/placement/setTranslation      -2. 0. 0. cm<br>/gate/collimator/setMaterial                              Lead<br>/gate/collimator/vis/setColor                            black<br>#/gate/collimator/vis/forceWireframe<br><br># Insert the first hole of air in the collimator<br>/gate/collimator/daughters/name                     hole<br>/gate/collimator/daughters/insert                     hexagone<br><b>/gate/hole/geometry/setHeight              2.405 cm</b><br><b>/gate/hole/geometry/setRadius              0.0555 cm</b><br>/gate/hole/placement/setRotationAxis             0 1 0<br>/gate/hole/placement/setRotationAngle          90 deg<br>/gate/hole/setMaterial                                      Air<br><br># Repeat the hole in an array<br>/gate/hole/repeaters/insert                                 cubicArray<br>/gate/hole/cubicArray/setRepeatNumberX         1<br>/gate/hole/cubicArray/setRepeatNumberY         164<br>/gate/hole/cubicArray/setRepeatNumberZ         71<br>/gate/hole/cubicArray/setRepeatVector             0. 0.36  0.624 cm<br><br># Repeat these holes in a linear<br>/gate/hole/repeaters/insert                       linear<br>/gate/hole/linear/setRepeatNumber         2<br>/gate/hole/linear/setRepeatVector            0. 0.18 0.312 cm<br><br>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!!<br>You know if RADIUS is half the diameter of the hole?<br><br>Thanks in advance<br><br clear="all"></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><pre>LORENA VÁZQUEZ CANELAS<br><span>PhD student, </span><span lang="en">Research staff</span><br>LRI - Datación (Despacho Técnicos)<br>Calle Espejo, s/n, 37008 - Salamanca (Spain)<br>TEL: +34 923 294930 ext:4932 email: <a href="mailto:lorena.vazquez@usal.es" target="_blank">lorena.vazquez@usal.es</a><br>web: <a href="http://campus.usal.es/%7Elri/" target="_blank">http://campus.usal.es/~lri/</a></pre></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div>