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

[Elnur Sadihov]

In the previous letter I forgot to write one more point:
In the block where you describe dose actor, you should try to write
/gate/actor/doseDistribution/saveEveryNSeconds 1 instead of
/gate/actor/doseDistribution/saveEveryNSeconds 10.

Because when you save every 10 seconds you miss all information which is
accumulated before tenth second.

Best regards,
Elnur

2017-01-05 17:32 GMT+03:00 Theodora Kostou <theokost at upatras.gr>:

> 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
>
>
>
>
>
>
>
>
> _______________________________________________
> Gate-users mailing list
> 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/20170106/ab23c700/attachment-0001.html>