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

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.