Two Higgs doublet model predictions for (Formula presented) in NLO QCD

The decay (Formula presented) is studied at next to leading order in QCD in a class of models containing at least two Higgs doublets and with only one charged Higgs boson nondecoupled at low energy. The two-loop matching condition is calculated and it is found to agree with existing results. The complete dependence of the Wilson coefficients on the matching scale is given. The size of the next to leading order corrections is extensively discussed. Results for branching ratios, possible CP asymmetries, and lower bounds on the charged Higgs boson mass are presented when the convergence of the perturbative series appears fast enough to yield reliable predictions. Regions in the parameter space of these models where the next to leading order calculation is still not a good approximation of the final result for these observables are singled out.