Observation of prompt J/$\psi$ meson elliptic flow in high-multiplicity pPb collisions at $\sqrt{s_\mathrm{NN}} =$ 8.16 TeV

A measurement of the elliptic flow ($v_2$) of prompt J/$\psi$ mesons in high-multiplicity pPb collisions is reported using data collected by the CMS experiment at a nucleon-nucleon center-of-mass energy $\sqrt{s_\mathrm{NN}} =$ 8.16 TeV. Prompt J/$\psi$ mesons decaying into two muons are reconstructed in the rapidity region in the nucleon-nucleon center-of-mass frame ($y_\mathrm{cm}$), corresponding to either $-$2.86 $\lt y_{\mathrm{cm}}<$ $-$1.86 or 0.94 $\lt y_{\mathrm{cm}}<$ 1.94. The average $v_2$ result from the two rapidity ranges is reported over the transverse momentum (p$_\mathrm{T}$) range from 0.2 to 10 GeV. Positive $v_2$ values are observed for the prompt J/$\psi$ meson, as extracted from long-range two-particle correlations with charged hadrons, for 2 $\lt p_{\mathrm{T}}<$ 8 GeV. The prompt J/$\psi$ results are compared with previous CMS measurements of elliptic flow for open charm mesons (D$^0$) and strange hadrons. From these measurements, constraints can be obtained on the collective dynamics of charm quarks produced in high-multiplicity events arising from small systems.


Introduction
Strong collective behavior is found in the azimuthal correlations of particles emitted in relativistic nucleus-nucleus (AA) collisions at the BNL RHIC [1-4] and at the CERN LHC [5-10]. These correlations, which are long-range in pseudorapidity (η), suggest the formation of a strongly interacting quark-gluon plasma (QGP) that exhibits nearly ideal hydrodynamic behavior [11][12][13]. The azimuthal correlation structure of emitted particles is typically characterized by its Fourier components [14]. In particular, within a hydrodynamic picture, the second and third Fourier anisotropy components are known as elliptic (v 2 ) and triangular (v 3 ) flow, respectively, and reflect the QGP medium response to the initial collision geometry and its fluctuations [15][16][17]. In recent years, similar long-range collective azimuthal correlations have also been observed in events with high final-state particle multiplicity in proton-proton (pp) [18][19][20][21], proton-nucleus (pA) [22][23][24][25][26][27][28][29][30], and lighter AA collisions [31][32][33], raising the question of whether a fluid-like QGP is created in these much smaller systems. While experimental measurements in these small systems are consistent with the hydrodynamic expansion of a tiny QGP droplet, alternative scenarios based on gluon saturation in the initial state also claim to capture the main features of the correlation data (recent reviews are provided in Refs. [34,35]).
Heavy-flavor quarks (charm and bottom) are primarily produced at a very early stage via initial hard scattering because of their large masses. As such, they are largely decoupled from the bulk production of soft gluons and light-flavor quarks at a later stage in AA collisions, and thereby probe the properties and dynamics of the QGP through its entire evolution [36]. A strong elliptic flow (v 2 ) signal has been observed for open heavy-flavor D 0 mesons in both AuAu collisions at RHIC [37] and PbPb collisions at the LHC [38][39][40], suggesting that charm quarks may develop strong collective flow behavior. Furthermore, a recent measurement of the elliptic flow of J/ψ mesons in PbPb collisions at √ s NN = 5.02 TeV [41] has provided additional evidence for the collective behavior of charm quarks in the QGP.
In the study of collectivity in small systems, such as that occurring in pp or pPb collisions, a key open question is whether the strong collective behavior observed for bulk constituents in high-multiplicity events also extends to charm and bottom quarks. Long-range correlations involving inclusive muons at high transverse momentum (p T ) reveal a hint of heavy-flavor quark collectivity in pPb collisions [42]. Furthermore, the recent observation of a significant elliptic flow signal for prompt D 0 mesons in pPb collisions has provided evidence for charm quark collectivity in a small system [43]. The v 2 signal for D 0 mesons is found to be smaller than that of light-flavor hadrons at a given p T , indicating that in these small systems there is a weaker collective motion for charm quarks, as compared to that of the bulk medium, than found in large AA systems. However, as the D 0 meson carries both a light and a charm quark, the relative contribution of these different flavor quarks to the observed v 2 signal is not fully constrained. Without detailed theoretical modeling, a scenario is not excluded where the D 0 meson v 2 signal is entirely carried by the light-flavor quark. The observation of an elliptic flow signal for J/ψ mesons in a small system could provide more direct evidence of charm quark collectivity and could impose new constraints on the collective dynamics of heavy-quark production in such collisions. Furthermore, heavy-quark collectivity may also provide a hint of how, in small systems, hard probes interact with the QGP [36], assuming this is formed. A measurement of inclusive J/ψ (combined charmonia and J/ψ mesons from decay of open beauty hadrons) v 2 in pPb collisions was reported in Ref. [44], where positive v 2 coefficients were found in the range of 3 < p T < 6 GeV with center-of-mass rapidities −4.46 < y cm < −2.96 or 2.03 < y cm < 3.53. A recent model calculation of J/ψ v 2 in pPb collisions suggests little v 2 signal arising from final-state interactions between charm quarks and the QGP medium [45]. This Letter presents the first measurement of prompt J/ψ meson elliptic flow (excluding contributions from B hadron decays) from long-range two-particle correlations in very high multiplicity pPb collisions at √ s NN = 8.16 TeV. The v 2 harmonics for prompt J/ψ mesons in the ranges −2.86 < y cm < −1.86 and 0.94 < y cm < 1.94 are determined over a wide p T range from 0.2 to 10 GeV. To estimate the possible residual contribution from back-to-back jet-like correlations, the v 2 values are also presented after subtracting correlations obtained from lowmultiplicity pPb events (denoted as v sub 2 ), where jet-like correlations are assumed to dominate. The results are compared to those of the light strange-flavor K 0 S and Λ hadrons, and the open heavy-flavor prompt D 0 meson, which were previously reported by CMS [43] in the same p T range but in a different rapidity range of −1.46 < y cm < 0.54. In order to explore possible collectivity at the partonic level, a comparison is also presented in terms of the transverse kinetic energy per constituent quark (KE T /n q , where KE T = √ m 2 + p 2 T − m, and n q is the number of constituent quarks).

The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume, there are four primary subdetectors including a silicon pixel and strip tracker detector, a lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron calorimeter, each composed of a barrel and two endcap sections. Iron and quartz-fiber Cherenkov hadron forward (HF) calorimeters cover the range 3.0 < |η| < 5.2. Muons are measured in the range |η| < 2.4 in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid, with detection planes made using three technologies: drift tubes, cathode strip chambers, and resistive-plate chambers. The silicon tracker measures charged particles within the range |η| < 2.5. For charged particles with 1 < p T < 10 GeV and |η| < 1.4, the track resolutions are typically 1.5% in p T and 25-90 (45-150) µm in the transverse (longitudinal) impact parameter [46]. A detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [47].

Data selection and J/ψ meson reconstruction
The pPb data at √ s NN = 8.16 TeV used in this analysis were collected in 2016, and correspond to an integrated luminosity of 186 nb −1 . The beam energies are 6.5 TeV for the protons and 2.56 TeV per nucleon for the lead nuclei. Because of the asymmetric beam conditions, particles selected in the laboratory rapidity range of 1.4 < y lab < 2.4 (−2.4 < y lab < −1.4) have a corresponding nucleon-nucleon center-of-mass frame rapidity range of 0.94 < y cm < 1.94 (−2.86 < y cm < −1.86), with positive rapidity defined in the proton beam direction. To minimize statistical uncertainties, the quoted J/ψ meson v 2 results combine the individual values obtained for the proton and lead beam directions.
The pPb data are analyzed in different ranges of N offline trk , where N offline trk is the number of primary charged particle tracks [46] with |η| < 2.4 and p T > 0.4 GeV. The main results are obtained with events in the high-multiplicity range 185 ≤ N offline trk < 250. To select these events, dedicated triggers were developed, as discussed in Refs. [48,49]. Events with N offline trk < 35 are also used to estimate the possible contribution of residual back-to-back jet-like correlations. These lower-multiplicity events are selected online with a hardware-based trigger requiring two muon candidates in the muon detectors with no explicit momentum or rapidity threshold [50]. In the offline analysis, hadronic collisions are selected by requiring at least one HF calorimeter tower with more than 3 GeV of total energy in each of the two HF detectors. Events must contain a primary vertex close to the nominal interaction point of the beams, within 15 cm along the beam direction, and 0.2 cm in the plane transverse to beam direction. The N offline trk range limits correspond to fractional inelastic cross sections from 100 to 57% for N offline trk < 35, and from 0.33 to 0.01% for 185 ≤ N offline trk < 250, respectively.
The offline muon reconstruction algorithm starts either by finding tracks in the muon detectors, which are then fitted together with tracks reconstructed in the silicon tracker (global muons), or by extrapolating tracks from the silicon tracker to match a hit on at least one segment of the muon detectors (tracker muons). The muon candidates are required to pass the identification criteria of the particle-flow algorithm [51], which suppresses contamination of "punchthrough" hadrons misidentified as muons, based on energy deposition in the calorimeters. The soft muon selection criteria are also imposed, as defined in Ref. [52], to further improve the purity of muons.
The J/ψ meson candidates are formed from pairs of oppositely charged muons, originating from a common vertex. Based on the vertex probability distributions for signal and background candidates, the probability that the dimuon pair shares a common vertex is required to be larger than 1%, lowering the background from random combinations as well as from semileptonic decays of bottom and charm hadrons. Because of the long lifetime of b hadrons compared to that of J/ψ mesons, the nonprompt J/ψ meson component can be reduced by placing constraints on the pseudo-proper decay length [53]. This is defined by 3D J/ψ = L xyz m J/ψ /|p µµ |, where L xyz is the distance between the primary and dimuon vertices, m J/ψ is the Particle Data Group [54] world average value of the J/ψ meson mass (assumed for all dimuon candidates), and p µµ is the dimuon momentum. The upper limit (decreasing as a function of p T ) imposed on the 3D J/ψ value is based on Monte Carlo (MC) studies with simulated event samples of PYTHIA 8.209 [55,56], and found to reject 75-90% (from low to high p T ) of nonprompt J/ψ mesons, largely independent of multiplicity. The residual nonprompt J/ψ meson fraction in the data is estimated to be approximately 5% across the full p T range, and its effect on the v 2 measurement is propagated as a systematic uncertainty, as described in Section 5.

Analysis technique
The azimuthal anisotropy of J/ψ mesons is extracted from the long-range (|∆η| > 1) twoparticle azimuthal correlations, following an identical procedure to that described in Refs. [21,27,43]. A two-dimensional (2D) correlation function is constructed by pairing each J/ψ candidate with reference primary charged-particle tracks with 0.3 < p T < 3 GeV and |η| < 2.4 (denoted "ref" particles), and calculating where ∆η and ∆φ are the differences in η and in the azimuthal angle (φ) of the pair. The same-event pair distribution, S(∆η, ∆φ), represents the yield of particle pairs normalized by the number of J/ψ candidates from the same event. The mixed-event pair yield distribution, B(∆η, ∆φ), is constructed by pairing J/ψ candidates in each event with the reference primary charged-particle tracks from 20 different randomly selected events, from the same N offline trk range and having a primary vertex falling in the same 2 cm wide range of reconstructed z coordinate. The analysis procedure is performed in each p T and invariant mass (m inv ) range of J/ψ candi-dates. A correction for the acceptance and efficiency of the J/ψ meson yields is applied, but found to have a negligible effect on the measurements. The ∆φ correlation functions averaged over |∆η| > 1 (to remove short-range correlations, such as jet fragmentation) are then obtained from the 2D distributions and fitted by the first three terms of a Fourier series (including additional terms has a negligible effect on the fit results): Here, V n∆ are the Fourier coefficients and N assoc represents the total number of same-event pairs per J/ψ candidate for a given invariant mass interval. By assuming that V n∆ is the product of single-particle anisotropies of J/ψ mesons and reference charged particles [57], , the v n anisotropy harmonics for J/ψ candidates can be extracted as a function of invariant mass, v n (J/ψ)  To extract the genuine v 2 values of the J/ψ meson signal (v S 2 ), the contribution from background candidates (v B 2 ) has to be subtracted from the v 2 values of all J/ψ meson candidates, as obtained in the previous step. The procedure is to first fit the dimuon mass spectrum with a function composed of three components: two Crystal Ball functions [58] with different widths but common mean and tail parameters for the J/ψ signal (the tail parameters are fixed to the values obtained from simulation), S(m inv ), and an exponential function to model the combinatorial background, B(m inv ). Then, the signal plus background v S+B where Here, v B 2 (m inv ) for the background J/ψ candidates is modeled as an exponential function of the invariant mass, and α(m inv ) is the J/ψ signal fraction obtained from the mass spectrum fit. An example of fits to the mass spectrum and v S+B 2 (m inv ) in the p T interval 6.0-8.0 GeV for the multiplicity range 185 ≤ N offline trk < 250 is shown in Fig. 1. The residual contribution of back-to-back dijets to the measured v 2 results is estimated from low-multiplicity pPb events and is removed from the signal after accounting for the jet yield ratio of the selected events, following a jet subtraction procedure similar to that established in Refs. [21,43,57]. The Fourier coefficients, V n∆ , extracted from Eq. (2) for N offline trk < 35, are subtracted from the V n∆ coefficients obtained in the high-multiplicity region, with Here, Y jet represents the jet yield obtained by integrating the difference of the short-range (|∆η| < 1) and long-range event-normalized associated yields for each multiplicity class. The ratio, Y jet /Y jet (N offline trk < 35), is introduced to account for the enhanced jet correlations resulting from the selection of higher-multiplicity events. For p T (J/ψ) < 4.5 GeV, the jet yield ratio cannot be directly estimated from the two-particle azimuthal correlations, as the J/ψ candidates tend to have larger η values than the acceptance for charged particles. Therefore, the value is assumed to be the same as that for the high-p T region, where no p T dependence has been observed. It was also previously observed that the values of jet yield ratio for D 0 and strange particle species show little dependence on p T over the full p T range [43].

Systematic uncertainties
Sources of systematic uncertainties on the prompt J/ψ meson v 2 measurement include the J/ψ meson yield correction (acceptance and efficiency correction derived from PYTHIA simulation), the nonprompt J/ψ meson contamination, the background v B 2 (m inv ) function form, the signal and background invariant mass PDF, the jet subtraction procedure, the contamination of events containing more than one pPb interaction (pileup), and the trigger bias. In this Letter, the quoted uncertainties in v 2 are absolute values, and are found to have no dependence on p T , except those for the jet subtraction procedure. Systematic uncertainties originating from different sources are added in quadrature to obtain the overall systematic uncertainty shown as boxes in the figures.
To evaluate the uncertainties arising from the efficiency correction to the J/ψ meson yield, the v 2 values are compared to the uncorrected ones, yielding an uncertainty of 0.008. The effect on the measured v 2 due to the residual contribution from nonprompt J/ψ mesons is evaluated by varying the 3D J/ψ requirement such that the nonprompt J/ψ meson yield is doubled. The v 2 values are found not to change by more than ±0.004, which is assigned as the systematic uncertainty due to the J/ψ meson yield correction. Possible differences in the rejection efficiency of nonprompt J/ψ mesons between data and simulation are investigated and found to be negligible. The systematic uncertainties from the background v 2 function form are evaluated by comparing v B 2 (m inv ) values based on first-, second-, and third-order polynomial fits to the background distribution. The resulting J/ψ signal v 2 values are found to vary by less than 0.009. Systematic effects related to signal invariant mass PDF are found to be negligible by releasing, one at a time, the fixed tail parameters of the Crystal Ball functions. The variation of v 2 , while changing the background invariant mass PDF to a second-or third-order polynomial function is also found to be negligible. In the jet subtraction procedure, the statistical precision of the jet yield ratio is limited. The v sub 2 results are found to be consistent within ±0.002 to ±0.014 (increasing with p T ) when varying the jet yield ratio by its statistical uncertainty. The systematic uncertainties from the potential pileup effect and the trigger bias are taken to be the same as for inclusive charged particles in Ref.
[49], where they can be established with good statistical precision. The pileup and trigger bias uncertainties are negligible compared to the other sources of systematic uncertainties, as the fraction of residual pileup events is only a few % and the trigger efficiency is close to 100%. Figure 2 shows the v 2 results of prompt J/ψ mesons at forward rapidities (−2.86 < y cm < −1.86 or 0.94 < y cm < 1.94) for high-multiplicity (185 ≤ N offline trk < 250) pPb collisions, covering a p T range from 0.2 to 10 GeV. Results obtained separately for J/ψ meson rapidity in the Pb-and pgoing direction are compared, and found to be consistent within statistical uncertainties. Thus, as mentioned earlier, combined v 2 values are presented for the best statistical precision. The  Positive prompt J/ψ meson v 2 values are observed over a wide p T range from about 2 to 8 GeV. The prompt J/ψ meson v 2 results show a trend of first increasing up to p T ∼ 4 GeV and then decreasing toward higher p T . This observed trend appears to be in common with the other hadron species shown. In the p T range below 5 GeV, the v 2 values for J/ψ and D 0 mesons are consistent with each other within the uncertainties, while an indication of smaller v 2 values for J/ψ mesons than that for D 0 mesons is seen for p T > 5 GeV. Over the full p T range, the v 2 signal values for both J/ψ and D 0 hadrons are smaller than those for K 0 S and Λ hadrons. This observation is consistent with the earlier conclusion that charm quarks develop a weaker collective dynamics than light quarks in small systems [43], unlike what is seen in AA collisions. Because of experimental limitation, v 2 values for the prompt J/ψ meson and the other meson species are not compared within the same rapidity range, possibly affecting their comparison. The rapidity dependence of v 2 values for charged particles in pPb collisions has been measured [59, 60], suggesting up to around 15% variation from |y lab | ∼ 0 to 2.4.

Results
To better study the elliptic flow signal coming purely from long-range collective correlations, the J/ψ v 2 results are corrected for residual jet correlations. The resulting (v sub 2 ) values are shown in Fig. 3 (upper) for prompt J/ψ mesons as a function of p T with 185 ≤ N offline trk < 250, and compared to similarly corrected K 0 S , Λ, and D 0 hadron results [43]. The effect of the correction for all particle species is most noticeable at very high p T , while the overall p T dependence of the v 2 data remains unchanged. The K 0 S mesons have a larger correction applied to their v 2 values (possibly because K 0 S mesons are more correlated with the bulk multiplicity, and thus are biased toward stronger jet correlations due to the selection of high multiplicities) and their v sub 2 values after the correction tend to converge to those of the prompt J/ψ and D 0 mesons at high p T .
A recent model calculation of J/ψ v 2 in pPb collisions, based on final-state interactions between produced charm quarks and a QGP medium, suggests a very small v 2 signal of less than 0.01 [45]. This calculation indicates that additional contributions, e.g., those from initialstate interactions, may be needed to account for the observed v 2 signal of prompt J/ψ mesons reported in this Letter. values is observed between the K 0 S meson and Λ baryon, shown in Fig. 3 (lower). This scaling between light baryon and meson species systems produced in the collision (known as the number-ofconstituent-quark or NCQ scaling) was first discovered in AA colliding systems [61, 62, 66], indicating that collectivity is first developed among the partons, which later recombine into final-state hadrons. The elliptic flow signal per quark (v sub 2 /n q ) for prompt J/ψ mesons at low KE T /n q range is consistent with those of K 0 S , Λ, and prompt D 0 hadrons within large statistical uncertainties for the current data. There is a hint that the prompt J/ψ meson data tend to fall on the same trend as those of K 0 S and Λ baryons, all of which are above the prompt D 0 meson data. A more definitive conclusion could be drawn with future high precision data. For KE T /n q > 1 GeV, the v sub 2 /n q for prompt D 0 and J/ψ mesons are consistently below that of the K 0 S meson. An indication of smaller v sub 2 /n q values for J/ψ mesons than for D 0 mesons is seen for KE T /n q > 1 GeV. As J/ψ mesons contain two charm quarks, while D 0 mesons contain a charm and a light-flavor quark, this observation would be consistent with a weaker collective behavior of heavy-flavor quarks than light quarks, possibly a consequence of the much smaller size of the collision system. Future data with improved precision will provide crucial insights to fully constrain the collective behavior of light-and heavy-flavor quarks in high-multiplicity, small systems.

Summary
In summary, the elliptic flow harmonic (v 2 ) for prompt J/ψ mesons in high-multiplicity protonlead (pPb) collisions at √ s NN = 8.16 TeV is presented as a function of transverse momentum (p T ). Positive v 2 values are observed for prompt J/ψ mesons at forward rapidity (−2.86 < y cm < −1.86 or 0.94 < y cm < 1.94) over a wide p T range (2 < p T < 8 GeV). This observation provides evidence for charm quark collectivity in high-multiplicity pPb collisions, similar to that first observed for light-flavor hadrons. The observed ordering of v 2 among light-flavor, open and hidden heavy-flavor hadrons at intermediate and high-p T regions (e.g., above 4 GeV) adds support to the earlier conclusion that heavy quarks exhibit weaker collective behavior than light quarks or gluons in small systems, unlike what is found in AA collisions. For particle transverse kinetic energy per constituent quark values less than 1 GeV, the v 2 of prompt J/ψ mesons is consistent with prompt D 0 , K 0 S and Λ hadrons, within current uncertainties. A model calculation based on final-state interactions between charm quarks and a QGP medium in pPb collisions significantly underestimates the measured prompt J/ψ v 2 signal. The new prompt J/ψ meson results, together with previous results for light-flavor and open heavy-flavor hadrons, provide novel insights into the dynamics of the heavy quarks produced in small systems that lead to high final-state multiplicities.       [27] CMS Collaboration, "Long-range two-particle correlations of strange hadrons with charged particles in pPb and PbPb collisions at LHC energies", Phys. [35] J. L. Nagle and W. A. Zajc, "Small system collectivity in relativistic hadron and nuclear collisions", (2018). arXiv:1801.03477.