• Nie Znaleziono Wyników

Measurement of long-range pseudorapidity correlations and azimuthal harmonics in $\sqrt{s_{NN}}=5.02$ TeV proton-lead collisions with the ATLAS detector

N/A
N/A
Protected

Academic year: 2022

Share "Measurement of long-range pseudorapidity correlations and azimuthal harmonics in $\sqrt{s_{NN}}=5.02$ TeV proton-lead collisions with the ATLAS detector"

Copied!
29
0
0

Pełen tekst

(1)

Measurement of long-range pseudorapidity correlations and azimuthal harmonics in

s

N N

= 5.02 TeV proton-lead collisions with the ATLAS detector

G. Aad et al. (ATLAS Collaboration)

(Received 5 September 2014; published 9 October 2014)

Measurements of two-particle correlation functions and the first five azimuthal harmonics, v1 to v5, are presented, using 28 nb−1of p + Pb collisions at a nucleon-nucleon center-of-mass energy of

sNN = 5.02 TeV measured with the ATLAS detector at the LHC. Significant long-range “ridgelike” correlations are observed for pairs with small relative azimuthal angle (|φ| < π/3) and back-to-back pairs (|φ| > 2π/3) over the transverse momentum range 0.4 < pT< 12 GeV and in different intervals of event activity. The event activity is defined by either the number of reconstructed tracks or the total transverse energy on the Pb-fragmentation side. The azimuthal structure of such long-range correlations is Fourier decomposed to obtain the harmonics vn

as a function of pT and event activity. The extracted vnvalues for n = 2 to 5 decrease with n. The v2 and v3

values are found to be positive in the measured pTrange. The v1is also measured as a function of pTand is observed to change sign around pT≈ 1.5–2.0 GeV and then increase to about 0.1 for pT> 4 GeV. The v2(pT), v3(pT), and v4(pT) are compared to the vncoefficients in Pb+ Pb collisions at√

sNN = 2.76 TeV with similar event multiplicities. Reasonable agreement is observed after accounting for the difference in the average pTof particles produced in the two collision systems.

DOI:10.1103/PhysRevC.90.044906 PACS number(s): 25.75.Dw

I. INTRODUCTION

One striking observation in high-energy nucleus-nucleus (A + A) collisions is the large anisotropy of particle pro- duction in the azimuthal angle φ [1,2]. This anisotropy is often studied via a two-particle correlation of particle pairs in relative pseudorapidity (η) and azimuthal angle (φ) [3,4].

The anisotropy manifests itself as a strong excess of pairs at

φ ∼ 0 and π, and the magnitude of the excess is relatively constant out to large |η| [5–9]. The azimuthal structure of this “ridgelike” correlation is commonly characterized by its Fourier harmonics, dNpairs/dφ ∼ 1 + 2

nv2ncos nφ.

While the elliptic flow, v2, and triangular flow, v3, are the dominant harmonics in A + A collisions, significant v1, v4, v5, and v6harmonics have also been measured [8–13]. These vnvalues are commonly interpreted as the collective hydrody- namic response of the created matter to the collision geometry and its fluctuations in the initial state [14]. The success of hydrodynamic models in describing the anisotropy of particle production in heavy-ion collisions at BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC) places important constraints on the transport properties of the produced matter [15–20].

For a small collision system, such as proton-proton (p + p) or proton-nucleus (p + A) collisions, it was assumed that the transverse size of the produced system is too small for the hydrodynamic flow description to be applicable. Thus, it came

Full author list given at the end of the article.

Published by the American Physical Society under the terms of the Creative Commons Attribution 3.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI.

as a surprise that ridge-like structures were also observed in two-particle correlations in high-multiplicity p + p [21]

and proton-lead (p + Pb) [22–24] collisions at the LHC and later in deuteron-gold collisions [25] at RHIC. A Fourier decomposition technique has been employed to study the azimuthal distribution of the ridge in p + Pb collisions. The transverse momentum (pT) dependence of the extracted v2

and v3 [23,24], and the particle-mass dependence of v2 [26]

are found to be similar to those measured in A + A collisions.

Large v2 coefficients are also measured via the four-particle cumulant method [27–29], suggesting that the ridge reflects genuine multiparticle correlations.

The interpretation of the long-range correlations in high- multiplicity p + p and p + Pb collisions is currently a subject of intense study. References [30–33] argue that the produced matter in these collisions is sufficiently voluminous and dense that the hydrodynamic model framework may still apply.

However, models based on gluon saturation and color con- nections suggest that the long-range correlations are an initial- state effect, intrinsic to QCD at high gluon density [34–38].

Recently, a hybrid model that takes into account both the initial- and the final-state effects has been proposed [39].

All these approaches can describe, qualitatively and even quantitatively, the v2and v3data in the p + Pb collisions.

To provide more insights into the nature of the ridge correlation and to discriminate between different theoretical interpretations, this paper provides a detailed measurement of the two-particle correlation and vn coefficients in p + Pb collisions at a nucleon-nucleon center-of-mass energy of

sNN = 5.02 TeV. The data correspond to an integrated luminosity of approximately 28 nb−1, recorded in 2013 by the ATLAS experiment at the LHC. This measurement benefits from a dedicated high-multiplicity trigger (see Sec. II B) that provides a large sample of high-multiplicity events, not only extending the previous v2 and v3 results to higher pT,

(2)

but also enabling the first measurement of v1, v4, and v5. The results are extracted independently for two different event-activity definitions: the total transverse energy in the forward calorimeter on the Pb-fragmentation side1 (−4.9 <

η < −3.2), ETPb, or the number of reconstructed tracks in

|η| < 2.5, Nchrec. The results are also compared to the Pb+ Pb data with similar multiplicity. The analysis technique follows closely the previous ATLAS study of v2 and v3 based on a much smaller dataset from a short p + Pb run in 2012 [24].

II. EXPERIMENTAL SETUP A. Detector and dataset

The ATLAS detector [40] provides nearly full solid- angle coverage of the collision point with tracking detectors, calorimeters, and muon chambers. The measurements pre- sented in this paper are performed using the ATLAS inner detector (ID), forward calorimeters (FCals), minimum-bias trigger scintillators (MBTSs), zero-degree calorimeter (ZDC), and the trigger and data acquisition systems. The ID measures charged particles within the pseudorapidity region|η| < 2.5 using a combination of silicon pixel detector, silicon microstrip detector (SCT), and a straw-tube transition radiation tracker, all immersed in a 2-T axial magnetic field. The MBTS detects charged particles over 2.1 < |η| < 3.9 using two hodoscopes of 16 counters positioned at z = ±3.6 m. The FCal consists of two sections that cover 3.2 < |η| < 4.9. The FCal modules are composed of tungsten and copper absorbers with liquid argon as the active medium, which provide ten interaction lengths of material. The ZDC, situated at≈140 m from the collision vertex, detects neutral particles, mostly neutrons and photons, with|η| > 8.3.

This analysis uses approximately 28 nb−1of p + Pb data recorded by the ATLAS experiment at the LHC in 2013.

The LHC was configured with a 4-TeV proton beam and a 1.57-TeV-per-nucleon Pb beam that together produced collisions at √

sNN = 5.02 TeV. The beam directions were reversed approximately half-way through the running period.

The higher energy of the proton beam results in a net rapidity shift of the nucleon-nucleon center-of-mass frame relative to the ATLAS rest frame. This rapidity shift is 0.47 towards the proton beam direction.

B. Trigger

The minimum-bias (MB) level-1 (L1) trigger [41] used for this analysis requires a signal in at least one MBTS counter on

1ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the center of the detector and the z axis along the beam pipe. The x axis points from the IP towards the center of the LHC ring, and the y axis completes the right-handed system. Cylindrical coordinates (r,φ) are used in the transverse plane, φ being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2). During 2013 p + Pb data taking, the beam directions were reversed approximately half-way through the running period, but in presenting results the direction of the proton beam is always chosen to point to positive η.

TABLE I. A list of thresholds in ETL1and NtrkHLTfor the HMTs used in this analysis.

NtrkHLT 100 130 150 180 200 225

EL1T (GeV) 10 10 50 50 65 65

each side, or a signal in the ZDC on the Pb-fragmentation side with the trigger threshold set just below the peak corresponding to a single neutron. A timing requirement based on signals from each side of the MBTS is imposed to suppress beam backgrounds. Owing to the high event rate during the run, only a small fraction of the MB events (∼1/1000) were recorded.

To enhance the number of events with high multiplicity, a dedicated high-multiplicity trigger (HMT) was implemented, which uses the ATLAS L1 and high-level trigger (HLT) systems [42]. At L1, the total transverse energy EL1T in the FCal rapidity interval is required to be above a certain threshold.

In the HLT, the charged-particle tracks are reconstructed by requiring at least two hits in the pixel detector and three hits in the SCT. For each event, the collision vertex reconstructed with the highest number of online tracks is selected, and the number of tracks (NtrkHLT) associated with this vertex with pT>

0.4 GeV and a distance of closest approach of less than 4 mm is calculated.

The HMT triggers are implemented by requiring different thresholds on the values of ETL1 and NtrkHLTwith prescale factors adjusted to the instantaneous luminosity provided by the LHC [42]. This analysis uses the six pairs of thresholds on EL1T and NtrkHLTlisted in TableI. The NtrkHLT 225 trigger was not prescaled throughout the entire running period.

III. DATA ANALYSIS A. Event and track selections

In the off-line analysis, p + Pb events are required to have a reconstructed vertex containing at least two associated off- line tracks, with its z position satisfying |zvtx| < 150 mm.

Noncollision backgrounds and photonuclear interactions are suppressed by requiring at least one hit in a MBTS counter on each side of the IP and the difference between times measured on the two sides to be less than 10 ns. In the 2013 p + Pb run, the luminosity conditions provided by the LHC result in an average probability of 3% that an event contains two or more p + Pb collisions (pileup). The pileup events are suppressed by rejecting events containing more than one good reconstructed vertex. The remaining pileup events are further suppressed based on the signal in the ZDC on the Pb-fragmentation side.

This signal is calibrated to the number of detected neutrons (Nn) based on the location of the peak corresponding to a single neutron. The distribution of Nnin events with pileup is broader than that for the events without pileup. Hence, a simple cut on the high tail end of the ZDC signal distribution further suppresses the pileup, while retaining more than 98% of the events without pileup. After this pileup rejection procedure, the residual pileup fraction is estimated to be10−2 in the event class with the highest track multiplicity studied in this analysis. About 57× 106 MB-selected events and 15× 106 HMT-selected events are included in this analysis.

(3)

Charged-particle tracks are reconstructed in the ID using an algorithm optimized for p + p MB measurements [43]: The tracks are required to have pT > 0.3 GeV and |η| < 2.5, at least seven hits in the pixel detector and the SCT, and a hit in the first pixel layer when one is expected. In addition, the transverse (d0) and longitudinal (z0 sin θ) impact parameters of the track relative to the vertex are required to be less than 1.5 mm. They are also required to satisfy|d0d0| < 3 and |z0sin θ/σz| < 3, respectively, where σd0 and σz are uncertainties on d0 and z0sin θ obtained from the track-fit covariance matrix.

The efficiency, (pT,η), for track reconstruction and track selection cuts is obtained using p + Pb Monte Carlo events produced with version 1.38b of theHIJINGevent generator [44]

with a center-of-mass boost matching the beam conditions. The response of the detector is simulated usingGEANT4 [45,46] and the resulting events are reconstructed with the same algorithms as applied to the data. The efficiency increases with pTby 6%

between 0.3 and 0.5 GeV and varies only weakly for pT >

0.5 GeV, where it ranges from 82% at η = 0 to 70% at |η| = 2 and 60% at|η| > 2.4. The efficiency is also found to vary by less than 2% over the multiplicity range used in the analysis.

The extracted efficiency function (pT,η) is used in the corre- lation analysis, as well as to estimate the average efficiency- corrected charged-particle multiplicity in the collisions.

B. Characterization of the event activity

The two-particle correlation (2PC) analyses are performed in event classes with different overall activity. The event activity is characterized by either ETPb, the sum of transverse energy measured on the Pb-fragmentation side of the FCal with

−4.9 < η < −3.2, or Nchrec, the off-line-reconstructed track multiplicity in the ID with|η| < 2.5 and pT > 0.4 GeV. These event-activity definitions have been used in previous p + Pb analyses [21,22,24,27,28]. Events with larger activity have on average a larger number of participating nucleons in the Pb nucleus and a smaller impact parameter. Hence, the term

“centrality,” familiar in A + A collisions, is used to refer to the event activity. The terms “central” and “peripheral” are used to refer to events with large activity and small activity, respectively.

Owing to the wide range of trigger thresholds and the prescale values required by the HMT triggers, the ETPb and Nchrec distributions are very different for the HMT events and the MB events. To properly include the HMT events in the event-activity classification, an event-by-event weight, w = 1/P , is utilized. The combined probability, P , for a given event to be accepted by the MB trigger or any of the HMT triggers is calculated via the inclusion-exclusion principle as

P = 

1iN

pi− 

1i<jN

pipj+ 

1i<j<kN

pipjpk− · · · , (1) where N is the total number of triggers and piis the probability for the ith trigger to accept the event, defined as zero if the event does not fire the trigger and otherwise as the inverse of the prescale factor of the trigger. The higher-order terms in Eq. (1)

account for the probability of more than one trigger being fired.

The weight factor, w, is calculated and applied event by event.

The distribution for all events after reweighting has the same shape as the distribution for MB events, as should be the case if the reweighting is done correctly.

Figure1shows the distribution of Nchrec(left panels) and ETPb (right panels) for the MB and MB+ HMT events before (top panels) and after (bottom panels) the reweighting procedure.

For MB-selected events, the reweighted distribution differs from the original distribution by a constant factor, reflecting the average prescale. The multiple steps in the Nchrecdistribution (top-left panel) reflect the rapid turn-on behavior of individual HMT triggers in Nchrec. The broad shoulder in the ETPbdistribu- tion (top-right panel) is attributable to the finite width of the Nchrecvs ETPbcorrelation, which smears the contributions from different HMT triggers in ETPb. All these structures disappear after the reweighting procedure. The results of this analysis are obtained using the MB+ HMT combined dataset with event reweighting.

Owing to the relatively slow turn on of the HMT triggers as a function of ETPb [Fig. 1(b)], the events selected in a given ETPb range typically receive contributions from several HMT triggers with very different weights. Hence, the effective increase in the number of events from the HMT triggers in the large ETPbregion is much smaller than the increase in the large Nchrecregion.

Figure 2(a) shows the correlation between ETPb and Nchrec from MB+ HMT p + Pb events. This distribution is similar to that obtained for the MB events, except that the HMT triggers greatly extend the reach in both quantities. The ETPbvalue grows with increasing Nchrec, suggesting that, on average, EPbT in the nucleus direction correlates well with the particle production at midrapidity. However, the broad distribution of EPbT at fixed Nchrecalso implies significant fluctuations. To study the relation between ETPband Nchrec, events are divided into narrow bins in Nchrec, and the mean and root-mean-square values of the ETPb

distribution are calculated for each bin. The results are shown in Fig.2(b). A nearly linear relation betweenETPb and Nchrecis observed. This relationship is used to match a given Nchrecevent class to the corresponding ETPbevent class. This approximately linear relation can also be parameterized [indicated by the solid line in Fig.2(b)] as

ETPb

/GeV ≈ 0.60Nchrec. (2)

The 2PC analysis is performed in different intervals of the event activity defined by either ETPb or Nchrec. TableIIgives a list of representative event-activity classes, together with the fraction of MB+ HMT events [after reweighting as shown in Fig.2(a)] contained in each event class. The table also provides the average Nchrecand ETPbvalues for each event-activity class, as well as the efficiency-corrected number of charged particles within|η| < 2.5 and pT> 0.4 GeV, Nch. The event classes defined in narrow ETPb or Nchrec ranges are used for detailed studies of the centrality dependence of the 2PC, while the event classes in broad ETPbor Nchrecranges are optimized for the studies of the pTdependence. As the number of events at large ETPbis smaller than at large Nchrec, the main results in this paper are presented for event classes defined in Nchrec.

(4)

Events

1 10 102

103

104

105

106

w/o weighting MB MB+HMT

(a)

Events

1 10 102

103

104

105

106

w/o weighting

(b)

rec

Nch

0 100 200 300

Events

1 102

104

106

108

109 re-weighted

(c)

ATLAS p+Pb

= 5.02 TeV sNN

28 nb-1 int >65 GeV L

L1

>225, ET HLT

Ntrk

>65 GeV

L1

>200, ET HLT

Ntrk

>50 GeV

L1

>180, ET HLT

Ntrk

>50 GeV

L1

>150, ET HLT

Ntrk

>10 GeV

L1

>130, ET HLT

Ntrk

>10 GeV

L1

>100, ET HLT

Ntrk

[GeV]

Pb

ET

0 100 200

Events

1 102

104

106

108

109 re-weighted

(d)

ATLAS p+Pb

= 5.02 TeV sNN

28 nb-1 int L

FIG. 1. (Color online) The distributions of Nchrec(left panels) and ETPb(right panels) for MB and MB+ HMT events before (top panels) and after (bottom panels) applying an event-by-event weight (see text). The smaller symbols in the top panels represent the distributions from the six HMT triggers listed in TableI.

C. Two-particle correlation

For a given event class, the 2PCs are measured as functions of relative azimuthal angle, φ = φa− φb, and relative pseudorapidity, η = ηa− ηb, with |η|  ηmax= 5. The labels a and b denote the two particles in the pair, which may be selected from different pTintervals. The particles a and b are conventionally referred to as the “trigger” and “associated”

particles, respectively. The correlation strength, expressed in terms of the number of pairs per trigger particle, is defined as [4–6]

Y (φ,η) =

 B(φ,η)dφdη πηmax

S(φ,η) B(φ,η)

 , Y (φ) =

 B(φ)dφ π

S(φ) B(φ)



, (3)

where S and B represent pair distributions constructed from the same event and from “mixed events” [4], respectively,

which are then normalized by the number of trigger particles in the event. These distributions are also referred to as per-trigger yield distributions. The mixed-event distribution, B(φ,η), measures the distribution of uncorrelated pairs.

The B(φ,η) distribution is constructed by choosing the two particles in the pair from different events of similar Nchrec(match to|Nchrec| < 10 tracks), ETPb(match to|ETPb| < 10 GeV), and zvtx(match to|zvtx| < 10 mm), so that B(φ,η) properly reflects the known detector effects in S(φ,η). The one- dimensional (1D) distributions S(φ) and B(φ) are obtained by integrating S(φ,η) and B(φ,η), respectively, over a η range. The region |η| < 1 is chosen to focus on the short-range features of the correlation functions, while the region|η| > 2 is chosen to focus on the long-range features of the correlation functions. These two regions are hence referred to as the “short-range region” and the “long-range region,” respectively. The normalization factors in front of the S/B ratio are chosen such that the (φ,η)-averaged

(5)

1 102

104

106

108

109

rec

Nch

0 100 200 300 400

[GeV]Pb TE

0 100 200

300 ATLAS p+Pb = 5.02 TeV sNN

28 nb-1

Lint

MB+HMT, re-weighted (a)

rec

Nch

0 100 200 300

[GeV]Pb TEσ or 〉Pb TE〈

0 50 100 150 200

ATLAS p+Pb = 5.02 TeV sNN

28 nb-1 int L

Pb ET

Pb ET

σ

0.6 Linear fit, slope (b)

FIG. 2. (Color online) (a) Correlation between EPbT and Nchrecin MB+ HMT events. (b) The mean ETPb and root-mean-square σEPbT of the ETPbdistributions for slices of Nchrec. The line is a linear fit to all points.

value of B(φ,η) and φ-averaged value of B(φ) are both unity. When measuring S and B, pairs are filled in one quadrant of the (φ,η) space and then reflected to the other quadrants [24]. To correct S(φ,η) and B(φ,η) for the individual inefficiencies of particles a and b, the pairs are

weighted by the inverse product of their tracking efficiencies 1/(ab). Remaining detector distortions not accounted for in the efficiency largely cancel in the S/B ratio.

Examples of two-dimensional (2D) correlation functions are shown in Fig. 3for charged particles with 1 < pa,bT < 3

TABLE II. A list of event-activity classes defined in Nchrec (left) and ETPb(right) ranges, where the notation [a,b) implies a  Nchrecor ETPb< b. For each event class, the fraction of MB + HMT events after trigger reweighting [Fig.2(a)], the average values ofETPb and Nchrec and the efficiency-corrected average number of charged particles with pT> 0.4 GeV and |η| < 2.5, Nch, are also listed.

Event-activity classes based on Nchrec Event-activity classes based on ETPb

Nchrecrange Fraction ETPb Nchrec Nch ETPbrange Fraction EPbT Nchrec Nch

(GeV) (GeV) (GeV)

<20 0.31 7.3 10.3 12.6 ± 0.6 <10 0.28 4.8 12.4 15.4 ± 0.7

[20,40) 0.27 18.6 29.1 37.9 ± 1.7 [10,23) 0.26 16.1 29.2 38.1 ± 1.7

[40,60) 0.19 30.8 48.8 64.3 ± 2.9 [23,37) 0.19 29.5 47.3 62.3 ± 2.8

[60,80) 0.12 42.8 68.6 90.7 ± 4.1 [37,52) 0.12 43.8 64.0 84.7 ± 3.8

[80,100) 0.064 54.9 88.3 117± 5 [52,68) 0.067 58.8 80.4 107± 5

[100,120) 0.029 66.4 108 144± 7 [68,83) 0.028 74.2 96.1 128± 6

[120,140) 0.011 78.4 128 170± 8 [83,99) 0.012 89.7 111 147± 7

[140,160) 0.0040 90.3 148 196± 9 [99,116) 0.0043 106 126 168± 8

[160,180) 0.0013 102 168 223± 10 [116,132) 0.0012 122 141 187± 8

[180,200) 3.6 × 10−4 113 187 249± 11 [132,148) 3.6 × 10−4 138 155 206± 9

[200,220) 9.4 × 10−5 125 207 276± 12 [148,165) 1.0 × 10−4 155 169 225± 10

[220,240) 2.1 × 10−5 134 227 303± 14 [165,182) 2.2 × 10−5 171 184 244± 11

[240,260) 4.6 × 10−6 145 247 329± 15 [182,198) 4.6 × 10−6 188 196 261± 12

[260,290) 1.1 × 10−6 157 269 358± 16 [198,223) 1.1 × 10−6 206 211 281± 13

[290,370) 8.9 × 10−8 174 301 393± 18 [223,300) 9.6 × 10−8 232 230 306± 14

<40 0.58 12.5 19.0 24.4 ± 1.1 <25 0.59 10.2 21.7 28.0 ± 1.3

[40,80) 0.32 35.3 56.4 74.4 ± 3.3 [25,50) 0.27 35.1 54.7 72.2 ± 3.3

[80,110) 0.081 56.8 91.7 122± 6 [50,75) 0.096 61.5 81.4 108± 5

[110,140) 0.023 74.2 121 161± 7 [75,100) 0.025 84.5 106 141± 6

[140,180) 0.0053 93.0 153 203± 9 [100,130) 0.0051 110 130 173± 8

[180,220) 4.6 × 10−4 116 192 255± 12 [130,165) 5.6 × 10−4 141 156 208± 9

[220,260) 2.6 × 10−5 136 231 307± 14 [165,200) 2.7 × 10−5 174 186 248± 11

[260,370) 1.2 × 10−6 158 271 361± 16 [200,300) 1.0 × 10−6 208 214 284± 13

(6)

Δφ 0 2 4

Δη

)ηΔ, φΔY(

0.12 0.14

-4 -2 0 2 4 < 10 GeV

Pb

ET

(a)

Δφ 0 2 4

Δη

)ηΔ, φΔY( 0.1

0.12 0.14

-4 -2 0 2 4 < 20

rec

Nch

(b)

Δφ 0 2 4

Δη

)ηΔ, φΔY( 0.95

1

-4 -2 0 2 4 100 GeV

Pb≥ ET

(c)

Δφ 0 2 4

Δη

)ηΔ, φΔY( 1.65

1.7 1.75

-4 -2 0 2 4

≥ 220

rec

Nch

(d)

ATLAS

p+Pb

28 nb-1 int

= 5.02 TeV, L sNN

< 3 GeV

a,b

1 < pT

FIG. 3. (Color online) The 2D correlation function in φ and η for the peripheral event class selected by either (a) ETPb< 10 GeV or (b) Nchrec< 20 and the central event class selected by either (c) EPbT  100 GeV or (d) Nchrec 220.

GeV in low-activity events, EPbT < 10 GeV or Nchrec< 20, in the top panels, and high-activity events, ETPb> 100 GeV or Nchrec>

220, in the bottom panels. The correlation for low-activity events shows a sharp peak centered at (φ,η) = (0,0) owing to short-range correlations for pairs resulting from jets, high-pT resonance decays, and Bose-Einstein correlations.

The correlation function also shows a broad structure at

φ ∼ π from low-pT resonances, dijets, and momentum conservation that is collectively referred to as “recoil” [24]

in the remainder of this paper. In the high-activity events, the correlation reveals a flat ridgelike structure at φ ∼ 0 (the “near side”) that extends over the full measured η range. This η independence is quantified by integrating the 2D correlation functions over |φ| < 1 to obtain Y (η) =



|φ|<1Y (φ,η)φ. The yield associated with the near-side short-range correlation peak centered at (φ,η) = (0,0) can then be estimated as

Yn-peak =

|η|<1Y (η)dη − 1 5− ηmin

×

ηmin <|η|<5

Y (η)dη, (4)

where the second term accounts for the contribution of uncorrelated pairs and the ridge component under the near-side peak. The default value of Yn-peakis obtained with a lower-end of the integration range of ηmin = 2, but the value of ηmin is varied from 2 to 4 to check the stability of Yn-peak. The

distribution at φ ∼ π (the “away side”) is also broadened in high-activity events, consistent with the presence of a long- range component in addition to the recoil component [24].

This recoil component can be estimated from the low-activity events and subtracted from the high-activity events using the procedure detailed in the next section.

D. Recoil subtraction

The correlated structure above a flat pedestal in the correlation functions is calculated using a zero-yield-at- minimum (ZYAM) method [4,47] following previous mea- surements [22–24],

Ycorr(φ,η) =

B(φ,η)dφdη πηmax

×

S(φ,η) B(φ,η)− bZYAM

 ,

(5) Ycorr(φ) =

B(φ)dφ π

×

S(φ) B(φ)− bZYAM

 ,

where the parameter bZYAM represents the pedestal formed by uncorrelated pairs. A second-order polynomial fit to the 1D Y (φ) distribution in the long-range region is used to find the location of the minimum point, φZYAM, and from this the value

(7)

of bZYAM is determined and subtracted from the 2D correlation function. The Ycorr(φ,η) functions differ, therefore, by a constant from the Y (φ,η) functions, such as those in Fig.3.

In low-activity events, Ycorr(φ,η) contains mainly the short-range correlation component and the recoil component.

In high-activity events, the contribution from the long-range

“ridge” correlation also becomes important. This long-range component of the correlation function in a given event class is obtained by estimating the short-range correlation component using the peripheral events and is then subtracted,

Ysub(φ,η) = Y (φ,η) − αYpericorr(φ,η), (6) Ysub(φ) = Y (φ) − αYpericorr(φ),

where the Ycorr in a low-activity or peripheral event class, denoted by Ypericorr, is used to estimate and subtract [hence, the superscript “sub” in Eq. (6)] the short-range correlation at the near side and the recoil at the away side. The parameter α is chosen to adjust the near-side short-range correlation yield in the peripheral events to match that in the given event class for each pTa and pbT combination, α = Yn-peak/Yperin-peak. This scaling procedure is necessary to account for enhanced short-range correlations and away-side recoil in higher-activity events, under the assumption that the relative contribution of

the near-side short-range correlation and away-side recoil is independent of the event activity. A similar rescaling procedure has also been used by the CMS Collaboration [28]. The default peripheral event class is chosen to be ETPb< ET0= 10 GeV.

However, the results have also been checked with other E0T values, as well as with a peripheral event class defined by Nchrec< 20. In the events with the highest multiplicity, the value of α determined with the default peripheral event class varies from ∼2 at pT ≈ 0.5 GeV to ∼1 for pT> 3 GeV, with a pT-dependent uncertainty of 3%–5%.

The uncertainty on bZYAMonly affects the recoil-subtracted correlation functions through the Ypericorr term in Eq. (6). This uncertainty is usually very small in high-activity p + Pb collisions, owing to their much larger pedestal level than for the peripheral event class.

Figures4(a)and4(b)show, respectively, the 2-D correlation functions before and after the subtraction procedure given by Eq. (6). Most of the short-range peak and away-side recoil structures are removed by the subtraction, and the remaining distributions exhibit a φ-symmetric double ridge that is al- most independent of η. Figure4(c)shows the corresponding 1D correlation functions before and after recoil subtraction in the long-range region of|η| > 2. The distribution at the near-side is not affected because the near-side short-range peak

Δφ 0 2 4

Δη

)ηΔ, φΔY( 1.65

1.7 1.75

-4 -2 0 2 4

≥ 220

rec

Nch (a)

Δφ 0 2 4

Δη

)ηΔ, φΔ(sub Y 1.65

1.7 1.75

-4 -2 0 2 4

≥ 220

rec

Nch (b)

φ Δ

0 1 2 3

Per-trigger yield

11 11.2 11.4 11.6

φ) Δ Y(

φ) Δ

sub( Y ATLAS p+Pb

= 5.02 TeV sNN

28 nb-1 int≈ L

| < 5 η Δ 2 < |

<3 GeV

a,b

1< pT

(c)

η| Δ

|

0 2 4

n

v

0 0.05 0.1 0.15

unsub

v2

× 2 unsub

v3

× 3 unsub

v4

v2

× 2

v3

× 3

v4

(d)

FIG. 4. (Color online) The 2D correlation function in φ and η for events with Nchrec 220 (a) before and (b) after subtraction of the peripheral yield. Panel (c) shows the corresponding 1D correlation functions in φ for pairs integrated over 2 < |η| < 5 from panels (a) and (b), together with Fourier fits including the first five harmonics. Panel (d) shows the second-, third-, and fourth-order Fourier coefficients as functions of|η| calculated from the 2D distributions in panel (a) or panel (b), represented by the open or solid symbols, respectively. The error bars and shaded boxes are statistical and systematic uncertainties, respectively.

(8)

is narrow in η [Fig.4(a)], while the away-side distribution is reduced owing to the removal of the recoil component.

E. Extraction of the azimuthal harmonics associated with long-range correlation

The azimuthal structure of the long-range correlation is studied via a Fourier decomposition similar to the approach used in the analysis of Pb+ Pb collisions [7,9],

Ysub(φ) =

 Ysub(φ)dφ π

1+

n

2vn,ncos(nφ)

, (7) where vn,nare the Fourier coefficients calculated via a discrete Fourier transformation,

vn,n=

N

m=1cos(nφm)Ysub(φm)

N

m=1Ysub(φm) , (8)

where N = 24 is the number of φ bins from 0 to π. The first five Fourier coefficients are calculated as functions of pTaand pbTfor each event-activity class.

The azimuthal anisotropy coefficients for single particles, vn, can be obtained via the factorization relation commonly used for heavy-ion collisions [7,9,48],

vn,n

pTa,pbT

= vn

paT

vn

pbT

. (9)

From this the pTdependence of vnfor n = 2–5 are calculated as

vn

paT

= vn,n

paT,pTb

vn,n

pTb,pbT

, (10)

where the default transverse momentum range for the asso- ciated particle (b) is chosen to be 1 < pTb< 3 GeV, and the Fourier coefficient as a function of the transverse momentum of the trigger particle is denoted by vn(pTa) or simply vn(pT) where appropriate. The extraction of v1requires a slight modification and is discussed separately in Sec. IV C. The factorization behavior is checked by comparing the vn(pTa) obtained for different pTbranges, as discussed in Sec.IV B.

A similar Fourier decomposition procedure is also carried out for correlation functions without peripheral subtraction, i.e., Y (φ). The harmonics obtained in this way are denoted by vn,nunsuband vnunsub, respectively.

Figure 4(d) shows the azimuthal harmonics obtained by Fourier decomposition of the Y (φ,η) and Ysub(φ,η) distributions in Figs.4(a)and4(b)for different, narrow slices of η. The resulting vnunsub and vn values are plotted as functions of η for n = 2, 3, and 4. The vnvalues are much smaller than vunsubn for|η| < 1, reflecting the removal of the short-range correlations at the near side. The v2values are also systematically smaller than v2unsubfor|η| > 1, reflecting the removal of the away-side recoil contribution.

F. Systematic uncertainties

The systematic uncertainties in this analysis arise from pair acceptance, the ZYAM procedure, tracking efficiency, Monte Carlo consistency, residual pileup, and the recoil subtraction.

Each source is discussed separately below.

The correlation functions rely on the pair acceptance functions, B(φ,η) and B(φ) in Eq. (3), to reproduce detector acceptance effects in the signal distribution. A natural way of quantifying the influence of detector effects on vn,n

and vn is to express the single-particle and pair acceptance functions as Fourier series, similar to Eq. (7). The resulting coefficients for pair acceptance vn,ndetare the product of those for the two single-particle acceptances vndet,aand vndet,b. In general, the pair acceptance function in φ is quite flat: The maximum fractional variation from its average value is observed to be less than 0.001 for pairs integrated in 2 < |η| < 5, and the corresponding |vn,ndet| values are found to be less than 2× 10−4. These vn,ndet values are expected to mostly cancel in the correlation function, and only a small fraction contributes to the uncertainties in the pair acceptance function. Possible residual effects on the pair acceptance are evaluated following Ref. [9] by varying the criteria for matching in Nchrec, ETPb, and zvtx. In each case, the residual vn,ndet values are evaluated by a Fourier expansion of the ratio of the pair acceptances before and after the variation. This uncertainty varies in the range of (5–8)× 10−6. It is negligible for v2and v3, but becomes sizable for higher-order harmonics, particularly at low pT, where the vnvalues are small.

As discussed in Sec.III D, the value of bZYAM is determined by a second-order polynomial fit of the Y (φ) distribution.

The stability of the fit is studied by varying the φ range in the fit. The uncertainty in bZYAM depends on the local curvature around φZYAMand is estimated to be 0.0003–0.001 of the minimum value of Y (φ). This uncertainty contributes directly to Ycorr(φ), but contributes to Ysub(φ) and vn

indirectly through the peripheral subtraction [see Eq. (6)]. The resulting uncertainty on vn is found to be less than 2%, for all n.

The values of per-trigger yields, Y (φ), Ycorr(φ), and Ysub(φ), are sensitive to the uncertainty on the tracking efficiency correction for the associated particles. This un- certainty is estimated by varying the track quality cuts and the detector material in the simulation, re-analyzing the data using corresponding Monte Carlo efficiencies, and evaluating the change in the extracted yields. The resulting uncertainty is estimated to be 2.5% owing to the track selection and 2%–3% related to our limited knowledge of the detector material. The vn,n and vn values depend only on the shape of the Ysub(φ) distribution and hence are not sensitive to the tracking efficiency.

The analysis procedure is also validated by measuring vn

values in fully simulatedHIJINGevents [45,46] and comparing them to those measured using the generated particles. A small but systematic difference between the two results are included in the systematic uncertainties.

Nearly all of the events containing pileup are removed by the procedure described in Sec. III A. The influence of the residual pileup is evaluated by relaxing the pileup rejection criteria and then calculating the change in the per-trigger yields and vnvalues. The differences are taken as an estimate of the uncertainty and are found to be negligible in low event-activity classes and increase to 2% for events with ETPb> 200 GeV or Nchrec> 300.

(9)

TABLE III. Summary of relative systematic uncertainties for Y (φ), Ycorr(φ), and Ysub(φ).

Residual pair acceptance (%) 1–2

ZYAM procedure (%) 0.2–1.5

Tracking efficiency and material (%) 4.2

Residual pileup (%) 0–2

According to Table II, the low-activity events used in the peripheral subtraction (EPbT < E0T= 10 GeV) correspond to 28% of the MB-triggered events. The pair distributions for these events may contain a small genuine long-range component, leading to a reduction of the long-range correlation signal in a high-activity class via the peripheral subtraction procedure. The influence of this oversubtraction is evaluated by varying the definition of the low-activity events in the range of ET0 = 5 GeV to E0T= 20 GeV. The Ysub(φ) and vnvalues are calculated for each variation. The vn values are found to decrease approximately linearly with increasing E0T. The amount of oversubtraction can be estimated by extrapolating ET0 to zero. The estimated changes of vn and Ysub(φ) vary from less than 1% for ETPb> 100 GeV or Nchrec> 150 and increase for lower event-activity classes approximately as 1.5/Nchrec. The relative change of vn is also found to be independent of pT. As a cross check, the analysis is repeated by defining peripheral events as Nchrec< 20. The variation of vn

values is found to be consistent with the variation from varying ET0.

The stability of the scale factor, α, is evaluated by varying the η window of the long-range region in Eq. (4). A 3%–

5% uncertainty is quoted for α from these variations. The resulting uncertainty on vnfor n = 2–5 is within 1% at low pT

(<4 GeV) and increases to ∼10% at the highest pT. However, the v1extraction is directly affected by the subtraction of the recoil component, and hence the v1value is very sensitive to the uncertainty in α. The estimated uncertainty is 8%–12% for pT< 1 GeV and is about 20%–30% for pT > 3 GeV.

The different sources of the systematic uncertainties de- scribed above are added in quadrature to give the total systematic uncertainties for per-trigger yields and vn, which are summarized in Tables III and IV, respectively. The systematic uncertainty quoted for each source usually covers the maxmium uncertainty over the measured pT range and event-activity range. However, because v1(pT) changes sign

within 1.5–2.0 GeV (see Fig. 15), the relative uncertainties are quoted for pT< 1 GeV and pT> 3 GeV. The uncertainty of pair acceptance, which is less than 8× 10−6 for vn,n, was converted to percent uncertainties. This uncertainty can be significant at high pT.

IV. RESULTS

A. Correlation functions and integrated yields Figure5shows the 1D correlation functions after the ZYAM procedure, Ycorr(φ), in various ranges of pTa for a fixed pbT

range of 1–3 GeV. The correlation functions are obtained in the long-range region (|η| > 2) and are shown for events selected by Nchrec 220. This event class contains a small fraction (3× 10−5) of the MB p + Pb events with highest multiplicity. The correlation functions are compared to the distributions of the recoil component, αYpericorr(φ) in Eq. (6), estimated from the peripheral event class defined by EPbT <

10 GeV. The scale factor α is chosen such that the near-side short-range yield matches between the two event classes [see Eq. (6) and discussion around it]. Figure 5 shows a clear near-side excess in the full paTrange studied in this analysis. An excess above the estimated recoil contribution is also observed on the away side over the same pTrange.

To further quantify the properties of the long-range components, the Ycorr(φ) distributions are integrated over

|φ| < π/3 and |φ| > 2π/3, similar to the procedure used in previous analyses [23,24]. The integrated yields, Yint, are obtained in several event classes and are plotted as functions of paTin Fig.6. The near-side yields increase with trigger pT, reach a maximum at pT ∼ 3 GeV, and then decrease to a value close to zero at pT> 10 GeV. This trend is characteristic of the pTdependence of the Fourier harmonics in A + A collisions.

In contrast, the away-side yields show a continuous increase across the full pTrange, owing to the contribution of the recoil component that mostly results from dijets.

Figure7shows the centrality dependence of the long-range integrated yields for the event-activity based on Nchrec (left) and ETPb(right) for particles in 1 < pa,bT < 3 GeV range. The near-side yield is close to zero in low-activity events and increases with ETPbor Nchrec. The away-side yield shows a similar increase as a function of ETPbor Nchrec, but it starts at a value significantly above zero. The yield difference between these two regions is found to vary slowly with ETPbor Nchrec, indicating

TABLE IV. Summary of relative systematic uncertainties on vn, for n = 1 to 5.

n = 1 n = 2 n = 3 n = 4 n = 5

Residual pair acceptance (%) 1.0–5.0 <0.5 1.0–4.0 7.0–12 7.0–20

ZYAM procedure (%) 0.6 0.3 0.3 0.5 0.6

Tracking efficiency and material (%) 1.0 0.4 0.8 1.2 2.4

Monte Carlo consistency (%) 4.0 1.0 2.0 4.0 8.0

Residual pileup (%) 0–2.0 0–2.0 0–2.0 0–2.0 0–2.0

Uncertainty on scale factor α (%) 8.0–30 0.2–10 0.2–12 0.2–14 1.0–14

Choice of peripheral events

for Nchrec> 160 or ETPb>100 GeV (%) 4.0 1.0 1.0 2.0 4.0

(10)

Per-trigger yield

0 0.2 0.4

< 1 GeV

a

0.5 < pT

< 3 GeV

b

1 < pT

| < 5 η Δ 2 < |

0 0.2 0.4 0.6 0.8

< 3 GeV

a

1 < pT

220

rec

), Nch

φ Δ

corr( Y

φ) Δ

peri( Ycorr

α

recoil Y

0 0.5 1

< 4 GeV

a

3 < pT

0 0.5 1

< 5 GeV

a

4 < pT

φ Δ

0 1 2 3

Per-trigger yield

0 0.5 1

< 7 GeV

a

5 < pT

ATLAS p+Pb = 5.02 TeV sNN

28 nb-1

Lint

φ Δ

0 1 2 3

0 0.5 1

1.5 7 < paT < 9 GeV

φ Δ

0 1 2 3

0 0.5 1 1.5

< 12 GeV

a

9 < pT

FIG. 5. The per-trigger yield distributions Ycorr(φ) and Yrecoil(φ) for events with Nchrec 220 in the long-range region |η| > 2. The distributions are shown for 1 < pbT< 3 GeV in various paTranges. They are compared to the recoil contribution estimated from a peripheral event class defined by ETPb< 10 GeV using a rescaling procedure [see Eq. (6) and discussion around it]. The curves are Fourier fits including the first five harmonics.

that the growth in the integrated yield with increasing event activity is similar on the near side and the away side. This behavior suggests the existence of an away-side long-range component that has a magnitude similar to the near-side long-range component.

Figure 7 also shows (solid lines) the recoil component estimated from the low event-activity class (ETPb< 10 GeV)

via the rescaling procedure discussed in Sec.III D. The yield difference between the away side and the near side in this pT

range is reproduced by this estimate of the recoil component.

In other pT ranges, a systematic difference between the recoil component and the yield difference is observed and is attributed to the contribution of a genuine dipolar flow, v1,1, to the correlation function (see discussion in Sec.IV C).

[GeV]

a

pT

0 5 10

intY

0 0.2 0.4 0.6

/3 (a) π

|<

φ Δ

|

< 3 GeV

b

1 < pT

| < 5 η Δ 2 < |

260

rec

Nch

< 260

rec

Nch

220

< 220

rec

Nch

180

< 180

rec

Nch

140

< 140

rec

Nch

110 ATLAS p+Pb

28nb-1 int = 5.02 TeV, L sNN

[GeV]

a

pT

0 5 10

intY

0 0.5 1

/3 (b) π

|>2 φ Δ

|

FIG. 6. (Color online) Integrated per-trigger yields Yintas a function of paTfor 1 < pbT< 3 GeV, for events in various Nchrecranges on (a) the near side and (b) the away side. The errors bars and shaded bands represent the statistical and systematic uncertainties, respectively.

Cytaty

Powiązane dokumenty

51b High Energy Physics Institute, Tbilisi State University, Tbilisi, Georgia. 52 II Physikalisches Institut, Justus-Liebig-Universita¨t Giessen,

51b High Energy Physics Institute, Tbilisi State University, Tbilisi, Georgia. 52 II Physikalisches Institut, Justus-Liebig-Universität Giessen,

51b High Energy Physics Institute, Tbilisi State University, Tbilisi, Georgia. 52 II Physikalisches Institut, Justus-Liebig-Universit¨at Giessen,

51b High Energy Physics Institute, Tbilisi State University, Tbilisi, Georgia. 52 II Physikalisches Institut, Justus-Liebig-Universität Giessen,

51b High Energy Physics Institute, Tbilisi State University, Tbilisi, Georgia. 52 II Physikalisches Institut, Justus-Liebig-Universität Giessen,

51b High Energy Physics Institute, Tbilisi State University, Tbilisi, Georgia. 52 II Physikalisches Institut, Justus-Liebig-Universität Giessen,

51b High Energy Physics Institute, Tbilisi State University, Tbilisi, Georgia. 52 II Physikalisches Institut, Justus-Liebig-Universita¨t Giessen,

51b High Energy Physics Institute, Tbilisi State University, Tbilisi, Georgia. 52 II Physikalisches Institut, Justus-Liebig-Universit¨at Giessen,