<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:13px"><div id="yui_3_16_0_1_1419600534486_48074">Hey,</div><div id="yui_3_16_0_1_1419600534486_48192"><br></div><div dir="ltr" id="yui_3_16_0_1_1419600534486_48016"><span id="yui_3_16_0_1_1419600534486_48349">When you attach a volume to an actor, you should assure that the "setposition" is matched to your given volume (if you want to score the dose in the whole volume). If this is the case, so I recommend you just to attach your volume, set the resolution, and comment the followings: setposition, voxel size and setsize.</span></div><div id="yui_3_16_0_1_1419600534486_49704" dir="ltr"><span id="yui_3_16_0_1_1419600534486_48349">Good luck</span></div><div id="yui_3_16_0_1_1419600534486_49707" dir="ltr"><span id="yui_3_16_0_1_1419600534486_48349">- Soman</span></div> <div class="qtdSeparateBR"><br><br></div><div style="display: block;" class="yahoo_quoted"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 13px;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div dir="ltr"> <font size="2" face="Arial"> On Thursday, December 25, 2014 2:55 PM, Leonid Ivanov <leonid78.lenya@mail.ru> wrote:<br> </font> </div> <br><br> <div class="y_msg_container"><div id="yiv9377917922">
<div>Hello, Gaters! I am simulating CATANA facility (proton therapy of ocular tumours). It is necessary to find dose distribution in water phantom (box). I saw in manual and <br>examples (radiotherapy) for understanding how does it work. But in my macros adding dose actor doesn't result in something good. There is zero dose in whole phantom. I can't understand what am I doing wrong. Please give me advice how can I solve this problem. Thanks in advance.<br><br>Here is a part of my macros which is responsible for Physics, Dose Actor and Source.<br><br><div>#=====================================================<br># PHYSICS<br>#=====================================================</div><div>/gate/physics/processList Available<br>/gate/physics/processList Enabled<br>/gate/physics/processList Initialized</div><div>/gate/geometry/setIonisationPotential Water 75 eV<br>/gate/geometry/setIonisationPotential Air 85.7 eV<br>/gate/geometry/setIonisationPotential PMMA 74 eV</div><div>/control/execute physicslist_EM_std.mac<br>/control/execute physicslist_HAD_std.mac</div><div>/gate/physics/init</div><div>/gate/physics/Gamma/SetCutInRegion world 1 mm<br>/gate/physics/Electron/SetCutInRegion world 1 mm<br>/gate/physics/Positron/SetCutInRegion world 1 mm</div><div>/gate/physics/Gamma/SetCutInRegion Phantom 0.1 mm<br>/gate/physics/Electron/SetCutInRegion Phantom 0.1 mm<br>/gate/physics/Positron/SetCutInRegion Phantom 0.1 mm</div><div>#/gate/physics/SetMaxStepSizeInRegion world 1 mm<br>/gate/physics/SetMaxStepSizeInRegion Phantom 0.01 mm<br>/gate/physics/ActivateStepLimiter proton</div><div>/gate/physics/displayCuts<br>#/gate/physics/print physics.txt</div><div>#=====================================================<br># DETECTORS<br>#=====================================================</div><div>/gate/actor/addActor DoseActor doseDistribution<br>/gate/actor/doseDistribution/save IDD-proton.txt<br>/gate/actor/doseDistribution/attachTo Phantom<br>/gate/actor/doseDistribution/stepHitType random<br>/gate/actor/doseDistribution/setPosition 0 0 +2.262 m<br>#/gate/actor/doseDistribution/setResolution 1 1 400 <br>#/gate/actor/doseDistribution/setSize 10 10 10 cm<br>/gate/actor/doseDistribution/setVoxelSize 2 2 2 mm<br>/gate/actor/doseDistribution/saveEveryNSeconds 60<br>/gate/actor/doseDistribution/enableEdep false<br>/gate/actor/doseDistribution/enableUncertaintyEdep false<br>/gate/actor/doseDistribution/enableDose true<br>/gate/actor/doseDistribution/enableUncertaintyDose false<br>/gate/actor/doseDistribution/enableNumberOfHits false</div><div>#/gate/actor/addActor SimulationStatisticActor stat<br>#/gate/actor/stat/save stat-proton.txt<br>#/gate/actor/stat/saveEveryNSeconds 60<br><br></div><div># SOURCE PBS<br>/gate/source/addSource PBS PencilBeam</div><div>#activity<br>/gate/source/PBS/setActivity 1000000 becquerel</div><div>#particle type<br>/gate/source/PBS/setParticleType proton</div><div>#energy spectrum<br>/gate/source/PBS/setEnergy 62 MeV<br>/gate/source/PBS/setSigmaEnergy 0.4 MeV</div><div>#position<br>/gate/source/PBS/setPosition 0 0 -3.2 m</div><div># spot size<br>/gate/source/PBS/setSigmaX 4 mm<br>/gate/source/PBS/setSigmaY 4 mm</div><div>#divergence<br>/gate/source/PBS/setSigmaTheta 3 mrad<br>/gate/source/PBS/setSigmaPhi 3 mrad</div><div>#emittance<br>/gate/source/PBS/setEllipseXThetaEmittance 30 mm*mrad<br>/gate/source/PBS/setEllipseYPhiEmittance 30 mm*mrad</div><div>#positive-> converging beam; negative diverging beam<br>/gate/source/PBS/setEllipseXThetaRotationNorm negative<br>/gate/source/PBS/setEllipseYPhiRotationNorm negative</div><div>#beam direction rotation<br>#/gate/source/PBS/setRotationAxis 0 0 1<br>#/gate/source/PBS/setRotationAngle 180 deg</div>And here is a result of Dose Actor's working:<br><br><div>#################################### <br># Matrix Size= (600,600,500)<br># Resol = (300,300,250)<br># VoxelSize = (2,2,2)<br># nbVal = 22500000<br>#################################### <br>0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0</div>and so on (many zeros)<br>-- <br>Best regards, Leonid</div>
</div><br>_______________________________________________<br>Gate-users mailing list<br><a href="" class="removed-link" ymailto="mailto:Gate-users@lists.opengatecollaboration.org">Gate-users@lists.opengatecollaboration.org</a><br><a href="" class="removed-link" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br><br></div> </div> </div> </div> </div></body></html>