<div dir="ltr"><div dir="auto"><div dir="ltr"><div>Hello Gate users,<br></div><div><br></div><div>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.</div><div><br></div><div><br></div><div>#################################################### Actors ####################################<br>/gate/actor/addActor  SimulationStatisticActor stat<br>/gate/actor/stat/save output/stats.txt<br><br>/gate/actor/addActor ProductionAndStoppingActor range_actor_a<br>/gate/actor/range_actor_a/attachTo Xcat<br>/gate/actor/range_actor_a/setResolution                    256 256 101<br>/gate/actor/range_actor_a/stepHitType                      post<br>/gate/actor/range_actor_a/save  output/Xcat256_256_101_0.1_no_filter.mhd<br><br>/gate/actor/addActor ProductionAndStoppingActor range_actor<br>/gate/actor/range_actor/attachTo Xcat<br>/gate/actor/range_actor/addFilter                        particleFilter<br>/gate/actor/range_actor/particleFilter/addParticle       e+ # beta<br>/gate/actor/range_actor/addFilter                        energyFilter<br>/gate/actor/range_actor/energyFilter/setEmin             0.0 MeV<br>/gate/actor/range_actor/energyFilter/setEmax             0.0 MeV<br>/gate/actor/range_actor/setResolution                    256 256 101<br>/gate/actor/range_actor/stepHitType                      post<br>/gate/actor/range_actor/save  output/Xcat256_256_101_0.1_filter_positrons.mhd<br><br>/gate/actor/addActor ProductionAndStoppingActor range_actor_b<br>/gate/actor/range_actor_b/attachTo Xcat<br>/gate/actor/range_actor_b/addFilter                        particleFilter<br>/gate/actor/range_actor_b/particleFilter/addParticle       e+ # beta<br>/gate/actor/range_actor_b/setResolution                    256 256 101<br>/gate/actor/range_actor_b/stepHitType                      post<br>/gate/actor/range_actor_b/save  output/Xcat256_256_101_0.1_positrons.mhd<br><br>/gate/actor/addActor ProductionAndStoppingActor range_actor_c<br>/gate/actor/range_actor_c/attachTo Xcat<br>/gate/actor/range_actor_c/addFilter                        particleFilter<br>/gate/actor/range_actor_c/particleFilter/addParticle       gamma # beta<br>/gate/actor/range_actor_c/addFilter                        energyFilter<br>/gate/actor/range_actor_c/energyFilter/setEmin             0.511 MeV<br>/gate/actor/range_actor_c/energyFilter/setEmax             0.511 MeV<br>/gate/actor/range_actor_c/setResolution                    256 256 101<br>/gate/actor/range_actor_c/stepHitType                      post<br>/gate/actor/range_actor_c/save  output/Xcat256_256_101_0.1_filter_gammas.mhd<br><br>/gate/actor/addActor ProductionAndStoppingActor range_actor_d<br>/gate/actor/range_actor_d/attachTo Xcat<br>/gate/actor/range_actor_d/addFilter                        particleFilter<br>/gate/actor/range_actor_d/particleFilter/addParticle       gamma # beta<br>/gate/actor/range_actor_d/setResolution                    256 256 101<br>/gate/actor/range_actor_d/stepHitType                      post<br>/gate/actor/range_actor_d/save  output/Xcat256_256_101_0.1_gammas.mhd</div><div><br></div><div><br></div><div>/gate/physics/addPhysicsList emstandard <br><br>/gate/physics/addProcess PhotoElectric<br>/gate/physics/processes/PhotoElectric/setModel StandardModel<br><br>/gate/physics/addProcess Compton<br>/gate/physics/processes/Compton/setModel StandardModel<br><br>/gate/physics/addProcess RayleighScattering <br>/gate/physics/processes/RayleighScattering/setModel LivermoreModel<br><br>/gate/physics/addProcess ElectronIonisation<br>/gate/physics/processes/ElectronIonisation/setModel StandardModel e-<br>/gate/physics/processes/ElectronIonisation/setModel StandardModel e+<br><br>/gate/physics/addProcess Bremsstrahlung<br>/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-<br>/gate/physics/processes/Bremsstrahlung/setModel StandardModel e+<br><br>/gate/physics/addProcess PositronAnnihilation e+<br><br>#/gate/physics/addProcess MultipleScattering e-<br>#/gate/physics/addProcess MultipleScattering e+<br>/gate/physics/addProcess RadioactiveDecay</div></div></div>
</div>