<html><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8" /></head><body style='font-size: 10pt'>
<p>Hello, </p>
<p>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. </p>
<p>From the code below I can visualize the rotation of the xcat phantom but I take zero counts as output. </p>
<p><br />/control/execute mac/verbose.mac<br />#/control/execute mac/visu.mac<br />/gate/geometry/setMaterialDatabase data/GateMaterials.db<br />/gate/world/setMaterial Vacuum<br />/gate/world/geometry/setXLength 1000. cm<br />/gate/world/geometry/setYLength 1000. cm<br />/gate/world/geometry/setZLength 1000. cm<br /># Container to rotate the phantom<br />/gate/world/daughters/name CT<br />/gate/world/daughters/insert box<br />/gate/CT/setMaterial Vacuum<br />/gate/CT/geometry/setXLength 640. mm<br />/gate/CT/geometry/setYLength 640. mm<br />/gate/CT/geometry/setZLength 1825. mm<br />/gate/CT/vis/forceWireframe<br />/gate/CT/moves/insert rotation<br />/gate/CT/rotation/setSpeed 45 deg/s<br />/gate/CT/rotation/setAxis 0 0 1<br /># Patient phantom (128x128x128 voxels, 5 mm voxelize)<br />/gate/CT/daughters/name xcat<br />/gate/CT/daughters/insert ImageRegularParametrisedVolume<br />/gate/xcat/geometry/setImage data/xcat.h33<br />/gate/xcat/geometry/setRangeToMaterialFile data/Range.dat<br />/gate/xcat/geometry/TranslateTheImageAtThisIsoCenter 0 0 0 mm<br />/gate/xcat/placement/setTranslation -320 -320 -912.5 mm<br />#/gate/xcat/placement/setRotationAxis 1 0 0<br />#/gate/xcat/placement/setRotationAngle 0 deg<br />/gate/xcat/vis/forceWireframe<br />#/gate/xcat/setSkipEqualMaterials 0<br />/gate/CT/attachPhantomSD<br />/gate/xcat/attachPhantomSD</p>
<p>/gate/physics/addProcess PhotoElectric<br />/gate/physics/processes/PhotoElectric/setModel StandardModel<br />/gate/physics/addProcess Compton<br />/gate/physics/processes/Compton/setModel StandardModel<br />/gate/physics/addProcess RayleighScattering<br />/gate/physics/processes/RayleighScattering/setModel PenelopeModel<br />/gate/physics/addProcess ElectronIonisation<br />/gate/physics/processes/ElectronIonisation/setModel StandardModel e-<br />/gate/physics/addProcess Bremsstrahlung<br />/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-<br />/gate/physics/addProcess MultipleScattering e-<br />/gate/physics/processList Enabled<br />/gate/physics/processList Initialized</p>
<p>/gate/actor/addActor DoseActor doseDistribution<br />/gate/actor/doseDistribution/save output/myactor.mhd<br />/gate/actor/doseDistribution/attachTo xcat<br />/gate/actor/doseDistribution/stepHitType random<br />/gate/actor/doseDistribution/setPosition 0 0 0 mm<br />/gate/actor/doseDistribution/setVoxelSize 5 5 5 mm<br />/gate/actor/doseDistribution/enableEdep true<br />/gate/actor/doseDistribution/enableUncertaintyEdep true<br />/gate/actor/doseDistribution/enableDose true<br />/gate/actor/doseDistribution/enableSquaredEdep false<br />/gate/actor/doseDistribution/enableNumberOfHits false<br />/gate/actor/doseDistribution/saveEveryNSeconds 10</p>
<p>/gate/source/addSource xraygun gps<br />/gate/source/xraygun/gps/verbose 1<br />/gate/source/xraygun/gps/particle gamma<br />/gate/source/xraygun/gps/energytype Arb<br />/gate/source/xraygun/gps/histname arb<br />/gate/source/xraygun/gps/emin 12.00 keV<br />/gate/source/xraygun/gps/emax 120.00 keV<br />/control/execute mac/histogram.mac <br />/gate/source/xraygun/gps/arbint Lin<br />/gate/source/xraygun/gps/type Plane<br />/gate/source/xraygun/gps/shape Rectangle<br />/gate/source/xraygun/gps/halfx 0.9 mm<br />/gate/source/xraygun/gps/halfy 0.8 mm<br />/gate/source/xraygun/gps/centre -60. 0. 0. cm<br />/gate/source/xraygun/gps/angtype iso<br />/gate/source/xraygun/gps/pos/rot1 1 0 0<br />/gate/source/xraygun/gps/pos/rot2 0 1 0<br />#cone angle <br />/gate/source/xraygun/gps/mintheta 88.1 deg<br />/gate/source/xraygun/gps/maxtheta 91.9 deg<br />#fan angle <br />/gate/source/xraygun/gps/minphi 155. deg<br />/gate/source/xraygun/gps/maxphi 205.0 deg<br />/gate/source/list</p>
<p>/gate/run/initialize<br />/gate/random/setEngineName MersenneTwister<br />/gate/random/setEngineSeed auto <br />/gate/application/setTimeStart 0 s<br />/gate/application/setTimeStop 10. s<br />/gate/application/setTimeSlice 1. s<br />/gate/application/setTotalNumberOfPrimaries 10<br />/gate/application/startDAQ</p>
<p>Is something wrong or missing in the code? </p>
<p>I would appreciate any kind of advice. </p>
<p>Thank you </p>
<p>Theodora</p>
<p> </p>
<p> </p>
<p> </p>
<div>
<div class="pre" style="margin: 0; padding: 0; font-family: monospace"> </div>
</div>
</body></html>