[Gate-users] Dose actor
Dimitrios Thanasas
dthanas at uoa.gr
Mon Feb 10 16:36:05 CET 2020
Yes I think +18.5 is near the top surface.
So you are right.
Dimitris
> Thank you very much.
> my source is placed at 0 0 100 cm
> the waterbox is placed at 0 0 -20 cm
> and i want to calculate the profil to 1.5 cm depth in water
> the code become:
> /gate/actor/addActor DoseActor doseprofile
>> /gate/actor/doseprofile/save
>> output/electron-profile.txt
>> /gate/actor/doseprofile/attachTo waterbox
>> /gate/actor/doseprofile/stepHitType random
> /gate/actor/doseprofile/setPosition 0 0 -18.5 cm
> /gate/actor/doseprofile/setSize 40 1 1 cm
> /gate/actor/doseprofile/setVoxelSize 0.4 1 1 cm
> --
> Am I right?
>
> Le lun. 10 févr. 2020 à 13:56, Dimitrios Thanasas <dthanas at uoa.gr> a écrit
> :
>
>> Hi Mohamed,
>>
>> One way to calculate profile using dose actor is the following:
>>
>> > /gate/actor/addActor DoseActor doseprofile
>> > /gate/actor/doseprofile/save
>> output/electron-profile.txt
>> > /gate/actor/doseprofile/attachTo waterbox
>> > /gate/actor/doseprofile/stepHitType random
>> > /gate/actor/doseprofile/setPosition 0 0 18.5 cm (be sure
>> that
>> +18.5 is on the side of the source. You may need to put -18.5 cm)
>> > #/gate/actor/doseprofile/setResolution 100 1 1 (do not use
>> setResolution. I put this command in comment.)
>> > /gate/actor/doseprofile/setSize 40 1 1 cm (this is a
>> line
>> along x axis)
>> > /gate/actor/doseprofile/setVoxelSize 0.4 1 1 cm (0.4x1x1 cm
>> each voxel, along the line in order to have 100 voxels)
>> > /gate/actor/doseprofile/saveEveryNSeconds 60
>> > /gate/actor/doseprofile/enableEdep true
>> > /gate/actor/doseprofile/enableUncertaintyEdep true
>> > /gate/actor/doseprofile/enableDose false
>> > /gate/actor/doseprofile/enableUncertaintyDose false
>> > /gate/actor/doseprofile/enableNumberOfHits false
>>
>> I think this will help you.
>>
>> Dimitris Thanasas
>>
>> > Hi Mohammed,
>> >
>> >
>> > Please keep the conversation on the mailing list.
>> >
>> > The next few days I'll be very busy and won't have time to answer
>> myself.
>> > Maybe someone else can have a look at your macro and tell if it's fine
>> for
>> > your purpose.
>> >
>> >
>> > Good luck!
>> >
>> > /David
>> >
>> >
>> > ________________________________
>> > Von: Mohammed REZZOUG <nmedrezz at gmail.com>
>> > Gesendet: Montag, 10. Februar 2020 10:39:25
>> > An: David Boersma
>> > Betreff: Re: [Gate-users] Dose actor
>> >
>> > Hi david
>> > thank you for your reply
>> > i want to validate a phase space source ( condition DSP= 100 cm
>> waterbox
>> > 40*40*40 cm)
>> > in this order, I use this cod macro.
>> >
>> --------------------------------------------------------------------------
>> >
>> > # WORLD
>> > /gate/world/setMaterial Air
>> > /gate/world/geometry/setXLength 3.0 m
>> > /gate/world/geometry/setYLength 3.0 m
>> > /gate/world/geometry/setZLength 3.0 m
>> > /gate/world/vis/setVisible 1
>> >
>> > # WATER BOX
>> >
>> > /gate/world/daughters/name waterbox
>> > /gate/world/daughters/insert box
>> > /gate/waterbox/setMaterial Water
>> > /gate/waterbox/placement/setTranslation 0 0 -200 mm
>> > /gate/waterbox/geometry/setXLength 400 mm
>> > /gate/waterbox/geometry/setYLength 400 mm
>> > /gate/waterbox/geometry/setZLength 400 mm
>> > /gate/waterbox/vis/setColor blue
>> >
>> >
>> >
>> >
>> >
>> >
>> > # VIRTUAL PLANE FOR PHASE SPACE
>> > # It should be the same geometry than PhS generation Here we add a
>> > # "gantry" rotation. To do that easily, we define two box, on for the
>> > # rotation (according to the isocenter=world origin), then one another
>> > # to translation the phase-space plane to the right distance
>> >
>> > /gate/world/daughters/name linacGantry
>> > /gate/world/daughters/insert box
>> > /gate/linacGantry/setMaterial Air
>> > /gate/linacGantry/geometry/setXLength 1.0 mm
>> > /gate/linacGantry/geometry/setYLength 1.0 mm
>> > /gate/linacGantry/geometry/setZLength 1.0 mm
>> > /gate/linacGantry/placement/setRotationAxis 0 1 0
>> > /gate/linacGantry/placement/setRotationAngle 0 deg
>> > /gate/linacGantry/vis/setVisible 1
>> >
>> > /gate/linacGantry/daughters/name linac
>> > /gate/linacGantry/daughters/insert box
>> > /gate/linac/setMaterial Air
>> > /gate/linac/geometry/setXLength 1.0 mm
>> > /gate/linac/geometry/setYLength 1.0 mm
>> > /gate/linac/geometry/setZLength 1.0 mm
>> > /gate/linac/placement/setTranslation 0 0 1000 mm
>> > /gate/linac/vis/setVisible 1
>> >
>> > # /gate/geometry/setMaterialDatabase data/PreciseMaterials.db
>> > # /control/execute mac/linac_head.mac
>> > /gate/linac/placement/setTranslation 0 0 1000 mm
>> >
>> > /gate/linac/daughters/name PhS_Plane
>> > /gate/linac/daughters/insert cylinder
>> > /gate/PhS_Plane/setMaterial Vacuum
>> > /gate/PhS_Plane/geometry/setRmin 0 mm
>> > /gate/PhS_Plane/geometry/setRmax 70 mm
>> > /gate/PhS_Plane/geometry/setHeight 1 nm
>> > /gate/PhS_Plane/placement/setTranslation 0.0 0.0 -297 mm
>> > /gate/PhS_Plane/vis/setColor white
>> > /gate/PhS_Plane/vis/setVisible 1
>> >
>> > # =====================================================
>> > # PHYSICS
>> > # =====================================================
>> >
>> > /gate/physics/addPhysicsList emstandard_opt1
>> >
>> > /gate/physics/Gamma/SetCutInRegion world 1 mm
>> > /gate/physics/Electron/SetCutInRegion world 1 mm
>> > /gate/physics/Positron/SetCutInRegion world 1 mm
>> >
>> > /gate/physics/Gamma/SetCutInRegion waterbox 0.5 mm
>> > /gate/physics/Electron/SetCutInRegion waterbox 0.5 mm
>> > /gate/physics/Positron/SetCutInRegion waterbox 0.5 mm
>> > ------------------------------------------------------------------------
>> >
>> >
>> >
>> >
>> > I have to calculate the profile dose at equal depth 1.5 cm
>> >
>> > /gate/actor/addActor DoseActor doseprofile
>> > /gate/actor/doseprofile/save
>> output/electron-profile.txt
>> > /gate/actor/doseprofile/attachTo waterbox
>> > /gate/actor/doseprofile/stepHitType random
>> > /gate/actor/doseprofile/setPosition 0 0 18.5 cm
>> > /gate/actor/doseprofile/setResolution 100 1 1
>> > /gate/actor/doseprofile/setSize 40 40 40 cm
>> > /gate/actor/doseprofile/saveEveryNSeconds 60
>> > /gate/actor/doseprofile/enableEdep true
>> > /gate/actor/doseprofile/enableUncertaintyEdep true
>> > /gate/actor/doseprofile/enableDose false
>> > /gate/actor/doseprofile/enableUncertaintyDose false
>> > /gate/actor/doseprofile/enableNumberOfHits false
>> >
>> >
>> > --
>> >
>> > This method is valid or not?(I have to calculate the profile dose at
>> equal
>> > depth 1.5 cm)
>> >
>> > --
>> >
>> >
>> >
>> >
>> > #=====================================================
>> > # INITIALISATION
>> > #=====================================================
>> >
>> > /gate/run/initialize
>> >
>> > #=====================================================
>> > # BEAMS
>> > #=====================================================
>> >
>> > # Read from Phase Space files. Here we indicate the type of particle
>> > # (e-) because this information is not stored in the PhS to gain
>> > # space. Moreover, we set 'useRandomSymmetry' to generate more
>> > # particle than the number of particle in the PhS, using a Z axis
>> > # symmetry.
>> >
>> > /gate/source/addSource beam_g phaseSpace
>> > /gate/source/beam_g/addPhaseSpaceFile
>> output/output-PhS-g.root
>> > /gate/source/beam_g/setParticleType e-
>> > /gate/source/beam_g/attachTo PhS_Plane
>> > /gate/source/beam_g/useRandomSymmetry true
>> >
>> >
>> > #=====================================================
>> > # START BEAMS
>> > #=====================================================
>> >
>> > /gate/random/setEngineName MersenneTwister
>> > /gate/random/setEngineSeed auto
>> > /gate/application/noGlobalOutput
>> >
>> > #/tracking/verbose 2
>> >
>> > /gate/application/setTotalNumberOfPrimaries 5000
>> > /gate/application/start
>> >
>> >
>> > Thank you for any help
>> >
>> > Mohammed.
>> >
>> >
>> > Le lun. 10 févr. 2020 à 03:30, David Boersma
>> > <david.boersma at acmit.at<mailto:david.boersma at acmit.at>> a écrit :
>> >
>> > Hi Mohammed,
>> >
>> >
>> > Your question is a bit short and unspecific, we need more information
>> in
>> > order to have a chance to give you a useful answer. Can you describe
>> what
>> > you have tried so far to configure the dose actor in your simulation,
>> and
>> > what are the problems you encountered? Which Gate release are you
>> using?
>> >
>> >
>> > The documentation of the dose actor is here:
>> >
>> >
>> https://opengate.readthedocs.io/en/latest/tools_to_interact_with_the_simulation_actors.html#dose-measurement-doseactor
>> >
>> >
>> > Examples for how to use the dose actor can be found in the GateContrib
>> > project, for instance here:
>> >
>> >
>> https://github.com/OpenGATE/GateContrib/tree/master/dosimetry/dosimetry/electron_radiotherapy
>> >
>> >
>> > Hope that helps,
>> >
>> > David
>> >
>> > ________________________________
>> > Von: Gate-users
>> > <gate-users-bounces at lists.opengatecollaboration.org<mailto:
>> gate-users-bounces at lists.opengatecollaboration.org>>
>> > im Auftrag von Mohammed REZZOUG
>> > <nmedrezz at gmail.com<mailto:nmedrezz at gmail.com>>
>> > Gesendet: Sonntag, 9. Februar 2020 19:53:53
>> > An:
>> > gate-users at lists.opengatecollaboration.org<mailto:
>> gate-users at lists.opengatecollaboration.org>
>> > Betreff: [Gate-users] Dose actor
>> >
>> > I want to simulate the electron accelerator on Gate in order to
>> calculate
>> > the pdd and the dose profile in a water tank, I want to know how to
>> use
>> > the "Dose actor"
>> > Thanks for any help
>> > Cordially.
>> > Mohammed
>> >
>> > _______________________________________________
>> > 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