[Gate-users] Dose Actor Problem
Alex Vergara Gil
alexvergaragil at gmail.com
Mon Apr 13 17:29:42 CEST 2015
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