<HTML><BODY>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><p>#=====================================================<br># PHYSICS<br>#=====================================================</p><p>/gate/physics/processList Available<br>/gate/physics/processList Enabled<br>/gate/physics/processList Initialized</p><p>/gate/geometry/setIonisationPotential Water 75 eV<br>/gate/geometry/setIonisationPotential Air 85.7 eV<br>/gate/geometry/setIonisationPotential PMMA 74 eV</p><p>/control/execute physicslist_EM_std.mac<br>/control/execute physicslist_HAD_std.mac</p><p>/gate/physics/init</p><p>/gate/physics/Gamma/SetCutInRegion world 1 mm<br>/gate/physics/Electron/SetCutInRegion world 1 mm<br>/gate/physics/Positron/SetCutInRegion world 1 mm</p><p>/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</p><p>#/gate/physics/SetMaxStepSizeInRegion world 1 mm<br>/gate/physics/SetMaxStepSizeInRegion Phantom 0.01 mm<br>/gate/physics/ActivateStepLimiter proton</p><p>/gate/physics/displayCuts<br>#/gate/physics/print physics.txt</p><p>#=====================================================<br># DETECTORS<br>#=====================================================</p><p>/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</p><p>#/gate/actor/addActor SimulationStatisticActor stat<br>#/gate/actor/stat/save stat-proton.txt<br>#/gate/actor/stat/saveEveryNSeconds 60<br><br></p><p># SOURCE PBS<br>/gate/source/addSource PBS PencilBeam</p><p>#activity<br>/gate/source/PBS/setActivity 1000000 becquerel</p><p>#particle type<br>/gate/source/PBS/setParticleType proton</p><p>#energy spectrum<br>/gate/source/PBS/setEnergy 62 MeV<br>/gate/source/PBS/setSigmaEnergy 0.4 MeV</p><p>#position<br>/gate/source/PBS/setPosition 0 0 -3.2 m</p><p># spot size<br>/gate/source/PBS/setSigmaX 4 mm<br>/gate/source/PBS/setSigmaY 4 mm</p><p>#divergence<br>/gate/source/PBS/setSigmaTheta 3 mrad<br>/gate/source/PBS/setSigmaPhi 3 mrad</p><p>#emittance<br>/gate/source/PBS/setEllipseXThetaEmittance 30 mm*mrad<br>/gate/source/PBS/setEllipseYPhiEmittance 30 mm*mrad</p><p>#positive-> converging beam; negative diverging beam<br>/gate/source/PBS/setEllipseXThetaRotationNorm negative<br>/gate/source/PBS/setEllipseYPhiRotationNorm negative</p><p>#beam direction rotation<br>#/gate/source/PBS/setRotationAxis 0 0 1<br>#/gate/source/PBS/setRotationAngle 180 deg</p>And here is a result of Dose Actor's working:<br><br><p>#################################### <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</p>and so on (many zeros)<br>-- <br>Best regards, Leonid</BODY></HTML>