Centrality dependence of dihadron correlations and azimuthal anisotropy harmonics in PbPb collisions at sqrt(s[NN]) = 2.76 TeV

Measurements from the CMS experiment at the LHC of dihadron correlations for charged particles produced in PbPb collisions at a nucleon-nucleon centre-of-mass energy of 2.76 TeV are presented. The results are reported as a function of the particle transverse momenta (pt) and collision centrality over a broad range in relative pseudorapidity [Delta(eta)] and the full range of relative azimuthal angle [Delta(phi)]. The observed two-dimensional correlation structure in Delta(eta) and Delta(phi) is characterised by a narrow peak at (Delta(eta), Delta(phi)) approximately (0, 0) from jet-like correlations and a long-range structure that persists up to at least |Delta(eta)| = 4. An enhancement of the magnitude of the short-range jet peak is observed with increasing centrality, especially for particles of pt around 1-2 GeV/c. The long-range azimuthal dihadron correlations are extensively studied using a Fourier decomposition analysis. The extracted Fourier coefficients are found to factorise into a product of single-particle azimuthal anisotropies up to pt approximately 3-3.5 GeV/c for at least one particle from each pair, except for the second-order harmonics in the most central PbPb events. Various orders of the single-particle azimuthal anisotropy harmonics are extracted for associated particle pt of 1-3 GeV/c, as a function of the trigger particle pt up to 20 GeV/c and over the full centrality range.


Introduction
Measurements of dihadron correlations are a well established technique for studying the properties of particle production in the high density medium created in heavy ion collisions. Early results from PbPb collisions at the Large Hadron Collider (LHC) [1,2] extended these studies into a regime of much higher beam energies as compared to those from the Relativistic Heavy Ion Collider (RHIC) [3][4][5][6][7][8][9][10]. These results complement other LHC measurements of medium properties, including a large deficit of charged particles at high-p T [11] and the observations of an enhanced fraction of dijets with very asymmetric energies [12,13].
The Compact Muon Solenoid (CMS) experiment at the LHC has studied dihadron correlations over a broad range of relative azimuthal angles (|∆φ|) and pseudorapidity (|∆η|, where η = − ln[tan(θ/2)] and θ is the polar angle relative to the counterclockwise beam axis) in the most central PbPb collisions at a nucleon-nucleon centre-of-mass energy ( √ s NN ) of 2.76 TeV [1].
Concentrating on large |∆η|, previous measurements at RHIC established some of the properties of the so-called "ridge" [4,6,9], an enhancement of pairs with |∆φ| ≈ 0. While a variety of theoretical models have been proposed to interpret the ridge phenomena as a consequence of jet-medium interactions [14][15][16][17][18][19], recent theoretical developments indicate that, because of event-by-event fluctuations in the initial shape of the interacting region, sizeable higher-order hydrodynamic flow terms could also be induced, e.g., triangular flow [20][21][22][23][24][25][26][27][28]. The triangular flow effect will contribute to the dihadron correlations in the form of a cos(3∆φ) component, which also gives a maximum near-side correlation at ∆φ ≈ 0, similarly to the elliptic flow contribution. It has been proposed that by taking into account various higher-order terms, the ridge structure could be described entirely by hydrodynamic flow effects [23]. To investigate this possibility, a Fourier decomposition of the CMS data at large |∆η| was performed, finding a strong dependence on p T [1]. Similar results, although with a smaller |∆η| gap, have been reported by ALICE [2,29]. The observations by CMS of a ridge-like structure in very high multiplicity proton-proton (pp) collisions at a centre-of-mass energy of 7 TeV [30], where no medium effect is expected, may also challenge the interpretations of these long-range correlations.
This paper, expanding on previous CMS results [1], presents dihadron correlation measurements from PbPb collisions at √ s NN = 2.76 TeV acquired in 2010, for all collision centralities and over a broader range of hadron p T . As in Refs. [1,30], the yield of particles (binned in p T ) associated with a trigger particle (also binned in p T ) is extracted as a function of their relative pseudorapidity and azimuthal angle. Such a study of hadron pairs in either the same or different p T ranges can reveal important information about the production of particles and their propagation through the medium. A Fourier decomposition technique is used to quantify the long-range azimuthal correlations. The potential connection between the extracted Fourier coefficients from the correlation data and the azimuthal anisotropy harmonics for single particles is investigated. This measurement provides a comprehensive examination of the centrality and transverse momentum (1 < p T < 20 GeV/c) dependencies of the short-range (|∆η| < 1) and long-range (2 < |∆η| < 4) dihadron correlations in PbPb collisions at LHC energies, as well as the relationship between these two-particle correlations and single-particle angular distributions. These results provide extensive input to the interpretation of these observables in terms of broad theoretical concepts such as hydrodynamic flow and quantitative models of particle production and propagation in the high-density medium.
The detector, the event selection and the extraction of the correlation functions are described in Sections 2 and 3, while the extracted results are described in Sections 4 and 5.

CMS Detector
The ability of the CMS detector (all components of which are described in Ref. [31]) to extract the properties of charged particles over a large solid angle is particularly important in the study of dihadron correlations. This study is based primarily on data from the inner tracker contained within the 3.8 T axial magnetic field of the large superconducting solenoid. The tracker consists of silicon pixel and strip detectors. The former includes 1 440 modules arranged in 3 layers, while the latter consists of 15 148 modules arranged in 10 (11) layers in the barrel (endcap) region. The trajectories of charged particles can be reconstructed for p T > 100 MeV/c and within |η| < 2.5.
The field volume of the solenoid also contains crystal electromagnetic and brass/scintillator hadron calorimeters. Although not included in the present results, muons are detected using gas-ionisation counters embedded in the steel return yoke. In addition to these components in and around the barrel and endcap of the solenoid, CMS also has extensive forward calorimetry. In the right-handed coordinate system used by CMS, the x-, y-, and z-axes are aligned with the radius of the LHC ring, the vertical direction, and the counterclockwise beam direction, respectively, with the origin located at the centre of the nominal interaction region.
For PbPb collisions, the primary minimum-bias trigger uses signals from either the beam scintillator counters (BSC, 3.23 < |η| < 4.65) or the steel/quartz-fibre Cherenkov forward hadron calorimeters (HF, 2.9 < |η| < 5.2). Coincident signals from detectors located at both ends of the detector (i.e., a pair of BSC or a pair of HF modules) are required. Events due to noise, cosmic-ray muons, double-firing triggers, and beam backgrounds are suppressed by further requiring the presence of colliding beam bunches. The fraction of inelastic hadronic PbPb collisions accepted by this primary trigger is (97 ± 3)% [13].

Data and Analysis
The procedure used in the present analysis follows that described in the previous CMS correlation paper [1]. Offline event selection requires a reconstructed vertex with at least two tracks (i.e., at least one pair of charged particles). This vertex must be within 15 cm along the beam axis relative to the centre of the nominal collision region and within 0.02 cm in the transverse plane relative to the average position of all vertices in a given data sample. In addition, various background events (for example beam-gas and beam-halo collisions, cosmic muons, and large-impact-parameter electromagnetic collisions) are suppressed by requiring at least three signals in the HF calorimeters at both positive and negative η, with at least 3 GeV of energy in each signal.
The analysis is based on a data sample of PbPb collisions corresponding to an integrated luminosity of approximately 3.9 µb −1 [32,33], which contains 30 million minimum-bias collisions after all event selections are applied. The pp data at √ s = 2.76 TeV, the reference for comparison to the PbPb data, were collected during a short low-energy LHC run at the end of March 2011. Minimum-bias-triggered pp events corresponding to an integrated luminosity of 520 µb −1 are selected for this analysis.
The energy released in the collisions is related to the centrality of the heavy ion interactions, i.e., the geometrical overlap of the incoming nuclei. The event centrality is defined as the fraction of the total cross section, starting at 0% with the most central collisions (i.e., smallest impact parameter). This fraction is determined from the distribution of total energy measured in both HF calorimeters.
The event centrality can be correlated with the total number of nucleons in the two Pb nuclei that experienced at least one inelastic collision, N part . The average values of N part for the various centrality bins used in this analysis are given in Table 1. The N part values are obtained using a Glauber Monte Carlo (MC) simulation [34,35] with the same parameters as in Ref. [13]. These calculations are translated into reconstructed centrality bins using correlations between N part and the measured total energy in the HF calorimeters, obtained from fully simulated MC events. The systematic uncertainties on the N part values in Table 1 are derived from propagation of the uncertainties in the parameters of the Glauber model. More details on the determination of centrality and N part can be found in Refs. [13,36,37]. The reconstruction of charged particles in PbPb collisions is based on signals in the silicon pixel and strip detectors, similarly to the reconstruction for pp collisions [38]. However, a number of settings are adjusted to cope with the challenges presented by the much higher signal density in central PbPb collisions. A set of tight quality selections are imposed on the collection of fully reconstructed tracks to minimise the contamination from misidentified tracks. These include requirements of at least 13 signals on the track, a relative momentum uncertainty of less than 5%, a normalised χ 2 of less than 0.15 times the number of signals, and transverse and longitudinal impact parameters of less than three times the sum in quadrature of the uncertainties on the impact parameter and the primary vertex position. Studies with simulated MC events show that the combined geometrical acceptance and reconstruction efficiency for the primary-track reconstruction reaches about 60% for the 0-5% most central PbPb collisions at p T > 2 GeV/c over the full CMS tracker acceptance (|η| < 2.4) and 66% for |η| < 1.0. The fraction of misidentified tracks is about 1-2% for |η| < 1.0, but increases to 10% at |η| ≈ 2.4 for the 5% most central PbPb collisions. For the peripheral PbPb events (70-80%), the overall tracking efficiency improves by up to 5%, with a much lower fraction of misidentified tracks.
The analysis of dihadron angular correlations in this paper follows exactly the procedure established in Ref. [1]. Any charged particle associated with the primary vertex and in the range |η| < 2.4 can be used as a trigger particle. A variety of bins of trigger transverse momentum, denoted by p trig T , are considered. There can be more than one such trigger particle in a single event and their total multiplicity in a particular data sample is denoted by N trig . Within each event, every trigger particle is then paired with all of the remaining particles (again within |η| < 2.4). As for the trigger particles, these associated particles are binned in transverse momentum (p assoc T ). The differential yield of associated particles per trigger particle is given by where N pair is the total number of correlated hadron pairs. The functions S(∆η, ∆φ) and B(∆η, ∆φ) are called the signal and background distributions, respectively. The value of the latter at ∆η = 0 and ∆φ = 0 (B(0, 0)) is a normalisation factor.
The signal distribution is the per-trigger-particle yield of pairs found in the same event, where N same is the number of such pairs within a (∆η,∆φ) bin. The background distribution is found using a mixed-event technique, wherein trigger particles from one event are combined (mixed) with all of the associated particles from a different event. In the analysis, associated particles from 10 randomly chosen events are used. The result is given by where N mix denotes the number of mixed-event pairs. This background distribution represents the expected correlation if the only effects present were random combinatorics and pairacceptance.
The value of B(∆η, ∆φ) at ∆η = 0 and ∆φ = 0 (with a bin width of 0.3 in ∆η and π/16 in ∆φ) is used to find the normalisation factor B(0, 0). In this case, the two particles have the maximum possible geometric pair acceptance since they are travelling in essentially the same direction. The effect of two tracks merging into a single reconstructed track is negligible. The extent to which the background distribution at larger angular separation is smaller than this value (more specifically the ratio B(0, 0)/B(∆η, ∆φ)) can be used to determine the pair acceptance correction factor. Multiplying the signal distribution by this ratio gives the acceptance-corrected pertrigger-particle associated yield. Since the distributions should, in principle, be symmetric, the statistical precision is maximised by filling only one quadrant using the absolute values of ∆η and ∆φ. For illustration purposes only (for example, see Fig. 1), the other three quadrants are filled by reflection, giving distributions that are symmetric about (∆η, ∆φ) = (0, 0) by construction. The pair acceptance decreases rapidly with ∆η and so, to avoid large fluctuations due to statistical limitations, the distributions are truncated at |∆η| = 4. The analysis is performed in twelve centrality classes of PbPb collisions ranging from the most central 0-5% to the most peripheral 70-80%. Within each centrality range, the yield described in Eq. 1 is calculated in 0.5 cm wide bins of the vertex position (z vtx ) along the beam direction and then averaged over the range |z vtx | < 15 cm.
When filling the signal and background distributions, each pair is weighted by the product of correction factors for the two particles. These factors are the inverse of an efficiency that is a function of each particle's pseudorapidity and transverse momentum, where A(η, p T ) is the geometrical acceptance, E(η, p T ) is the reconstruction efficiency, and F(η, p T ) is the fraction of misidentified tracks. The effect of this weighting factor only changes the overall scale but not the shape of the associated yield distribution, which is determined by the signal-to-background ratio.
As described in Ref. [1], the track-weighting procedure is tested using MC events generated with HYDJET [39] (version 1.6) propagated through a full detector simulation. The tracking efficiencies themselves are checked using simulated tracks embedded into actual data events. Systematic uncertainties due to variations of the track reconstruction efficiency as a function of vertex location and also the procedure used to generate the background events are evaluated. The individual contributions are added in quadrature to find the final systematic uncertainties of 7.3-7.6%.

4 Correlation Functions and Near-Side Yields
The two-dimensional (2D) per-trigger-particle associated yield distribution of charged hadrons as a function of |∆η| and |∆φ| is measured for each p trig T and p assoc T interval, and in different centrality classes of PbPb collisions. An example for trigger particles with 3 < p trig T < 3.5 GeV/c and associated particles with 1 < p assoc T < 1.5 GeV/c is shown in Fig. 1, for centralities ranging from the 0-5% most central collisions, to the most peripheral (70-80%) events. The 2D correlations are rich in structure, and evolve with centrality. The p trig T and p assoc T ranges shown in this figure were chosen as an example because they demonstrate a good balance of the following features. For the most central PbPb collisions, a clear and significant ridge-like structure mostly flat in ∆η, and extending to the limit of |∆η| = 4, is observed at ∆φ ≈ 0. At mid-peripheral events, a pronounced cos(2∆φ) component emerges, originating predominantly from elliptic flow [10]. Lastly, in the most peripheral collisions, the near-side ridge structure has largely diminished, while the away-side back-to-back jet correlations can be clearly seen at ∆φ ≈ π, but spread out in ∆η.  : Two-dimensional (2D) per-trigger-particle associated yield of charged hadrons as a function of |∆η| and |∆φ| for 3 < p trig T < 3.5 GeV/c and 1 < p assoc T < 1.5 GeV/c, for twelve centrality ranges of PbPb collisions at √ s NN = 2.76 TeV. The near-side peak is truncated in the two most peripheral distributions to better display the surrounding structure.
As was done in Ref. [1], to quantitatively examine the features of short-range and long-range azimuthal correlations, one-dimensional (1D) ∆φ correlation functions are calculated by aver-aging the 2D distributions over a limited region in ∆η from ∆η min to ∆η max : The results of extracting the 1D ∆φ correlations in the short-range (0 < |∆η| < 1) and longrange (2 < |∆η| < 4) regions are shown in Fig. 2. The associated yield distribution per trigger particle is extracted for the same p trig T and p assoc T ranges as in Fig. 1. In order to study the short-range ∆φ correlations in the absence of the flat background in ∆η, the 1D ∆φ distribution in the long-range region is subtracted from that in the short-range region. The resulting difference of the distributions is shown in Fig. 3. The near-side peak (∆φ ≈ 0) represents mainly the correlations from jet fragmentation, whereas the away-side region (∆φ ≈ π) is mostly flat and close to zero due to the weak ∆η dependence of the away-side jet peak. A comparison to the pp data at √ s = 2.76 TeV is also presented, showing a similar structure to that in the very peripheral 70-80% PbPb data. However, the magnitude of the near-side peak is significantly enhanced in the most central PbPb collisions as compared to pp. Most of the systematic uncertainties manifest themselves as an overall change in the scale of the correlation functions, with little dependence on ∆φ and ∆η. Therefore, they largely cancel when the difference between the short-range and long-range regions is taken. The difference between short-range (0 < |∆η| < 1) and long-range (2 < |∆η| < 4) per-trigger-particle associated yields of charged hadrons as a function of |∆φ| for 3 < p trig T  and the values one would expect from a trivial extrapolation of the pp results, the latter are also represented using horizontal lines covering the full N part range. The yield of the near-side peak increases by a factor of 1.7 in going from the very peripheral 70-80% to the most central 0-5% PbPb events, for the lowest p assoc T interval of 1-1.5 GeV/c. As p assoc T increases, the centrality dependence of the near-side yield becomes less prominent. An increase by a factor of only 1.3 is observed for the highest p assoc T interval of 2.5-3 GeV/c. This is of particular interest because at RHIC energies for p assoc T down to 2 GeV/c and similar p trig T ranges and methodology, although for a lower density system (AuAu at √ s NN = 0.2 TeV), there is almost no centrality dependence observed [40]. On both near and away sides, the yield in PbPb matches that in pp for the most peripheral events. On the away side, the yield in PbPb decreases with centrality, becoming negative for the most central events. Variations in the event-mixing procedure can cause large fluctuations but only at the very edge of the acceptance around |∆η| = 4.8. However, the correlation function is only studied up to |∆η| < 4 in this paper so these fluctuations do not affect the results. The negative values of the yields in Fig. 4 are caused by a slightly concave structure on the away-side region (1.18 < |∆φ| < π), i.e. the yields near |∆η| ≈ 0 are smaller than those at higher |∆η|. The effect is more prominent for central PbPb events. However, this concavity is seen only for |∆η| < 2. Beyond that region, the ∆η distribution is found to be largely flat up to |∆η| = 4. Similar behaviour was also observed at RHIC for AuAu collisions at √ s NN = 200 GeV [41]. This deviation from pp may be related to the jet quenching phenomena, which leads to a modification in the back-to-back jet correlations in PbPb. Any effect that modifies the kinematics of dijet production could also result in a modification of away-side distributions in ∆η. Additionally, any slight dependence of the flow effect on η could also play a role. More detailed theoretical models will be required to fully understand the origin of this small effect.

Fourier Decomposition Analysis of the PbPb Data
The first Fourier decomposition analysis of long-range dihadron azimuthal correlations for PbPb collisions at √ s NN = 2.76 TeV was presented in Ref. [1]. This analysis was motivated by the goal of determining whether the long-range ridge effect was caused by higher-order hydrodynamic flow harmonics induced by the initial geometric fluctuations. This decomposition involves fitting the 1D ∆φ-projected distribution for 2 < |∆η| < 4 (to avoid the jet peak) with a Fourier series given by where V n∆ are the Fourier coefficients and N assoc represents the total number of hadron pairs per trigger particle for the given |∆η| range and (p trig T , p assoc T ) bin. The first five Fourier terms (N max = 5) are included in both the current fits and those in Ref. [1]. In this paper, the analysis of the Fourier decomposition is extended to the full centrality range, and is performed as a function of both p trig T and p assoc T .
The Fourier decomposition results have several systematic uncertainties. Because the trackingcorrection-related systematic uncertainties only change the overall scale of the correlation functions, instead of the shape, they have only a ±0.8% uncertainty on the extracted Fourier coefficients (V n∆ ), largely independent of n and collision centrality. In addition, the results are insensitive to looser or tighter track selections to within ±0.5%. By comparing the Fourier coefficients derived for two different z vtx ranges, |z vtx | < 15 cm and |z vtx | < 5 cm, the systematic uncertainties due to the dependence on the vertex position are estimated to be less than ±0.5%. Variations from the finite bin width of the ∆φ histograms contribute the largest systematic uncertainty to the analysis, especially for the higher-order components, which are more sensitive to the fine structure of the distributions. Reducing the binning of the ∆φ histograms by factors of 2, 4, and 8, the extracted Fourier coefficients vary by ±0.3-2.2%. The effect of including additional higher-order Fourier terms in the fit using Eq. 6 results in changes of at most ±1.0% (for n = 5), with Fourier terms up to n = 10 included (N max = 10). The values of additional higher-order Fourier terms included in the fit are all consistent with zero. Table 2 summarises the different sources of uncertainty for the first five Fourier coefficients. These uncertainties are added in quadrature to obtain the total systematic uncertainties, also given in Table 2. The fitted Fourier coefficients (V n∆ ) up to n = 5, for two representative low-p trig T ranges of 1 < p trig T < 1.5 GeV/c and 3 < p trig T < 3.5 GeV/c, with p assoc T fixed at 1-1.5 GeV/c, are presented in Fig. 5 for various centrality ranges. The values of V n∆ peak at n = 2 and then drop dramatically toward larger n values at all centralities, although this behaviour is less pronounced for the 0-5% centrality bin. The error bars show the statistical uncertainties only, while the systematic uncertainties are indicated in Table 2.  Table 2.

Factorisation of Fourier Coefficients
If the observed azimuthal dihadron correlations at large ∆η are driven only by the singleparticle azimuthal anisotropy with respect to a particular direction in the event, the extracted Fourier coefficients (V n∆ ) from long-range azimuthal dihadron correlations can be factorised into a product of the single-particle azimuthal anisotropy harmonics, v n , via where v n (p trig T ) and v n (p assoc T ) are the harmonics for the trigger and associated particles [2] averaged over all the events, respectively. One source of v n is the collective-flow harmonics arising from hydrodynamic expansion of the medium (e.g., anisotropic elliptic flow contribution to v 2 ) [42], particularly in the low-p T regime where hadron production in heavy ion collisions is thought to be mainly from the bulk medium [43]. In addition, for very high p T particles that are predominantly produced by the fragmentation of energetic jets, v n could also be induced by the path-length dependence of the jet-quenching effect inside the medium [44][45][46][47][48][49]. This path difference can lead to a stronger suppression of the high-p T hadron yield along the long axis of the elliptically-shaped system than along its short axis, resulting in an azimuthal anisotropy characterised by the v 2 harmonic. Both scenarios satisfy the factorisation relation of Eq. 7. However, note that the p T dependent event-by-event fluctuations of v n could break the factorisation in general, even though v n may be still related to the single-particle azimuthal anisotropy. This possibility is not investigated in this paper.
This relation (Eq. 7) is a necessary ingredient for the extraction of single-particle azimuthal anisotropy harmonics using the dihadron correlation data, since a Fourier series can be used to decompose any functional form by construction. The relation can be tested by first assuming that factorisation is valid for pairs including one particle in a fixed, low p assoc T range, denoted by p low T , correlated with a second particle of any p T . The range of p low T is chosen to be 1-1.5 GeV/c, where particle production is expected to be predominantly driven by hydrodynamics. The value of v n (p low T ) is first calculated as the square root of V n∆ (p low This is effectively equivalent to the two-particle cumulant method of flow measurement [50,51]. Next, the ratio of V n∆ (p  Figures 6, 7, and 8 show the ratios for n = 2, 3, and 4, over five p assoc T ranges as a function of p trig T , for both short-range (0 < |∆η| < 1) and long-range (2 < |∆η| < 4) regions. Three different centrality intervals 0-5%, 15-20% and 35-40% are presented. The ratio for n = 2 in pp data is also shown in the last column of Fig. 6. The first point of each panel equals 1.0 by construction. The error bars correspond to the statistical uncertainties. The total systematic uncertainties are estimated to be 1.5% (n = 2)-3.6% (n = 5), approximately √ 2 times the systematic uncertainties of V n∆ shown in Table 2.
First of all, no evidence of factorisation is found in the pp data and for the short-range region (0 < |∆η| < 1) in any of the centrality ranges of the PbPb data, where dijet production is expected to be the dominant source of correlations. In Ref. [2], it has been shown that V n∆ factorises for jet-like correlations at very high-p T (e.g., p T > 5 GeV/c) as the direction of the dijet forms a special axis, to which produced particles are strongly correlated. This is similar to the elliptic flow effect, where particles are preferentially produced along the short axis of the elliptically-shaped overlapping region. However, the lack of factorisation in pp and the shortrange region of PbPb observed for the p T range of 1 < p assoc T < 3.5 GeV/c primarily investigated in this paper may suggest a complicated interplay of different particle production mechanisms between low-p T (hydrodynamic flow for PbPb and underlying event for pp) and high-p T (dijet production) particles.
In contrast, the long-range region (2 < |∆η| < 4) for mid-peripheral 15-20% and 35-40% events does show evidence of factorisation for V 2∆ to V 4∆ with p assoc T up to 3-3.5 GeV/c and p trig T up to approximately 8 GeV/c. The data are also consistent with factorisation for even higher p trig T (> 8 GeV/c) combined with low p assoc T , but the current event sample is not large enough to draw a firm conclusion. Note that V n∆ varies by almost 60% in the p trig T range from 1 to 3.5 GeV/c as shown in Fig. 5, whereas factorisation is found to hold to better than 5%. This suggests a potential connection between the extracted Fourier coefficients from long-range dihadron correlations and the single-particle azimuthal anisotropy harmonics. For the most central 0-5% collisions, the ratio for n = 2 deviates significantly from unity, while V 3∆ and V 4∆ still show a similar level of factorisation to that of the 15-20% and 35-40% mid-peripheral data. This may indicate the existence of other sources of long-range correlations for the most central collisions that violate the factorisation relation. The breakdown of factorisation for p assoc T > 4 GeV/c in the long-range region is likely due to dijet correlations. Higher-order Fourier coefficients and a wider p T range can be investigated once larger samples of PbPb data are collected. Factorisation of V 1∆ is not discussed in this paper as it contains an additional negative contribution from momentum conservation [52], which is not related to the collective flow effect and needs to be accounted for in further studies.

Elliptic and Higher-Order Single-Particle Azimuthal Anisotropy Harmonics
As discussed in Section 5.1, except for v 2 in the very central PbPb events, the factorisation relation given by Eq. 7 for long-range (2 < |∆η| < 4) azimuthal dihadron correlations is found to be valid for sufficiently low p assoc T , combined with low p trig T , and possibly high p trig T as well. Therefore, the single-particle azimuthal anisotropy harmonics v n (p trig T ) can be extracted using Eq. 8 with 1 < p assoc T < 3 GeV/c. Values are found for centralities ranging from 0-5% to 70-80%, and presented in Fig. 9. The 1-3 GeV/c p assoc T range is chosen in order to reduce the statistical uncertainty by utilising as many associated particles as possible over the p assoc T range where factorisation is valid. Data for the most central and most peripheral events are included for completeness, although the results for v 2 in those events are clearly demonstrated by Fig. 6 to be more complicated in nature. The value of v 2 is extracted up to p trig T ∼ 20 GeV/c for all but the 2 most peripheral centralities, whereas higher-order v n are truncated at p trig T ∼ 10 GeV/c or less for the peripheral data due to statistical limitations. As mentioned previously, factorisation is not demonstrated conclusively at very high p trig T . For the most central 0-5% events, all the harmonics are of similar magnitude across the entire p trig T range. The p T dependence of all v n shows the same trend of a fast rise to a maximum around p T ≈ 3 GeV/c, followed by a slower fall, independent of centrality up to 50-60%. The magnitude of v 2 increases when moving away from the most central events. At very high p trig T , sizeable v 2 signals are observed, which exhibit an almost flat p T dependence from 10 to 20 GeV/c for most of the centrality ranges. This is not the case for the higher-order harmonics. In order to explicitly investigate the centrality (GeV/c) trig T p 2 4 6 8 1.0 1.0 1.0  Table 2.
A strong centrality dependence of v 2 is observed in Figs. 9 and 10 for all p trig T ranges, while the higher-order harmonics v 3 -v 5 do not vary significantly over the range of N part . This behaviour is expected in the context of both the hydrodynamic flow phenomena for lower-p T particles [43] and the path-length dependence of the parton energy-loss scenario for high-p T particles [44].
The v 2 harmonics are sensitive to the eccentricity of the almond-shaped initial collision region that becomes larger for the peripheral events, whereas the higher-order harmonics are driven by fluctuations of the initial geometry that have little dependence on the collision centrality [53]. In the most peripheral events, the p T dependence of v 2 is found to be very different from that in the central events, as shown in Fig. 9. In the high-p T interval 8-20 GeV/c (third panel of Fig. 10), v 2 increases rapidly at low N part (very peripheral). A possible explanation is the presence of non-flow effects due to back-to-back jets. Based on the factorised Fourier coefficients from long-range dihadron correlations, the single-particle azimuthal anisotropy harmonics extracted over a wide range of p T and centrality allow a detailed comparison to theoretical calculations of the hydrodynamics and path-length dependence of in-medium parton energy loss.  Table 2.

Summary
The previous CMS analysis of angular correlations between charged particles has been expanded to cover a wide centrality range of PbPb collisions at √ s NN = 2.76 TeV. As was seen previously for central PbPb collisions, the associated yields differ significantly from those observed in pp interactions. Correlations with both small (0 < |∆η| < 1) and large (2 < |∆η| < 4) relative pseudorapidities were again studied as a function of the transverse momentum of the trigger and associated particle p T . The integrated yield of the near-side region shows an increasing enhancement towards more central PbPb collisions, especially for low-p T associated particles.
To further characterise the dependence of the correlations on relative azimuthal angle, a Fourier decomposition of the distributions projected onto ∆φ was performed, as a function of both centrality and particle p T . Evidence of a factorisation relation was observed between the Fourier coefficients (V n∆ ) from dihadron correlations and single-particle azimuthal anisotropy harmonics (v n ). This holds for p assoc T 3-3.5 GeV/c over the p trig T range up to at least 8 GeV/c in central and mid-peripheral PbPb collisions, except for v 2 in the most central events. The observed factorisation is absent in pp and very peripheral PbPb data, indicating a strong connection of the observed long-range azimuthal dihadron correlations to the single-particle azimuthal anisotropy in heavy ion collisions, such as the one driven by the hydrodynamic expansion of the system. The single-particle azimuthal anisotropy harmonics v 2 through v 5 were extracted over a wide range in both p T and collision centrality, profiting from the broad solid-angle coverage of the CMS detector. The comprehensive correlation data presented in this paper are very useful for studies of the path-length dependence of in-medium parton energy-loss, and provide valuable inputs to a variety of theoretical models, including hydrodynamic calculations of higher-order Fourier components.