[Gate-users] Time-activity curves for voxelized phantom/source

[Ida Häggström]

Hi!

I'm trying to implement a dynamic voxelized source, but have run into some
problems. I've followed the users guide, section 7.4 regarding how to set up
my source/phantom, but I'm not successful. =(
First of all, I have two time-activity curves (TACs), one for tissue and one
for the heart. I want a simple phantom/source, a 4x4x4 voxel cube with the
central 2x2x2  voxels corresponding to heart and the remaining outer voxels
being tissue. The actual phantom/source (to the header file) is a 32bit
binary image of the 4x4x4 voxels, voxel value 2 for the inner heart voxels
and voxel value 1 for the outer tissue voxels.

One thing that confuses me in the users guide is that in the source macro,
is that they call for both* activityRange.dat* and *TimeActivity_Tables.dat*.
The first translates the voxel value to an activity, and the second contains
the TACs for the different voxel values. Shouldn't you just use one of them?
And further more, they add the "RTVPhantom". I cannot find anything about
what this is in the guide, and no hits when I google it...

Anyway, when I run I get the error:

 creating a RTV Phantom at Address 0x8fb6ac0
 The Object RTPhantom RTVPhantom is enabled.
GateRTPhantom::AttachToGeometry
 GateRTPhantom::AttachToGeometry : The RTPhantom RTVPhantom is now ATTACHED
to Inserter Object phantom : OK !!!
 Please set it before setting the base file name. Aborting.

*** G4Exception: Aborting execution ***
Gatedied with signal 6


All macros are found below. Thanks in advance for any suggestions! Any
suggestions at all regarding implementing TACs into sources are very
welcome!
Ida


-------------------------------------------------------------------
                   phantom.mac
-------------------------------------------------------------------
#       P H A N T O M
#

/gate/world/daughters/name phantom
/gate/world/daughters/insert compressedMatrix

/gate/phantom/geometry/insertReader interfile
/gate/RTPhantom/insert RTVPhantom
/gate/RTVPhantom/AttachTo phantom

/gate/RTVPhantom/setBaseFileName dynamicPhantom/phantom
/gate/RTVPhantom/setHeaderFileName dynamicPhantom/phantom_header.h33

/gate/RTVPhantom/SetNumberOfFrames 11
/gate/RTVPhantom/SetTimePerFrame 1 s

/gate/phantom/interfileReader/insertTranslator range
/gate/phantom/interfileReader/rangeTranslator/readTable
dynamicPhantom/range.dat
/gate/phantom/interfileReader/rangeTranslator/describe 1
/gate/phantom/attachPhantomSD
/gate/phantom/placement/setTranslation 0. 0. 0. mm
/gate/phantom/interfileReader/describe 1
-------------------------------------------------------------------
               source.mac
-------------------------------------------------------------------
#     S O U R C E
#
/gate/source/addSource voxelSource voxel
/gate/RTVPhantom/AttachToSource voxel

/gate/source/voxelSource/reader/insert interfile
/gate/source/voxelSource/interfileReader/translator/insert range
/gate/source/voxelSource/interfileReader/rangeTranslator/readTable
dynamicPhantom/actRange.dat
/gate/source/voxelSource/interfileReader/SetTimeActivityTablesFrom
dynamicPhantom/TAC_table.dat
/gate/source/voxelSource/interfileReader/SetTimeSampling 1 s

# The deafult position of the voxellized source is in the 1^{st}
# quarter. So the voxellized source has to be shifted over half its
# dimension in the negative direction on each axis
/gate/source/voxelSource/setPosition -40. -40. -40. mm

/gate/source/voxelSource/setType backtoback
/gate/source/voxelSource/gps/particle gamma
/gate/source/voxelSource/gps/energytype Mono
/gate/source/voxelSource/setForcedUnstableFlag true
/gate/source/voxelSource/setForcedHalfLife 1223 s #Carbon-11
/gate/source/voxelSource/gps/monoenergy 511. keV
/gate/source/voxelSource/gps/angtype iso
/gate/source/voxelSource/gps/mintheta 90. deg
/gate/source/voxelSource/gps/maxtheta 90. deg
/gate/source/voxelSource/gps/minphi     0. deg
/gate/source/voxelSource/gps/maxphi   360. deg
/gate/source/voxelSource/gps/confine NULL #voxelPhantom_P

/gate/source/voxelSource/dump 1
/gate/source/list
------------------------------------------------------------------
             phantom_header.h33
------------------------------------------------------------------
!INTERFILE  :=
!GENERAL DATA :=
!GENERAL IMAGE DATA :=
!type of data := PET
imagedata byte order := LITTLEENDIAN
!PET STUDY (General) :=
!PET data type := Image
process status := Reconstructed
!number format := float
!number of bytes per pixel := 4
number of dimensions := 3
matrix axis label [1] := x
!matrix size [1] := 4
scaling factor (mm/pixel) [1] := 20
matrix axis label [2] := y
!matrix size [2] := 4
scaling factor (mm/pixel) [2] := 20
matrix axis label [3] := z
!matrix size [3] := 4
scaling factor (mm/pixel) [3] := 20
image scaling factor[1] := 1
data offset in bytes[1] := 0
quantification units := 1
!END OF INTERFILE :=
------------------------------------------------------------------
             TACs.dat
------------------------------------------------------------------
2
1 1   tissue.dat
2 2   heart.dat
------------------------------------------------------------------
            actRange.dat
------------------------------------------------------------------
2
1 1  0
2 2 10
------------------------------------------------------------------
             tissue.dat
------------------------------------------------------------------
11
0.0 0
1. 100
2. 200
3. 300
4. 500
5. 800
6. 900
7. 910
8. 920
9. 920
10 920
------------------------------------------------------------------
             heart.dat
------------------------------------------------------------------
11
0.0 0.0
1. 100
2. 200
3. 300
4. 500
5. 800
6. 1000
7. 700
8. 400
9. 200
10 100
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20091012/e19bb858/attachment.htm>