Pseudorapidity dependence of the bulk properties of hadronic medium in pp collisions at 7 TeV

The measured charged particle \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_T$$\end{document}pT spectra in proton-proton collisions obtained by the CMS experiment at CERN is compared with the simulation results of EPOS–LHC and Pythia8.24 models at 7 TeV center-of-mass energy. The Pythia8.24 model describes the experimental data very well, particularly in the high \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_T$$\end{document}pT region. The model also predicts the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_T$$\end{document}pT spectra for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|$$\end{document}| \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document}η \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|$$\end{document}| < 2.4 at 0 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\le $$\end{document}≤ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_T$$\end{document}pT \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\le $$\end{document}≤ 6 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {GeV/}c$$\end{document}GeV/c. The EPOS–LHC model underpredicts the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_T$$\end{document}pT spectra from 0.1 to 2 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {GeV/}c$$\end{document}GeV/c in all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document}η bins for about 20% and the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_T$$\end{document}pT spectrum from 0.1 to 4.2 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {GeV/}c$$\end{document}GeV/c for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|$$\end{document}| \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document}η \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|$$\end{document}| < 2.4 by about 15% while reasonably predicts well for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_T$$\end{document}pT > 4.2 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {GeV/}c$$\end{document}GeV/c within the experimental errors. Furthermore, to get information about collective properties of the hadronic matter, modified Hagedorn function with embedded transverse flow velocity and thermodynamically consistent Tsallis distribution functions are used to fit the experimental data and simulated results. The values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2/ndf$$\end{document}χ2/ndf show that the functions fit the data and simulation results well. The parameter extracted by the functions: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _T$$\end{document}βT, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_0$$\end{document}T0, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{eff}$$\end{document}Teff decreases with increasing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document}η. The decrease in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _T$$\end{document}βT with increasing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document}η is due to the large energy deposition in lower rapidity bins producing rapid expansion due to large pressure gradient resulting quick expansion of the fireball. Similarly, large energy transfer in the lower pseudo-rapidity bin results in higher degree of excitation of the system which results larger values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_0$$\end{document}T0 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{eff}$$\end{document}Teff. The values of the fit constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_0$$\end{document}N0 increase with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document}η where the values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_0$$\end{document}N0 extracted from Pythia8.24 are closer to the data than the EPOS–LHC model. The Pythia8.24 model has better prediction than the EPOS–LHC model which might be connected to its flow-like features and color re-connections resulting from different Parton interactions in the initial and final state.

High energy particle collision is a very complex topic that makes it hard to measure and understand the geometry, quantities, and global properties of the collisions. Therefore, in addition to the measurements of some quantities by experiments, one must rely on theoretical models to explain some other characteristics based on experimental results. Experiments using high energy pp collisions such as the Relativistic Heavy Ion Collider (RHIC) 1,2 and the Large Hadron Collider (LHC) have been very useful to study the characteristics of these collisions such as the flow effects, the pseudorapidity, and the transverse momentum distributions. On the other hand, different models have been presented in the literature which tried to describe the collision process at different stages 3,4 .
A measurement of pseudorapidity density provides constraints to the modeling of the characteristics of pp collisions such as the direction of particle emission. The pseudorapidity ( η ) is given by η = − ln(tan(θ/2)) , where θ is the polar angle that the charged particles make with the anticlockwise beam direction. The presented analysis are of dN ch / d η and dN ch /dp T in a pseudorapidity range of η < 2.4 in steps of 0.2. In the current analysis, the charged hadrons ( N ch ) includes decay products of particles with a lifetime of < 1 cm, while the products of secondary particles are excluded and a correction factor is applied for prompt leptons.
The p T and η distributions of charged particles, pions, kaons, protons, and anti-protons produced in nucleusnucleus, hadron-nucleus, and pp collisions are very important observables because they provide very crucial information about the anisotropy and dynamics of the final state particles produced in collisions at high energies. In addition, p T spectra are different for different particles in the range of 0 to 100 GeV/c , which makes it difficult

Method and models
Pseudorapidity (η ) distributions in pp collisions at 7 TeV are obtained from the inner tracking system of the Compact Muon Solenoid experiment (CMS) at the Large Hadron Collider (LHC) 26 on March 30, 2010. The spectra obtained by the experiment were normalized to all non single diffractive events with corrections for trigger and selection efficiency, acceptance, and branching ratios. Simulations of two Monte Carlo models, EPOS 24 and Pythia 25 at the same energy are performed and are then compared with the measurements from the experimental data. These models are based on the simple Parton model using theory of Reggie-Gribov 27 , which is a QCD based effective field theory that accounts for multiple interactions in parallel with Pomerons considered to represent the partons interactions 28 . A brief description about each of these models is provided below.
The generator EPOS 24 is a hadronic interaction Monte Carlo simulation package that is well known to simulate the hadronic interactions in high energy cosmic ray simulations. In addition, it is used to describe the minimum bias interactions and centrality dependence of heavy-ion collisions. EPOS is an acronym of Energy conserving quantum mechanical multiple scattering approach based on Partons Off-shell remnants and Splitting of parton ladders. Soft, semi-soft and hard Pomeron exchanges are used to describe the interaction between particles in EPOS, where the particle production originates from two kind of sources, cut Pomeron and remnant decays 29 . In its last improvement, EPOS was updated to a new version, called EPOS-LHC 24 , to cover the energies of the LHC. The new version can reproduce all minimum bias results of all particles with transverse momentum from 0 to a few GeV/c . In addition, in case of very dense system in a small volume in pp collisions, a different parametrization of flow has been introduced in EPOS, compared to the large volume produced in heavy-ion collisions. In this paper, we have used the EPOS-LHC version of the EPOS model but for simplicity only EPOS will be used throughout the manuscript.
Pythia is one of the most widely used Monte Carlo event generator of particle collisions in high energy physics with emphasis on pp interactions 25 . The hadronization process of transforming the final outgoing coloured Partons into colourless particles is based on the Lund string fragmentation model 30 . The main event in a pp collision can be represented by a large number of processes, such as elastic and diffractive processes, electroweak processes, QCD hard and soft processes and top quark production. The model also implements the initial and final state radiations, multi Partonic interactions, beam remnants which are formed after the extraction of multi Partonic interactions 31 . We used Phythia8.24 version of the Pythia model but for simplicity we will be using Pythia throughout the manuscript. It is also pertinent to mention that one million events are simulated in case of both the models.
In a standard analysis, it is always important to compare the results of the simulation generated with the Monte Carlo models to experimental data. In addition to the comparison of models' prediction with the experimental data, we use the theoretical Hagedorn function 32 to fit the experimental data as well as the models' simulations to extract the freeze-out parameters. The Hagedorn function, a QCD -inspired inverse power law, which can produce the transverse mass m T distribution of hadrons in pp and AA collisions. This Hagedorn function describes the bulk spectra in the low transverse mass region as well as the particles produced in QCD hard scatterings. The function is given as: Here, m t = p 2 T + m 2 0 is the transverse mass of hadrons, m 0 the rest mass while C is the normalization constant. The n and p 0 in the equation are two free parameters of the function. Moreover, the Tsallis function [15][16][17][18] can excellently describe the p T and m T invariant distribution measured in pp collisions at high energies. There are several version of the Tsallis function that can give good fit results to the p T spectra, however, the following (1) d 2 N 2πN ev p T dp T dy www.nature.com/scientificreports/ expression is a simple version of the Tsallis function 22 that describes the invariant spectra of particles in terms of the effective temperature T eff and non-extensivity parameter q which accounts for the deviation of the p T spectra from the usual Boltzmann-Gibbs exponential distribution function. This version of the Tsallis function is a consistent version of thermodynamics for the particle number, energy, density and pressure at mid-rapidity ( y ≈ 0 ), and is given by the following expression: where C q is constant of the fit function, T eff is the effective temperature, q is the non-extensive parameter and can also be considered as a measure of the non-thermalization 23 . When the parameter q is close to one, the thermalization degree of a system is larger and the Tsallis distribution approaches the normal Boltzmann-Gibbs exponential distribution function. On the other hand, the effective temperature parameter T eff represents the contribution from the thermal motion of particles as well as the collective flow of expanding matter 22 .
which is called the thermodynamically consistent Tsallis function throughout the paper or Tsallis function for shortness.
Assuming that n = (q − 1) −1 and p 0 = nT , we notice that Eqs. (1) and (2) are mathematically equivalent. One can see that the parameters n and q are inversely proportional to each other, as the value of the parameter n increases, the parameter q decreases.
Finally to include the transverse flow in Eq. (4), m T is equated with These modifications have been discussed in 18,33,34 and have successfully calculated the average transverse flow velocity < β T > and T 0 , and the equation is now known as the Hagedorn equation with embedded transverse flow velocity 35 : where C is the normalization constant and is to be normalized to 1. < γ t > = 1/ 1− < β 2 t > and < β t > is the average transverse flow velocity. This simple form of the equation having a few parameters, is very powerful tool to compare different collisions with a small number of parameters. The Hagedorn function has reproduced the spectra described in 18,32,33 with physical parameters. We can re-arrange these terms in our current analysis as follows, The p T spectra of charged hadrons from experiment as well as model simulations are fitted with Hagedorn equations (Eq. 7) and thermodynamically consistent Tsallis function (Eq. 8) to extract the fit parameters as explained above. We have used the method of least square for the extraction of the related parameters. The Tsallis function is very effective in providing an excellent description of the transverse momentum spectra of particles, for nucleus-nucleus (AA) as well as hadron-hadron (hh) collisions at high energies. The function has three free parameters including normalization/fitting constant indicating volume information, the non-extensivity parameter used to see the deviation of the distribution from the exponential Boltzman and Gibbs distribution, and the effective temperature which includes the flow information along with the thermal motion of the particles. Furthermore, to describe p T distribution of particles in AA and hh collisions, different models with flow definitions are included in the Tsallis distribution function. Hagedorn function is one of them which among other parameters gives direct access to the transverse expansion (flow velocity) and kinetic freeze-out temperature with an excellent description of the p T spectra at high energies. In addition, the modified Hagedorn function is very close to the ideal gas model. Both the Hagedorn and Tsallis models have the advantage that they cover a wide range of p T due to the entropy parameter q(n) where q = 1/(n − 1). . .

Results and discussion
The transverse momentum ( p T ) spectra of charged hadrons simulated with the EPOS and Pythia models are compared with the pp collisions at √ s = 7 TeV measured by the CMS experiment 26 . The initial conditions used for simulations are similar to that of experimental conditions. The p T spectra from 0.1 to 2 GeV/c in different η bins from 0 to 2.4 in steps of 0.2, and the p T spectrum from 0.1 to 6 GeV for 0 ≤ η < 2.4 as a single η bin are studied using the aforementioned models. Figure 1 shows the p T distributions of charged hadrons in different η bins starting from 0.1 (corresponding to η = 0.0−0.2 ) from top left to η = 2.3 ( 2.2−2.4 ) to the bottom right. The experimental data is shown by solid black markers while lines of different colors represent the two models calculations. A red solid line is used to show the results of the EPOS model, while the blue solid line is for Pythia results. The quadrature sum of the systematic and statistical errors are considered in the experimental results, while model calculations include the statistical errors only. The horizontal error bars show the bin size shown in the experimental results. The two event generators have the same horizontal bin size in p T as the experimental data. It has been observed that for p T above 0.8 GeV/c , Pythia measurements fully describe the p T spectra at all the pseudorapidity regions. Below 0.3 GeV/c , Pythia model under-predicts the experimental data, while for 0.3 ≤ p T ≤ 0.8 GeV/c , the model slightly overshoot the experimental data which is about 10%. There is no significant η dependence observed in the model's predictions and have similar results for all pseudorapidity regions. The EPOS model under-predicts the experimental data over the entire p T range up to about 20% and is also observed to be independent of all η regions understudy. The data to Monte Carlo ratio is shown in the lower panel of each plot which supports the above statement. Figure 2 shows the predictions of the p T spectra of the two models in comparison to the experimental data for η < 2.4 GeV/c . Again the Pythia model reproduced the experimental data over most of the p T region within the experimental errors. For p T < 1 GeV/c , the model slightly overshoot the data in a narrow region of the p T for about 10%. The EPOS model again underestimates the experimental data over the entire p T region with decreasing discrepancy with increasing p T . The ratio of the Monte Carlo predictions to the data is also shown in the lower panel of the plot which supports our statements. Figure 3 shows the η of charged hadrons integrated over the p T for − 2.5 ≤ η ≤ 2.5 in comparison with the models' predictions. Pythia model reproduced the experimental results of dN ch /dη very well for the whole η region, − 2.5 ≤ η ≤ 2.5. The EPOS model does not reproduce the data and hence underpredicts the charged particles η integrated over p T over the entire region of − 2.5 ≤ η ≤ 2.5 from 10% to about 15%. The lower panel of the graph shows the ratio of the Monte Carlo predictions to the data. The description mentioned can easily be inferred from the ratio plot as well.
To extract the bulk properties of the hadronic matter created as a result of pp collisions at √ s = 7 TeV, the CMS experimental data 36 and models' simulations are fitted with two statistical functions. The fit results of the data and two models' predictions are presented in Fig. 4 which show that the fit functions fit the data very well. The experimental data and model simulations are fitted by the modified Hagedorn function with embedded transverse flow velocity (Eq. 7) and thermodynamically consistent Tsallis distribution function (Eq. 8). Figure 4a Table 1. The values of χ 2 /ndf also show that the two functions fit the data and models' predictions very well given in the last column of the table. The values of transverse flow velocity ( β T ), n, and Kinetic freeze-out temperature ( T 0 ) are extracted by fitting the spectra of the experimental data and models' simulations tabulated as a third, fourth, and fifth column in the Table 1. These parameters are directly connected to the scattering centers involved in the interaction process.
The values of different parameters are extracted by fitting the experimental data and models' simulations in each η bin for better comparison of predictions with the measurements. The values of N 0 , β T , and T 0 decrease with η in data and both the models. The values of β T and T 0 extracted by the fit functions from the EPOS model are closer to the experimental data than the Pythia model while opposite in the case for N 0 where the latter has similar values in Pythia with the data than the EPOS model.
In case of experimental data, the function described in Eq. (7) gives the value of β T = (0.3836 ± 0.0002 ) c for η = 0.1 , which corresponds to the η range from 0.0 to 0.2. The extracted value of β T = (0.3505 ± 0.0002 ) c for the η = 2.3 . In the case of Pythia model, the value of β T varies from ( 0.6178 ± 0.0002 ) c for η = 0.1 to ( 0.593 ± 0.0002 ) c for η = 2.3 , whereas it varies from ( 0.482 ± 0.0002 ) c to ( 0.347 ± 0.0002 ) c in case of EPOS model for the first and last regions of η . It has been observed that the variation in the value of β T is monotonic in all the cases. The kinetic freeze-out temperature T 0 , extracted by fitting the experimental data with Hagedorn function also shows a decreasing trend with increasing η starting from ( 76.37 ± 0.01 ) MeV at η = 0.1 to ( 71.86 ± 0.01 ) MeV at η = 2.3. The reason behind this decreasing trend of T 0 is the decrease of the energy transfer in the system due to the large penetration between participants particles as the system goes from mid-rapidity to the forward-rapidity region and this result is in agreement with our recent work 37 . A similar decreasing trend of the T 0 is observed by applying the fit function on the two models' predictions. In both of the cases, the value of  The relation of q with n can be established by comparing the two tables resulting from two different functions. A slight variation in q yields an appreciable inverse variation in the value of n. It has resulted from this study that the n parameter decreases while the q increases with η.
An increasing trend is observed in the values of N 0 with increasing η extracted by the fit using Tsallis function. The value varies from 93.0 ± 0.2 to 97.1 ± 0.2 in case of experimental data while for Pythia it yields 88.8 ± 0.2 for η = 0.1 while 98.4 ± 0.2 for η = 2.3 . The function yield a value of 77.2 ± 0.2 for the lower η bin while 86.3 ± 0.2 for the higher η bin. The values of N 0 in case of the Pythia model is closer to the experimental data than the EPOS model. Since N 0 is proportional to the multiplicity of particles and Pythia reproduced similar multiplicity hence predicted better results than the EPOS model. This might be connected to the effects that are incorporated in the Pythia model such as flow-like effect and color reconnection effect which is produced from simultaneous hard sub-collisions forming color strings between initial and final state partons from separate hard scatterings due to which the model predicts the data well.
Before going to conclusion section we would like to clarify that the T 0 and T e eff are not different. The former does not include the flow effect, however the later includes the flow effect. We observed that both of them decrease with increasing pseudo-rapidity due to the large energy transfer in lower pseudo-rapidity intervals. The physics behind this is that larger energy transfer in the lower pseudo-rapidity bin results in higher degree of excitation of the system which results in larger T 0 and T eff 10,38-41 and this indicates that the system in lower pseudo-rapidity www.nature.com/scientificreports/ intervals comes to equilibrium quickly. In addition we also observed that q (n) is increasing (decreasing) with the increase of rapidity. As discussed above, q and n are reciprocal. A system with smaller q (more closer to unity) and larger n indicates to be more closer to equilibrium state [42][43][44][45] . In the present work the increase (decrease) of q (n) with increasing η claims that the system goes far from the equilibrium state as the psuedo-rapidity intervals is larger. The β T is also observed to decrease with increasing η due to large energy deposition in lower rapidity bins. Actually, the large deposition of energy in a system indicates its rapid expansion. In the present work, the larger β T in lower η bins shows that there is large pressure gradient due to the fact that large amount of energy is transferred to the system which results in a quick expansion of the fireball.

Summary and conclusion
Simulations of EPOS and Pythia models are performed and then compared with the measurements of the CMS experimental data in pp collisions at √ s = 7 TeV. The transverse-momentum ( p T ) spectra are presented from 0.1 to 2 GeV/c in twelve pseudorapidity ( η ) intervals from 0 to 2.4 in steps of 0.2 and the p T spectrum from 0.1 to 6 GeV/c for a wider η interval, 0 ≤ η ≤ 2.4 as a single bin. The Pythia model predicts the p T spectra at all the pseudorapidity bins very well particularly at the higher value of p T , while over predicts the distribution slightly at the lower p T region. The model also predicts the p T spectra for a wider range of η , η < 2.4 and over a wide range of p T , 0 ≤ p T ≤ 6 GeV/c . However, within the experimental error a slight bump is predicted in a narrow region at lower p T interval. The EPOS model under-predicts the p T spectra over the entire p T range and in all η intervals. The model also underestimates the experimental results of p T spectra for all η intervals for about 15% at η < 2.4 while a reasonable agreement is shown for p T > 4.2 GeV/c within the experimental errors. The Pythia model reproduce the charged particle pseudorapidity integrated over p T for the whole region of η , − 2.5 ≤ η ≤ 2.5 very well whereas the EPOS model under-predicts from 10% to about 15%.
Furthermore, the models' predictions and the experimental data are fitted by two statistical functions to get information about some collective properties of the hadronic matter. The measured experimental data and models simulations are fitted by modified Hagedorn function with embedded transverse flow velocity and thermodynamically consistent Tsallis distribution function. The two function fitted the data and models prediction very well. The values of transverse flow velocity ( β T ), and Kinetic freeze-out temperature ( T 0 ) extracted from the Hagedorn distribution function while effective temperature ( T eff ) is obtained from the Tsallis function. The value of β T , T 0 , and T eff decrease with increasing η because large energy transfer in the lower pseudo-rapidity bin results in higher degree of excitation of the system which results larger values of T 0 and/or T eff . Furthermore, www.nature.com/scientificreports/ the decrease in β T with increasing η is due to the large energy deposition in lower rapidity bins producing rapid expansion. In the present work, the larger β T in lower η bins shows that there is large pressure gradient due to the fact that large amount of energy is transferred to the system resulting in a quick expansion of the fireball. It is concluded that the Pythia model describes the experimental data at most of the p T regions for all the η bins, while the EPOS underpredicts mostly. The Pythia model also predicted the p T distribution over a wider p T range well which might be connected to the color re-connection and flow-like feature of the Pythia model, whereas the EPOS reproduced the distribution at higher values of the p T only and largely underpredicts the distribution. Although, the models under-study reproduced the p T distribution of the charged particles in some region of p T in different η regions presented from 0 to 2.4 in steps of 0.2, but none of them completely describe all the distribution over the entire p T range.