Charged Particles Transverse Momentum and Pseudorapidity Distribution in Hadronic Collisions at LHC Energies

We present an analysis of the pseudorapidity η and transverse momentum pT distributions of charged hadrons in pp collisions for the kinematic range of 0<pT<4 GeV/c and |η|<2.4 at 0.9, 2.36, and 7 TeV. Charged particles are produced in pp collision using several Monte Carlo event generators (Pythia Simple, Vincia, Dire showers, Sibyll2.3d, QGSJETII-04, EPOS-LHC) and compared with CMS data at LHC. It is observed that the Simple parton showers can explain the CMS data very well for pT>1 GeV/c at 0.9 and 2.36 TeV within the experimental errors, while Dire overshoots and Vicia undershoots the data by 50% each. At 7 TeV, the Dire module presents a good prediction, whereas the Simple and Vincia modules underestimate the data within 30% and 50%. Comparing the Simple module of the Pythia model and the predictions of the CRMC models with the experimental data shows that at 0.9 TeV, EPOS-LHC has better results than the others. At 2.36 GeV, the cosmic rays Monte Carlo (CRMC) models have better prediction than the Simple module of Pythia at low pT, while QGSJETII-04 predicts well at high pT. QGSJETII-04 and EPOS-LHC have closer results than the Pythia-Simple and Sibyll2.3d at 7 TeV. In the case of the pseudorapidity distributions, only the Pythia-Simple reproduced the experimental measurements at all energies. The Dire module overestimates, while Vincia underestimates the data in decreasing order of discrepancy (20%, 12%, 5%) with energy. All CRMC models underestimate the data over the entire η range at all energies by 20%. The angular ordering of partons and the parton fragmentation could be possible reasons for this deviation. Furthermore, we used the two-component standard distribution to fit the pT spectra to the experimental data and extracted the effective temperature (Teff) and the multiplicity parameter (N0). It is observed that Teff increases with the increase in the center of mass energy. The fit yielded 0.20368±0.01, 0.22348±0.011, and 0.24128±0.012 GeV for 0.9, 2.36, and 7 TeV, respectively. This shows that the system at higher energies freezes out earlier than lower ones because they quickly attain the equilibrium state.


Introduction
It is essential to measure particle yields and kinematic distributions to explore new energy regimes of hadronic collisions. Such studies provide insights into the soft and hard scattering contributions and a deep understanding of the hadron production, which in its turn provides a better platform for new investigations [1][2][3]. For the pp collisions at LHC, an intensive knowledge of inclusive particle production will make it easy to distinguish the background from the signal in soft hadronic interactions. Furthermore, the pp study provides a baseline for the heavy systems by giving information on nuclear medium effects [4,5]. In pp collisions, soft interactions produce most particles and can be modeled phonologically. The most widely used models and event generators can be tuned with the help of comparison to experimental data. These soft collisions include the inelastic single diffractive (ND), non-single diffractive (NSD), double diffractive (DD), and elastic scatterings [6].
The current study aims to contrast the prediction of event generators using the transverse momentum (p T ) and pseudorapidity (η) distributions of charged particles at several LHC energies. Such parameters are key observables to distinguish general properties of the collision that are more significant for higher collision energies [7,8]. Similar studies showing models' prediction compared to experimental data at different energies are given in the references [9][10][11][12][13][14][15][16][17][18][19][20]. The data at 0.9 and 2.36 TeV are taken from [4] and 7 TeV from [5].
The rest of the paper is organized as follows: Section 2 explains the methods adopted for this work. Section 3 contains a detailed discussion of the obtained results, and Section 4 contains concluding remarks.

The Method and Formalism
For the production of charged hadrons, mini bias Monte Carlo event generators, including PYTHIA8 and parton showers, EPOS-LHC, QGSJETII-04, and Sibyll2.3d, have been used, and the obtained results are contrasted against the CMS data at LHC. For the simulation purpose, the definitions of the kinematic observable have been used in the models. In the momentum space, the transverse momentum is defined as follows: x is the momentum of the particles along x-axis and p 2 y is the momentum of the particles along y-axis. For the pseudorapidity measure, the definition of η = −lntan(θ/2) is used. A total of five hundred thousand events are generated in each case.
PTHIA8 is a general-purpose Monte Carlo event generator for heavy ion hadronic collisions. PYTHIA8 rests on the physics concepts of soft and hard processes, multiparton interactions, parton showers, and string fragmentation. It also links the parton showers, such as Dire showers, Vincia showers, and Simple showers (TimeShowers and SpaceShowers). All three shower models perform ordinary QCD and QED radiation, whereas parton showers provide the baseline for linking external showers [21,22]. In order to tune PYTHIA8 to different parton showers, the selection of mode from 1 to 3 is chosen as PartonShowers : modelmode = 1, 2, 3. For mode = 1, the Simple parton shower is activated; for mode = 2, the Vincia shower is activated; for mode = 3, the Dire shower is activated [23][24][25]. A description of the models' prediction and their comparison to the experimental data in this study is given below.
Cosmic-ray Monte Carlo (CRMC) models (EPOS-LHC, Sibyll2.3d, and QGSJETII-04) are minimum-biased event generators used for the hadron production in hadronic and nuclear interactions at high energies. Moreover, cosmic ray models are very well suited and flexible to describe the hadronic bulk particle production; therefore, they can be used at LHC energies for efficiency corrections. The production of events is very trivial With CRMC models [7,26]. These models are based on Gribov's Reggeon Field Theory (RFT) [27]. It works on the production of particles, and studies the cross-section via the exchange of Pomerons and Reggeons. At higher energies, the contribution of Pomeron dominates through the exchange of multi-gluons and those from secondary Reggeons. Notably, the exchange of soft Pomerons, which are of colorless nature, is mainly responsible for the production of particles in non-central hadronic collisions and diffractive dissociations. The latter is responsible for the total inelastic cross-section up to asymptotic energies.
All the CRMC models are based on soft and (semi)-hard processes of particle production by considering the multi-Pomeron interactions. As a result of high density of the partons at small values of x, the use of the remnants hadrons is different in different generators. The Sibyll2.3d model makes use of an energy-dependent transverse momentum cut-off for the production of mini-jet and relies on the less number of gluons in a hadron for given transverse area, While in EPOS-LHC, the elastic and inelastic fusion of parton-ladders is included to address the nonlinear phenomenological effects [28].
Besides comparing the experimental measurements with the predictions of models, we used the standard distribution to reproduce the experimental data and extract the effective temperature T e f f of the produced medium. The standard distribution is chosen considering the temperature has the closest concept as in the ideal gas model, which includes the chemical potential (µ) and the index of standard distribution (S). This means that the p T -dependent probability density function has the form [29]: In Equation (1), the transverse mass (m T ) is given by where m 0 is used to represent the rest-mass, and N shows the particles' multiplicity. The minimum and maximum rapidity are represented by y minimum and y maximum , respectively. S has a value of 1, while the −ve and +ve values are for bosons and fermions, respectively. Finally, C is used as a constant of normalization, which gives the integral of Equation (1) to be normalized to unity. The parameter µ in Equation (1) is particle-dependent, where i is a particle type. The following equation can show the chemical potential (µ i ) of these particles by [30][31][32]: In Equation (2), z i represents the negative to positive particles' ratio, while T ch in the last equation is the temperature at the chemical freeze-out stage in the statistical thermal model [33][34][35][36][37] and is given by the following expression: where T limit = 0.1580 GeV is the saturation or limiting temperature [38], and √ s NN has the usual meaning of the cms energy in GeV units. One may need to use standard distributions with two-component to fit the p T spectra. A special case is when there is a contribution from the decay of resonances; then, a two-component distribution for a good fit is required. Similarly, if hard scattering has a considerable contribution in the p T region under consideration, a two-component distribution will be a good choice. In the former case, where the decay of resonances contributes is covering a very low-p T region of the spectra (<0.20∼0.30 GeV/c), the total probability density function of the p T of the two components can be written in the form of the superposition principle as follows: In Equation (4), z is the contribution fraction coming from the first component ( . The integral of Equation (4) is normalized to one. Hence, Equation (5) is the weighted average of the two components. The temperature (T) given in Equation (5) is the T e f f , assuming both components stay in equilibrium.

Results and Discussions
We investigate the transverse momentum (p T ) and pseudorapidity (η) distributions of primary charged hadrons in the kinematic range of 0 < p T < 4 GeV/c and |η| < 2.4 at 0.9, 2.36, and 7 TeV. The minimum biased Monte Carlo (MC) event generators, i.e., PYTHIA8.307, EPOS-LHC, QGSJETII-04, and Sibyll2.3d, are used to produce primary charged hadrons, and then the predictions are compared to the CMS data at 0.9 TeV, 2.36 TeV [4], and 7 TeV [5] at the LHC. Since there are three different modules in the PYTHIA8 for parton showers, Simple (Timelike, Spacelike), Vincia, and Dire showers, the predictions of these three modules are first compared to the experimental data [4,5]. The best of the three modules is then used with other MC generators for comparison with the experimental data. A description of the models' prediction compared to the data is given in the following section.

Comparison of Models' Prediction with Data
A detailed description of the model's predictions compared to the CMS data at LHC is outlined here. Figure 1 shows the distributions of the p T of all charged hadrons at the three CMS energies within the kinematic range mentioned above. It is worth mentioning that filled black circles show the experimental data, while lines of different colors represent the model results. The left panels of the figure compare the three modules of Pythia with the experimental data, where red, blue, and green colors are used to display the predictions of the Pythia-Simple, Pythia-Vincia, and Pythia-Dire, respectively. The right panels compare CRMC models and Pythia (Simple module) with the experimental data, where red, blue, green, and yellow show the results of Sibyll2.3d, QGSJETII-04, EPOS-LHC, and Pythia-Simple, respectively. Figure 1a compares parton showers predictions to the CMS data at 0.9 TeV. It is clear that the Simple showers very well describe the data over the entire p T range, while Dire shower predictions are higher than experimental data, which could be due to the kinematics involved in Dire showers and may be due to the nPQCD effects. The Vincia shower predictions are 10% lower than the data at low p T values, but apparent underestimation is seen at high p T values, which goes up to 50%, possibly due to the angular ordering of partons. Figure 1c compares the predictions of the parton showers with the CMS data at 2.36 TeV. It is observed that the behavior of the distributions is the same as in the case of 0.9 TeV energy, while for the right side of the panel, the models' predictions come close to experimental data. At 0 < p T < 1 GeV/c, all the models describe the CMS very well, while for the late p T region, deviation from data is observed, although it is less compared to 0.9 TeV energy. Figure 1e shows that the Dire shower explains the experimental data well over the entire p T range. In contrast, the Simple shower agrees well for up to p T = 2.6 GeV/c. In the later p T region, the models' predictions deviate slightly from experimental data, which could be due to the multi-parton interactions.
While comparing the different modules of the Pythia8 model, we found that although the Dire model predicts better than the other two at 7 TeV, overall, the Simple module reproduced good results. Therefore, in the right panel of Figure 1, we used the predictions of the Simple module of the Pythia8 with Sibyll2.3d, QGSJETII-04, and EPOS-LHC for comparison to the experimental data. Figure 1b compares different models' predictions with the CMS data at 0.9 TeV. One can see that all the models give good descriptions for 0 < p T < 1 GeV/c, while in later p T regions, all models underestimate the experimental data. The deviation is more in the case of the Sibyll2.3d event generator, while the EPOS-LHC and Pythia yielded a good comparison, which is within 10% in most of the cases throughout the p T range. Figure 1d compares the CRMC models and PYTHIA8 (Simple module) predictions with the experimental data. It is clear from the figure that the PYTHIA8 model describes the experimental data with 10% over the entire p T region, while all other models deviate from the experimental data by 30%. Figure 1f shows that the QGSJETII-04 model reproduces the data within 15 % over the entire p T range up to p T = 1 GeV/c. The Sibyll2.3d model has better prediction than the Pythia8 and EPOS-LHC model. The former has a deviation up to 30% while the latter has one up to 40%, mostly in the region of p T ; 1.5 < p T < 5, which may be due to parton fragmentation.  Figure 2 is the pseudorapidity distribution of primary charged hadrons within the −2.4 < η < 2.4 range integrated over transverse momentum at LHC energies. In Figure 2a,b (upper panels), at 0.9 TeV, the Simple showers agree very well with the experimental data, while Dire overestimates and Vincia underestimates the data. The observed deviations may be due to the angular ordering of partons. The PYTHIA8 predictions can also describe the experimental data over the entire η region. On the other hand, Sibyll2.3d, EPOS-LHC, and QGSJETII-04 underestimate the CMS data. Figure 2c,d (middle panels), at 2.36 TeV, shows that non of the parton showers, PYTHIA8, and CRMC models could explain the experimental data. Figure 2e,f (lower panel), at 7 TeV, shows that Dire and Vincia showers do not agree well with the experimental data. On the other hand, Simple shower data are slightly above the experimental data. It is clear from the figure that EPOS-LHC, QGSJETII-04, and Sibyll2.3d are lower than the experimental data, and the PYTHIA8 predictions are slightly higher than the experimental data.

Fit Procedure by Standard Distribution
We used the two-components standard distribution (Equation (4)) to fit the distributions of charged particles' p T of the experimental data in the range of p T , 0 < p T < 4 GeV/c and for pseudorapidity range of |η| < 2.4 at 0.9, 2.36, and 7 TeV. The solid curves resulted from the fit procedure given in Figure 3. The markers with different colors show the experimental data, while lines with the same color as the corresponding data points show the fit results. We used a blue square, red circle, and orange triangle to show the data and a line with blue, red, and orange colors to show the fit curve at 0.9, 2.36, and 7 TeV. The parameters extracted from the fit function are given in Table 1. The χ 2 /nd f values, given in the second last column of the table, show that the function fits the data well. The T 1 values are 0.163 ± 0.008, 0.169 ± 0.008, and 0.182 ± 0.009 at 0.9, 2.36, and 7 TeV, respectively, while the respective values for T 2 are 0.389 ± 0.012, 0.396 ± 0.013, and 0.494 ± 0.013. The value of T, calculated using Equation (5), varies as 0.20368 ± 0.01, 0.22348 ± 0.011, and 0.24128 ± 0.012 for 0.9, 2.36, and 7 TeV, respectively, which shows a growth of T with increasing the center of mass energy. This is due to the fact that more energy is stored in the system as a result of a very violent reaction at higher energies, and correspondingly the system is excited to a high degree. This observation is harmonious with some of our previous results, such as [39][40][41][42], but it disagrees with other results, such as [43,44]. Finally, the parameter N has the values 360 ± 18, 460 ± 23, and 620 ± 31 from the fit procedure for the data at 0.9, 2.36, and 7 TeV, hence showing an increasing trend with the rising energy. N 0 is a free parameter with a special significance in that it reveals the multiplicity. At higher energies, there is a harsh squeeze in the colliding system, where copious particles are produced. Therefore, larger N 0 at higher energies evinces the production of copious particles. It is pertinent to mention that N 0 is different from C. The constant C is used to normalize the integral of Equation (1), while the constant N 0 is used for the comparison of the fit function f (pT) and experimental data. Table 1. Result obtained from the fit function using Equation (1) for the free parameters T 1 , T 2 , T, (N), (χ 2 ), corresponding to the curves in Figure 1. The fit curves resulting from the fit function are given in Figure 3.

Summary and Conclusions
The pseudorapidity and transverse momentum distributions of charged hadrons are presented for the kinematic range of 0 < p T < 4 GeV/c and |η| < 2.4 at 0.9, 2.36, and 7 TeV in pp collisions. The data of the CMS measurements are compared to simulation results using Monte Carlo event generators, i.e., PYTHIA8 and CRMC models. The there modules for parton showers in PYTHIA8 (Simple Showers, Vincia showers, and Dire showers) are analyzed. First, the predictions of the three modules are compared to experimental data, and we found that the Simple showers model yields the best results. In the second step, The prediction of the Simple showers module along with Sibyll2.3d, QGSJETII-04, and EPOS-LHC are compared to experimental measurements. Such studies are very useful to soft parton interactions and hadronization processes and tune event generators, including but not limited to, parton interaction, correlation, hadronization, spin, final state effects, etc. We record our findings of the comparison in the following:

•
The Simple parton showers explain the CMS data very well for p T > 1 GeV/c at 0.9 and 2.36 TeV within the experimental errors, while Dire overshoots and Vicia undershoots the data by 50% each; • At 7 TeV, the Dire module has a good prediction, whereas the Simple and Vincia modules underestimate the data within 30% and 50%; • The comparison of the Simple module of the Pythia model, along with the predictions of the CRMC models with the experimental data, show that at 0.9 TeV, EPOS-LHC has better results than the others; • At 2.36 GeV, CRMC models have better prediction than the Simple module of Pythia at low p T , while QGSJETII-04 predicts well at high p T . Both QGSJETII-04 and EPOS-LHC have closer results than the Pythia-Simple and Sibyll2.3d at 7 TeV; • In the case of the pseudorapidity distributions, only the Pythia-Simple reproduced the experimental measurements at all energies. The Dire module overestimates while Vincia underestimates in decreasing order of discrepancy (20%, 12%, 5%) with energy; • All CRMC models underestimate over the entire η range and for all three energies by 20%; • Furthermore, A fit procedure using the standard distribution resulting in the T e f f for experimental data shows that the T e f f of the hadronic matter increases with an increase in the center of mass energy; • The observed inconsistencies may be linked to the kinematics involved in the simulation of charged hadrons; • In the case of parton showers, the angular ordering of partons could be the possible reason for the deviation; • In contrast, the multi-parton interaction and parton fragmentation could be a possible reason in the case of CRMC models.