Resummation and renormalization of kinematical effects in inclusive $P$-wave quarkonium production

We investigate the renormalization properties of the shape function formalism for inclusive production of $P$-wave heavy quarkonia, which arises from resumming a class of corrections coming from kinematical effects associated with the motion of the heavy quark and antiquark pair relative to the quarkonium. Such kinematical effects are encoded in the nonperturbative shape functions, which are normalized to the corresponding nonrelativistic QCD long-distance matrix elements. By using the known ultraviolet divergences in the matrix elements, we derive the large-momentum asymptotic behavior of the shape functions. This strongly constrains the form of the shape functions and significantly reduces the dependence on the nonperturbative model. Based on these results we show that the shape function formalism at loop level can be useful in taming the threshold logarithms at large transverse momentum, and at small transverse momentum the kinematical corrections reduce the sizes of $\chi_c$ and $\chi_b$ cross sections which may improve agreement with measurements.