Recording the positions of stopping protons

Dear Experts and Users
My simulation consists of monitoring the positions and energies of the stopping protons through the target(this is in continuation of my earlier queries!). I want to record the positions and deposited energies at the moment when protons stop in the target(here exists different materials in the target, Ta and H2O).
And I set the transport threshold for protons in the PARA-THR at 1e-14 GeV. I have modified a mgdraw routine. I am using ENDRAW entry. Is this the correct way to go or is there an another way? I am attaching the input file and mgdraw file.
target.inp (7.5 KB)
mgdraw_end.f (3.5 KB)

Thanks in advance.

Best regards,

Qi Ding

1 Like