<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body>
<div dir="ltr" style="">
<div></div>
<div style="">
<div>
<div dir="ltr" style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
Hi Dimitris, </div>
<div dir="ltr" style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
<br>
</div>
<div dir="ltr" style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
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.</div>
</div>
<div dir="ltr" style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
<br>
</div>
<div dir="ltr" style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
You can see how to do it here: </div>
<div dir="ltr" style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
<a rel="noreferrer noopener" href="https://opengate.readthedocs.io/en/latest/cut_and_variance_reduction_technics.html">https://opengate.readthedocs.io/en/latest/cut_and_variance_reduction_technics.html</a><br>
</div>
<div dir="ltr" style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
<br>
</div>
<div dir="ltr" style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
Thanks,</div>
<div dir="ltr" style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
Mahdjoub</div>
<div id="ms-outlook-mobile-signature">
<div><br>
</div>
Get <a href="https://aka.ms/o0ukef" style="">Outlook for iOS</a></div>
</div>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Gate-users <gate-users-bounces@lists.opengatecollaboration.org> on behalf of Dimitris Anestis <janest755@gmail.com><br>
<b>Sent:</b> Tuesday, April 12, 2022 4:42:25 AM<br>
<b>To:</b> gate-users@lists.opengatecollaboration.org <gate-users@lists.opengatecollaboration.org><br>
<b>Subject:</b> [Gate-users] How to store the stopping position of positrons of 511keV</font>
<div> </div>
</div>
<div>
<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>
</div>
</body>
</html>