[Gate-users] [Gate-devel] problem with parametrised collimator setup ...

[Alex Vergara Gil]

Dear Uwe

the first thing I have noticed is the header of GateSPECTHeadSystem
class which says

/*! \class  GateSPECTHeadSystem
    \brief  The GateSPECTHeadSystem is a model for SPECT scanners and
gamma cameras

    - GateSPECTHeadSystem - by Daniel.Strul at iphe.unil.ch

    - A GateSPECTHeadSystem is a model for SPECT scanners and gamma
cameras. The component tree
      is a linear hierarchy of 3 components, the base, the crystal and the pixel

    - The base is the SPECT head itself (including casing, colli,
shielding...). You may
      have multiple heads by repeating the base with a
'ring-repeater'. If you do so,
      the heads can be at 180, 90, 120 degrees (be careful when you
choose the angle: it
      should be a multiple of the angular step of the scanner). You
can also have the heads
      rotating using a 'rotation-move'. Note that, if you have
multiple heads, the rotation
      axis and the repetition axis (i.e. the scanner axis) must be the same.

    - The level below is the crystal, which can be monoblock or
pixellated. It must be square
      in shape (round crystals not handled)

    - The level below is optional, and is meant to be used for
pixellated gamma-cameras.

    - Note that this component tree does not include a component for
the collimator, as I
      had the feeling reconstruction did not have an urgent need of
collimator information.
      If needed, such a component could be added to this model.

    - When this system is used, it will automatically insert into the
vector of output modules
      2 modules dedicated to SPECT, GateToProjectionSet and
GateToInterfile, to build
      projections and save them respectively.

      \sa GateProjectionSet, GateToProjectionSet, GateToInterfile
*/

The most relevant part here is

- Note that this component tree does not include a component for the
collimator, as I
      had the feeling reconstruction did not have an urgent need of
collimator information.
      If needed, such a component could be added to this model.

So no collimator is included yet to the SPECThead model, I think here
lies the issue. Another thing is that this class needs to be updated
ASAP, because it also included a GateToInterfile when a more general
GateImage is required for every other image format.

So we need to contact Daniel Strul (the original creator) to clarify this class

Regards

Alex Vergara



2015-07-07 11:40 GMT-04:00, Alex Vergara Gil <alexvergaragil at gmail.com>:
> Dear Uwe
>
> I have tried your example with the last update on Gate's git branch
> develop and I can reproduce your bug
>
> GATE object:        'systems/SPECThead'
> Components:
>
> GATE object:        'systems/SPECThead/base'
> Attached to volume: SPECThead
> Nb of children:       1
>
>         GATE object:        'systems/SPECThead/crystal'
>         Attached to volume: crystal
>         Nb of children:       1
>
>                 GATE object:        'systems/SPECThead/pixel'
>                 Attached to volume: ---
> [Physic-0] The following Geant4's physic-list is enabled :emstandard
> [Core-0] Initialization of geometry
> GateVVolume.cc (l.189): The material of the volume hole is not defined.
>                 Nb of children:       0
> tems/SPECThead/pixel'
>                 Attached to volume: ---
>
> It seems something related to the SPECThead system since it is declared in
>
> alex at alex-kubuntu:~/Gate/Correr/fanbeam
> collimator/FanBeam_Gamma_Camera/FanBeam_Gamma_Camera$ grep -r 'hole'
> colli_parallelbeam.mac:#specify the hole radius
> colli_parallelbeam.mac:#distance between the holes
> camera_FB.mac:#set the septal thickness to the required distance
> between the holes
> camera_FB.mac:#specify the hole radius
> camera_PB.mac:#specify the hole radius
> camera_PB.mac:#distance between the holes
> colli_fanbeam.mac:#set the septal thickness to the required distance
> between the holes
> colli_fanbeam.mac:#specify the hole radius
> camera_FB_save.mac:#set the septal thickness to the required distance
> between the holes
> camera_FB_save.mac:#specify the hole radius
> camera_FB_save.mac:#set the septal thickness to the required distance
> between the holes
> camera_FB_save.mac:#specify the hole radius
> camera_Tc.mac:# Insert the first hole of air in the collimator
> camera_Tc.mac:/gate/collimator/daughters/name
> hole
> camera_Tc.mac:/gate/hole/geometry/setHeight
>          3. cm
> camera_Tc.mac:/gate/hole/geometry/setRadius
>          .15 cm
> camera_Tc.mac:/gate/hole/placement/setRotationAxis                      0 1
> 0
> camera_Tc.mac:/gate/hole/placement/setRotationAngle                     90
> deg
> camera_Tc.mac:/gate/hole/setMaterial
>                  Air
> camera_Tc.mac:# Repeat the hole in an array
> camera_Tc.mac:/gate/hole/repeaters/insert
>          cubicArray
> camera_Tc.mac:/gate/hole/cubicArray/setRepeatNumberX                    1
> camera_Tc.mac:/gate/hole/cubicArray/setRepeatNumberY                    52
> camera_Tc.mac:/gate/hole/cubicArray/setRepeatNumberZ                    44
> camera_Tc.mac:/gate/hole/cubicArray/setRepeatVector
>  0. 0.36  0.624 cm
> camera_Tc.mac:# Repeat these holes in a linear
> camera_Tc.mac:/gate/hole/repeaters/insert
>          linear
> camera_Tc.mac:/gate/hole/linear/setRepeatNumber
>          2
> camera_Tc.mac:/gate/hole/linear/setRepeatVector
>          0. 0.18 0.312 cm
> camera_I_131.mac:# Insert the first hole of air in the collimator
> camera_I_131.mac:/gate/collimator/daughters/name
>          hole
> camera_I_131.mac:/gate/hole/geometry/setHeight
>          5. cm
> camera_I_131.mac:/gate/hole/geometry/setRadius
>          .15 cm
> camera_I_131.mac:/gate/hole/placement/setRotationAxis                   0 1
> 0
> camera_I_131.mac:/gate/hole/placement/setRotationAngle                  90
> deg
> camera_I_131.mac:/gate/hole/setMaterial
>                          Air
> camera_I_131.mac:# Repeat the hole in an array
> camera_I_131.mac:/gate/hole/repeaters/insert
>          cubicArray
> camera_I_131.mac:/gate/hole/cubicArray/setRepeatNumberX
>          1
> camera_I_131.mac:/gate/hole/cubicArray/setRepeatNumberY
>          52
> camera_I_131.mac:/gate/hole/cubicArray/setRepeatNumberZ
>          44
> camera_I_131.mac:/gate/hole/cubicArray/setRepeatVector
>  0. 0.36  0.624 cm
> camera_I_131.mac:# Repeat these holes in a linear
> camera_I_131.mac:/gate/hole/repeaters/insert
>