2020 059 15 April

Production of light-flavor hadrons in pp collisions at $\sqrt{s}~{}=~{}7\text{
and }\sqrt{s}=13\text{\,TeV}$

ALICE Collaboration††thanks: See Appendix A for the list of collaboration
members ALICE Collaboration

The production of $\pi^{\pm}$, $\rm{K}^{\pm}$, $\rm{K}^{0}_{S}$,
$\rm{K}^{*}(892)^{0}$, $\rm{p}$, $\phi(1020)$, $\Lambda$, $\Xi^{-}$,
$\Omega^{-}$, and their antiparticles was measured in inelastic proton-proton
(pp) collisions at a center-of-mass energy of $\sqrt{s}$ = 13 TeV at
midrapidity ($|y|<0.5$) as a function of transverse momentum ($p_{\rm{T}}$)
using the ALICE detector at the CERN LHC. Furthermore, the single-particle
$p_{\rm{T}}$ distributions of $\rm{K}^{0}_{S}$, $\Lambda$, and
$\overline{\Lambda}$ in inelastic pp collisions at $\sqrt{s}$ = 7 TeV are
reported here for the first time. The $p_{\rm{T}}$ distributions are studied
at midrapidity within the transverse momentum range $0\leq p_{\rm{T}}\leq 20$
GeV/$c$, depending on the particle species. The $p_{\rm{T}}$ spectra,
integrated yields, and particle yield ratios are discussed as a function of
collision energy and compared with measurements at lower $\sqrt{s}$ and with
results from various general-purpose QCD-inspired Monte Carlo models. A
hardening of the spectra at high $p_{\rm{T}}$ with increasing collision energy
is observed, which is similar for all particle species under study. The
transverse mass and $x_{\rm{T}}\equiv 2p_{\rm{T}}/\sqrt{s}$ scaling properties
of hadron production are also studied. As the collision energy increases from
$\sqrt{s}$ = 7 to 13 TeV, the yields of non- and single-strange hadrons
normalized to the pion yields remain approximately constant as a function of
$\sqrt{s}$, while ratios for multi-strange hadrons indicate enhancements. The
$p_{\rm{T}}$-differential cross sections of $\pi^{\pm}$, $\rm{K}^{\pm}$ and
$\rm{p}$ ($\overline{\rm{p}}$) are compared with next-to-leading order
perturbative QCD calculations, which are found to overestimate the cross
sections for $\pi^{\pm}$ and $\rm{p}$ ($\overline{\rm{p}}$) at high
$p_{\rm{T}}$.

##  1 Introduction

Identified particle spectra and yields, which are among the most fundamental
physical observables in high-energy hadronic collisions, have been intensively
studied in hadron-collider and cosmic-ray physics for many decades [1]. Hadron
production at collider energies originates from soft and hard scattering
processes at the partonic level. Hard scatterings, where two partons interact
with a large momentum transfer, are responsible for the production of
particles with high transverse momenta. This process is theoretically
described by perturbative Quantum Chromodynamics (pQCD) calculations based on
the factorization theorem [2]. In this approach, the cross section is a
convolution of the parton distribution function (PDF), the partonic QCD matrix
elements, and the fragmentation function (FF). The PDFs describe the
probability densities of finding a parton with a specific flavor carrying
fraction $x$ of the proton momentum, whereas the FFs encode the probability
densities that the parton with a specific flavor fragments into a hadron
carrying a fraction of the parton’s longitudinal momentum; both considered at
a given energy scale. At the LHC, with increasing center-of-mass collision
energy ($\sqrt{s}$), the lower $x$ regime is probed and contributions from
hard-scattering processes increase. In the kinematic region probed by these
measurements, high-$p_{\mathrm{T}}$ particles dominantly originate from the
fragmentation of gluons [3, 4]. Parameterizations of both the PDFs and FFs are
derived from global analyses [5, 6] based on fits to experimental data at
various $\sqrt{s}$ with next-to-leading order (NLO) accuracy. These include
single-inclusive hadron production in semi-inclusive electron-positron
annihilation data, semi-inclusive deep-inelastic scattering, and single
inclusive hadron spectra at high $p_{\mathrm{T}}$, notably including results
at LHC energies. Results presented in this paper can be used as further input
for these studies. In particular, identified particle spectra provide new
constraints on the gluon-to-pion and, especially, gluon-to-kaon fragmentation
functions [7, 8, 9, 10]. While particle production at high $p_{\mathrm{T}}$ is
expected to be calculable with pQCD, the LHC results are in general not well
reproduced by pQCD calculations, see Ref. [6] and references therein. Charged
particle production at high $p_{\mathrm{T}}$ is known to scale with
$x_{\mathrm{T}}\equiv 2p_{\mathrm{T}}/\sqrt{s}$, as observed in a wide energy
range up to $\sqrt{s}=7\text{\,TeV}$. This has been observed by the CDF
Collaboration in p$\bar{\text{p}}$ collisions at the Tevatron [11, 12], by the
UA1 Collaboration at the CERN SPS [13], by the STAR Collaboration in pp
collisions at RHIC [14], and by the CMS Collaboration [15] at the CERN LHC.
Above $x_{\mathrm{T}}\simeq 10^{-2}$, significant deviations from the leading-
twist NLO pQCD predictions have been reported in Ref. [16] and are
investigated in this paper.

The bulk of particles produced at low transverse momenta
($p_{\mathrm{T}}<2\text{\,GeV/$c$}$) originate from soft scattering processes
involving small momentum transfers. In this regime, particle production cannot
be calculated from first principles. Instead, calculations rely on QCD
inspired phenomenological models, which are tuned to reproduce previous
measurements. Hence, measurements at low $p_{\mathrm{T}}$ provide further
important constraints on such models. The universal transverse mass,
$m_{\mathrm{T}}\equiv\sqrt{m^{2}+p_{\mathrm{T}}^{2}}$, scaling originally
proposed by R. Hagedorn [17] was first seen to hold approximately at ISR
energies [18]. It was then observed by the PHENIX [19, 20] and STAR [21]
collaborations to hold only separately for mesons and baryons at RHIC
energies, by applying the approximate $m_{\mathrm{T}}$ scaling relation
respectively for pions and protons. At $\sqrt{s}=900\text{\,GeV}$ a
disagreement was observed for charged kaons and $\phi(1020)$ mesons, which
indicated a breaking of the generalized scaling behavior [22]. Moreover,
recent studies, based on identified particle spectra measured in pp collisions
at $\sqrt{s}=7\text{\,TeV}$ by ALICE, indicate that $m_{\mathrm{T}}$ scaling
also breaks in the low-$p_{\mathrm{T}}$ region [23]. These observations
motivate studies of the applicability of $m_{\mathrm{T}}$ scaling of particle
production through the precise measurement of identified hadrons at
$\sqrt{s}=13\text{\,TeV}$.

The results reported in this paper are the measurements of the production of
$\pi^{\pm}$, $\mathrm{K}^{\pm}$, $\mathrm{K_{S}^{0}}$,
$\mathrm{K^{*}\kern-0.89996pt(892)^{0}}$,
$\mathrm{\overline{K}^{*}(892)^{0}}$, ${\rm p}$, ${\rm\overline{p}}$,
$\phi(1020)$, $\Lambda$, $\overline{\Lambda}$, $\Xi^{-}$,
$\overline{\Xi}^{+}$, $\Omega^{-}$ and $\overline{\Omega}^{+}$ at the highest
collision energies, and therefore extend the studies of the energy dependence
of the production of light-flavor hadrons into new territory. The study of the
production of the $\mathrm{K^{*}\kern-0.89996pt(892)^{0}}$ and $\phi(1020)$
resonances, containing respectively one and two strange valence quarks,
contributes to a better understanding of strange particle production
mechanisms. Because of their short lifetimes (${\sim}4\text{\,fm/$c$}$ for
$\mathrm{K^{*}\kern-0.89996pt(892)^{0}}$ and ${\sim}46\text{\,fm/$c$}$ for
$\phi(1020)$), their decay daughters may undergo re-scattering and/or
regeneration processes that affect their yields and the shapes of their
$p_{\mathrm{T}}$ distributions. In addition, multi-strange baryons,
$\Omega^{-}$ ($\overline{\Omega}^{+}$) and $\Xi^{-}$ ($\overline{\Xi}^{+}$),
are of crucial importance due to their dominant strange (s) quark content.
Furthermore, the production of $\mathrm{K_{S}^{0}}$, $\Lambda$, and
$\overline{\Lambda}$ in pp collisions at $\sqrt{s}=7\text{\,TeV}$ is reported
here for the first time, completing the set of reference measurements at that
energy [24, 25, 26, 27].

The present measurements serve as important baselines for studies of particle
production as a function of the charged-particle multiplicity [28] or event
shape (e.g. spherocity) [29] and also provide input to tune the modeling of
several contributions in Monte Carlo (MC) event generators such as PYTHIA [30,
31] and EPOS-LHC [32]. In addition, measurements in minimum bias pp collisions
reported in this paper serve as reference data to study nuclear effects in
proton-lead (p–Pb) and lead-lead (Pb–Pb) collisions.

The paper is organized as follows. In Section 2 the ALICE experimental
apparatus and the analyzed data samples are described, focusing on the
detectors which are relevant for the presented measurements. In Section 3 the
details of the event and track selection criteria and of the Particle
IDentification (PID) techniques are discussed. The results are given in
Section 4, in which the $p_{\mathrm{T}}$ spectra and the extraction procedures
for the $p_{\mathrm{T}}$-integrated yield and average $p_{\mathrm{T}}$ are
presented. Section 5 discusses the results, followed by a summary in Section
6. For the remainder of this paper, the masses will be omitted from the
symbols of the strongly decaying particles, which will be denoted as
$\mathrm{K^{*0}}$, $\mathrm{\overline{K}^{*0}}$, and $\phi$.

##  2 Experimental setup

A detailed description of the ALICE detector and its performance can be found
in Refs. [33, 34]. The main subsystems of the ALICE detector used in this
analysis are the V0 detector, the Inner Tracking System (ITS), the Time
Projection Chamber (TPC), the Time of Flight (TOF) detector, and the High-
Momentum Particle Identification Detector (HMPID).

The V0 detector [35] is used for triggering and beam background suppression.
It is made up of two scintillator arrays placed along the beam axis on each
side of the interaction point (IP) at $z=340$ cm and $z=-90$ cm, covering the
pseudorapidity regions $2.8<\eta<5.1$ (V0A) and $-3.7<\eta<-1.7$ (V0C),
respectively.

In the measurements of light-flavor hadrons, primary charged particles are
considered. Primary particles are defined as particles with a mean proper
lifetime $\tau$ that is larger than 1 cm/$c$, which are either produced
directly in the interaction or from decays of particles produced at the
interaction vertex with $\tau$ shorter than 1 cm/$c$. This excludes particles
produced in interactions with the detector material [36]. Primary charged-
hadron tracks are reconstructed by the ITS and TPC detectors, which have full
azimuthal acceptance within $|\eta|<0.8$ for full-length tracks. They are
located inside a solenoidal magnet providing a magnetic field of $B=0.5$ T.

The ITS [33, 37] is a silicon tracking detector made up of six concentric
cylindrically-shaped layers, measuring high-resolution space points near the
collision vertex. The two innermost layers consist of Silicon Pixel Detectors
(SPD) used to reconstruct the primary vertex of the collision and short track
segments called “tracklets". The four outer layers are equipped with silicon
drift (SDD) and strip (SSD) detectors and allow measurement of the specific
energy loss ($\text{d}E/\text{d}x$) with a relative resolution of about 10%.
The ITS is also used as a stand-alone tracking detector to reconstruct charged
particles with momenta below $200\text{\,MeV/$c$}$ that are deflected or decay
before reaching the TPC.

The TPC [38] is the main tracking detector of ALICE. It is a large volume
cylindrical drift detector spanning the approximate radial and longitudinal
ranges $85<r<250$ cm and $-250<z<250$ cm, respectively. The endcaps of the TPC
are equipped with multiwire proportional chambers (MWPCs) segmented radially
into pad rows. Together with the measurement of the drift time, the TPC
provides three dimensional space point information, with up to 159 tracking
points. Charged tracks originating from the primary vertex can be
reconstructed down to $p{\sim}100\text{\,MeV/$c$}$ [34], albeit with a lower
tracking efficiency for identified charged hadrons with
$p_{\mathrm{T}}<200\text{\,MeV/$c$}$. Combining information from the ITS and
TPC allows the momenta of charged particles to be measured for momenta from
0.05 to $100\text{\,GeV/$c$}$ with a resolution of 1–10%, depending on
$p_{\mathrm{T}}$. The TPC provides charged-hadron identification via
measurement of the specific energy loss $\text{d}E/\text{d}x$ in the fill gas,
with a resolution of ${\sim}5\%$ [38].

The Time of Flight detector (TOF) [39, 40, 41] is a cylindrical array of
multi-gap resistive plate chambers which sits outside the TPC. It covers the
pseudorapidity range $|\eta|<0.9$ with (almost) full azimuthal acceptance. The
total time-of-flight resolution, including the resolution on the collision
time, is about 90 ps in pp collisions.

The HMPID consists of seven proximity focusing Ring Imaging Cherenkov (RICH)
counters. Primary charged particles penetrate the radiator volume, filled with
liquid $\text{C}_{6}\text{F}_{14}$, and generate Cherenkov photons that are
converted into photoelectrons in thin CsI-coated photocathodes. Photo-electron
clusters, together with pad clusters (also called “MIP” clusters) associated
with the primary ionization of a particle, form Cherenkov rings. The amplified
signal is read out by MWPCs, filled with $\text{CH}_{4}$. The detector covers
$|\eta|<0.5$ and $1.2^{\circ}<\varphi<58.5^{\circ}$, which corresponds to
${\sim}5\%$ of the TPC geometrical acceptance.

##  3 Event and track selection

###  3.1 Event selection

The measurements at $\sqrt{s}=13\text{\,TeV}$ are obtained from a minimum bias
data sample of pp collisions collected in June 2015 during a period of low
pileup in LHC Run 2. The minimum bias trigger required at least one hit in
both of the V0 scintillator arrays in coincidence with the arrival of proton
bunches from both directions along the beam axis. The mean number of inelastic
proton-proton interactions per single bunch crossing ranges between 2% and
14%. A requirement of a coincidence of signals in both V0A and V0C detectors
removes contamination from single-diffractive and electromagnetic events.
Contamination arising from beam-induced background events, produced outside
the interaction region, is removed offline by using timing information from
the V0 detector, which has a time resolution better than 1 ns. Background
events are further rejected by exploiting the correlation between the number
of clusters and the multiplicity of tracklets in the SPD. From the triggered
events, only events with a reconstructed primary vertex are considered for the
analyses. Additionally, the position of the primary vertex along the beam axis
is required to be within $\pm 10$ cm with respect to the nominal interaction
point (center of the ALICE barrel). This requirement ensures that the vast
majority of reconstructed tracks are within the central barrel acceptance
($|\eta|<0.8$) and it reduces background events by removing unwanted
collisions from satellite bunches. Contamination from pileup events, which
have more than one pp collision per bunch crossing, were rejected offline by
excluding events with multiple primary vertices reconstructed in the SPD [34].
The pileup-rejected events are less than 1% of the total sample of minimum
bias events. The size of the analyzed sample after selections ranges between
40 and 60 million events (corresponding to an integrated luminosity
$0.74-1.1\,\rm{nb}^{-1}$), depending on the requirements of the analyses of
the different particle species.

The measurements of $\mathrm{K_{S}^{0}}$, $\Lambda$, and $\overline{\Lambda}$
at $\sqrt{s}=7\text{\,TeV}$ are obtained by analyzing a sample of about 150
million events (corresponding to an integrated luminosity
$2.41\,\rm{nb}^{-1}$) collected in 2010 during the LHC Run 1 data taking
period. The corresponding trigger and event selection criteria applied were
very similar to those used for the measurements at $\sqrt{s}=13\text{\,TeV}$;
see Refs. [24, 42] for details on the triggering and event selection for these
periods.

All corrections are calculated using Monte Carlo events from PYTHIA 6 and
PYTHIA 8. PYTHIA 6.425 (Perugia 2011 tune) and PYTHIA 8.210 (Monash 2013 tune)
event generators were used for $\sqrt{s}=13\text{\,TeV}$. The PYTHIA 6.421
(Perugia 0 tune) was used for $\mathrm{K_{S}^{0}}$ and $\Lambda$ at
$\sqrt{s}=7\text{\,TeV}$ because that production was used for correcting the
other $\sqrt{s}=7\text{\,TeV}$ analyses. The particles produced using these
event generators were propagated through a simulation of the ALICE detector
using GEANT3 [43].

###  3.2 Track selection

Tracks from charged particles are reconstructed in the TPC and ITS detectors
and then propagated to the outer detectors and matched with reconstructed
points in the TOF and HMPID. Additionally, in the analysis of $\pi^{\pm}$,
$\mathrm{K}^{\pm}$ and ${\rm p}({\rm\overline{p}})$, a dedicated tracking
algorithm based only on ITS information (ITS stand-alone, ITS-sa) was used to
reconstruct low momentum tracks. In the measurements, global tracks, which are
reconstructed using the combined ITS and TPC information, are distinguished
from ITS-sa tracks.

For analyses using global tracks, track selection criteria are applied to
limit the contamination due to secondary particles, to maximize tracking
efficiency and improve the $\text{d}E/\text{d}x$ and momentum resolution for
primary charged particles, and to guarantee an optimal PID quality. The number
of crossed pad rows in the TPC is required to be at least 70 (out of a maximum
possible of 159); the ratio of the number of crossed pad rows to the number of
findable clusters (that is the number of geometrically possible clusters which
can be assigned to a track) is restricted to be greater than 0.8, see Ref.
[34] for the details. The goodness-of-fit values $\chi^{2}$ per cluster
($\chi^{2}/N_{\rm clusters}$) of the track fit in the TPC must be less than 4.
Tracks must be associated with at least one cluster in the SPD and the
$\chi^{2}$ values per cluster in the ITS are restricted in order to select
high-quality tracks. The distance of closest approach (DCA) to the primary
vertex in the plane perpendicular to the beam axis (DCAxy) is required to be
less than 7 times the resolution of this quantity; this selection is
$p_{\mathrm{T}}$ dependent, i.e.
$\text{DCA}_{xy}<7\times(0.0015+0.05\times(p_{\mathrm{T}}/(\text{\,GeV/$c$}))^{-1.01})$
cm. A loose selection criterion is also applied on the the DCA in the beam
direction (DCAz), by rejecting tracks with DCAz larger than 2 cm, to remove
tracks from possible residual pileup events. The transverse momentum of each
track must be greater than 150 $\mathrm{MeV}\kern-1.79993pt/\kern-0.89996ptc$
and the pseudorapidity is restricted to the range $|\eta|<0.8$ to avoid edge
effects in the TPC acceptance. Additionally, tracks produced by the
reconstructed weak decays of pions and kaons (the “kink" decay topology) are
rejected.

For the topological reconstruction of weakly decaying particles, the selected
global tracks are combined using specific algorithms, as described in Sec.
4.3. Track selection criteria are the same applied for global tracks with a
few exceptions: for tracks used in the reconstruction of $\mathrm{K_{S}^{0}}$,
$\Lambda$, $\overline{\Lambda}$, $\Xi^{-}$, $\overline{\Xi}^{+}$,
$\Omega^{-}$, and $\overline{\Omega}^{+}$, no ITS information is required and
special selection criteria are applied on the DCA to the collision vertex, as
shown in Table 3. The kink topology tracks that are used to reconstruct the
weak decays of $\mathrm{K}^{\pm}$ do not have ITS information. For the latter,
removal of contributions from pileup collisions outside the trigger proton
bunch (“out-of-bunch pileup”) is achieved by requiring that at least one
charged decay track matches a hit in a “fast” detector (either the ITS or the
TOF detector).

ITS stand-alone tracking uses similar selection criteria to those mentioned
above. Tracks are required to have at least four ITS clusters, with at least
one in the SPD, three in the SSD and SDD and $\chi^{2}/N_{\rm clusters}<2.5$.
This further reduces contamination from secondary tracks and provides high
resolution for the track impact parameter and optimal resolution for
$\text{d}E/\text{d}x$. Similar to global tracks, a $p_{\mathrm{T}}$-dependent
parameterization of the DCAxy selection is used, but with different parameters
to account for the different resolution. For the $p_{\mathrm{T}}$ ranges used
in this analysis, the selected ITS-sa tracks have the same $p_{\mathrm{T}}$
resolution as those measured in pp collisions at $\sqrt{s}=7\text{\,TeV}$: 6%
for pions, 8% for kaons, and 10% for protons [24].

##  4 Data analysis techniques

Table 1 lists the basic characteristics of the particles studied in this
paper. This section describes the techniques used to measure the yields of the
various hadron species. In Sec. 4.1, aspects common to all analyses are
described, including the correction and normalization procedure and the common
sources of systematic uncertainties. Next, the analysis of each hadron species
is described in detail. The measurements of charged pions, charged kaons, and
(anti)protons, which are performed using several different PID techniques, are
described in Sec. 4.2. It is worth noting that charged kaons are also
identified using the kink topology of their two-body decays. The measurements
of weakly decaying strange hadrons ($\mathrm{K_{S}^{0}}$, $\Lambda$,
$\Xi^{-}$, $\Omega^{-}$ and their antiparticles) are reported in Sec. 4.3,
followed by the strongly decaying resonances ($\mathrm{K^{*0}}$,
$\mathrm{\overline{K}^{*0}}$, and $\phi$) in Sec. 4.4.

Table 1: Main characteristics of the reconstructed particles: valence quark
content, mass, proper decay length ($c\tau$), the decay channel studied in
this paper, and the corresponding branching ratio (B.R.) [44].

| Particle |  | Valence |  | Mass |  | $c\tau$ |  | Decay |  | B.R. (%) |
---|---|---|---|---|---|---|---|---|---|---|---|---
|  | Quark Content |  | (MeV/$c^{2}$) |  |  |  |
| $\pi^{+}$ |  | $\mathrm{u\overline{d}}$ |  | 139.57 |  | 7.8 m |  | — |  | — |
| $\mathrm{K}^{+}$ |  | $\mathrm{u\overline{s}}$ |  | 493.68 |  | 3.7 m |  | $\mathrm{K}^{+}\rightarrow\mu^{+}\nu_{\mu}$ |  | $63.56\pm 0.11$ |
| $\mathrm{K_{S}^{0}}$ |  | $\frac{1}{\sqrt{2}}(\mathrm{d\overline{s}}-\mathrm{\overline{d}s})$ |  | 497.61 |  | 2.68 cm |  | $\mathrm{K_{S}^{0}}\rightarrow\pi^{+}+\pi^{-}$ |  | $69.20\pm 0.05$ |
Mesons | $\mathrm{K^{*0}}$ |  | $\mathrm{d\overline{s}}$ |  | 895.55 |  | 4.16 fm |  | $\mathrm{K^{*0}}\rightarrow\pi^{-}+\mathrm{K}^{+}$ |  | $\sim 66.6$ |
| $\phi$ |  | $\mathrm{s\overline{s}}$ |  | 1019.46 |  | 46.2 fm |  | $\phi\rightarrow\mathrm{K}^{+}+\mathrm{K}^{-}$ |  | $49.2\pm 0.5$ |
| ${\rm p}$ |  | uud |  | 938.27 |  | — |  | — |  | — |
| $\Lambda$ |  | uds |  | 1115.68 |  | 7.89 cm |  | $\Lambda\rightarrow{\rm p}+\pi^{-}$ |  | $63.9\pm 0.5$ |
| $\Xi^{-}$ |  | dss |  | 1321.71 |  | 4.91 cm |  | $\Xi^{-}\rightarrow\Lambda+\pi^{-}$ |  | $99.887\pm 0.035$ |
Baryons | $\Omega^{-}$ |  | sss |  | 1672.45 |  | 2.46 cm |  | $\Omega^{-}\rightarrow\Lambda+\mathrm{K}^{-}$ |  | $67.8\pm 0.7$ |

###  4.1 Common aspects of all analyses

In several of the analyses presented below, the measured PID signal is
compared to the expected value based on various particle mass hypotheses. The
difference between the measured and expected values is expressed in terms of
$\sigma$, the standard deviation of the corresponding measured signal
distribution. The size of this difference, in multiples of $\sigma$, is
denoted $n_{\sigma}$. In the following, the $\sigma$ values accounting for the
resolution of the PID signals measured in the TPC and TOF detectors are
denoted as $\sigma_{\mathrm{TPC}}$ and $\sigma_{\mathrm{TOF}}$, respectively.

The corrected yield of each hadron species as a function of $p_{\mathrm{T}}$
is

$Y_{\mathrm{corr}}=\frac{Y_{\mathrm{raw}}}{\Delta
p_{\mathrm{T}}\,\Delta{y}}\times\frac{f_{\mathrm{SL}}}{A\times\varepsilon}\times(1-f_{\mathrm{cont}})\times
f_{\mathrm{cross.sec.}}~{}.$ (1)

$Y_{\mathrm{corr}}$ is obtained by following the procedure described in
previous publications. Here, $Y_{\mathrm{raw}}$ is the number of particles
measured in each $p_{\mathrm{T}}$ bin and $A\times\varepsilon$ is the product
of the acceptance and the efficiency (including PID efficiency, matching
efficiency, detector acceptance, reconstruction, and selection efficiencies).
Monte Carlo simulations are used to evaluate $A\times\varepsilon$, which takes
on similar values to those found in our previous analyses. The factor
$f_{\mathrm{SL}}$, also known as the “signal-loss” correction, accounts for
reductions in the measured particle yields due to event triggering and primary
vertex reconstruction. Such losses are more important at low $p_{\mathrm{T}}$,
since events that fail the trigger conditions or fail to have a
reconstructible primary vertex tend to have softer particle $p_{\mathrm{T}}$
spectra than the average inelastic collision. For $\sqrt{s}=13\text{\,TeV}$,
$f_{\mathrm{SL}}$ deviates from unity by a few percent at low $p_{\mathrm{T}}$
to less than one percent for $p_{\mathrm{T}}\gtrsim 2\text{\,GeV/$c$}$. The
trigger configuration used for $\sqrt{s}=7\text{\,TeV}$ resulted in negligible
signal loss, thus $f_{\mathrm{SL}}$ is set to unity for this energy. The
factor $(1-f_{\mathrm{cont}})$ is used to correct for contamination from
secondary and misidentified particles; $f_{\mathrm{cont}}$ is non-zero only
for the measurements of $\pi^{\pm}$, $\mathrm{K}^{\pm}$, ${\rm
p}({\rm\overline{p}})$, $\Lambda$, and $\overline{\Lambda}$, and it is more
important at low $p_{\mathrm{T}}$. The computation of $f_{\mathrm{cont}}$ for
those species is described further in the relevant sections below. The factor
$f_{\mathrm{cross.sec.}}$ corrects for inaccuracies in the hadronic production
cross sections in GEANT3, which is used in the calculation of
$A\times\varepsilon$ to describe the interactions of hadrons with the detector
material of ALICE. GEANT4 and FLUKA [45], which have more accurate
descriptions of the hadronic cross sections, are used to calculate the
correction factor, which can be different from unity by up to a few percent.
The correction $f_{\mathrm{cross.sec.}}$ is applied only for the analyses of
$\mathrm{K}^{-}$, ${\rm\overline{p}}$, $\overline{\Lambda}$,
$\overline{\Xi}^{+}$, and $\overline{\Omega}^{+}$.

After correction, the yields are normalized to the number of inelastic pp
collisions using the ratio of the ALICE visible cross section to the total
inelastic cross section. This ratio is $0.852^{+0.062}_{-0.030}$ for
$\sqrt{s}=7\text{\,TeV}$ [46] and $0.7448\pm 0.0190$ for
$\sqrt{s}=13\text{\,TeV}$ [47, 48].

The procedures for the estimation of systematic uncertainties strictly follow
those applied in our measurements from LHC Run 1. All described uncertainties
are assumed to be strongly correlated among adjacent $p_{\mathrm{T}}$ bins.
For the evaluation of the total systematic uncertainty in every analysis, all
contributions originating from different sources are considered to be
uncorrelated and summed in quadrature. Components of uncertainties related to
the ITS-TPC matching efficiency correction and to the event selection are
considered correlated among different measurements. The systematic uncertainty
due to the normalization to the number of inelastic collisions is
$\pm~{}2.6\%$ for $\sqrt{s}=13\text{\,TeV}$ and ${}^{+7.3\%}_{-3.5\%}$ for
$\sqrt{s}=7\text{\,TeV}$ independent of $p_{\mathrm{T}}$. This uncertainty is
common to all measured $p_{\mathrm{T}}$ spectra and $\text{d}N/\text{d}y$
values (see Sec. 5.1) at a given energy. The systematic uncertainty associated
to possible residual contamination from pileup events was estimated varying
pileup rejection criteria and was found to be of 1%. The singal loss
correction has a small dependence on the Monte Carlo event generator used to
calculate it. These variations result in $p_{\mathrm{T}}$-dependent
uncertainties that are largest at low $p_{\mathrm{T}}$, where they have values
of $0.2\%$ for $\Omega$, ${\sim}1\%$ for $\pi^{\pm}$, $\mathrm{K}^{\pm}$,
${\rm p}({\rm\overline{p}})$, and $\Xi$, and ${\sim}2\%$ for
$\mathrm{K_{S}^{0}}$, $\overline{\Lambda}$, $\mathrm{K^{*0}}$, and $\phi$.

The systematic uncertainty accounting for the limited knowledge of the
material budget is estimated by varying the amount of detector material in the
MC simulations within its expected uncertainties [34]. For the analysis of
$\pi^{\pm}$, $\mathrm{K}^{\pm}$, ${\rm p}({\rm\overline{p}})$,
$\mathrm{K^{*0}}$, and $\phi$, the values are taken from the studies reported
in Refs. [49] and [50]. This uncertainty is estimated to be around $3.3\%$ for
$\mathrm{K}^{\pm}$, $1.1\%$ for $\pi^{\pm}$, $1.8\%$ for ${\rm
p}({\rm\overline{p}})$, $3\%$ for $\mathrm{K^{*0}}$, and $2\%$ for $\phi$; it
is largest at low momenta and tends to be negligible towards higher momenta.
For the measurement of $\mathrm{K_{S}^{0}}$ and $\Lambda$ at
$\sqrt{s}=7\text{\,TeV}$, the material budget uncertainty is estimated to be
$4\%$, independent of $p_{\mathrm{T}}$. For the measurements of
$\mathrm{K_{S}^{0}}$, $\Lambda$, $\Xi$ and $\Omega$ at
$\sqrt{s}=13\text{\,TeV}$, the material budget uncertainty is $p_{\mathrm{T}}$
dependent for low $p_{\mathrm{T}}$ ($\lesssim 2\text{\,GeV/$c$}$) and constant
at higher $p_{\mathrm{T}}$. For low $p_{\mathrm{T}}$, the uncertainty reaches
maximum values of about 4.7% for $\mathrm{K_{S}^{0}}$, 6.7% for $\Lambda$, 6%
for $\Xi$, and 3.5% for $\Omega$; at high $p_{\mathrm{T}}$, the uncertainty is
less than 1% for $\mathrm{K_{S}^{0}}$, $\Lambda$, and $\Xi$, and about 1.5%
for $\Omega$.

The systematic uncertainty due to the limited description of the hadronic
interaction cross sections in the transport code is evaluated using GEANT4 and
FLUKA. This leads to uncertainties of up to $2.8\%$ for $\pi^{\pm}$, $2.5\%$
for $\mathrm{K}^{\pm}$, $0.8\%$ for ${\rm p}$, and $5\%$ for
${\rm\overline{p}}$ [49]. It is at most $3\%$ for $\mathrm{K^{*0}}$, $2\%$ for
$\phi$ and $1-2\%$ for the strange baryons. It is negligible for
$\mathrm{K_{S}^{0}}$ at both reported collision energies. In the following
sections, details are given on the contributions (specific to each analysis)
related to track or topological selections and signal extraction methods, as
well as those related to feed-down.

###  4.2 Identification of primary charged pions, charged kaons, and
(anti)protons

To measure the production of primary charged pions, kaons, and (anti)protons
over a wide range of $p_{\mathrm{T}}$, five analyses using distinct PID
techniques were carried out. The individual analyses follow the techniques
adopted in previous measurements based on data collected at lower center-of-
mass energies and for different collision systems during LHC Run 1 [51, 52,
24, 25, 53]. The $p_{\mathrm{T}}$ spectra have been measured from
$p_{\mathrm{T}}=0.1\text{\,GeV/$c$}$ for pions,
$p_{\mathrm{T}}=0.2\text{\,GeV/$c$}$ for kaons, and
$p_{\mathrm{T}}=0.3\text{\,GeV/$c$}$ for protons, up to $20\text{\,GeV/$c$}$
for all three species. The individual analyses with their respective
$p_{\mathrm{T}}$ reaches are summarized in Tab. 2. All the analysis techniques
are extensively described in Refs. [51, 49, 24, 7]. Each procedure is
discussed separately in Secs. 4.2.1–4.2.5, with special emphasis on those
aspects that are relevant for the current measurements. The results for the
different analyses are then combined as described in Sec. 4.2.6.

The calculation of $f_{\mathrm{cont}}$ in Eq. 1 at low $p_{\mathrm{T}}$ is
performed by subtracting the secondary $\pi^{\pm}$, $\mathrm{K}^{\pm}$, and
${\rm p}({\rm\overline{p}})$ from the primary particle sample. This method is
data-driven and it is based on the measured distance of closest approach to
the primary vertex in the plane transverse to the beam direction
($\text{DCA}_{xy}$), following the same procedure adopted in Ref. [24]. The
$\text{DCA}_{xy}$ distribution of the selected tracks was fitted in every
$p_{\mathrm{T}}$ bin with Monte Carlo templates composed of three ingredients:
primary particles, secondaries from material and secondaries from weak decays,
each accounting for the expected shapes of the distribution. Because of the
different track and PID selection criteria, the contributions are different
for each analysis. The resulting corrections are significant at low
$p_{\mathrm{T}}$ and decrease towards higher $p_{\mathrm{T}}$ due to decay
kinematics. Up to $p_{\mathrm{T}}=2\text{\,GeV/$c$}$, the contamination is
$2-10\%$ for pions, up to $20\%$ for kaons (in the narrow momentum range where
the $\text{d}E/\text{d}x$ response for kaons and secondary electrons overlap),
and $15-20\%$ for protons.

Table 2: Summary of the kinematic ranges ($p_{\rm T}$ (GeV/$c$) and $\eta$ or
$y$) covered by the individual analyses for the measurement of $\pi^{\pm}$,
$\text{K}^{\pm}$ and $(\bar{\text{p}})\text{p}$ in pp collisions at
$\sqrt{s}=13\text{\,TeV}$.

Analysis technique |  |  $p_{\mathrm{T}}$ range (GeV/$c$) | $\eta$ or $y$ range
---|---|---|---
|  | $\pi^{+}+\pi^{-}$ |  | ${\rm K}^{+}+{\rm K}^{-}$ |  | ${\rm p}+\bar{\rm p}$ |  |
ITS-sa |  | $0.1-0.7$ |  | $0.2-0.6$ |  | $0.3-0.65$ |  | $|y|<0.5$
TPC-TOF fits |  | $0.3-3.0$ |  | $0.3-3.0$ |  | $0.4-3.0$ |  | $|y|<0.5$
|  |  |  | $|\eta|<0.4$
HMPID |  | $1.5-4$ |  | $1.5-4$ |  | $1.5-6$ |  | $|y|<0.5$
Kinks |  | $-$ |  | $0.2-7$ |  | $-$ |  | $|y|<0.5$
TPC rel. rise |  | $2-20$ |  | $3-20$ |  | $3-20$ |  | $|\eta|<0.8$

The main sources of systematic uncertainties for each analysis are summarized
in Tab. LABEL:tab:syst:pikp, including contributions common to all analyses.
The systematic uncertainty due to the subtraction of secondary particles is
estimated by changing the fit range of the $\text{DCA}_{xy}$ distribution,
resulting in uncertainties of up to $4\%$ for protons and $1\%$ for pions,
with negligible uncertainties for kaons. The uncertainty due to the matching
of TPC tracks with ITS hits is estimated to be in the range ${\sim}1-5\%$ for
$p_{\mathrm{T}}\lesssim 3\text{\,GeV/$c$}$ depending on $p_{\mathrm{T}}$,
while it takes values around $6\%$ at higher $p_{\mathrm{T}}$. This
uncertainty together with that resulting from the variation of the track
quality selection criteria lead to the systematic uncertainty of the global
tracking efficiency that varies from 2.2 to 7.3% from low to high
$p_{\mathrm{T}}$, independent of particle species.

####  4.2.1 ITS stand-alone

In the “ITS stand-alone” analysis, both tracking and PID are performed based
on information from the ITS detector only. For the present data sample, the
contribution of tracks with wrongly assigned clusters in the ITS is negligible
due to the low pseudorapidity density of charged particles,
$\langle\text{d}N_{\text{ch}}/\text{d}\eta\rangle_{|\eta|<0.5}=5.31\pm 0.18$,
measured in the pseudorapidity region $|\eta|<0.5$ [54]. The average
$\text{d}E/\text{d}x$ signal in the four outer ITS layers used for PID is
estimated by means of the truncated mean method [51]. The measured
$\text{d}E/\text{d}x$ for the sample of ITS-sa tracks is shown in the top left
panel of Fig. LABEL:fig:PidPerformance, along with the Bethe–Bloch
parametrization of the most probable values, which is the same as the one used
for the LHC Run 1 analyses [24]. Two identification strategies were used. In
the main analysis, a unique identity is assigned to the ITS-sa track according
to the mass hypothesis for which the expected specific energy-loss value is
the closest to the measured $\text{d}E/\text{d}x$ for a track with momentum
$p$. The second analysis strategy uses the Bayesian PID approach [55], based
on likelihood parametrization with a set of iterative prior probabilities. The
identification is based on the maximal probability method in which the species
with the highest probability is assigned to a track. For
$p_{\mathrm{T}}<160\text{\,MeV/$c$}$, where the e$/\pi$ separation power in
the ITS allows high-purity identification of electrons, four mass hypotheses
(e/$\pi$/K/p) are considered. For $p_{\mathrm{T}}>160\text{\,MeV/$c$}$,
electrons and pions cannot be separated using their $\text{d}E/\text{d}x$ in
the ITS detector, and the Bayesian approach is based on the $\pi$/K/p mass
hypotheses only.

For the ITS-sa analysis, the systematic uncertainties related to the PID
procedure originate from the different techniques that are used (the truncated
mean method and the Bayesian PID approach). These range from about 1 to 5%
depending on particle species and $p_{\mathrm{T}}$. The Lorentz force causes
the migration of the cluster position in the ITS by driving the charge in
opposite directions depending on the polarity of the magnetic field of the
experiment ($\mathbf{E}\times\mathbf{B}$ effect). The uncertainty related to
this effect is estimated by analyzing data samples with opposite magnetic
field polarities, for which a difference at the level of ${\sim}2\%$ is
observed.

####  4.2.2 TPC-TOF fits

In the so-called “TPC-TOF fits” analysis, the distributions of the specific
energy loss $\text{d}E/\text{d}x$ measured in the TPC and the velocity $\beta$
measured in the TOF detector are fitted with functions that describe the PID
signals for different track momentum ($p$) intervals. The TPC provides a
$3\sigma_{\mathrm{TPC}}$ separation between pions and kaons up to
$p_{\mathrm{T}}{\sim}600\text{\,MeV/$c$}$ and between kaons and protons up to
$p_{\mathrm{T}}{\sim}800\text{\,MeV/$c$}$.

Particle identification using this technique is possible in the
$p_{\mathrm{T}}$ ranges $0.3<p_{\mathrm{T}}<0.5\text{\,GeV/$c$}$,
$0.3<p_{\mathrm{T}}<0.6\text{\,GeV/$c$}$, and
$0.4<p_{\mathrm{T}}<0.8\text{\,GeV/$c$}$ for $\pi^{\pm}$, $\mathrm{K}^{\pm}$,
and ${\rm p}({\rm\overline{p}})$, respectively. The extraction of the raw
yield for a given species is done by integrating the
$\mathrm{d}^{2}N/\mathrm{d}p_{\mathrm{T}}\mathrm{d}n_{\sigma}$ distribution in
these $p_{\mathrm{T}}$ intervals. In addition, in the transverse momentum
ranges $p_{\mathrm{T}}<0.4\text{\,GeV/$c$}$ for $\pi^{\pm}$,
$p_{\mathrm{T}}>0.45\text{\,GeV/$c$}$ for $\mathrm{K}^{\pm}$, and
$p_{\mathrm{T}}>0.6\text{\,GeV/$c$}$ for ${\rm p}({\rm\overline{p}})$, a
Gaussian fit is used to remove the background contribution ($\text{e}^{\pm}$
for pions, $\pi^{\pm}$ for kaons, and $\mathrm{K}^{\pm}+\text{e}^{\pm}$ for
(anti)protons). The background contribution is small ($<1\%$) in all cases,
except for $\mathrm{K}^{\pm}$ at $p_{\mathrm{T}}>0.55\text{\,GeV/$c$}$, where
it reaches $\sim 13\%$. The TOF analysis uses the sub-sample of global tracks
for which the time-of-flight measurement is available. The procedure is
performed in narrow regions of pseudorapidity, $|\eta|<0.2$ and
$0.2<|\eta|<0.4$, in order to achieve a sufficient level of separation and to
strengthen the correlation between the total momentum and the transverse
momentum. The TOF matching efficiency for the presented data sample in the
pseudorapidity region $|\eta|<0.2$ ($0.2<|\eta|<0.4$) increases rapidly with
$p_{\mathrm{T}}$ up to around $50\%$ $(60\%)$ for pions at
$p_{\mathrm{T}}{\sim}700\text{\,MeV/$c$}$, $45\%$ $(55\%)$ for kaons at
$p_{\mathrm{T}}{\sim}1\text{\,GeV/$c$}$, and $50\%$ $(65\%)$ for protons at
$p_{\mathrm{T}}{\sim}800\text{\,MeV/$c$}$; it saturates at higher momenta
[56]. The PID performance of the TOF detector is shown in the bottom right
panel of Fig. LABEL:fig:PidPerformance, where the velocity $\beta$ of the
particles is reported as a function of momentum $p$. The raw particle yields
are then obtained by fitting the measured $\beta$ distributions with a TOF
response function, where the contribution for a given species is centered
around $\beta=p/E$. An additional template is added to account for wrongly
associated (mismatched) hits in the TOF detector.

The TPC-TOF fits analysis includes uncertainties related to the PID procedure
from several sources. For the TPC part, uncertainties are estimated by
integrating the ${\rm d}^{2}N/{\rm d}p_{\rm T}{\rm d}n_{\sigma}$ distribution
for each particle species in $|n_{\sigma}|<2.5$ and $|n_{\sigma}|<3.5$, and
comparing to the integral of the underlying distribution performed within
$|n_{\sigma}|<3$; the larger resulting uncertainty is used. Similarly, for TOF
fits for $p_{\mathrm{T}}<1.5\text{\,GeV/$c$}$, the uncertainty related to PID
is estimated by integrating the ${\rm d}^{2}N/{\rm d}p_{\rm T}{\rm
d}n_{\sigma}$ in the range of $|n_{\sigma}|<3$. This results in uncertainties
of up to $1\%$ for $\pi^{\pm}$, $5\%$ for $\mathrm{K}^{\pm}$, and $2\%$ for
(anti)protons. At higher $p_{\mathrm{T}}$, where the separation between
particle species becomes small, a better estimate of the uncertainty of the
method can be achieved by varying simultaneously the resolution
$\sigma_{\mathrm{TOF}}$ and the tail parameter [24, 56] of the fit function
used to describe the PID signal around their nominal values. An additional
uncertainty is included to account for the TOF miscalibration, which becomes
significant for $p_{\mathrm{T}}>1.5\text{\,GeV/$c$}$. This results in
uncertainties of up to $8\%$ for $\pi^{\pm}$, $14\%$ for $\mathrm{K}^{\pm}$,
and $7\%$ for (anti)protons. An uncertainty related to the matching of tracks
to TOF hits arises from differences in the TPC-TOF matching efficiency in real
data and simulations. This uncertainty is $3\%$ for $\pi^{\pm}$, $6\%$ for
$\mathrm{K}^{\pm}$, and $4\%$ for (anti)protons, independent of
$p_{\mathrm{T}}$.

####  4.2.3 HMPID

The HMPID analysis extends charged hadron identification into the
intermediate-$p_{\mathrm{T}}$ region ($2\text{\,GeV/$c$}\lesssim
p_{\mathrm{T}}\lesssim 10\text{\,GeV/$c$}$), combining the measurement of the
emission angle of the Cherenkov photons $\theta_{\rm Ch}$ and the momentum
information of the particle under study. The Cherenkov photons are selected
using the Hough Transform Method [57]. The measurement of the single photon
$\theta_{\rm Ch}$ angle in the HMPID requires the determination of the track
parameters, which are calculated for tracks propagated from the central
tracking detectors to the radiator volume where the Cherenkov photons are
emitted. Each track is extrapolated to the HMPID cathode plane and matched to
the closest primary ionization (MIP) cluster. The distance within the cathode
plane between the extrapolated track and the matched cluster (denoted
$d_{\rm{MIP-track}}$) is restricted to be less than 5 cm to reduce false
matches. A mean Cherenkov angle $\left<\theta_{\rm Ch}\right>$, computed as
the weighted average of single photon angles, is associated to each particle
track. The PID performance is shown in the bottom left panel of Fig.
LABEL:fig:PidPerformance, where the correlation between the reconstructed
Cherenkov angle and the track momentum is shown, indicating good agreement
with the theoretically expected values. For yield extraction, a statistical
unfolding technique is applied by fitting the reconstructed Cherenkov angle
distribution in a given momentum interval, which requires a precise knowledge
of the detector response function. Yields are evaluated from the integral of
each of the three Gaussian functions, corresponding to the signals from pions,
kaons, and protons. The HMPID allows pion and kaon identification in the
momentum range $1.5<p<4\text{\,GeV/$c$}$, while (anti)protons can be
distinguished from pions and kaons up to $p=6\text{\,GeV/$c$}$, with
separation powers larger than $2\sigma$.

Additionally, in the HMPID analysis, a data-driven correction for the
selection criterion on the distance $d_{\rm{MIP-track}}$ has been evaluated by
taking the ratio between the number of tracks that pass the selection
criterion on $d_{\rm{MIP-track}}$ and all the tracks in the detector
acceptance. This matching efficiency correction is $p_{\mathrm{T}}$ dependent
and it is about $20-40\%$; it is lower for particles with velocity
$\beta{\sim}1$ [52]. Negatively charged particle tracks have a distance
correction ${\sim}3\%$ lower than positive ones due to a radial residual
misalignment of the HMPID chambers and an imperfect estimation of the energy
loss in the material traversed by the track.

In the HMPID analysis, the systematic uncertainty has contributions from
tracking, PID, and track association [24, 52]. The PID uncertainties are
estimated by varying the parameters of the fit function used to extract the
raw particle yield. This uncertainty is $p_{\mathrm{T}}$-dependent and
increases with $p_{\mathrm{T}}$ to a maximum value of $12\%$ for $\pi^{\pm}$
and $\mathrm{K}^{\pm}$, and $11\%$ for (anti)protons. Furthermore, the
uncertainty of the association of the global track to the charged particle
signal in the HMPID is obtained by varying the default value of the
$d_{\rm{MIP-track}}$ distance criterion required for the matching. The
resulting uncertainty is $p_{\mathrm{T}}$-dependent with a maximum value of
about 4% for (anti)protons at $p_{\mathrm{T}}=1.5\text{\,GeV/$c$}$.

####  4.2.4 TPC relativistic rise

In the TPC $\text{d}E/\text{d}x$ relativistic rise analysis charged pions,
charged kaons, and (anti)protons can be identified up to
$p_{\mathrm{T}}=20\text{\,GeV/$c$}$. The identification is achieved by
measuring the specific energy loss $\text{d}E/\text{d}x$ in the TPC in the
relativistic rise regime of the Bethe–Bloch curve. The $\text{d}E/\text{d}x$
as a function of momentum $p$ is shown in the top right panel of Fig.
LABEL:fig:PidPerformance, indicating the $\left<\text{d}E/\text{d}x\right>$
response for charge-summed $\pi^{\pm}$, $\mathrm{K}^{\pm}$, ${\rm
p}({\rm\overline{p}})$, and $e^{\pm}$. The separation power between particle
species is about $4.5\sigma_{\mathrm{TPC}}$ ($1.5\sigma_{\mathrm{TPC}}$) for
$\pi-\text{p}$ ($\text{K}-\text{p}$) at $p_{\mathrm{T}}=10\text{\,GeV/$c$}$,
and it is nearly constant at similar values for larger momenta. The results
presented in this paper were obtained using the method detailed in Ref. [52].
As discussed in Refs. [7, 52, 24, 25], $\text{d}E/\text{d}x$ is calibrated
taking into account chamber gain variations, track curvature and diffusion to
obtain the best possible overall performance, which results in a response that
essentially only depends on $\beta\gamma$ ($=p/m$). The resolution is better
at larger rapidities for the same $\left<\text{d}E/\text{d}x\right>$ because
of the longer integrated track lengths. Hence, to analyze homogeneous samples,
the analysis is performed in four equal-width intervals within $|\eta|<0.8$
($|\eta|<0.2$, $0.2<|\eta|<0.4$, $0.4<|\eta|<0.6$ and $0.6<|\eta|<0.8$).
Samples of topologically identified pions from $\mathrm{K_{S}^{0}}$ decays,
protons from $\Lambda$ decays, and electrons from photon conversions were used
to parameterize the Bethe–Bloch response $\left<\text{d}E/\text{d}x\right>$ as
a function of $\beta\gamma$ and the relative resolution
$\sigma_{\mathrm{TPC}}/\left<\text{d}E/\text{d}x\right>$ as a function of
$\left<\text{d}E/\text{d}x\right>$. The relative yields of pions, kaons,
protons, and electrons are obtained as the $\pi^{+}+\pi^{-}$, ${\rm
K}^{+}+{\rm K}^{-}$, ${\rm p}+{\rm\overline{p}}$, and ${\rm e}^{+}+{\rm
e}^{-}$ yields normalized to that for inclusive charged particles. They are
obtained using four-Gaussian fits to $\text{d}E/\text{d}x$ distributions
differentially in $p$ and $|\eta|$ intervals. The parameters (mean and width)
of the fits are fixed using the parameterized Bethe–Bloch and resolution
curves. The relative yields as a function of $p_{\mathrm{T}}$ are found to be
independent of $\eta$ and therefore averaged. Particle yields are constructed
using the corrected relative yields and the corrected charged particle yields
[54]. A Jacobian correction is applied to account for the pseudorapidity-to-
rapidity conversion.

In the TPC $\text{d}E/\text{d}x$ relativistic rise analysis, the pion and
(anti)proton yields are corrected for secondary particles from weak decays
using MC simulations for the relative fraction of secondaries. The obtained
fraction of secondary pions and (anti)protons are scaled to those extracted
from $\text{DCA}_{xy}$ template fits to data. For $p_{\mathrm{T}}\gtrsim
3\text{\,GeV/$c$}$, the correction is negligible for pions. It is $\sim 2\%$
for (anti)protons at $p_{\mathrm{T}}=3\text{\,GeV/$c$}$, decreases to $\sim
1\%$ at $p_{\mathrm{T}}=10\text{\,GeV/$c$}$, and stays constant from that
$p_{\mathrm{T}}$ onward. Moreover, at high $p_{\mathrm{T}}$, there is a small
contamination of primary muons in the pion yields. Due to the similar muon and
pion masses, the electron (fractional) yield is subtracted from the pion yield
to correct for the muon contamination. This procedure gives a $<1\%$
correction to the pion yield in the entire $p_{\mathrm{T}}$ range considered
in the analysis. Furthermore, above $p_{\mathrm{T}}=3\text{\,GeV/$c$}$, both
the contamination of kaons and the contamination of (anti)deuterons in the
(anti)proton sample are negligible.

The tracking efficiency component is calculated as a relative correction
factor. It is the ratio of the inclusive to identified charged particle
efficiencies, and is applied to the relative yields. At high $p_{\mathrm{T}}$
this correction is nearly constant, of the order of $3-6\%$, depending on the
particle species.

In the TPC $\text{d}E/\text{d}x$ relativistic rise analysis, the systematic
uncertainties mainly originate from event and track selection and the PID
procedure. The first component is based on the study of inclusive charged
particles [54], and it was recalculated to meet the event selection condition
for inelastic events. Its value is estimated to be $7.3\%$ at high
$p_{\mathrm{T}}$ (at which the value is largest). It is the main contribution
for pions. The second component was measured following the procedure explained
in Ref. [52]. Here, the largest contribution is due to the uncertainties in
the parameterization of the $\text{d}E/\text{d}x$ response, resulting in
uncertainties of $6.5-15.0\%$ for $\mathrm{K}^{\pm}$ and $13.0-17.0\%$ for
(anti)protons, depending on $p_{\mathrm{T}}$.

####  4.2.5 Topological reconstruction of $\mathrm{K}^{\pm}$ kink decays

Charged kaons are also measured by reconstructing the vertex of their weak
decay in the TPC. The procedure extends the $p_{\mathrm{T}}$ reach of the
identification of charged kaons on a track-by-track basis from
$4\text{\,GeV/$c$}$ (available with the HMPID) up to $7\text{\,GeV/$c$}$. This
method exploits the characteristic kink topology defined by the decay of a
charged mother particle to a daughter with the same charge and a neutral
daughter [24]. Thanks to the two-body kinematics of the kink topology, it is
possible to separate kaon decays from the background mainly caused by pion
decays. For this purpose, a topological selection is applied by imposing a
selection criterion on the daughter track’s momentum with respect to that of
the mother track, and on the decay angle, defined as the angle between the
momenta of the mother and the charged daughter track. Furthermore, mother
tracks are selected inside a $3.5\sigma_{\mathrm{TPC}}$ band of the expected
$\text{d}E/\text{d}x$ for kaons to enhance the purity of the sample. With the
assumption that the charged daughter track is a muon and the undetected
neutral daughter particle is a neutrino, the reconstructed invariant mass
$M_{\mu\nu}$ is calculated and is shown in Fig. 1. The raw yield of the
topologically selected kaons in a given $p_{\mathrm{T}}$ bin are obtained from
the integral of the invariant mass distribution after the topological
selection criteria. The contamination due to fake kinks increases with
$p_{\mathrm{T}}$ and saturates at $p_{\mathrm{T}}{\sim}1\text{\,GeV/$c$}$,
reaching a maximum value of about 5%.

For the topological identification of charged kaons, the size of the
correction related to contamination (arising from background or
misidentification) was assigned as a $p_{\mathrm{T}}$-dependent uncertainty on
the purity. The kink identification uncertainty is $p_{\mathrm{T}}$ dependent
and it ranges from 2.5% at low $p_{\mathrm{T}}$ to 2.2% at high
$p_{\mathrm{T}}$. The systematic uncertainty on the efficiency for findable
kink vertices was estimated to be $3\%$, independent of $p_{\mathrm{T}}$. The
uncertainty due to contamination from fake kinks was estimated to be at most
about 5% around $p_{\mathrm{T}}=2\text{\,GeV/$c$}$, decreasing towards higher
$p_{\mathrm{T}}$.

Figure 1: Kink invariant mass distribution for charge summed particles in the
mother particle transverse momentum interval
$0.2\,\text{GeV}/c<p_{\mathrm{T}}<6\,\text{GeV}/c$, before (open circle) and
after (full circle) topological selections.

####  4.2.6 Combination of $\pi^{\pm}$, $\mathrm{K}^{\pm}$, and ${\rm
p}({\rm\overline{p}})$ spectra from different analyses

The charged pion, charged kaon and (anti)proton transverse momentum spectra
were measured via several independent analyses as described in the preceding
sections. To ensure the maximal $p_{\mathrm{T}}$ coverage, the final
$p_{\mathrm{T}}$ spectra were calculated as the average of all analyses
weighted by the systematic uncertainties that are not shared between analyses,
i.e. uncorrelated. The uncertainties related to the ITS-TPC matching
efficiency and the global tracking efficiency are largely correlated and were
summed in quadrature with the uncorrelated part of the systematic
uncertainties obtained after the averaging. Only the TPC relativistic rise
analysis is used above $p_{\mathrm{T}}=4,7$, and $6\text{\,GeV/$c$}$
respectively for pions, kaons, and protons. To verify the validity of the
procedure, the spectra obtained from the individual analyses were compared to
the final combined ones. The left panel of Fig. 2 shows the $\pi^{\pm}$,
$\mathrm{K}^{\pm}$, and ${\rm p}({\rm\overline{p}})$ spectra obtained from the
five analyses discussed above, which are normalized to the number of inelastic
collisions ($N_{\rm INEL}$). The right panel of Fig. 2 shows the ratios of the
individual spectra to the combined spectra, which illustrates an excellent
agreement in the overlapping $p_{\mathrm{T}}$ regions for every particle
species within the uncorrelated part of the systematic uncertainties.

Figure 2: Left panel: $p_{\mathrm{T}}$ spectra of $\pi^{\pm}$,
$\mathrm{K}^{\pm}$, and ${\rm p}({\rm\overline{p}})$ measured at midrapidity
($|y|<0.5$) in pp collisions at $\sqrt{s}=13\text{\,TeV}$ using different PID
techniques. The spectra are normalized to the number of inelastic collisions.
Right panel: The ratios of individual spectra to the combined spectra as a
function of $p_{\mathrm{T}}$ for $\pi^{\pm}$ (top), $\mathrm{K}^{\pm}$
(middle), and ${\rm p}({\rm\overline{p}})$ (bottom). Only the
$p_{\mathrm{T}}$-range where the analyses overlap is shown. The vertical bars
indicate statistical uncertainties while the bands show the uncorrelated
systematic uncertainties.

###  4.3 Topological identification of weakly-decaying strange hadrons

Primary strange hadrons $\mathrm{K_{S}^{0}}$, $\Lambda$, $\overline{\Lambda}$,
$\Xi^{-}$, $\overline{\Xi}^{+}$, $\Omega^{-}$, and $\overline{\Omega}^{+}$ are
reconstructed at midrapidity ($|y|<$0.5) via their characteristic weak decay
topologies in the channels presented in Tab. 1. Single-strange hadrons
$\mathrm{K_{S}^{0}}$, $\Lambda$ and $\overline{\Lambda}$ decay into two
oppositely charged daughter particles ($V^{0}$ decay). Multi-strange hadrons
($\Xi$ and $\Omega$) decay into a charged meson (bachelor) plus a $V^{0}$
decaying particle, giving the two-step process known as a cascade. The
identification methods for the $V^{0}$ ($\mathrm{K_{S}^{0}}$ and $\Lambda$)
and cascade-like ($\Xi$ and $\Omega$) candidates strictly follow those
presented in earlier works [58, 27]. An additional selection criterion on
(anti)proton momentum ($p>0.31\text{\,GeV/$c$}$) measured at the inner wall of
TPC was introduced because of the observed instability of track reconstruction
for the lowest $p_{\mathrm{T}}$ bin of $\Lambda$ and $\overline{\Lambda}$
spectra. Several track, PID, and topological selection criteria are applied in
order to find $V^{0}$ and cascade decay candidates. Charged tracks are
selected using the standard criteria described in Sec. 3.2. The identity of
these daughter tracks is established with the requirement that the specific
energy loss $\text{d}E/\text{d}x$ measured in the TPC is compatible with the
expected mass hypothesis within $5\sigma_{\mathrm{TPC}}$
($4\sigma_{\mathrm{TPC}}$) for the analysis of $\mathrm{K_{S}^{0}}$ and
$\Lambda$ ($\Xi$ and $\Omega$). These identified tracks are then combined to
form invariant mass distributions, and fake combinations are reduced by
applying selection criteria on topological variables.

Table 3: Selection criteria for secondary and bachelor tracks as well as for $V^{0}$ and cascade candidates applied in the presented work. Hadron species |  | $\mathrm{K_{S}^{0}}$ | $\Lambda$ | $\Xi$ | $\Omega$
---|---|---|---|---|---
Secondary track selections |  |  |  |
Pseudorapidity range $|\eta|$ |  $<$ 0.8 |  $<$ 0.8 |  $<$ 0.8 |  $<$ 0.8
DCAxy of $V^{0}$ daughter track |  |  |  |
to primary vertex (cm) |  $>$ 0.06 |  $>$ 0.06 |  $>$ 0.04 |  $>$ 0.03
DCAxy of bachelor track |  |  |  |
to primary vertex (cm) | — | — |  $>$ 0.05 |  $>$ 0.05
TPC $\text{d}E/\text{d}x$ PID selection ($n_{\sigma}$) |  $<$ 5 |  $<$ 5 |  $<$ 4 |  $<$ 4
$V^{0}$ selections |  |  |  |
Rapidity range $|y|$ |  $<$ 0.5 |  $<$ 0.5 | — | —
Transverse decay radius (cm) |  $>$ 0.5 |  $>$ 0.5 |  $>$ 1.4 |  $>$ 1.4
DCAxy of $V^{0}$ to primary vertex (cm) | — | — |  $>$ 0.07 |  $>$ 0.07
DCA between $V^{0}$ daughter tracks ($n_{\sigma}$) |  |  $<$ 1 |  $<$ 1 |  $<$ 1.5 |  $<$ 1.5
Cosine of Pointing Angles |  $>$ 0.97 |  $>$ 0.995 |  $>$ 0.97 |  $>$ 0.97
$|\Delta m|$ around nominal $\Lambda$ mass (MeV/$c^{2}$) | — | — |  $<$ 6 |  $<$ 6
Cascade selections |  |  |  |
Rapidity range $|y|$ | — | — |  $<$ 0.5 |  $<$ 0.5
Transverse decay radius (cm) | — | — |  $>$ 0.8 |  $>$ 0.6
DCA between $V^{0}$ and bachelor track (cm) | — | — |  $<$ 1.6 |  $<$ 1
Cosine of Pointing Angle | — | — |  $>$ 0.97 |  $>$ 0.97
$|\Delta m|$ around nominal $\Xi$ mass (MeV/$c^{2}$) | — | — |  $<$ 8 |  $>$ 8

Values for these selection criteria are summarized in Tab. 3 and a detailed
description can be found in Ref. [58]. $\mathrm{K_{S}^{0}}$ ($\Lambda$)
candidates compatible with the alternative $V^{0}$ hypothesis, obtained by
changing the mass assumption for the daughter tracks accordingly, are rejected
if they lie within a fiducial window around the nominal $\Lambda$
($\mathrm{K_{S}^{0}}$) mass. A similar selection is applied for the $\Omega$
analysis where candidates are rejected if the corresponding invariant mass,
obtained by assuming the pion mass hypothesis for the bachelor track, is
compatible within $\pm 8\ \mathrm{MeV}/c^{2}$ with the nominal $\Xi$ mass. A
selection is also made on the proper lifetime $c\tau=mL/p$, where $m$ is the
particle mass and $L$ is the distance from the primary vertex to the decay
vertex; $c\tau$ is required to be less than $20$ cm/$c$. A further selection
is applied on the pointing angle $\Theta$, the angle between the strange
hadron’s momentum vector and the position vector of its decay point with
respect to the primary collision vertex. With requirements of
$\cos\Theta>0.97$ for $\mathrm{K_{S}^{0}}$ and $\cos\Theta>0.995$ for
$\Lambda$ and $\overline{\Lambda}$, about 1% of secondary $\Lambda$ and
$\overline{\Lambda}$ generated in the detector material is removed.

The particle yields are obtained as a function of $p_{\mathrm{T}}$ by
extracting the signals from the relevant invariant mass distributions.
Examples of the invariant mass peaks at $\sqrt{s}=13\text{\,TeV}$ are shown in
Fig. LABEL:fig:InvMass; the distributions of $\mathrm{K_{S}^{0}}$ and
$\Lambda$ are very similar to those at $\sqrt{s}=7\text{\,TeV}$. The mean
($\mu$) and the width ($\sigma$) values of the distributions are found by
fitting the distribution with a Gaussian for the signal plus a linear function
describing the background. The extracted $\mu$ values of the distributions
both for $V^{0}$s and cascades are in good agreement with the accepted values
[44] and are well reproduced by MC simulations at $\sqrt{s}=7$ and
$\sqrt{s}=13\text{\,TeV}$ in all measured $p_{\mathrm{T}}$ bins. The widths of
the distributions evolve with $p_{\mathrm{T}}$ at $\sqrt{s}=13$ (7) TeV by
about 7 (14) MeV/$c^{2}$ for $\mathrm{K_{S}^{0}}$ and 2 (4) MeV/$c^{2}$ for
$\Lambda$ and $\overline{\Lambda}$, which agrees with MC simulations within
$15-20\%$ at both reported energies. For the $V^{0}$s (cascades) a region
containing all the signal (signal region) is defined around the mean within
$\pm~{}6\sigma$ ($3\sigma$), while a region to estimate the background
(background region) is defined as side-bands from -12$\sigma$ to -6$\sigma$
(-12$\sigma$ to -6$\sigma$) and from 6$\sigma$ to 12$\sigma$ (6$\sigma$ to
19$\sigma$). Given the flatness of the invariant mass distribution in the
background region, the estimate of the background in the signal region is
obtained rescaling, by the ratio of the widths in the two regions, the sum of
the entries of all the bins in the background region. An alternative method,
used to estimate a possible systematic uncertainty, uses the fit of the
invariant mass distribution in the background region to estimate the
background contribution inside the peak region. In both cases, the signal is
obtained subtracting the estimated background in the peak region from the
integrated counts in the peak region.

The yields for $\Lambda$ ($\overline{\Lambda}$) are significantly affected by
secondary particles coming from the decays of $\Xi^{-}$ ($\overline{\Xi}^{+}$)
and $\Xi^{0}$ ($\overline{\Xi}^{0}$). The feed-down fraction is computed for
each $p_{\mathrm{T}}$ bin as the detection efficiency of $\Lambda$
($\overline{\Lambda}$) from $\Xi$ decays multiplied by the measured
$\overline{\Xi}^{+}$ ($\Xi^{-}$) spectra, thereby assuming that the production
rates of charged and neutral $\Xi$ are equal. The ratio of secondary $\Lambda$
to the measured primary yield is about $6-22\%$, depending on
$p_{\mathrm{T}}$.

For the $V^{0}$ and cascade analyses, the main sources of uncertainties are
listed in Tab. LABEL:tab:syst:lambdak0s and Tab. LABEL:tab:syst:xiomega,
respectively. The values of the selection criteria on the topological
variables were varied around their nominal values and the observed deviations
for each component were summed in quadrature. Uncertainties related to signal
extraction were estimated by varying the values of the width of the signal and
background sampling regions with respect to the default one, and adopting an
alternative method to estimate the background in the peak region. The
procedure resulted in a $p_{\mathrm{T}}$-dependent uncertainty which rages
from 0.2% to 6.8% for $V^{0}$s and from 0.9% to 2.4% for the cascades.

For the measurement of $\mathrm{K_{S}^{0}}$, $\Lambda$ and
$\overline{\Lambda}$ at $\sqrt{s}=13\text{\,TeV}$, contributions from out-of-
bunch pileup are removed as stated in Sec. 3.2. The applied correction reaches
a maximum value of about $2\%$ for $\mathrm{K_{S}^{0}}$ and about $3\%$ for
$\Lambda+\overline{\Lambda}$ at high $p_{\mathrm{T}}$. For
$\overline{\Xi}^{+}$ and $\Xi^{-}$ at $\sqrt{s}=13\text{\,TeV}$, a
$p_{\mathrm{T}}$-dependent correction factor, taken from Ref. [28] is applied
to remove the out-of-bunch pileup contribution; this correction is $0.5\%$ at
low $p_{\mathrm{T}}$ and rises to $2.1\%$ at high $p_{\mathrm{T}}$. A similar
correction for $\Omega^{-}$ and $\overline{\Omega}^{+}$ is negligible and
hence not applied. The out-of-bunch pileup contribution is found to be
negligible for the measurement of $\mathrm{K_{S}^{0}}$, $\Lambda$, and
$\overline{\Lambda}$ at $\sqrt{s}=7\text{\,TeV}$. A systematic uncertainty to
account for any residual effect is $0.3-1.1\%$ for the
$\overline{\Xi}^{+}+\Xi^{-}$, $1.2-4.6\%$ for $\Lambda+\overline{\Lambda}$,
depending on $p_{\mathrm{T}}$, while it is negligible for
$\mathrm{K_{S}^{0}}$.

The $V^{0}$ (topological selection) and track selection criteria are varied
around their default values, producing up to $9.5\%$ ($5.8\%$) uncertainties
in the $\mathrm{K_{S}^{0}}$ and $\Lambda$ yields at $\sqrt{s}=13\text{\,TeV}$
($\sqrt{s}=7\text{\,TeV}$). The feed-down corrections for $\Lambda$ and
$\overline{\Lambda}$ carry uncertainties associated with the method and the
uncertainty of the measured $\Xi$ spectrum. The values are estimated to be
around $2\%$ ($4-6\%$) for 13 (7) TeV, depending on $p_{\mathrm{T}}$. The
applied selection criteria for PID in the TPC, used for better discrimination
between the combinatorial background and the signal for strange baryons, are
varied in the range of $4-7\sigma_{\mathrm{TPC}}$. For $\Lambda$ and
$\overline{\Lambda}$, this uncertainty is $\sim$0.8–6% at $\sqrt{s}=7$ and
0.2–3.3% at $\sqrt{s}=13\text{\,TeV}$. For $\Xi$ and $\Omega$, it is at most
1% for both collision energies. For $\mathrm{K_{S}^{0}}$ the difference was
found to be negligible at 13 TeV and at most 2.6% at 7 TeV. The resulting
total uncertainties from low to high $p_{\mathrm{T}}$ vary in the ranges
3.6–9.5% for $\mathrm{K_{S}^{0}}$, 6.1–7.8% for $\Lambda$, 6.6–3.7% for $\Xi$,
and 6.0–6.7% for $\Omega$.

In the invariant mass distribution (shown in Fig. LABEL:fig:InvMass), one
could observe for $V^{0}$s and cascades an imperfect fit of the tails close to
the peak; the excess of candidates comes from misidentified daughters. The
removal of this effect though application of further PID selection criteria
would reduce significantly the number of candidates while not improving the
signal extraction.

###  4.4 Reconstruction of resonances

Using analysis techniques similar to those described in Ref. [59],
$\mathrm{K^{*0}}$ and $\mathrm{\overline{K}^{*0}}$ mesons are identified
through reconstruction of their decays to charged pions and kaons, while
$\phi$ mesons are identified via their decays to pairs of charged kaons.
Primary charged tracks are selected using the standard criteria described in
Sec. 3.2. Pion and kaon candidates are identified using their specific energy
loss $\text{d}E/\text{d}x$ measured in the TPC and their velocity $\beta$
measured with the TOF. The specific energy loss for each pion (kaon) candidate
with $p>0.4\text{\,GeV/$c$}$ is required to be within 2$\sigma_{\mathrm{TPC}}$
of the expected mean value for pions (kaons); a less restrictive selection is
applied for lower momenta. In addition, if the charged track is matched to a
hit in the TOF, $\beta$ must be within 3$\sigma_{\mathrm{TOF}}$ of the
expected mean value. For the analysis of the
$\mathrm{K^{*0}}\kern-1.79993pt+\kern-0.89996pt\mathrm{\overline{K}^{*0}}$
($\phi$), all kaon candidates are paired with all oppositely charged pion
(kaon) candidates from the same event and the pair invariant mass is
calculated. The resulting distributions are shown in Fig. LABEL:fig:InvMass.
The combinatorial background is estimated by calculating the invariant mass
distribution of like-charge $\pi$K or KK pairs from the same event, by
parameterizing the combinatorial background with a simple function (for $\phi$
only), or by pairing tracks from two different events (the “mixed-event"
technique). In order to ensure that the mixed events have similar
characteristics, the $z$ positions of their primary vertices are required to
be separated by less than 1 cm and their charged particle multiplicities are
required to differ by no more than 5 reconstructed tracks. The background-
subtracted invariant mass distributions are then fitted with a peak function
added to a function that parameterizes the residual background contribution
from correlated pairs. The $\phi$ peak is described with a Voigtian function,
the convolution of a Breit–Wigner function and a Gaussian. The mass resolution
for the $\mathrm{K^{*0}}$ is much smaller than the width of that resonance and
a Breit–Wigner function is used to describe the $\mathrm{K^{*0}}$ peak. The
yield of the
$\mathrm{K^{*0}}\kern-1.79993pt+\kern-0.89996pt\mathrm{\overline{K}^{*0}}$
($\phi$) is then calculated by integrating the invariant mass distribution
within $0.8<m_{\pi\mathrm{K}}<1.0$
$\mathrm{GeV}\kern-1.79993pt/\kern-0.89996ptc^{2}$
$(1.01<m_{\mathrm{KK}}<1.03~{}\mathrm{GeV}\kern-1.79993pt/\kern-0.89996ptc^{2})$,
subtracting the integral of the residual background in the same region, and
adding the yield from outside the peak region obtained from the peak fit
function. The variations in the
$\mathrm{K^{*0}}\kern-1.79993pt+\kern-0.89996pt\mathrm{\overline{K}^{*0}}$ and
$\phi$ yields due to alternate treatments of the combinatorial background,
residual background parameterizations, and peak fitting functions are
incorporated into the systematic uncertainties.

The systematic uncertainties in the
$\mathrm{K^{*0}}\kern-1.79993pt+\kern-0.89996pt\mathrm{\overline{K}^{*0}}$ and
$\phi$ yields include the contributions listed in Tab.
LABEL:tab:syst:phikstar. The contribution due to variations in the PID, track,
and event selection criteria (“Track Selection" in the table) is 2–8% for
$\mathrm{K^{*0}}\kern-1.79993pt+\kern-0.89996pt\mathrm{\overline{K}^{*0}}$ and
1-5% for $\phi$. Variations in the combinatorial background construction,
residual background parameterization, peak parameterization, and fit range
(“Signal Extraction" in the table) combine to give an uncertainty of 5–11% for
$\mathrm{K^{*0}}\kern-1.79993pt+\kern-0.89996pt\mathrm{\overline{K}^{*0}}$ and
2–6% for $\phi$, depending on $p_{\mathrm{T}}$. The ITS-TPC matching
uncertainty for single particles (pions and kaons) is 2% for all
$p_{\mathrm{T}}$ intervals. The uncertainty in the branching ratio is $1\%$
for $\phi$ and negligible for
$\mathrm{K^{*0}}\kern-1.79993pt+\kern-0.89996pt\mathrm{\overline{K}^{*0}}$.
The total systematic uncertainty for
$\mathrm{K^{*0}}\kern-1.79993pt+\kern-0.89996pt\mathrm{\overline{K}^{*0}}$ and
$\phi$ is estimated to be about $10-18\%$ and $8-13\%$, depending on
$p_{\mathrm{T}}$.

##  5 Results

For the light-flavor hadrons discussed in this paper, the ratios of yields for
particles and antiparticles are around one within the uncertainties, as
expected at these collision energies in the midrapidity region. Therefore, all
the $p_{\mathrm{T}}$ spectra shown in the following are reported after summing
particles and antiparticles, when a distinct antiparticle state exists. Unless
explicitly stated, the sums of particles and antiparticles, $\pi^{+}+\pi^{-}$,
$\mathrm{K}^{+}+\mathrm{K}^{-}$, $\mathrm{K^{*0}}+\mathrm{\overline{K}^{*0}}$,
${\rm p}+{\rm\overline{p}}$, $\Lambda+\overline{\Lambda}$,
$\overline{\Xi}^{+}+\Xi^{-}$, and $\overline{\Omega}^{+}+\Omega^{-}$ are
denoted as $\pi^{\pm}$, $\mathrm{K}^{\pm}$,
$\mathrm{K^{*0}}(\mathrm{\overline{K}^{*0}})$, ${\rm p}({\rm\overline{p}})$,
$\Lambda(\overline{\Lambda})$, $\Xi^{-}(\overline{\Xi}^{+})$,
$\Omega^{-}(\overline{\Omega}^{+})$ or simply as $\pi$, K, $\mathrm{K^{*0}}$,
p, $\Lambda$, $\Xi$, and $\Omega$, respectively, unless explicitly written.

The uncertainty related to the overall normalization to inelastic (INEL)
events is fully correlated between particle species and is not shown
explicitly when plotting the results.

###  5.1 Transverse momentum distributions, integrated yields, and average
transverse momenta

Figure 3: Transverse momentum spectra of light-flavor hadrons measured at
midrapidity ($|y|<0.5$) in inelastic pp collisions at
$\sqrt{s}=13\text{\,TeV}$ (filled symbols) and $\sqrt{s}=7\text{\,TeV}$ (open
symbols, scaled by a factor of $1/2$) [25, 27, 60]. Statistical and systematic
uncertainties are shown as vertical error bars and boxes, respectively. The
data points are fitted using a Lévy–Tsallis function. The normalization
uncertainty of ${}^{+7.3}_{-3.5}\%$ $(\pm 2.6\%)$ for pp collisions at
$\sqrt{s}=7$ (13) TeV is common to all particle species and is excluded from
the plotted uncertainties.

The $p_{\mathrm{T}}$ spectra of light-flavor hadrons measured at midrapidity
in inelastic pp collisions at $\sqrt{s}=7$ and $13$ TeV are given in Fig. 3.
The $p_{\mathrm{T}}$ spectra of $\mathrm{K_{S}^{0}}$ and $\Lambda$ measured in
this paper at $\sqrt{s}=7\text{\,TeV}$ are shown together with other particle
species from previous ALICE measurements at the same center-of-mass energy
[27, 25, 60]. For $\sqrt{s}=7\text{\,TeV}$ INEL pp collisions, the reported
$p_{\mathrm{T}}$ distributions of $\pi$, K, p [27, 25], $\mathrm{K^{*0}}$, and
$\phi$ [60] are from the updated measurements of ALICE, with extended
$p_{\mathrm{T}}$ reach, and, for resonances, additionally with an improved
estimate of the systematic uncertainties. For clarity, some of the spectra
have been scaled with the factors indicated in the legends.

The $p_{\mathrm{T}}$ distributions are fitted with Lévy–Tsallis functions [61,
62] in order to extrapolate the spectra to the unmeasured $p_{\mathrm{T}}$
regions, i.e. down to zero and up to high $p_{\mathrm{T}}$, similar to what
was done in previous measurements [24, 58, 27, 26, 63]. This procedure allows
the $p_{\mathrm{T}}$-integrated yields $\text{d}N/\text{d}y$ and the average
transverse momenta $\langle p_{\mathrm{T}}\rangle$ to be extracted, for which
the measured as well as the extrapolated distributions are used. The obtained
values are given in Tab. 4. The fit function describes well both the
low-$p_{\mathrm{T}}$ exponential and high-$p_{\mathrm{T}}$ power-law nature of
the $p_{\mathrm{T}}$ distribution, with $\chi^{2}/\rm{ndf}$ values in the
range $0.2-2.3$. The slope parameter of the Lévy–Tsallis functions decreases
for all particle species going from $\sqrt{s}=7$ to $13$ TeV. For example, for
charged pions and (anti)protons it changes from $6.65\pm 0.03$ to $6.42\pm
0.03$ and from $7.99\pm 0.13$ to $7.71\pm 0.11$. No extrapolation is needed
for $\mathrm{K_{S}^{0}}$ and $\mathrm{K^{*0}}$, as their yields are measured
down to $p_{\mathrm{T}}=0\text{\,GeV/$c$}$. The fractions of extrapolated
particle yields outside the measured $p_{\mathrm{T}}$ range at low
$p_{\mathrm{T}}$ are given in the last column of Tab. 4.

Other fit ranges and other parameterizations ($m_{\mathrm{T}}$ exponential,
Boltzmann distribution, Bose–Einstein distribution, Fermi–Dirac distribution,
Boltzmann–Gibbs blast-wave function [64]) are also used and the resulting
variations in the $\text{d}N/\text{d}y$ and $\langle p_{\mathrm{T}}\rangle$
values are incorporated into the systematic uncertainties. The systematic
uncertainties are 4–10% for $\text{d}N/\text{d}y$ and 1–3% for $\langle
p_{\mathrm{T}}\rangle$, similar to those estimated for measurements during Run
1. For the present measurements, systematic uncertainties are dominant. The
average yields for all particles species increase with collision energy.
Compared to pp collisions at $\sqrt{s}=7\text{\,TeV}$ [24, 27], the average
increase of the $\langle p_{\mathrm{T}}\rangle$ and the yield per inelastic
collision for all measured particle species is about $8\%$ and $11\%$,
respectively. This is in agreement with the ${\approx}15\,\%$ increase of the
average pseudorapidity density of charged particles produced in $|\eta|<0.5$
as the collision energy increases from $\sqrt{s}=7\text{\,TeV}$ to
$\sqrt{s}=13\text{\,TeV}$ [54].

Table 4: Average transverse momentum $\langle p_{\mathrm{T}}\rangle$ and
$p_{\mathrm{T}}$-integrated yield $\text{d}N/\text{d}y$ values for light-
flavor hadrons measured in pp collisions at $\sqrt{s}=13\text{\,TeV}$ and 7
TeV. The first uncertainty is statistical and the second is systematic
(including the low- and high-$p_{\mathrm{T}}$ extrapolation). The
uncertainties due to the normalization to the number of inelastic events,
which are ${}^{+7.3}_{-3.5}\%$ and $\pm 2.6\%$ for pp collisions at
$\sqrt{s}=7$ and 13 TeV, respectively, are not included. The last column
represents the fraction of extrapolated yield at low transverse momenta.

Hadron species | $\text{d}N/\text{d}y$ | $\langle p_{\mathrm{T}}\rangle$ (GeV/$c$) | Extrapolated
---|---|---|---
fraction (%)
Collision energy: $\sqrt{s}=13\text{\,TeV}$
$\pi^{+}+\pi^{-}$ | $4.775\pm 0.001\pm 0.243$ | $(4.915\pm 0.001\pm 0.099)\times 10^{-1}$ | $8\pm 1$
${\rm K}^{+}+{\rm K}^{-}$ | $(6.205\pm 0.004\pm 0.303)\times 10^{-1}$ | $(8.099\pm 0.007\pm 0.099)\times 10^{-1}$ | $9\pm 0.1$
$\text{K}^{0}_{\text{S}}$ | $(3.192\pm 0.004\pm 0.111)\times 10^{-1}$ | $(8.207\pm 0.008\pm 0.087)\times 10^{-1}$ | negl.
$\mathrm{K^{*0}}\kern-1.79993pt+\kern-0.89996pt\mathrm{\overline{K}^{*0}}$ | $(2.098\pm 0.016\pm 0.200)\times 10^{-1}$ | $1.121\pm 0.005\pm 0.030$ | negl.
${\rm p}+\overline{\rm p}$ | $(2.750\pm 0.002\pm 0.188)\times 10^{-1}$ | $(9.659\pm 0.008\pm 0.144)\times 10^{-1}$ | $11\pm 0.5$
$\phi$ | $(3.734\pm 0.040\pm 0.213)\times 10^{-2}$ | $1.236\pm 0.009\pm 0.027$ | $13\pm 1$
$\Lambda+\overline{\Lambda}$ | $(1.807\pm 0.005\pm 0.102)\times 10^{-1}$ | $1.078\pm 0.002\pm 0.030$ | $22\pm 12$
$\Xi^{-}+\overline{\Xi}^{+}$ | $(1.980\pm 0.012\pm 0.082)\times 10^{-2}$ | $1.296\pm 0.004\pm 0.015$ | $20\pm 4$
$\Omega^{-}+\overline{\Omega}^{+}$ | $(1.846\pm 0.046\pm 0.122)\times 10^{-3}$ | $1.527\pm 0.023\pm 0.027$ | $34\pm 7$
Collision energy: $\sqrt{s}=7\text{\,TeV}$
$\text{K}^{0}_{\text{S}}$ | $(2.802\pm 0.002\pm 0.149)\times 10^{-1}$ | $(7.731\pm 0.006\pm 0.100)\times 10^{-1}$ | negl.
$\Lambda+\overline{\Lambda}$ | $(1.523\pm 0.002\pm 0.110)\times 10^{-1}$ | $1.028\pm 0.001\pm 0.019$ | $16\pm 1$

##  6 Discussion

Figure 4: Ratios of the transverse-momentum spectra of light-flavor hadrons in
inelastic pp collisions at $\sqrt{s}=13\text{\,TeV}$ to those for
$\sqrt{s}=7\text{\,TeV}$ [25, 27, 60]. The ratio for $\pi^{\pm}$ is shown in
each panel with grey crosses and boxes. Statistical and systematic
uncertainties are shown as vertical bars and boxes, respectively. The
normalization uncertainty (${}^{+10.8}_{-6.3}\%$) is excluded from the plotted
uncertainties.

The transverse-momentum spectra reported in Fig. 3 indicate a progressive and
significant evolution of the spectral shapes at high $p_{\mathrm{T}}$ with
increasing collision energy, which is similar for all particle species under
study. This behavior is better visualized in Fig. 4 which shows the
corresponding ratios of $p_{\mathrm{T}}$ spectra at $\sqrt{s}=13\text{\,TeV}$
to those at $\sqrt{s}=7\text{\,TeV}$ [25, 27, 60]. The systematic
uncertainties at both collision energies are largely uncorrelated and
therefore their sum in quadrature is taken as the systematic uncertainty on
the ratios. The uncertainty on the ratio due to normalization is
${}^{+10.8}_{-6.3}\%$.

The ratios for all hadron species are above unity, which is consistent with
the observed increase in the pseudorapidity density of inclusive charged
particles with increasing collision energy [54]. Furthermore, all the ratios
exhibit a clear increase as a function of $p_{\mathrm{T}}$, indicating that
hard processes become dominant in the production of high-$p_{\mathrm{T}}$
particles. The $p_{\mathrm{T}}$ dependence demonstrates that the spectral
shapes are significantly harder at $\sqrt{s}=13\text{\,TeV}$ than at
$\sqrt{s}=7\text{\,TeV}$, which is also evident in the reduction of the slope
parameter of the Lévy–Tsallis functions mentioned above. A universal shape —
independent of $p_{\mathrm{T}}$ within uncertainties — can be observed for
most species (excluding $\Xi$ and $\Omega$) in the soft regime,
$p_{\mathrm{T}}\lesssim 1\text{\,GeV/$c$}$. There is a hint that the ratio for
${\rm p}({\rm\overline{p}})$ may be enhanced above the one for $\pi^{\pm}$ in
the $p_{\mathrm{T}}$ region ${\sim}3-6\text{\,GeV/$c$}$, although the
enhancement is barely significant given the uncertainties. Such an enhancement
would be consistent with the appearance of the baryon anomaly, an increased
baryon-to-meson production ratio at intermediate transverse momenta
($2\text{\,GeV/$c$}\lesssim p_{\mathrm{T}}\lesssim 10\text{\,GeV/$c$}$),
observed in previous ALICE measurements of light-flavor hadron production [7,
65, 53, 25, 66]. It is worth noting that the hardening of the $p_{\mathrm{T}}$
spectra with increasing collision energy has been reported in our earlier work
for inclusive charged particles [54], although with different event selection
criteria. There, a requirement of at least one charged particle with
$p_{\mathrm{T}}>0\text{\,GeV/$c$}$ in $|\eta|<1$ was imposed, selecting events
corresponding to 75% of the total inelastic cross section. In Ref. [54], the
observed trend was found to be well captured by the PYTHIA and EPOS-LHC MC
generators. In Sec. 6.5 the ratios of $p_{\mathrm{T}}$ spectra of light-flavor
hadrons at $\sqrt{s}=13\text{\,TeV}$ to those at $\sqrt{s}=7\text{\,TeV}$ are
compared to results from these common event generators.

###  6.1 Scaling properties of hadron production

Two kinds of universal scaling of identified particle production have been
observed in high energy pp collisions: transverse mass ($m_{\mathrm{T}}$)
scaling, which was originally seen in the lower $p_{\mathrm{T}}$ region of
hadron spectra, and $x_{\mathrm{T}}$ scaling [67, 68, 69, 16, 70], observed in
the higher $p_{\mathrm{T}}$ region. New studies of the $m_{\mathrm{T}}$ and
$x_{\mathrm{T}}$ scaling properties of light-flavor hadrons in pp collisions
at $\sqrt{s}=7$ and 13 TeV are discussed below.

####  6.1.1 Transverse mass ($m_{\mathrm{T}}$) scaling

At ISR energies [71, 72] it has been observed that hadron $m_{\mathrm{T}}$
spectra in pp collisions seem to follow an approximately universal curve after
scaling with arbitrary normalization factors, an effect known as
$m_{\mathrm{T}}$ scaling. The STAR and PHENIX collaborations observed the
breaking of $m_{\mathrm{T}}$-scaling in pp collisions at
$\sqrt{s}=200\text{\,GeV}$ [21, 19], where a clear separation between baryon
and meson spectra was observed for $m_{\mathrm{T}}\geq 2$
$\mathrm{GeV}\kern-1.79993pt/\kern-0.89996ptc^{2}$. In pp collisions at
$\sqrt{s}=200\text{\,GeV}$, the separation between the baryon and meson
spectra seems to increase over the measured $m_{\mathrm{T}}$ range. We, the
ALICE collaboration, recently reported a similar breaking in pp collisions at
$\sqrt{s}=8\text{\,TeV}$ [10]. The measurements in this paper allow this study
to be extended to higher $m_{\mathrm{T}}$ and to the highest LHC energies.

The charged kaon and (anti)proton $m_{\mathrm{T}}$ spectra are fitted
separately with a modified Hagedorn function of the form
$A\times(e^{-am_{\mathrm{T}}}+m_{\mathrm{T}}/b)^{-n}$ [22]. The meson
$m_{\mathrm{T}}$ spectra are then scaled by multiplicative factors so that
their integrals over the measured $m_{\mathrm{T}}$ ranges (or
$m_{\mathrm{T}}>4$ $\mathrm{GeV}\kern-1.79993pt/\kern-0.89996ptc^{2}$ for the
pions) match the integral of the kaon fit function. Similarly, the baryon
$m_{\mathrm{T}}$ spectra are scaled so that their integrals match the integral
of the proton fit function over their measured $m_{\mathrm{T}}$ ranges. The
baryon spectra are further scaled so that all fits to the $m_{\mathrm{T}}$
spectra have the same value at $m_{\mathrm{T}}=1$
$\mathrm{GeV}\kern-1.79993pt/\kern-0.89996ptc^{2}$. Figure 5 shows the result
of this study for pp collisions at $\sqrt{s}=7\text{\,TeV}$ [25, 27, 60] and
13 TeV, with the lower panels showing the ratios of the various scaled
$m_{\mathrm{T}}$ spectra to the $\mathrm{K}^{\pm}$ fit function.

The pion $m_{\mathrm{T}}$ spectra deviate from the trend followed by the other
mesons for $m_{\mathrm{T}}\lesssim 2$
$\mathrm{GeV}\kern-1.79993pt/\kern-0.89996ptc^{2}$, which is likely due to
feed-down from resonance decays [73]. The measured primary $\pi^{\pm}$ yield
contains a significant contribution mostly from $\rho$ and $\omega$ decays,
which, according to a recent study [23], affects the low-$p_{\mathrm{T}}$
($\lesssim 1\text{\,GeV/$c$}$) part of the spectrum with increasing importance
towards higher collision energies. This is the reason why the
$\mathrm{K}^{\pm}$ $m_{\mathrm{T}}$ spectrum is used as the reference for the
other mesons. A clear difference in the slope is observed between the baryon
and meson spectra for $m_{\mathrm{T}}\gtrsim 2$
$\mathrm{GeV}\kern-1.79993pt/\kern-0.89996ptc^{2}$. The separation between the
meson and baryon $m_{\mathrm{T}}$ spectra may be a reflection of the fact
that, according to the Lund model of hadronization, meson formation via the
fragmentation of strings requires the breakup to only a (quark, anti-quark)
pair, while a baryon–antibaryon pair can be formed by the (diquark, anti-
diquark) break up of the string [74, 75]. The separation between the
(anti)proton and meson $m_{\mathrm{T}}$ spectra becomes approximately constant
for $m_{\mathrm{T}}>10$ $\mathrm{GeV}\kern-1.79993pt/\kern-0.89996ptc^{2}$.

Figure 5: Upper panels: Scaled $m_{\mathrm{T}}$ spectra for identified
particles in pp collisions at $\sqrt{s}=$ 7 TeV (left) and 13 TeV (right).
Lower panels: Ratios of the scaled $m_{\mathrm{T}}$ spectra to the
$\mathrm{K}^{\pm}$ fit functions. The reference data for pp at
$\sqrt{s}=7\text{\,TeV}$ are from Refs. [25, 27, 60]

The breaking of $m_{\mathrm{T}}$ scaling for pions at low $m_{\mathrm{T}}$ as
shown in Fig. 5 serves as a motivation to quantify the effect for different
mesons. Following Refs. [23, 76], instead of using the $m_{\mathrm{T}}$
spectra, the scaling law can be better studied practically as a function of
$p_{\mathrm{T}}$, (note that the invariant yields are equal in terms of these
variables). This requires changing the functional form of the invariant yield
parameterization through the substitution
$p_{\mathrm{T}}\to\sqrt{m_{\mathrm{T}}^{2}-m^{2}}$, where $m$ is the particle
mass. In doing so, the $p_{\mathrm{T}}$-differential invariant yield
$Y_{s\textprime}$ of a particle species $s\textprime$ can be obtained by
scaling the parameterization of the yield $Y^{\rm ref}_{s}$ of reference
particle species $s$. When both spectra are evaluated at the same transverse
mass,
$(p_{\mathrm{T},s}^{2}+m_{s}^{2})=(p_{\mathrm{T},s\textprime}^{2}+m_{s\textprime}^{2})$,
the yield for species $s\textprime$ is
$Y_{s\textprime}(p_{\mathrm{T},s\textprime})=c\times Y^{\rm
ref}_{s}\left(\sqrt{p^{2}_{\mathrm{T},s\textprime}+m^{2}_{s\textprime}-m^{2}_{s}}\right)$.
Here, $c$ is a constant offset between the yields for species $s$ and
$s\textprime$, determined in the $p_{\mathrm{T}}$ region where the spectral
shapes for the two species are the same (i.e. at high $p_{\mathrm{T}}$).

Charged pions are used as the reference species to test the $m_{\mathrm{T}}$
scaling of the $\mathrm{K}^{\pm}$/$\pi^{\pm}$,
$\mathrm{K_{S}^{0}}$/$\pi^{\pm}$, $\mathrm{K^{*0}}$/$\pi^{\pm}$ and
$\phi$/$\pi^{\pm}$ particle ratios. The reference $p_{\mathrm{T}}$ spectrum of
$\pi^{\pm}$ is parameterized with a Lévy–Tsallis function, which is also used
for the extraction of the $p_{\mathrm{T}}$-integrated yields and is found to
describe the data within $15\%$ over its entire $p_{\mathrm{T}}$ range. The
$m_{\mathrm{T}}$ scaling relation is then applied to the parametrized pion
yield. The appropriate offset parameter $c$ is determined by fitting the
measured particle ratios separately in the high-$p_{\mathrm{T}}$ region where
they saturate. The measured $\mathrm{K}^{\pm}$, $\mathrm{K_{S}^{0}}$, and
$\mathrm{K^{*0}}$ ($\phi$) yield in pp collisions at $\sqrt{s}=$ 13 (7) TeV
are fitted with a constant linear function in the range
$p_{\mathrm{T}}>6\text{\,GeV/$c$}$ ($p_{\mathrm{T}}>4\text{\,GeV/$c$}$) with
statistical and systematic uncertainties summed in quadrature. The resulting
$c$ values for the various particle ratios are given in the legends of Fig. 6
and are shown as shaded bands in the $p_{\mathrm{T}}$ region from which they
are extracted. The reported widths of the bands correspond to the statistical
uncertainties on $c$ obtained from the fits.

Along with the measured particle yield ratios, Fig. 6 shows the ratios of the
$m_{\mathrm{T}}$-scaled parameterizations to the reference parametrization
(solid blue lines) for the particle species in question. A significant
deviation between the parameterized curves and the measured data is observed
in the low-to-mid $p_{\mathrm{T}}$ region ($p_{\mathrm{T}}\lesssim
5-6\text{\,GeV/$c$}$) for the $\mathrm{K}^{\pm}$/$\pi^{\pm}$ and
$\mathrm{K_{S}^{0}}$/$\pi^{\pm}$ ratios, which indicates the breaking of
empirical $m_{\mathrm{T}}$ scaling. In contrast, the $m_{\mathrm{T}}$-scaling
predictions for the $\mathrm{K^{*0}}/\pi^{\pm}$ and $\phi/\pi^{\pm}$ ratios
are notably closer to the measurements. In order to quantify the level of the
scale breaking, the double ratios between the measured yield ratios and those
based on the $m_{\mathrm{T}}$-scaling relation are evaluated. The double ratio
values for $p_{\mathrm{T}}$ below the threshold of $6\text{\,GeV/$c$}$ (i.e.,
below the region used for the determination of the offset parameter $c$)
decrease towards lower $p_{\mathrm{T}}$, deviating beyond $16\%$ for
$p_{\mathrm{T}}\lesssim 2\text{\,GeV/$c$}$. The significance of the deviation
from the $m_{\mathrm{T}}$ scaling hypothesis is $4.7\sigma$ at
$p_{\mathrm{T}}=1\text{\,GeV/$c$}$, far from the threshold (note that adjacent
$p_{\mathrm{T}}$ bins have fully uncorrelated uncertainties). This further
confirms the breaking of the empirical $m_{\mathrm{T}}$ scaling relation for
the quoted $p_{\mathrm{T}}$ region.

Figure 6: Particle ratios $\mathrm{K}^{\pm}$/$\pi^{\pm}$,
$\mathrm{K_{S}^{0}}$/$\pi^{\pm}$ (top) and $\mathrm{K^{*0}}$/$\pi^{\pm}$,
$\phi$/$\pi^{\pm}$ (bottom) as a function of $p_{\mathrm{T}}$ measured in pp
collisions at $\sqrt{s}=13\text{\,TeV}$. The measured ratio is reported
together with that obtained from transverse mass ($m_{\mathrm{T}}$) scaling of
charged pions shown as a solid line. The shaded red band indicates the
constant fit, in the $p_{\mathrm{T}}$ region where it was performed, which is
used to determine the constant offset parameter $c$. See the text for details.
The $\phi$ spectrum in pp at $\sqrt{s}=7\text{\,TeV}$ is from Ref. [60]

####  6.1.2 $x_{\mathrm{T}}$ scaling

The validity of empirical $x_{\mathrm{T}}$ scaling is tested using the
$\sqrt{s}=13\text{\,TeV}$ pp measurements reported here and those obtained at
$\sqrt{s}=2.76\text{\,TeV}$ [7, 63] and $\sqrt{s}=7\text{\,TeV}$ [25, 60]. Due
to the lack of a complete set of measurements of hadron yields at high
$p_{\mathrm{T}}$ at these three collision energies, only $\pi^{\pm}$,
$\mathrm{K}^{\pm}$, $\mathrm{K^{*0}}$, and ${\rm p}({\rm\overline{p}})$ are
used.

The invariant cross sections are determined from the measured invariant yields
as $E\,\mathrm{d}^{3}\sigma/\mathrm{d}p^{3}=\sigma_{\rm
inel}\times\,E\,\mathrm{d}^{3}N/\mathrm{d}p^{3}$, where $\sigma_{\rm inel}$ is
the inelastic cross section in pp collisions at $\sqrt{s}=13\text{\,TeV}$
[48]. The logarithm of the ratio of the invariant cross sections at two
different collision energies, scaled by the logarithm of the ratio of the two
collision energies is calculated. This quantity, denoted as $n$ in the
following, depends on $x_{\mathrm{T}}$ and $\sqrt{s}$ [77]. It increases with
$x_{\mathrm{T}}$ in the low $x_{\mathrm{T}}$ region, where particle production
is dominated by soft processes, and appears to saturate in the high
$x_{\mathrm{T}}$ region.

The $x_{\mathrm{T}}$ scaling exponents $n$ as a function of $x_{\mathrm{T}}$
are extensively studied in Ref. [16] for unidentified hadrons and pions using
the data from the E706 [78], PHENIX/ISR [79, 80], PHENIX [79, 81], UA1 [13]
and CDF [12, 82] collaborations. The value of $n$ decreases towards the
high-$x_{\mathrm{T}}$ region. The measured mean exponent value was about $8.2$
in the range $x_{\mathrm{T}}=10^{-1}-4\times 10^{-1}$ and reaches a value of
$\approx 5.2$ in the range $x_{\mathrm{T}}=10^{-2}-7\times 10^{-3}$. In our
measurements at $\sqrt{s}=13\text{\,TeV}$, we are able to probe an even lower
$x_{\mathrm{T}}$ region.

Each $n(x_{\mathrm{T}},\sqrt{s})$ distribution is fitted with a constant in
the range $2\times 10^{-3}\leq x_{\mathrm{T}}\leq 6\times 10^{-3}$ to obtain
the respective $n$ values for different energy combinations; these are then
averaged to obtain the mean value $\langle n\rangle$ for each particle
species. The $x_{\mathrm{T}}$ spectra for different particle species are
scaled by $(\sqrt{s}/\text{GeV})^{\langle n\rangle}$. Within the quoted
$x_{\mathrm{T}}$ range, the best scaling is achieved with the exponents
$\langle n\rangle=5.04\pm 0.02$ for $\pi^{\pm}$, $\langle
n\rangle=5.02^{+0.21}_{-0.25}$ for $\mathrm{K}^{\pm}$, $\langle
n\rangle=5.83^{+0.13}_{-0.21}$ for ${\rm p}({\rm\overline{p}})$, and $\langle
n\rangle=5.23\pm 0.15$ for $\mathrm{K^{*0}}$. The uncertainties on the
$\langle n\rangle$ values are the maximum observed deviations from the mean.
The scaling exponents for the meson spectra are found to be consistent within
systematic uncertainties, indicating that ratios of meson spectra at high
$p_{\mathrm{T}}$ are constant and attain similar values for all beam energies.
Figure 7 shows the $x_{\mathrm{T}}$-scaled spectra for $\pi^{\pm}$,
$\mathrm{K}^{\pm}$, ${\rm p}({\rm\overline{p}})$, and $\mathrm{K^{*0}}$ at
collision energies of $\sqrt{s}=2.76\text{\,TeV}$ [7, 63], 7 TeV [25, 60] and
13 TeV.

These results suggest that identified particle yields at LHC energies follow
$x_{\mathrm{T}}$ scaling above $x_{\mathrm{T}}{\sim}10^{-3}$. The measured
$x_{\mathrm{T}}$ scaling exponent at LHC energies is consistent with the
$x_{\mathrm{T}}$-dependent trend of the exponent values [16]. The value of the
parameter $n$ is found to be reduced at LHC energies compared to lower
collision energies [79, 80, 78, 81, 12, 82, 13, 83, 84, 85, 86, 14], which is
attributed to the increasing importance of hard scattering processes at higher
$\sqrt{s}$, as suggested by Fig. 4.

It is also interesting to note that the exponent $n$ takes on larger values
for baryons than for mesons in the investigated $x_{\mathrm{T}}$ range. On the
one hand, this is connected to the decrease of the ${\rm p}/\pi$ ratio with
increasing $p_{\mathrm{T}}$ (see Fig. 11), as opposed to the constant behavior
of the ${\rm K}/\pi$ and ${\rm K^{*0}}/\pi$ ratios (see Fig. 6), which
suggests that meson spectra are harder than baryon spectra in the
corresponding $p_{\mathrm{T}}$ range. On the other hand, as discussed in Ref.
[16], the NLO pQCD predictions including higher-twist processes, i.e. in which
the detected hadron can be exclusively produced in the hard subprocess
reaction, there is evidence for the larger value of the exponent for baryons
than for mesons. This is in contrast to the observations based on the leading-
twist processes, where the exponent $n$ has only a weak dependence on hadron
species.

Figure 7: Scaled invariant yields of $\pi^{\pm}$, $\mathrm{K}^{\pm}$, ${\rm
p}({\rm\overline{p}})$, and $\mathrm{K^{*0}}$ as a function of
$x_{\mathrm{T}}=2p_{\mathrm{T}}/\sqrt{s}$ at different collision energies of
$\sqrt{s}=2.76\text{\,TeV}$ [7, 63], $\sqrt{s}=7\text{\,TeV}$ [25, 60], and
$\sqrt{s}=13\text{\,TeV}$. The solid line represents a combined power-law fit
in the high-$x_{\mathrm{T}}$ region where the distributions show a scaling
behavior.

The quality of the scaling behavior is verified by combined fitting the
differential cross sections with a power-law function of the form $a\times
x_{\mathrm{T}}^{b}\times(1+x_{\mathrm{T}})^{c}$. Here, $a$, $b$, and $c$ are
free parameters and the region below $x_{\mathrm{T}}=1.5\times 10^{-3}$
($x_{\mathrm{T}}=2\times 10^{-3}$ for protons) is excluded to avoid the
dominant contribution from soft particle production, which does not follow
$x_{\mathrm{T}}$ scaling. The fits are of good quality with
$\chi^{2}/\rm{ndf}$ values in the range $0.4-1.5$. In spite of the naive
assumption of a power law function and the expected non-scaling behaviors
discussed in Ref. [87], the measurements agree with the global power law fits
within the region of overlap ($2\times 10^{-3}\lesssim x_{\mathrm{T}}\lesssim
6\times 10^{-3}$) within roughly $40\%$, depending on particle species. The
measurements from ALICE performed at $\sqrt{s}=13\text{\,TeV}$ are consistent
over the accessible $x_{\mathrm{T}}$ range ($2\times 10^{-3}\leq
x_{\mathrm{T}}\leq 6\times 10^{-3}$) with empirical $x_{\mathrm{T}}$ scaling
and with measurements from pp collisions at $\sqrt{s}=$ 2.76 and 7 TeV.

Figure 8:  Average transverse momentum $\langle p_{\mathrm{T}}\rangle$ as a
function of the center-of-mass energy. Open boxes indicate the statistical and
systematic uncertainties (when available) summed in quadrature. Results from
ALICE [51, 7, 24, 58, 63, 27, 26, 60] are compared with those from STAR
measured at $\sqrt{s}=200\text{\,GeV}$ [88]. Some data points are slightly
offset from their true energy for better visibility. Dashed curves show linear
fits in $\ln{s}$. Note that the data points of $\mathrm{K_{S}^{0}}$ are not
fitted due to their very similar values of $\langle p_{\mathrm{T}}\rangle$ to
those of charged kaons. Red data points are for $\sqrt{s}=13\text{\,TeV}$.

###  6.2 Excitation functions

Figure 8 compiles the excitation functions of the average transverse momenta
$\langle p_{\mathrm{T}}\rangle$ for light-flavor hadrons in inelastic pp
collisions; the focus is on the LHC energy regime, but some RHIC results are
also shown. As discussed in Ref. [89], the measurement of $\langle
p_{\mathrm{T}}\rangle$ as a function of the collision energy is particularly
useful in probing the saturation scale of the gluons inside the proton.

Results at midrapidity are presented from $\sqrt{s}=200\text{\,GeV}$ up to the
top LHC energy $\sqrt{s}=13\text{\,TeV}$, spanning nearly two orders of
magnitude in center-of-mass energy. The average $p_{\mathrm{T}}$ increases
with $\sqrt{s}$; it rises steeper for heavier particles, as seen in our
earlier measurements at lower collision energies [24, 27, 26, 60]. The growth
of $\langle p_{\mathrm{T}}\rangle$ with $\sqrt{s}$ is attributed to the
increasing importance of hard processes for higher collision energies. For
single- and multi-strange hadrons, this observation is equivalent to the
hardening of $\langle p_{\mathrm{T}}\rangle$ as the collision energy increases
from $\sqrt{s}=$ 7 to 13 TeV for event classes with a similar
$\left<\text{d}N_{\rm{ch}}/\text{d}\eta\right>$, as reported in Ref. [28].

Figure 9: Ratio of particle yields to $\pi^{\pm}$ and $\mathrm{K}^{\pm}$ as a
function of the collision energy $\sqrt{s}$ (GeV). Results from ALICE [51, 7,
24, 58, 63, 27, 26, 60] are compared with those from CMS [90] at the same
energies, and with those from STAR [88, 91, 92] at $\sqrt{s}=200\text{\,GeV}$
and PHENIX [20, 19] at $\sqrt{s}=62.4\text{\,GeV}$ obtained at BNL RHIC. The
CMS data points are slightly shifted horizontally for clarity. Open boxes
represent statistical and systematic uncertainties summed in quadrature. Red
data points are ALICE measurements for $\sqrt{s}=13\text{\,TeV}$.

It is worth noting that the proton and $\Lambda$ have different $\langle
p_{\mathrm{T}}\rangle$ values than the $\mathrm{K^{*0}}$ and $\phi$, despite
the similar masses of these particles. In particular, regardless of the
considered $\sqrt{s}$, the average transverse momenta are greater for those
resonances than for protons or $\Lambda$, with the $\phi$ having the largest
value. This clearly indicates a violation of mass ordering among these
particles. It is noteworthy that in our measurements at
$\sqrt{s}=7\text{\,TeV}$ as a function of charged-particle multiplicity, the
difference in $\langle p_{\mathrm{T}}\rangle$ between ${\rm
p}({\rm\overline{p}})$ and $\mathrm{K^{*0}}$ (and $\phi$) increases with
multiplicity [53].

It is of particular importance to study the $p_{\mathrm{T}}$-integrated
particle ratios as a function of collision energy, which might provide more
information on hadron production mechanisms and their dependence on the
collision energy or charged particle multiplicity. Given the observation of an
increase in (multi-)strange hadron production as a function of multiplicity in
pp collisions at $\sqrt{s}=7\text{\,TeV}$ [93], an increase in these yields
with collision energy would also be expected. Figure 9 shows the ratios of
particle yields to the yields of pions and kaons as a function of the
collision energy; results from ALICE are compared with those from CMS [90] and
with STAR [88, 91, 92] and PHENIX [20, 19] results from RHIC measured at
$\sqrt{s}=$ 200 and 62.4 GeV. All the ratios except $\Omega/\pi$ and $\Xi/\pi$
appear to saturate in the LHC energy regime. For the $\Omega/\pi$ and
$\Xi/\pi$ ratios, the relative increases are ${\sim}29\%$ and ${\sim}18\%$,
respectively; these are larger than the ${\sim}12\%$ increase of the
$\Lambda/\pi$ ratio, indicating that the strangeness content may control the
magnitude of the increase in the yield ratios as the energy changes from
$\sqrt{s}=$ 7 to 13 TeV.

The ratios of the $\phi$ and $\mathrm{K^{*0}}$ yields to those of charged
kaons ($\phi/\text{K}$ and $\mathrm{K^{*0}}/\text{K}$) do not exhibit any
dependence on the collision energy. It is worth noting that there may be a
decrease in the $\mathrm{K^{*0}}/\text{K}$ ratio in high-multiplicity pp
collisions at $\sqrt{s}=7\text{\,TeV}$ [53]. In nuclear collisions, such a
modification of resonance yields is often described as a consequence of
scattering processes during the hadron gas phase of the collision system
evolution.

###  6.3 Mass and baryon number effects on $\langle p_{\mathrm{T}}\rangle$

Figure 10: Average transverse momenta of light-flavor hadrons as a function of
hadron mass (left) and as a function of hadron mass normalized to the number
of constituent quarks (right) for different collision energies ($\sqrt{s}$).
Results reported from ALICE [51, 7, 24, 58, 63, 27, 26] are compared with
those measured by STAR at RHIC at $\sqrt{s}=200\text{\,GeV}$ [88, 91, 92].
Uncertainties that are not visible are smaller than the size of the symbol.
Solid and dashed lines are drawn as visual aids and represent separate linear
fits to meson and baryon data at $\sqrt{s}=13\text{\,TeV}$. The data points
for $\mathrm{K_{S}^{0}}$ are slightly shifted horizontally for clarity. Red
data points are for $\sqrt{s}=13\text{\,TeV}$.

The average transverse momenta $\langle p_{\mathrm{T}}\rangle$ as a function
of the particle mass are reported in the left panel of Fig. 10 for all light-
flavor hadrons under study. The $\sqrt{s}=13\text{\,TeV}$ pp results are
compared to our earlier measurements reported at lower collision energies [51,
24, 63, 58, 27, 7, 26], and to those measured at $\sqrt{s}=200\text{\,GeV}$ by
the STAR Collaboration at RHIC [88, 91, 92]. At $\sqrt{s}=13\text{\,TeV}$, two
different linear trends can be observed for mesons and baryons separately,
reflecting the violation of mass ordering in $\langle p_{\mathrm{T}}\rangle$.

The observation of scale breaking between mesons and baryons as a function of
the transverse mass, as discussed in Sec. 6.1.1, might lead to a violation of
mass ordering in the average transverse momenta of the produced particles. For
particles with similar masses (like $\mathrm{K^{*0}}$, p, $\phi$, and
$\Lambda$), the meson spectra will be harder and the $\langle
p_{\mathrm{T}}\rangle$ values greater in comparison to the baryons.
Furthermore, the softer pion $m_{\mathrm{T}}$ spectra could lead to a
deviation of the pion $\langle p_{\mathrm{T}}\rangle$ values from the trend
observed for the other mesons.

For $\sqrt{s}<1$ TeV all the considered hadrons appear to exhibit scaling of
$\langle p_{\mathrm{T}}\rangle$ with the reduced hadron mass $m/n_{q}$, i.e.
the mass normalized by the number of constituent quarks $n_{q}$. This is
observed in pp collisions at $\sqrt{s}=900\text{\,GeV}$ at ALICE and is shown
in the right panel of Fig. 10. Reference [94] suggests that such scaling holds
even for multi-strange baryons ($\Xi$, $\Omega$) when the average
pseudorapidity density of charged particles measured at mid-pseudorapidity is
small ($\left<\text{d}N_{\rm{ch}}/\text{d}\eta\right>=3.81\pm
0.01~{}(\text{stat.})\pm 0.07~{}(\text{syst.})$ [95]). At higher collision
energies, the scaling behavior is obviously broken. Two separate trends are
observed at $\sqrt{s}=13\text{\,TeV}$: one for mesons (slope $1.71\pm 0.06$)
and one for baryons (slope $2.42\pm 0.15$).

###  6.4 Yield ratios

The ratios of hadron yields are investigated as a function of transverse
momentum. This allows the $p_{\mathrm{T}}$ spectra of different particle
species, characterized by their unique mass and quark content, to be compared.
Additionally, measurements at different $\sqrt{s}$ are included, which helps
to quantify any change in spectral shapes with $\sqrt{s}$. The uncertainties
related to normalization cancel in these ratios. Figure 11 shows the ${\rm
p}/\pi$, $\Lambda/\mathrm{K_{S}^{0}}$ (left panel), $\Xi/\phi$, and
$\Omega/\phi$ (right panel) baryon-to-meson ratios as a function of
$p_{\mathrm{T}}$ at $\sqrt{s}=7$ [25, 66, 27] (open symbols) and
$\sqrt{s}=13\text{\,TeV}$ (full symbols). The left panel includes particle
ratios with baryons containing zero (${\rm p}/\pi$) and one
($\Lambda/\mathrm{K_{S}^{0}}$) strange valence quark, whereas the right panel
collects ratios for hadrons with two ($\Xi/\phi$) and three ($\Omega/\phi$)
strange quarks. The $\Omega/\phi$ ratio compares hadrons that consist entirely
of strange valence (anti)quarks: three for $\Omega$ in and two for $\phi$.

At low $p_{\mathrm{T}}$, all of the ratios increase with $p_{\mathrm{T}}$ as
expected from the higher $\langle p_{\mathrm{T}}\rangle$ observed for higher
mass particles. In this $p_{\mathrm{T}}$ regime, all of the ratios at
$\sqrt{s}=13\text{\,TeV}$ show good agreement with those at
$\sqrt{s}=7\text{\,TeV}$ within their systematic uncertainties, suggesting
that the collision energy has no observable effect on the magnitude or shape
of these yield ratios. This observation remains valid in the higher
$p_{\mathrm{T}}$ region. Within the systematic uncertainties the ${\rm
p}/\pi$, $\Lambda/\mathrm{K_{S}^{0}}$, $\Xi/\phi$, and $\Omega/\phi$ ratios
are consistent for the two collision energies. It is noteworthy that the ${\rm
p}/\pi$ ratio at $\sqrt{s}=13\text{\,TeV}$ has a hint of enhancement at
intermediate $p_{\mathrm{T}}$ with respect to that at
$\sqrt{s}=7\text{\,TeV}$, however this is barely significant given the quoted
uncertainties. From $p_{\mathrm{T}}>10\text{\,GeV/$c$}$ onward, the ${\rm
p}/\pi$ ratio becomes fairly constant for both collision energies.

Figure 11: Baryon-to-meson ratios as a function of transverse momentum
($p_{\mathrm{T}}$) measured in pp collisions at $\sqrt{s}=$ 13 (full symbols)
and 7 TeV [25, 27] (open symbols). Note the different scale of the horizontal
axis in the right panel.

###  6.5 Comparison to theoretical models

The high-precision measurements of the $p_{\mathrm{T}}$ spectra reported in
this paper are crucial inputs for the further tuning of Monte Carlo event
generators and to improve the understanding of particle production mechanisms
at the investigated collision energies. The measurements of the light-flavor
hadron species presented above are compared with Monte Carlo model predictions
based on general-purpose event generators: PYTHIA 6, PYTHIA 8, and EPOS-LHC.

The PYTHIA MC model contains a rigorous description of hard scatterings
through pQCD, combined with phenomenological models for semi-hard/soft
processes. The Perugia-2011 tune [96] for PYTHIA 6 and the Monash 2013 tune
[31] for PYTHIA 8, which have different sets of parameters, are used. The
PYTHIA 8 model uses an updated parameter set for Lund hadronization for light
(and heavy) quarks. The widely-used Monash 2013 tune has an improved
description of diffractive processes with respect to PYTHIA 6. Both sets of
parameters have been obtained from recent (2011 and 2013) analyses of minimum
bias, underlying event, and/or Drell–Yan measurements in pp collisions at
$\sqrt{s}=7\text{\,TeV}$; the Monash 2013 tune was optimized to describe early
data collected by the LHC experiments as well as lower energy data. Moreover,
both versions of the model have strong final-state parton interactions
implemented through different color reconnection models [97, 98]. As a
consequence of the different energy evolution of the cutoff for multiple
parton interactions, the Monash tune has larger multiple parton interaction
activity at a given collision energy than the Perugia tune.

Conversely, the EPOS-LHC event generator (used with the CRMC package, version
1.5.4) invokes Gribov’s Reggeon Field Theory [99] for multiple scatterings;
this formalism features collective hadronization with the core-corona
mechanism [100]. After multiple scattering, the final-state partonic system
consists mainly of longitudinal flux tubes that fragment into string segments.
If the energy density from string segments is high enough, they fuse into the
so-called “core” region, which then evolves hydrodynamically and eventually
hadronizes to form the bulk part of the system. On the other hand, in the low-
density region the strings expand and eventually break via the production of
quark-antiquark pairs, which hadronize using the unmodified string
fragmentation and form the “corona” region. The EPOS-LHC model uses recent
data available from LHC, which helps in reproducing minimum bias results with
transverse momenta up to a few GeV/$c$.

Figure LABEL:fig:models shows the ratios of $p_{\mathrm{T}}$ spectra extracted
from the PYTHIA 8, PYTHIA 6, and EPOS-LHC models to the spectra measured at
$\sqrt{s}=$ 7 and 13 TeV; the models give similar descriptions of the data at
both energies. For the comparison of the measured data to the MC generators,
the total fractional uncertainties of the data are shown, i.e. the statistical
and systematic uncertainties of the measurement have been summed in
quadrature. Both PYTHIA versions can generally describe the shapes of particle
$p_{\mathrm{T}}$ spectra with reasonable accuracy at intermediate
$p_{\mathrm{T}}$ ($2\text{\,GeV/$c$}\lesssim p_{\mathrm{T}}\lesssim
10\text{\,GeV/$c$}$), but generally give softer $p_{\mathrm{T}}$ spectra than
observed at low $p_{\mathrm{T}}$ ($\lesssim 2\text{\,GeV/$c$}$), although
pions show the opposite trend, with a harder spectrum. This
low-$p_{\mathrm{T}}$ behavior is due to the fact that the MC generators are
known to have difficulties at describing diffractive processes that play a
role at very low $p_{\mathrm{T}}$ [46]. In contrast, at high $p_{\mathrm{T}}$
PYTHIA predicts harder spectra.

For the most abundant pions and charged kaons, the spectral shapes above
${\sim}1\text{\,GeV/$c$}$ are described significantly better by PYTHIA 8 than
PYTHIA 6. For both $\mathrm{K}^{\pm}$ and $\mathrm{K_{S}^{0}}$, the model-to-
measured ratios for both PYTHIA versions are qualitatively similar, although
PYTHIA 8 provides a better agreement with the data at mid-to-high
$p_{\mathrm{T}}$ than PYTHIA 6. These PYTHIA tunes describe the (anti)proton
$p_{\mathrm{T}}$ spectra best for $0.8<p_{\mathrm{T}}<4\text{\,GeV/$c$}$. For
single-strange ($\Lambda$) and multi-strange ($\Xi$, $\Omega$) baryons, both
of the PYTHIA tunes underestimate the amount of produced particles in almost
the entire $p_{\mathrm{T}}$ range. A similar discrepancy is observed for the
case of PYTHIA 6 with the Perugia 2011 tune at $\sqrt{s}=7\text{\,TeV}$ [27]
and with older tunes (D6T, ATLAS-CSC, Perugia 0) for the $\Xi$ as reported for
$\sqrt{s}=900\text{\,GeV}$ in Ref. [58]. It is worth noting that a discrepancy
is also observed when the yields of $\Lambda$, $\Xi$, and $\Omega$ measured as
a function of the charged-particle multiplicity [28] are compared to results
from PYTHIA and EPOS-LHC. In all cases, the observed deviation becomes larger
for hadrons with higher strange valence quark content. In Fig.
LABEL:fig:models, the model underestimates the measured $\Xi$ $(\Omega)$
yields by a factor of two (four), although the discrepancy weakens for the
$\Xi$ towards higher $p_{\mathrm{T}}$ values. For the $\phi$ meson (with zero
net strangeness content), PYTHIA predicts yields within $10-20\%$ for
$p_{\mathrm{T}}{\sim}8\text{\,GeV/$c$}$, while at higher $p_{\mathrm{T}}$ the
measured yields and the PYTHIA values agree within uncertainties.

For $\pi^{\pm}$, $\mathrm{K}^{\pm}$, and ${\rm p}({\rm\overline{p}})$, the
EPOS-LHC model predicts the spectral shape and normalization better than both
of the PYTHIA tunes in the entire $p_{\mathrm{T}}$ range. For the resonances
and multi-strange baryons, EPOS-LHC gives harder $p_{\mathrm{T}}$ spectra than
the measured ones and performs better than PYTHIA at describing the yields of
the multi-strange baryons. The EPOS-LHC model gives an accurate description
for the $\Xi$ baryon at $p_{\mathrm{T}}\lesssim 2\text{\,GeV/$c$}$ both in
shape and normalization, but deviates at higher $p_{\mathrm{T}}$. The model
also describes the shape of the $p_{\mathrm{T}}$ spectrum of $\Omega$ baryons,
which have higher strangeness content, but does not reproduce the yield. For
the $\mathrm{K_{S}^{0}}$ and $\Lambda$ (reported here for the first time at
both $\sqrt{s}=$ 7 and 13 TeV), as well as $\mathrm{K^{*0}}$ and $\phi$, the
models give similar descriptions at both collision energies. However, at
$\sqrt{s}=7\text{\,TeV}$ the agreement in normalization worsens for
$\mathrm{K_{S}^{0}}$ at intermediate $p_{\mathrm{T}}$, while EPOS-LHC predicts
harder spectra at higher $p_{\mathrm{T}}$. For the $\phi$ meson, the EPOS-LHC
model shows a marginal agreement with the data at low $p_{\mathrm{T}}$, and
monotonically deviates from the measured spectrum as the $p_{\mathrm{T}}$
increases.

Generally, the deviations of these models from the ALICE measurements are
similar to those observed at lower $\sqrt{s}$, which were reported in Refs.
[63, 26] for $\mathrm{K^{*0}}$ and $\phi$ and in Ref. [24] for $\pi^{\pm}$,
$\mathrm{K}^{\pm}$, and ${\rm p}({\rm\overline{p}})$, although with a
restricted $p_{\mathrm{T}}$ reach for the latter three particle species. To
study how the models follow the changes in spectral shapes and normalization
as a function of $p_{\mathrm{T}}$ as the energy increases from $\sqrt{s}=$ 7
to 13 TeV, Fig. 12 shows the double ratios: the ratios of the measured
$p_{\mathrm{T}}$ spectra at $\sqrt{s}=13\text{\,TeV}$ to
$\sqrt{s}=7\text{\,TeV}$ (see Fig. 4) are divided by the same ratios obtained
from MC models. The comparisons indicate that the models capture well the
increase of the yields with $\sqrt{s}$ for the mesons shown in panels (a) –
(e); for $\pi^{\pm}$, $\mathrm{K}^{\pm}$, and $\mathrm{K_{S}^{0}}$, such an
observation between two distinct collision energies in this $p_{\mathrm{T}}$
regime is reported here for the first time. For charged kaons, the models
predict a more pronounced hardening with collision energy at high
$p_{\mathrm{T}}$. The tension seen in the range $2\lesssim
p_{\mathrm{T}}\lesssim 6\text{\,GeV/$c$}$ originates from the different
analysis techniques which were used to obtain the combined spectra for
$\sqrt{s}=7\text{\,TeV}$ and at $\sqrt{s}=13\text{\,TeV}$. Note that the
normalization uncertainties at both energies are not included in the reported
fractional uncertainties in the figure. For both $\mathrm{K^{*0}}$ and $\phi$
mesons, the similarity of the deviations of the PYTHIA model (though with a
different tune) from ALICE measurements was seen in our earlier measurements
at $\sqrt{s}=2.76\text{\,TeV}$ and $\sqrt{s}=7\text{\,TeV}$, as reported in
Ref. [63]. Panels (f) and (g) demonstrate that the PYTHIA 6 and EPOS-LHC
models agree with the measurements of ${\rm p}({\rm\overline{p}})$ and
$\Lambda$, in contrast, PYTHIA 8 systematically underestimates the measured
data for $p_{\mathrm{T}}\lesssim 8\text{\,GeV/$c$}$. For the multi-strange
baryons (panels (h) and (i)), EPOS-LHC predicts the $\sqrt{s}$ evolution of
the $p_{\mathrm{T}}$ spectra above $p_{\mathrm{T}}=1\text{\,GeV/$c$}$ the
best, however the PYTHIA models also agree with the measured data for the
$\Omega$ baryons above $p_{\mathrm{T}}{\sim}3\text{\,GeV/$c$}$. For baryons,
the model description for PYTHIA 8 improves as the strangeness content
increases.

Figure 12: The measured 13-to-7 TeV ratios of hadron $p_{\mathrm{T}}$ spectra
are compared with model calculations: where $Y$ is the particle yield, the
plotted double ratio is
$(Y^{13\;\mathrm{TeV}}_{\mathrm{model}}/Y^{7\;\mathrm{TeV}}_{\mathrm{model}})/(Y^{13\;\mathrm{TeV}}_{\mathrm{measured}}/Y^{7\;\mathrm{TeV}}_{\mathrm{measured}})$.
The measured $p_{\mathrm{T}}$ spectra at $\sqrt{s}=7\text{\,TeV}$ are from
Refs. [25, 27, 60]. The total fractional uncertainties of the measured yield
ratios are shown as shaded boxes.  Figure 13: Baryon-to-meson particle ratios
as a function of transverse momentum measured in pp collisions at
$\sqrt{s}=13\text{\,TeV}$. The data are compared with several general-purpose
Monte Carlo models. Statistical and systematic uncertainties are summed in
quadrature, and are shown as vertical error bars.

The relative impact of the hardening observed in the $p_{\mathrm{T}}$ spectra
in Fig. 4 for the various light-flavor hadrons is better seen in terms of the
$p_{\mathrm{T}}$-dependent yield ratios of different particle species. The
comparison of these ratios to MC models allows the different hadronization
mechanisms implemented in the event generators to be tested. In Fig. 13 the
$p_{\mathrm{T}}$-dependent baryon-to-meson ratios measured in pp collisions at
$\sqrt{s}=13\text{\,TeV}$ and discussed in the previous section are shown as a
function of $p_{\mathrm{T}}$ and compared with the same MC models discussed
above. The ratios for multi-strange baryons are better approximated by PYTHIA
6 than PYTHIA 8. As discussed for multi-strange baryon measurements performed
by ALICE in pp collisions at $\sqrt{s}=7\text{\,TeV}$ [27], this is a
consequence of the removal of the “popcorn” mechanism [101] in the Perugia
2011 tune which in turn suppresses baryon production by favoring soft quark-
antiquark pairing.

Although the Perugia 2011 tune reproduces the shape of the $\Omega$ baryon
$p_{\mathrm{T}}$ spectrum within roughly $10\%$ (as seen in Fig.
LABEL:fig:models), the $\Omega/\phi$ ratio given by this model has large
deviations from the data for $p_{\mathrm{T}}\lesssim 4\text{\,GeV/$c$}$, the
same $p_{\mathrm{T}}$ region where the model $\phi$ meson spectrum has its
largest disagreement with the measured data (about $40\%$). The Perugia 2011
tune better characterizes the $p_{\mathrm{T}}$ evolution of the $\Xi/\phi$
ratio than the Monash 2013 tune of PYTHIA 8 for
$p_{\mathrm{T}}>2\text{\,GeV/$c$}$. The EPOS-LHC model gives a good
description of the shapes of the multi-strange baryon-to-meson ratios in the
entire $p_{\mathrm{T}}$ range of the measurements. The model predicts
quantitatively the magnitude of the $\Xi/\phi$ ($\Omega/\phi$) ratio with
reasonable accuracy above (below) $p_{\mathrm{T}}=2\text{\,GeV/$c$}$. All
three models approximate the $p_{\mathrm{T}}$ evolution of the
$\Lambda/\mathrm{K_{S}^{0}}$ ratio qualitatively, but systematically over-
predict its magnitude for the full $p_{\mathrm{T}}$ range. The PYTHIA 8 and
EPOS-LHC models predict the maximum of the ratio at higher $p_{\mathrm{T}}$
values than measured. EPOS-LHC fails to reproduce the magnitude of the
$\Lambda/\mathrm{K_{S}^{0}}$ ratio, deviating from the measured ratio by about
a factor of three at the maximum of the peak. The ${\rm p}/\pi$ ratio is well
described by the models in the low-$p_{\mathrm{T}}$ region ($\lesssim
2\text{\,GeV/$c$}$), however they fail to follow the enhancement and depletion
behavior seen in the $p_{\mathrm{T}}$ region between $2\text{\,GeV/$c$}$ and
$6\text{\,GeV/$c$}$. All three models indicate a flattening behavior above
${\sim}10\text{\,GeV/$c$}$, similar to what is seen in the measured data.

###  6.6 Comparison to pQCD calculations

The measured invariant cross sections are compared to next-to-leading order
(NLO) perturbative QCD calculations using CT10NLO proton PDFs [102] with the
DSS (de Florian, Sassot, and Stratmann) FF set [103, 3]. For charged pions a
new version of the DSS FFs is available: the DSS14 FF [5] set. The NLO
calculations are based on Ref. [104] which applies the same factorization
scale value, $\mu=p_{\mathrm{T}}$ for the factorization, renormalization and
fragmentation scales. The variation of the scales to $\mu=p_{\mathrm{T}}/2$
and $\mu=2\,p_{\mathrm{T}}$ gives an estimate of the theoretical uncertainty;
the PDF uncertainties are negligible in comparison to the scale uncertainty.
The rather large scale uncertainty observed at lower $p_{\mathrm{T}}$
($2<p_{\mathrm{T}}<10\text{\,GeV/$c$}$) stabilizes at $\pm 20-30\%$ for
$p_{\mathrm{T}}\simeq 10\text{\,GeV/$c$}$, which is the region where the NLO
calculations are trustworthy and free from non-perturbative effects.

The production of $\pi^{\pm}$, $\mathrm{K}^{\pm}$, and ${\rm
p}({\rm\overline{p}})$ from hard scattering becomes dominated by gluon
fragmentation with increasing collision energy in the $p_{\mathrm{T}}$ range
of the measurement [3]. The presented identified charged-hadron spectra can
therefore help to constrain the gluon-to-charged-hadron fragmentation function
[4] which is of crucial importance to a better description of the LHC charged-
hadron data with NLO pQCD (see e.g. Ref. [105]). The invariant differential
cross sections for $\pi^{\pm}$, $\mathrm{K}^{\pm}$, and ${\rm
p}({\rm\overline{p}})$ are shown in the left panel of Fig. 14 in comparison to
NLO pQCD calculations. Since, to date, no calculation exists for the scale
uncertainties of the DSS14, only the DSS FFs with the corresponding scale
uncertainties are reported in the figure. In the right panel of Fig. 14 the
ratios of the measured data and the NLO pQCD calculations to the Lévy–Tsallis
fits of the $\pi^{\pm}$, $\mathrm{K}^{\pm}$, and ${\rm p}({\rm\overline{p}})$
cross sections are shown. For $p_{\mathrm{T}}>10\text{\,GeV/$c$}$ the NLO pQCD
calculations, employing the DSS14 FFs and using the DSS scale uncertainties,
over-predict the measured pion cross section by up to a factor of
approximately two, but describe the shape of the $p_{\mathrm{T}}$ spectrum
rather well. Similar discrepancies between NLO pQCD calculations and the
measured cross sections have also been reported for the measurements of
neutral pions ($\pi^{0}$) at $\sqrt{s}~{}=~{}7$ and 8 TeV [106, 10] from
ALICE, leaving room for future improvements in the calculations. It is worth
noting that the published $\pi^{0}$ measurement at $\sqrt{s}=7\text{\,TeV}$
[106] adds important constraints for gluon FFs, which would help reduce the FF
uncertainties. The NLO calculations describe charged kaons better than pions,
which is reflected in the better agreement between the calculated cross
section and with the measured data-to-fit ratio within the quoted
uncertainties. The deviations between the NLO calculations and the data
generally increase with $p_{\mathrm{T}}$. This increase is significantly
stronger for protons, which deviate the most from the measured values for all
$\mu$ scale choices. The NLO calculations significantly overestimate the
measured data at high $p_{\mathrm{T}}$.

Figure 14: (Left) Invariant differential cross sections for $\pi^{\pm}$,
$\mathrm{K}^{\pm}$, and ${\rm p}({\rm\overline{p}})$ production, compared with
NLO pQCD calculations using CT10NLO PDFs [102] with DSS14 FFs [5] for
$\pi^{\pm}$ and DSS FFs [103, 3] for $\mathrm{K}^{\pm}$ and ${\rm
p}({\rm\overline{p}})$. The shaded bands around the NLO calculations
correspond to the scale uncertainty. (Right) Ratio of measured data (points)
or NLO calculations (lines) to a Lévy–Tsallis function that is fitted to the
data. The NLO-to-fit ratio is shown for the scale $\mu=p_{\mathrm{T}}$ and the
variations $\mu=p_{\mathrm{T}}/2$ and $\mu=2\,p_{\mathrm{T}}$. The fully
correlated normalization uncertainty is indicated in the legend.  Figure 15:
Ratios of transverse-momentum spectra of $\pi^{\pm}$, $\mathrm{K}^{\pm}$, and
${\rm p}({\rm\overline{p}})$ in minimum bias inelastic pp collisions at
$\sqrt{s}=13\text{\,TeV}$ to those at $\sqrt{s}=7\text{\,TeV}$ [25]. Solid
lines show predictions from a next-to-leading order (NLO) pQCD calculation
using DSS fragmentation functions [104]. Statistical and systematic
uncertainties are shown as vertical error bars and boxes, respectively.

At lower $p_{\mathrm{T}}$ ($<10\text{\,GeV/$c$}$) the discrepancy between the
NLO calculations and the measured data is reduced for all particle species,
but at the cost of an increase of the scale uncertainties. In this
$p_{\mathrm{T}}$ regime, soft parton interactions and resonance decays
dominate particle production, which cannot be described within the framework
of pQCD. This is also reflected in the shape of the distributions, which is
not described by the calculation. The presented results show that independent
fragmentation works rather well for mesons, in particular for
$p_{\mathrm{T}}>6\text{\,GeV/$c$}$, however for (anti)protons the spectral
shape deviates towards the high-$p_{\mathrm{T}}$ region.

Though the $p_{\mathrm{T}}$ dependence of the cross sections at a given
collision energy is not described well by the NLO calculations, the
$p_{\mathrm{T}}$ dependence of the ratios of the $p_{\mathrm{T}}$ spectra at
$\sqrt{s}=$ 7 [25] and 13 TeV is described better. Such a comparison is shown
in Fig. 15, where the ratio of invariant yields in inelastic pp collisions at
$\sqrt{s}=$ 13 TeV to those at 7 TeV is compared to the same ratio calculated
using NLO pQCD. The agreement between the data and NLO calculations is notably
improved compared to the NLO results for the spectra themselves. Taking a
double ratio, where the ratios of the measured spectra at $\sqrt{s}=$ 13 and 7
TeV are divided by the ratios from the NLO pQCD calculations, the observed
difference is below $10\%$ ($20\%$) for pions (kaons and protons).

##  7 Summary

The production of light-flavor hadrons at midrapidity was measured in
inelastic pp collisions at $\sqrt{s}=13\text{\,TeV}$. Additionally, single-
particle $p_{\mathrm{T}}$ distributions of $\mathrm{K_{S}^{0}}$, $\Lambda$,
and $\overline{\Lambda}$ were also measured in inelastic pp collisions at
$\sqrt{s}=7\text{\,TeV}$. The presented measurements complement the existing
ones at lower collision energies, allowing particle production to be studied
over a wide range of $\sqrt{s}$. All $p_{\mathrm{T}}$ distributions are
observed to become harder for $p_{\mathrm{T}}>2\text{\,GeV/$c$}$ with
increasing collision energy.

The $p_{\mathrm{T}}$-integrated hadron yields normalized to the pion yields
saturate as a function of $\sqrt{s}$ at LHC energies for the studied non-
strange and single-strange hadrons. In contrast, a hint of an increase in the
multi-strange hadron yields is apparent as $\sqrt{s}$ increases from
$\sqrt{s}=$ 7 to 13 TeV, the increase being more pronounced for $\Omega$
baryons, which have the largest strange quark content among the studied
hadrons. This observation is similar to our recent measurements performed as a
function of charged-particle multiplicity in pp collisions at $\sqrt{s}=$ 7
and 13 TeV [53, 28].

A significant deviation from the empirical transverse mass scaling relation of
the production cross sections between $\pi^{\pm}$ and $\mathrm{K}^{\pm}$ and
$\mathrm{K_{S}^{0}}$ is observed for $p_{\mathrm{T}}<2\text{\,GeV/$c$}$.
Empirical $x_{\mathrm{T}}$ scaling for $\pi^{\pm}$, $\mathrm{K}^{\pm}$,
$\mathrm{K^{*0}}$, and ${\rm p}({\rm\overline{p}})$ is well satisfied (within
roughly 20–40%) in the hard scattering region of particle production.