Multiplicity dependence of π , K, and p production in pp collisions at s √ =13 TeV
EEUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH
CERN-EP-2020-02403 March 2020c (cid:13)
Multiplicity dependence of π , K, and p production in pp collisions at √ s = TeV
ALICE Collaboration ∗ Abstract
This paper presents the measurements of π ± , K ± , p and p transverse momentum ( p T ) spectra as afunction of charged-particle multiplicity density in proton-proton (pp) collisions at √ s =
13 TeVwith the ALICE detector at the LHC. Such study allows us to isolate the center-of-mass energy de-pendence of light-flavour particle production. The measurements reported here cover a p T range from0.1 GeV / c to 20 GeV / c and are done in the rapidity interval | y | < .
5. The p T -differential particleratios exhibit an evolution with multiplicity, similar to that observed in pp collisions at √ s = p T -integrated hadron-to-pion yield ratios measured in ppcollisions at two different center-of-mass energies are consistent when compared at similar multi-plicities. This also extends to strange and multi-strange hadrons, suggesting that, at LHC energies,particle hadrochemistry scales with particle multiplicity the same way under different collision ener-gies and colliding systems. ∗ See Appendix A for the list of collaboration members a r X i v : . [ nu c l - e x ] M a r LICE Collaboration
The unprecedented energies available at the Large Hadron Collider (LHC) provide unique opportuni-ties to investigate the properties of strongly-interacting matter. Particle production at large transversemomenta ( p T ) is well-described by perturbative Quantum Chromodynamics (pQCD). The soft regime( p T (cid:46) c ), in which several collective phenomena are observed in proton-proton (pp), proton-lead(p–Pb), and heavy-ion (A–A) collisions, is not calculable from first principles of QCD. Instead, in orderto describe bulk particle production in A–A collisions, one usually relies on hydrodynamic and ther-modynamic modelling, which assumes the system to be in kinetic and chemical equilibrium [1, 2]. Onthe other hand, the description of low- p T particle spectra in smaller systems such as pp collisions is of-ten based on phenomenological modelling of multi-partonic interactions (MPI) and color reconnection(CR) [3, 4] or overlapping strings [5].Recent reports on the enhancement of (multi-)strange hadrons [6], double-ridge structure [7, 8], non-zero v coefficients [9], mass ordering in hadron p T spectra, and characteristic modifications of baryon-to-meson ratios [10] suggest that collective phenomena are present at the LHC energies also in p–Pbcollisions. This is further extended to even smaller systems, such as pp collisions at √ s = N ch / d η from pp to p–Pband then to Pb–Pb collisions was found [11–13]. The observed similarities suggest the existence of acommon underlying mechanism determining the chemical composition of particles produced in thesethree collision systems.Results from pp [11] and p–Pb [10] collisions indicate that particle production scales with d N ch / d η de-spite the colliding system. However, measurements reported in previous multiplicity-dependent studieshave considered different colliding systems, each at a different center-of-mass energy. In this work, weextend the existing observations by performing a detailed study of pp collisions at √ s =
13 TeV. Thanksto the availability of Run 2 data from the LHC, for the first time, in pp collisions, we can disentangle theeffect of center-of-mass energy from the multiplicity dependence of π ± , K ± and p (p) production.In this paper, we report on the multiplicity dependence of the production of primary π ± , K ± and p (p)at √ s =
13 TeV. Particles are considered as primary if they are created in the collision (includingproducts of strong and electromagnetic decays), but not from a weak decay of other light-flavor hadronsor muons [14]. The reported particle spectra are measured in the rapidity region | y | < . p T = . / c to several tens of GeV / c [16]. As particles and anti-particles are produced in equal amounts at LHCenergies [17], we adopt a notation where π , K, and p refer to ( π + + π − ) , ( K + + K − ) , and ( p + p) unlessstated otherwise. This paper is organized as follows. In Sec. 2, the details on particle identificationtechniques, systematic uncertainties, spectra corrections and normalization are provided. The results arepresented and discussed in Sec. 3, together with comparisons to Monte Carlo model predictions. Finally,the most important findings are summarized in Sec. 4. The dataset used for this study was recorded by the ALICE Experiment during the 2016 LHC pp datataking period. Overall ∼ .
47 nb − considering the visible cross-section measured with the V0 detector [18]. A detailed descrip-tion of the ALICE detector and its performance is provided in [15, 16]. Measurements of identified par-ticle spectra have been performed by using the central barrel detectors: the Inner Tracking System (ITS)(Sec. 3.1 of [15]), the Time Projection Chamber (TPC) [19] and the Time-of-Flight detector (TOF) [20].The charged-particle multiplicity estimation is done by the V0 detector (Sec. 5.4 of [15]), which consists2LICE Collaborationof two arrays of 32 scintillators each, positioned in the forward (V0A, 2 . < η < .
1) and backward(V0C, − . < η < − .
7) rapidity regions. In addition, the V0 is also used for triggering purposes aswell as background rejection. The determination of the event collision time [21] is performed by the T0detector as well as the TOF detector. The former consists of two arrays of Cherenkov counters, posi-tioned on both sides of the interaction region, and covering a pseudorapidity range of − . < η < − . . < η < | η | < .
9. The twoinnermost layers form the Silicon Pixel Detector (SPD), which features binary readout and is also usedas a trigger detector. The Silicon Drift Detector (SDD) and the Silicon Strip Detector (SSD), which formthe four outer layers of the ITS, provide the amplitude of the charge signal, which is used for particleidentification through the measurement of specific energy loss at low transverse momenta ( p T (cid:38) < r <
247 cm and −
250 cm < z <
250 cm, respectively. The TPC is read out by multi-wire proportional chambers (MWPC) with cathodepad readout, located at its endplates. With the measurement of drift time, the TPC provides three-dimensional space-point information for each charged track in pseudorapidity range | η | < . .
2% [16].The TOF is a large-area array of multigap resistive plate chambers (MRPC), formed into a ∼ | η | < . . (cid:46) p T (cid:46) / c .The arrival time is measured by the MRPCs with an intrinsic time resolution of 50 ps, while the eventcollision time is determined by combining the T0 detector measurement with the estimate using theparticle arrival times at the TOF [21]. The analysed data were recorded using a minimum-bias trigger requiring signals in both V0A and V0Cscintillators in coincidence with the arrival of the proton bunches from both directions. The backgroundevents produced outside the interaction region are rejected using the correlation between the SPD clus-ters and the tracklets reconstructed in SPD. The out-of-bunch pileup was rejected offline using the timinginformation from the V0 counter. The primary vertex was reconstructed either using global tracks (re-constructed using ITS and TPC information) or SPD tracklets (reconstructed using only the SPD infor-mation) with | z vtx | <
10 cm along the beam axis. Events with in-bunch pileup were removed if a secondvertex was reconstructed within 8 mm of the primary vertex in the beam direction. The typical interactionrate of pp collisions in the 2016 data taking periods was around 120 kHz while beam-gas interactionsoccurred at a rate of 1.2 kHz.In the analysis presented in this paper, we consider the event class INEL > | η | <
1, which corresponds to ∼
75% of the total inelasticscattering cross-section [22]. To avoid auto-correlation biases [11, 22], the events are classified using thetotal charge collected in the V0 detector (V0M amplitude), which scales linearly with the total numberof the corresponding charged particles in its acceptance. For each event class, the corresponding meancharged-particle multiplicity density (cid:104) d N ch / d η (cid:105) is measured at mid-rapidity ( | η | < . ) as summarisedin Table 1. 3LICE Collaboration Table 1:
Mean charged-particle multiplicity density (cid:104) d N ch / d η (cid:105) measured in different event multiplicity classes. V0M mult. class I II III IV V σ / σ INEL > (%) 0–0.92 0.92–4.6 4.6–9.2 9.2–13.8 13.8–18.4 (cid:104) d N ch / d η (cid:105) . ± .
35 20 . ± .
27 16 . ± .
22 13 . ± .
19 12 . ± . σ / σ INEL > (%) 18.4–27.6 27.6–36.8 36.8–46.0 46.0–64.5 64.5–100 (cid:104) d N ch / d η (cid:105) . ± .
14 7 . ± .
11 6 . ± .
09 4 . ± .
07 2 . ± . Table 2:
Different p T ranges used for the identification of pions, kaons and protons. The final p T spectra havebeen obtained by combining the results of the various PID techniques. Analysis PID Technique p T ranges (GeV / c ) (pseudo)rapidity π ± K ± p (p) rangeITSsa n σ integral 0.1 − − − | y | < σ fits to TPC, β fits to TOF 0.25 − − − | y | < | η | < ∆ t − − − | y | < − − − | y | < E /d x fits 2 −
20 3 −
20 3 − | η | < In order to measure particle spectra in a wide p T range, several sub-analyses employing different detec-tors and particle identification (PID) techniques were performed and combined. As a result, the combinedspectra cover transverse momenta ranges from 0.1/0.2/0.3 GeV / c to 20 GeV / c for π /K/p. The p T and(pseudo)rapidity ranges covered by each analysis for different particle species are summarized in Table 2.At low p T , hadron spectra were measured by the ITS stand-alone (ITSsa) analysis. The dynamic rangeof the analogue readout of SDD and SSD allows for d E / d x measurements of highly ionizing particles,which otherwise do not reach the outer detectors. Hadron identification in the ITS is carried out bycalculating the truncated mean of d E / d x and comparing it to the expected energy loss under differentmass hypotheses. The difference between measured and expected d E / d x is then estimated in terms of thestandard deviation σ and the particle mass hypothesis with the lowest score is assigned. This is feasibleeven for pp collisions with the highest multiplicities, as the number of charge clusters wrongly assignedto the reconstructed tracks is negligible. A detailed description of the method is provided in [11].Hadrons at intermediate p T enter the fiducial volume of the TPC where they can be identified by measur-ing the charge generated in the gas. The truncated mean of d E / d x is calculated for the global tracks andcompared to the expected energy loss under a given mass hypothesis. At low transverse momenta wherethe separation between different species is sufficiently large, tracks within three standard deviations fromthe expected d E / d x are assigned to a given hypothesis. In the regions where signals from several speciesoverlap ( p T < . / c for π , p T > .
45 GeV / c for K, and p T > . / c for p), statistical unfoldingtechniques are used to remove the background contributions [11].4LICE CollaborationIn the p T region where the statistical unfolding of the TPC signal becomes unfeasible, particle iden-tification is performed using the time-of-flight measurements. The results presented in this paper wereobtained by combining the particle spectra estimated with two separate TOF analyses, taking into accountthe non-common part of the respective systematic uncertainties. In the “TOF template fits”, the PID isbased on a statistical unfolding method, where the distribution of the difference between measured andexpected time-of-flight (i.e. ∆ t ) is fitted with templates for pions, kaons and protons in each p T and mul-tiplicity bin [23]. An additional template is needed to take into account the background due to wronglyassociated tracks with hits in the TOF detector. The template for each particle is built from data, consid-ering the measured TOF time response function (Gaussian with an additional exponential tail for largerarrival times). The fits are repeated separately for each particle hypothesis in | y | < .
5. In contrast to this,in the “TOF fits” analysis, the velocity β distribution is simultaneously fitted for all three particle types.For this purpose, four analytic functions, three for π , K and p, and one for mismatches, are employed.The analysis is performed in two narrow pseudorapidity slices ( | η | < . . < | η | < .
4) and inmomentum bins, which are then unfolded to transverse momenta. The corresponding rapidity intervalis determined under the assumption of a flat d N ch / d η distribution in the aforementioned pseudorapiditybins [24].Charged kaons can also be identified via the kink decay topology, where a charged particle decays intoa charged and a neutral daughter ( K ± → µ ± ν µ or K ± → π ± π ). This secondary vertex where bothdecaying particle and the charged decay product have the same charge is reconstructed inside the ALICETPC detector. This technique extends the charged kaon identification up to 6 GeV/ c on a track-by-trackbasis. The algorithm for selecting kaons via their kink decay is used in a fiducial volume inside the TPCcorresponding to a radial distance of 120 < R <
210 cm. This selection allows for an adequate numberof TPC clusters to be associated with the decaying particle and its products. The track of the decayingparticle is required to fulfil all the criteria of the global tracks except for the minimum number of clusters,which in this case is 30.The topological selection of the kaon candidates and their separation from the pion decays ( π ± → µ ± ν µ )is based on the two-body decay kinematics. The transverse momentum of the decay product with respectto the decaying particle’s direction ( q T ) has an upper limit of 236 MeV/ c for kaons and 30 MeV/ c forpions for the two-body decay to µ ± ν µ . Similarly, for kaons decaying to pions, this limit is 205 MeV/ c .Thus, a selection of q T <
120 MeV / c rejects the majority (85%) of pion decays. In addition, the anglebetween the mother and the daughter tracks is selected to be above the maximum allowed decay anglefor pions and below the maximum allowed decay angle for kaons [25]. The invariant mass for the decay µ ± ν µ , M µν is calculated by assuming the daughter track to be a muon and the undetected track to be aneutrino. These selection criteria lead to a kaon sample with a purity of 97%.The strategy employed to measure particle production in the region of the relativistic rise of the TPC wasreported in [26]. The d E / d x signal in the relativistic rise ( < β γ (cid:0) = pm (cid:1) < ) follows the functionalform ln ( β γ ) . In addition to the logarithmic growth, the separation in number of standard deviationsbetween pions and protons, pions and kaons, and kaons and protons as a function of momentum isnearly constant, which allows identification of charged pions, kaons, and (anti)protons with a statisticaldeconvolution approach from p T ≈ − / c up to p T =
20 GeV / c . In order to describe the TPCresponse in the relativistic rise, clean external samples of secondary particles were used to parametrizethe Bethe-Bloch and resolution curves. These correspond to pions (protons) from weak decays: K S → π + + π − ( Λ → p + π − ) and electrons from photon conversion. Moreover, primary pions measured withthe TOF detector were used. The parametrization is done in four different pseudorapidity intervals. Forshort (long) tracks, i.e tracks within | η | < . ( . < | η | < . ) , the resolution for protons is ≈ . ( ≈ . ) , while for pions it is ≈ . ( ≈ . ) . To extract the fraction of charged pions, kaons, and protonsin the four different pseudorapidity intervals a 4-Gaussian fit (three for π , K, p and one to remove theunwanted electron contribution) to the d E / d x distribution in momentum bins is performed. The only5LICE Collaborationfree parameter in each of the Gaussian functions is the normalization, while the (cid:104) d E / d x (cid:105) and σ (cid:104) d E / d x (cid:105) are obtained and fixed using the Bethe-Bloch and resolution parametrizations, respectively. A weightedaverage of the four different measurements is calculated to obtain the particle fractions in | η | < .
8. Theyields are obtained by multiplying the particle fractions by the measured unidentified charged particlespectrum.
The raw particle distributions are normalized to the total number of events analysed in each multiplicityclass. To obtain the p T distributions of primary π , K, and p, the raw particle distributions obtained fromthe different PID approaches need to be also corrected for the detector efficiency and acceptance, theITS-TPC, and TPC-TOF matching efficiency, the PID efficiency, the trigger efficiency and the contami-nation from secondary particles.Secondary particles are either produced in weak decays or from the interaction of particles with the de-tector material. The estimation of secondary particle contribution is based on the Monte Carlo (MC)templates of the distance of closest approach of the track to the primary vertex in the transverse planewith respect to the beam axis (DCA xy ), as carried out in previous works [11, 23]. The DCA xy distri-butions of the tracks in data are fitted with three MC templates corresponding to the expected shapesof primary particles, secondaries from material and secondaries from weak decays to obtain the correctfraction of primary particles in the data. This procedure is repeated in each p T and multiplicity bin andthus takes into account the possible differences in the feed-down corrections due to the change in theabundances and spectral shapes of the weakly decaying particles. The contamination is different in eachPID analysis due to different track selection criteria and PID techniques and hence it is estimated sepa-rately for each analysis. The contribution of secondary particles was found to be significant for π (up to2%) and p (up to 15%) whereas the contribution for K is negligible.The spectra are corrected for the detector acceptance and track reconstruction efficiencies based on asimulation using the P YTHIA − with the detector material, an additional correction factor to the efficiency for these two par-ticles is estimated with GEANT4 [29] and FLUKA [30], respectively, where the interaction processesare known to be better reproduced [23]. Additional corrections to the efficiency are applied when TPCor TOF information is used to take into account the track matching between ITS and TPC, and betweenTPC and TOF.Signal losses due to the trigger selection are extracted from P YTHIA p T in the V0M class X (the lowestmultiplicity), and reduces to ∼ ∼
2% in classes IX and VIII, respectively. The correction is negligiblein higher multiplicity pp collisions and for p T (cid:38) c in all multiplicity bins except in class X. In thelatter, the correction reaches ∼
2% at p T = 7 GeV/ c . Finally, an additional correction is applied to passfrom triggered INEL > > | η true | < | V truez | <
10 cm. The correction is independentof particle species and is found to be negligible from V0M I (the highest multiplicity) to V0M VI, whileit ranges from 1% in class VII to 11% in class X. The correction is about 8% for multiplicity-integratedINEL > Events that passed all the selection criteria.
The systematic uncertainties are divided into two categories, those common to all analyses and thosewhich are analysis specific. The common systematic uncertainties are those due to tracking, which in-cludes track quality criteria and the p T -dependent ITS-TPC matching efficiency (except for the ITSsaanalysis), the TPC-TOF matching efficiency (for TPC-TOF and TOF analyses), and the signal loss cor-rection. In addition, the systematic source related to the effect of the material budget on the globaltracking ( p T dependent) is also added. The uncertainties on global tracking and TPC-TOF matchingdue to material budget are calculated by varying the material budget in the simulation by ± xy , on the χ of the track, and on the number of clusters required in the ITS layers.The uncertainty related to the particle identification is calculated by using a Bayesian technique andcomparing the results obtained with the standard n σ method as already performed in [31]. Due to theLorentz force, the positions of ITS clusters are shifted depending on the magnetic field polarity, givingrise to a 3% uncertainty. Finally, the energy-independent uncertainty related to the ITS material budget isestimated with a simulation of pp collisions at √ s = 900 GeV by varying the material budget of the ITS by ± p T (below 500, 600, and 800 MeV / c for π , K, and p,respectively), the systematic uncertainty associated with the PID technique is calculated by integratingthe measured d E / d x of charged tracks in the ranges of ± . σ and ± . σ , where σ represents onestandard deviation from the (cid:104) d E / d x (cid:105) under given mass hypothesis. At higher p T values, where only thetime-of-flight information is used, the associated uncertainties are calculated by simultaneously varyingthe width and tail parameters by 10%. An additional uncertainty is calculated by fixing the centralvalues of the fit functions to the β calculated for each particle species in a given momentum range. Thiswas found to be the dominant source of systematic uncertainty for π and K at the highest p T values( (cid:38) . / c ). For the TOF template fits analysis, PID uncertainties are estimated by simultaneouslyvarying the spread and tail slope of each ∆ t template by 10%. In addition to this, for both the TPC-TOF and TOF template fits analyses, systematic uncertainties associated with tracking are calculatedby varying the track selection criteria: the number of crossed rows in the TPC, the distance of closestapproach in beam and transverse directions, and the quality of the global track fit χ . For the kink analysisthe sources of systematic uncertainties are: the kink vertex finding efficiency (3% constant in p T ), thekink PID efficiency (calculated by taking into account the position of the kink vertex), the number ofTPC clusters of the decaying particle track and the q T of the decay product, and the uncertainty relatedto the purity of the selected sample. The contamination due to the random association of tracks wronglyattributed to kaon decays is of the order of 2.3% at low transverse momenta and reaches the value of 3.4%above 4 GeV/ c . The largest component of the systematic uncertainties in the analysis of the relativisticrise of the TPC arises from the imprecise parametrization of both the Bethe-Bloch and resolution curves.To quantify this uncertainty, the variations of the Bethe-Bloch resolution parametrizations with respectto the measured (cid:104) d E / d x (cid:105) ( σ (cid:104) d E / d x (cid:105) ) are used to vary the values of the mean and σ in the 4-Gaussianfit [26]. The largest relative deviation between the nominal particle ratios and the ones obtained after thevariations are assigned as a systematic uncertainty. 7LICE Collaboration Table 3:
Sources of the relative systematic uncertainties of the p T -differential yields of π , K, and p. The uncer-tainties are split into two categories, the common and the individual-analysis specific for low, intermediate andhigh p T . Numbers in parenthesis in the p column refer to p uncertainties. In the last rows, the maximum (amongmultiplicity classes) total systematic uncertainty is reported. Uncertainty (%)
Common source π K p (p) p T (GeV / c ) 0.1 3.0 20.0 0.2 2.5 20.0 0.3 4.0 20.0Correction for secondaries 1 1 1 negl. 4 1 1Hadronic interactions 2 2.4 2.4 2.7 1.8 1.8 1(3.6) 1(3.6) 1(3.6)ITS-TPC matching efficiency 0.7 1.5 2.9 0.7 1.5 2.9 0.7 1.5 2.9Global tracking efficiency 0.7 0.5 1.5TOF matching efficiency(TPC-TOF fits,TOF template fits) 3 6 4Signal-loss correction 0.2 1 3.3 p T (GeV / c ) 0.3 3.0 20.0 0.3 2.5 20.0 0.4 4.0 20.0Material budget(TPC-TOF) 0.5 1.0 0.2 1.5 1.0 0.4 2.9 1.7 0.1 Specific source π K p (p)
ITSsa , p T (GeV / c ) 0.10 0.70 0.20 0.60 0.30 0.65Tracking 1.4 1.4 1.5 1.5 1.1 1.1Material budget 4.8 0.3 2.3 0.6 5.0 0.9 E × B effect 3.0 3.0 3.0PID 0.4 0.4 3.9 3.9 4.2 4.2 TOF templates , p T (GeV / c ) 0.7 4.0 0.6 3.0 0.9 4.0PID 1 9.4 1 12 1 21 TPC-TOF fits , p T (GeV / c ) 0.3 3 0.3 3 0.4 3PID 1.4 7 3 16 1 4.3 rTPC , p T (GeV / c ) 2.0 20.0 2.0 20.0 2.0 20.0Bethe–Bloch parameterization 8.1 4 14.8 8.0 14 6.0Feed-Down 0.5 0.5 - - 2.4 2.0 Kinks , p T (GeV / c ) 0.6 6.0 0.6 6.0 0.6 6.0PID - - 0.75 5.3 - -Kink vertex finding efficiency - - 3 3 - -Contamination - - 2.3 3.4 - - Total π K p (p) p T (GeV / c ) 0.3 3.0 20.0 0.3 3.0 20.0 0.3 3.0 20.0Total 6.6 5.1 4.3 6.6 6.8 7.8 9.6 9.7 12.4 Particle ratios K/ π p/ π p T (GeV / c ) 0.2 3.0 20.0 0.3 3.0 20.0Total 7.2 14.6 7.7 10.2 12.2 11.58LICE Collaboration The p T -differential spectra of π , K, and p measured as a function of the charged-particle multiplicitydensity in pp collisions at √ s =
13 TeV are shown in Fig. 1. For each V0M class, charged-particlemultiplicity density has been measured in the central region ( | η | < . > p T spectra become harder with increasing (cid:104) d N ch / d η (cid:105) , and the effect is morepronounced for protons. The hardening of the inclusive charged-hadron spectra with (cid:104) d N ch / d η (cid:105) has beenalso recently reported in [33], where different MPI models were shown to describe such effect. On theother hand, the mass dependence of spectral shape modifications is also observed in Pb–Pb collisions at √ s NN = .
76 TeV [26], where it is usually associated with the hydrodynamical evolution of the system.At higher p T ( (cid:38) / c ), we find that slopes of particle spectra become independent of the multiplicityclass considered, as expected from pQCD calculations [34].The p T -differential K/ π and p/ π ratios as a function of (cid:104) d N ch / d η (cid:105) measured at low and intermedi-ate transverse momenta are shown in Fig. 2 together with those measured in pp collisions at √ s = √ s =
13 TeV. The mea-sured K/ π ratio shows no evident sign of evolution with multiplicity in all p T ranges considered, whilethe p/ π ratio shows depletion at low p T , an increase at intermediate p T , and constant behavior at high p T .In addition, the measured K/ π and p/ π ratios are consistent between the two center-of-mass energies [11].For MC predictions, the event classification is based on the number of charged tracks simulated at for-ward and backward pseudorapidities covered by the V0 detector, in a way similar to the way the eventclassification is done for the data. The mean charged-particle multiplicity density is then calculated inthe central pseudorapidity region, | η | < .
5. HERWIG 7 [35, 36], where a clustering approach is usedfor hadronization, provides a good description of the evolution of the K/ π and p/ π ratios with (cid:104) d N ch / d η (cid:105) in the low and intermediate p T ranges and is consistent with the measured ratios within 1-2 standarddeviations. P YTHIA π and p/ π ratios.The CR scheme, which has been shown to capture the modifications of the baryon-to-meson ratios [3],provides only a qualitative description of the evolution of the p/ π ratio with (cid:104) d N ch / d η (cid:105) and overesti-mates the absolute values of the ratio at low and high p T . The implementation of color ropes [5, 38, 39]in P YTHIA
8, which results in higher effective string tension and thus enhances strange- and di-quarkproduction, provides a qualitative description K/ π (p/ π ) ratio only at low (intermediate) p T and overes-timates the p/ π ratio at low p T . This could be understood considering that larger effective string tensionis mostly translated to hadronic mass and thus feeds down the low p T part of the spectrum.In large collision systems such as Pb–Pb, multiplicity-dependent modifications of hadron p T spectracan be interpreted as the hydrodynamical radial expansion of the system and studied in the context ofthe Boltzmann-Gibbs Blast-Wave model [40]. In this model, a thermalized medium expands radiallyand undergoes an instantaneous kinematic freeze-out. The average expansion velocity (cid:104) β T (cid:105) , the kineticfreeze-out temperature T kin , and the velocity profile exponent n can be extracted from simultaneousmodel fits to hadron spectra. As the trends observed in the evolution of particle spectra measured in ppcollisions are highly reminiscent to those in p–Pb and Pb–Pb, it is interesting to check whether the Blast-Wave model can be extended to describe pp collisions. Such study has been previously reported in [11],where pp, p–Pb, and Pb–Pb collisions at √ s NN =
7, 5.02, and 2.76 TeV were considered. Now, for thefirst time, we can study the evolution of (cid:104) β T (cid:105) , T kin and n in pp collisions as a function of the collisionenergy.At low transverse momenta ( p T (cid:46)
500 MeV / c ), the dominant mechanism of π production is from res-onance decays. To account for this in the Blast-Wave model fits, spectral measurements of all stronglydecaying hadrons are required. Alternatively, one can choose to omit the low- p T pions. Noting that thereis a strong dependence of Blast-Wave parameters on the fitting range [23], it is important to consider9LICE Collaboration ) c (GeV/ T p - ) c ) ( G e V / T p d y / ( d N d e v t N / - - - - - - -
10 110 V0M multiplicity classes) · I ( ) · II () · III ( ) · IV () · V ( ) · VI () · VII ( ) · VIII () · IX ( ) · X ( | < 0.5 y | = 13 TeV s pp, ALICE - p + + p =26.02 æh /d ch N d Æ =2.55 æh /d ch N d Æ ) c (GeV/ T p R a t i o t o I N E L > -
10 110 ) c (GeV/ T p - ) c ) ( G e V / T p d y / ( d N d e v t N / - - - - - - -
10 110 - +K + K =26.02 æh /d ch N d Æ =2.55 æh /d ch N d Æ ) c (GeV/ T p R a t i o t o I N E L > -
10 110 ) c (GeV/ T p - ) c ) ( G e V / T p d y / ( d N d e v t N / - - - - - - - -
10 110 pp+ =26.02 æh /d ch N d Æ =2.55 æh /d ch N d Æ ) c (GeV/ T p R a t i o t o I N E L > -
10 110
Figure 1:
Transverse momentum spectra of π , K, and p for different multiplicity event classes. Spectra are scaledby powers of 2 in order to improve visibility. The corresponding ratios to INEL > ) - p + + p ) / ( - + K + ( K ALICE |<0.5 y = 13 TeV, | s pp |<0.5 y = 7 TeV, | s pp c < 0.55 GeV/ T p ) - p + + p ) / ( p ( p + (x 4) ) - p + + p ) / ( - + K + ( K PYTHIA8 + color ropesHERWIG7PYTHIA8 MonashPYTHIA8 Monash, NoCR c < 2.60 GeV/ T p ) - p + + p ) / ( p ( p + |<0.5 h | æh /d ch N d Æ ) - p + + p ) / ( - + K + ( K c < 20.00 GeV/ T p ) - p + + p ) / ( p ( p + (x 1.5) Figure 2:
Multiplicity dependence of p T -differential K/ π (upper panels) and p/ π (lower panels) ratios measured inpp collisions at √ s = 7 TeV [11] and 13 TeV (blue and red, respectively). Lines represent different MC generatorpredictions for pp collisions at √ s =
13 TeV. Left to right: low-, intermediate-, and high-transverse momenta.Vertical bars, open, and shaded bands represent statistical, total systematic, and multiplicity uncorrelated system-atic uncertainties, respectively. the same p T range in the fitting procedure in order to obtain a consistent comparison between differ-ent colliding systems. The comparison of the (cid:104) β T (cid:105) - T kin correlations measured in different systems andcenter-of-mass energies is shown in Fig. 3. In this paper we consider three different approaches to theBlast-Wave model fits to particle spectra measured in pp collisions at √ s =
13 TeV: a) traditionalfits as done in [10, 11, 23], where π , K, and p spectra are fitted and resonance feed-down is neglected(represented by markers in Fig. 3), b) simultaneously fitting K, p, and Λ spectra [22] noting that Λ arenot significantly affected by resonance decays (represented by a solid line in Fig. 3), and c) a methodproposed in [41, 42], where the resonance feed-down is calculated before the Cooper–Frye freeze-outusing a statistical hadronization model (represented by a dashed line in Fig. 3). We find that the T kin - (cid:104) β T (cid:105) correlation in pp collisions at √ s =
13 TeV follows similar trends as seen at lower energies. When Λ ’s are considered instead of pions, the trends seen in (cid:104) β T (cid:105) - T kin correlation do not change significantlyand only at highest multiplicities we find a larger T kin . On the other hand, when a proper treatment ofresonance decays is used, we find a significantly lower T kin of around 135 MeV at the lowest multiplic-ities, which then grows with increasing (cid:104) d N ch / d η (cid:105) and approaches the pseudocritical QCD temperature T c = ± . (cid:104) β T (cid:105) , T kin ,and n with (cid:104) d N ch / d η (cid:105) is shown in Fig. 4 for different colliding systems. From the lowest multiplicities, T kin grows with (cid:104) β T (cid:105) until it saturates at around 180 MeV. At larger multiplicities ( (cid:104) d N ch / d η (cid:105) (cid:38) T kin decreases and becomes similar to that measured in p–Pb collisions at √ s = .
02 TeV, suggesting thatthe system decouples at lower temperature and thus is longer-lived. The average expansion velocity (cid:104) β T (cid:105) increases with (cid:104) d N ch / d η (cid:105) and its values are consistent for pp collisions at different √ s as well as withthe corresponding values for p–Pb collisions, indicating that small systems become more explosive atlarger multiplicities. In contrast to this, (cid:104) β T (cid:105) measured in Pb–Pb collisions is lower than that in smallersystems for the common (cid:104) d N ch / d η (cid:105) range, see Fig. 4. This indicates that the size of the colliding system11LICE Collaboration æ T bÆ ( G e V ) k i n T ) c ), p (0.3-3.0 GeV/ c ), K (0.2-1.5 GeV/ c (0.5-1.0 GeV/ p Global Blast-Wave fit range:ALICE = 7 TeV s pp, = 13 TeV s pp, = 5.02 TeV NN s Pb, - p = 5.02 TeV NN s Pb, - Pb p ), no c (0.6-3.0 GeV/ L = 13 TeV, s pp, = 13 TeV, FastReso s pp, Figure 3:
Correlation of kinetic freeze-out temperature T kin and average expansion velocity (cid:104) β T (cid:105) , extracted fromsimultaneous Blast-Wave fits to π , K, and p spectra measured in pp, p–Pb, and Pb–Pb collisions. Contours repre-sent 1 σ uncertainty. The shaded ellipses represent the T kin - (cid:104) β T (cid:105) correlation calculated from Blast-Wave fit to K, p,and Λ spectra [22] measured in pp collisions at √ s =
13 TeV. The empty circles represent Blast-Wave fits withresonance decays [42]. References from [10–12, 42]. might have significant effects on the final state particle dynamics. This is also reflected in the expansionvelocity profile power n shown in Fig. 4: in pp and p–Pb collisions, large n suggests high pressure gra-dients which lead to larger (cid:104) β T (cid:105) , while in Pb–Pb collisions, n ∼ p T -integrated particle yields with (cid:104) d N ch / d η (cid:105) [11], which extends across different colliding systems andcollision energies. Now for the first time we can isolate the center-of-mass energy dependence of thisscaling for π , K, and p in pp collisions. The p T -integrated particle yields (d N / d y ) and average transversemomenta ( (cid:104) p T (cid:105) ) are calculated by integrating the measured transverse momentum spectra and using theLévy-Tsallis parametrization [45–47] to extrapolate to the low p T regions not covered by the measure-ments. The extrapolated fractions of the yields at low p T are 8% (10%) for π , 6% (13%) for K, and7% (20%) for p for the highest (lowest) multiplicities. For systematic uncertainties on the extrapolation,Bylinkin, Bose-Einstein, Fermi-Dirac, m T -exponential and Hagedorn functions are used to fit particlespectra. The statistical uncertainties of d N / d y and (cid:104) p T (cid:105) are calculated by coherently shifting the cen-tral values of each spectra point by a fraction of its statistical uncertainty. The fraction is randomlydrawn from Gaussian distribution and new values of integrated yields and mean transverse momentaare calculated. The procedure is repeated 1000 times to calculate the standard deviations of d N / d y and (cid:104) p T (cid:105) , which are then used as the statistical uncertainties. To estimate the systematic uncertainty on theintegrated yields, the spectra points are moved to maximal/minimal values allowed by their respectivesystematic uncertainties before repeating the fit procedure. For (cid:104) p T (cid:105) , each point of the spectra is shiftedto the upper/lower edge of the corresponding p T bin to obtain the hardest/softest particle distribution.The largest differences to the nominal yield and (cid:104) p T (cid:105) values are combined with the extrapolation un-certainties to calculate the total systematic uncertainties. The kaon- and proton-to-pion integrated yieldratios measured in pp collisions at √ s =
13 TeV are found to be in a good agreement within systematicuncertainties with those measured in pp, p–Pb, and Pb–Pb collisions at √ s NN =
7, 5.02, and 2.76 TeV,12LICE Collaboration |< 0.5 h | æh /d ch N d Æ æ T b Æ ALICE = 7 TeV s pp = 5.02 TeV NN s p-Pb = 2.76 TeV NN s Pb-Pb = 5.02 TeV NN s Pb-Pb = 13 TeV s pp Global Blast-Wave fit) c (0.5-1 GeV/ p ) c K (0.2-1.5 GeV/ ) c p (0.3-3.0 GeV/ |< 0.5 h | æh /d ch N d Æ ) c ( G e V / k i n T |< 0.5 h | æh /d ch N d Æ n Figure 4:
Evolution of (cid:104) β T (cid:105) , T kin , and n with (cid:104) d N ch / d η (cid:105) . (cid:104) β T (cid:105) , T kin and n are extracted from Blast-Wave fits to π ,K, and p p T spectra measured in pp, p–Pb, and Pb–Pb collisions at different √ s . respectively, as shown in Fig. 5. This indicates that hadrochemistry at LHC energies scales with charged-particle multiplicity density in a uniform way, despite the colliding system or collision energy.The description of hadron-to-pion ratio factorization with multiplicity at lower center-of-mass energiesin MC generators has been previously shown to be qualitative at best [11]. In fact, both P YTHIA (cid:104) d N ch / d η (cid:105) . In thispaper, we consider more recent versions of the two MC generators. In particular, the hadronizationin P YTHIA
YTHIA Ω / π ratio at the lowest multiplicity in pp collisions. This might pointto the need of a further refinement of MC generator tuning, as similar discrepancies are already observedfor e + e − data [48]. The integrated K/ π yield ratio shown in Fig. 5 at high multiplicity pp collisions arecaptured by P YTHIA (cid:104) d N ch / d η (cid:105) which is not observed in the data. The predictions from P YTHIA π ratios in pp collisions at √ s =
13 TeV, whether color reconnection is considered ornot. The quantitative description of p/ π ratio is given only by HERWIG 7, while all considered versionsof P YTHIA
YTHIA π ratio with (cid:104) d N ch / d η (cid:105) , which could be attributed to the enhanced production of strange- and di-quark inthe rope fragmentation. Overall, we conclude that none of the models considered provide a consistentdescription of the data.The average transverse momenta of identified particles are found to increase with multiplicity in ppcollisions at √ s = (cid:104) p T (cid:105) . Similar observations have been previouslyreported in pp [49] and p–Pb [10] collisions at lower energies and for strange hadrons in pp collisions at √ s =
13 TeV [22]. The solid red line in Fig. 7 represents a fit of the form a − b / ( c − (cid:104) d N ch / d η (cid:105) ) to the √ s =
13 TeV data, which is then used for a better comparison of (cid:104) p T (cid:105) between the two center-of-massenergies, see lower panels of the same figure. We find a small hint of an increase with √ s for similarmultiplicities for π , while the (cid:104) p T (cid:105) of protons is similar at the two center-of-mass energies. Note thatsimilar observations have been already reported in [22], where spectra of K were found to become harderwith √ s at similar multiplicities. In addition, we find that P YTHIA
YTHIA π (cid:104) p T (cid:105) evolution with (cid:104) d N ch / d η (cid:105) .This is expected as pions are the most abundant particles produced in collisions, and the three generatorsare tuned to explicitly to describe the (cid:104) p T (cid:105) of charged-particles. On the other hand, we observe thatthe (cid:104) p T (cid:105) of K and p are well described only by HERWIG 7, while P YTHIA (cid:104) p T (cid:105) in the whole (cid:104) d N ch / d η (cid:105) range considered. This could be understood consideringthat the additional energy available during the rope fragmentation predominantly enhances the productionof heavier hadrons at low p T . We have studied π , K, and p production as a function of multiplicity in pp collisions at √ s =
13 TeV.To avoid auto-correlation biases, the event classification has been based on multiplicity measurements atforward (backward) pseudorapidity, while event activity (cid:104) d N ch / d η (cid:105) has been correspondingly estimatedat central pseudorapidities, | η | < .
5. We find that hadron p T spectra become harder with multiplicity,and the effect is more pronounced for heavier particles. The hardening of the spectra is predicted byP YTHIA
YTHIA π (cid:104) p T (cid:105) , while K and p aredescribed qualitatively only by HERWIG7. At high p T ( (cid:38) / c ) we find that spectral shapes becomeindependent of (cid:104) d N ch / d η (cid:105) as predicted by pQCD calculations [34].The measured p T -differential K/ π ratios show no evolution with multiplicity in the p T range considered.In contrast to this, a depletion (enhancement, saturation) is visible for the p/ π ratios at low (intermediate,high) p T . In addition, we find that the ratios measured in pp collisions at √ s =
13 TeV are consistentwith those measured at √ s = p T -differential K/ π and p/ π ratios, none of them provides a consistent description of the data.The saturation at high p T is captured by P YTHIA
YTHIA π and p/ π ratios at high p T .The study of hadron p T spectra in the context of the Blast-Wave model reveals that the kinetic freeze-outtemperature T kin , average expansion velocity (cid:104) β T (cid:105) , and the velocity profile exponent n show little or nodependence on the center-of-mass energy and are consistent within uncertainties with those extractedfrom particle spectra measured in pp collisions at √ s = (cid:104) d N ch / d η (cid:105) .14LICE Collaboration |<0.5 h | æh /d ch N d Æ ) - p + + p ) / ( - + K + ( K = 13 TeV s pp PYTHIA8 + color ropes = 7 TeV s pp HERWIG7 = 5.02 TeV NN s Pb - p PYTHIA8 Monash, NoCR = 5.02 TeV NN s Pb - Pb PYTHIA8 Monash = 2.76 TeV NN s Pb - Pb ALICE |<0.5 h | æh /d ch N d Æ ) - p + + p ) / ( p ( p + = 13 TeV s pp PYTHIA8 + color ropes = 7 TeV s pp HERWIG7 = 5.02 TeV NN s Pb - p PYTHIA8 Monash, NoCR = 5.02 TeV NN s Pb - Pb PYTHIA8 Monash = 2.76 TeV NN s Pb - Pb ALICE
Figure 5:
Integrated K/ π (top) and p/ π (bottom) yield ratios as a function of charged-particle multiplicity densitymeasured in pp, p–Pb, and Pb–Pb collisions at different center-of-mass energies. Empty (shaded) boxes representtotal (multiplicity uncorrelated) systematic uncertainties. Black lines represent predictions from different MCgenerators for pp collisions at √ s =
13 TeV. References from [11, 23, 31, 50]. |< 0.5 h | æh /d ch N d Æ ) -p + + p R a t i o o f y i e l d s t o ( - - - · ( + W + - W · ( + X + - X · ( L + L S0 ALICE = 13 TeV s pp, = 7 TeV s pp, = 5.02 TeV NN s p-Pb, = 2.76 TeV NN s Pb-Pb, PYTHIA8 + color ropesHERWIG7PYTHIA8 MonashPYTHIA8 Monash, NoCR
Figure 6:
Integrated strange hadron-to-pion ratios as a function of (cid:104) d N ch / d η (cid:105) measured in pp, p–Pb, and Pb–Pbcollisions. The open (shaded) boxes around markers represent full (multiplicity uncorrelated) systematic uncer-tainties. Different lines represent predictions from different MC generators for pp collisions at √ s =
13 TeV.References from [6, 11, 12, 22].
10 20 ( G e V / c ) æ T p Æ - p + + p ALICE = 13 TeV s pp, = 7 TeV s pp,
10 20 ( G e V / c ) æ T p Æ - +K + K
10 20 ( G e V / c ) æ T p Æ pp+
10 20 D a t a / F i t |<0.5 h | æh /d ch N d Æ
10 20 ( G e V / c ) æ T p Æ ( G e V / c ) æ T p Æ Figure 7:
Upper panels: average transverse momenta of π , K, and p as a function of charged-particle multiplicitydensity measured in pp collisions at √ s = a − b ( c − (cid:104) d N ch / d η (cid:105) ) − fit to the 13 TeV data to guide the eye. Open (shaded) boxes represent total (multiplicity uncorrelated) systematicuncertainties. Black lines represent predictions from different MC generators for pp collisions at √ s =
13 TeV.Bottom panels: ratios of (cid:104) p T (cid:105) to the fits. Data at √ s = The p T -integrated hadron-to-pion ratios as a function of multiplicity show no center-of-mass dependenceand the measurement in pp collisions at √ s =
13 TeV are compatible to those in pp, p–Pb, and Pb–Pbcollisions at √ s NN =
7, 5.02, and 2.76 TeV, respectively. This suggests that, at the LHC energies, thechemical composition of primary hadrons scales with charged-particle multiplicity density in a uniformway, despite the colliding system and collision energy. Comparisons of the integrated hadron-to-pionratios to the predictions from MC generators show that P
YTHIA π ratio. HERWIG7 also capturesthe evolution of the ratios with (cid:104) d N ch / d η (cid:105) , but underestimates the absolute values of Ξ / π and Ω / π .Overall, none of the generators are able to provide a consistent quantitative description of the measuredhadron-to-pion ratios. Acknowledgements
The ALICE Collaboration would like to thank all its engineers and technicians for their invaluable con-tributions to the construction of the experiment and the CERN accelerator teams for the outstandingperformance of the LHC complex. The ALICE Collaboration gratefully acknowledges the resources andsupport provided by all Grid centres and the Worldwide LHC Computing Grid (WLCG) collaboration.The ALICE Collaboration acknowledges the following funding agencies for their support in buildingand running the ALICE detector: A. I. Alikhanyan National Science Laboratory (Yerevan Physics In-stitute) Foundation (ANSL), State Committee of Science and World Federation of Scientists (WFS),Armenia; Austrian Academy of Sciences, Austrian Science Fund (FWF): [M 2467-N36] and National-stiftung für Forschung, Technologie und Entwicklung, Austria; Ministry of Communications and HighTechnologies, National Nuclear Research Center, Azerbaijan; Conselho Nacional de DesenvolvimentoCientífico e Tecnológico (CNPq), Financiadora de Estudos e Projetos (Finep), Fundação de Amparo àPesquisa do Estado de São Paulo (FAPESP) and Universidade Federal do Rio Grande do Sul (UFRGS),Brazil; Ministry of Education of China (MOEC) , Ministry of Science & Technology of China (MSTC)17LICE Collaborationand National Natural Science Foundation of China (NSFC), China; Ministry of Science and Educationand Croatian Science Foundation, Croatia; Centro de Aplicaciones Tecnológicas y Desarrollo Nuclear(CEADEN), Cubaenergía, Cuba; Ministry of Education, Youth and Sports of the Czech Republic, CzechRepublic; The Danish Council for Independent Research | Natural Sciences, the VILLUM FONDEN andDanish National Research Foundation (DNRF), Denmark; Helsinki Institute of Physics (HIP), Finland;Commissariat à l’Energie Atomique (CEA), Institut National de Physique Nucléaire et de Physique desParticules (IN2P3) and Centre National de la Recherche Scientifique (CNRS) and Région des Pays dela Loire, France; Bundesministerium für Bildung und Forschung (BMBF) and GSI Helmholtzzentrumfür Schwerionenforschung GmbH, Germany; General Secretariat for Research and Technology, Ministryof Education, Research and Religions, Greece; National Research, Development and Innovation Office,Hungary; Department of Atomic Energy Government of India (DAE), Department of Science and Tech-nology, Government of India (DST), University Grants Commission, Government of India (UGC) andCouncil of Scientific and Industrial Research (CSIR), India; Indonesian Institute of Science, Indonesia;Centro Fermi - Museo Storico della Fisica e Centro Studi e Ricerche Enrico Fermi and Istituto Nazionaledi Fisica Nucleare (INFN), Italy; Institute for Innovative Science and Technology , Nagasaki Instituteof Applied Science (IIST), Japanese Ministry of Education, Culture, Sports, Science and Technology(MEXT) and Japan Society for the Promotion of Science (JSPS) KAKENHI, Japan; Consejo Nacionalde Ciencia (CONACYT) y Tecnología, through Fondo de Cooperación Internacional en Ciencia y Tec-nología (FONCICYT) and Dirección General de Asuntos del Personal Academico (DGAPA), Mexico;Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO), Netherlands; The Research Coun-cil of Norway, Norway; Commission on Science and Technology for Sustainable Development in theSouth (COMSATS), Pakistan; Pontificia Universidad Católica del Perú, Peru; Ministry of Science andHigher Education and National Science Centre, Poland; Korea Institute of Science and Technology In-formation and National Research Foundation of Korea (NRF), Republic of Korea; Ministry of Educationand Scientific Research, Institute of Atomic Physics and Ministry of Research and Innovation and Insti-tute of Atomic Physics, Romania; Joint Institute for Nuclear Research (JINR), Ministry of Education andScience of the Russian Federation, National Research Centre Kurchatov Institute, Russian Science Foun-dation and Russian Foundation for Basic Research, Russia; Ministry of Education, Science, Research andSport of the Slovak Republic, Slovakia; National Research Foundation of South Africa, South Africa;Swedish Research Council (VR) and Knut & Alice Wallenberg Foundation (KAW), Sweden; EuropeanOrganization for Nuclear Research, Switzerland; Suranaree University of Technology (SUT), NationalScience and Technology Development Agency (NSDTA) and Office of the Higher Education Commis-sion under NRU project of Thailand, Thailand; Turkish Atomic Energy Agency (TAEK), Turkey; Na-tional Academy of Sciences of Ukraine, Ukraine; Science and Technology Facilities Council (STFC),United Kingdom; National Science Foundation of the United States of America (NSF) and United StatesDepartment of Energy, Office of Nuclear Physics (DOE NP), United States of America.
References [1] R. Derradi de Souza, T. Koide, and T. Kodama, “Hydrodynamic Approaches in Relativistic HeavyIon Reactions,”
Prog. Part. Nucl. Phys. (2016) 35–85, arXiv:1506.03863 [nucl-th] .[2] A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel, “Decoding the phase structure ofQCD via particle production at high energy,” Nature no. 7723, (2018) 321–330, arXiv:1710.09425 [nucl-th] .[3] A. Ortiz Velasquez, P. Christiansen, E. Cuautle Flores, I. Maldonado Cervantes, and G. PaiÄ ˘G,“Color Reconnection and Flowlike Patterns in pp Collisions,”
Phys. Rev. Lett. no. 4, (2013)042001, arXiv:1303.6326 [hep-ph] .[4] T. SjÃ˝ustrand, S. Ask, J. R. Christiansen, R. Corke, N. Desai, P. Ilten, S. Mrenna, S. Prestel, C. O.18LICE CollaborationRasmussen, and P. Z. Skands, “An Introduction to PYTHIA 8.2,”
Comput. Phys. Commun. (2015) 159–177, arXiv:1410.3012 [hep-ph] .[5] C. Bierlich, G. Gustafson, L. LÃ˝unnblad, and A. Tarasov, “Effects of Overlapping Strings in ppCollisions,”
JHEP (2015) 148, arXiv:1412.6259 [hep-ph] .[6] ALICE
Collaboration, J. Adam et al. , “Multi-strange baryon production in p–Pb collisions at √ s NN = .
02 TeV,”
Phys. Lett.
B758 (2016) 389–401, arXiv:1512.07227 [nucl-ex] .[7]
CMS
Collaboration, S. Chatrchyan et al. , “Observation of long-range near-side angularcorrelations in proton-lead collisions at the LHC,”
Phys. Lett.
B718 (2013) 795–814, arXiv:1210.5482 [nucl-ex] .[8]
ALICE
Collaboration, B. Abelev et al. , “Long-range angular correlations on the near and awayside in p–Pb collisions at √ s NN = .
02 TeV,”
Phys. Lett.
B719 (2013) 29–41, arXiv:1212.2001[nucl-ex] .[9]
ALICE
Collaboration, B. Abelev et al. , “Multiparticle azimuthal correlations in p–Pb and Pb–Pbcollisions at the CERN Large Hadron Collider,”
Phys. Rev.
C90 no. 5, (2014) 054901, arXiv:1406.2474 [nucl-ex] .[10]
ALICE
Collaboration, B. Abelev et al. , “Multiplicity Dependence of Pion, Kaon, Proton andLambda Production in p–Pb Collisions at √ s NN = 5.02 TeV,” Phys. Lett.
B728 (2014) 25–38, arXiv:1307.6796 [nucl-ex] .[11]
ALICE
Collaboration, S. Acharya et al. , “Multiplicity dependence of light-flavor hadronproduction in pp collisions at √ s = 7 TeV,” Phys. Rev.
C99 no. 2, (2019) 024906, arXiv:1807.11321 [nucl-ex] .[12]
ALICE
Collaboration, J. Adam et al. , “Enhanced production of multi-strange hadrons inhigh-multiplicity proton-proton collisions,”
Nature Phys. (2017) 535–539, arXiv:1606.07424[nucl-ex] .[13] ATLAS
Collaboration, G. Aad et al. , “Observation of Long-Range Elliptic AzimuthalAnisotropies in √ s =
13 and 2.76 TeV pp Collisions with the ATLAS Detector,”
Phys. Rev. Lett. no. 17, (2016) 172301, arXiv:1509.04776 [hep-ex] .[14]
ALICE
Collaboration, “The ALICE definition of primary particles,”. https://cds.cern.ch/record/2270008 .[15]
ALICE
Collaboration, K. Aamodt et al. , “The ALICE experiment at the CERN LHC,”
JINST (2008) S08002.[16] ALICE
Collaboration, B. Abelev et al. , “Performance of the ALICE Experiment at the CERNLHC,”
Int. J. Mod. Phys.
A29 (2014) 1430044, arXiv:1402.4476 [nucl-ex] .[17]
ALICE
Collaboration, E. Abbas et al. , “Mid-rapidity anti-baryon to baryon ratios in pp collisionsat √ s = 0.9, 2.76 and 7 TeV measured by ALICE,” Eur. Phys. J.
C73 (2013) 2496, arXiv:1305.1562 [nucl-ex] .[18]
ALICE
Collaboration, “ALICE luminosity determination for pp collisions at √ s =
13 TeV,”. http://cds.cern.ch/record/2160174 .[19] J. Alme et al. , “The ALICE TPC, a large 3-dimensional tracking device with fast readout forultra-high multiplicity events,”
Nucl. Instrum. Meth.
A622 (2010) 316–367, arXiv:1001.1950[physics.ins-det] . 19LICE Collaboration[20]
ALICE
Collaboration, D. De Gruttola, “Particle IDentification with the ALICE Time-Of-Flightdetector at the LHC,”
JINST no. 10, (2014) C10019.[21] ALICE
Collaboration, J. Adam et al. , “Determination of the event collision time with the ALICEdetector at the LHC,”
Eur. Phys. J. Plus no. 2, (2017) 99, arXiv:1610.03055[physics.ins-det] .[22]
ALICE
Collaboration, S. Acharya et al. , “Multiplicity dependence of (multi-)strange hadronproduction in proton-proton collisions at √ s = 13 TeV,” arXiv:1908.01861 [nucl-ex] .[23] ALICE
Collaboration, B. Abelev et al. , “Centrality dependence of π , K, p production in Pb–Pbcollisions at √ s NN = 2.76 TeV,” Phys. Rev.
C88 (2013) 044910, arXiv:1303.0737 [hep-ex] .[24]
ALICE
Collaboration, J. Adam et al. , “Charged-particle multiplicities in protonâ ˘A ¸Sprotoncollisions at √ s = . Eur. Phys. J.
C77 no. 1, (2017) 33, arXiv:1509.07541[nucl-ex] .[25]
ALICE
Collaboration, J. Adam et al. , “Measurement of pion, kaon and proton production inprotonâ ˘A ¸Sproton collisions at √ s = Eur. Phys. J.
C75 no. 5, (2015) 226, arXiv:1504.00024 [nucl-ex] .[26]
ALICE
Collaboration, J. Adam et al. , “Centrality dependence of the nuclear modification factorof charged pions, kaons, and protons in Pb–Pb collisions at √ s NN = .
76 TeV,”
Phys. Rev.
C93 no. 3, (2016) 034913, arXiv:1506.07287 [nucl-ex] .[27] P. Skands, S. Carrazza, and J. Rojo, “Tuning PYTHIA 8.1: the Monash 2013 Tune,”
Eur. Phys. J.
C74 no. 8, (2014) 3024, arXiv:1404.5630 [hep-ph] .[28] R. Brun, F. Bruyant, M. Maire, A. C. McPherson, and P. Zanarini,
GEANT 3: user’s guide Geant3.10, Geant 3.11; rev. version . CERN, Geneva, 1987. https://cds.cern.ch/record/1119728 .[29]
GEANT4
Collaboration, S. Agostinelli et al. , “GEANT4: A Simulation toolkit,”
Nucl. Instrum.Meth.
A506 (2003) 250–303.[30] A. Ferrari, P. R. Sala, A. Fasso, and J. Ranft, “FLUKA: A multi-particle transport code (Programversion 2005),”.[31]
ALICE
Collaboration, J. Adam et al. , “Multiplicity dependence of charged pion, kaon, and(anti)proton production at large transverse momentum in p–Pb collisions at √ s NN = 5.02 TeV,” Phys. Lett.
B760 (2016) 720–735, arXiv:1601.03658 [nucl-ex] .[32]
ALICE
Collaboration, K. Aamodt et al. , “Production of pions, kaons and protons in pp collisionsat √ s =
900 GeV with ALICE at the LHC,”
Eur. Phys. J.
C71 (2011) 1655, arXiv:1101.4110[hep-ex] .[33]
ALICE
Collaboration, S. Acharya et al. , “Charged-particle production as a function ofmultiplicity and transverse spherocity in pp collisions at √ s = .
02 and 13 TeV,” arXiv:1905.07208 [nucl-ex] .[34] S. Kretzer, “Fragmentation functions from flavour inclusive and flavour tagged e + e − annihilations,” Phys. Rev.
D62 (2000) 054001, arXiv:hep-ph/0003177 [hep-ph] .[35] J. Bellm et al. , “Herwig 7.0/Herwig++ 3.0 release note,”
Eur. Phys. J.
C76 no. 4, (2016) 196, arXiv:1512.01178 [hep-ph] . 20LICE Collaboration[36] M. Bahr et al. , “Herwig++ Physics and Manual,”
Eur. Phys. J.
C58 (2008) 639–707, arXiv:0803.0883 [hep-ph] .[37] T. Sjostrand, S. Mrenna, and P. Z. Skands, “A Brief Introduction to PYTHIA 8.1,”
Comput. Phys.Commun. (2008) 852–867, arXiv:0710.3820 [hep-ph] .[38] C. Flensburg, G. Gustafson, and L. Lönnblad, “Inclusive and Exclusive Observables from Dipolesin High Energy Collisions,”
JHEP (2011) 103, arXiv:1103.4321 [hep-ph] .[39] C. Bierlich, G. Gustafson, L. LÃ˝unnblad, and H. Shah, “The Angantyr model for Heavy-IonCollisions in PYTHIA8,” JHEP (2018) 134, arXiv:1806.10820 [hep-ph] .[40] E. Schnedermann, J. Sollfrank, and U. W. Heinz, “Thermal phenomenology of hadrons from200-AGeV S+S collisions,” Phys. Rev.
C48 (1993) 2462–2475, arXiv:nucl-th/9307020[nucl-th] .[41] A. Mazeliauskas, S. Floerchinger, E. Grossi, and D. Teaney, “Fast resonance decays in nuclearcollisions,” arXiv:1809.11049 [nucl-th] .[42] A. Mazeliauskas and V. Vislavicius, “Temperature and fluid velocity on the freeze-out surfacefrom π , K , p spectra in pp, p–Pb and Pb–Pb collisions,” arXiv:1907.11059 [hep-ph] .[43] HotQCD
Collaboration, A. Bazavov et al. , “Chiral crossover in QCD at zero and non-zerochemical potentials,”
Phys. Lett.
B795 (2019) 15–21, arXiv:1812.08235 [hep-lat] .[44] E. Shuryak and I. Zahed, “High-multiplicity pp and pA collisions: Hydrodynamics at its edge,”
Phys. Rev.
C88 no. 4, (2013) 044915, arXiv:1301.4470 [hep-ph] .[45] C. Tsallis, “Possible Generalization of Boltzmann-Gibbs Statistics,”
J. Statist. Phys. (1988)479–487.[46] G. Wilk and Z. Wlodarczyk, “On the interpretation of nonextensive parameter q in Tsallis statisticsand Levy distributions,” Phys. Rev. Lett. (2000) 2770, arXiv:hep-ph/9908459 [hep-ph] .[47] STAR
Collaboration, J. Adams et al. , “K(892)* resonance production in Au+Au and p+pcollisions at √ s NN =
200 GeV at STAR,”
Phys. Rev.
C71 (2005) 064902, arXiv:nucl-ex/0412019 [nucl-ex] .[48] C. Bierlich and J. R. Christiansen, “Effects of color reconnection on hadron flavor observables,”
Phys. Rev.
D92 no. 9, (2015) 094010, arXiv:1507.02091 [hep-ph] .[49]
CMS
Collaboration, S. Chatrchyan et al. , “Study of the inclusive production of charged pions,kaons, and protons in pp collisions at √ s = .
9, 2.76, and 7 TeV,”
Eur. Phys. J.
C72 (2012) 2164, arXiv:1207.4724 [hep-ex] .[50]
ALICE
Collaboration, S. Acharya et al. , “Production of charged pions, kaons and (anti-)protonsin Pb-Pb and inelastic pp collisions at √ s NN = 5.02 TeV,” arXiv:1910.07678 [nucl-ex] .21LICE Collaboration A The ALICE Collaboration
S. Acharya , D. Adamová , A. Adler , J. Adolfsson , M.M. Aggarwal , G. Aglieri Rinella ,M. Agnello , N. Agrawal
10 ,53 , Z. Ahammed , S. Ahmad , S.U. Ahn , A. Akindinov , M. Al-Turany ,S.N. Alam , D.S.D. Albuquerque , D. Aleksandrov , B. Alessandro , H.M. Alfanda , R. AlfaroMolina , B. Ali , Y. Ali , A. Alici
10 ,26 ,53 , A. Alkin , J. Alme , T. Alt , L. Altenkamper ,I. Altsybeev , M.N. Anaam , C. Andrei , D. Andreou , H.A. Andrews , A. Andronic , M. Angeletti ,V. Anguelov , C. Anson , T. Antiˇci´c , F. Antinori , P. Antonioli , N. Apadula , L. Aphecetche ,H. Appelshäuser , S. Arcelli , R. Arnaldi , M. Arratia , I.C. Arsene , M. Arslandok , A. Augustinus ,R. Averbeck , S. Aziz , M.D. Azmi , A. Badalà , Y.W. Baek , S. Bagnasco , X. Bai ,R. Bailhache , R. Bala , A. Balbino , A. Baldisseri , M. Ball , S. Balouza , D. Banerjee ,R. Barbera , L. Barioglio , G.G. Barnaföldi , L.S. Barnby , V. Barret , P. Bartalini , K. Barth ,E. Bartsch , F. Baruffaldi , N. Bastid , S. Basu , G. Batigne , B. Batyunya , D. Bauri ,J.L. Bazo Alba , I.G. Bearden , C. Beattie , C. Bedda , N.K. Behera , I. Belikov , A.D.C. BellHechavarria , F. Bellini , R. Bellwied , V. Belyaev , G. Bencedi , S. Beole , A. Bercuci ,Y. Berdnikov , D. Berenyi , R.A. Bertens , D. Berzano , M.G. Besoiu , L. Betev , A. Bhasin ,I.R. Bhat , M.A. Bhat , H. Bhatt , B. Bhattacharjee , A. Bianchi , L. Bianchi , N. Bianchi ,J. Bielˇcík , J. Bielˇcíková , A. Bilandzic
104 ,117 , G. Biro , R. Biswas , S. Biswas , J.T. Blair , D. Blau ,C. Blume , G. Boca , F. Bock
33 ,95 , A. Bogdanov , S. Boi , L. Boldizsár , A. Bolozdynya ,M. Bombara , G. Bonomi , H. Borel , A. Borissov , H. Bossi , E. Botta , L. Bratrud ,P. Braun-Munzinger , M. Bregant , M. Broz , E. Bruna , G.E. Bruno , M.D. Buckland ,D. Budnikov , H. Buesching , S. Bufalino , O. Bugnon , P. Buhler , P. Buncic , Z. Buthelezi
72 ,131 ,J.B. Butt , J.T. Buxton , S.A. Bysiak , D. Caffarri , A. Caliva , E. Calvo Villar , R.S. Camacho ,P. Camerini , A.A. Capon , F. Carnesecchi
10 ,26 , R. Caron , J. Castillo Castellanos , A.J. Castro ,E.A.R. Casula , F. Catalano , C. Ceballos Sanchez , P. Chakraborty , S. Chandra , W. Chang ,S. Chapeland , M. Chartier , S. Chattopadhyay , S. Chattopadhyay , A. Chauvin , C. Cheshkov ,B. Cheynis , V. Chibante Barroso , D.D. Chinellato , S. Cho , P. Chochula , T. Chowdhury ,P. Christakoglou , C.H. Christensen , P. Christiansen , T. Chujo , C. Cicalo , L. Cifarelli
10 ,26 ,F. Cindolo , G. Clai
53 ,ii , J. Cleymans , F. Colamaria , D. Colella , A. Collu , M. Colocci ,M. Concas
58 ,iii , G. Conesa Balbastre , Z. Conesa del Valle , G. Contin
24 ,59 , J.G. Contreras ,T.M. Cormier , Y. Corrales Morales , P. Cortese , M.R. Cosentino , F. Costa , S. Costanza ,P. Crochet , E. Cuautle , P. Cui , L. Cunqueiro , D. Dabrowski , T. Dahms
104 ,117 , A. Dainese ,F.P.A. Damas
114 ,137 , M.C. Danisch , A. Danu , D. Das , I. Das , P. Das , P. Das , S. Das ,A. Dash , S. Dash , S. De , A. De Caro , G. de Cataldo , J. de Cuveland , A. De Falco , D. DeGruttola , N. De Marco , S. De Pasquale , S. Deb , H.F. Degenhardt , K.R. Deja , A. Deloff ,S. Delsanto
25 ,131 , W. Deng , D. Devetak , P. Dhankher , D. Di Bari , A. Di Mauro , R.A. Diaz ,T. Dietel , P. Dillenseger , Y. Ding , R. Divià , D.U. Dixit , Ø. Djuvsland , U. Dmitrieva ,A. Dobrin , B. Dönigus , O. Dordic , A.K. Dubey , A. Dubla , S. Dudi , M. Dukhishyam ,P. Dupieux , R.J. Ehlers
95 ,146 , V.N. Eikeland , D. Elia , E. Epple , B. Erazmus , F. Erhardt ,A. Erokhin , M.R. Ersdal , B. Espagnon , G. Eulisse , D. Evans , S. Evdokimov , L. Fabbietti
104 ,117 ,M. Faggin , J. Faivre , F. Fan , A. Fantoni , M. Fasel , P. Fecchio , A. Feliciello , G. Feofilov ,A. Fernández Téllez , A. Ferrero , A. Ferretti , A. Festanti , V.J.G. Feuillard , J. Figiel ,S. Filchagin , D. Finogeev , F.M. Fionda , G. Fiorenza , F. Flor , S. Foertsch , P. Foka , S. Fokin ,E. Fragiacomo , U. Frankenfeld , U. Fuchs , C. Furget , A. Furs , M. Fusco Girard , J.J. Gaardhøje ,M. Gagliardi , A.M. Gago , A. Gal , C.D. Galvan , P. Ganoti , C. Garabatos , E. Garcia-Solis ,K. Garg , C. Gargiulo , A. Garibli , K. Garner , P. Gasik
104 ,117 , E.F. Gauger , M.B. Gay Ducati ,M. Germain , J. Ghosh , P. Ghosh , S.K. Ghosh , M. Giacalone , P. Gianotti , P. Giubellino
58 ,106 ,P. Giubilato , P. Glässel , A. Gomez Ramirez , V. Gonzalez , L.H. González-Trueba , S. Gorbunov ,L. Görlich , A. Goswami , S. Gotovac , V. Grabski , L.K. Graczykowski , K.L. Graham ,L. Greiner , A. Grelli , C. Grigoras , V. Grigoriev , A. Grigoryan , S. Grigoryan , O.S. Groettvik ,F. Grosa , J.F. Grosse-Oetringhaus , R. Grosso , R. Guernane , M. Guittiere , K. Gulbrandsen ,T. Gunji , A. Gupta , R. Gupta , I.B. Guzman , R. Haake , M.K. Habib , C. Hadjidakis ,H. Hamagaki , G. Hamar , M. Hamid , R. Hannigan , M.R. Haque
63 ,85 , A. Harlenderova ,J.W. Harris , A. Harton , J.A. Hasenbichler , H. Hassan , D. Hatzifotiadou
10 ,53 , P. Hauer ,S. Hayashi , S.T. Heckel
68 ,104 , E. Hellbär , H. Helstrup , A. Herghelegiu , T. Herman ,E.G. Hernandez , G. Herrera Corral , F. Herrmann , K.F. Hetland , H. Hillemanns , C. Hills ,B. Hippolyte , B. Hohlweger , J. Honermann , D. Horak , A. Hornung , S. Hornung , R. Hosokawa , P. Hristov , C. Huang , C. Hughes , P. Huhn , T.J. Humanic , H. Hushnud ,L.A. Husova , N. Hussain , S.A. Hussain , D. Hutter , J.P. Iddon
33 ,127 , R. Ilkaev , H. Ilyas ,M. Inaba , G.M. Innocenti , M. Ippolitov , A. Isakov , M.S. Islam , M. Ivanov , V. Ivanov ,V. Izucheev , B. Jacak , N. Jacazio , P.M. Jacobs , S. Jadlovska , J. Jadlovsky , S. Jaelani ,C. Jahnke , M.J. Jakubowska , M.A. Janik , T. Janson , M. Jercic , O. Jevons , M. Jin ,F. Jonas
95 ,144 , P.G. Jones , J. Jung , M. Jung , A. Jusko , P. Kalinak , A. Kalweit , V. Kaplin ,S. Kar , A. Karasu Uysal , O. Karavichev , T. Karavicheva , P. Karczmarczyk , E. Karpechev ,U. Kebschull , R. Keidel , M. Keil , B. Ketzer , Z. Khabanova , A.M. Khan , S. Khan , S.A. Khan ,A. Khanzadeev , Y. Kharlov , A. Khatun , A. Khuntia , B. Kileng , B. Kim , B. Kim , D. Kim ,D.J. Kim , E.J. Kim , H. Kim
17 ,147 , J. Kim , J.S. Kim , J. Kim , J. Kim , J. Kim , M. Kim ,S. Kim , T. Kim , T. Kim , S. Kirsch
38 ,68 , I. Kisel , S. Kiselev , A. Kisiel , J.L. Klay , C. Klein ,J. Klein
33 ,58 , S. Klein , C. Klein-Bösing , M. Kleiner , A. Kluge , M.L. Knichel , A.G. Knospe ,C. Kobdaj , M.K. Köhler , T. Kollegger , A. Kondratyev , N. Kondratyeva , E. Kondratyuk ,J. Konig , P.J. Konopka , L. Koska , O. Kovalenko , V. Kovalenko , M. Kowalski , I. Králik ,A. Kravˇcáková , L. Kreis , M. Krivda
64 ,110 , F. Krizek , K. Krizkova Gajdosova , M. Krüger ,E. Kryshen , M. Krzewicki , A.M. Kubera , V. Kuˇcera
33 ,60 , C. Kuhn , P.G. Kuijer , L. Kumar ,S. Kundu , P. Kurashvili , A. Kurepin , A.B. Kurepin , A. Kuryakin , S. Kushpil , J. Kvapil ,M.J. Kweon , J.Y. Kwon , Y. Kwon , S.L. La Pointe , P. La Rocca , Y.S. Lai , R. Langoy ,K. Lapidus , A. Lardeux , P. Larionov , E. Laudi , R. Lavicka , T. Lazareva , R. Lea , L. Leardini ,J. Lee , S. Lee , F. Lehas , S. Lehner , J. Lehrbach , R.C. Lemmon , I. León Monzón ,E.D. Lesser , M. Lettrich , P. Lévai , X. Li , X.L. Li , J. Lien , R. Lietava , B. Lim ,V. Lindenstruth , S.W. Lindsay , C. Lippmann , M.A. Lisa , A. Liu , J. Liu , S. Liu , W.J. Llope ,I.M. Lofnes , V. Loginov , C. Loizides , P. Loncar , J.A. Lopez , X. Lopez , E. López Torres ,J.R. Luhder , M. Lunardon , G. Luparello , Y.G. Ma , A. Maevskaya , M. Mager , S.M. Mahmood ,T. Mahmoud , A. Maire , R.D. Majka , M. Malaev , Q.W. Malik , L. Malinina
75 ,iv , D. Mal’Kevich ,P. Malzacher , G. Mandaglio , V. Manko , F. Manso , V. Manzari , Y. Mao , M. Marchisone ,J. Mareš , G.V. Margagliotti , A. Margotti , J. Margutti , A. Marín , C. Markert , M. Marquard ,C.D. Martin , N.A. Martin , P. Martinengo , J.L. Martinez , M.I. Martínez , G. Martínez García ,S. Masciocchi , M. Masera , A. Masoni , L. Massacrier , E. Masson , A. Mastroserio
52 ,138 ,A.M. Mathis
104 ,117 , O. Matonoha , P.F.T. Matuoka , A. Matyja , C. Mayer , F. Mazzaschi ,M. Mazzilli , M.A. Mazzoni , A.F. Mechler , F. Meddi , Y. Melikyan
62 ,92 , A. Menchaca-Rocha ,C. Mengke , E. Meninno
29 ,113 , M. Meres , S. Mhlanga , Y. Miake , L. Micheletti , D.L. Mihaylov ,K. Mikhaylov
75 ,91 , A.N. Mishra , D. Mi´skowiec , A. Modak , N. Mohammadi , A.P. Mohanty ,B. Mohanty , M. Mohisin Khan
16 ,v , Z. Moravcova , C. Mordasini , D.A. Moreira De Godoy ,L.A.P. Moreno , I. Morozov , A. Morsch , T. Mrnjavac , V. Muccifora , E. Mudnic , D. Mühlheim ,S. Muhuri , J.D. Mulligan , M.G. Munhoz , R.H. Munzer , H. Murakami , S. Murray , L. Musa ,J. Musinsky , C.J. Myers , J.W. Myrcha , B. Naik , R. Nair , B.K. Nandi , R. Nania
10 ,53 , E. Nappi ,M.U. Naru , A.F. Nassirpour , C. Nattrass , R. Nayak , T.K. Nayak , S. Nazarenko , A. Neagu ,R.A. Negrao De Oliveira , L. Nellen , S.V. Nesbo , G. Neskovic , D. Nesterov , L.T. Neumann ,B.S. Nielsen , S. Nikolaev , S. Nikulin , V. Nikulin , F. Noferini
10 ,53 , P. Nomokonov , J. Norman
78 ,127 ,N. Novitzky , P. Nowakowski , A. Nyanin , J. Nystrand , M. Ogino , A. Ohlson
80 ,103 , J. Oleniacz ,A.C. Oliveira Da Silva , M.H. Oliver , C. Oppedisano , A. Ortiz Velasquez , A. Oskarsson ,J. Otwinowski , K. Oyama , Y. Pachmayer , V. Pacik , D. Pagano , G. Pai´c , J. Pan ,S. Panebianco , P. Pareek
49 ,141 , J. Park , J.E. Parkkila , S. Parmar , S.P. Pathak , B. Paul , H. Pei ,T. Peitzmann , X. Peng , L.G. Pereira , H. Pereira Da Costa , D. Peresunko , G.M. Perez , Y. Pestov ,V. Petráˇcek , M. Petrovici , R.P. Pezzi , S. Piano , M. Pikna , P. Pillot , O. Pinazza
33 ,53 , L. Pinsky ,C. Pinto , S. Pisano
10 ,51 , D. Pistone , M. Płosko´n , M. Planinic , F. Pliquett , S. Pochybova
145 ,i ,M.G. Poghosyan , B. Polichtchouk , N. Poljak , A. Pop , S. Porteboeuf-Houssais , V. Pozdniakov ,S.K. Prasad , R. Preghenella , F. Prino , C.A. Pruneau , I. Pshenichnov , M. Puccio , J. Putschke ,L. Quaglia , R.E. Quishpe , S. Ragoni , S. Raha , S. Rajput , J. Rak , A. Rakotozafindrabe ,L. Ramello , F. Rami , S.A.R. Ramirez , R. Raniwala , S. Raniwala , S.S. Räsänen , R. Rath ,V. Ratza , I. Ravasenga , K.F. Read
95 ,130 , A.R. Redelbach , K. Redlich
84 ,vi , A. Rehman , P. Reichelt ,F. Reidt , X. Ren , R. Renfordt , Z. Rescakova , K. Reygers , V. Riabov , T. Richert
80 ,88 , M. Richter ,P. Riedler , W. Riegler , F. Riggi , C. Ristea , S.P. Rode , M. Rodríguez Cahuantzi , K. Røed ,R. Rogalev , E. Rogochaya , D. Rohr , D. Röhrich , P.S. Rokita , F. Ronchetti , E.D. Rosas ,K. Roslon , A. Rossi
28 ,56 , A. Rotondi , A. Roy , P. Roy , O.V. Rueda , R. Rui , B. Rumyantsev , A. Rustamov , E. Ryabinkin , Y. Ryabov , A. Rybicki , H. Rytkonen , O.A.M. Saarimaki ,S. Sadhu , S. Sadovsky , K. Šafaˇrík , S.K. Saha , B. Sahoo , P. Sahoo , R. Sahoo , S. Sahoo ,P.K. Sahu , J. Saini , S. Sakai , S. Sambyal , V. Samsonov
92 ,97 , D. Sarkar , N. Sarkar ,P. Sarma , V.M. Sarti , M.H.P. Sas , E. Scapparone , J. Schambach , H.S. Scheid , C. Schiaua ,R. Schicker , A. Schmah , C. Schmidt , H.R. Schmidt , M.O. Schmidt , M. Schmidt ,N.V. Schmidt
68 ,95 , A.R. Schmier , J. Schukraft , Y. Schutz
33 ,136 , K. Schwarz , K. Schweda ,G. Scioli , E. Scomparin , M. Šefˇcík , J.E. Seger , Y. Sekiguchi , D. Sekihata , I. Selyuzhenkov
92 ,106 ,S. Senyukov , D. Serebryakov , E. Serradilla , A. Sevcenco , A. Shabanov , A. Shabetai ,R. Shahoyan , W. Shaikh , A. Shangaraev , A. Sharma , A. Sharma , H. Sharma , M. Sharma ,N. Sharma , S. Sharma , A.I. Sheikh , K. Shigaki , M. Shimomura , S. Shirinkin , Q. Shou ,Y. Sibiriak , S. Siddhanta , T. Siemiarczuk , D. Silvermyr , G. Simatovic , G. Simonetti , B. Singh ,R. Singh , R. Singh , R. Singh , V.K. Singh , V. Singhal , T. Sinha , B. Sitar , M. Sitta ,T.B. Skaali , M. Slupecki , N. Smirnov , R.J.M. Snellings , C. Soncco , J. Song ,A. Songmoolnak , F. Soramel , S. Sorensen , I. Sputowska , J. Stachel , I. Stan , P. Stankus ,P.J. Steffanic , E. Stenlund , D. Stocco , M.M. Storetvedt , L.D. Stritto , A.A.P. Suaide ,T. Sugitate , C. Suire , M. Suleymanov , M. Suljic , R. Sultanov , M. Šumbera , V. Sumberia ,S. Sumowidagdo , S. Swain , A. Szabo , I. Szarka , U. Tabassam , S.F. Taghavi , G. Taillepied ,J. Takahashi , G.J. Tambave , S. Tang , M. Tarhini , M.G. Tarzila , A. Tauro , G. Tejeda Muñoz ,A. Telesca , L. Terlizzi , C. Terrevoli , D. Thakur , S. Thakur , D. Thomas , F. Thoresen ,R. Tieulent , A. Tikhonov , A.R. Timmins , A. Toia , N. Topilskaya , M. Toppi , F. Torales-Acosta ,S.R. Torres
36 ,120 , A. Trifiro , S. Tripathy
49 ,69 , T. Tripathy , S. Trogolo , G. Trombetta , L. Tropp ,V. Trubnikov , W.H. Trzaska , T.P. Trzcinski , B.A. Trzeciak
36 ,63 , T. Tsuji , A. Tumkin , R. Turrisi ,T.S. Tveter , K. Ullaland , E.N. Umaka , A. Uras , G.L. Usai , M. Vala , N. Valle , S. Vallero ,N. van der Kolk , L.V.R. van Doremalen , M. van Leeuwen , P. Vande Vyvre , D. Varga , Z. Varga ,M. Varga-Kofarago , A. Vargas , M. Vasileiou , A. Vasiliev , O. Vázquez Doce
104 ,117 , V. Vechernin ,E. Vercellin , S. Vergara Limón , L. Vermunt , R. Vernet , R. Vértesi , L. Vickovic , Z. Vilakazi ,O. Villalobos Baillie , G. Vino , A. Vinogradov , T. Virgili , V. Vislavicius , A. Vodopyanov ,B. Volkel , M.A. Völkl , K. Voloshin , S.A. Voloshin , G. Volpe , B. von Haller , I. Vorobyev ,D. Voscek , J. Vrláková , B. Wagner , M. Weber , A. Wegrzynek , S.C. Wenzel , J.P. Wessels ,J. Wiechula , J. Wikne , G. Wilk , J. Wilkinson
10 ,53 , G.A. Willems , E. Willsher , B. Windelband ,M. Winn , W.E. Witt , Y. Wu , R. Xu , S. Yalcin , Y. Yamaguchi , K. Yamakawa , S. Yang ,S. Yano , Z. Yin , H. Yokoyama , I.-K. Yoo , J.H. Yoon , S. Yuan , A. Yuncu , V. Yurchenko ,V. Zaccolo , A. Zaman , C. Zampolli , H.J.C. Zanoli , N. Zardoshti , A. Zarochentsev , P. Závada ,N. Zaviyalov , H. Zbroszczyk , M. Zhalov , S. Zhang , X. Zhang , Z. Zhang , V. Zherebchevskii ,D. Zhou , Y. Zhou , Z. Zhou , J. Zhu , Y. Zhu , A. Zichichi
10 ,26 , G. Zinovjev , N. Zurlo , Affiliation notes i Deceased ii Italian National Agency for New Technologies, Energy and Sustainable Economic Development (ENEA),Bologna, Italy iii
Dipartimento DET del Politecnico di Torino, Turin, Italy iv M.V. Lomonosov Moscow State University, D.V. Skobeltsyn Institute of Nuclear, Physics, Moscow, Russia v Department of Applied Physics, Aligarh Muslim University, Aligarh, India vi Institute of Theoretical Physics, University of Wroclaw, Poland
Collaboration Institutes A.I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation, Yerevan, Armenia Bogolyubov Institute for Theoretical Physics, National Academy of Sciences of Ukraine, Kiev, Ukraine Bose Institute, Department of Physics and Centre for Astroparticle Physics and Space Science (CAPSS),Kolkata, India Budker Institute for Nuclear Physics, Novosibirsk, Russia California Polytechnic State University, San Luis Obispo, California, United States Central China Normal University, Wuhan, China Centre de Calcul de l’IN2P3, Villeurbanne, Lyon, France Centro de Aplicaciones Tecnológicas y Desarrollo Nuclear (CEADEN), Havana, Cuba Centro de Investigación y de Estudios Avanzados (CINVESTAV), Mexico City and Mérida, Mexico Centro Fermi - Museo Storico della Fisica e Centro Studi e Ricerche “Enrico Fermi’, Rome, Italy Chicago State University, Chicago, Illinois, United States China Institute of Atomic Energy, Beijing, China Comenius University Bratislava, Faculty of Mathematics, Physics and Informatics, Bratislava, Slovakia COMSATS University Islamabad, Islamabad, Pakistan Creighton University, Omaha, Nebraska, United States Department of Physics, Aligarh Muslim University, Aligarh, India Department of Physics, Pusan National University, Pusan, Republic of Korea Department of Physics, Sejong University, Seoul, Republic of Korea Department of Physics, University of California, Berkeley, California, United States Department of Physics, University of Oslo, Oslo, Norway Department of Physics and Technology, University of Bergen, Bergen, Norway Dipartimento di Fisica dell’Università ’La Sapienza’ and Sezione INFN, Rome, Italy Dipartimento di Fisica dell’Università and Sezione INFN, Cagliari, Italy Dipartimento di Fisica dell’Università and Sezione INFN, Trieste, Italy Dipartimento di Fisica dell’Università and Sezione INFN, Turin, Italy Dipartimento di Fisica e Astronomia dell’Università and Sezione INFN, Bologna, Italy Dipartimento di Fisica e Astronomia dell’Università and Sezione INFN, Catania, Italy Dipartimento di Fisica e Astronomia dell’Università and Sezione INFN, Padova, Italy Dipartimento di Fisica ‘E.R. Caianiello’ dell’Università and Gruppo Collegato INFN, Salerno, Italy Dipartimento DISAT del Politecnico and Sezione INFN, Turin, Italy Dipartimento di Scienze e Innovazione Tecnologica dell’Università del Piemonte Orientale and INFNSezione di Torino, Alessandria, Italy Dipartimento Interateneo di Fisica ‘M. Merlin’ and Sezione INFN, Bari, Italy European Organization for Nuclear Research (CERN), Geneva, Switzerland Faculty of Electrical Engineering, Mechanical Engineering and Naval Architecture, University of Split,Split, Croatia Faculty of Engineering and Science, Western Norway University of Applied Sciences, Bergen, Norway Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague, Prague,Czech Republic Faculty of Science, P.J. Šafárik University, Košice, Slovakia Frankfurt Institute for Advanced Studies, Johann Wolfgang Goethe-Universität Frankfurt, Frankfurt,Germany Fudan University, Shanghai, China Gangneung-Wonju National University, Gangneung, Republic of Korea Gauhati University, Department of Physics, Guwahati, India Helmholtz-Institut für Strahlen- und Kernphysik, Rheinische Friedrich-Wilhelms-Universität Bonn, Bonn,Germany Helsinki Institute of Physics (HIP), Helsinki, Finland High Energy Physics Group, Universidad Autónoma de Puebla, Puebla, Mexico Hiroshima University, Hiroshima, Japan Hochschule Worms, Zentrum für Technologietransfer und Telekommunikation (ZTT), Worms, Germany Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest, Romania Indian Institute of Technology Bombay (IIT), Mumbai, India Indian Institute of Technology Indore, Indore, India Indonesian Institute of Sciences, Jakarta, Indonesia INFN, Laboratori Nazionali di Frascati, Frascati, Italy INFN, Sezione di Bari, Bari, Italy INFN, Sezione di Bologna, Bologna, Italy INFN, Sezione di Cagliari, Cagliari, Italy INFN, Sezione di Catania, Catania, Italy INFN, Sezione di Padova, Padova, Italy INFN, Sezione di Roma, Rome, Italy INFN, Sezione di Torino, Turin, Italy INFN, Sezione di Trieste, Trieste, Italy Inha University, Incheon, Republic of Korea Institut de Physique Nucléaire d’Orsay (IPNO), Institut National de Physique Nucléaire et de Physique desParticules (IN2P3/CNRS), Université de Paris-Sud, Université Paris-Saclay, Orsay, France Institute for Nuclear Research, Academy of Sciences, Moscow, Russia Institute for Subatomic Physics, Utrecht University/Nikhef, Utrecht, Netherlands Institute of Experimental Physics, Slovak Academy of Sciences, Košice, Slovakia Institute of Physics, Homi Bhabha National Institute, Bhubaneswar, India Institute of Physics of the Czech Academy of Sciences, Prague, Czech Republic Institute of Space Science (ISS), Bucharest, Romania Institut für Kernphysik, Johann Wolfgang Goethe-Universität Frankfurt, Frankfurt, Germany Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, Mexico City, Mexico Instituto de Física, Universidade Federal do Rio Grande do Sul (UFRGS), Porto Alegre, Brazil Instituto de Física, Universidad Nacional Autónoma de México, Mexico City, Mexico iThemba LABS, National Research Foundation, Somerset West, South Africa Jeonbuk National University, Jeonju, Republic of Korea Johann-Wolfgang-Goethe Universität Frankfurt Institut für Informatik, Fachbereich Informatik undMathematik, Frankfurt, Germany Joint Institute for Nuclear Research (JINR), Dubna, Russia Korea Institute of Science and Technology Information, Daejeon, Republic of Korea KTO Karatay University, Konya, Turkey Laboratoire de Physique Subatomique et de Cosmologie, Université Grenoble-Alpes, CNRS-IN2P3,Grenoble, France Lawrence Berkeley National Laboratory, Berkeley, California, United States Lund University Department of Physics, Division of Particle Physics, Lund, Sweden Nagasaki Institute of Applied Science, Nagasaki, Japan Nara Women’s University (NWU), Nara, Japan National and Kapodistrian University of Athens, School of Science, Department of Physics , Athens,Greece National Centre for Nuclear Research, Warsaw, Poland National Institute of Science Education and Research, Homi Bhabha National Institute, Jatni, India National Nuclear Research Center, Baku, Azerbaijan National Research Centre Kurchatov Institute, Moscow, Russia Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark Nikhef, National institute for subatomic physics, Amsterdam, Netherlands NRC Kurchatov Institute IHEP, Protvino, Russia NRC ´nKurchatov InstituteÂ˙z - ITEP, Moscow, Russia NRNU Moscow Engineering Physics Institute, Moscow, Russia Nuclear Physics Group, STFC Daresbury Laboratory, Daresbury, United Kingdom Nuclear Physics Institute of the Czech Academy of Sciences, ˇRež u Prahy, Czech Republic Oak Ridge National Laboratory, Oak Ridge, Tennessee, United States Ohio State University, Columbus, Ohio, United States Petersburg Nuclear Physics Institute, Gatchina, Russia Physics department, Faculty of science, University of Zagreb, Zagreb, Croatia Physics Department, Panjab University, Chandigarh, India
Physics Department, University of Jammu, Jammu, India
Physics Department, University of Rajasthan, Jaipur, India
Physikalisches Institut, Eberhard-Karls-Universität Tübingen, Tübingen, Germany
Physikalisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany
Physik Department, Technische Universität München, Munich, Germany
Politecnico di Bari, Bari, Italy
Research Division and ExtreMe Matter Institute EMMI, GSI Helmholtzzentrum fürSchwerionenforschung GmbH, Darmstadt, Germany
Rudjer Boškovi´c Institute, Zagreb, Croatia
Russian Federal Nuclear Center (VNIIEF), Sarov, Russia
Saha Institute of Nuclear Physics, Homi Bhabha National Institute, Kolkata, India
School of Physics and Astronomy, University of Birmingham, Birmingham, United Kingdom
Sección Física, Departamento de Ciencias, Pontificia Universidad Católica del Perú, Lima, Peru
St. Petersburg State University, St. Petersburg, Russia
Stefan Meyer Institut für Subatomare Physik (SMI), Vienna, Austria
SUBATECH, IMT Atlantique, Université de Nantes, CNRS-IN2P3, Nantes, France
Suranaree University of Technology, Nakhon Ratchasima, Thailand
Technical University of Košice, Košice, Slovakia
Technische Universität München, Excellence Cluster ’Universe’, Munich, Germany
The Henryk Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, Cracow, Poland
The University of Texas at Austin, Austin, Texas, United States
Universidad Autónoma de Sinaloa, Culiacán, Mexico
Universidade de São Paulo (USP), São Paulo, Brazil
Universidade Estadual de Campinas (UNICAMP), Campinas, Brazil
Universidade Federal do ABC, Santo Andre, Brazil
University of Cape Town, Cape Town, South Africa
University of Houston, Houston, Texas, United States
University of Jyväskylä, Jyväskylä, Finland
University of Liverpool, Liverpool, United Kingdom
University of Science and Technology of China, Hefei, China
University of South-Eastern Norway, Tonsberg, Norway
University of Tennessee, Knoxville, Tennessee, United States
University of the Witwatersrand, Johannesburg, South Africa
University of Tokyo, Tokyo, Japan
University of Tsukuba, Tsukuba, Japan
Université Clermont Auvergne, CNRS/IN2P3, LPC, Clermont-Ferrand, France
Université de Lyon, Université Lyon 1, CNRS/IN2P3, IPN-Lyon, Villeurbanne, Lyon, France
Université de Strasbourg, CNRS, IPHC UMR 7178, F-67000 Strasbourg, France, Strasbourg, France
Université Paris-Saclay Centre d’Etudes de Saclay (CEA), IRFU, Départment de Physique Nucléaire(DPhN), Saclay, France
Università degli Studi di Foggia, Foggia, Italy
Università degli Studi di Pavia, Pavia, Italy
Università di Brescia, Brescia, Italy
Variable Energy Cyclotron Centre, Homi Bhabha National Institute, Kolkata, India
Warsaw University of Technology, Warsaw, Poland
Wayne State University, Detroit, Michigan, United States
Westfälische Wilhelms-Universität Münster, Institut für Kernphysik, Münster, Germany
Wigner Research Centre for Physics, Budapest, Hungary
Yale University, New Haven, Connecticut, United States