[Gate-users] NCAT

Sat May 3 23:07:40 CEST 2008

Hi,

My guess from what you have posted here is that you are  not inserting the activity range file before loading your source file.  

Try inserting the range file first  and then read the interfile if you are not doing that already 
eg-
/gate/source/voxel/interfileReader/translator/insert range
/gate/source/voxel/interfileReader/rangeTranslator/readTable activityRange.dat
/gate/source/voxel/interfileReader/readFile hof.h33

Bosky

________________________________

From: gate-users-bounces at lists.healthgrid.org on behalf of vieira.lina at sapo.pt
Sent: Sat 5/3/2008 4:45 PM
To: gate-users at lists.healthgrid.org
Subject: Fwd: [Gate-users] NCAT

        Hello,

        I`m beginning to use the NCAT phantom in my simulations. But I have a problem...

I used the parameter file general.samp.par to make my file
visible_phantom_act.bin.
Then, I decoded the file visible_phantom_act.bin to UNSIGNED INTEGER, and I got
the file ncat.i33.
Afterwords, I made the file ncat.h33 according to the parameter file titled
"general.samp.par".

My problem appeared when I made the file activityRange.dat.

I don`t know how I should proceed.

I saw the activity for each organ, and in each frame, I added different activity
for different organs.

For example:

Frame 61/202 =  152 Bq

Body = 2 Bq
Liver = 75 Bq
Heart = 75 Bq

My file activityRange.dat is:
 15

1     3        152
4      11      227
12    12      181
13     22     302
23     24      227
25     25      272
26     27      287
28     37      212
38     50      152
51      51     115
52      52       77
53      62     152
63       63    116
64       89       77
90       202       2

But when I ran the simulations I got the following message:

Range   0 : imageValue in [    1.000e+00 , 3.000e+00 ]  ---> activity (Bq)
1.520e+02
        Range   1 : imageValue in [    4.000e+00 , 1.100e+01 ]  ---> activity
(Bq) 2.270e+02
        Range   2 : imageValue in [    1.200e+01 , 1.200e+01 ]  ---> activity
(Bq) 1.810e+02
        Range   3 : imageValue in [    1.300e+01 , 2.200e+01 ]  ---> activity
(Bq) 3.020e+02
        Range   4 : imageValue in [    2.300e+01 , 2.400e+01 ]  ---> activity
(Bq) 2.270e+02
        Range   5 : imageValue in [    2.500e+01 , 2.500e+01 ]  ---> activity
(Bq) 2.720e+02
        Range   6 : imageValue in [    2.600e+01 , 2.700e+01 ]  ---> activity
(Bq) 2.870e+02
        Range   7 : imageValue in [    2.800e+01 , 3.700e+01 ]  ---> activity
(Bq) 2.120e+02
        Range   8 : imageValue in [    3.800e+01 , 5.000e+01 ]  ---> activity
(Bq) 1.520e+02
        Range   9 : imageValue in [    5.100e+01 , 5.100e+01 ]  ---> activity
(Bq) 1.150e+02
        Range  10 : imageValue in [    5.200e+01 , 5.200e+01 ]  ---> activity
(Bq) 7.700e+01
        Range  11 : imageValue in [    5.300e+01 , 6.200e+01 ]  ---> activity
(Bq) 1.520e+02
        Range  12 : imageValue in [    6.300e+01 , 6.300e+01 ]  ---> activity
(Bq) 1.160e+02
        Range  13 : imageValue in [    6.400e+01 , 8.900e+01 ]  ---> activity
(Bq) 7.700e+01
        Range  14 : imageValue in [    9.000e+01 , 2.020e+02 ]  ---> activity
(Bq) 2.000e+00
GateSourceVoxelImageReader::ReadFile : fileName: ncat.h33
 Header read from       'ncat.h33'
 Data file name         'ncat.i33'
 Nb of planes:           202
 Nb of pixels per plane: 128 128
 Pixel size:             3.125e+00 3.125e+00
 Slice thickness:        3.125e+00
 Matrix size:            4.000e+02 4.000e+02
 Data type:              UNSIGNED INTEGER

 Ignoring confine condition
Source ---------------> voxel_NCAT
  ID                  : 0
  type                :
  startTime (s)       : 0.000e+00
  time (s)            : 0.000e+00
  forcedUnstable      : 0
  verboseLevel        : 0
-----------------------
  Voxel reader ---------->   number of voxels       : 0

total activity (Bq)    : 0.000e+00

position  (mm)         : 0.000e+00 0.000e+00 0.000e+00
 voxel size  (mm)    : 3.125e+00 3.125e+00 3.125e+00
    GPS info ---------------->
    particle                 : gamma
    SourcePosType            : Point
    Shape                    : NULL
    halfx,halfy,halfz  (mm)  : 0.000e+00 0.000e+00 0.000e+00
    Radius             (mm)  : 0.000e+00
    Radius0            (mm)  : 0.000e+00
    CentreCoords       (mm)  : 0.000e+00 0.000e+00 0.000e+00
    Direction         : (1.000e+00,0.000e+00,0.000e+00)
    Confine                  : 0
    VolName                  : NULL
    EnergyDisType            : Mono
    MonoEnergy         (MeV) : 1.400e-01
    MinTheta, MaxTheta (deg) : 0.000e+00 1.800e+02
    MinPhi, MaxPhi     (deg) : 0.000e+00 3.600e+02
    --------------------------

Time set to (s) 0.000e+00

      .....
========= Table of registered couples ==============================

....
====================================================================

----- HepJamesRandom engine status -----
 Initial seed = 9876
 u[] = 0.174685 ....... 0.485254 0.222273

 c = 0.0216029, cd = 0.456233, cm = 1
 i97 = 96, u[i97] = 0.222273
 j97 = 32, u[j97] = 0.370831
----------------------------------------
[GateSource::Update] Source name: voxel_NCAT
GateSource::GetNextTime : next time (s) 1.798e+299
GateSourceMgr::PrepareNextEvent : m_time > m_timeLimit. No vertex generated

----- HepJamesRandom engine status -----
 Initial seed = 9876
 u[] = 0.211593 .............0.475358 0.406137

 c = 0.665798, cd = 0.456233, cm = 1
 i97 = 78, u[i97] = 0.940329
 j97 = 14, u[j97] = 0.360893
----------------------------------------
Time set to (s) 5.000e+01
 Current Output Root File: 0xc8bcac0 named benchSPECT.root
GateToRoot: ROOT: files writing...
GateToRoot: ROOT: files closing...
Idle>

Can any help me please?

Thank you very much.

Best regards,

Lina Vieira

----- Finalizar mensagem reenviada -----