Neutron and isotopes dump with mgdraw.f

Advanced Features and User Routines
5

Dear @hanhechong,

  1. The mother particle for JTRACK is incoming one, and you would be interested in the target as well from what I understand. For your process, in addition to the proton, you could maybe have nuclei that you seem to be interested in.
    From a separate question of yours, I see that you are using the stuprf.f routine, but this will store the information from ALL the subsequent interactions, which might be too much. Alternatively, you could use the mdstck.f subroutine in a similar way.
    To get the target nuclei, you would have to add INCLUDE 'resnuc.inc' to access the parameters that you would be interested in.

  2. The isotopes will not have corresponding particle codes, but you can nevertheless retrieve their A and Z numbers, as suggested in this post:

      IARES  = MOD ( MOD ( IAZTRK , 1000000), 1000  ) 
      IZRES  = MOD ( IAZTRK, 1000000) / 1000
      ISRES  = IAZTRK/ 1000000

Hope this helps.

Best,
Daniel

2 Likes