[Gate-users] Problem with Dose Actor

[Leonid Ivanov]

 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 
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.

Here is a part of my macros which is responsible for Physics, Dose Actor and Source.

#=====================================================
# PHYSICS
#=====================================================
/gate/physics/processList Available
/gate/physics/processList Enabled
/gate/physics/processList Initialized
/gate/geometry/setIonisationPotential Water 75 eV
/gate/geometry/setIonisationPotential Air 85.7 eV
/gate/geometry/setIonisationPotential PMMA 74 eV
/control/execute physicslist_EM_std.mac
/control/execute physicslist_HAD_std.mac
/gate/physics/init
/gate/physics/Gamma/SetCutInRegion world 1 mm
/gate/physics/Electron/SetCutInRegion world 1 mm
/gate/physics/Positron/SetCutInRegion world 1 mm
/gate/physics/Gamma/SetCutInRegion Phantom 0.1 mm
/gate/physics/Electron/SetCutInRegion Phantom 0.1 mm
/gate/physics/Positron/SetCutInRegion Phantom 0.1 mm
#/gate/physics/SetMaxStepSizeInRegion world 1 mm
/gate/physics/SetMaxStepSizeInRegion Phantom 0.01 mm
/gate/physics/ActivateStepLimiter proton
/gate/physics/displayCuts
#/gate/physics/print physics.txt
#=====================================================
# DETECTORS
#=====================================================
/gate/actor/addActor DoseActor doseDistribution
/gate/actor/doseDistribution/save IDD-proton.txt
/gate/actor/doseDistribution/attachTo Phantom
/gate/actor/doseDistribution/stepHitType random
/gate/actor/doseDistribution/setPosition 0 0 +2.262 m
#/gate/actor/doseDistribution/setResolution 1 1 400 
#/gate/actor/doseDistribution/setSize 10 10 10 cm
/gate/actor/doseDistribution/setVoxelSize 2 2 2 mm
/gate/actor/doseDistribution/saveEveryNSeconds 60
/gate/actor/doseDistribution/enableEdep false
/gate/actor/doseDistribution/enableUncertaintyEdep false
/gate/actor/doseDistribution/enableDose true
/gate/actor/doseDistribution/enableUncertaintyDose false
/gate/actor/doseDistribution/enableNumberOfHits false
#/gate/actor/addActor SimulationStatisticActor stat
#/gate/actor/stat/save stat-proton.txt
#/gate/actor/stat/saveEveryNSeconds 60

# SOURCE PBS
/gate/source/addSource PBS PencilBeam
#activity
/gate/source/PBS/setActivity 1000000 becquerel
#particle type
/gate/source/PBS/setParticleType proton
#energy spectrum
/gate/source/PBS/setEnergy 62 MeV
/gate/source/PBS/setSigmaEnergy 0.4 MeV
#position
/gate/source/PBS/setPosition 0 0 -3.2 m
# spot size
/gate/source/PBS/setSigmaX 4 mm
/gate/source/PBS/setSigmaY 4 mm
#divergence
/gate/source/PBS/setSigmaTheta 3 mrad
/gate/source/PBS/setSigmaPhi 3 mrad
#emittance
/gate/source/PBS/setEllipseXThetaEmittance 30 mm*mrad
/gate/source/PBS/setEllipseYPhiEmittance 30 mm*mrad
#positive-> converging beam; negative diverging beam
/gate/source/PBS/setEllipseXThetaRotationNorm negative
/gate/source/PBS/setEllipseYPhiRotationNorm negative
#beam direction rotation
#/gate/source/PBS/setRotationAxis 0 0 1
#/gate/source/PBS/setRotationAngle 180 deg And here is a result of Dose Actor's working:

#################################### 
# Matrix Size= (600,600,500)
# Resol = (300,300,250)
# VoxelSize = (2,2,2)
# nbVal = 22500000
#################################### 
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 and so on (many zeros)
-- 
Best regards, Leonid
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20141225/880b964d/attachment.html>