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

[Dimitris Anestis]

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/970bd3ca/attachment.html>