Measurement of flow harmonics correlations with mean transverse momentum in lead–lead and proton–lead collisions at sNN=5.02TeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{\mathrm{NN}}} = 5.02~\hbox {TeV}$$\end{document} with the ATLAS detector

To assess the properties of the quark–gluon plasma formed in ultrarelativistic ion collisions, the ATLAS experiment at the LHC measures a correlation between the mean transverse momentum and the flow harmonics. The analysis uses data samples of lead–lead and proton–lead collisions obtained at the centre-of-mass energy per nucleon pair of 5.02 TeV, corresponding to total integrated luminosities of 22μb-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$22~\upmu \text {b}^{-1}$$\end{document} and 28nb-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$28~\text {nb}^{-1}$$\end{document}, respectively. The measurement is performed using a modified Pearson correlation coefficient with the charged-particle tracks on an event-by-event basis. The modified Pearson correlation coefficients for the 2nd-, 3rd-, and 4th-order flow harmonics are measured in the lead–lead collisions as a function of event centrality quantified as the number of charged particles or the number of nucleons participating in the collision. The measurements are performed for several intervals of the charged-particle transverse momentum. The correlation coefficients for all studied harmonics exhibit a strong centrality evolution, which only weakly depends on the charged-particle momentum range. In the proton–lead collisions, the modified Pearson correlation coefficient measured for the 2nd-order flow harmonics shows only weak centrality dependence. The lead-lead data is qualitatively described by the predictions based on the hydrodynamical model.


Introduction
The large azimuthal anisotropy observed for particles produced in heavy-ion collisions at RHIC [1][2][3][4] and the LHC [5][6][7][8] is one of the main signatures of the formation of strongly interacting matter called quark-gluon plasma (QGP). A standard picture of an ultrarelativistic heavy-ion collision is that the initial, asymmetric 'almond' shape of the colliding nuclei's overlap region leads to the formation of pressure gradients in the QGP. These pressure gradients transform the e-mail: atlas.publications@cern.ch initial shape into an azimuthal anisotropy of the final-state particle distributions through a nearly ideal hydrodynamic evolution and subsequent QGP hadronisation process [9]. The azimuthal anisotropy is customarily decomposed into Fourier components with the amplitude of the nth term denoted by v n and known as a flow harmonic [10]. Theoretical hydrodynamical models successfully describe observed flow phenomena at low particle transverse momenta [11]. The properties of QGP were recently studied with measurements of correlations between flow harmonics of different order [12][13][14][15][16] as well as with analyses of event shapes [16][17][18][19][20]. It is expected that in lead-lead (Pb+Pb) collisions the magnitudes of the azimuthal flow harmonics [6,7] should be correlated with the mean transverse momentum [ p T ] of the particles on an event-by-event basis [21]. In this paper, that correlation is called the v n -[ p T ] correlation. In proton-lead ( p+Pb) collisions, the measurements of multiparticle correlations [22] show evidence of collective phenomena. The spectra of identified particles in p+Pb collisions are consistent with a presence of the radial flow [23] while the nuclear modification factor at high p T approaches unity [24]. Despite intensive studies, the mechanism responsible for the collective behaviour in small collision systems still remains unknown [9]. In p+Pb collisions the v n -[ p T ] correlation could provide constraints on the initial geometry of the particle source, thereby reducing the overall modelling uncertainty. According to the hydrodynamical model predictions [25], in p+Pb collisions the v n -[ p T ] correlation is sensitive to the distribution of energy deposition in the first stage of the collision. For a larger source a positive v 2 -[ p T ] correlation is expected while for a compact source the negative correlation is obtained. Simultaneous measurements of v n -[ p T ] correlations in small and large systems may help disentangle the role of initial conditions and subsequent dynamical QGP evolution in final-state particle distributions.
To measure the strength of the v n -[ p T ] correlation, the Pearson correlation coefficient (PCC) R [25] is used where The term v n {2} 2 is the square of the nth-order flow harmonic obtained using the two-particle correlation method [26], cov(v n {2} 2 , [ p T ]) is the covariance between v n {2} 2 and [ p T ], and Var(v n {2} 2 ) and Var([ p T ]) are the variances of the v n {2} 2 and [ p T ] distributions, respectively. Experimentally, however, the finite event-by-event charged-particle track multiplicity results in an additional broadening of the v n {2} 2 and [ p T ] distributions due to statistical fluctuations. Thus, the values of the respective variances are increased, especially for [ p T ]. The magnitude of this broadening depends on the choice of kinematic region and on detector performance, making direct comparisons between experimental results and with theoretical calculations difficult. To overcome this problem, a modified correlation coefficient ρ, less sensitive to the charged-particle multiplicity than R, was suggested in Ref. [25]. To reduce the auto-correlation effects and those due to the finite charged-particle multiplicity in an event, the variances of the v n {2} 2 and [ p T ] distributions are replaced by corresponding dynamical variables, which are more sensitive to intrinsic initial-state fluctuations. The variance of v n {2} 2 is replaced by its dynamical counterpart [27] Var where corr n {2} and corr n {4} are the two-and four-particle correlations [26] and where angular brackets denote that they are averaged over events. These correlations are described in detail in Sect. 4. The variance of [ p T ] is replaced by the dynamical p T fluctuation magnitude [28,29] c k defined as where [ p T ] is the average [ p T ] over the all analysed events. The modified PCC ρ is thus defined as It was demonstrated in Ref.
[25] that the ρ coefficient calculated using realistic and finite multiplicities provides a reliable estimate of the true value of R found in the limit of infinite multiplicity, whereas the coefficient R, calculated using Eq. (1) for finite multiplicity underestimates the true value.
The ALICE experiment measured [20] that the chargedparticle p T spectrum is correlated with the magnitude of the elliptic flow. It is measured to be harder in collisions with the higher second flow harmonics and softer in collisions where the elliptic flow is smaller. The magnitude of spectra modification is observed to increase with p T , starting to be significant at around 1 GeV and reaching a few percent at around 5 GeV. The modification is found to be most significant in the mid-central collisions, decreasing in the most central ones. The ALICE results suggest that the value of the correlation coefficient should be significant in mid-central and central collisions and that its magnitude and centrality dependence should be sensitive to the scale of intrinsic fluctuations of v 2 and p T . Including particles of higher p T in the measurement is expected to result in increased values of the Pb+Pb collisions,v 3 and v 4 in wide centrality range as well as for the v 2 in high multiplicity p+Pb are unexplored by measurements.
This paper reports on the first measurement of the ρ coefficient with the ATLAS detector in Pb+Pb and p+Pb collisions at a centre-of-mass energy per nucleon pair of 5.02 TeV. The Pb+Pb data sample with a total integrated luminosity of 22 μb −1 was collected in 2015, and the p+Pb sample with 28 nb −1 in 2013.
This paper is organised as follows. Section 2 gives a brief description of the ATLAS detector. Details of the event selection and charged-particle reconstruction are provided in Sect. 3. Section 4 describes the analysis procedure for calculating the ρ coefficient. Systematic uncertainties are described in Sect. 5 and Appendix A. Results are presented in Sect. 6, followed by a summary in Sect. 7.

Experimental setup
The ATLAS experiment [30] at the LHC is a multipurpose particle detector with a forward-backward symmetric cylindrical geometry and a near 4π solid angle coverage. The inner detector (ID) covers the pseudorapidity 1 range |η| < 2.5 and is surrounded by a thin superconducting solenoid providing a 2 T axial magnetic field. The ID consists of silicon pixel, silicon microstrip (SCT), and straw tube tracking detectors. After the 2013 p+Pb run, an additional pixel silicon layer, the insertable B-layer [31,31,32], was installed prior to the 5.02 TeV Pb+Pb data-taking to attain more precise track-ing. Lead/liquid-argon (LAr) sampling calorimeters provide electromagnetic (EM) energy measurements with high granularity. A steel/scintillator tile hadronic calorimeter covers the central pseudorapidity range (|η| < 1.7). The endcap and forward regions are instrumented with LAr calorimeters for EM and hadronic energy measurements up to |η| = 4.9. The forward calorimeter (FCal) covers 3.2 < |η| < 4.9 and is used for centrality estimation [10]. The minimumbias trigger scintillators (MBTS) are located on each side of the detector at z = ±3.6 m and detect charged particles with 2.07 < |η| < 3.86. The zero-degree calorimeter (ZDC), located in the LHC tunnel and covering |η| > 8.3, is used for triggering on collision events and pile-up event rejection. It is calibrated to resolve an individual neutron originating from the collision spectators.
A two-level trigger system selects events [33,34]. The level-1 trigger is implemented in hardware and preselects up to 10 5 events per second for further decisions by the high-level trigger (HLT). The software-based HLT tuned for Pb+Pb collision data selects up to 1000 events per second for recording. This analysis primarily uses charged-particle tracks in the ID, but information from the central calorimeters and the ZDC is also used for triggering, event selection, and analysis.

Event and track selection
The Pb+Pb data in this analysis were selected using two mutually exclusive minimum-bias triggers. Events with semi-central and central collisions were selected if the scalar sum of transverse energy in the entire ATLAS calorimeter system exceeded 50 GeV. Peripheral events, i.e. those with large impact parameter of the colliding Pb nuclei, fail the 50 GeV selection and were instead selected by requiring a deposition in the ZDC corresponding to at least one neutron and by requiring at least one track reconstructed in the HLT. Data in this analysis are required to come from periods when the entire detector was functioning normally. The events are required to have a reconstructed vertex within 100 mm of the nominal interaction point. The contribution from events containing more than one inelastic interaction (pile-up) is studied by exploiting correlations between the transverse energy measured in the FCal ( E FCal T ) with the estimated number of neutrons in the ZDC, and with the number of tracks associated with a primary vertex [27,35]. The distribution of E FCal T and the distribution of the number of neutrons in events with more than one collision are broader than the corresponding distributions in events with only one collision. Pile-up events are suppressed by rejecting events with abnormally large values of either E FCal T or the number of neutrons in the ZDC compared with the charged-particle multiplicity in the event. Approximately 0.2% of the events are rejected by these requirements.
The p+Pb data in this analysis were selected using minimum-bias triggers and high-multiplicity triggers (HMT). The minimum-bias trigger required signals in both sides of the MBTS system with a timing difference of less than 10 ns to eliminate non-collision backgrounds. The HMT required the total transverse energy in the calorimeter at level-one and the number of ID track candidates reconstructed in the HLT to be above predefined thresholds. Six combinations of thresholds were used to optimise data-taking during periods with different luminosities. Samples of events collected by these triggers are combined by applying event weights to reproduce the charged-particle multiplicity distribution of the minimum-bias trigger. Further details of the data selection are given in Refs. [22,36]. The average pile-up probability in the p+Pb dataset is approximately 3% but can be significantly larger in high-multiplicity events. Events with more than one reconstructed vertex are removed from the sample. Similarly to the Pb+Pb dataset, to remove events where the two interaction vertices are too close to resolve as independent ones, the ZDC signal on the Pb fragmentation side is used. The distribution of the number of neutrons, which is broader in events with pile-up than that for the events without pile-up is exploited for that purpose [36]. The fraction of rejected events varies with the event activity and reaches a maximum of 10% for events with the highest multiplicities.
The analysis for both collision systems is performed in narrow bins of event activity defined by the charged-particle multiplicity N ch (described in Sect. 4), which estimates the collision centrality. In addition, the Pb+Pb results are presented as a function of collision centrality expressed by the average number of nucleons participating in the collision, N part , to allow comparison with theoretical predictions [37]. The centrality is estimated from the E FCal T distribution [6,10] using the Glauber model [38]. The number of events passing the selection requirements is 1.3 × 10 8 for Pb+Pb within the 0-80% centrality interval. For the p+Pb system, about 0.64 × 10 8 events enter the analysis.
The charged-particle tracks reconstructed in the ID are required to satisfy selection criteria in order to suppress the contribution of incorrectly reconstructed tracks and secondary products of particle decays. The selection criteria include the requirement that the number of hits in the pixel and SCT detectors should be greater than two and eight, respectively, for the Pb+Pb data and greater than one and six for the p+Pb data. The track impact parameters relative to the collision vertex in the transverse direction, |d 0 |, and longitudinal direction, |z 0 sin θ |, are required to be less than 1 mm for tracks in the Pb+Pb data sample and less than 1.5 mm in the p+Pb sample. In addition, in p+Pb collisions, the track impact parameter significances must satisfy |d 0 /σ d 0 | < 3 and |z 0 sin θ/σ z | < 3, where σ d 0 and σ z are the uncertainties in d 0 and z 0 sin θ determined from the covariance matrix of the track fit. The different selection criteria for Pb+Pb and p+Pb optimise the performance of the track reconstruction in differing running conditions.
Corrections needed due to track reconstruction effects are evaluated using 4×10 6 Pb+Pb and 10 7 p+Pb minimumbias Monte Carlo (MC) events generated by the HIJING v1.38b [39] event generator. After the generation, an azimuthal flow is implemented using the afterburner technique [40], and the p T spectrum is reweighted to match the data. Generated events were simulated in the detector by the Geant 4-based [41] ATLAS detector simulation programs [42] and reconstructed using the same procedures and detector conditions as the data. Track reconstruction corrections are applied to each selected track using weights to account for the tracking efficiency and the fake-track fraction f . The efficiency is defined as the fraction of primary MC charged particles that are matched to reconstructed tracks, and f is the fraction of tracks that are not matched to primary MC particles or are produced from random combinations of hits in the ID. A similar analysis procedure is described in Refs. [10,16]. The fake-track fraction and tracking efficiency are determined as functions of the track p T and η and of the track multiplicity in the event. Tracks included in the analysis are weighted with the factor (1− f )/ . An additional multiplicative weight evaluated from data is applied to the data to correct for detector non-uniformity in the azimuthal angle. These weights are obtained by requiring the tracks to be distributed uniformly in azimuth in all pseudorapidity slices of width 0.1.
In the Pb+Pb data, the contribution of fake tracks is largest in central collisions at the lowest analysed track p T of 0.5 GeV and at the largest |η|, reaching up to 20%. The faketrack rate is below 1% for tracks with p T above 2 GeV and |η| < 1.5. The tracking efficiency depends weakly on centrality, and in the most central events it is about 3% less than in more peripheral events. The efficiency increases with the track p T from about 50% at the lowest analysed p T to 70% above 2 GeV. It is highest at mid-rapidity and drops by about 15% for |η| > 1. For p+Pb collisions, with p T increasing from 0.3 to 1 GeV the efficiency increases from about 75% (60%) to 83% (70%) at η ≈ 0 (|η| > 2). The p+Pb tracking efficiency is independent of the event's multiplicity for N ch ≥ 10, i.e. in the multiplicity range used in the analysis. The fake rate in p+Pb collisions is very low, below 1% (3%) at η ≈ 0 (|η| > 2).

Correlation coefficient ρ
In each event, charged-particle tracks are grouped into three regions of subevents based on their pseudorapidity: region A with −2.5 < η < −0.75, central region B with |η| < 0.5 and region C with 0.75 < η < 2.5. The v 2 n for the n = 2-4 harmonics are calculated by correlating chargedparticle tracks from subevents A and C, which are separated in pseudorapidity to suppress non-flow contributions. Tracks in central region B are used to obtain the mean value of the charged-particle transverse momentum in the event, [ p T ], defined as where the summation is over tracks in region B, labelled by index b. The variable c k (Eq. (3)) is also calculated using tracks from region B. Here, and in following formulas, the weights w include the fake-track fraction, efficiency, and azimuthal non-uniformity corrections, as discussed in Sect. 3. The covariance term from the numerator of Eq. (4) is defined as where φ is the azimuthal angle and indices a and c span the tracks in regions A and C, respectively. The two-and four-particle correlations used to define the dynamical variance in Eq. (2), which enters the denominator of Eq. (4), are calculated as in Ref. [26] where the q a and q c are the complex flow vectors of subevent A and subevent C, respectively, and the asterisk denotes the complex conjugate. The flow vectors are q n,a = 1 a w a a w a e inφ a and q n,c = 1 The four-particle correlation is obtained from the expression where for subevent A and similarly for subevent C. Equation (7) represents the sum e in(φ a 1 +φ a 2 −φ c 3 −φ c 4 ) over all particles from subevents A and C normalised by the number of quadruplets without autocorrelations in each subevent.
The second factor in the denominator of Eq. (4), the mean p T fluctuation in the event class c k , is defined by Eq. (3) and in this analysis it is calculated as The summation indices b and b run over all charged particles in region B.
The correlation coefficient expressed by Eq. (4) is evaluated for the range 0.5 < p T < 2 GeV in Pb+Pb collisions and 0.3 < p T < 2 GeV in p+Pb collisions. These intervals, called 'main', contain a large number of soft particles and constitute the main result of the analysis which can be compared with hydrodynamical models. For each system, two additional p T ranges are considered: 0.5 < p T < 5 GeV and 1 < p T < 2 GeV in the analysis of Pb+Pb collisions, and 0.3 < p T < 5 GeV and 0.5 < p T < 2 GeV in p+Pb collisions. These ranges facilitate the study of the sensitivity of ρ(v n {2} 2 , [ p T ]) to the high p T part of the particle spectrum and to the lower charged-particle multiplicity from the higher minimum p T value. The charged-particle p T range 0.5 < p T < 2 GeV is common to both systems and can be used to compare the ρ(v 2 {2} 2 , [ p T ]) results from Pb+Pb and p+Pb collisions.
, are determined in bins of reconstructed track multiplicity M AC measured in the combination of regions A and C. This is done to avoid a negative correlation between the multiplicity in subevents A+C and B that occurs if the analysis is binned in multiplicity in the entire ID. Narrow M AC bins are also chosen due to the sensitivity to multiplicity fluctuations of the multi-particle correlations that are used to obtain the Var v n {2} 2 dyn [27]. The events are grouped in fine bins with a width of ten in M AC for 0.5 < p T < 5 GeV in the Pb+Pb analysis and 0.3 < p T < 5 GeV in the p+Pb analysis. It was crosschecked that the variables of interest obtained with a finer binning in M AC are consistent with the measurement with the nominal binning.
To enable comparisons with the theoretical predictions and with future experimental results, measurements obtained in M AC are presented as a function of the ATLAS ID multiplicity N ch of 0.5 < p T < 5 GeV and |η| < 2.5. They are projected from the M AC values taking into account tracking efficiency and fake-track production as described in the pre-vious section. A similar analysis procedure is described in Ref. [22]. For the N part dependencies in the Pb+Pb system, the results measured in M AC multiplicity intervals are averaged, with weights equal to the probabilities to find any given M AC value in the centrality intervals.
The formulation of the modified PCC ρ(v n {2} 2 , [ p T ]) requires that there should be at least two tracks in each region (A, B, and C). Further, Var v n {2} 2 dyn calculated according to Eq. (6) can be negative at low multiplicities due to statistical fluctuations, which renders Eq. (4) invalid because of the √ Var v n {2} 2 dyn term. For each M AC bin, p T interval, and harmonic, a criterion is applied that Var v n {2} 2 dyn needs to be positive at a level of at least one standard deviation of its statistical uncertainty. Results presented as a function of N ch are produced only for those M AC intervals. For the N part dependencies in the Pb+Pb system, it is additionally required for each centrality interval that the fraction of rejected events due to this criterion does not exceed 1%.

Systematic uncertainties
The systematic uncertainty is estimated by varying individual aspects of the analysis. The systematic uncertainties for the main p T interval are discussed for each collision system. Systematic uncertainties for the other p T intervals behave consistently with the ones for the main p T interval. Since the modified PCC ρ(v n {2} 2 , [ p T ]) is a ratio of quantities which are calculated using tracks, many variations largely cancel out and the resulting systematic uncertainties are small. To suppress the statistical fluctuations and to get more robust estimation of systematic uncertainties, they are averaged over several, wide ranges of the charged-particle multiplicity. For each uncertainty source and for each measurement point, the maximum variation from the baseline measurement is used. The total resulting uncertainty is the sum of the individual contributions combined in quadrature. The following sources of systematic uncertainties are considered.

Track selection
The tracking performance has a relatively small impact on v n {2}, but it directly affects the [ p T ] and c k via the admixture of the fake tracks, especially at low p T . To assess the impact on ρ(v n {2} 2 , [ p T ]), the measurement is repeated with tracks selected with looser and tighter track quality criteria, thus increasing and decreasing the fake-track rate, respectively. The weights used in the evaluation of measured quantities take the modified selection into account. The loose track selection in the Pb+Pb analysis relaxes requirements on the number of pixel and SCT hits to at least one and six, respectively. Additionally, the requirements on the transverse and longitudinal impact parameters of the track are relaxed to 1.5 mm. The tighter selection in the Pb+Pb analysis tightens the requirement on the transverse and lon-gitudinal impact parameters of the track to 0.5 mm. For the p+Pb analysis, the loose selection relaxes the requirements on the transverse and longitudinal impact parameters of the track to 2 mm and on the impact parameter significances to less than 4. In the tight selection, the impact parameter values and their significances must be less than 1 mm and 2, respectively. For each of the two track selections the absolute difference is calculated with respect to the baseline measurement: The largest difference is taken as a systematic uncertainty.
Detector material Since the tracks that are used in the calculation of ρ(v n {2} 2 , [ p T ]) are weighted by the inverse of the tracking efficiency, a bias in its estimation due to inaccurate modelling of the material in the detector may change the balance between low-and highp T tracks in the sums. Based on simulations, the estimated uncertainty in the detector description is obtained [43,44]. The resulting p T -and η-dependent uncertainties in the track efficiency of up to 4% are used to determine the systematic uncertainty.
Tracking azimuthal uniformity In this analysis, the weighting factors w correct for any non-uniformity in the azimuthal angle distribution of reconstructed tracks. The weights are obtained from the data by requiring azimuthal uniformity over the two-dimensional distribution of reconstructed tracks in the η-φ plane. The effect of that correction on the result is conservatively estimated by comparing the baseline measurement and the measurement obtained without applying this weight. The uncertainty is small, and it envelopes potential effects of imperfections in the weighting factors determination, including their dependence on the transverse momentum, collision centrality, run-by-run differences, on dead module maps or the vertex position.
Residual pile-up events The selection criteria discussed in Sect. 3 suppress the fraction of pile-up events accepted for analysis to almost zero in central Pb+Pb collisions. To estimate the systematic uncertainty related to pile-up, the measurement is conservatively repeated without this event rejection, resulting in at most a 1% difference in the most central Pb+Pb events for the ρ(v 2 {2} 2 , [ p T ]) coefficient. The p+Pb data were taken with higher pile-up than the Pb+Pb data. To estimate the impact of contamination by residual pile-up events, p+Pb results were obtained with only the vertex criteria applied. The variation covers the estimated residual pileup fraction in events of the highest track multiplicity [36].

Centrality selection
The minimum-bias trigger is fully efficient for the 0-85% centrality interval. However, the total fraction of inelastic Pb+Pb events selected is known only to 1% accuracy due to trigger inefficiency and possible sample contamination in more peripheral interactions. The centrality is estimated using the E FCal T distribution [6,10] and the Glauber model [38] to obtain the mapping from the observed E FCal T to the number of nucleons participating in the collision, N part . The modified PCC uncertainty is evaluated by repeating the analysis with the altered centrality selections on the E FCal T distribution, which results in ±1% uncertainty in the total fraction of inelastic Pb+Pb events. The centrality selection contributes mainly to uncertainties for peripheral collisions. Figure 1 shows the magnitude of the systematic uncertainties δρ(v n {2} 2 , [ p T ]) for n = 2 − 4 in Pb+Pb collisions as a function of N ch . In Pb+Pb collisions, the systematic uncertainty of the measured correlation coefficients across different order harmonics and centralities is not dominated by a single source. One of the largest uncertainties comes from restoring the azimuthal uniformity, and dominates for the second order harmonic in the most central collisions and for the third and fourth order harmonics almost over the full centrality range. A sizeable contribution to the uncertainty for all three harmonics is due to the track selection. The impact of the detector material is rather small except for a significant contribution for the forth order harmonic in the most central events. The residual pile-up in Pb+Pb collisions gives a negligible contribution. Figure 1d shows systematic uncertainties for ρ(v 2 {2} 2 , [ p T ]) coefficients in p+Pb collisions for the main interval of 0.3 < p T < 2 GeV as a function of event activity. In p+Pb interactions the largest uncertainty in the most active collisions (N ch > 150) originates from pile-up. The track selection is a source of sizeable uncertainty for this collision system, while the azimuthal uniformity correction procedure and the detector material have a small impact.
Details on the contributions to systematic uncertainties from different sources of c k , Var v n {2} 2 dyn and cov(v n {2} 2 , [ p T ]) are included in the Appendix.

The constituents of the modified PCC
The constituents of the modified PCC, c k , Var v n {2} 2 dyn and cov(v n {2} 2 , [ p T ]) and are combined, using Eq. (4), to obtain ρ. Figure 2 shows the dynamical p T fluctuation coefficient c k as a function of charged-particle multiplicity in Pb+Pb and p+Pb collision systems for tracks in three different p T intervals. A strong decrease of c k with increasing N ch is observed in all measured results. A similar decrease was seen for c k in Au+Au and Pb+Pb data at lower centreof-mass energies [28,29], evaluated for lower p T range, 0.15 < p T < 2 GeV, not accessible with the ATLAS detector. For the same N ch , the c k values differ by an order of magnitude for different p T ranges of tracks used in the analysis. For the intervals with the same lower p T limit, the c k values are higher for the interval with the larger upper p T limit. Figure 3 shows Var v n {2} 2 dyn for n = 2 − 4 as function of N ch for Pb+Pb collisions. For low multiplicities, Var v n {2} 2 dyn increases with increasing N ch , reaching a maximum at N ch of approximately 500 (1000) for n = 2 (n = 3), respectively. At higher N ch values the variances decrease with multiplicity. The dynamical variance for n = 4, measured for N ch 500, decreases with increas-  observed. The largest variances are observed for the p T intervals with an increased lower limit. This is expected as the v n {2} value increases strongly with p T below 3 GeV [10]. Additionally, the interval in which the upper limit on p T is set to 5 GeV integrates the region with the highest values of v n {2} (which occur around 3 GeV) and thus the values of the variance are expected to be larger than that for the main p T range.
In Fig. 4, the covariances cov(v n {2} 2 , [ p T ]) are shown for the 2nd-, 3rd-, and 4th-order harmonics in Pb+Pb collisions and for the second-order harmonics in p+Pb collisions. They are presented as a function of N ch for three p T intervals. Significant positive correlations between v n {2} and [ p T ] are observed in the Pb+Pb events. The measured covariances depend on the charged-particle multiplicity and the p T range of the charged particles. In Pb+Pb collisions, a strong dependence on the multiplicity is observed for n = 2 and 4. The is similar for all p T intervals. It starts at negative values for N ch < 200 and rapidly increases with multiplicity up to ∼ 1500 particles where the increase slows down and reaches the maximum at N ch ≈ 4500 of 0.24-0.3, depending on the p T interval. At even higher N ch , the ρ(v 2 {2} 2 , [ p T ]) value decreases rapidly. The significant correlation observed for mid-central events suggests a connection between anisotropic and radial [46] flows which might be attributed to stronger hydrodynamic response (larger pressure gradients) to the large initial-state eccentricities [47]. The modified PCC multiplicity dependence could reflect a balance between stronger radial flow observed in central collision and the larger initial eccentricity seen in peripheral interactions. The decrease observed in central collisions, for N ch 5000, might be related to the increased role of initial-state fluctuations in anisotropic flow [27]. However, a complete understanding of this effect would require a more precise modelling of heavy ion collisions. The correlation coefficients calculated with the upper p T limit of 2 GeV are 10-20% smaller than the values obtained with a p T limit of 5 GeV. The correlation coefficient ρ(v 3 {2} 2 , [ p T ]) is evaluated in Pb+Pb collisions for the same three p T ranges. The magnitudes measured for ρ(v 3 {2} 2 , [ p T ]) are significantly smaller than those measured for ρ(v 2 {2} 2 , [ p T ]) and similar to the magnitudes of ρ(v 4 {2} 2 , [ p T ]). All three curves increase with N ch in the range of 1000 < N ch < 5000. At low values of N ch , a flattening of the trend can be noticed. In the most central collisions, a breakdown of the rise is seen, similarly to the ρ(v 2 {2} 2 , [ p T ]). Above N ch ∼ 1500, the curves for the two intervals with the same maximum p T are consistent with each other and are below the curve for the interval which uses tracks with p T up to 5 GeV. The largest values of ρ(v 4 {2} 2 , [ p T ]) are observed at N ch ≈ 1000. For high N ch , ρ(v 4 {2} 2 , [ p T ]) decreases with N ch up to about N ch ≈ 4000 and rises slowly at higher values. The trends obtained for p T intervals with the same minimum value are consistent above N ch ∼ 1500 as is the case for ρ(v 3 {2} 2 , [ p T ]). The decrease for N ch < 4000 might be due to a contribution to v 4 from a non-linear term containing v 2 2 , decreasing with increasing centrality [13]. However, a theoretical modelling of the initial state and its subsequent evolution would be required to support this interpretation. Similarly to the ρ( correlations measured with the larger upper p T limit have larger magnitudes. The results for the larger upper p T limit show the sensitivity of the ρ(v n {2} 2 , [ p T ]) coefficients to the high p T part of the particle spectrum contaminated with non-flow correlations from jets. On the other hand, the correlations measured for the intervals with fixed upper p T limit (2 GeV) and varied lower p T limits are similar, demonstrating insensitivity of the modified PCC coefficients to a significant change of the event charged-particle multiplicity as expected [25]. The fourth-order correlations are weaker than those for the second-order flow harmonic and for N ch > 4000 are comparable to ρ(v 3 {2} 2 , [ p T ]). The results for all harmonics indicate a change in the trend in events with high N ch around 4500, which suggests a change in the nature of the correlations in those events [47].
In p+Pb collisions, ρ(v 2 {2} 2 , [ p T ]) exhibits much weaker N ch dependence than that in Pb+Pb collisions. For the main p T interval, the modified PCC assumes a negative value of approximately −0.1 and is almost constant within uncertain- ties. Values for different lower p T limits are similar, and the ρ(v 2 {2} 2 , [ p T ]) magnitudes for the larger upper p T limit are smaller. The magnitude (and sign) of the modified PCC in p+Pb collisions is expected to be related to the distribution of the energy deposition in the initial state, as predicted by the hydrodynamic model [25]. In hydrodynamics, in p+Pb collision, for small sources a higher initial pressure gradients and smaller eccentricities are expected to be generated. This mechanism could lead to the negative correlation of the final state observables, this is the mean transverse momentum and higher order flow harmonics. Thus, the negative value of the modified PCC for v 2 {2} in p+Pb and peripheral Pb+Pb that is measured should provide valuable constraints for models describing the collectivity in small systems.

(d) (c)
response to the initial geometry of the source, characterized by the positive value of the modified PCC.

Comparison to theoretical predictions
To compare the Pb+Pb results with a theoretical prediction in Ref. [25], the ρ(v n {2} 2 , [ p T ]) coefficients for 0.5 < p T < 2 GeV are obtained as a function of centrality intervals expressed by N part using the procedure described in Sect. 4. Figure 7 shows the N part dependence of ρ(v n {2} 2 , [ p T ]) for n = 2 − 4 in Pb+Pb collisions. It resembles the trends observed in Fig. 5, which show the modified PCC as a function of N ch , a measure of event activity. The theoretical predictions of the ρ(v n {2} 2 , [ p T ]) coefficient are based on a model in which the initial conditions were generated with nucleon positions by a MC Glauber model [48]. These initial conditions are then evolved using the pressure-driven 3+1D hydrodynamical simulations with viscous effects followed by the statistical particle emission to match multiplicities observed experimentally [37]. The modified Pearson correlation coefficient is then extracted from the final-state particles. The predictions for all harmonics are consistent with the data within the large model uncertainties except for the most central collisions where the predictions underestimate the measured ρ(v 2 {2} 2 , [ p T ]) and for the semi-peripheral collisions, for N part ∼ 130, where the predictions overestimate the ρ(v 2 {2} 2 , [ p T ]) and underestimate ρ(v 4 {2} 2 , [ p T ]).

Summary
The first measurement of the modified PCC ρ(v n {2} 2 , [ p T ]), which quantifies the correlation between the flow harmonics and the mean transverse momentum, is performed by ATLAS experiment at the LHC. The measurement uses 22 μb −1 of Pb+Pb data and 28 nb −1 of p+Pb data at the same centre-ofmass energy per nucleon pair of 5.02 TeV.
The correlation coefficient for several charged-particle p T ranges is measured as a function of the number of charged particles N ch and, in Pb+Pb collisions, the average number of nucleons participating in the collision, N part . For the 2 nd -, 3 rd -, and 4 th -order harmonics, the measured quantities exhibit a dependence on the choice of charged-particle p T range. Measurements with an upper limit of 5 GeV on p T indicate a stronger correlation than those with an upper limit  [20]. In p+Pb collisions, the value of The modified PCC is a valuable tool for studying the dynamics of heavy-ion collisions. It provides a reliable estimate of the magnitude of correlations calculated using finite multiplicities. In comparison with existing results, it allows quantitative comparisons between the experimental data and theoretical models. The precise measurements of this observable, presented in this paper, provide useful insights into the interplay of the azimuthal anisotropies (azimuthal flow) and the mean event p T (radial flow), providing input for a better understanding of QGP dynamics and for constraining the theoretical models. The obtained ρ(v n {2} 2 , [ p T ]) coefficients for 0.5 < p T < 2 GeV were compared with a theoretical prediction based on the pressure-driven 3+1D hydrodynamical simulations with viscous effects. The predictions for all harmonics are consistent with the data within the large model uncertainties. The only exception are the most central collisions, where the predictions underestimate the measured ρ(v 2 {2} 2 , [ p T ]) and the semi-peripheral collisions, where the predictions overestimate the ρ(v 2 {2} 2 , [ p T ]) and underestimate ρ(v 4 {2} 2 , [ p T ]). Sizeable positive correlations observed for non-peripheral Pb+Pb collisions support a qualitatively expected scenario in which the azimuthal flow originates from the pressure gradients.
In small system collisions the magnitude of the transverse flow is expected to be very sensitive to the size of the initial source in the hydrodynamic model. In particular, in the compact source scenario in p+Pb collisions, the smaller source sizes are expected to yield larger transverse flow and smaller initial eccentricities. The negative sign of the modified PCC measured in p+Pb collisions seems to support the compact source scenario, and indicates the role of the initial conditions in these systems.
Acknowledgements We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently. We acknowledge the support of ANPCyT

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: "All ATLAS scientific output is published in journals, and preliminary results are made available in Conference Notes. All are openly available, without restriction on use by external parties beyond copyright law and the standard conditions agreed by CERN. Data associated with journal publications are also made available: tables and data from plots (e.g. cross section values, likelihood profiles, selection efficiencies, cross section limits, ...) are stored in appropriate repositories such as HEPDATA (http:// hepdata.cedar.ac.uk/). ATLAS also strives to make additional material related to the paper available that allows a reinterpretation of the data in the context of new theoretical models. For example, an extended encapsulation of the analysis is often provided for measurements in the framework of RIVET (http://rivet.hepforge.org/)." This information is taken from the ATLAS Data Access Policy, which is a public document that can be downloaded from http://opendata.cern.ch/record/413.] Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .