[Gate-users] Regarding ct dosimetry with voxelized moving phantom

[Theodora Kostou]

Hello, 

I am running Gate version 7.2 and I am trying to extract the dose from a
moving phantom in order to simulate a CT scan. 

>From the code below I can visualize the rotation of the xcat phantom but
I take zero counts as output. 

/control/execute mac/verbose.mac
#/control/execute mac/visu.mac
/gate/geometry/setMaterialDatabase data/GateMaterials.db
/gate/world/setMaterial Vacuum
/gate/world/geometry/setXLength 1000. cm
/gate/world/geometry/setYLength 1000. cm
/gate/world/geometry/setZLength 1000. cm
# Container to rotate the phantom
/gate/world/daughters/name CT
/gate/world/daughters/insert box
/gate/CT/setMaterial Vacuum
/gate/CT/geometry/setXLength 640. mm
/gate/CT/geometry/setYLength 640. mm
/gate/CT/geometry/setZLength 1825. mm
/gate/CT/vis/forceWireframe
/gate/CT/moves/insert rotation
/gate/CT/rotation/setSpeed 45 deg/s
/gate/CT/rotation/setAxis 0 0 1
# Patient phantom (128x128x128 voxels, 5 mm voxelize)
/gate/CT/daughters/name xcat
/gate/CT/daughters/insert ImageRegularParametrisedVolume
/gate/xcat/geometry/setImage data/xcat.h33
/gate/xcat/geometry/setRangeToMaterialFile data/Range.dat
/gate/xcat/geometry/TranslateTheImageAtThisIsoCenter 0 0 0 mm
/gate/xcat/placement/setTranslation -320 -320 -912.5 mm
#/gate/xcat/placement/setRotationAxis 1 0 0
#/gate/xcat/placement/setRotationAngle 0 deg
/gate/xcat/vis/forceWireframe
#/gate/xcat/setSkipEqualMaterials 0
/gate/CT/attachPhantomSD
/gate/xcat/attachPhantomSD 

/gate/physics/addProcess PhotoElectric
/gate/physics/processes/PhotoElectric/setModel StandardModel
/gate/physics/addProcess Compton
/gate/physics/processes/Compton/setModel StandardModel
/gate/physics/addProcess RayleighScattering
/gate/physics/processes/RayleighScattering/setModel PenelopeModel
/gate/physics/addProcess ElectronIonisation
/gate/physics/processes/ElectronIonisation/setModel StandardModel e-
/gate/physics/addProcess Bremsstrahlung
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-
/gate/physics/addProcess MultipleScattering e-
/gate/physics/processList Enabled
/gate/physics/processList Initialized 

/gate/actor/addActor DoseActor doseDistribution
/gate/actor/doseDistribution/save output/myactor.mhd
/gate/actor/doseDistribution/attachTo xcat
/gate/actor/doseDistribution/stepHitType random
/gate/actor/doseDistribution/setPosition 0 0 0 mm
/gate/actor/doseDistribution/setVoxelSize 5 5 5 mm
/gate/actor/doseDistribution/enableEdep true
/gate/actor/doseDistribution/enableUncertaintyEdep true
/gate/actor/doseDistribution/enableDose true
/gate/actor/doseDistribution/enableSquaredEdep false
/gate/actor/doseDistribution/enableNumberOfHits false
/gate/actor/doseDistribution/saveEveryNSeconds 10 

/gate/source/addSource xraygun gps
/gate/source/xraygun/gps/verbose 1
/gate/source/xraygun/gps/particle gamma
/gate/source/xraygun/gps/energytype Arb
/gate/source/xraygun/gps/histname arb
/gate/source/xraygun/gps/emin 12.00 keV
/gate/source/xraygun/gps/emax 120.00 keV
/control/execute mac/histogram.mac 
/gate/source/xraygun/gps/arbint Lin
/gate/source/xraygun/gps/type Plane
/gate/source/xraygun/gps/shape Rectangle
/gate/source/xraygun/gps/halfx 0.9 mm
/gate/source/xraygun/gps/halfy 0.8 mm
/gate/source/xraygun/gps/centre -60. 0. 0. cm
/gate/source/xraygun/gps/angtype iso
/gate/source/xraygun/gps/pos/rot1 1 0 0
/gate/source/xraygun/gps/pos/rot2 0 1 0
#cone angle 
/gate/source/xraygun/gps/mintheta 88.1 deg
/gate/source/xraygun/gps/maxtheta 91.9 deg
#fan angle 
/gate/source/xraygun/gps/minphi 155. deg
/gate/source/xraygun/gps/maxphi 205.0 deg
/gate/source/list 

/gate/run/initialize
/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto 
/gate/application/setTimeStart 0 s
/gate/application/setTimeStop 10. s
/gate/application/setTimeSlice 1. s
/gate/application/setTotalNumberOfPrimaries 10
/gate/application/startDAQ 

Is something wrong or missing in the code? 

I would appreciate any kind of advice. 

Thank you 

Theodora 

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