[Gate-users] Dose Actor Problem
AYCIN BORA
aycinbora at gmail.com
Mon Apr 13 18:21:37 CEST 2015
Dear Alex
There is still no output I achieved.
Changed whole physic.mac content to:
/gate/physics/addPhysicsList emstandard_opt3
Changed source.mac content to:
/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.0 cm
/gate/source/list
Finally changed the order of lines in Actor.mac file.
Thanks for your attention
Best Regards
Aycin BORA
On 13/04/2015, Alex Vergara Gil <alexvergaragil at gmail.com> wrote:
> 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