How to separate pKa by reaction channel (n,p), (n,α)… in mgdraw.f?

Dear FLUKA experts,

At the moment I retrieve pKa simply by setting icode = 300.
https://fluka-forum.web.cern.ch/t/primary-knock-on-atom-pka-using-mgdraw/7752?u=131415248
I would like to go one step further and distinguish the contributions coming from different nuclear reaction channels, e.g. (n,p), (n,α), etc., before merging them for subsequent analysis.

  1. Is it possible to implement such a channel-wise selection inside mgdraw.f?

  2. If yes, which variables / common blocks should I use to identify the channel reliably?

Any hints or minimal code snippets are highly appreciated.

mgdraw_empty (1).f (10.3 KB)

Thank you in advance.

Cheers,
Yingming

You should gather all reaction products that are split between genstk.inc and fheavy.inc, as already explained. Also, you can cross-check that nothing was forgotten, out of the neutron induced reaction, by looking at the target nucleus’ mass and atomic numbers (IBTAR and ICHTAR in resnuc.inc).

Dear Francesco,

Could you kindly clarify whether, for low-energy neutron reactions (icode = 3xx), it would be possible to further distinguish the specific reaction channels—such as elastic scattering, inelastic scattering, and neutron capture?

Thank you in advance.

Best regards

Yes. Note that for low-energy neutron interactions the only available ICODE value is 300.
Still, if you look at the reaction products as earlier mentioned, you can naturally figure out the specific reaction channel: if you find just one proton and possible gammas (in genstk.inc) and a residue (in fheavy.inc), it was (n,p); if you find just one neutron and some de-excitation gamma (in genstk.inc) and a residue matching the target nucleus (in fheavy.inc), it was (n,n’); if you find nothing in genstk.inc and a residue matching the target nucleus in fheavy.inc, it was (n,el); if you find just gammas (in genstk.inc) and a residue (in fheavy.inc), it was capture, namely (n,gamma); and so on for any possible channel.
Therefore, gathering all the interaction products, you have the exact picture of the reaction channel.