Saturation models of HERA DIS data and inclusive hadron distributions in p+p collisions at the LHC