[Gate-users] Dose Actor Problem
Alex Vergara Gil
alexvergaragil at gmail.com
Mon Apr 13 18:12:10 CEST 2015
MMM
Try putting this line
/gate/actor/USMdoseDistribution/save Sphere.hdr
after the line
/gate/actor/USMdoseDistribution/enableUncertaintyDoseToWater true
Let me know what happens
Alex
2015-04-13 12:06 GMT-04:00, AYCIN BORA <aycinbora at gmail.com>:
> Dear Alex
>
> I made changes that you advised but dose actor acquired no output again.
>
> Best Regards
> Aycin BORA
>
> On 13/04/2015, Alex Vergara Gil <alexvergaragil at gmail.com> wrote:
>> First of all, change
>>
>> #PHYSICS
>>
>> #=================================================================================================
>> # Electromagnetic processes : physics list standard option 3
>> (issue de egammaStandardPhys.mac)
>> #=================================================================================================
>>
>> /gate/physics/addProcess PhotoElectric
>> /gate/physics/processes/PhotoElectric/setModel StandardModel
>>
>> /gate/physics/addProcess Compton
>> /gate/physics/processes/Compton/setModel StandardModel
>>
>> /gate/physics/addProcess RayleighScattering
>> /gate/physics/processes/RayleighScattering/setModel PenelopeModel
>>
>> /gate/physics/addProcess ElectronIonisation
>> /gate/physics/processes/ElectronIonisation/setModel StandardModel e-
>> /gate/physics/processes/ElectronIonisation/setModel StandardModel e+
>>
>> /gate/physics/addProcess Bremsstrahlung e-
>> /gate/physics/addProcess Bremsstrahlung e+
>> /gate/physics/processes/Bremsstrahlung/setModel StandardModel e-
>> /gate/physics/processes/Bremsstrahlung/setModel StandardModel e+
>>
>> /gate/physics/addProcess PositronAnnihilation e+
>>
>> /gate/physics/addProcess eMultipleScattering e+
>> /gate/physics/addProcess eMultipleScattering e-
>>
>> #=================================================================================================
>> # Additional processes not included in the list option 3
>> #=================================================================================================
>>
>> #pair production
>> /gate/physics/addProcess GammaConversion
>> /gate/physics/processes/GammaConversion/setModel StandardModel
>>
>> By
>>
>> #=====================================================
>> # Physic List Mechanism
>> #=====================================================
>>
>> /gate/physics/addPhysicsList emstandard_opt3
>>
>>
>> Second
>>
>> For fluor 18 there is already a fast ion source that should be used like
>>
>> #-------------------oooooOOOOO00000OOOOOooooo---------------------#
>> # #
>> # D E F I N I T I O N O F #
>> # Y O U R S O U R C E S #
>> # #
>> #-------------------oooooOOOOO00000OOOOOooooo---------------------#
>>
>>
>> # F A S T P R O D U C T I O N M O D U L E F O R F - 1 8 I S
>> O T O P E S
>> /gate/source/addSource F18CylinderSource
>> /gate/source/F18CylinderSource/setActivity 10000. becquerel
>> /gate/source/F18CylinderSource/gps/particle e+
>> /gate/source/F18CylinderSource/setForcedUnstableFlag true
>> /gate/source/F18CylinderSource/setForcedHalfLife 6586 s
>> /gate/source/F18CylinderSource/gps/energytype Fluor18
>> /gate/source/F18CylinderSource/gps/type Volume
>> /gate/source/F18CylinderSource/gps/shape Cylinder
>> /gate/source/F18CylinderSource/gps/radius 0.5 mm
>> /gate/source/F18CylinderSource/gps/halfz 0.5 mm
>> /gate/source/F18CylinderSource/gps/angtype iso
>> /gate/source/F18CylinderSource/gps/centre 0.0 0.0 0. cm
>>
>>
>> /gate/source/list
>>
>> Do this and tell me what happens
>>
>> Regards
>> Alex
>>
>> 2015-04-13 11:11 GMT-04:00, AYCIN BORA <aycinbora at gmail.com>:
>>> Thank you very much Alex for your attention.I can send my macro files.
>>> =====================================================
>>> # DETECTORS
>>> #=====================================================
>>>
>>> /gate/actor/addActor DoseActor USMdoseDistribution
>>> /gate/actor/USMdoseDistribution/save Sphere.hdr
>>> /gate/actor/USMdoseDistribution/attachTo unitsphere
>>> #/gate/actor/USMdoseDistribution/setSize 0.26 0.26 0.26 cm
>>> /gate/actor/USMdoseDistribution/setPosition 0 0 0 cm
>>> /gate/actor/USMdoseDistribution/setResolution 3 3 2600
>>> /gate/actor/USMdoseDistribution/stepHitType random
>>> #/gate/actor/USMdoseDistribution/saveEveryNEvents 100
>>>
>>>
>>> /gate/actor/USMdoseDistribution/enableEdep true
>>> /gate/actor/USMdoseDistribution/enableUncertaintyEdep true
>>> /gate/actor/USMdoseDistribution/enableSquaredEdep true
>>> /gate/actor/USMdoseDistribution/enableDose true
>>> /gate/actor/USMdoseDistribution/enableUncertaintyDose true
>>> /gate/actor/USMdoseDistribution/enableSquaredDose true
>>> /gate/actor/USMdoseDistribution/enableNumberOfHits true
>>> /gate/actor/USMdoseDistribution/enableDoseToWater true
>>> /gate/actor/USMdoseDistribution/enableUncertaintyDoseToWater true
>>>
>>> /gate/actor/addActor SimulationStatisticActor USMstat
>>> /gate/actor/USMstat/save USMstat.txt
>>>
>>> #GEOMETRY
>>> /gate/geometry/setMaterialDatabase GateMaterials.db
>>>
>>>
>>> # WORLD
>>> /gate/world/geometry/setXLength 0.3 cm
>>> /gate/world/geometry/setYLength 0.3 cm
>>> /gate/world/geometry/setZLength 0.3 cm
>>> /gate/world/setMaterial Air
>>> /gate/world/vis/setVisible 1
>>> /gate/world/vis/forceWireframe
>>>
>>>
>>>
>>>
>>> #SPHERE
>>> /gate/world/daughters/name unitsphere
>>> /gate/world/daughters/insert sphere
>>> /gate/unitsphere/setMaterial Water
>>> /gate/unitsphere/geometry/setRmin 0 cm
>>> /gate/unitsphere/geometry/setRmax 0.13 cm
>>> /gate/unitsphere/vis/setColor yellow
>>> /gate/unitsphere/vis/forceWireframe
>>> /gate/unitsphere/setTranslation 0. 0. 0. cm
>>>
>>> #PHYSICS
>>>
>>> #=================================================================================================
>>> # Electromagnetic processes : physics list standard option 3
>>> (issue de egammaStandardPhys.mac)
>>> #=================================================================================================
>>>
>>> /gate/physics/addProcess PhotoElectric
>>> /gate/physics/processes/PhotoElectric/setModel StandardModel
>>>
>>> /gate/physics/addProcess Compton
>>> /gate/physics/processes/Compton/setModel StandardModel
>>>
>>> /gate/physics/addProcess RayleighScattering
>>> /gate/physics/processes/RayleighScattering/setModel PenelopeModel
>>>
>>> /gate/physics/addProcess ElectronIonisation
>>> /gate/physics/processes/ElectronIonisation/setModel StandardModel e-
>>> /gate/physics/processes/ElectronIonisation/setModel StandardModel e+
>>>
>>> /gate/physics/addProcess Bremsstrahlung e-
>>> /gate/physics/addProcess Bremsstrahlung e+
>>> /gate/physics/processes/Bremsstrahlung/setModel StandardModel e-
>>> /gate/physics/processes/Bremsstrahlung/setModel StandardModel e+
>>>
>>> /gate/physics/addProcess PositronAnnihilation e+
>>>
>>> /gate/physics/addProcess eMultipleScattering e+
>>> /gate/physics/addProcess eMultipleScattering e-
>>>
>>> #=================================================================================================
>>> # Additional processes not included in the list option 3
>>> #=================================================================================================
>>>
>>> #pair production
>>> /gate/physics/addProcess GammaConversion
>>> /gate/physics/processes/GammaConversion/setModel StandardModel
>>>
>>> #==================================================================================================
>>> # customized CUTS
>>> #==================================================================================================
>>>
>>>
>>> /gate/physics/Gamma/SetCutInRegion world 1.0 mm
>>> /gate/physics/Electron/SetCutInRegion world 1.0 mm
>>> /gate/physics/Positron/SetCutInRegion world 1.0 mm
>>>
>>> /gate/physics/Gamma/SetCutInRegion unitsphere 0.1 mm
>>> /gate/physics/Electron/SetCutInRegion unitsphere 0.1 mm
>>> /gate/physics/Positron/SetCutInRegion unitsphere 0.1 mm
>>>
>>> #1/1 of dosel size
>>> /gate/physics/processes/ElectronIonisation/setStepFunction e+ 0.1 0.001
>>> mm
>>> /gate/physics/processes/ElectronIonisation/setStepFunction e- 0.1 0.001
>>> mm
>>>
>>>
>>> /gate/physics/displayCuts
>>>
>>>
>>> /gate/physics/setEMin 0.1 keV
>>> /gate/physics/setEMax 10 GeV
>>>
>>> /gate/physics/setDEDXBinning 220 $ pouvoir d'arret
>>> /gate/physics/setLambdaBinning 220 $ mean free path
>>>
>>> #===================================================================================================
>>> #additional features
>>> #===================================================================================================
>>>
>>>
>>> /gate/physics/processes/eMultipleScattering/setGeometricalStepLimiterType
>>> e- distanceToBoundary
>>> /gate/physics/processes/eMultipleScattering/setGeometricalStepLimiterType
>>> e+ distanceToBoundary
>>>
>>> #SOURCE
>>>
>>> # 18F
>>> #taken from (MIRD radionuclide data and decay schemes, 2nd edition, KF
>>> Eckerman & A Endo Eds, The Society of Nuclear Medicine, Reston, VA,
>>> SNM, 2008.)
>>>
>>> #=====================================================
>>> # Positron emission spectrum
>>> #=====================================================
>>> #Distrbution is automatically normalized to 1. Number of particle per
>>> nt must be adjusted post-treatment
>>>
>>>
>>> /gate/source/addSource positron
>>> /gate/source/positron/gps/type Volume
>>> /gate/source/positron/gps/shape Sphere
>>> /gate/source/positron/gps/radius 0.13 cm
>>> #/gate/source/positron/gps/centre 0. 0. 0. cm
>>> /gate/source/positron/gps/confine unitsphere
>>>
>>> /gate/source/positron/gps/particle e+
>>> /gate/source/positron/gps/angtype iso
>>> /gate/source/positron/gps/energytype Arb
>>> /gate/source/positron/gps/histname arb
>>>
>>>
>>> # Warning : probability density (or point wise) and not bin values
>>> /gate/source/positron/gps/histpoint 0.0001 1.781E-08
>>> /gate/source/positron/gps/histpoint 0.00011 4.220E-08
>>> /gate/source/positron/gps/histpoint 0.00012 8.955E-08
>>> /gate/source/positron/gps/histpoint 0.00013 1.739E-07
>>> /gate/source/positron/gps/histpoint 0.00014 3.142E-07
>>> /gate/source/positron/gps/histpoint 0.00015 5.345E-07
>>> /gate/source/positron/gps/histpoint 0.00016 8.643E-07
>>> /gate/source/positron/gps/histpoint 0.00018 1.998E-06
>>> /gate/source/positron/gps/histpoint 0.0002 4.060E-06
>>> /gate/source/positron/gps/histpoint 0.00022 7.472E-06
>>> /gate/source/positron/gps/histpoint 0.00024 1.272E-05
>>> /gate/source/positron/gps/histpoint 0.00026 2.034E-05
>>> /gate/source/positron/gps/histpoint 0.00028 3.090E-05
>>> /gate/source/positron/gps/histpoint 0.0003 4.499E-05
>>> /gate/source/positron/gps/histpoint 0.00032 6.322E-05
>>> /gate/source/positron/gps/histpoint 0.00036 1.145E-04
>>> /gate/source/positron/gps/histpoint 0.0004 1.893E-04
>>> /gate/source/positron/gps/histpoint 0.00045 3.230E-04
>>> /gate/source/positron/gps/histpoint 0.0005 5.090E-04
>>> /gate/source/positron/gps/histpoint 0.00055 7.551E-04
>>> /gate/source/positron/gps/histpoint 0.0006 1.068E-03
>>> /gate/source/positron/gps/histpoint 0.00065 1.455E-03
>>> /gate/source/positron/gps/histpoint 0.0007 1.920E-03
>>> /gate/source/positron/gps/histpoint 0.00075 2.468E-03
>>> /gate/source/positron/gps/histpoint 0.0008 3.102E-03
>>> /gate/source/positron/gps/histpoint 0.00085 3.824E-03
>>> /gate/source/positron/gps/histpoint 0.0009 4.634E-03
>>> /gate/source/positron/gps/histpoint 0.001 6.518E-03
>>> /gate/source/positron/gps/histpoint 0.0011 8.743E-03
>>> /gate/source/positron/gps/histpoint 0.0012 1.129E-02
>>> /gate/source/positron/gps/histpoint 0.0013 1.412E-02
>>> /gate/source/positron/gps/histpoint 0.0014 1.722E-02
>>> /gate/source/positron/gps/histpoint 0.0015 2.055E-02
>>> /gate/source/positron/gps/histpoint 0.0016 2.409E-02
>>> /gate/source/positron/gps/histpoint 0.0018 3.167E-02
>>> /gate/source/positron/gps/histpoint 0.002 3.977E-02
>>> /gate/source/positron/gps/histpoint 0.0022 4.827E-02
>>> /gate/source/positron/gps/histpoint 0.0024 5.705E-02
>>> /gate/source/positron/gps/histpoint 0.0026 6.603E-02
>>> /gate/source/positron/gps/histpoint 0.0028 7.515E-02
>>> /gate/source/positron/gps/histpoint 0.003 8.436E-02
>>> /gate/source/positron/gps/histpoint 0.0032 9.363E-02
>>> /gate/source/positron/gps/histpoint 0.0036 1.123E-01
>>> /gate/source/positron/gps/histpoint 0.004 1.309E-01
>>> /gate/source/positron/gps/histpoint 0.0045 1.540E-01
>>> /gate/source/positron/gps/histpoint 0.005 1.769E-01
>>> /gate/source/positron/gps/histpoint 0.0055 1.995E-01
>>> /gate/source/positron/gps/histpoint 0.006 2.218E-01
>>> /gate/source/positron/gps/histpoint 0.0065 2.438E-01
>>> /gate/source/positron/gps/histpoint 0.007 2.655E-01
>>> /gate/source/positron/gps/histpoint 0.0075 2.868E-01
>>> /gate/source/positron/gps/histpoint 0.008 3.078E-01
>>> /gate/source/positron/gps/histpoint 0.0085 3.285E-01
>>> /gate/source/positron/gps/histpoint 0.009 3.488E-01
>>> /gate/source/positron/gps/histpoint 0.01 3.887E-01
>>> /gate/source/positron/gps/histpoint 0.011 4.274E-01
>>> /gate/source/positron/gps/histpoint 0.012 4.650E-01
>>> /gate/source/positron/gps/histpoint 0.013 5.016E-01
>>> /gate/source/positron/gps/histpoint 0.014 5.372E-01
>>> /gate/source/positron/gps/histpoint 0.015 5.720E-01
>>> /gate/source/positron/gps/histpoint 0.016 6.059E-01
>>> /gate/source/positron/gps/histpoint 0.018 6.713E-01
>>> /gate/source/positron/gps/histpoint 0.02 7.338E-01
>>> /gate/source/positron/gps/histpoint 0.022 7.937E-01
>>> /gate/source/positron/gps/histpoint 0.024 8.513E-01
>>> /gate/source/positron/gps/histpoint 0.026 9.066E-01
>>> /gate/source/positron/gps/histpoint 0.028 9.599E-01
>>> /gate/source/positron/gps/histpoint 0.03 1.011E+00
>>> /gate/source/positron/gps/histpoint 0.032 1.061E+00
>>> /gate/source/positron/gps/histpoint 0.036 1.156E+00
>>> /gate/source/positron/gps/histpoint 0.04 1.245E+00
>>> /gate/source/positron/gps/histpoint 0.045 1.349E+00
>>> /gate/source/positron/gps/histpoint 0.05 1.446E+00
>>> /gate/source/positron/gps/histpoint 0.055 1.537E+00
>>> /gate/source/positron/gps/histpoint 0.06 1.622E+00
>>> /gate/source/positron/gps/histpoint 0.065 1.702E+00
>>> /gate/source/positron/gps/histpoint 0.07 1.777E+00
>>> /gate/source/positron/gps/histpoint 0.075 1.848E+00
>>> /gate/source/positron/gps/histpoint 0.08 1.914E+00
>>> /gate/source/positron/gps/histpoint 0.085 1.976E+00
>>> /gate/source/positron/gps/histpoint 0.09 2.035E+00
>>> /gate/source/positron/gps/histpoint 0.1 2.142E+00
>>> /gate/source/positron/gps/histpoint 0.11 2.236E+00
>>> /gate/source/positron/gps/histpoint 0.12 2.318E+00
>>> /gate/source/positron/gps/histpoint 0.13 2.389E+00
>>> /gate/source/positron/gps/histpoint 0.14 2.449E+00
>>> /gate/source/positron/gps/histpoint 0.15 2.499E+00
>>> /gate/source/positron/gps/histpoint 0.16 2.540E+00
>>> /gate/source/positron/gps/histpoint 0.18 2.594E+00
>>> /gate/source/positron/gps/histpoint 0.2 2.616E+00
>>> /gate/source/positron/gps/histpoint 0.22 2.607E+00
>>> /gate/source/positron/gps/histpoint 0.24 2.570E+00
>>> /gate/source/positron/gps/histpoint 0.26 2.507E+00
>>> /gate/source/positron/gps/histpoint 0.28 2.421E+00
>>> /gate/source/positron/gps/histpoint 0.3 2.314E+00
>>> /gate/source/positron/gps/histpoint 0.32 2.188E+00
>>> /gate/source/positron/gps/histpoint 0.36 1.891E+00
>>> /gate/source/positron/gps/histpoint 0.4 1.549E+00
>>> /gate/source/positron/gps/histpoint 0.45 1.094E+00
>>> /gate/source/positron/gps/histpoint 0.5 6.557E-01
>>> /gate/source/positron/gps/histpoint 0.55 2.879E-01
>>> /gate/source/positron/gps/histpoint 0.6 5.162E-02
>>> /gate/source/positron/gps/histpoint 0.6335 0.000E+00
>>> /gate/source/positron/gps/arbint Spline
>>>
>>>
>>> /gate/source/list
>>>
>>>
>>> On 13/04/2015, Alex Vergara Gil <alexvergaragil at gmail.com> wrote:
>>>> Please Start with an example that is related to your needs, learn from
>>>> it and then adapt it according to your necessities. Another thing is
>>>> that you send a minimal example of your macro to see what is wrong.
>>>>
>>>> Regards
>>>> Alex
>>>>
>>>> 2015-04-13 10:57 GMT-04:00, AYCIN BORA <aycinbora at gmail.com>:
>>>>> Dear Gate Users;
>>>>>
>>>>> I'm new at Gate.I'm trying to set up a simulation about radiotheraphy
>>>>> application.I confine my source distribution to my volume-unitvolume-
>>>>> and attach the dose actor to that volume. But I get no output. Would
>>>>> you please help me?
>>>>>
>>>>> Kind Regards
>>>>>
>>>>> Aycin BORA
>>>>> _______________________________________________
>>>>> 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