J/ ψ production as a function of charged-particle multiplicity in p-Pb collisions at s NN − − − √ = 8.16 TeV
aa r X i v : . [ nu c l - e x ] A p r EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH
CERN-EP-2020-05815 April 2020c (cid:13) J / ψ production as a function of charged-particle multiplicity in p–Pbcollisions at √ s NN = .
16 TeV
ALICE Collaboration ∗ Abstract
Inclusive J / ψ yields and average transverse momenta in p–Pb collisions at a center-of-mass energyper nucleon pair √ s NN = .
16 TeV are measured as a function of the charged-particle pseudora-pidity density with ALICE. The J / ψ mesons are reconstructed at forward (2 . < y cms < .
53) andbackward ( − . < y cms < − .
96) center-of-mass rapidity in their dimuon decay channel while thecharged-particle pseudorapidity density is measured around midrapidity. The J / ψ yields at forwardand backward rapidity normalized to their respective average values increase with the normalizedcharged-particle pseudorapidity density, the former showing a weaker increase than the latter. Thenormalized average transverse momenta at forward and backward rapidity manifest a steady increasefrom low to high charged-particle pseudorapidity density with a saturation beyond the average value. ∗ See Appendix A for the list of collaboration members ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration
Quarkonium states have long been considered as probes of the Quark–Gluon Plasma (QGP) produced inultra-relativistic heavy-ion collisions [1]. The large color-charge density in the plasma prevents the for-mation of bound states, in an analogous process to the Debye screening for electromagnetic processes [2].The suppression of J / ψ production in nucleus–nucleus (AA) with respect to proton–proton (pp) colli-sions was observed by several experiments [3–11]. To determine whether the origin of this suppression isthe influence of the QGP or of Cold Nuclear Matter (CNM), data on proton(deuteron)–nucleus collisionsare also scrutinized.The measurements in p–Pb collisions at the LHC show a suppression of J / ψ production [12–14] atlow transverse momentum ( p T ) and forward center-of-mass rapidity (p-going direction, positive y cms ),consistent with various combinations of CNM effects: modification of the parton distribution functions(PDFs) in nuclei, i.e. shadowing [15, 16], the Color-Glass Condensate (CGC) [17, 18], or coherentparton energy loss [19]. The measurement of ψ ( ) production in p–Pb collisions [20] exhibits a largersuppression than the one measured for J / ψ , both at forward and backward rapidity, which was notexpected from CNM predictions. This effect is reproduced by models which consider the break-up ofthe bound quark–anti-quark pair via interactions with the final-state comoving particles [21, 22].The p–Pb data at the center-of-mass energy per nucleon–nucleon collision of √ s NN = .
02 TeV [23, 24]showed that these effects depend on the centrality of the collision, as estimated from the energy depositedin the Zero Degree Calorimeter in the Pb-going direction [25], and/or the produced charged-particlemultiplicity [26]. An increase of the relative J / ψ and ϒ [26–28] yields with the relative charged-particlemultiplicity is observed, similarly to the results from pp collisions [27–29]. The increase of the J / ψ normalized yields was observed to be similar to the increase for D mesons [30, 31], suggesting thata common mechanism may be at its origin. The excited-to-ground state ratios, ϒ ( nS ) / ϒ ( ) , werefound to decrease with increasing charged-particle multiplicity, which was not expected from CNMpredictions [27, 28].The measurements of two-particle angular correlations in small systems have shown interesting struc-tures in the angular correlation function. A near-side ridge, located at ( ∆ϕ ) ≈
0, is observed in high-multiplicity pp [32] and p–Pb [33] collisions, accompanied by an away-side structure, located at ∆ϕ ≈ π and exceeding the away-side jet contribution, in p–Pb collisions [34, 35]. These structures are remi-niscent of those in Pb–Pb data [36], interpreted as signatures of the collective motion of the particlesduring the hydrodynamic evolution of the hot and dense medium. Correlations of J / ψ (at large rapidity)and charged particles (at midrapidity) in p–Pb collisions [37, 38] revealed persisting long-range corre-lation structures at high p T , similar to those observed with charged hadrons. The corresponding ellipticflow coefficients are found to be positive and of comparable magnitude to those measured in Pb–Pbcollisions [39–41], indicating that the mechanism at its origin could be similar in both collision systems.This letter reports the measurement of the multiplicity-differential inclusive J / ψ yield and average trans-verse momentum in p–Pb collisions at √ s NN = .
16 TeV. The J / ψ mesons are reconstructed at for-ward and backward center-of-mass rapidities in their dimuon decay channel. The charged-particle pseu-dorapidity density is measured around midrapidity. It complements and extends previous J / ψ mea-surements performed as a function of the collision centrality and the charged-particle multiplicity at √ s NN = .
02 TeV [23, 26]. The classification of events as a function of their charged-particle pseudora-pidity density enables the scrutiny of rare events, corresponding to the 0.01–0.04% highest multiplicitiesin the collision. This allows p–Pb events to be studied from low multiplicities, similar to those of ppcollisions, up to very large multiplicities corresponding to ∼
100 produced charged particles per rapidityunit, similar to those of peripheral Pb–Pb collisions, which exhibit collective-like effects.2ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration
In this section, the detector subsystems relevant for this analysis are presented. A complete descriptionof the ALICE detector and its performance can be found in [42, 43].The muon spectrometer [42, 43] covers the pseudorapidity window of − . < η < − . λ int ), starting at90 cm from the nominal interaction point, ten layers of muon tracking chambers (MCH), coupled to adipole magnet with a 3 Tm field integral, and four layers of muon trigger chambers (MTR). The MCHand MTR systems are separated by an additional iron wall of about 7 . λ int that absorbs the remaininghadronic and low-momentum particle contamination. A rear absorber positioned downstream of theMTR filters out the background from beam-gas interactions. A conical absorber surrounds the beam pipeand protects the spectrometer against secondary particles produced mainly by large- η primary particlesinteracting with the beam pipe.The Silicon Pixel Detector (SPD) [44] is the innermost part of the Inner Tracking System (ITS). Itconsists of two cylindrical silicon pixel layers at radial distances of 3 . . | η | < | η | < .
4. The SPD is used toreconstruct the primary vertex and to measure the charged-particle pseudorapidity density at midrapidity.The V0 scintillator arrays [45] are located at each side of the interaction point, covering the pseudora-pidity ranges of − . < η < − . . < η < .
1. In this analysis, the V0 provides an online triggerand helps to reject contamination from beam-gas events.The neutron Zero Degree Calorimeter (ZDC) [42] located at about 112.5 m on either side from theinteraction point are used to reject electromagnetic interactions and beam-induced background.The results presented in this letter are obtained with data recorded during the p–Pb run at √ s NN = .
16 TeVin 2016. The J / ψ are reconstructed in the dimuon channel with data taken in two different beam con-figurations. Due to the asymmetry of the beam energy per nucleon in p–Pb collisions at the LHC, thenucleon–nucleon center-of-mass rapidity frame is shifted by ∆ y = .
465 in the direction of the protonbeam. As a consequence, the J / ψ are measured in the forward rapidity range of 2 . < y cms < .
53 (withprotons going in the direction of the muon spectrometer, p-going direction) and in the backward rapidityregion − . < y cms < − .
96 (Pb-going direction). Events used in this analysis were collected with adedicated dimuon trigger which requires the coincidence of signals in both V0 arrays (minimum biastrigger, MB) with at least two opposite-sign muons registered in the MTR. The trigger has an adjustableonline threshold, which for this data sample was set to only accept muons with transverse momenta p T > . c ( p T for which an efficiency of 50% is reached). The p T differential single-muon triggerefficiency reaches a plateau of ∼
96% at p T ∼ . c . In this data-taking period, the maximum pile-up probability was about 4%. A dedicated event-selection strategy—exploiting the signals in the V0 andthe ZDC, the correlation of the number of clusters and track segments reconstructed in the SPD, as wellas an algorithm to tag events with multiple vertices—allowed us to keep the pile-up below 0.5% for theanalysed events, even at large multiplicities. The charged-particle pseudorapidity density (d N ch / d η ) is measured at midrapidity exploiting the in-formation provided by the SPD detector [46, 47]. It is evaluated by counting the number of tracklets( N tracklet ), i.e. track segments joining pairs of hits in the two layers of the SPD pointing to the primaryvertex. The primary vertex is also computed with the SPD information. To minimize non-uniformities inthe SPD acceptance, only events with a z -vertex position determined within | z vtx | <
10 cm are considered,and tracklets are counted within | η | <
1. 3ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE CollaborationThe raw N tracklet counts are corrected ( N corrtracklet ) for the variation of the detector conditions with time (frac-tion of active SPD channels) and its limited acceptance as a function of z vtx using a data-driven event-by-event correction [29, 30]. In this analysis, the correction is done by renormalising the N tracklet ( z vtx ) distributions to the overall maximum with a Poissonian smearing to account for the fluctuations. Theevents are sliced in N corrtracklet intervals. Monte Carlo (MC) simulations using the DPMJET [48] event gen-erator and the GEANT3 transport code [49] are used to estimate d N ch / d η from N corrtracklet . A second orderpolynomial correlation is assumed between these two quantities for the full N corrtracklet interval. Severalsources of systematic uncertainty were taken into account. Possible deviations from the second orderpolynomial correlation were estimated by using other functions to quantify the correlation or MC aver-ages in each interval, with values ranging from 0.1% at intermediate multiplicities to 6.9% (5.8%) at thelowest (highest) multiplicity intervals. The systematic uncertainty on the residual z vtx dependence due todifferences between data and MC amounts to 3%. Finally, the event generator influence was consideredand evaluated by comparing the DPMJET simulations with events generated in EPOS [50], resulting ina 2% uncertainty.The average charged-particle pseudorapidity density, h d N ch / d η i , in non-single diffractive (NSD) eventswas obtained from an independent analysis and amounts to h d N ch / d η i = . ± .
83 (20 . ± .
83) inp–Pb (Pb–p) collisions for | η | < Table 1:
Sources of systematic uncertainties on the normalized charged-particle multiplicity.
Source | η | < N corrtracklet to d N ch / d η correlation 0.1 – 6.9(5.8)% z -vertex dependence 3%Monte Carlo event generator 2% h d N ch / d η i J / ψ measurement The normalized J / ψ yield, i.e. the yield in each multiplicity interval i normalized to the multiplicity-integrated value, is evaluated as d N i / d y h d N / d y i = N i J / ψ N J / ψ N eqMB N i , eqMB ( A ε ) J / ψ ( A ε ) i J / ψ ε i MB ε MB , (1)from the reconstructed number of J / ψ , N J / ψ , the number of minimum bias (MB) events equivalent tothe analysed dimuon sample, N eqMB , the J / ψ acceptance and efficiency correction, ( A ε ) J / ψ , and the NSDevent selection efficiency in the minimum bias sample, ε MB .The J / ψ are reconstructed for each multiplicity interval by combining opposite-sign muons and comput-ing the invariant mass of the pairs. The muon identification is ensured by requiring that the track candi-dates reconstructed in the MCH have a matching track segment in the MTR. Furthermore, the individualtracks must fulfill the following criteria to make sure they are within the acceptance of the spectrometer:their radial distance from the beam axis at the end of the front absorber is within 17 . < R abs < . − < η < − . / ψ acceptance timesefficiency ( A ε ), differentially in p T and y . The resulting distributions are then fitted with a superpo-sition of J / ψ and ψ ( ) signals and a background lineshape. Various combinations of lineshapes areused in order to evaluate the signal counts and their uncertainties. The two charmonium resonancesare parametrized by a sum of either two Crystal Ball or two pseudo-Gaussian functions with power-law4ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaborationtails [51]. The tail parametrizations are fixed to the values determined from either fits of the J / ψ signalfrom MC simulations or to values taken from fits to the multiplicity-integrated distribution in p–Pb dataat √ s NN = .
16 TeV [13] and in pp data at √ s =
13 TeV [52]. The tails obtained from fitting themultiplicity-integrated distributions using the Crystal Ball function are also considered, and fixed in thebinned fits. The J / ψ peak mean position and width are left free in the multiplicity-integrated fit, whilstthe ψ ( ) ones are bound to those of the J / ψ following the same procedure as in [53]. Note that the ψ ( ) yields obtained are not physical values, as the invariant-mass spectrum is corrected by the A ε correction for the J / ψ . In the multiplicity-differential fits, the mass and width of the J / ψ peak are fixedto the integrated values to ensure the convergence of the fits in the few cases where statistical significanceis low. The background is parameterized by either a sum of two exponentials or the product of an expo-nential and a fourth-order polynomial. Two fit mass ranges are taken into account when computing theaverage number of J / ψ and its uncertainty: 1 . < m µµ < . / c and 2 . < m µµ < . / c . Ex-amples of fits at low, intermediate, and high multiplicity for data in the rapidity range 2 . < y cms < . c (GeV/ µµ m c c oun t s pe r M e V / = 8.16 TeV NN s Pb, − ALICE, p < 3.53 cms y trackletscorr N ≤
1 423 ± = 18328 ψ J/ N c ± = 3095 ψ J/ µ c ± = 75 ψ J/ σ c (GeV/ µµ m c c oun t s pe r M e V / = 8.16 TeV NN s Pb, − ALICE, p < 3.53 cms y trackletscorr N ≤
49 762 ± = 27843 ψ J/ N c ± = 3096 ψ J/ µ c ± = 76 ψ J/ σ c (GeV/ µµ m c c oun t s pe r M e V / = 8.16 TeV NN s Pb, − ALICE, p < 3.53 cms y trackletscorr N ≤
111 520 ± = 14032 ψ J/ N c ± = 3096 ψ J/ µ c ± = 72 ψ J/ σ Figure 1:
Opposite-sign muon pair invariant mass distributions for selected multiplicity intervals, corrected forthe J / ψ acceptance and efficiency, at forward rapidity. The distributions are shown together with a typical fitfunction (solid line, see text for details). The J / ψ signal contribution is also depicted by a dot-dashed red line, andthe background by a dotted line. / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaborationdoes change with multiplicity. Therefore, in order to minimize the uncertainty on the signal extraction,the same signal lineshape is used in the fit function for both the numerator and denominator in Eq. 1.The number of equivalent MB events N eqMB is computed from the number of dimuon triggered events, N µµ ,and the normalization factor of dimuon triggered to MB events (calculated as explained in next section)as N eqMB = F norm · N µµ . The number needs to be corrected for by the NSD event selection efficiency, ε MB = ( ± ) % [47], to take into account the fraction of events without a reconstructed SPD vertexthat are rejected. This factor ε MB is found to be independent of the charged-particle multiplicity in allthe intervals studied, with the exception of the lowest multiplicity interval, where it decreases by 1%.The J / ψ acceptance and efficiency correction is obtained from MC simulations as a function of p T and y cms . The J / ψ are generated using p T and y cms distributions tuned to data [13]. They are simulated todecay into a muon pair using EvtGen [54]. The final state radiation is described with PHOTOS [55].The acceptance and efficiency correction is independent of multiplicity in the measurement intervals.Therefore, when estimating the uncertainty on the MC input, only the possible variation of the input p T and y cms distributions is taken into account by using as input a subsample of the lower/higher multiplicityevents.To extract the J / ψ mean transverse momentum h p J / ψ T i , the A ε -corrected transverse momentum of thedimuon pair is fitted with the following function [26]: h p µµ T i ( m µµ ) = α J / ψ ( m µµ ) h p J / ψ T i + α ψ ′ ( m µµ ) h p ψ ′ T i + (cid:16) − α J / ψ ( m µµ ) − α ψ ′ ( m µµ ) (cid:17) h p bkgdT i ( m µµ ) , (2)where the ratios of signal over the sum of signal and background of the two charmonium states α J / ψ = S J / ψ / ( S J / ψ + S ψ ′ + B ) and α ψ ′ = S ψ ′ / ( S J / ψ + S ψ ′ + B ) are fixed to the value extracted from fittingthe invariant-mass spectrum corrected by the J / ψ A ε . The background is described by a function h p bkgdT i ( m µµ ) . Two functional forms are used: either a sum of two exponentials or the product of anexponential and a fourth-order polynomial. Note that the h p ψ ′ T i does not represent a physical mean trans-verse momentum of the ψ ( ) as the spectra are corrected by the A ε for J / ψ . Figure 2 illustrates typical h p µµ T i distributions for selected multiplicity intervals. c (GeV/ µµ m ) c ( G e V / 〉 µµ T p 〈 = 8.16 TeV NN s Pb, − ALICE, p < 3.53 cms y c ± = 2.45 〉 ψ J/T p 〈 < 10 trackletscorr N ≤ c (GeV/ µµ m ) c ( G e V / 〉 µµ T p 〈 = 8.16 TeV NN s Pb, − ALICE, p < 3.53 cms y c ± = 2.79 〉 ψ J/T p 〈 < 155 trackletscorr N ≤ Figure 2:
Average transverse momentum of opposite-sign muon pairs for selected multiplicity intervals, correctedfor the J / ψ acceptance and efficiency. The distributions are shown together with a typical fit function (see text fordetails). / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration
The following sources of systematic uncertainty on the J / ψ yields in multiplicity classes are considered:(i) the signal extraction, (ii) the normalisation, (iii) the effect of resolution and pile-up, (iv) the even-t-by-event N tracklet to N corrtracklet correction, and (v) the event selection efficiency of the NSD event class.For the measurement of the yields in each multiplicity interval normalized to the event average, the sys-tematic uncertainties are estimated directly for this ratio. Details on the signal extraction uncertaintywere addressed in the previous section. The values are estimated by varying the signal and backgroundshapes of the fit function, as well as by varying the invariant-mass range of the fit. The systematic un-certainty is computed as the root-mean-square of the uncertainties on the ratio for each of these fits,ranging between 0.8–2.3% (0.5–1.9%) at forward (backward) rapidity, being larger at large multiplic-ities where the number of events is smaller. The normalisation factor of the dimuon triggered to MBevents F norm is studied using three alternative methods [13]. The first method evaluates the probabilityof a coincidence of a dimuon- and a MB-triggered event in a MB-triggered data set. The second methodexploits the higher probability of occurrence of a single-muon trigger by looking at the product of theprobability of coincidence of a single-muon- and MB-triggered event and of the probability of findinga dimuon event in the single-muon triggered data. The third method is based on information from thetrigger scalers. The run-by-run spread of the F norm / F i norm values computed for these three methods de-termines a 2.5% systematic uncertainty, independent of multiplicity. The effect of the method of choicefor the event-by-event correction from N tracklet to N corrtracklet on the J / ψ yield is also studied [29, 30]. Boththe randomisation function (Poisson or binomial) and the reference normalisation of the correction arevaried. The Poissonian smearing is applied when the maximum is selected as normalisation reference,while the binomial correction should be used when considering all other possible reference values (in ourcase the minimum). The influence of these modifications on the yield ranges from 0.1% to 2.6% (4.3%)at forward (backward) rapidity, as a function of multiplicity. The uncertainty coming from pile up andmultiplicity axis resolution is estimated as a single contribution by repeating the analysis multiple timeswith a different randomisation seed for the event-by-event correction, or introducing a small shift of the N corrtracklet intervals, or varying the pile-up rejection criteria. The uncertainty amounts to 2%, independentof multiplicity. The uncertainty on the event selection efficiency for the NSD event class is estimatedas in Ref. [47]. The uncertainty amounts to 1% and is correlated in all multiplicity intervals. Table 2summarizes all contributions to the systematic uncertainty on the normalized yield.For the h p T i , the effects of the uncertainty on the h p T i extraction procedure and of the A ε are considered.Similar to the yields, the signal extraction uncertainty is estimated by varying the fit function and itsrange. In addition, as the S / ( S + B ) terms in Eq. 2 are fixed in the fit to the h p T i invariant-mass spectrum,the influence of the statistical uncertainty on the J / ψ signal S is introduced via a Gaussian smearing of S (with respect to its statistical uncertainty) to prevent artificially minimising the uncertainty. It ranges from0.2% to 3.0% (1.2%) at forward (backward) rapidity, increasing with multiplicity as a consequence of thesmaller number of events. The uncertainty on the absolute h p T i also takes into account the uncertaintyon: (i) the MC input shapes as a function of p T and y cms , ranging from < . < . N corrtracklet distribution for each rapidity interval. For each of these bins, the h p T i is estimated using amodified A ε correction, which was re-weighted to better describe the p T - and y -dependent distributionsof J / ψ in given bin. The systematic uncertainty is taken as the difference of the original h p T i value,computed with the initial A ε correction, and the new h p T i estimated with re-weighted correction. Theuncertainty on all the measured multiplicity intervals is extrapolated from these two values assuming thatin each class the uncertainty is proportional to the h p T i . The contributions of the tracking, the triggerand their matching to the uncertainty are correlated between multiplicity intervals. The normalized h p T i values are only affected by the uncertainty on the signal extraction procedure and the MC input,7ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaborationwhich is partly correlated in multiplicity and ranges from < .
1% to 2% ( < .
1% to 4%) in the forward(backward) rapidity interval. Table 3 summarizes all contributions to the average and normalized average p T measurements. Table 2:
Sources of systematic uncertainties on the normalized yield. The contributions marked with an asteriskare correlated in multiplicity.
Source 2 . < y cms < . − . < y cms < − . . .
3% 0 . . F norm ) 2 .
5% 2 . N corrtracklet . .
6% 0 . . ∗ ∗ Table 3:
Systematic uncertainty sources on the average and normalized average p T . The values in parenthesescorrespond to the multiplicity-integrated uncertainties related to the signal extraction. The contributions markedwith an asterisk are correlated in multiplicity. The uncertainty on MC input, marked with a diamond, is partiallycorrelated in multiplicity. . < y cms < . − . < y cms < − . h p T i h p T i (cid:14) h p intT i h p T i h p T i (cid:14) h p intT i Signal extraction 0.2–3.0% (0.2%) 0.3–3.0% 0.2–1.2% (0.2%) 0.3–1.3%Tracking efficiency 1%* – 1%* –Trigger efficiency 2.6%* – 3.1%* –Track–trigger matching 1%* – 1%* –Monte Carlo input < . − ⋄ < . − ⋄ < . − ⋄ < . − ⋄ The normalized J / ψ yield, at forward and backward rapidities, is presented in Fig. 3 as a function ofthe normalized charged-particle pseudorapidity density, measured at midrapidity ( | η | < x regime ( x Pb ∼ − in a naive 2-body calculation for p T = x Pb ∼ − ). The observed suppression of the p T - andmultiplicity-integrated J / ψ yield at forward rapidity is described by different cold nuclear matter modelsconsidering the probed shadowing/saturation domain [13]. The centrality-differential measurements at √ s NN = .
02 TeV [23] of the nuclear modification factor, h p T i and h p i , corresponding to relativemultiplicities of at most 2.5 times the average one, can also be described by these models.The normalized J / ψ yield is compared with the one measured in p–Pb collisions at √ s NN = .
02 TeV [26]in Fig. 4. Good agreement is found for both rapidity intervals. These results extend the probed charged-particle pseudorapidity density interval, both at low and high multiplicity, examining events of up toalmost six times the average value. The more precise √ s NN = .
16 TeV data evidence a continuousincrease with multiplicity up to the largest multiplicities attained. The similarities at √ s NN = .
16 TeVand √ s NN = .
02 TeV suggest a common origin of the multiplicity trend, with a mechanism whoseeffect varies with rapidity, but might have a small dependence on the collision energy. This is consistentwith the large variation of the probed x Pb with rapidity and its relative slow evolution on the collisionenergy (typically a factor of 2 in the simplified 2-body picture).8ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration |<1 η | NSD 〉η / d ch N d 〈 η / d ch N d02468 N S D 〉 y / d N d 〈 y / d N d = 8.16 TeV NN s Pb, − ALICE, p - µ + µ → ψ J/ < 3.53 (p-going) cms y − < cms y −
1% norm. unc. not shown ± Figure 3:
Normalized yield of inclusive J / ψ , at forward and backward rapidities, as a function of the normalizedcharged-particle pseudorapidity density, measured at midrapidity, in p–Pb collisions at √ s NN = .
16 TeV. Thevertical bars represent the statistical uncertainties. The vertical and horizontal widths of the boxes represent therespective systematic uncertainties for the J / ψ yields and the multiplicities. The dashed line indicates the diagonalline, to guide the eye. / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration |<1 η | NSD 〉η / d ch N d 〈 η / d ch N d02468 N S D 〉 y / d N d 〈 y / d N d Pb − ALICE, p < 3.53 (p-going) cms y - µ + µ → ψ J/ = 8.16 TeV NN s = 5.02 TeV NN s ±
1% norm. unc. not shown at 8.16 TeV ± |<1 η | NSD 〉η / d ch N d 〈 η / d ch N d02468 N S D 〉 y / d N d 〈 y / d N d Pb − ALICE, p − < cms y − - µ + µ → ψ J/ = 8.16 TeV NN s = 5.02 TeV NN s ±
1% norm. unc. not shown at 8.16 TeV ± Figure 4:
Normalized yield of inclusive J / ψ as a function of the normalized charged-particle pseudorapiditydensity, measured at midrapidity, in p–Pb collisions at √ s NN = .
16 TeV and √ s NN = .
02 TeV [26]. The top(bottom) panel presents the measurement at forward (backward) rapidity. The vertical bars represent the statisticaluncertainties, the boxes the systematic ones. The dashed line indicates the diagonal line, to guide the eye. / ψ production in p–Pb at √ s NN = .
16 TeV ALICE CollaborationFigure 5 presents a comparison of the normalized J / ψ p–Pb yields with results from pp collisions at √ s NN = ( . < y cms < . ) and Pb–Pb collisions at √ s NN = .
02 TeV [56] ( . < y cms < . ) .The ratio of the yields over the corresponding charged-particle multiplicity is also shown in Fig. 5.The trend exhibited by the pp data is similar to the one observed in the backward (Pb-going) direction. |<1 η | 〉η / d ch N d 〈 η / d ch N d02468 〉 y / d N d 〈 y / d N d ALICE - µ + µ → ψ J/ = 8.16 TeV NN s Pb, − p < 3.53 (p-going) cms y − < cms y − < 4.0 cms y s pp, = 5.02 TeV NN s Pb, − Pb ±
1% norm. unc. not shown at 8.16 TeV ± |<1 η | 〉η / d ch N d 〈 η / d ch N d00.511.52 〉 η / d c h N d 〈 η / d c h N d / 〉 y / d N d 〈 y / d N d ALICE - µ + µ → ψ J/ = 8.16 TeV NN s Pb, − p < 3.53 (p-going) cms y − < cms y − < 4.0 cms y s pp, = 5.02 TeV NN s Pb, − Pb ±
1% norm. unc. not shown at 8.16 TeV ± Figure 5:
Top: Normalized yield of inclusive J / ψ as a function of the normalized charged-particle pseudorapid-ity density, measured at midrapidity, in various collision systems. Bottom: Ratio of the normalized yields to thecorresponding normalized charged-particle pseudorapidity density. The pp results are normalized to INEL colli-sions [29], whereas p–Pb ones refer to the NSD event class; all for p T >
0. The Pb–Pb data points include J / ψ with 0 . < p T <
12 GeV / c to reduce the low- p T contribution from photoproduction, which is significant only inmore peripheral collisions [56–58]. The vertical bars represent the statistical uncertainties, the boxes the systematicones. The dashed line indicates the diagonal line, to guide the eye. It should be noted that the pp results are normalized to the inelastic ‘INEL’ event class, whereas thep–Pb measurements are normalized to the non-single-diffractive ‘NSD’ one. In p–Pb collisions thesetwo event classes mostly overlap when comparing MB results [47]. The Pb–Pb data also show a faster-than-linear increase with the normalized charged-particle pseudorapidity density. They are compatiblewithin uncertainties with the p–Pb backward rapidity result in the restricted multiplicity interval of themeasurement. Whereas the pp and p–Pb data include J / ψ with p T >
0, the Pb–Pb data points include J / ψ / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaborationwith 0 . < p T <
12 GeV / c to reduce the low- p T contribution from photoproduction, which is significantonly in more peripheral collisions [56–58].The measured yield in p–Pb collisions can be described with the EPOS 3 event generator [59, 60] (seeFig. 6) based on a combination of Gribov-Regge theory and pQCD: where the individual scatteringsare identified with parton ladders emerging as flux tubes, the existence of multiple nucleon–nucleoncollisions in pPb collisions is accounted for, the initial conditions of the collision are modified due toCNM effects including parton saturation, and slow string segments (far from the surface) can be furthermapped to fluid dynamic fields using a core-corona description. The J / ψ bound-state formation inEPOS 3 assumes a color-evaporation approach, i.e. it is associated to a charm quark–anti-quark pairin a given mass range. The influence of the viscous hydrodynamic evolution of the bulk in the EPOS 3calculation is small (see Fig. 6). However the number of simulated events at large multiplicities is limitedand does not allow us to elucidate possible hydrodynamic effects. |<1 η | NSD 〉η / d ch N d 〈 η / d ch N d02468 N S D 〉 y / d N d 〈 y / d N d = 8.16 TeV NN s Pb, − ALICE, p < 3.53 (p-going) cms y - µ + µ → ψ J/ dataEPOS without hydroEPOS with hydro
1% norm. unc. not shown ± |<1 η | NSD 〉η / d ch N d 〈 η / d ch N d02468 N S D 〉 y / d N d 〈 y / d N d = 8.16 TeV NN s Pb, − ALICE, p − < cms y − - µ + µ → ψ J/ dataEPOS without hydroEPOS with hydro
1% norm. unc. not shown ± Figure 6:
Normalized yield of inclusive J / ψ as a function of the normalized charged-particle pseudorapiditydensity, measured at midrapidity, in p–Pb collisions at √ s NN = .
16 TeV compared with EPOS 3 [59, 60] cal-culations. Left (right) panel presents the measurement at forward (backward) rapidity. The vertical bars representthe statistical uncertainties, the boxes the systematic ones. The dashed line indicates the diagonal line, to guide theeye. The shaded areas represent the statistical uncertainties on the EPOS 3 calculations.
Figure 7 presents h p T i as a function of the relative charged-particle pseudorapidity density. The mea-sured h p T i is systematically smaller at backward than at forward rapidity. This is also true for themultiplicity-integrated value, which is consistent with the observed decrease of h p T i with increasing | y cms | in pp collisions [61]. The h p T i increases steadily for multiplicities below the average, and saturatesabove the average multiplicity. Two naive scenarios are typically considered to explain high-multiplicityevents: the incoherent superposition of multiple parton–parton collisions, or single parton interactionswith higher energy transfer. One would expect the latter to be characterized by a higher h p T i of the pro-duced J / ψ . Reality is probably somewhere in between these two simplified scenarios. The simultaneousincrease of the yield together with the saturation of h p T i may point to J / ψ production from an incoherentsuperposition of parton–parton collisions, as suggested by two-particle correlation studies [62].A comparison of the normalized h p T i at √ s NN = .
16 TeV and √ s NN = .
02 TeV [26] as a functionof the normalized charged-particle pseudorapidity density is shown in Fig. 8. The measurements arein remarkable agreement, within the uncertainties, at both energies and rapidities confirming that themechanism governing J / ψ production and its variation with the charged-particle pseudorapidity densityhas small or no variation with the collision energy, in the explored interval.12ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration |<1 η | NSD 〉η / d ch N d 〈 η / d ch N d1.522.533.5 ) c ( G e V / N S D 〉 T p 〈 = 8.16 TeV NN s Pb, − ALICE, p - µ + µ → ψ J/ < 3.53 (p-going) cms y − < cms y − for the p(Pb)-going 3.0 (3.4)% norm. unc. not shown ± Figure 7:
Average transverse momentum of inclusive J / ψ at forward and backward rapidities as a func-tion of the normalized charged-particle pseudorapidity density, measured at midrapidity, in p–Pb collisions at √ s NN = .
16 TeV. The vertical bars represent the statistical uncertainties, the boxes the systematic ones.
The production of inclusive J / ψ at large rapidities in p–Pb collisions at √ s NN = .
16 TeV is reported asa function of the charged-particle pseudorapidity density at midrapidity. The normalised J / ψ yield showsan increase with increasing normalised charged-particle pseudorapidity density. The yield at backwardrapidity grows faster than the forward rapidity one, reaching values above those of the linear (with slopeunity) increase estimate at large normalised multiplicity, whereas the values at forward rapidity show aslower-than-linear increase. The trends of the normalised yield are reproduced by the EPOS 3 [59, 60]event generator. The h p T i is smaller at backward than at forward rapidity, consistent with the expectedsoftening of the spectra with increasing | y cms | . The h p T i increases steadily for multiplicities below theaverage, and saturates above the average multiplicity. The simultaneous increase of the yield togetherwith the saturation of h p T i may point to J / ψ production from an incoherent superposition of parton–parton collisions, as suggested by two-particle correlation studies [62]. These measurements show trendscompatible with those observed at √ s NN = .
02 TeV [26], but in this work an improved precisionand extended multiplicity coverage were reached. The similarities suggest a common origin, with amechanism whose effect varies with rapidity, but with only a small dependence (if any) on the collisionenergy. 13ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration |<1 η | NSD 〉η / d ch N d 〈 η / d ch N d0.70.80.911.11.2 N S D 〉 i n t T p 〈 / 〉 T p 〈 Pb − ALICE, p < 3.53 (p-going) cms y - µ + µ → ψ J/ = 8.16 TeV NN s = 5.02 TeV NN s |<1 η | NSD 〉η / d ch N d 〈 η / d ch N d0.70.80.911.11.2 N S D 〉 i n t T p 〈 / 〉 T p 〈 Pb − ALICE, p − < cms y − - µ + µ → ψ J/ = 8.16 TeV NN s = 5.02 TeV NN s Figure 8:
Normalized average transverse momentum of inclusive J / ψ as a function of the normalizedcharged-particle pseudorapidity density, measured at midrapidity, in p–Pb collisions at √ s NN = .
16 TeV and √ s NN = .
02 TeV [26]. Top (bottom) panel presents the measurement at forward (backward) rapidity. Thevertical bars represent the statistical uncertainties, the boxes the systematic ones. / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration
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)and 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) and Institut National de Physique Nucléaire et de Physiquedes Particules (IN2P3) and Centre National de la Recherche Scientifique (CNRS), France; Bundesmin-isterium für Bildung und Forschung (BMBF) and GSI Helmholtzzentrum für SchwerionenforschungGmbH, Germany; General Secretariat for Research and Technology, Ministry of Education, Researchand Religions, Greece; National Research, Development and Innovation Office, Hungary; Departmentof Atomic Energy Government of India (DAE), Department of Science and Technology, Governmentof India (DST), University Grants Commission, Government of India (UGC) and Council of Scientificand Industrial Research (CSIR), India; Indonesian Institute of Science, Indonesia; Centro Fermi - MuseoStorico della Fisica e Centro Studi e Ricerche Enrico Fermi and Istituto Nazionale di Fisica Nucleare(INFN), Italy; Institute for Innovative Science and Technology , Nagasaki Institute of Applied Science(IIST), Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT) and Japan So-ciety for the Promotion of Science (JSPS) KAKENHI, Japan; Consejo Nacional de Ciencia (CONACYT)y Tecnología, through Fondo de Cooperación Internacional en Ciencia y Tecnología (FONCICYT) andDirección General de Asuntos del Personal Academico (DGAPA), Mexico; Nederlandse Organisatievoor Wetenschappelijk Onderzoek (NWO), Netherlands; The Research Council of Norway, Norway;Commission on Science and Technology for Sustainable Development in the South (COMSATS), Pak-istan; Pontificia Universidad Católica del Perú, Peru; Ministry of Science and Higher Education, NationalScience Centre and WUT ID-UB, Poland; Korea Institute of Science and Technology Information andNational Research Foundation of Korea (NRF), Republic of Korea; Ministry of Education and ScientificResearch, Institute of Atomic Physics and Ministry of Research and Innovation and Institute of AtomicPhysics, Romania; Joint Institute for Nuclear Research (JINR), Ministry of Education and Science ofthe Russian Federation, National Research Centre Kurchatov Institute, Russian Science Foundation andRussian Foundation for Basic Research, Russia; Ministry of Education, Science, Research and Sport ofthe Slovak Republic, Slovakia; National Research Foundation of South Africa, South Africa; SwedishResearch Council (VR) and Knut & Alice Wallenberg Foundation (KAW), Sweden; European Organi-zation for Nuclear Research, Switzerland; Suranaree University of Technology (SUT), National Scienceand Technology Development Agency (NSDTA) and Office of the Higher Education Commission underNRU project of Thailand, Thailand; Turkish Atomic Energy Agency (TAEK), Turkey; National Academyof Sciences of Ukraine, Ukraine; Science and Technology Facilities Council (STFC), United Kingdom;National Science Foundation of the United States of America (NSF) and United States Department of15ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE CollaborationEnergy, Office of Nuclear Physics (DOE NP), United States of America.
References [1] P. Braun-Munzinger and J. Stachel, “The quest for the quark-gluon plasma”,
Nature (2007) 302–309.[2] T. Matsui and H. Satz, “J / ψ Suppression by Quark-Gluon Plasma Formation”,
Phys. Lett.
B178 (1986) 416–422.[3]
NA50
Collaboration, B. Alessandro et al. , “A New measurement of J / ψ suppression in Pb-Pbcollisions at 158-GeV per nucleon”, Eur. Phys. J.
C39 (2005) 335–345, arXiv:hep-ex/0412036 [hep-ex] .[4]
NA60
Collaboration, R. Arnaldi et al. , “J / ψ production in indium-indium collisions at158-GeV/nucleon”, Phys. Rev. Lett. (2007) 132302.[5] PHENIX
Collaboration, A. Adare et al. , “ J / ψ suppression at forward rapidity in Au+Aucollisions at √ s NN =
200 GeV”,
Phys. Rev.
C84 (2011) 054912, arXiv:1103.6269 [nucl-ex] .[6]
ALICE
Collaboration, B. Abelev et al. , “ J / ψ suppression at forward rapidity in Pb-Pb collisionsat √ s NN = .
76 TeV”,
Phys. Rev. Lett. (2012) 072301, arXiv:1202.1383 [hep-ex] .[7]
ALICE
Collaboration, B. Abelev et al. , “Centrality, rapidity and transverse momentumdependence of J / ψ suppression in Pb-Pb collisions at √ s NN = 2.76 TeV”, Phys. Lett.
B734 (2014) 314–327, arXiv:1311.0214 [nucl-ex] .[8]
STAR
Collaboration, L. Adamczyk et al. , “ J / ψ production at low p T in Au + Au and Cu + Cucollisions at √ s NN =
200 GeV with the STAR detector”,
Phys. Rev.
C90 no. 2, (2014) 024906, arXiv:1310.3563 [nucl-ex] .[9]
CMS
Collaboration, S. Chatrchyan et al. , “Suppression of non-prompt J / ψ , prompt J / ψ , andY(1S) in PbPb collisions at √ s NN = .
76 TeV”,
JHEP (2012) 063, arXiv:1201.5069 [nucl-ex] .[10] ATLAS
Collaboration, G. Aad et al. , “Measurement of the centrality dependence of J / ψ yieldsand observation of Z production in lead-lead collisions with the ATLAS detector at the LHC”, Phys. Lett.
B697 (2011) 294–312, arXiv:1012.5419 [hep-ex] .[11]
ALICE
Collaboration, J. Adam et al. , “J/ ψ suppression at forward rapidity in Pb-Pb collisions at √ s NN = .
02 TeV”,
Phys. Lett.
B766 (2017) 212–224, arXiv:1606.08197 [nucl-ex] .[12]
ALICE
Collaboration, B. Abelev et al. , “ J / ψ production and nuclear effects in p-Pb collisions at √ s NN = 5.02 TeV”, JHEP (2014) 073, arXiv:1308.6726 [nucl-ex] .[13] ALICE
Collaboration, S. Acharya et al. , “Inclusive J/ ψ production at forward and backwardrapidity in p-Pb collisions at √ s NN = 8.16 TeV”, JHEP (2018) 160, arXiv:1805.04381 [nucl-ex] .[14] LHCb
Collaboration, R. Aaij et al. , “Study of J / ψ production and cold nuclear matter effects in pPb collisions at √ s NN = JHEP (2014) 072, arXiv:1308.6729 [nucl-ex] .[15] J. L. Albacete et al. , “Predictions for p + Pb Collisions at √ s NN = Int. J. Mod. Phys.
E25 no. 9, (2016) 1630005, arXiv:1605.09479 [hep-ph] .16ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration[16] A. Kusina, J.-P. Lansberg, I. Schienbein, and H.-S. Shao, “Gluon Shadowing in Heavy-FlavorProduction at the LHC”,
Phys. Rev. Lett. no. 5, (2018) 052004, arXiv:1712.07024 [hep-ph] .[17] Y.-Q. Ma, R. Venugopalan, K. Watanabe, and H.-F. Zhang, “ ψ ( S ) versus J / ψ suppression inproton-nucleus collisions from factorization violating soft color exchanges”, Phys. Rev.
C97 no. 1, (2018) 014909, arXiv:1707.07266 [hep-ph] .[18] B. Ducloué, T. Lappi, and H. Mäntysaari, “Forward J / ψ and D meson nuclear suppression at theLHC”, Nucl. Part. Phys. Proc. (2017) 309–312, arXiv:1612.04585 [hep-ph] .[19] F. Arleo and S. Peigné, “Quarkonium suppression in heavy-ion collisions from coherent energyloss in cold nuclear matter”,
JHEP (2014) 073, arXiv:1407.5054 [hep-ph] .[20] ALICE
Collaboration, B. Abelev et al. , “Suppression of ψ (2S) production in p-Pb collisions at √ s NN = 5.02 TeV”, JHEP (2014) 073, arXiv:1405.3796 [nucl-ex] .[21] E. G. Ferreiro, “Excited charmonium suppression in proton-nucleus collisions as a consequence ofcomovers”, Phys. Lett.
B749 (2015) 98–103, arXiv:1411.0549 [hep-ph] .[22] B. Chen, T. Guo, Y. Liu, and P. Zhuang, “Cold and Hot Nuclear Matter Effects on CharmoniumProduction in p+Pb Collisions at LHC Energy”,
Phys. Lett.
B765 (2017) 323–327, arXiv:1607.07927 [nucl-th] .[23]
ALICE
Collaboration, J. Adam et al. , “Centrality dependence of inclusive J / ψ production in p-Pbcollisions at √ s NN = .
02 TeV”,
JHEP (2015) 127, arXiv:1506.08808 [nucl-ex] .[24] ALICE
Collaboration, J. Adam et al. , “Centrality dependence of ψ (2S) suppression in p-Pbcollisions at √ s NN = 5.02 TeV”, JHEP (2016) 050, arXiv:1603.02816 [nucl-ex] .[25] ALICE
Collaboration, J. Adam et al. , “Centrality dependence of particle production in p-Pbcollisions at √ s NN = 5.02 TeV”, Phys. Rev.
C91 no. 6, (2015) 064905, arXiv:1412.6828 [nucl-ex] .[26]
ALICE
Collaboration, D. Adamová et al. , “J/ ψ production as a function of charged-particlepseudorapidity density in p-Pb collisions at √ s NN = .
02 TeV”,
Phys. Lett.
B776 (2018) 91–104, arXiv:1704.00274 [nucl-ex] .[27]
CMS
Collaboration, S. Chatrchyan et al. , “Event activity dependence of Y(nS) production in √ s NN =5.02 TeV pPb and √ s =2.76 TeV pp collisions”, JHEP (2014) 103, arXiv:1312.6300 [nucl-ex] .[28] ATLAS
Collaboration, M. Aaboud et al. , “Measurement of quarkonium production inprotonâ ˘A ¸Slead and protonâ ˘A ¸Sproton collisions at 5 .
02 TeV with the ATLAS detector”,
Eur. Phys. J.
C78 no. 3, (2018) 171, arXiv:1709.03089 [nucl-ex] .[29]
ALICE
Collaboration, B. Abelev et al. , “ J / ψ production as a function of charged particlemultiplicity in pp collisions at √ s = Phys. Lett.
B712 (2012) 165–175, arXiv:1202.2816 [hep-ex] .[30]
ALICE
Collaboration, J. Adam et al. , “Measurement of charm and beauty production at centralrapidity versus charged-particle multiplicity in proton-proton collisions at √ s = JHEP (2015) 148, arXiv:1505.00664 [nucl-ex] .[31] ALICE
Collaboration, J. Adam et al. , “Measurement of D-meson production versus multiplicityin p-Pb collisions at √ s NN = .
02 TeV”,
JHEP (2016) 078, arXiv:1602.07240 [nucl-ex] .17ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration[32]
CMS
Collaboration, V. Khachatryan et al. , “Observation of Long-Range Near-Side AngularCorrelations in Proton-Proton Collisions at the LHC”,
JHEP (2010) 091, arXiv:1009.4122 [hep-ex] .[33] 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] .[34]
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] .[35]
ATLAS
Collaboration, G. Aad et al. , “Observation of Associated Near-Side and Away-SideLong-Range Correlations in √ s NN = 5.02 TeV Proton-Lead Collisions with the ATLAS Detector”, Phys. Rev. Lett. no. 18, (2013) 182302, arXiv:1212.5198 [hep-ex] .[36]
ALICE
Collaboration, K. Aamodt et al. , “Harmonic decomposition of two-particle angularcorrelations in Pb-Pb collisions at √ s NN = Phys. Lett.
B708 (2012) 249–264, arXiv:1109.2501 [nucl-ex] .[37]
ALICE
Collaboration, S. Acharya et al. , “Search for collectivity with azimuthal J/ ψ -hadroncorrelations in high multiplicity p-Pb collisions at √ s NN = 5.02 and 8.16 TeV”, Phys. Lett.
B780 (2018) 7–20, arXiv:1709.06807 [nucl-ex] .[38]
CMS
Collaboration, A. M. Sirunyan et al. , “Observation of prompt J/ ψ meson elliptic flow inhigh-multiplicity pPb collisions at √ s NN = Phys. Lett.
B791 (2019) 172–194, arXiv:1810.01473 [hep-ex] .[39]
ALICE
Collaboration, S. Acharya et al. , “J/ ψ elliptic flow in Pb-Pb collisions at √ s NN = . Phys. Rev. Lett. no. 24, (2017) 242301, arXiv:1709.05260 [nucl-ex] .[40]
ALICE
Collaboration, S. Acharya et al. , “Study of J/ ψ azimuthal anisotropy at forward rapidityin Pb-Pb collisions at √ s NN = .
02 TeV”,
JHEP (2019) 012, arXiv:1811.12727 [nucl-ex] .[41] ATLAS
Collaboration, M. Aaboud et al. , “Prompt and non-prompt J / ψ elliptic flow in Pb+Pbcollisions at √ s NN = .
02 Tev with the ATLAS detector”,
Eur. Phys. J.
C78 no. 9, (2018) 784, arXiv:1807.05198 [nucl-ex] .[42]
ALICE
Collaboration, K. Aamodt et al. , “The ALICE experiment at the CERN LHC”,
JINST (2008) S08002.[43] 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] .[44]
ALICE
Collaboration, K. Aamodt et al. , “Alignment of the ALICE Inner Tracking System withcosmic-ray tracks”,
JINST (2010) P03003, arXiv:1001.0502 [physics.ins-det] .[45] ALICE
Collaboration, E. Abbas et al. , “Performance of the ALICE VZERO system”,
JINST (2013) P10016, arXiv:1306.3130 [nucl-ex] .[46] ALICE
Collaboration, B. Abelev et al. , “Pseudorapidity density of charged particles in p + Pbcollisions at √ s NN = .
02 TeV”,
Phys. Rev. Lett. no. 3, (2013) 032301, arXiv:1210.3615 [nucl-ex] . 18ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration[47]
ALICE
Collaboration, S. Acharya et al. , “Charged-particle pseudorapidity density at mid-rapidityin p-Pb collisions at √ s NN = 8.16 TeV”, Eur. Phys. J.
C79 no. 4, (2019) 307, arXiv:1812.01312 [nucl-ex] .[48] S. Roesler, R. Engel, and J. Ranft, “The Monte Carlo event generator DPMJET-III”, in
AdvancedMonte Carlo for radiation physics, particle transport simulation and applications. Proceedings,Conference, MC2000, Lisbon, Portugal, October 23-26, 2000 , pp. 1033–1038. 2000. arXiv:hep-ph/0012252 [hep-ph] .[49] R. Brun, F. Bruyant, F. Carminati, S. Giani, M. Maire, A. McPherson, G. Patrick, and L. Urban,
GEANT: Detector Description and Simulation Tool; Oct 1994 . CERN Program Library. CERN,Geneva, 1993. https://cds.cern.ch/record/1082634 . Long Writeup W5013.[50] T. Pierog, I. Karpenko, J. M. Katzy, E. Yatsenko, and K. Werner, “EPOS LHC: Test of collectivehadronization with data measured at the CERN Large Hadron Collider”,
Phys. Rev.
C92 no. 3, (2015) 034906, arXiv:1306.0121 [hep-ph] .[51]
ALICE
Collaboration, J. Adam et al. , “Quarkonium signal extraction in ALICE”,ALICE-PUBLIC-2015-006, Oct, 2015.[52]
ALICE
Collaboration, S. Acharya et al. , “Energy dependence of forward-rapidity J / ψ and ψ ( ) production in pp collisions at the LHC”, Eur. Phys. J.
C77 no. 6, (2017) 392, arXiv:1702.00557 [hep-ex] .[53]
ALICE
Collaboration, J. Adam et al. , “Inclusive quarkonium production at forward rapidity in ppcollisions at √ s = Eur. Phys. J.
C76 no. 4, (2016) 184, arXiv:1509.08258 [hep-ex] .[54] D. J. Lange, “The EvtGen particle decay simulation package”,
Nucl. Instrum. Meth.
A462 (2001) 152–155.[55] E. Barberio and Z. Was, “PHOTOS: A Universal Monte Carlo for QED radiative corrections.Version 2.0”,
Comput. Phys. Commun. (1994) 291–308.[56] ALICE
Collaboration, S. Acharya et al. , “Studies of J/ ψ production at forward rapidity in Pb-Pbcollisions at √ s NN = 5.02 TeV”, JHEP (2020) 041, arXiv:1909.03158 [nucl-ex] .[57] ALICE
Collaboration, J. Adam et al. , “Measurement of an excess in the yield of J / ψ at very low p T in Pb-Pb collisions at √ s NN = 2.76 TeV”, Phys. Rev. Lett. no. 22, (2016) 222301, arXiv:1509.08802 [nucl-ex] .[58]
ALICE
Collaboration, S. Acharya et al. , “Coherent J/ ψ photoproduction at forward rapidity inultra-peripheral Pb-Pb collisions at √ s NN = .
02 TeV”,
Phys. Lett.
B798 (2019) 134926, arXiv:1904.06272 [nucl-ex] .[59] H. J. Drescher, M. Hladik, S. Ostapchenko, T. Pierog, and K. Werner, “Parton based Gribov-Reggetheory”,
Phys. Rept. (2001) 93–289, arXiv:hep-ph/0007198 [hep-ph] .[60] K. Werner, B. Guiot, I. Karpenko, and T. Pierog, “Analysing radial flow features in p-Pb and p-pcollisions at several TeV by studying identified particle production in EPOS3”,
Phys. Rev.
C89 no. 6, (2014) 064903, arXiv:1312.1233 [nucl-th] .[61]
LHCb
Collaboration, R. Aaij et al. , “Measurement of J / ψ production in pp collisions at √ s = Eur. Phys. J.
C71 (2011) 1645, arXiv:1103.0423 [hep-ex] .[62]
ALICE
Collaboration, B. Abelev et al. , “Multiplicity dependence of jet-like two-particlecorrelation structures in p–Pb collisions at √ s NN = 5.02 TeV”, Phys. Lett.
B741 (2015) 38–50, arXiv:1406.5463 [nucl-ex] . 19ultiplicity dependent J / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration
A The ALICE Collaboration
S. Acharya , D. Adamová , A. Adler , J. Adolfsson , M.M. Aggarwal , G. Aglieri Rinella ,M. Agnello , N. Agrawal
10 ,54 , Z. Ahammed , S. Ahmad , S.U. Ahn , Z. Akbar , A. Akindinov ,M. Al-Turany , S.N. Alam , D.S.D. Albuquerque , D. Aleksandrov , B. Alessandro , H.M. Alfanda ,R. Alfaro Molina , B. Ali , Y. Ali , A. Alici
10 ,26 ,54 , 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 , G. Biro , R. Biswas , S. Biswas , J.T. Blair , D. Blau ,C. Blume , G. Boca , F. Bock , A. Bogdanov , S. Boi , J. Bok , 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 , S.A. Bysiak , D. Caffarri , A. Caliva , E. Calvo Villar , R.S. Camacho , P. Camerini ,A.A. Capon , F. Carnesecchi , 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
54 ,ii ,J. Cleymans , F. Colamaria , D. Colella , A. Collu , M. Colocci , M. Concas
59 ,iii , G. ConesaBalbastre , Z. Conesa del Valle , G. Contin
24 ,60 , J.G. Contreras , T.M. Cormier , Y. Corrales Morales ,P. Cortese , M.R. Cosentino , F. Costa , S. Costanza , J. Crkovska , P. Crochet , E. Cuautle ,P. Cui , L. Cunqueiro , D. Dabrowski , T. Dahms , A. Dainese , F.P.A. Damas
115 ,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. De Gruttola , 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
90 ,107 , S. Dudi , M. Dukhishyam , P. Dupieux , R.J. Ehlers
96 ,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 , 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 ,A.N. Flores , 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
105 ,107 , E.F. Gauger , M.B. Gay Ducati , M. Germain , J. Ghosh , P. Ghosh ,S.K. Ghosh , M. Giacalone , P. Gianotti , P. Giubellino
59 ,107 , P. Giubilato , P. Glässel , A. GomezRamirez , V. Gonzalez
107 ,143 , 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 ,86 , A. Harlenderova , J.W. Harris , A. Harton , J.A. Hasenbichler ,H. Hassan , D. Hatzifotiadou
10 ,54 , P. Hauer , L.B. Havener , S. Hayashi , S.T. Heckel , 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 , / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration
P. Huhn , T.J. Humanic , H. Hushnud , L.A. Husova , N. Hussain , S.A. Hussain , D. Hutter ,J.P. Iddon
34 ,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
96 ,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 , J. Kim , J.S. Kim ,J. Kim , J. Kim , J. Kim , M. Kim , S. Kim , T. Kim , T. Kim , S. Kirsch , I. Kisel ,S. Kiselev , A. Kisiel , J.L. Klay , C. Klein , J. Klein
34 ,59 , 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 , S.A. Konigstorfer , P.J. Konopka ,G. Kornakov , L. Koska , O. Kovalenko , V. Kovalenko , M. Kowalski , I. Králik ,A. Kravˇcáková , L. Kreis , M. Krivda
64 ,111 , F. Krizek , K. Krizkova Gajdosova , M. Krüger ,E. Kryshen , M. Krzewicki , A.M. Kubera , V. Kuˇcera
34 ,61 , 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 , A. Lindner , 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ópezTorres , J.R. Luhder , M. Lunardon , G. Luparello , Y.G. Ma , A. Maevskaya , M. Mager ,S.M. Mahmood , T. Mahmoud , A. Maire , R.D. Majka
146 ,i , M. Malaev , Q.W. Malik , L. Malinina
75 ,iv ,D. Mal’Kevich , P. Malzacher , G. Mandaglio
32 ,56 , 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
53 ,138 , A.M. Mathis , 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 ,93 ,A. Menchaca-Rocha , C. Mengke , E. Meninno
29 ,114 , M. Meres , S. Mhlanga , Y. Miake ,L. Micheletti , L.C. Migliorin , D.L. Mihaylov , K. Mikhaylov
75 ,92 , 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 ,54 , 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 ,54 , P. Nomokonov , J. Norman
79 ,127 , N. Novitzky ,P. Nowakowski , A. Nyanin , J. Nystrand , M. Ogino , A. Ohlson
81 ,104 , J. Oleniacz , A.C. Oliveira DaSilva , 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
50 ,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
34 ,54 , L. Pinsky , C. Pinto ,S. Pisano
10 ,52 , D. Pistone , M. Płosko´n , M. Planinic , F. Pliquett , 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 , S. Qiu ,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
96 ,130 , A.R. Redelbach , K. Redlich
85 ,vi , A. Rehman , P. Reichelt ,F. Reidt , X. Ren , R. Renfordt , Z. Rescakova , K. Reygers , V. Riabov , T. Richert
81 ,89 , 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 , A. Rosano , / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration
E.D. Rosas , K. Roslon , A. Rossi
28 ,57 , 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
93 ,98 , 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 ,96 , A.R. Schmier , J. Schukraft , Y. Schutz
34 ,136 , K. Schwarz ,K. Schweda , G. Scioli , E. Scomparin , J.E. Seger , Y. Sekiguchi , D. Sekihata ,I. Selyuzhenkov
93 ,107 , S. Senyukov , D. Serebryakov , 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.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 , A. Trifiró
32 ,56 ,S. Tripathy
50 ,69 , T. Tripathy , S. Trogolo , G. Trombetta , L. Tropp , V. Trubnikov , W.H. Trzaska ,T.P. Trzcinski , B.A. Trzeciak
37 ,63 , 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. vanDoremalen , M. van Leeuwen , P. Vande Vyvre , D. Varga , Z. Varga , M. Varga-Kofarago ,A. Vargas , M. Vasileiou , A. Vasiliev , O. Vázquez Doce , V. Vechernin , E. Vercellin , S. VergaraLimó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 ,54 , G.A. Willems , E. Willsher , B. Windelband , M. Winn , W.E. Witt ,J.R. Wright , 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 / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration 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 di Scienze MIFT, Università di Messina, Messina, 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 / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration INFN, Sezione di Torino, Turin, Italy INFN, Sezione di Trieste, Trieste, Italy Inha University, Incheon, Republic of Korea 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 des 2 Infinis, Irôlne Joliot-Curie, Orsay, France 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 / ψ production in p–Pb at √ s NN = .
16 TeV ALICE Collaboration
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
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