[1] E. Gotsman, E. Levin, U. Maor and J.S. Miller, Eur. Phys. J. C57 (2008) 689 [arXiv:0805.2799 [hep- ph]]; [arXiv:0903.0247 [hep-ph]]. [2] Alice collaboration, Eur. Phys. J. C. 65 (2010) 111 [arXiv:0911.5430 [hep-ex]]. [3] CMS Collaboration, JHEP 02 (2010) 041 [arXiv:1002.0621 [hep-ex]]. [4] M. L. Good and W. D. Walker, Phys. Rev. 120 (1960) 1857. Acknowledgement: This research was supported in part by BSF grant #20004019. To summarize, we have presented a theoretical formulation for single inclusive hadron-hadron interactions based on our GLMM model. We have reproduced the p-p data[2, 3, 27] on single inclusive density as a function of the pseudo rapidity. Our results provide additional support for our proposed Pomeron approach. We have also presented predictions for the LHC en- ergy range. These predictions may soon be tested during the preliminary low luminosity LHC runs. 5 Despite this, we consider this fit to be acceptable, as the data points ”oscillate” about a uniform line with error bars which are much smaller than their deviation from a smooth average. The results of this fit are a$_{IPIP}$ = 75.7, a$_{IPIR}$ = 0.12 and Q = 3.8 GeV. In our procedure, the line for 546 GeV in Fig. 2 is calculated with the model parameters and is visually compatible with the experimental data points. Note that both the axes of Fig. 2 are linear, and that our calculation coincides with the LHC experimental results [2] and [3]. We have also made predictions for the higher energies at which the LHC is expected to run, see Fig. 2. The contributions of the secondary Regge trajectories are minimal. The experimental values for σ$_{in}$ = σ$_{tot}$ − σ$_{el}$ − σ$_{diff}$ were taken from Refs.[2, 3, 27]. For our predictions we have used the values of σ$_{in}$ calculated in our GLMM model. Our output over-estimates the few data points with η > 4 data at 546 and 900 GeV by up to 20%. This is to be expected, as we have not taken into account the parton correlations due to energy conservation, which are important in the fragmentation region, but difficult to include in the framework of Pomeron calculus. As noted, this calculation entails three additional parameters. The determination of these parameters from existing data [27] is not trivial. Comparing the numbers corresponding to the data shown in Fig. 2, it is evident that a conventional overall χ 2 analysis is impractical, owing to the quoted error bars of the 546 GeV data points, which are considerably smaller than the error bars quoted for the other energies. The full lines in Fig. 2 are the results derived from a χ 2 fit to the 200-1800 GeV data, excluding the 546 GeV points. This fit yields a seemingly poor χ $^{2}$/d.o.f = 3 . 2.