The presence of geometric scaling within the p T spectra of produced hadrons at high energy pp collisions using small-x k T-factorization is investigated. It is proposed a phenomenological parameterization for the unintegrated gluon distribution in the scaling range that reproduces the features of the differential cross section both in the saturated and dilute perturbative QCD regimes. As the saturation scale acts as an effective regulator of the infrared region, the extension of the model to quantities usually associated to soft physics is studied. The approach is applied to compute the average p T and the rapidity distribution of produced gluons at high energies.