强子的横向动量依赖的部分子分布函数由于包含了部分子横向动量的分布情况，对强子内部结构和非微扰现象进行了更完备的描述。在强子对撞过程中，重夸克偶素在小横向动量区域的产生给从实验数据中提取横向动量依赖的胶子分布函数提供了理想的平台。为达到这个目的我们需要对这些过程建立横向动量依赖的因子化定理，以保证所提取分布函数的普适性。本文研究了强子对撞中，S波重夸克偶素ηc，b和P波重夸克偶素在小横向动量区域产生的因子化。对于S波重夸克偶素的产生，通过计算单圈的横向动量依赖的胶子分布函数和微分散射截面，我们确认了因子化成立，并将微扰系数确定到了单圈水平。通过对比树图阶和单圈阶的因子化，我们发现需要引入一个软因子来抵消所有的软发散。这些结果可以用来在未来实验数据中提取横动量依赖的胶子分布函数。用类似的方法，我们也研究了P波重夸克偶素产生的因子化。我们发现因子化对于JPC=0++，2++重夸克偶素在单圈水平仍然成立，但在高圈阶会破坏。对于其它P波重夸克偶素，因子化在单圈阶就不再成立。因此P波重夸克偶素的产生不适合用来从实验数据中提取横向动量依赖的胶子分布函数。