[Gate-users] Dose Actor Problem
AYCIN BORA
aycinbora at gmail.com
Mon Apr 13 18:06:25 CEST 2015
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