We calculate inclusive hadron productions in pA collisions in the small-x saturation formalism at one-loop order. The differential cross section is written into a factorization form in the coordinate space at the next-to-leading order, while the naiv
e form of the convolution in the transverse momentum space does not hold. The rapidity divergence with small-x dipole gluon distribution of the nucleus is factorized into the energy evolution of the dipole gluon distribution function, which is known as the Balitsky-Kovchegov equation. Furthermore, the collinear divergences associated with the incoming parton distribution of the nucleon and the outgoing fragmentation function of the final state hadron are factorized into the splittings of the associated parton distribution and fragmentation functions, which allows us to reproduce the well-known DGLAP equation. The hard coefficient function, which is finite and free of divergence of any kind, is evaluated at one-loop order.
We demonstrate the QCD factorization for inclusive hadron production in $pA$ collisions in the saturation formalism at one-loop order, with explicit calculation of both real and virtual gluon radiation diagrams. The collinear divergences associated w
ith the incoming parton distribution of the nucleon and the outgoing fragmentation function of the final state hadron, as well as the rapidity divergence with small-$x$ dipole gluon distribution of the nucleus are factorized into the splittings of the associated parton distribution and fragmentation functions and the energy evolution of the dipole gluon distribution function. The hard coefficient function is evaluated at one-loop order, and contains no divergence.
