[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