  We study the impact of initial-state multiple parton radiation on transverse
momentum $(q_T)$ distribution of Higgs boson produced via bottom quark fusion
at hadron colliders. The shape of the resulting $q_T$ distribution is affected
by the bottom-quark mass corrections and by the strong kinematical behavior of
the bottom-quark parton density. We account for both features in the full range
of $q_T$. To do this, we formulate the resummation calculation in a
general-mass factorization (S-ACOT) scheme and introduce a correction in the
resummed-term to account for the effect from large-$q_T$ kinematics of Higgs
boson. The results of this resummation are compared to fixed-order and PYTHIA
predictions.
