<div dir="ltr">In the previous letter I forgot to write one more point:<div>In the block where you describe dose actor, you should try to write</div><div><span style="font-size:13.3333px">/gate/actor/doseDistribution/</span><wbr style="font-size:13.3333px"><span style="font-size:13.3333px">saveEveryNSeconds 1 instead of</span><br></div><div><span style="font-size:13.3333px">/gate/actor/doseDistribution/</span><wbr style="font-size:13.3333px"><span style="font-size:13.3333px">saveEveryNSeconds 10.</span><span style="font-size:13.3333px"><br></span></div><div><span style="font-size:13.3333px"><br></span></div><div><span style="font-size:13.3333px">Because when you save every 10 seconds you miss all information which is accumulated before tenth second.</span></div><div><span style="font-size:13.3333px"><br></span></div><div class="gmail_extra">Best regards,</div><div class="gmail_extra">Elnur<br>
<br><div class="gmail_quote">2017-01-05 17:32 GMT+03:00 Theodora Kostou <span dir="ltr"><<a href="mailto:theokost@upatras.gr" target="_blank">theokost@upatras.gr</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div 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/<wbr>setMaterialDatabase data/GateMaterials.db<br>/gate/world/setMaterial Vacuum<br>/gate/world/geometry/<wbr>setXLength 1000. cm<br>/gate/world/geometry/<wbr>setYLength 1000. cm<br>/gate/world/geometry/<wbr>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/<wbr>setRangeToMaterialFile data/Range.dat<br>/gate/xcat/geometry/<wbr>TranslateTheImageAtThisIsoCent<wbr>er 0 0 0 mm<br>/gate/xcat/placement/<wbr>setTranslation -320 -320 -912.5 mm<br>#/gate/xcat/placement/<wbr>setRotationAxis 1 0 0<br>#/gate/xcat/placement/<wbr>setRotationAngle 0 deg<br>/gate/xcat/vis/forceWireframe<br>#/gate/xcat/<wbr>setSkipEqualMaterials 0<br>/gate/CT/attachPhantomSD<br>/gate/xcat/attachPhantomSD</p>
<p>/gate/physics/addProcess PhotoElectric<br>/gate/physics/processes/<wbr>PhotoElectric/setModel StandardModel<br>/gate/physics/addProcess Compton<br>/gate/physics/processes/<wbr>Compton/setModel StandardModel<br>/gate/physics/addProcess RayleighScattering<br>/gate/physics/processes/<wbr>RayleighScattering/setModel PenelopeModel<br>/gate/physics/addProcess ElectronIonisation<br>/gate/physics/processes/<wbr>ElectronIonisation/setModel StandardModel e-<br>/gate/physics/addProcess Bremsstrahlung<br>/gate/physics/processes/<wbr>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/<wbr>save output/myactor.mhd<br>/gate/actor/doseDistribution/<wbr>attachTo xcat<br>/gate/actor/doseDistribution/<wbr>stepHitType random<br>/gate/actor/doseDistribution/<wbr>setPosition 0 0 0 mm<br>/gate/actor/doseDistribution/<wbr>setVoxelSize 5 5 5 mm<br>/gate/actor/doseDistribution/<wbr>enableEdep true<br>/gate/actor/doseDistribution/<wbr>enableUncertaintyEdep true<br>/gate/actor/doseDistribution/<wbr>enableDose true<br>/gate/actor/doseDistribution/<wbr>enableSquaredEdep false<br>/gate/actor/doseDistribution/<wbr>enableNumberOfHits false<br>/gate/actor/doseDistribution/<wbr>saveEveryNSeconds 10</p>
<p>/gate/source/addSource xraygun gps<br>/gate/source/xraygun/gps/<wbr>verbose 1<br>/gate/source/xraygun/gps/<wbr>particle gamma<br>/gate/source/xraygun/gps/<wbr>energytype Arb<br>/gate/source/xraygun/gps/<wbr>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/<wbr>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/<wbr>centre -60. 0. 0. cm<br>/gate/source/xraygun/gps/<wbr>angtype iso<br>/gate/source/xraygun/gps/pos/<wbr>rot1 1 0 0<br>/gate/source/xraygun/gps/pos/<wbr>rot2 0 1 0<br>#cone angle <br>/gate/source/xraygun/gps/<wbr>mintheta 88.1 deg<br>/gate/source/xraygun/gps/<wbr>maxtheta 91.9 deg<br>#fan angle <br>/gate/source/xraygun/gps/<wbr>minphi 155. deg<br>/gate/source/xraygun/gps/<wbr>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/<wbr>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><span class="HOEnZb"><font color="#888888">
<p>Theodora</p>
<p> </p>
<p> </p>
<p> </p>
<div>
<div class="m_7631602644511425504pre" style="margin:0;padding:0;font-family:monospace"> </div>
</div>
</font></span></div>
<br>______________________________<wbr>_________________<br>
Gate-users mailing list<br>
<a href="mailto:Gate-users@lists.opengatecollaboration.org">Gate-users@lists.<wbr>opengatecollaboration.org</a><br>
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" rel="noreferrer" target="_blank">http://lists.<wbr>opengatecollaboration.org/<wbr>mailman/listinfo/gate-users</a><br></blockquote></div><br></div></div>