Transverse momentum and process dependent azimuthal anisotropies in √ s NN = 8 . 16 TeV p + Pb collisions with the ATLAS detector The ATLAS Collaboration

The azimuthal anisotropy of charged particles produced in $$\sqrt{s_{\mathrm {NN}}}=8.16$$sNN=8.16 TeV p+Pb collisions is measured with the ATLAS detector at the LHC. The data correspond to an integrated luminosity of 165 $$\mathrm {nb}^{-1}$$nb-1 that was collected in 2016. Azimuthal anisotropy coefficients, elliptic $$v_2$$v2 and triangular $$v_3$$v3, extracted using two-particle correlations with a non-flow template fit procedure, are presented as a function of particle transverse momentum ($$p_\mathrm {T}$$pT) between 0.5 and 50 GeV. The $$v_2$$v2 results are also reported as a function of centrality in three different particle $$p_\mathrm {T}$$pT intervals. The results are reported from minimum-bias events and jet-triggered events, where two jet $$p_\mathrm {T}$$pT thresholds are used. The anisotropies for particles with $$p_\mathrm {T}$$pT less than about 2 GeV are consistent with hydrodynamic flow expectations, while the significant non-zero anisotropies for $$p_\mathrm {T}$$pT in the range 9–50 GeV are not explained within current theoretical frameworks. In the $$p_\mathrm {T}$$pT range 2–9 GeV, the anisotropies are larger in minimum-bias than in jet-triggered events. Possible origins of these effects, such as the changing admixture of particles from hard scattering and the underlying event, are discussed.


Introduction
The collisions of heavy nuclei at relativistic speeds generate hot and dense droplets of matter composed of deconfined quarks and gluons known as the quark-gluon plasma (QGP) [1,2]. Studies of the QGP at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) have yielded a wealth of surprising results that reveal a complex set of QGP-related phenomena. Bulk hadron production, occurring mainly at low transverse momentum (p T 3 GeV), exhibits significant azimuthal anisotropies that are well described in terms of nearly inviscid hydrodynamic flow of the QGP [1]. The final hadron momentum anisotropies arise from inhomogeneities in the initial spatial distribution of the QGP translated to momentum space via strong differential pressure gradients. These anisotropies are characterised in terms of a Fourier decomposition: where v n are the anisotropy coefficients, Ψ n is the n th -order orientation of the anisotropy, and the normalization, G, is set by the integral of the distribution. In particular, v 2 and v 3 are referred to as the elliptic and triangular coefficients, respectively.
In addition, the production of high transverse momentum hadrons (p T 10 GeV) is highly suppressed relative to the yields expected from nuclear thickness scaling of proton-proton collision yields [3]. This suppression is understood to result from high momentum transfer parton-parton interactions followed by the outgoing partons losing energy via radiative and collisional processes in the QGP -processes referred to as jet quenching [4][5][6]. These high-p T hadrons and associated jets are also observed to have a non-zero azimuthal anisotropy [7][8][9], despite being well outside the nominal domain where the anisotropies are interpreted in terms of hydrodynamic flow. Instead, these anisotropies are understood to also arise from inhomogeneities in the initial spatial distribution of the QGP, but in this case, where the jet quenching effect is stronger for partons traversing longer paths through the QGP and weaker for partons traversing shorter paths [10]. In this way, low-and high-p T hadrons have a common orientation of their azimuthal anisotropy in a given event, because both are correlated with the orientation of the initial geometry of the colliding nucleons. It is notable that, for more than a decade, an outstanding and challenging theoretical puzzle has been how to quantitatively describe both high-p T hadron suppression and azimuthal anisotropy in Pb+Pb collisions [11]. There are a number of proposed explanations for resolving this puzzle in heavy-ion collisions -see Refs. [12][13][14][15][16][17] for examples.
Measurements in smaller collision systems, pp and p+Pb collisions at the LHC [18][19][20][21][22][23] and p+Au, d+Au, and 3 He+Au at RHIC [24], indicate significant azimuthal anisotropies for low-p T hadrons with patterns quite similar to those observed in the larger heavy-ion collision systems. For a recent review see Ref. [25]. These observations have raised the question of whether smaller and shorter-lived droplets of QGP are formed in these smaller collision systems. Indeed, models employing nearly inviscid hydrodynamics for the QGP provide a quantitative description of this large body of data at low p T [26].
In contrast, measurements aimed at observing signatures of jet quenching in small collision systems have found no such effect. Measurements of hadron and jet p T spectra at high p T indicate production yields consistent with those in pp collisions scaled up by the expected nuclear thickness in p+Pb [27-29] and d+Au collisions [30], and that the p T -balance between dijets or hadron-jet pairs is unmodified in p+Pb collisions within uncertainties [31,32]. The ATLAS experiment has also published results for the hadron azimuthal anisotropy up to p T ≈ 12 GeV that hint at a non-zero anisotropy extending into the region beyond The correlation functions, S(∆φ), are defined as the yields of particle pairs passing the above event and pair selection, binned in ∆φ = φ A − φ B , and normalised by the total number of A-particles. Corrections for the imperfect trigger and tracking efficiencies are applied as weights to the entries of the correlation functions. A mixed-event correction, M(∆φ), is generated by correlating A-particles from one event with B-particles from a different event with a vertex z-position differing by less than 10 mm and a number of reconstructed charged particles (N rec ch ) differing by less than 10 for N rec ch < 100 and less than 20 for N rec ch > 100. Thus, the mixed events contain only trivial detector acceptance effects and no physical correlations. To reduce the statistical uncertainty introduced by the correction, each event is mixed with five others meeting the above vertex z and N rec ch conditions. The corrected correlation is, then, Y (∆φ) = S(∆φ)/M(∆φ), where M(∆φ) is normalised such that the ratio preserves the overall integral of S(∆φ). Jet events are mixed with other jet events, and the |∆η jB | condition is applied with respect to the jets in the A-particle event only. Thus, the B-particle acceptance is consistent between the same-and and mixed-event correlations.
To extract the anisotropy coefficients while accounting for residual non-flow, the ATLAS template fitting procedure, as used for previous results [19,20], is applied to Y (∆φ). In this procedure, Y (∆φ) is found for two different selections of event activity quantified by centrality: a central selection, Y cent (∆φ), and a peripheral selection, Y peri (∆φ). In this analysis, the peripheral selection corresponds to the 60-90% centrality interval. Assuming that the shape of the non-flow correlation is independent of centrality, where F and each v n,n are parameters of a global χ 2 fit, and G is fixed by the requirement that the integral of the fit model is that of Y cent . The parameter F allows for a linear scaling of the non-flow between the two centrality classes. The fit includes the fourth harmonic, v 4,4 , but it is not presented in the results because it is statistically insignificant. The fit χ 2 function incorporates the statistical uncertainties from both Y cent and the peripheral template, Y peri , although the examples shown in Figures 1 and 2 do not show the uncertainties of Y peri for readability. The statistical uncertainties of the extracted v n,n parameters are returned from the MINUIT χ 2 minimiser [55], accounting for correlations between parameters. Figure 1 shows an example of two template fits using jet-triggered events with jet p T > 100 GeV. The left plot shows the fit for correlations made without the B-particle jet rejection condition, and the right plot shows the same correlation, but with the condition |∆η jB | > 1 applied. In this figure, Y ridge N represents the N th -order harmonic component of the fit. The left plot has a dominant non-flow contribution, and a distortion in the resulting subtracted distribution is observed near ∆φ ≈ π. Removing much of the jet correlation in this way reduces the overall sensitivity to the template method assumption that the shape of the non-flow contribution is the same for the central and peripheral selections. However, violation of this assumption will introduce distortions that could potentially bias the harmonic coefficients. This is explicitly tested by varying the centrality selection of the peripheral template, as discussed further in Section 6. Two additional examples of template fits from the jet-triggered events with jet p T > 100 GeV and with the B-particle jet rejection are plotted in Figure 2. These show the behaviour of the template fits for high A-particle p T .
If the particle momentum correlations originate from a global field, as is the case for collective expansion, the v n,n will factorise such that v n,n (p A T , p B T ) = v n (p A T ) · v n (p B T ). By assuming this relation and making specific p T selections on A-and B-particles, the single-particle v n (p A T ) can be obtained from   Figure 1: Template fitting output for events with jet p T > 100 GeV. Both require 3.5 < p A T < 4.0 GeV and are made with 60-90% peripheral selection and 0-5% central selection. The left plot is made with no selection on the B-particles and the right plot is made requiring the B-particles to have |∆η jB | > 1 relative to all jets with p jet T > 15 GeV in the event. In the upper panels, the open circles show the scaled and shifted peripheral template with uncertainties omitted, the closed circles show the central data, and the red line shows the fit (template and harmonic functions). The blue dashed line shows the second-order harmonic component, Y (the n = 2 and n = 3 contributions to the sum in Eq. 1, respectively). The lower panels show the difference between the central data and the peripheral template along with the second and third harmonic functions. The resulting v 2,2 , v 3,3 , and global fit χ 2 /NDF values are reported in the legends. In these fits, NDF= 35.
where v n,n (p A T , p B T ) is determined with A-and B-particles having p T in range p A T and p B T , respectively, and v n,n (p B T , p B T ) is determined with A-and B-particles both having p T in range p B T . In this analysis, this range is nominally p B T > 0.4 GeV, although the dependence of the extracted anisotropy on this choice is explored in Section 7. The relative yield of particle pairs entering the correlation functions is estimated assuming a simple, two-component model of particle production. Particles are assumed to be produced either by hard scattering (HS) processes, such as jet production, or by soft underlying event (UE) processes. With this assumption, the correlation functions are constructed from pairs pulled from a mixture of the two sources. Particle pairs can be formed in the following four A-B combinations: UE-UE, UE-HS, HS-UE, and HS-HS. The event-by-event yields of the UE and HS processes are estimated by classifying the charged particles according to their azimuthal orientation relative to the leading jet or, in the case of MBT events that contain no jets with p T > 15 GeV and |η| < 4.9, relative to the leading hadron. The following regions are defined relative to this leading vector:   Then, assuming that HS particles are completely contained in the towards region and the UE particles are distributed uniformly in azimuth, the following relations are inferred: where N UE and N HS are the single-particle yields from UE and HS processes, respectively, and N trans and N toward are the particle yields in the transverse and toward regions, respectively. The quantities N UE and N HS are statistically determined from the event averaged N trans and N toward and are, thus, insensitive to event-by-event fluctuations. However, it is not possible to classify individual particles. It should be noted that the assumptions used in this derivation are likely not perfect; for example, the UE is not uniformly distributed in φ, event by event, due to the presence of azimuthal anisotropy. The leading object may be more likely to be oriented with the anisotropy, in which case the UE yield would be underestimated and the HS yield overestimated. However, the analysis proceeds with the assumptions as given and includes no additional uncertainty for this potential effect.
The ∆η AB and ∆η jB rejections produce a geometric correlation between the yields of A-and B-particles and, thus, the number of pairs is not the simple product of the two individual yields. Accounting for the dependencies, the total yield of particle pairs can be expressed in the following way where X-Z could be any pairwise combination of UE and HS. In the case of jet events, the ∆η jB condition is enforced when filling the η A and ∆η AB distributions so it's effects are taken into account.

Systematic uncertainties
The systematic uncertainties fall into two categories: those associated with both the MBT and jet events and those associated with only the jet events. The uncertainties are determined by assessing the difference between the nominal value of v 2 or v 3 and the value after a given variation. Unless otherwise stated, the uncertainties are defined as asymmetric one-standard-deviation errors. The final uncertainty is the quadrature sum of the uncertainty from each individual source. The relative downward and upward systematic uncertainties from different sources and all sources combined are shown for v 2 in Table 1 and for v 3 in Table 2. The rest of this section focuses primarily on the systematic uncertainties of v 2 . While the absolute uncertainties in v 2 and v 3 are of similar magnitude, this represents larger relative uncertainties in the v 3 values since they are generally smaller than the v 2 values at any given p T .
For both the MBT and jet events, the sensitivity to the trigger and tracking efficiency corrections was assessed by removing each. This variation ('Track/trig Eff.' in Table 1) yields a 0-2% (2-4%) difference for MBT (jet) events, depending on track p T , and is subdominant. In the construction of the correlation functions, the uncertainty in the mixed-event correction was again found by removing it from the analysis. This variation results in an uncertainty that vanishes at low p T but grows to 20% (10%) at high p T for MBT (jet) events, but remains subdominant to statistical uncertainties over the whole p T range. Regarding the template fitting procedure, the centrality range for the peripheral reference selection was varied from the nominal 60-90% to 50-70% and 70-90%. This variation ('Peri. reference') results in an uncertainty of about 2% at low p T and increasing to about 10% or 18% in the mid p T range between 2 and 10 GeV depending on the event trigger. This last uncertainty is dominant in this category for most of the p T range probed in the measurement. At high p T , the sensitivity of the measurements in MBT events to the mixed event correction and reference selection is significantly higher than in jet events; this is particularly noticeable for the v 3 values, where the relative uncertainties in the MBT events for p T > 10 GeV are 5-10 times larger than in the jet-triggered events.
The following set of uncertainties is associated with jet events only. To assess the sensitivity to the uncertainty in the jet energy scale and the impact of imperfect trigger efficiency, the jet p T thresholds used to select events were varied from 75 GeV and 100 GeV to 80 GeV and 105 GeV, respectively. This variation ('Trig jet p T ') results in a symmetric uncertainty that is smaller than 1% at low particle p T and that increases to about 10% with increasing p T . It is subdominant to other sources in this category. The jets used in the B-particle jet rejection were varied to include only jets with p T greater than 20 GeV instead of the nominal 15 GeV ('Reject jet p T '). The 2% and 10% differences at low and high p T are incorporated as a symmetric uncertainty that is subdominant to other sources in this category. The ∆η jB rejection allows jets to be composed of only a single particle that may originate in the tail of the UE particle p T spectrum. Thus, the jets used in this rejection were varied to require at least three tracks in a ∆R = 0.4 cone around the jet axis ('Reject jet mult.'). The uncertainty associated with this variation is about 5% Source and subdominant to the others in this category. An additional uncertainty is used to cover the impact of a sector of the hadronic endcap calorimeter (HEC) being disabled for the running period. The disabled sector was in the range 1.5 < η < 3.2 and −π < φ < −π/2. This uncertainty ('Disabled HEC sector') was assessed by requiring all B-particles to be outside the pseudorapidity region of the disabled HEC. The difference was found to be less than 1% at low p T and about 20% at high p T , where it is the dominant positive uncertainty. Finally, an uncertainty is assigned to account for the potential of the UE to bias the event selection. The azimuthal modulation of the UE increases the reconstructed p T of jets aligned with the flow orientation, and, thus, the event-wise jet-p T threshold will bias the events to have more jets correlated with the flow plane. The impact of this effect on the measured results was assessed in simulation by mixing jet events with a realistic UE containing azimuthal anisotropy. The resulting uncertainty only affects v 2 , is the dominant negative uncertainty for track p T greater than 3 GeV, and is about 30% (20%) for jet-triggered events with jet p T > 75 GeV (100 GeV). The effect is larger for lower-p T jets because the UE energy contribution is independent of jet energy. For a power-law spectrum, a given threshold change has a greater fractional effect on the yield for smaller values of the threshold.
In summary, the uncertainty in v 2 from the peripheral reference selection was found to be dominant for p T less than 10 GeV for MBT events, above which, the mixed event correction uncertainty is dominant, and between 2 and 5 GeV for jet events. The uncertainties associated with the jet selection were found to be dominant for p T 10 GeV in jet events. The total uncertainty in MBT events ranges from (−1%, +2%) at low p T to about (−10%, +25%) at high p T . For jet events, the total uncertainty ranges from about (−5%, +4%) at low p T to about (−35%, +50%) and (−30%, +25%) at high p T for events with jet p T > 75 GeV and jet p T > 100 GeV respectively.
The uncertainties associated with the measurement of particle pair yields are generated from some of the variations discussed above, namely the track and trigger efficiency variation, the trigger jet p T threshold variation, and each B-particle jet rejection variation. An additional variation was made to test the assumption that the toward region contains all HS particles. The two transverse region sides were tagged as having the minimum and maximum number of tracks out of the two. The pair yields were, then, calculated using the minimum and maximum sides only, as separate variations. This variation produces the dominant uncertainty in the relative pair yields, defined in Eq. (2), for all particle combinations at all p T . Source Figure 3 shows the extracted second-(v 2 ) and third-order (v 3 ) anisotropy coefficients for the MBT events compared to those from both selections of jet events plotted as a function of A-particle p T in the range 0.5 < p T < 100 GeV. Each set of values is from events with the same 0-5% centrality selection. Points are located on the horizontal axis at the mean p T of tracks within any given bin. The v 2 and v 3 coefficients increase as a function of p T in the low p T region (p T < 2-3 GeV), then decrease (2-3 < p T < 9 GeV), and finally plateau for high p T (p T > 9 GeV). The v 2 coefficients are consistent with being independent of p T for p T > 9 GeV, while the larger uncertainties in the values of v 3 preclude any strong conclusion.  The v 2 results show agreement within uncertainties between the MBT and jet events for the low p T (p T 2 GeV) and high p T (p T 9 GeV) regions. For the intermediate p T region, the MBT events yield a higher v 2 value than jet events, although the trends are qualitatively similar. Similarly to v 2 , the v 3 results show agreement between the MBT and jet events for p T < 2 GeV, and higher values from MBT events for p T > 2 GeV. The points with different B-particle p T selections are offset slightly from the nominal horizontal-axis positions to make the uncertainties visible. For clarity, systematic uncertainties are omitted from the three sets of restricted B-particle p T ranges; they are, however, consistent with those from the inclusive results and are highly correlated between the selections.

Results
As mentioned in Section 5, if the measured anisotropy originates from a global momentum field, the v 2 and v 3 values, extracted for a given p A T range, will be independent of B-particle selection. This assumption of factorisation is explicitly tested by carrying out the analysis for different selections of p B T . Figure 4 shows the v 2 values, from each event trigger, for the nominal results using p B T > 0.4 GeV overlaid with results using 0.4 < p B T < 1 GeV, 1 < p B T < 2 GeV, and 2 < p B T < 3 GeV. The test shows factorisation breaking at the level of 5% for p A T < 5 GeV in MBT events. However, at higher p A T , the differences grow with p A T to be 10-100% from the nominal values. For jet events, factorisation holds within about 10-20% for all values of p A T , except for 4 < p A T < 9 GeV in p jet T > 100 GeV events, where it is within about 30-40%. Although the large uncertainties prevent strong conclusions from being drawn, there is a hint of a difference in behaviour at high p A T where the factorisation breaking is greater for MBT events than for jet events. This result could be due to the B-particle jet rejection scheme used for the jet events. Correlations resulting from hard-process, e.g. from back-to-back jets, specifically violate the factorisation assumption, and the B-particle jet rejection dramatically limits the contribution from these processes from entering the correlation functions in jet events. However, the correlations from MBT events have no such rejection, and could, therefore, be more susceptible to hard-process correlations at high p A T .  Figure 5: Distribution of v 2 plotted as a function of centrality for MBT events (black squares), events with jet p T > 75 GeV (blue circles), and events with jet p T > 100 GeV (orange diamonds). The results are obtained in three different selections of the A-particle p T : 0.5 < p T < 2 GeV (top left), 2 < p T < 9 GeV (top right), and 9 < p T < 100 GeV (bottom). Statistical uncertainties are shown as narrow vertical lines on each point, and systematic uncertainties are presented as coloured boxes behind the points. Figure 5 shows v 2 plotted as a function of centrality for MBT events and both classes of jet events. The results are divided into three regions in A-particle p T : 0.5 < p T < 2 GeV, 2 < p T < 9 GeV, and 9 < p T < 100 GeV. The v 2 results show agreement, within uncertainties, between the MBT and jet events for p T selections 0.5 < p T < 2 GeV and p T > 9 GeV for all centralities and are found to be nearly independent of centrality. For 2 < p T < 9 GeV, the MBT events give a higher v 2 value than the jet events, and all three sets show a trend to lower values of v 2 as the collisions become more peripheral.
Focusing on the overall p T dependence of the anisotropies, Figure 6 (left panel) shows v 2 and v 3 coefficients from events with jet p T > 100 GeV compared with theoretical calculations from Ref. [14]. This theoretical calculation, within the jet quenching paradigm, invokes a stronger parton coupling to the QGP near the  transition temperature, which helps to reduce the tension in simultaneously matching the nucleus-nucleus high-p T hadron spectrum suppression and the azimuthal anisotropy v 2 . The calculation tests two different initial p+Pb geometries referred to as 'size a' and 'size b', where the latter has a smaller initial QGP volume. The predictions are slightly lower than the data for both v 2 and v 3 , and the 'size a' curve is within two standard deviations of all points. However, in the right panel of Figure 6, the same calculation predicts a substantial suppression of high-p T hadrons, as expressed by the quantity R pPb = d 2 N pPb /dp T dy/(T pPb ×d 2 σ pp /dp T dy) where T pPb represents the nuclear thickness of the Pb nucleus, as determined via a Monte Carlo Glauber calculation [56]. Shown in comparison are published experimental results from ATLAS and ALICE for R pPb in central events that are consistent with no nuclear suppression, i.e. R pPb = 1 [41, 57]. The ALICE experiment uses the notation Q pPb for the same quantity to describe a bias that may exist due to the centrality categorisation. There are uncertainties in the experimental measurements related to the centrality or multiplicity selection in p+Pb collisions, particularly in determining the nuclear thickness value T pPb . However, there is no indication of the large R pPb suppression predicted by the jet quenching calculation. Thus, the jet quenching calculation is disfavoured as it cannot simultaneously describe the non-zero high-p T azimuthal anisotropy and the lack of yield suppression. The calculations are derived from two opposite limits of kinetic theory. The low momentum bands represent zeroth-order hydrodynamic calculations for high-multiplicity p+Pb events that give quantitative agreement with v 2 up to p T = 2 GeV while predicting values of v 3 that are too high. Above some high p T threshold, hadrons are expected to result, not from hydrodynamics, but instead from jets where the resulting partons have the opposite limit than in hydrodynamics, i.e. a large mean free path. To model this region, a non-hydrodynamic 'eremitic' expansion calculation (see Ref. [26] for the detailed calculation), shown as the bands at high p T , indicates slowly declining v 2 and v 3 coefficients. The dashed lines are a simple Padé-type fit connecting the two regimes [26]. The trends are qualitatively similar to those in the data, although there is not quantitative agreement. In particular, the calculation predicts values of v 2 and v 3 substantially below the experimental results for p T = 4-15 GeV. It should be noted that calculations presented in Ref. [26] are performed, consistently between the hydrodynamic and eremitic components, only for massless partons and with an ideal equation of state. Thus, one does not expect quantitative agreement and is looking for rather qualitative trends. More sophisticated treatments in the hydrodynamic regime result in better quantitative agreement with the anisotropy coefficients at low p T [58,59]. It is worth highlighting that traditional parton energy-loss calculations connect the high-p T v 2 with a suppression in the overall yield of high-p T particles. The same is true with this eremitic calculation, and thus, it should also be in contradistinction to p+Pb high-p T experimental data indicating almost no suppression, i.e. jet quenching.
Another possible source of the high-p T anisotropies could lie in an initial-state effect, potentially encoded in a model such as P 8. Shown in Figure 8 is a P 8 calculation with hard2 pp events overlaid on minimum-bias p+Pb events generated in the default Angantyr framework [60]. It is emphasised that this version of P does not include the recently developed string-string, or so-called string shoving, implementation [61]. The generator-level charged particles are then processed with the entire analysis procedure, including the non-flow template fit. The result is a negative v 2,2 for all momenta, in contradistinction to the experimental data. Further investigation reveals that P 8 run in 'hard' scattering mode has correlations with large pseudorapidity separation between particle pairs as a result of the specific implementation of initial-state radiation. This correlation is reduced in high-multiplicity 2 The term 'hard' refers to P 8 run with the following settings: HardQCD:all=on, PartonLevel:MPI=off, and containing a jet with p T > 100 GeV.  events because of the large number of uncorrelated UE particles, and thus results in a negative v 2,2 after subtracting the non-flow contribution. Figure 9 shows the published Pb+Pb results for v 2 as a function of p T in the 20-30% centrality selection [62] Figure 8: Predictions of azimuthal anisotropy from P 8 using the same two-particle formalism used for the data results. The events combine minimum-bias p+Pb underlying events generated in the Angantyr framework with hard pp events that require the presence of a jet with p T > 100 GeV. The two top plots show example correlation functions, with template fits, from a low particle-p T selection (top left) and a high particle-p T selection (top right). In the upper panels of the two top plots, the open circles show the scaled and shifted peripheral template with uncertainties omitted, the closed circles show the central data, and the red histogram shows the fit (template and harmonic functions). The blue dashed line shows the second-order harmonic component, Y compared to the v 2 from both the MBT p+Pb data and p+Pb containing a jet with p T > 100 GeV. This Pb+Pb centrality range is selected because the spatial elliptic eccentricity is approximately the same as in 0-5% centrality p+Pb collisions [63], despite having a much larger total particle multiplicity. The overall trends for Pb+Pb v 2 as a function of p T are qualitatively similar to those presented here for p+Pb from MBT events and the jet events with jet p T > 100 GeV. Both sets of the p+Pb values are scaled by a single multiplicative factor (1.5) to match the Pb+Pb rise at low p T . After scaling, the MBT p+Pb results quantitatively agree with those from the Pb+Pb system for 0.5 < p T < 8 GeV, except for a slight difference in the peak value near p T ≈ 3 GeV. For p T above about 8 GeV, the Pb+Pb results indicate a slow decline of v 2 values with increasing p T , while the p+Pb results exhibit more of a plateau. Strikingly, the overall behaviour of the v 2 values are quite similar.
As described above, the physics interpretations of the Pb+Pb elliptic anisotropies are hydrodynamic flow at low p T , differential jet quenching at high p T , and a transition between the two in the intermediate region of approximately 2 < p T < 10 GeV. Since these effects all relate to the initial QGP geometric inhomogeneities, a common shape with a single scaling factor for p+Pb could indicate a common physics interpretation. This scaling factor of 1.5, as empirically determined, may be the result of slightly different initial spatial deformations, or from the much larger Pb+Pb overall multiplicity, which enables a stronger translation of spatial deformations into momentum space. For the high p T region, this presents a conundrum in that it is difficult for differential jet quenching to cause the v 2 anisotropy in p+Pb collisions when there is no evidence for jet quenching overall. These measurements showing non-zero high p T v 2 in p+Pb collisions in the absence the jet quenching observed in Pb+Pb collisions suggest there might be additional contributions to v 2 at high p T in Pb+Pb collisions.  Figure 9: Scaled p+Pb v 2 values plotted as a function of the A-particle p T overlaid with v 2 from 20-30% central Pb+Pb data at √ s NN = 5.02 TeV [62]. Results from MBT and jet p T > 100 GeV p+Pb events are plotted as black squares and orange diamonds, respectively, and those from Pb+Pb are plotted as green circles. Statistical uncertainties are shown as narrow vertical lines on each point, and systematic uncertainties are presented as coloured boxes behind the points.
Returning to the issue of the difference in the intermediate p T region between the p+Pb MBT and jet event results, the source of hadrons in this region should be considered. As detailed previously, in a highly simplified picture one can classify hadrons as originating from hard scatterings (HS) or from the underlying event (UE). Thus, pairs of particles of A and B types can come from the combinations HS-HS, HS-UE, UE-HS, and UE-UE. Figure 10 presents the measured pair fractions for both MBT and jet, 0-5% central events plotted as a function of the A-particle p T . UE-UE pairs dominate the correlation functions at low p T in each case, and HS-UE combinations dominate at high p T . Combinations with HS B-particles are sub-dominant, because there are fewer jet particles than UE particles in central events; for the jet selected events, these combinations are further suppressed by the B-particle jet rejection condition. Figure 11 shows the dominant contributions from the MBT and jet events overlaid. Although the same qualitative behaviour is found in each case, the point at which the HS-UE pairs become dominant over the other combinations is at a lower p T for jet events than for MBT events.
This behaviour can also be seen in Figure 12, in which the pair fractions are plotted as a function of centrality, and again, the values for MBT and jet events are overlaid. The centrality-dependent results are plotted for low, medium, and high A-particle p T ranges in the same way as in Figure 5. At low p T , pair fractions from MBT and jet events agree, and in the mid-p T transition region, MBT events have a larger UE-UE contribution and smaller HS-UE contribution compared to jet events. At high p T , central events show a difference between UE-UE and HS-UE that is reduced in more-peripheral events and absent for more peripheral than 25% centrality. The overall trend of the pair fractions with centrality is quite similar to that of v 2 shown in Figure 5; little centrality dependence for low and high p T , and significant centrality dependence in addition to MBT-jet event ordering in the mid-p T transition region.
Thus, a potential explanation for the lower v 2 and v 3 in the intermediate p T region is simply that, in that region, the HS particles have lower anisotropy coefficients than UE particles, and MBT events have a larger fraction of UE-UE pairs than jet-triggered events. In the low and high p T regions, the same types of pairs dominate in both the MBT and jet-triggered events, namely UE-UE and HS-UE respectively, and hence the anisotropy coefficients agree between the event samples. If this explanation is correct, it also aids in understanding Figure 9 in which there is a significant difference between the p+Pb jet event v 2 and the Pb+Pb v 2 in the intermediate p T region, because the relative pair fractions are potentially significantly different.
This particle mixing picture is attractive in that it naturally explains the general shape of the v 2 (p T ) and v 3 (p T ) distributions as well as the ordering of the different event samples. However, it is noted that the correspondence between the differences in the flow coefficients and pair fractions is not quantitative; the differences in the flow coefficients are fractionally much larger than the differences in the pair fractions. Thus, there are either additional sources of correlation or our assumptions are violated in some way (e.g. the two assumed HS and UE sources are too simplistic or the measured pair fractions do not accurately represent the sources, as is discussed in Section 5). That said, for particle p T > 20 GeV, where particle production in any model is thought to arise mainly from jet fragmentation, the non-zero v 2 demonstrates that a positive correlation exists between hard (high p T ) and soft (low p T ) particles, irrespective of the pair fractions.

Conclusion
This paper presents Fourier coefficients of the azimuthal distribution of unidentified charged particles from 165 nb −1 of √ s NN = 8.16 TeV p+Pb collisions at the LHC and measured with the ATLAS detector. Results are presented separately for minimum-bias and jet events, with jet p T thresholds of 75 and 100 GeV, as a function of particle p T and centrality. The results are extracted using two-particle azimuthal correlations combined with a non-flow template fit procedure. The charged particle p T dependence of v 2 and v 3 is found assuming the factorisation of v n,n . The v 2 results are presented for charged-particle transverse momentum p T = 0.5-20 GeV for minimum-bias events and p T = 0.5-50 GeV for jet-triggered events, and the v 3 results are for p T = 0.5-20 GeV in both cases.
For charged particles with p T between 0.5 and 2 GeV, the v n results from each event selection are quantitatively consistent with each other, rising steadily with p T , and the v 2 coefficients are roughly independent of centrality. The v 2 values at 0-5% centrality agree with those predicted by hydrodynamic calculations.
Between charged particle p T of 2 and 9 GeV, the v n values drop in each case, but are ordered with minimum-bias events yielding the highest v n values and the jet events with jet p T > 100 GeV the lowest. Charged particles in this p T range exhibit a significant centrality-dependent v 2 , monotonically decreasing from central to peripheral events. This behaviour can be qualitatively explained within a simplified two-component model of particle production, in which the magnitude of the correlation in this region is determined by the admixture of charged particles originating from soft and hard processes in the given event selection. The measured particle pair yields support this qualitative argument.
For charged particles with p T above 9 GeV, the v n results are again consistent between MBT and jet events. Although the uncertainties in the v 3 values make any quantitative statement difficult, v 2 plateaus at a value of 0.025 up to a p T of 50 GeV. This result cannot be explained in the theoretical context of jet quenching or eremitic expansion calculations while simultaneously describing the observed lack of suppression of high-p T hadron and jet yields in p+Pb collisions.    [57] ATLAS Collaboration, Transverse momentum, rapidity, and centrality dependence of inclusive charged-particle production in √ s NN = 5.02    [61] C. Bierlich, [64] ATLAS Collaboration, ATLAS Computing Acknowledgements, ATL-GEN-PUB-2016-002, : https://cds.cern.ch/record/2202407.