Measurement of the centrality and pseudorapidity dependence of the integrated elliptic flow in lead–lead collisions at \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}}}=2.76$$\end{document}sNN=2.76 TeV with the ATLAS detector

The integrated elliptic flow of charged particles produced in Pb+Pb collisions at \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}}}=2.76$$\end{document}sNN=2.76 TeV has been measured with the ATLAS detector using data collected at the Large Hadron Collider. The anisotropy parameter, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v_2$$\end{document}v2, was measured in the pseudorapidity range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|\eta |\le 2.5$$\end{document}|η|≤2.5 with the event-plane method. In order to include tracks with very low transverse momentum \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_{\mathrm {T}}$$\end{document}pT, thus reducing the uncertainty in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v_2$$\end{document}v2 integrated over \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_{\mathrm {T}}$$\end{document}pT, a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1~\mu \hbox {b}^{-1}$$\end{document}1μb-1 data sample recorded without a magnetic field in the tracking detectors is used. The centrality dependence of the integrated \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v_2$$\end{document}v2 is compared to other measurements obtained with higher \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_{\mathrm {T}}$$\end{document}pT thresholds. The integrated elliptic flow is weakly decreasing with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|\eta |$$\end{document}|η|. The integrated \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v_2$$\end{document}v2 transformed to the rest frame of one of the colliding nuclei is compared to the lower-energy RHIC data.


Introduction
The anisotropy in the azimuthal angle distribution of particles produced in heavy-ion collisions has been studied extensively due to its sensitivity to the properties of the produced hadronic medium [1,2]. The final-state anisotropy arises from the initial spatial asymmetry of the overlap zone in the collision of two nuclei with non-zero impact parameter. The initial spatial asymmetry induces asymmetric pressure gradients that are stronger in the direction of the reaction plane and, due to the collective expansion, lead to an azimuthally asymmetric distribution of the ejected particles. The finalstate anisotropy is customarily characterized by the coefficients v n of the Fourier decomposition of the azimuthal angle distribution of the emitted particles [3]. The second Fourier coefficient v 2 is related to the elliptical shape of the overlap region in non-central heavy-ion collisions, and the higher flow harmonics reflect fluctuations in the initial collision geometry [4]. The first observation of elliptic flow, quantified by measurements of v 2 , at RHIC [5][6][7][8] were found to be well e-mail: atlas.publications@cern.ch described by predictions based on relativistic hydrodynamics [9][10][11], providing compelling evidence that the created matter is strongly coupled and behaves like an almost perfect, non-viscous, fluid. Later studies show small deviations from ideal hydrodynamics, described in terms of the ratio of shear viscosity to entropy density [12][13][14][15].
First results from Pb+Pb collisions at √ s NN = 2.76 TeV [16][17][18][19][20][21] from the Large Hadron Collider (LHC) showed no change in the transverse momentum, p T , dependence of elliptic flow from that measured at the highest RHIC energy, while the elliptic flow integrated over p T [16,20]  This increase in the integrated elliptic flow with energy is therefore driven mostly by the increase in the mean p T of the produced particles. The dependence of elliptic flow on the geometry of the collision (the collision centrality) is of particular importance, since the flow is thought to depend strongly on the initial spatial anisotropy. Hydrodynamical models are used to quantitatively relate the initial geometry to the experimentally measured distributions. Furthermore, recent hydrodynamical calculations [22,23] also include a longitudinal dependence in the source shape, which can be deduced from flow measurements over a wide pseudorapidity range. This article presents measurements of the centrality and pseudorapidity dependence of the elliptic flow integrated over the p T of charged particles produced in Pb+Pb collisions at √ s NN = 2.76 TeV recorded in 2010 by the ATLAS detector.
In order to reduce the uncertainty in the p T -integrated coefficient v 2 by including tracks with p T lower than in the measurements reported by the ALICE [16] and CMS [20] experiments, a special track reconstruction procedure was applied to "field-off" data taken without the solenoid's magnetic field in the tracking detectors. Other track recon-struction methods, applicable at higher p T , were exploited in cross-checks using "field-on" data taken with the solenoid's magnetic field.

The ATLAS detector
The ATLAS detector is a multi-purpose particle physics apparatus and is described in detail elsewhere [24]. This analysis uses the three-level trigger system to select Pb+Pb collision events, the forward calorimeters (FCal) to measure the collision centrality, and the inner detector (ID) to measure charged-particle tracks. The ID tracking system comprises silicon pixel and microstrip detectors and a transition radiation tracker. It provides complete azimuthal coverage and spans the pseudorapidity region |η| < 2.5. 2 The pixel detector consists of a three-layer barrel section and three discs in each of the forward regions. The semiconductor tracker has four double layers of microstrip sensors in its barrel section and nine discs in each of the forward regions. The ID is surrounded by a thin superconducting solenoid, which produces a 2 T axial magnetic field for the field-on data. The FCal measures both electromagnetic and hadronic energy, using copper-tungsten/liquid-argon technology, and provides complete azimuthal coverage for 3.2 < |η| < 4.9. The hardware-based Level-1 trigger selected minimum-bias Pb+Pb collisions by requiring either a coincidence of signals recorded in the zero-degree calorimeters (ZDC) located symmetrically at z = ±140 m (|η| > 8.3) or a signal in at least one side of the minimum-bias trigger scintillators (MBTS) at z = ±3.6 m (2.1 < |η| < 3.9). To suppress beam backgrounds, the Level-2 trigger demanded MBTS signals from opposite sides of the interaction region and imposed a timing requirement on them.
With these trigger conditions, ATLAS recorded a sample of Pb+Pb collisions corresponding to an integrated luminosity of approximately 1 μb −1 taken with the field provided by the solenoid turned off. In addition, approximately 0.5 μb −1 of field-on data was used in studies of track reconstruction performance.

Event selection and centrality definition
The offline event selection required each event to have a vertex formed by at least three charged-particle tracks recon-2 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upward. 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). structed in the ID. The data were recorded at low instantaneous luminosity where the probability of multiple collisions per bunch crossing (pile-up) was negligible. The track reconstruction algorithms therefore allowed only one collision vertex (called the primary vertex) in each event, thereby reducing the processing time while maintaining efficiency. The time difference between the MBTS signals from the opposite sides of the interaction region was required to be less than 3 ns, and a coincidence of ZDC signals was also required. These additional selection criteria efficiently remove beamgas and photo-nuclear interactions. As shown in previous studies [18], the applied trigger and offline requirements provide a minimum-bias event sample in which the fraction of inelastic Pb+Pb collisions is 98 ± 2 %.
Events satisfying the above criteria were also required to have a primary vertex within 50 mm (100 mm) in the zdirection of the nominal centre of the ATLAS detector for the field-off (field-on) data subsample. After requiring all relevant subdetectors to be performing normally, the subsamples used in the analysis of the field-off and field-on data contained approximately 1.6 million and 3 million minimumbias events, respectively.
Monte Carlo (MC) event samples were used to determine the tracking efficiencies and the rates of fake tracks. The HIJING event generator [25] was used to produce minimumbias Pb+Pb collisions. Events were generated with the default parameters except for jet quenching, which was turned off. The effect of elliptic flow was implemented after event generation. The azimuthal angles of final-state particles were redistributed at generator level to produce an elliptic flow consistent with previous ATLAS measurements [18,19]. The simulation of the ATLAS detector's response [26] to the generated events is based on the GEANT4 package [27] and included a detailed description of the detector geometry and material in the 2010 Pb+Pb run. Two samples of 0.5 million MC events were simulated, one with the solenoid field switched off and the other with it switched on. Additional MC samples consisting of 50,000 events simulated with 10-20 % extra detector material were used to study systematic uncertainties. The generated charged particles were reweighted with p T -and centrality-dependent functions so that the p T spectra in the MC samples matched the experimental ones [28].
The centrality of the Pb+Pb collisions was characterized by the summed transverse energy, Σ E FCal T , measured in the FCal [18]. The Σ E FCal T distribution was divided into ten centrality bins, each representing 10 % of the full distribution after accounting for 2 % inefficiency in recording the most peripheral collisions (the 0-10 % centrality interval corresponds to the most central 10 % of collisions: those with the largest Σ E FCal T ). A small change in the overall recording efficiency leads to large fluctuations in the population of the most peripheral collisions. To avoid resulting large sys-tematic uncertainties, the 20 % of events with the smallest Σ E FCal T were not included in the analysis.

Elliptic-flow measurement
The final-state azimuthal anisotropy is quantified by the coefficients in the Fourier expansion of the φ distribution of charged particles [3], where v n and Ψ n are the magnitude and the azimuthal direction (called the event-plane angle) of the n-th flow harmonic, respectively.
The second flow harmonic, v 2 , for a given collision centrality is a function of p T and η, and is determined by where the numerator denotes the average over chargedparticle tracks in a given η and p T range, and the denominator, averaged over events, is a correction accounting for the finite experimental resolution in the determination of the eventplane angle Ψ 2 . This resolution correction was obtained using the sub-event technique [3] as described in Refs. [18,19]. The two "sub-events" defined for each event cover two η ranges of the same width in the positive and negative η hemispheres (3.2 < |η| < 4.8) of the FCal detector. The sub-event-plane angles are determined by where the sums run over transverse energies, E tower T , as measured in calorimeter towers located at negative (N) and positive (P) η in the first sampling layer of the FCal. The FCal towers consist of cells grouped into projective regions in Δη × Δφ of 0.1 × 0.1. The weights, w i (Δη i , Δφ i ) are used to correct for any non-uniformity in the event-averaged azimuthal angle distribution of E tower T . They are determined from the data in narrow Δη i and Δφ i slices.
In the sub-event approach, potential non-flow correlations are minimized by using the reaction plane estimated from the η side opposite to the tracks used for the v 2 measurement; this provides a separation of Δη > 3.2. This method was previously applied [18] to measure v 2 as a function of p T using charged-particle tracks reconstructed in the ID tracking system with a minimum p T of 0.5 GeV.
In order to perform the integration over p T , the differential v 2 measurements are weighted by the number of chargedparticle tracks N corr and summed over bins in η (denoted by the index i) and p T (index k). The number of charged-particle tracks is calculated as N corr where the N i,k is the observed number of tracks in a given η and p T bin, (i, k) is the track reconstruction efficiency and f (i, k) is the estimated rate of fake tracks. In the following sections, the lower limit in the integration of v 2 over p T is denoted by p T,0 .

Track reconstruction
The ID was used to reconstruct charged-particle trajectories. Three track reconstruction methods were applied in order to exploit a large range in particle p T : -the tracklet (TKT) method used for the field-off data in order to reach charged-particle p T below 0.1 GeV [28], -the pixel track (PXT) method used to reconstruct tracks with p T ≥ 0.1 GeV using only the pixel detector in the field-on data sample, -the ID track (IDT) method for the field-on data sample, the default ATLAS reconstruction method, for which all ID sub-detectors are used and the track p T is limited to p T ≥ 0.5 GeV [29].
In the TKT method for field-off data, tracks are formed from the positions of hit clusters in the inner two layers of the pixel detector and the position of the primary vertex reconstructed using ID tracks. In the first step, the η 0 and φ 0 coordinates are defined using the event's vertex position and the hit recorded in the first pixel layer. Then a search for a hit in the second pixel layer (with η 1 and φ 1 coordinates defined with respect to the vertex position) is performed and its consistency with a straight-track hypothesis is checked. Candidate tracklets are required to satisfy the condition where Δη = η 1 − η 0 and Δφ = φ 1 − φ 0 , and σ η (η 0 ) and σ φ (η 0 ) are pseudorapidity-dependent widths of the Δη and Δφ distributions, respectively. In this analysis, N σ = 3 was used as the default condition. Clusters located close to each other in the second pixel layer are most likely to originate from the same particle. Therefore, if more than one cluster located in the second pixel layer fulfils the selection criteria, the resulting tracklets are merged into a single tracklet. The Δη and Δφ distributions in data and MC simulation are compared in Fig. 1. The data and MC distributions agree well.
Candidates fulfilling the criterion in Eq. (5) were accepted for further analysis with η = η 0 and φ = φ 0 . for tracklets measured within the pseudorapidity range |η| < 2, for events in the 0-80 % centrality interval and ΔR < 4σ, 3σ and 2σ (see Sect. 5 for details) as described in the legend This method does not provide information about the track's p T ; nevertheless, its performance can be checked as a function of generator-level particle p T by applying the same reconstruction procedure to the MC simulation and using the p T of the generated particle corresponding to the reconstructed tracklet whenever applicable. Figure 2 compares the p T spectra of stable charged particles at the MC-generator level, N primary , to the spectra of reconstructed tracklets matched to charged particles, N matched , for three representative centrality bins and for |η| < 1. A particle was considered to be primary if it originated directly from the collision or resulted from the decay of a particle with cτ < 1 mm. The matching criterion required that the two hits forming the tracklet be identical to the hits associated with a charged particle. The distributions show that the TKT method is able to reconstruct particles with transverse , to all reconstructed tracklets is also shown momenta ∼0.07 GeV with 50 % efficiency, and that a plateau at about 80 % is reached for p T > 0.1 GeV in all centrality bins. For low p T , the efficiency decreases sharply, but the particle density is small in this region, as is v 2 ; thus the contribution from this region to the integrated elliptic flow is expected to be small. Figure 2 also shows the reconstruction efficiency, N matched /N primary , as a function of η. Here, N primary denotes all primary charged particles with p T ≥ 0.07 GeV, which defines the low-p T limit for integrating v 2 over p T . The efficiency is found to be ∼80 % and depends weakly on η. The rate of fake tracklets, N fake , measured as the ratio of the number of tracklets not matched to charged particles to the total number of reconstructed tracklets, N fake /N reco , increases with centrality and |η|, reaching about 35 % for the most central collisions at |η| = 2. For field-on data, the PXT method allows the transverse momentum range p T > 0.1 GeV to be examined. Tracks were reconstructed within the full acceptance of the pixel detector (|η| < 2.5).  To improve the track reconstruction's performance in the heavy-ion collision environment, the track-quality requirements were made more stringent than those for proton-proton collisions [30]. Pixel tracks were required to have no missing hits in the pixel layers, and the transverse and longitudinal impact parameters, d 0 and z 0 , with respect to the vertex were required to have |d 0 | and |z 0 sin(θ )| less than 1 mm and significances |d 0 /σ d 0 | and |z 0 sin θ/σ z 0 sin θ | less than 3.0. Figure 3 shows good agreement between data and MC simulation in the distributions of |d 0 /σ d 0 | and |z 0 sin θ/σ z 0 sin θ |.
The pixel track method's reconstruction efficiency was evaluated in MC simulation by matching reconstructed tracks to the generated charged particles. A track is matched to a generated charged particle if it is reconstructed from at least 69 % of the pixel hits originating from the latter. The transverse momentum, p T , dependence of the pixel track reconstruction efficiency (left) and the fake rate (right) for three pseudorapidity ranges and three centrality intervals as indicated in the legend struction efficiency on p T in three pseudorapidity ranges and for three selected centrality bins. The efficiency decreases slightly from peripheral to central collisions and also deteriorates when moving away from mid-rapidity. A weak p T dependence is observed above p T > 0.5 GeV for all collision centralities. At lower p T , the efficiency decreases with decreasing p T and to about 20 % at the lowest accessible p T . The fraction of fake tracks, defined as the ratio of reconstructed tracks not matched to generated charged particles to all reconstructed pixel tracks, was evaluated using MC simulation. Figure 4 shows the fake-rate dependence on p T in three pseudorapidity ranges and for three centrality bins. The fake rate is below 10 % for p T above 0.4 GeV and depends very weakly on p T and η for peripheral collisions. In more central collisions, the fake rate increases at low p T and shows a similar increase with increasing |η|.
The performance of the PXT reconstruction method can be compared with that of the IDT method. The track reconstruction efficiency and rate of fake tracks from the IDT method are shown in Fig. 5 (for reconstruction details see Ref. [18]). The minimum p T reached is 0.5 GeV. A comparison of Figs. 4 and 5 shows that the extension towards lower The transverse momentum, p T , dependence of the ID track reconstruction efficiency (left) and the fake rate (right) for three pseudorapidity ranges and three centrality intervals as indicated in the legend p T values for the PTX method is achieved at the expense of much larger fake rates than observed for the IDT method, whereas the reconstruction efficiencies are similar. The two methods have different p T resolutions: it is very good for ID tracks, the root mean square of ( p reco T / p true T − 1) being, in |η| < 1, about 4 % and independent of the track p T in the used range, whereas for pixel tracks it is about 10 % at the lowest p T and increases to about 15 % at 5 GeV.
The performance of the MC simulation in describing the fake rates in the data was checked by comparing the Δη, Δφ, d 0 /σ d 0 and z 0 sin θ/σ z 0 sin θ distributions, like the ones shown in Figs. 1 and 3. Additionally, the distributions of the ratios of the number of tracklets and pixel tracks to the number of ID tracks in data and MC simulation were compared, as shown in Fig. 6. It can be concluded that the MC description of the TKT and PXT methods' performance is adequate.
The elliptic flow depends on the particle type [31] as does the reconstruction efficiency. Although the track reconstruction efficiency is averaged over all particle types in this analysis, the reconstruction efficiencies for simulated pions, kaons and protons are shown as a function of p T in the Appendix. At low transverse momenta, which are the focus of this analysis, the measured v 2 is dominated by pions with negligi-

Corrections to measured v 2
The event-plane method [3] is applied to measure the differential elliptic flow harmonic v 2 (η) in small η bins with the TKT method, and v 2 (η, p T ) in small η and p T bins with the PXT and IDT methods. The differential v 2 measurements are then corrected for detector-related effects. The first correction is associated with the variation in tracking efficiency induced by the flow itself. It is applied only to the PXT method, which is found to be sensitive to the detector occupancy. Such sensitivity is not observed for the IDT method. Since the flow phenomenon is a modulation of the multiplicity, it may induce a variation of the tracking efficiency in an event. Higher occupancy causes lower efficiency, and the number of tracks observed in the event plane is reduced more strongly than the number of tracks observed in other directions. As a consequence, the observed v 2 is smaller. In order to correct for this effect, an appropriate weight was applied to the tracks in the calculation of the numerator of Eq. (2). This weight, the inverted efficiency parameterized as a function of detector occupancy in the vicinity of the track, was derived from MC simulation. In the data, the occupancy was determined for each track from the number of hits near the track in the first layer of the pixel detector. The corrected v 2 ( p T ) was compared to the measurement obtained from the IDT method in the same data. In the MC simulation, the comparison was made to v 2 ( p T ) determined using generated particles. The relative increases in the value of v 2 ( p T ) in data and in simulation were found to be compatible for p T > 0.5 GeV, the range in which the comparison could be performed.
The occupancy correction results in an increase of about 12 % in the integrated v 2 for the 0-20 % centrality interval while it amounts to only 1 % for the most peripheral collisions, when using a lower p T integration limit of p T,0 = 0.1 GeV. For higher values of p T,0 the correction gradually becomes smaller. For p T,0 = 0.5 GeV it decreases to about 7 % for the most central collisions.
An additional correction, applied to the differential measurement of v 2 , accounts for the difference between v 2 measured only with fake tracks and v 2 measured with chargedparticle tracks from the primary vertex. The corrected v 2 is calculated as where v 2,m is the elliptic flow measured with all tracks, v 2, f is the flow of fake tracks, and f is the fake-track rate. This correction was applied to the differential v 2 measured with the TKT, PXT and IDT methods with the corresponding fake rates and v 2, f values. The rate and v 2, f of the fake tracks were derived from MC simulation and then cross-checked in the data with a sample, obtained with inverted track selection criteria, in which fake tracks dominate. Differences between the MC simulation and the data of up to 20 % were observed and included in the systematic uncertainties. The fake tracks reduce the values of v 2 integrated over the p T ranges considered in this analysis. The correction is a function of the fake-track rate and accordingly exhibits a dependence on centrality, p T and η. For |η| < 1, the largest correction, about 15 %, was obtained for the PXT method with p T,0 = 0.1 GeV. For peripheral collisions in the same kinematic range, it decreases to about 11 %. The correction is smaller for higher values of p T,0 . It decreases to about 2 % for p T,0 = 0.5 GeV for the 0-10 % centrality interval and gradually drops to zero for the most peripheral collisions. The fake-track flow correction for the integrated v 2 obtained with the IDT method ( p T,0 = 0.5 GeV) is less than 2 % for the most central collisions and even smaller for the more peripheral ones. For the TKT method, the correction is about 1 % for the most central collisions.
The limited p T resolution for tracks reconstructed in the pixel detector and the rapidly changing dN ch /d p T distribution lead to a significant bin-to-bin migration in p T . As a consequence of the variation of v 2 with p T , v 2 measured in a given p T bin is contaminated by v 2 values of particles from the neighbouring bins. In order to compensate for this effect, a correction derived from MC simulation was applied to the v 2 ( p T ) values. This correction was determined, using pixel tracks matched to generated particles, by comparing the v 2 ( p T ) distribution as a function of reconstructed p T to v 2 ( p T ) as a function of generated p T . In order to validate the correction derived from the MC simulation, the same procedure was applied in the data and in the simulation in the region of p T > 0.5 GeV , where the ID tracks were used instead of the generated particles. The ID tracks were matched by requiring an angular separation (Δη) 2 + (Δφ) 2 < 0.02. A comparison between the corrections obtained in the data and in the MC simulation, as a function of measured p T , showed a good agreement.
The correction for p T -bin migration of the reconstructed tracks was found to be small compared to the occupancy and fake-track flow corrections, and to depend only on the value of p T,0 . It increases the integrated v 2 value by 1 % (1.5 %) for p T,0 = 0.1 GeV (p T,0 = 0.5 GeV) independently of collision centrality.

Uncertainties in the v 2 determination
The systematic uncertainties include those common to different tracking methods, as well as method-specific ones.
The uncertainty which originates from the statistics of the MC samples is treated as a source of systematic uncertainty.
The v 2 values determined for samples enriched in fake tracks in data and MC simulation were compared and differences of up to 20 % for both the PXT and IDT methods were observed. For the PXT method, this difference resulted in a change of v 2 , integrated from p T,0 = 0.1 GeV, for the most central (0-10 %) collisions of 3 % at mid-rapidity and of 15 % at |η| ∼ 2. The impact on the integrated v 2 decreases with increasing centrality. For higher p T,0 values, the change was found to be negligible. For the IDT method, the uncertainty on the v 2 value of fake tracks induces a systematic uncertainty in the integrated v 2 for central collisions of less than 4 % at mid-rapidity and of 9 % at |η| ∼ 2; for peripheral collisions the uncertainty is smaller.
The variation of the fake tracklets' v 2 , at the level of 10 %, obtained from the comparison of data and MC simulation, results in an uncertainty at the level of 2 % in the integrated v 2 across the centrality range 0-40 %.  A comparison of v 2 values obtained with the TKT method for a MC sample with the nominal detector geometry to that with 10 % more active material and 15-20 % more inactive material shows agreement to better than 2 %. Therefore it was assumed that possible inaccuracies in the description of the detector material in the GEANT4 simulation have a negligible effect on the final results. The same holds for the measurements with the PXT and IDT methods.
An overall scale uncertainty on v 2 originates from the uncertainty on the fraction of the total inelastic cross section accepted by the trigger as well as from the event selection criteria, which affects the population of the centrality bins. It is negligibly small (below 1 %) for central collisions and increases to about 6 % for the most peripheral collisions for the TKT method and to about 5 % for both the PXT and IDT methods.
The influence of the detector nonuniformities on the measured v 2 was checked by comparing the v 2 values obtained for positive and negative η. This led to a typical uncertainty of 1 % except for the most peripheral collisions where it increased to about 2 %.
Deviations of sin 2[φ − Ψ 2 ] from zero point to detector non-uniformities and biases in the event-plane determination. The magnitude of the sine term relative to the cosine term is included in the systematic uncertainty of v 2 . For the TKT method, its contribution to the relative systematic uncertainty is negligibly small. For the PXT and IDT methods, it is small for most centrality bins, and increases to 2 % only for the most peripheral collisions.
The analysis procedure was checked with MC studies in which the generated elliptic flow signal was compared to the v 2 values obtained with the same analysis chain as used for the data. In this MC closure test, relative differences of up to 2 % in central collisions and of up to 5 % in peripheral collisions were observed for the TKT method. For the IDT method, the relative difference reaches 2 %; for the PXT method, it remains within 2 % except for the most peripheral collisions where it increases to 5 %. The relative difference between the expected and measured values is included in the total systematic uncertainty.
The ΔR parameter used in the tracklet reconstruction was varied by ±1σ from the nominal value. The resulting variation in the value of v 2 at the level of 1 % is included in the systematic uncertainty. For the PXT and IDT methods, differences between v 2 determined from tracks of negatively and positively charged particles as well as between the baseline v 2 and that obtained with tighter or looser tracking requirements (in which the transverse and longitudinal impact parameter significance criteria are changed by ±1) also contribute to the systematic uncertainty at the level of a few percent.
For the PXT method, the corrections due to the limited p T resolution were varied within their statistical uncertainties and the resulting variation was found to be at the level of 0.5 %, independently of the centrality.
The p T spectrum of charged particles in the MC simulation was reweighted so that the expected detector-level distribution agrees with that observed in the data. This changes the effective fake-track rate and therefore the weights used in the calculation of v 2 . A variation of these weights by up to 10 % has a negligible effect on the determination of v 2 .
The different contributions to the total systematic uncertainty on the integrated v 2 for |η| < 1 are shown in Fig. 7 and summarized in Table 1 for the three tracking methods. The total systematic uncertainties are determined by adding in quadrature all the individual contributions and are treated as ±1σ uncertainties. Figure 8 shows the centrality dependence of v 2 integrated over |η| < 1. For the TKT method, v 2 is integrated over p T > 0.07 GeV. For the PXT method, v 2 is integrated over p T,0 < p T < 5 GeV and p T,0 is varied from 0.1 to 0.5 GeV in steps of 0.1 GeV. Also shown is the v 2 value obtained from the IDT method integrated over 0.5 < p T < 5 GeV. The TKT method with p T,0 = 0.07 GeV gives results consistent with the v 2 values obtained with the PXT method with p T,0 = 0.1 GeV, as could be expected due to the very low charged-particle density and small v 2 signal in the momentum range below 0.1 GeV. This indicates that there is no need to extrapolate the measurements obtained with tracklets down to p T = 0 in order to obtain a reliable estimate of v 2 integrated over the whole kinematic range in p T . Furthermore, for the PXT method such an extrapolation would result in a very small correction to the measured integrated flow, well within the uncertainties of the measurement. This is in contrast to the integrated v 2 with p T,0 chosen at higher values, as also shown in Fig. 8. It can be seen that the integrated v 2 increases almost linearly with p T,0 for p T,0 > 0.1 GeV. Good agreement between the PXT and IDT methods is observed at p T,0 = 0.5 GeV. In Fig. 9, the results of this analysis are compared to the integrated v 2 measured by CMS [20] with p T,0 = 0.3 GeV. In this comparison, the sensitivity to p T,0 is clearly visible. A systematically larger v 2 is observed for the higher value of p T,0 as a consequence of the strong increase of v 2 with increasing p T . The η dependence of the p T -integrated v 2 provides useful constraints on the initial conditions of heavy-ion collisions used in model descriptions of the system's evolution (see, e.g., Refs. [1,2]). Figure 10 shows the η dependence of the p T -integrated v 2 . As already shown in Fig. 9, the differ- Centrality dependence of elliptic flow, v 2 , measured for |η| < 1 and integrated over transverse momenta, p T , for different chargedparticle reconstruction methods as described in the legend. Also shown are v 2 measurements by CMS integrated over 0.3 < p T < 5 GeV and |η| < 0.8 [20] (open crosses). Error bars show statistical and systematic uncertainties added in quadrature ence between the results obtained with p T,0 values of 0.07 and 0.1 GeV is very small and the two measurements agree within uncertainties. The results obtained using the PXT and IDT methods for the same p T,0 are also consistent. The η dependence of the integrated v 2 is very weak. A decrease with increasing |η| of about 5-10 % is seen. A comparison with the results from the CMS experiment [20] is shown in Fig. 11 for the 40-50 % centrality interval. The ATLAS measurements performed with the PXT method were integrated over p T for different p T,0 values, including one adjusted to match that used by CMS. The results agree, within uncertainties, provided the same p T,0 is used.

Results
The different upper limits in the p T integration, 3 GeV for CMS and 5 GeV for ATLAS, have negligible effect on the integrated v 2 value. A systematic decrease in v 2 with decreasing p T,0 is observed as expected from the linear dependence of v 2 on p T for p T ≈ 0. The decreasing value of p T,0 together with that of v 2 makes the integration over the full p T range less sensitive to the uncertainties in the extrapolation down to p T = 0. The large acceptance in η of the ATLAS and CMS experiments makes it possible to study whether the observation of the extended longitudinal scaling of v 2 , when viewed in the rest frame of one of the colliding nuclei, reported by the PHOBOS experiment at RHIC [6,32], holds at the much higher LHC energy. The PHOBOS measurements of elliptic flow over a range of Au+Au collision energies, √ s NN = 19.6, 62.4, 130 and 200 GeV, showed energy independence of the integrated v 2 as a function of |η| − y beam , where y beam = ln ( √ s NN /m) is the beam rapidity and m is the proton mass. A similar effect was also observed in chargedparticle densities [6] and is known as limiting fragmentation [33]. In Fig. 12, the integrated v 2 is plotted as a function of |η| − y beam and compared to the PHOBOS results for three centrality bins matching those used by PHOBOS. The PHO-BOS results are obtained with the event-plane method for charged particles with a low-p T limit of 0.035 GeV at midrapidity and of 0.004 GeV around the beam rapidity [34]. The CMS data [20] obtained with the event-plane method are also shown. The CMS measurement represents v 2 integrated over p T from 0 to 3 GeV. This measurement was obtained by extrapolating v 2 ( p T ) measured for p T > 0.3 GeV and the charged-particle spectra down to p T = 0 under the assumption that v 2 ( p T = 0) = 0 and with the charged-particle yield constrained by the measured dN ch /dη distribution [35]. The ATLAS and CMS results agree within the uncertainties, although the CMS v 2 is systematically smaller by about 5 % than the ATLAS measurement. This small systematic difference can be attributed to the uncertainty in the CMS extrapolation to p T = 0 or the p T threshold of 0.07 GeV for the ATLAS measurement, or the combination of both.
As can be seen from the figure, there is no overlap in |η| − y beam between the PHOBOS and LHC data, so a direct comparison with the low-energy data is not possible. Nevertheless, it can be concluded, keeping in mind the relatively large uncertainties in the low-energy results, that the extrapolation of the trend observed at RHIC to the LHC energy appears to be consistent with the LHC measurements, although the dependence on |η| − y beam may be weaker at the LHC energy.

Summary and conclusions
Measurements of the integrated elliptic flow of charged particles in Pb+Pb collisions at √ s NN = 2.76 TeV are presented by the ATLAS experiment at the LHC. The elliptic anisotropy parameter v 2 is measured with the event-plane method over a broad range of collision centralities (0-80 %). The kinematic range in pseudorapidity extends out to |η| = 2.5, and in p T down to 0.07 GeV. This low-p T region is reached by using a tracklet reconstruction algorithm to analyze about 1 μb −1 of data taken with the solenoid field turned off. Other The CMS result is obtained by averaging the v 2 ( p T ) with the charged particle spectra over the range 0 < p T < 3 GeV. Error bars show statistical and systematic uncertainties added in quadrature track reconstruction methods with low-p T thresholds of 0.1 and 0.5 GeV respectively, are exploited in order to verify the tracklet measurement and provide results that can be directly compared to other LHC measurements. The value of v 2 integrated from p T = 0.07 GeV provides a reliable estimate of the elliptic flow measured over the range p T ≥ 0. The p T -integrated elliptic flow as a function of collision centrality shows a clear dependence on p T,0 , both within the present measurements and in comparison to the CMS results obtained with higher low-p T limits. The integrated elliptic flow increases with centrality, reaching a maximum of 0.08 for mid-central collisions (40-50 %) and then decreases to about 0.02 for the most central collisions.
The pseudorapidity dependence of the p T -integrated v 2 is very weak, with a slight decrease in v 2 as |η| increases. The results are in agreement with the CMS measurements covering the same η range, provided the same low-p T cutoff is used. The integrated v 2 transformed to the rest frame of one of the colliding nuclei is compared to the lower-energy RHIC data. Although a direct comparison is not possible due to the non-overlapping kinematic regions, the general trend observed in the RHIC energy regime seems consistent with the LHC measurements, while the latter may have a weaker dependence on pseudorapidity.
In the low-p T region, the track reconstruction efficiency depends strongly on the particle type. This information is important for comparison of measurements with theory predictions in which the elliptic flow depends on the particle type.
The efficiency of the PXT and TKT methods in reconstructing tracks with |η| < 1 generated as π ± , K ± , p, and p in MC simulation is shown in Fig. 13 as a function of p T . Large differences in efficiency are observed for the PXT method at p T below about 1 GeV and for the TKT method at p T below about 0.4 GeV. Above these values, the reconstruction efficiency is independent of particle type. The efficiency is lowest for p andp. For the TKT method, which is most relevant at low p T , the efficiency for reconstructing protons drops to zero below 0.2 GeV. The transverse momentum, p T , dependence of the TKT (left) and PXT (right) track reconstruction efficiency for π ± , K ± and p ± in the pseudorapidity range |η| < 1 for three centrality intervals, as indicated in the legend 146 (a) Department of Physics, University of Cape Town, Cape Town, South Africa; (b) Department of Physics, University of m Also at CERN, Geneva, Switzerland