We study the role of the non-perturbative input to the transverse momentum dependent (TMD) gluon density in hard processes at the LHC. We derive the input TMD gluon distribution at low scale mu0^2 ~ 1 GeV^2 from the fit of the inclusive hadron spectra measured at low transverse momenta in pp collisions at the LHC and demonstrate that the best description of these spectra for larger hadron transverse momenta can be achieved by matching the derived TMD gluon distribution with the exact solution of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation obtained at low x and small gluon transverse momenta outside the saturation region. Then, we extend the input TMD gluon density to higher mu^2 numerically using the Catani-Ciafoloni-Fiorani-Marchesini (CCFM) gluon evolution equation. A special attention is put to the phenomenological applications of obtained TMD gluon density to some LHC processes, which are sensitive to the gluon content of a proton.