We show that in order to obtain a successful description of the transverse momenta distribution for charged particles in ion-ion collisions, one must include a thermal emission term. The temperature of this emission T th turns out to be proportional to the saturation scale, T th = 1.8/2π Qs. The formalism for the calculation of the transverse momenta spectra in CGC/saturation approach is developed, in which two stages of the process are seen: creation of the colour glass condensate, and hadronization of the gluon jets. Our calculations are based on the observation that even for small values of pT , the main contribution in the integration over the dipole sizes stems from the kinematic region in vicinity of the saturation momentum, where theoretically, we know the scattering amplitude. Non-perturbative corrections need to be included in the model of hadronization. This model incorporates the decay of a gluon jet with effective mass m 2 eff = 2Qsµ soft where µ soft denotes the soft scale, with the fragmentation functions at all values of the transverse momenta.We use the KLN model which, provides a simple way to estimate the cross sections for the different centrality classes. Comparing the results of this paper with the transverse distribution in the protonproton scattering, we see two major differences. First, a larger contribution of the thermal radiation term is needed, in accord with higher parton densities of the produced colour glass condensate. Second, even changing the model for the hadronization, without a thermal radiation term, we fail to describe the pT spectrum. Consequently, we conjecture that the existence of the thermal radiation term is independent of the model of confinement. PACS numbers: 12.38.Cy, 12.38g,24.85.+p,25.30.Hm * We have to mention, that actually the model of hadronization with given fragmentation function describes the experimental data at large values of p T (p T ≥ 3 GeV [7] and p T ≥ 5 GeV [38]). In our model of the hadronization we assume that we can describe the data with the fragmentation function at any value of p T . We will demonstrate that such a description is possible.