[Gate-users] How to store the stopping position of positrons of 511keV

[Mahdjoub Hamdi]

Hi Dimitris,

Have you tried the production threshold? I think you can stop some events, I.e., positions, based on their kinetic energy or the distance they are transported in the medium.

You can see how to do it here:
https://opengate.readthedocs.io/en/latest/cut_and_variance_reduction_technics.html

Thanks,
Mahdjoub

Get Outlook for iOS<https://aka.ms/o0ukef>
________________________________
From: Gate-users <gate-users-bounces at lists.opengatecollaboration.org> on behalf of Dimitris Anestis <janest755 at gmail.com>
Sent: Tuesday, April 12, 2022 4:42:25 AM
To: gate-users at lists.opengatecollaboration.org <gate-users at lists.opengatecollaboration.org>
Subject: [Gate-users] How to store the stopping position of positrons of 511keV

Hello Gate users,

I am trying to set up a simulation in order to produce an annihilation image. I am using a voxelized , voxelized source (Xcat) , and the ProductionAndStoppingActor actor with particleFilter for positrons in order to store the location of the positrons 3d images. However, some positrons stored in the image are not going to be annihilated so I am looking for a way to further filter the positrons, in terms of energy. I have already added an energy filter but i t seems it is not working properly.


#################################################### Actors ####################################
/gate/actor/addActor  SimulationStatisticActor stat
/gate/actor/stat/save output/stats.txt

/gate/actor/addActor ProductionAndStoppingActor range_actor_a
/gate/actor/range_actor_a/attachTo Xcat
/gate/actor/range_actor_a/setResolution                    256 256 101
/gate/actor/range_actor_a/stepHitType                      post
/gate/actor/range_actor_a/save  output/Xcat256_256_101_0.1_no_filter.mhd

/gate/actor/addActor ProductionAndStoppingActor range_actor
/gate/actor/range_actor/attachTo Xcat
/gate/actor/range_actor/addFilter                        particleFilter
/gate/actor/range_actor/particleFilter/addParticle       e+ # beta
/gate/actor/range_actor/addFilter                        energyFilter
/gate/actor/range_actor/energyFilter/setEmin             0.0 MeV
/gate/actor/range_actor/energyFilter/setEmax             0.0 MeV
/gate/actor/range_actor/setResolution                    256 256 101
/gate/actor/range_actor/stepHitType                      post
/gate/actor/range_actor/save  output/Xcat256_256_101_0.1_filter_positrons.mhd

/gate/actor/addActor ProductionAndStoppingActor range_actor_b
/gate/actor/range_actor_b/attachTo Xcat
/gate/actor/range_actor_b/addFilter                        particleFilter
/gate/actor/range_actor_b/particleFilter/addParticle       e+ # beta
/gate/actor/range_actor_b/setResolution                    256 256 101
/gate/actor/range_actor_b/stepHitType                      post
/gate/actor/range_actor_b/save  output/Xcat256_256_101_0.1_positrons.mhd

/gate/actor/addActor ProductionAndStoppingActor range_actor_c
/gate/actor/range_actor_c/attachTo Xcat
/gate/actor/range_actor_c/addFilter                        particleFilter
/gate/actor/range_actor_c/particleFilter/addParticle       gamma # beta
/gate/actor/range_actor_c/addFilter                        energyFilter
/gate/actor/range_actor_c/energyFilter/setEmin             0.511 MeV
/gate/actor/range_actor_c/energyFilter/setEmax             0.511 MeV
/gate/actor/range_actor_c/setResolution                    256 256 101
/gate/actor/range_actor_c/stepHitType                      post
/gate/actor/range_actor_c/save  output/Xcat256_256_101_0.1_filter_gammas.mhd

/gate/actor/addActor ProductionAndStoppingActor range_actor_d
/gate/actor/range_actor_d/attachTo Xcat
/gate/actor/range_actor_d/addFilter                        particleFilter
/gate/actor/range_actor_d/particleFilter/addParticle       gamma # beta
/gate/actor/range_actor_d/setResolution                    256 256 101
/gate/actor/range_actor_d/stepHitType                      post
/gate/actor/range_actor_d/save  output/Xcat256_256_101_0.1_gammas.mhd


/gate/physics/addPhysicsList emstandard

/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 LivermoreModel

/gate/physics/addProcess ElectronIonisation
/gate/physics/processes/ElectronIonisation/setModel StandardModel e-
/gate/physics/processes/ElectronIonisation/setModel StandardModel e+

/gate/physics/addProcess Bremsstrahlung
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e+

/gate/physics/addProcess PositronAnnihilation e+

#/gate/physics/addProcess MultipleScattering e-
#/gate/physics/addProcess MultipleScattering e+
/gate/physics/addProcess RadioactiveDecay
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20220412/df562d56/attachment-0001.html>