Charged-particle production as a function of multiplicity and transverse spherocity in pp collisions at √ s = 5 . 02 and 13 TeV

We present a study of the inclusive charged-particle transverse momentum ( p T ) spectra as a function of charged-particle multiplicity density at mid-pseudorapidity, d N ch / d η , in pp collisions at √ s = 5 . 02 and 13 TeV covering the kinematic range | η | < 0 . 8 and 0 . 15 < p T < 20 GeV/ c . The results are presented for events with at least one charged particle in | η | < 1 (INEL > 0). The p T spectra are reported for two multiplicity estimators covering different pseudorapidity regions. The p T spectra normalized to that for INEL > 0 show little energy dependence. Moreover, the high- p T yields of charged particles increase faster than the charged-particle multiplicity density. The average p T as a function of multiplicity and transverse spherocity is reported for pp collisions at √ s = 13 TeV. For low- (high-) spherocity events, corresponding to jet-like (isotropic) events, the average p T is higher (smaller) than that measured in INEL > 0 pp collisions. Within uncertainties, the functional form of (cid:104) p T (cid:105) ( N ch ) is not affected by the spherocity selection. While EPOS LHC gives a good description of many features of data, PYTHIA overestimates the average p T in jet-like events.


Introduction
Proton-proton collisions at the Large Hadron Collider (LHC) energies have unveiled features very similar to the ones observed in heavy-ion collisions [1]. The previous consensus of the heavy-ion community was that the partonic system created in nuclear collisions needs a large volume to thermalize and to lead to phenomena like collective flow. However, radial [2][3][4] and anisotropic flow [5], as well as strangeness enhancement [6], are also observed in pp and p-A collisions when they are studied as a function of event multiplicity. Surprisingly, with the same level of precision, microscopic and macroscopic approaches describe qualitatively well the observed features in pp collisions. While macroscopic models incorporate hydrodynamical evolution of the system [7], the others include overlapping strings [8], string percolation e-mail: alice-publications@cern.ch [9], multi-parton interactions and color reconnection [10,11]. The multiphase transport model [12], as well as the fragmentation of saturated gluon states [13,14], is able to describe some features of data.
The inclusive transverse momentum ( p T ) spectrum of charged particles carries information of the dynamics of soft and hard interactions. On one hand, the highp T ( p T > 10 GeV/c) particle production is quantitatively well described by perturbative QCD (pQCD) calculations; on the other hand, the understanding of particle production at lowp T has to resort to phenomenological QCD inspired models. Most of the new effects discovered in pp collisions have been unveiled in the low-( p T < 2 GeV/c) and intermediate-(2 ≤ p T < 10 GeV/c) p T domains [2][3][4][5][6]. The present paper reports a novel multi-differential analysis aimed at understanding charged-particle production associated to partonic scatterings with large momentum transfer and their possible correlations with soft particle production.
The transverse momentum distributions are reported for two multiplicity estimators which cover different pseudorapidity regions. The estimators are based on either the total charge deposited in the forward detector (covering the pseudorapidity regions 2.8 < η < 5.1 and −3.7 < η < −1.7) or on the number of tracks in the pseudorapidity region |η| < 0.8. The forward multiplicity estimator is commonly used by the ALICE collaboration to minimize the possible autocorrelations induced by the use of the mid-pseudorapidity estimator. One such examples is the "fragmentation bias" [15], which is the correlation between jet fragments and event multiplicity arising when the particle's p T and event multiplicity are both measured within the same pseudorapidity interval [16]. For each estimator, we defined different multiplicity classes based on either the number of tracks at mid-pseudorapidity (|η| < 0.8) or the signal in the forward detectors. It is worth mentioning that a similar study has been performed by ALICE using p-Pb data; the results showed different modifications of the spectral shapes depending on the multiplicity estimators which were used [17]. To disentangle the energy and multiplicity dependence, for a given 0123456789().: V,-vol 123 multiplicity class, the p T distributions are measured for pp collisions at √ s = 5.02 and 13 TeV. Particle production from intermediate to high p T (> 4 GeV/c) is studied by fitting a power-law function to the invariant yield, and studying the multiplicity and energy dependence of the exponent. This has been proposed in Ref.
[18] as a way to characterize the highp T tails of different systems and energies in a convenient way that may make the comparison for the different systems more straightforward.
Finally, we explore a new approach, which has been proposed to study multi-parton interaction effects in pp collisions. Transverse spherocity [19], hereinafter referred to as spherocity, has been proven to be a valuable tool to discriminate between jet-like and isotropic events [20] associated with an underlying event activity which is either suppressed or enhanced. The previous measurement of average transverse momentum of inclusive charged particles as a function of event multiplicity [21] is now explored adding a new dimension: the event shape characterized by spherocity. The aim of this study is to investigate the importance of jets in high-multiplicity pp collisions and their contribution to charged-particle production at low p T .
The paper is organized as follows: Sect. 2 describes the run conditions during the data taking and the main detectors used in the present analysis. Section 3 outlines the analysis details for the event and track selection, as well as the definitions of the different event classes. The correction procedures and the estimation of the systematic uncertainties are summarized in Sects. 4 and 5, respectively. Results and discussions are presented in Sect. 6. Finally, our summary and conclusions are reported in Sect. 7.

The ALICE apparatus
The main detectors used in the present work are the Inner Tracking System (ITS), the Time Projection Chamber (TPC) and the V0 detector. The ITS and TPC detectors are both used for primary vertex and track reconstruction. The V0 detector is used for triggering and for background rejection. More details concerning the full ALICE detector system can be found in Ref. [22].
The central barrel covers the pseudorapidity region |η| < 0.8 for full-length tracks. The main central-barrel tracking devices are the ITS and the TPC, which are located inside a solenoid magnet providing a 0.5 T magnetic field allowing the tracking of particles from 0.15 GeV/c. The ITS is composed of six cylindrical layers of high-resolution silicon tracking detectors. The innermost layers consist of two arrays of hybrid Silicon Pixel Detectors (SPD) located at an average radial distance (r ) of 3.9 and 7.6 cm from the beam axis and covering |η| < 2 and |η| < 1.4, respectively. The SPD is also used to reconstruct tracklets, which are track seg-ments built using the position of the reconstructed primary vertex and two hits, one on each SPD layer. The number of tracklets gives an excellent estimate of the charged-particle multiplicity at mid-pseudorapidity (N ch ). The outer layers of the ITS are composed of silicon strip and drift detectors, with the outermost layer sitting at r = 43 cm. The TPC is a large cylindrical drift detector of radial and longitudinal size of about 85 < r < 250 cm and −250 < z < 250 cm, respectively. It is segmented in radial "pad rows", providing up to 159 tracking points. The measurement of charged particles is based on "global tracks", reconstructed using the combined ITS and TPC information. The V0 detector consists of two forward scintillator arrays (V0-A and V0-C) employed for triggering, background suppression, and event-class determination. They are placed on either side of the interaction region at z = 3.3 m and z = −0.9 m, covering the pseudorapidity regions 2.8 < η < 5.1 and −3.7 < η < −1.7, respectively.
The data were collected using a minimum-bias trigger which required coincident signals in both V0-A and V0-C detectors . The events were recorded in coincidence with signals from two beam pick-up counters each positioned on either side of the interaction region to tag the arrival of proton bunches from both directions. Control triggers taken for various combinations of beam and empty buckets were used to measure beam-induced and accidental backgrounds. The contamination from background events was removed offline by using the timing information from the V0 detector, which has a time resolution better than 1 ns. Background events were also rejected by exploiting the correlation between the number of clusters of pixel hits and the number of tracklets in the SPD.

Analysis
The results presented here were obtained from the analysis of about 105 and 60 million minimum-bias pp events at √ s = 5.02 and 13 TeV, respectively. The interaction probability per single bunch crossing ranges between 2% and 14% for pp collisions at 13 TeV and from 0.3% to 6% for pp collisions at 5.02 TeV. The measurements have been obtained for events having at least one charged particle produced in the pseudorapidity interval |η| < 1 (INEL > 0). For the analysis, the events were furthermore required to have a reconstructed vertex located within |z| < 10 cm, where z is the position of the vertex along the beam axis, and z = 0 cm corresponds to the nominal center of the detector [22]. Events containing more than one distinct vertex were tagged as pileup and discarded from the analysis. The systematic uncertainty associated to pileup is between 3 -4% and is not the dominant source of uncertainty for the p T spectra reported here. The corrections are calculated using Monte Carlo events from the PYTHIA 6 [23] (tune Perugia 2011 [24]) event generator with particle transport performed via a GEANT 3 [25] simulation of the ALICE detector.
Only primary charged particles in the kinematic range |η| < 0.8 and 0.15 < p T < 20 GeV/c are considered in the transverse momentum analysis. A primary charged particle is defined as a charged particle with a mean proper lifetime τ larger than 1 cm/c, which is either produced directly in the interaction or from decays of particles with τ smaller than 1 cm/c, excluding particles produced in interactions with the detector material [26].
Transverse momentum distributions The measurement of the p T spectra follows the standard procedure already employed in several ALICE publications [27][28][29]. Tracks reconstructed using the information from the ITS and TPC detectors are used. The track selection criteria have been optimised for best track quality and minimal contamination from secondary particles. Tracks are required to have at least two hits in the ITS detector, of which at least one is in either of the two innermost SPD layers. The geometrical track length L (in cm) is calculated in the TPC readout plane, excluding the information from the pads at the sector boundaries (∼3 cm from the sector edges). The number of crossed TPC rows has to be larger than 0.85L. The number of TPC clusters has to be larger than 0.7L. The fit quality for the ITS and TPC track points must satisfy χ 2 ITS /N hits < 36 and χ 2 TPC /N clusters < 4, respectively, where N hits and N clusters are the numbers of hits in the ITS and the number of clusters in the TPC, respectively. Tracking information from the combined ITS and TPC track reconstruction algorithm is compared to that derived only from the TPC and constrained by the interaction vertex point. Then, the quantity χ 2 TPC−ITS is derived as described in Ref. [30]. Only tracks with χ 2 TPC−ITS < 36 are included in the analysis in order to improve the purity of primary track reconstruction at high p T . Tracks are rejected if their distance of closest approach to the reconstructed vertex in longitudinal and radial direction, d z and d xy , respectively, satisfies d z > 2 cm or d xy > 0.018 cm + 0.035 cm × p −1.01 T , with p T in GeV/c.

Multiplicity estimators
In order to study the multiplicity dependence of the inclusive charged particle p T distributions, the INEL > 0 sample is divided into event classes based on the total charge deposited in the V0 detector (V0M amplitude) and on the number of SPD tracklets (N SPD tracklets ) in the pseudorapidity region |η| < 0.8. The event classes used in the analysis and the corresponding mid-pseudorapidity charged particle densities are summarized in Tables 1 and 2. The average charged-particle multiplicity densities for INEL > 0 collisions and for the multiplicity classes are obtained by integrating the corresponding fully corrected p T spectra (measured using ITS and TPC information). To this end, the p T spectra were extrapolated to p T = 0 with a Hagedorn function [31]. Different functions were used and the differences with respect to the reference values were considered in the systematic uncertainties. For INEL > 0 pp collisions at √ s = 5.02 TeV the mid-pseudorapidity (|η| < 0.8) charged-particle density is dN ch /dη = 5.91 ± 0.45, while for √ s = 13 TeV the corresponding value is 7.60±0.50. The comparison of results obtained with these estimators allows to understand potential biases from measuring the multiplicity and p T distributions in overlapping η regions.
Spherocity For the data analysis we followed a strategy similar to that already reported in Ref. [32]. Spherocity, S 0 , originally proposed here [33] is defined for a unit vectorn s which minimizes the ratio: where the sum runs over all reconstructed ITS-TPC tracks. At least three tracks are required within |η| < 0.8 and p T > 0.15 GeV/c in order to achieve a good spherocity resolution. The spherocity resolution improves with the trackreconstruction efficiency, therefore the restrictions on the purity of primary charged particles can be relaxed. For spherocity we considered all tracks with at least 50 clusters in the TPC, which satisfy: d xy < 2.4 cm, d z < 3.2 cm, and χ 2 TPC /N clusters < 4. The exclusion of the ITS requirements guarantees a homogeneous azimuthal track-reconstruction efficiency.
It is worth mentioning some important features of spherocity: -The vector products are linear in particle momenta, therefore spherocity is a collinear safe quantity in pQCD. -The lower limit of spherocity (S 0 → 0) corresponds to event topologies where all transverse momentum vectors are (anti)parallel or the sum of the p T is dominated by a single track. -The upper limit of spherocity (S 0 → 1) corresponds to event topologies where transverse momentum vectors are "isotropically" distributed. S 0 = 1 can only be reached in the limit of an infinite amount of particles.
Since the goal of the present study is to separate jet events from isotropic ones, we study different spherocity classes for a given multiplicity value. The multiplicity is measured by counting the number of tracks within |η| < 0.8. As explained later, we adopted the procedure used in the analysis of average p T as a function of multiplicity to correct the number of tracks for detector effects [21]. The detector response is represented by a two-dimensional distribution: reconstructed spherocity as a function of generated spherocity, each bin of generated spherocity is normalized to unity. In this representation, the two-dimensional distribution gives 123 the normalized response matrix R (S 0 , S m ), which contains the probability that an event with spherocity S 0 is reconstructed with spherocity S m . Figure 1 shows the spherocity response matrices for two track multiplicity (N m ). Tracking efficiency effects on the spherocity resolution are relevant only for low-multiplicity events, therefore, the S 0 resolution improves with increasing multiplicity. In order to study the spherocity dependence of the particle production for a given track multiplicity value, the sample is divided into ten event sub-classes of equal size (percentiles), Page 5 of 22 857 The detector response for spherocity for two track multiplicity classes: dN m /dη = 3.12 (left) and dN m /dη = 25 (right). Protonproton collisions were simulated using PYTHIA 6. The simulations include the particle transport performed via a GEANT 3 simulation of the ALICE detector. The markers (boxes around the points) indicate the average (RMS) of the measured spherocity distributions for each bin of spherocity at generator level. The spherocity binning varies with dN m /dη, because in this way, it allows the analysis of ten event subclasses of equal size. The probability that an event with spherocity S 0 be reconstructed with spherocity S m is represented by P(S m ) based on the measured spherocity distribution. From now on, the most jet-like and isotropic events will be referred to as 0 -10% and 90 -100% spherocity event class, respectively.
It has been reported that the evolution of several observables as a function of center-of-mass energy can be factored out to be due to the changes in charged-particle multiplicities which in turn depend on the energy. For example, the particle production sensitive to the underlying event for different √ s exhibits approximate scaling properties connected to changes in N ch [34]. Moreover, within uncertainties, the average p T as a function of multiplicity exhibits a small energy dependence [21]. Therefore, the spherocity dependent average p T as a function of charged-particle multiplicity is only presented for pp collisions at √ s = 13 TeV. The physics message is valid for other center-of-mass energies, this was verified using data from pp collisions at √ s = 5.02 TeV.

Corrections
All the measurements presented in this paper are fully corrected for acceptance and tracking efficiency, contamination from secondary particles, event and signal loss, as well as multiplicity and spherocity resolution. Details of these corrections are presented below.

Transverse momentum distributions as a function of particle multiplicity
The transverse momentum spectrum for a specific event class is obtained by correcting the track yields N rec reconstructed in each ( η, p T ) interval for all detector effects that either influence the event reconstruction or the track reconstruction. The transverse momentum distribution is obtained as follows: The event selection (for a specific event class) and vertex reconstruction efficiencies are represented by ev. class and vz , respectively. The number of events of a given event class is represented by N rec ev . For the lowest multiplicity class selected using the V0M amplitude and for √ s = 5.02 TeV ( √ s = 13 TeV) they reach 66% and 95% (75% and 95%), respectively, while for the highest multiplicity class the detector is fully efficient. The track-level correction factors, C( η, p T ), are obtained for events which satisfy the selection criteria; they include acceptance, efficiency, purity, and p T resolution. The estimation of the four terms will be explained in detail in the following.
A data-driven method has been developed to reduce the systematic uncertainty related to incorrect description of the particle composition in Monte Carlo. The tracking efficiency is determined using the re-weighting procedure which is discussed for the first time in Ref. [29] and which is employed also in the present paper. The method uses the knowledge of the particle composition at LHC energies, i.e. the abundances of the different particle species within a specific interval of p T and for a specific event class.
To correct the distributions for secondary-particle contamination, i.e. the products of weak decays of kaons and baryons, and the particles originating from interactions in the detector material, we used the d xy distributions of particles in data and Monte Carlo simulations. Exploiting the differ-123 ences of the d xy distributions between primary and secondary particles, especially in the tails, the measured distributions were fitted by a linear combination of d xy distributions (templates) for primary and secondary particles obtained from Monte Carlo simulations in different p T bins. For INEL > 0 pp collisions at √ s = 13 TeV the contamination ranges from 8.5% at p T = 0.2 GeV/c to 1% for p T > 2 GeV/c. The contamination exhibits a small multiplicity dependence, which is below 2%. For pp collisions at √ s = 5.02 TeV, the correction factors reach similar values.
The transverse momentum spectra are also corrected for p T resolution; the correction factor is calculated using the covariance matrix of the Kalman fit [35]. The p T (multiplicity) dependence of the correction factor is negligible [29] (below 1%).
Finally, the p T spectra are corrected for the amount of signal which is missing from the yield due to the event selection (signal loss). This correction is negligible for highmultiplicity events and reaches 13% (4%) at p T = 0.2 ( p T = 1) GeV/c for the lowest multiplicity class based on N SPD tracklets .

Spherocity studies
The measurement of the average transverse momentum as a function of charged-particle multiplicity and spherocity is performed following a strategy close to that used in earlier publications [21,36]. The transverse momentum spectra for different multiplicity and spherocity classes are fully corrected as described in the previous section. The average transverse momentum is then calculated from the corrected spectra as the arithmetic mean in the kinematic range 0.15 < p T < 10 GeV/c and |η| < 0.8.
To extract the correlation between p T and the number of primary charged particles (N ch ) in |η| < 0.8 and for the spherocity class S 0 , the following re-weighting procedure is applied to account for the experimental resolution of the measured event multiplicity (N m ) and spherocity (S m ): This method is an extension to the one developed for the previous p T analysis [36]. It exploits the normalized response matrices R and R which encode the multiplicity, and spherocity detector resolutions, respectively. The average p T for the S 0 event class is encoded inside the inner sum, where the weights R (S 0 , S m ) are explicitly applied to p T values. The resulting p T (N m , S 0 ) is then corrected for multiplicity resolution. It is worth mentioning that the spherocity-integrated class (0 -100%) only requires the mul-tiplicity correction. The Monte Carlo non-closure, discussed in the next section, is assigned as systematic uncertainty.

Transverse momentum spectra
The relative systematic uncertainties on the p T spectra are summarized in Table 3. They include the effect of the event selection based on the vertex position, which is studied by comparing the fully corrected p T spectra obtained with alternative vertex selections: |z| < 5 cm and |z| < 20 cm. The corrections due to trigger and vertex selection were determined using the EPOS LHC [37] event generator and the deviations with respect to the nominal values, i.e. those obtained with PYTHIA 6, were assigned as systematic uncertainties. The same procedure was employed for the estimation of the systematic uncertainty associated to the signal loss correction. The systematic uncertainty related to the track selection was studied by varying the track cuts for which we used the variation intervals described in Ref. [29]. We also studied the systematic effects related to the uncertainty on the primary particle composition which is assumed for the efficiency correction. This uncertainty takes into account the extrapolation of the spectra to low p T , the relative particle abundances at high p T , the uncertainties of the measured particle spectra, and the Monte Carlo assumptions on the ± / spectra ratios. The systematic uncertainties of the correction for secondaries contamination is estimated by varying the fit model using two templates, i.e. for primaries and secondaries, or three templates, i.e. primaries, secondaries from interactions in the detector material, and secondaries from weak decays, as well as varying the fit momentum ranges. Since we are using the same event selection and track cuts as those used in Ref.
[29], the systematic uncertainties associated with matching efficiency, p T resolution and material budget, are identical.

Average transverse momentum
A summary of the systematic uncertainties for three multiplicity values and for different spherocity classes is shown in Table 4. In order to estimate the systematic uncertainties of p T , the results of the data analysis and of the evaluation of the corrections from Monte Carlo simulations were studied considering cut variations and Monte Carlo assumptions, within reasonable limits. The effect of the track cuts on p T was found to be spherocity independent and of the order of 1%. The efficiency correction is another spherocity independent contribution and it is found to be ∼1%. This contribution takes into account the different particle composition in data and models, as well as the multiplicity dependence of the correction. We also studied the multiplicity dependence of the purity correction; the effect was found to be smaller than 0.5%. The most relevant spherocity independent contribution is related to the re-weighting procedure to correct for the detector multiplicity resolution. This was quantified from the Monte Carlo non-closure, it amounts to ∼1.36%, ∼0.86% and ∼1.26% for dN ch /dη = 1.88, 6.25, 25.00, respectively. The set of track cuts used to measure spherocity was also varied compared to those used for the p T spectra analysis. The effect on the results amounted to 1%. The most relevant contribution to the systematic uncertainties originates from the re-weighting procedure method which is used to correct for the spherocity resolution. The Monte Carlo non-closure is assigned as a systematic uncertainty. For the lowest multiplicity value, dN ch /dη = 1.88, the uncertainty reaches 3.23%, 4.55%, and 7.06% for the 0 -10%, 40 -50%, and 90 -100% spherocity classes, respectively. For higher multiplicities, e.g. dN ch /dη = 25.0, the Monte Carlo nonclosure amounts to 0.57%, 1.07%, and 2.01% for the 0 -10%, 40 -50%, and 90 -100% spherocity classes, respec-123 tively. As expected from the detector response, the most relevant effects are observed for low-multiplicity events in particular for the isotropic classes. As will be seen later, in Monte Carlo jet-like events, the average p T shows a strong change with multiplicity at dN ch /dη ∼ 7. This effect increases the size of the uncertainty (Monte Carlo non-closure) in that multiplicity interval. This is the dominant contribution to the systematic uncertainties and covers the largest variations observed between data and PYTHIA 8 (version 8.212) [10] (tune Monash 2013 [38]).
The model dependence is also checked by using events simulated with PYTHIA 8 and EPOS LHC which include the particle transport through the detector. The corrections were calculated using these simulations and the maximum variation with respect to the nominal values (using PYTHIA 6 simulations) are below 1%.

Results
6.1 Transverse momentum spectra as a function of charged-particle multiplicity The p T distributions of charged particles, measured in |η| < 0.8 for pp collisions at √ s = 5.02 and 13 TeV, are shown in Fig. 2 for the different multiplicity classes selected using the estimator based on N SPD tracklets . The bottom panels depict the ratios to the p T distribution of the INEL > 0 event class. The features of the spectra, i.e. the change of the spectral shape going from low-to high-multiplicity values, are qualitatively the same for both energies. The only significant difference is the multiplicity reach which is higher at 13 TeV than that at 5.02 TeV. In the following we discuss the results for pp collisions at the highest energy. As shown in Fig. 2, the p T spectra become harder as the multiplicity increases, which contributes to the increase of the average transverse momentum with multiplicity. The ratios to the INEL> 0 p T distribution exhibit two distinct behavior. While at low p T (< 0.5 GeV/c) the ratios exhibit a modest p T dependence, for p T > 0.5 GeV/c they strongly depend on multiplicity and p T . Figure 3 shows the multiplicity dependent p T spectra using a multiplicity selection based on the V0M amplitude. Results for pp collisions at √ s = 5.02 and 13 TeV are shown. The average multiplicity values are significantly smaller than those reached with the mid-pseudorapidity estimator (based on N SPD tracklets ). For example, in pp collisions at √ s = 13 TeV, while the average charged-particle multiplicity density amounts to 56.55 for the highest N SPD tracklets class, it only reaches 27.61 for the highest V0M multiplicity class. We note that for similar average particle densities, e.g. the multiplicity classes II (V0M) and VII' (SPD tracklets) in pp collisions at √ s = 13 TeV, the ratios measured using Total syst. unc. ) 1 10 × X' ( ) 2 10 × IX' ( ) 3 10 × VIII' ( ) 4 10 × VII' ( ) 5 10 × VI' ( SPD tracklets mult. classes = 13 TeV s ) 6 10 × V' ( ) 7 10 × IV' ( ) 8 10 × III' ( ) 9 10 × II' ( ) 10 10 × I' ( Stat. unc.
Uncorr. syst. unc. 1 1 0 TeV are shown in the left and right panels, respectively. Statistical and total systematic uncertainties are shown as error bars and boxes around the data points, respectively. In the middle panels, ratios of multiplicity dependent spectra to INEL > 0 are shown in logarithmic scale. In the bottom panels we show the ratios in a linear scale to illustrate the dramatic behavior of the ratios. The systematic uncertainties on the ratios are obtained by considering only contributions uncorrelated across multiplicity. The spectra are scaled to improve the visibility the V0M amplitude and the N SPD tracklets are similar. The comparison of the p T spectra for these multiplicity classes is shown in Fig. 4. We observe that for transverse momentum below 0.5 GeV/c, the spectra exhibit the same shape. For transverse momenta within 0.5-3 GeV/c the spectra for the multiplicity class II is harder than that for the VII" class. At higher p T , the spectral shapes are the same, but the yield of the class II is ∼15% higher than that for the VII' class. Similar results are obtained if we compare the multiplicity classes I and VI' for pp collisions at 5.02 TeV.
Commonly, the particle production is characterized by quantities like integrated yields, or any fit parameter of the curve extracted from fits to the data, for example, the socalled inverse slope parameter reported by ALICE in Ref.
[39]. This facilitates the visualization of the evolution of the particle production as a function of multiplicity and the com- Total syst. unc. ) 1 10 × X ( ) 2 10 × IX ( ) 3 10 × VIII ( ) 4 10 × VII ( ) 5 10 × VI ( V0 mult. classes = 13 TeV s ) 6 10 × V ( ) 7 10 × IV ( ) 8 10 × III ( ) 9 10 × II ( ) 10 10 × I ( Stat. unc. Uncorr. syst. unc. 1 1 0 It is interesting and important to define a common quantity to compare the shape of the highp T part of the spectra of different particle species and collision systems. The natural choice is fitting a power-law function (α× p −n T ) to the invariant yield and studying the multiplicity dependence of the exponent (n) extracted from the fit. Figure 5 illustrates the results considering particles with transverse momentum within 6-20 GeV/c for pp at √ s = 13 TeV. It is worth mentioning that within uncertainties the power-law function describes rather well the data in that p T interval. Similarly, the p T spectra simulated with the different generators are well described (within 2%) by the power-law function.
Within uncertainties, going from low to high multiplicity n decreases taking values from 6 to 5, respectively. A similar behavior has been reported for heavy-ion collisions [41]. Moreover, the results using the two multiplicity estimators are consistent within the overlapping multiplicity interval. This result is consistent with that shown in Fig. 4. PYTHIA 6 and 8 simulations describe the trends very well, but a strong deviation between EPOS LHC and data is observed. In PYTHIA 8, it has been shown that the number of highp T jets increases with event multiplicity. Moreover, for a given event multiplicity and fixed jet p T , the highp T tails of the charged-particle spectra are very similar in low-and high-multiplicity events [16]. Therefore, based on PYTHIA 8 studies, the reduction of the power-law exponent with increasing multiplicity can be attributed to an increasing number of highp T jets.
As pointed out above, the ratios to the INEL > 0 p T distributions for dN ch /dη 25 exhibit a weak p T -dependence for p T > 4 GeV/c. This applies to both energies and to all multiplicity estimators. To illustrate better the behaviour of the yields at high momenta, we adopted a representation previously used for heavy-flavour hadrons [42] to point out to the similarities between the two results. The trend at highp T is highlighted in Fig. 6, which shows the integrated yields for three transverse momentum intervals (2 < p T < 10 GeV/c, 4 < p T < 10 GeV/c, and 6 < p T < 10 GeV/c) as a func-    Fig. 6 Self-normalized yields of charged particles integrated over different p T intervals: 2 < p T < 4 GeV/c, 4 < p T < 6 GeV/c, and 6 < p T < 10 GeV/c. The integrated yields for pp collisions at √ s = 13 TeV are shown as a function of charged-particle density at mid-pseudorapidity. Statistical and uncorrelated (across multiplicity) systematic uncertainties are shown as error bars and boxes around the data points, respectively. Statistical uncertainties are negligible compared to systematic uncertainties. Data are compared with PYTHIA 8 (dashed lines) and EPOS LHC (solid lines). The dotted line is drawn to see the differences between data and the linear dependence. Deviations of data from the linear trend are shown in the bottom panel tion of the average mid-pseudorapidity multiplicity. Both the charged-particle yields and the average multiplicity are selfnormalized, i.e. they are divided by their average value for the INEL > 0 sample. The highp T (> 4 GeV/c) yields of charged particles increase faster than the charged-particle multiplicity, while the increase is smaller when we consider lowerp T particles. The trend of the data is qualitatively well reproduced by PYTHIA 8, but for p T > 6 GeV/c the model significantly overestimates the ratio by a factor larger than 1.5. Although the shapes of the spectra (characterized by n) are not well reproduced by EPOS LHC, the model gives the best description of the self-normalized yields. Despite the large uncertainties, it is clear the data show a non-linear increase.

Double-differential study of the average transverse momentum
The spherocity-integrated average p T as a function of dN ch /dη for pp collisions at √ s = 13 TeV is shown in Fig. 7. In accordance with measurements at lower energies [21], the p T increases with dN ch /dη. In PYTHIA 8 the effect is enhanced by color reconnection, which allows the interaction among partons originating from multiple semi-hard scatterings via color strings. The minimum-bias data are compared with analogous measurements for the most jet-like structure (0 -10%) and isotropic (90 -100%) event classes. Studying observables as a function of spherocity reveals interesting features. On one hand, for isotropic events the average p T stays systematically below the spherocity-integrated p T over the full multiplicity range; on the other hand, for jet-like events the p T is higher than that for spherocity-integrated events. Moreover, within uncertainties the overall shape of the correlation, i.e. a steep linear rise below dN ch /dη = 10 followed by a less steep but still linear rise above, is not spherocity-dependent. Figure 8 shows that within uncertainties, PYTHIA 8 with color reconnection gives an adequate description of the spherocity-integrated event class. It is worth mentioning that color reconnection was originally introduced to explain the rise of p T with multiplicity [43]. However, PYTHIA 6 shows a steeper rise of p T with dN ch /dη than that seen in data. The Perugia 2011 tune relies on Tevatron and SPS minimum-bias data, while the Monash tune was constrained using the early LHC measurements [38]. The comparison of data with EPOS LHC is also shown. Clearly, the quantitative agreement is as good as that achieved by PYTHIA 8. The EPOS LHC model uses a different approach in order to simulate the hadronic interactions. Namely, the model considers a collective hadronization which depends only on the geometry and the density [37].
For the 0 -10% and 90 -100% spherocity classes, Fig. 8 also shows comparisons between data and Monte Carlo generators (PYTHIA 6, PYTHIA 8 and EPOS LHC). It is worth mentioning that we also used spherocity percentiles in all the Monte Carlo event generators reported in this paper because their spherocity distributions do not differ much from those measured in data. For further Monte Carlo comparisons the spherocity binning which was used in the analysis is provided as HEP data. In low-multiplicity events (dN ch /dη < 10), the deviations between data and PYTHIA 8 (without color reconnection) are smaller and larger respectively for the 0 - color reconnection, as well as PYTHIA 6 and EPOS LHC are displayed. Statistical uncertainties (error bars) are negligible compared to systematic uncertainties (shaded area around the data points). Data to model ratios are shown in the bottom panel. The color band around unity represents the systematic uncertainty 10% and 90 -100% spherocity classes than those seen for the 0 -100% spherocity class. The effect could be a consequence of the reduction of color reconnection contribution in events containing jets surrounded by a small underlying event activity. For isotropic events the three models quantitatively describe the correlation. Even for PYTHIA 6, the size of the discrepancy which was pointed out for the spherocityintegrated event class is reduced. On the contrary, for jet-like events both PYTHIA 6 and 8 exhibit a larger disagreement with the data. These models produce three distinct multiplicity regions, for dN ch /dη 7 the models give a steeper rise of p T than data. Within the intermediate multiplicity interval (7 dN ch /dη 25), the slope of p T given by models is more compatible with that seen in data, although the models overestimate the average p T . While in data the average p T increases at a constant rate with multiplicity for dN ch /dη 7, PYTHIA 6 and 8 shows a third change of the slope of p T , observed for dN ch /dη 25. The data to model ratio indicates a discrepancy larger than 10%, which is larger than the systematic uncertainties associated to p T in that multiplicity interval.
In order to study the details of the changes of the functional form of p T (N ch ) due to the spherocity selection, Fig. 9 shows the average p T of jet-like and isotropic events normalized to that for the spherocity-integrated event class. For jet-like events, the data exhibit a hint of a modest peak at dN ch /dη ∼ 7, which is not significant if we consider the size of the systematic uncertainties. Moreover, within uncertainties the ratio remains constant for dN ch /dη 25. EPOS LHC describes rather well the high-multiplicity behavior, however, it overestimates the peak. PYTHIA 6 and 8 show the worst agreement with the data. In this representation, the three distinct regions, which were described before are highlighted. In PYTHIA 8, the peak (at dN ch /dη ∼ 7) in jetlike events is caused by particles with transverse momentum above 2 GeV/c. The size of the peak is determined by particles with p T > 5 − 6 GeV/c. In contrast, data do not show a significant peak structure for any specific transverse momentum interval. We also varied the upper p T (0.15 < p T < p max T ) limit ( p max T = 10 GeV/c is the default) and studied the effect on the extracted p T . The p T remains constant within uncertainties for 4 < p max T < 10 GeV/c in data and for 6 < p max T < 10 GeV/c in PYTHIA 8. For p max T = 2 GeV/c the p T decreases by 23% (29%) in data (PYTHIA 8) compared to p max T = 10 GeV/c. The relative difference of p T between data and PYTHIA 8 amounts to 9% (4%) for p max T = 2 GeV/c ( p max T = 10 GeV/c). The results suggest that the power-law tail produces a smaller impact on data than in PYTHIA 8. A similar ratio for isotropic events shows a smaller structure at dN ch /dη ∼ 7. This effect is well reproduced by all models.
Finally, we also examined the evolution of p T (N ch ) going from the most jet-like to the most isotropic event classes. Figure 10 shows the spherocity-dependent p T (N ch ) in data and models, the data to model ratios are displayed in Fig. 11. The difference between the 0 -10% and 10 -20% spherocity classes is smaller for data and EPOS LHC than for PYTHIA 6 and 8. Moreover, within uncertainties PYTHIA 8 describes rather well the data for the 10 -20% spherocity class. This contrasts with the disagreement between the model and data for the 0 -10% spherocity class. Other features in PYTHIA 6 and 8 are the reduction of the bump at dN ch /dη ∼ 7 and the disappearance of a third rise of the p T for dN ch /dη 25 when one goes from the 0 -10% to the 10 -20% spherocity classes. The agreement among models and data for the 20 -100% spherocity classes is similar to that observed for the 10 -20% spherocity class. Within uncertainties, PYTHIA 8 and EPOS LHC qualitatively describe the data for dN ch /dη 10, while PYTHIA 6 overestimates the average p T .
From previous LHC studies we know that the production cross section of jets in high-multiplicity pp collisions is smaller in data than predicted from the Monte Carlo generators [32,44,45]. Therefore, a possible interpretation is that the low-momentum partons, color connected with higher momentum ones (jets), would produce an overall increase of the hadron transverse momentum. This would affect more the lowp T part of the spectrum associated with jet-enriched samples, which are achieved by requiring low-spherocity values. The incorporation of these new observables into the PYTHIA 8 tuning could be a challenge because, on one hand, the color reconnection has to be reduced to describe the low-S 0 data; on the other hand, the variation should not be too large because the good description of the spherocityintegrated and isotropic classes could be affected.

Summary and conclusions
In this paper, we have reported the transverse momentum spectra of inclusive charged particles in pp collisions at √ s = 5.02 and 13 TeV. The measurements were performed in the kinematic range of |η| < 0.8 and p T > 0.15 GeV/c. The particle production was studied as a function of event multiplicity quantified by two estimators, one based on the number of SPD tracklets within |η| < 0.8, and the second one based on the multiplicity in the V0 forward detector (V0M amplitude). For similar average charged-particle densities, the particle production above p T = 1 GeV/c is higher in pp collisions at √ s = 13 TeV than at √ s = 5.02 TeV. For a fixed center-of-mass energy, particle production above p T = 0.5 GeV/c exhibits a remarkable multiplicity dependence. Namely, for transverse momenta below 0.5 GeV/c, the ratio of the multiplicity dependent spectra to those for INEL > 0 pp collisions is rather constant, and for higher momenta, it shows a significant p T dependence. The behavior observed for each of the two multiplicity estimators are consistent within the dN ch /dη interval defined by the V0M multiplicity estimator, which gives a dN ch /dη reach of ∼25. For the highest V0M multiplicity class, the ratio increases going from p T = 0.5 GeV/c up to p T ≈ 4 GeV/c, then for higher p T , it shows a smaller increase.
The particle production at high transverse momenta is characterized by the exponent of a power-law function which is fitted to the invariant yield considering particles with 6 < p T < 20 GeV/c. Within that p T interval, the power-law function describes rather well the p T spectra. In concordance to the ratios discussed above, within uncertainties, the functional form of n as a function of dN ch /dη is the same for the two multiplicity estimators used in this analysis. Moreover, n is found to decrease with dN ch /dη . Within uncertainties, PYTHIA 8 (tune Monash 2013) and PYTHIA 6 (tune Perugia 2011) quantitatively reproduce the behavior of data, while EPOS LHC overestimates the value of the exponent. Nevertheless, all models describe the self-normalized highp T yields as a function of self-normalized charged-particle multiplicity. Finally, the measurement of the average transverse momentum as a function of event multiplicity at midpseudorapidity was presented. The results for the spherocityintegrated class (nearly identical to INEL > 0 pp collisions) at √ s = 13 TeV are consistent with previous measurements at lower energies. The increase of the average p T with increasing multiplicity is well captured by PYTHIA 8 and EPOS LHC. In order to get a better insight into the particle production mechanisms, the spherocity-integrated sample was separated into different sub-classes characterized by the event structure in the transverse plane. Jet-like and isotropic events were selected based on the spherocity of the events. Isotropic events are well described by the three models which were considered in this work. Interestingly, PYTHIA 6 reproduces these event classes better than the INEL > 0 sample. For jet-like events, the average p T is overestimated by PYTHIA 6 and 8 models in the full multiplicity interval reported. However, EPOS LHC gives the best description of the jet-like event subsample.
The results presented in this paper illustrate the difficulties for the models to describe different observables once they are differentially analyzed as a function of several variables. The measurements are important to better understand the similarities between heavy-ion and small collision systems, as well as for Monte Carlo tuning purposes.