[Gate-users] Problem with Dose Actor

Alex Vergara Gil alexvergaragil at gmail.com
Mon Jan 5 17:50:58 CET 2015


I first try to make a ROI for the volume, but it runs so fast that is
not necessary at all (1e8 primaries in 4 hours with just 1 processor).
So your suggestions work fine, except now for the results, I have send
another thread with this issue, please anwser me that
Regards
Alex

2014-12-26 15:22 GMT-05:00, S. Saghamanesh <s.saghamanesh at yahoo.com>:
> 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
>
>


More information about the Gate-users mailing list