[Gate-users] Simulating a Voxelized Phantom

[Marc Chamberland]

Hi Taro!

Look at the message it's giving you carefully, then look at your macro.

You gave an activity range data file for your phantom definition, which is wrong.

For a voxelized phantom, you need to give it an attenuation range data file. In the activity range file, there are no materials; that's why the message says it's looking for material " " in your material database.

Have a look at the Voxelized Phantom/Source example in the Gate folder. You need an activity range file and an attenuation range file.

Cheers!
Marc



__________________________

Marc Chamberland, MSc
PhD candidate
Department of Physics
Carleton University
Ottawa (ON)

Le 2012-09-09 à 10:14 PM, <dankemori at yahoo.co.jp<mailto:dankemori at yahoo.co.jp>>
 a écrit :

Hello dear GATE users.

I'm trying to simulate a Voxelized Phantom(Brain Phantom) in a PET system now.
I defined phantom and source according to Users Guide as follows

#  B R A I N   P H A N T O M
/gate/world/daughters/name hoffman
/gate/world/daughters/insert parameterizedBoxMatrix
/gate/hoffman/geometry/insertReader interfile
/gate/hoffman/interfileReader/insertTranslator range
/gate/hoffman/interfileReader/rangeTranslator/readTable activityRange.dat
/gate/hoffman/interfileReader/rangeTranslator/describe 1
/gate/hoffman/interfileReader/readFile hof3.h33
/gate/hoffman/placement/setTranslation 0. 0. 0. mm
/gate/hoffman/placement/setRotationAxis 1 0 0
/gate/hoffman/placement/setRotationAngle 0 deg
/gate/hoffman/attachVoxelPhantomSD

/gate/hoffman/addOutput doseOutput
/gate/output/doseOutput/saveUncertainty true
/gate/output/doseOutput/setFileName hoffman_paraDose.bin

# B R A I N
/gate/source/addSource voxel voxel
/gate/source/voxel/reader/insert interfile
/gate/source/voxel/interfileReader/translator/insert range
/gate/source/voxel/interfileReader/rangeTranslator/readTable activityRange.dat
/gate/source/voxel/interfileReader/rangeTranslator/describe 1
/gate/source/voxel/interfileReader/readFile hof.h33
/gate/source/voxel/interfileReader/verbose 1
/gate/source/voxel/setType backtoback
/gate/source/voxel/setAccolinearityFlag true
/gate/source/voxel/setAccoValue 0.5 deg
/gate/source/voxel/gps/particle gamma
/gate/source/voxel/gps/energytype Mono
/gate/source/voxel/gps/monoenergy 0.511 MeV
/gate/source/voxel/gps/angtype iso
/gate/source/voxel/gps/mintheta 0. deg
/gate/source/voxel/gps/maxtheta 180. deg
/gate/source/voxel/gps/minphi 0. deg
/gate/source/voxel/gps/maxphi 360. deg
/gate/source/voxel/gps/confine NULL
/gate/source/voxel/setPosition -128. -128. 0. mm
/gate/source/voxel/dump 1



But I had the following warnings:

-----------------------------------------------------------------------------------------------

PreInit> /gate/world/daughters/insert parameterizedBoxMatrix
/gate/world/daughters/insert parameterizedBoxMatrix
GateVoxelBoxParameterizedInserter::ConstructGeometry - Warning ! ConstructGeometry called without a reader
world_L is registered to the default region.
  Volume of the same name and copy number ("world_P", copy 0) still exists and is being used.
  WARNING: This does not necessarily guarantee it's the same
  volume you originally specified in /vis/scene/add/.
/vis/scene/notifyHandlers scene-0

PreInit> /gate/hoffman/geometry/insertReader interfile
/gate/hoffman/geometry/insertReader interfile

PreInit> /gate/hoffman/interfileReader/insertTranslator range
/gate/hoffman/interfileReader/insertTranslator range

PreInit> /gate/hoffman/interfileReader/rangeTranslator/readTable activityRange.dat
/gate/hoffman/interfileReader/rangeTranslator/readTable activityRange.dat
GateGeometryVoxelRangeTranslator::ReadFile : fileName: activityRange.dat
   GateMDBFile::ReadItem: I could not find the item '' in section [Materials] of the material database.
   GateMaterialDatabase: could not download the definition for ''
GateMaterialDatabase: failed to read the material '' in the database file!

PreInit> /gate/hoffman/interfileReader/rangeTranslator/describe 1
/gate/hoffman/interfileReader/rangeTranslator/describe 1
 Range Translator

!!! GateMDBFile::LookForElementIndbStream:
I could not find any section [Materials] in the material database file.
You should check the database file.
Computation aborted!!!


*** G4Exception: Aborting execution ***
Aborted

-----------------------------------------------------------------------------------------------


Any kind of hints or solutions to my problem are appreciated.


Thank you

*************************
Taro Kanamaru

Tohoku University Japan

email: dankemori at yahoo.co.jp<mailto:dankemori at yahoo.co.jp>
*************************
_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org<mailto:Gate-users at lists.opengatecollaboration.org>
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20120910/5fd9e0e3/attachment.html>