[Gate-users] Problem with Dose Actor

[Nada Fakhoury]

Hey Leonid,

If what you're trying to do is to measure the dose distribution in the  
entire phantom, you should uncomment the line for setting the  
resolution.

You can then for example set the resolution to 400 400 400, and get a  
3D dose map inside your phantom.

Would that help, or is your goal to calculate the dose in the entire  
phantom as a Sum?

Cheers,

Nada.

> Message: 1
> Date: Thu, 25 Dec 2014 14:24:23 +0300
> From: Leonid Ivanov <leonid78.lenya at mail.ru>
> To: gate-users at lists.opengatecollaboration.org
> Subject: [Gate-users] Problem with Dose Actor
> Message-ID: <1419506663.871301065 at f214.i.mail.ru>
> Content-Type: text/plain; charset="utf-8"
>
>  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