[Gate-users] Dose actor
Mohammed REZZOUG
nmedrezz at gmail.com
Mon Feb 10 16:04:22 CET 2020
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
>
>
>
-------------- section suivante --------------
Une pièce jointe HTML a été nettoyée...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200210/2de421bb/attachment-0001.html>
More information about the Gate-users
mailing list