[Gate-users] Problem with Dose Actor

[S. Saghamanesh]

Hey,
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.Good luck- Soman 

     On Thursday, December 25, 2014 2:55 PM, Leonid Ivanov <leonid78.lenya at mail.ru> wrote:
   

 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 degAnd 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 0and so on (many zeros)
-- 
Best regards, Leonid
_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users

   
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20141226/2e88c7ff/attachment.html>