Gluon saturation in pA collisions at energies available at the CERN Large Hadron Collider: Predictions for hadron multiplicities

The upcoming p+Pb run at the LHC will probe the nuclear gluon distribution at very small Bjorken x (from $x \sim 10^{-4}$ at mid-rapidity down to $x \sim 10^{-6}$ in the proton fragmentation region) and will allow to test approaches based on parton saturation. Here, we present the predictions of the KLN model for hadron multiplicities and multiplicity distributions in p+Pb collisions at a center-of-mass energy of 4.4 TeV. We also compare the model to the existing pp, dA and AA data from RHIC and LHC.