[Gate-users] How to create a Hoffman phantom?

Tsukahara Tatsuya Tatsuya.Tsukahara at mc7.seikyou.ne.jp
Fri Jan 7 06:27:41 CET 2011


Hello dear Gate-Users,

I have installed gate_v5.0.0_p01 and geant4.9.1.p03.
And I'm simulating a PET system.

Now, I'm trying to create a Hoffman phantom.
But I had the following warnings:
---------------------------------------------------------------------
GateMaterialDatabase: did not load properties table for material  'Vacuum'.
This is only a problem when OPTICAL PHOTONS are transported in this
material.'
---------------------------------------------------------------------
This kind of warning also appears with: Lung, Brest, Plastic.

In addition,
---------------------------------------------------------------------
GateGeometryVoxelInterfileReader::ReadFile: WARNING: voxel not added
(material translation not found); value: 1
---------------------------------------------------------------------
And finally,
---------------------------------------------------------------------
PANIC! - G4SmartVoxelHeader::BuildNodes()
         Daughter physical volume hof_brainVoxel_PVP
         is entirely outside mother logical volume hof_brain_L !!

*** G4Exception : InvalidSetup
      issued by : G4SmartVoxelHeader::BuildNodes()
Overlapping daughter with mother volume.
-------------------------------------------------------------------------

Below is the code that I used for the phantom and the source.
I read the user's guide but I didn't understand which parameter I
should change.
Any kind of hints or solutions to my problem are appreciated.

Additionally I have troubles with "point source".
I obtained correct data with Volume source but I couldn't get enough
counts with Point source.
Code for the point source is after Hoffman code.
Any kind of hints or solutions to this problem too are appreciated.


Thank you for your help.

************************************************
 Tohoku University, Japan
 Tatsuya Tsukahara
 E-Mail:Tatsuya.Tsukahara at mc7.seikyou.ne.jp
*************************************************


CODE
======================================================
#PHANTOM
/gate/world/daughters/name hof_brain
/gate/world/daughters/insert parameterizedBoxMatrix
/gate/hof_brain/geometry/insertReader interfile
/gate/hof_brain/interfileReader/insertTranslator range
/gate/hof_brain/interfileReader/rangeTranslator/readTable range.dat
/gate/hof_brain/interfileReader/rangeTranslator/describe 1
/gate/hof_brain/interfileReader/readFile hof3.h33
/gate/hof_brain/placement/setTranslation  0. 0. 0. mm
/gate/hof_brain/placement/setRotationAxis 1 0 0
/gate/hof_brain/placement/setRotationAngle 0 deg
/gate/hof_brain/attachVoxelPhantomSD

/gate/hof_brain/addOutput doseOutput
/gate/output/doseOutput/saveUncertainty true
/gate/output/doseOutput/setFileName hof_brain_paraDose.bin

#SOURCE
/gate/source/addSource voxel_brain voxel

/gate/source/voxel_brain/reader/insert interfile
/gate/source/voxel_brain/interfileReader/translator/insert range
/gate/source/voxel_brain/interfileReader/rangeTranslator/readTable
activityRange.dat
/gate/source/voxel_brain/interfileReader/rangeTranslator/describe 1
/gate/source/voxel_brain/interfileReader/readFile hof3.h33

/gate/source/voxel_brain/setType backtoback
/gate/source/voxel_brain/gps/particle gamma
/gate/source/voxel_brain/setForcedUnstableFlag true
/gate/source/voxel_brain/setForcedHalfLife 6586.2 s
/gate/source/voxel_brain/gps/energytype Mono
/gate/source/voxel_brain/gps/monoenergy 0.511 MeV
/gate/source/voxel_brain/setPosition   -128. -128. -55. mm
/gate/source/voxel_brain/gps/confine NULL
/gate/source/voxel_brain/gps/angtype iso
/gate/source/voxel_brain/dump 1

/gate/source/list

========================================================================
#POINT SOURCE
/gate/source/addSource twogamma
/gate/source/twogamma/setActivity 37000000. becquerel
/gate/source/twogamma/gps/particle gamma
/gate/source/twogamma/setType backtoback
/gate/source/twogamma/setAccolinearityFlag false
/gate/source/twogamma/setAccoValue 0.5 deg
/gate/source/twogamma/gps/centre 0.0 0.0 0.0 mm
/gate/source/twogamma/gps/energytype Mono
/gate/source/twogamma/gps/monoenergy 0.511 MeV
/gate/source/twogamma/gps/type Point
#/gate/source/twogamma/gps/type Volume
#/gate/source/twogamma/gps/shape Sphere
#/gate/source/twogamma/gps/shape Cylinder
#/gate/source/twogamma/gps/radius 0.3 mm
#/gate/source/twogamma/gps/halfz 20. mm
#/gate/source/twogamma/gps/confine NULL
/gate/source/twogamma/gps/angtype iso
/gate/source/twogamma/gps/mintheta 90. deg
/gate/source/twogamma/gps/maxtheta 90. deg
/gate/source/twogamma/gps/minphi 0. deg
/gate/source/twogamma/gps/maxphi 360. deg
/gate/source/list



_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users



More information about the Gate-users mailing list