Thank you very much for your helpful explanation!
Based on your suggestion, I added the RADDECAY card and performed the statistics for positrons and annihilation photons. I used USRBIN with the POSITRON PART and ANNIHRST PART scorings for this purpose. At the same time, I also recorded the positions of positrons and annihilation photons using mgdraw.f, with the conditions:
• KPART(KP) .EQ. 7 .AND. (ICODE .EQ. 214 .OR. ICODE .EQ. 215) for positrons
• KPART(KP) .EQ. 4 for annihilation photons
I would like to understand the difference between these two recording methods, and which one is more accurate or recommended in this context.
Additionally, I referred to https://fluka-forum.web.cern.ch/t/recording-position-of-particles-when-ke-goes-below-emfcut/1181 post and added the EMFCUT card in the setup with the following configuration:
EMFCUT -1E-06 1E-6 1 BLCKHOLE PMMA PROD-CUT
I’m not sure whether this threshold is low enough to properly simulate the behavior of low-energy annihilation photons and positrons.
I have also uploaded my input file for your reference.
Thank you again for your kind support!
Best regards
carbonposi.inp (2.1 KB)
mgdraw.f (15.8 KB)