[Gate-users] Dose actor
Dimitrios Thanasas
dthanas at uoa.gr
Mon Feb 10 21:45:55 CET 2020
No I think that as waterbox is from -40 to 0 cm and the position of the
source is +100 cm, the position +18.5 cm of the dose actor is the right
one (gate/actor/doseprofile/setPosition 0 0 +18.5 cm).
Dimitris
> I think for 1.5 cm in water /gate/actor/doseprofile/setPosition
> 0 0 -1,5 cm
> no ??
>
> Le lun. 10 févr. 2020 à 16:36, Dimitrios Thanasas <dthanas at uoa.gr> a écrit
> :
>
>> 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