How to generate process higgs+gamma for u/d/s quark

Dear experts,

I am currently trying to generate processes like u u~ > h a / d d~ > h a /s s~ > h a, but I found that I can not generate them using the loop_sm model. However, when I add the default restrict_c_mass.dat, I can generate c c~ > h a process correctly.

Actually, I found no vertices and couplings corresponding to u/d/s quark and higgs in this model, so it did not work when I only added the Yukawa values and masses for uds quark. and the calculations involve the goldstone couplings so it seems not easy to fix this model directly.

We did find some EFT models (like HEl_UFO [1] ) that can achieve these processes( like u u~ > h a ) at LO. But we would like to get SM results and be better at NLO. Could you give us some advice about this?

I really appreciate any help you can provide.

Best Regards,
Zhe Guan

[1] https://feynrules.irmp.ucl.ac.be/wiki/HEL#no1