[Gate-users] NCAT

Tue Apr 29 11:44:19 CEST 2008

	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

Body = 2
Liver = 75
Heart = 75

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
-------------- next part --------------
_______________________________________________
Gate-users mailing list
Gate-users at lists.healthgrid.org
http://lists.healthgrid.org/mailman/listinfo/gate-users