Azimuthally-differential pion femtoscopy relative to the third harmonic event plane in Pb-Pbcollisions at root ( NN )-N-S = 2 . 76

Acharya, S.; Acosta, F.T.; Adamova, D.; Adolfsson, J; Aggarwal, MM.; Rinella, G.A.; Agnello, Maria; Agrawal, N.; Ahammed, Z.; Ahn, S.U.; Aiola, S.; Akindinov, A.; Al-Turany, M.; Alam, S.; Albuquerque, DSD; Aleksandrov, D.; Alessandro, B; Molina, Rafael A.; Ali, Yusuf; Alici, A.; Alkin, A.; Alme, J.; Alt, T.; Bearden, Ian; Bilandzic, Ante; bsm989, bsm989; Gajdosova, Katarina; Bourjau, Christian Alexander; Gaardhøje, Jens Jørgen; Nielsen, Børge Svane; Thoresen, Freja; Zhou, You; Christensen, Christian Holm; Chojnacki, Marek; Ozelin De Lima Pimentel, Lais


Introduction
Heavy-ion collisions at LHC energies create a hot and dense medium known as the quark-gluon plasma (QGP) [1]. The QGP fireball first expands, cools, and then freezes out into a collection of final-state hadrons. Correlations among the particles carry information about the space-time extent of the emitting source, and are imprinted on the final-state spectra due to a quantummechanical interference effect [2]. Commonly known as intensity or Hanbury-Brown-Twiss (HBT) interferometry, the correlation of two identical particles at small relative momentum, is an effective tool to study the space-time ("femtoscopic") structure of the emitting source in relativistic heavy-ion collisions [3]. The initial state of a heavy-ion collision is characterized by spatial anisotropies that lead to anisotropies in pressure gradients, and consequently to azimuthal anisotropies in final particle distributions, commonly called anisotropic flow. Anisotropic flow is usually characterized by a Fourier decomposition of the particle azimuthal distribution and quantified by the flow coefficients v n and the corresponding symmetry plane angles n [4]. Elliptic flow is quantified by the second flow harmonic coefficient v 2 , whereas triangular flow [5] is quantified by v 3 . Due to the position-momentum correlations E-mail address: alice -publications @cern .ch. in particle emission [6], the particles emitted at a particular angle relative to the flow plane carry information about the source as seen from that corresponding direction; these correlations also lead to the HBT radii to be sensitive to the collective velocity fields, from which information about the dynamics of the system evolution can be extracted.
Azimuthally-differential femtoscopic measurements can be performed relative to the direction of different harmonics event planes [7,8]. The measurements of the HBT radii with respect to the first harmonic event plane (directed flow) at the AGS [9] revealed that the source was tilted relative to the beam direction [10]. The HBT radii variations relative to the second harmonic event plane angle ( 2 ) provide information on the pion source elliptic eccentricity at freeze-out. The recent ALICE measurements [11] indicate that due to the strong in-plane expansion the final-state source elliptic eccentricity is more than a factor 2-3 smaller compared to the initial-state. While the HBT radii modulations relative to 2 are defined mostly by the source geometry, the azimuthal dependence of the HBT radii relative to the third harmonic event plane ( 3 ) originate predominantly in the anisotropies of the collective velocity fields -for a triangular, but static source the radii do not exhibit any oscillations [12]. Models studies [13,14] show that the anisotropy in expansion velocity as well as the system geometrical shape can be strongly constrained by azimuthally differential femtoscopic measurements relative to https 3 . The HBT radii oscillations relative to the third harmonic event plane have been first observed in Au-Au collisions at RHIC energy by the PHENIX Collaboration [15]. Unfortunately, due to large uncertainties these measurements did not allow to obtain detailed information on the origin of the observed oscillations.
In this Letter, the first azimuthally-differential femtoscopic measurement relative to the third harmonic event plane in Pb-Pb collisions at √ s NN = 2.76 TeV from the ALICE experiment are presented. We compare our results to the toy-model calculations from [13] to get an insight on the role of the anisotropies in the velocity fields and the system shape. In addition, we compare our results to a 3 + 1D hydrodynamical calculations [14] and a Blast-Wave Model [16] for a quantitative characterization of the final source shape.

Data analysis
The analysis was performed over the data sample recorded in 2011 during the second Pb-Pb running period at the LHC. Approximately 2 million minimum bias events, 29.2 million central trigger events, and 34.1 million semi-central trigger events were used. The minimum bias, semi-central, and central triggers used all require a signal in both V0 detectors [17]. The V0 detector, also used for the centrality determination [18], is a small angle detector of scintillator arrays covering pseudorapidity ranges 2.8 < η < 5.1 and −3.7 < η < −1.7 for a collision vertex occurring at the center of the ALICE detector. The results of this analysis are reported for collision centrality classes expressed as ranges of the fraction of the inelastic Pb-Pb cross section: 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, and 40-50%. Events with the primary event vertex along the beam direction |V z | < 8 cm were used in this analysis to ensure a uniform pseudorapidity acceptance. A detailed description of the ALICE detector can be found in [19,20]. The Time Projection Chamber (TPC) has full azimuthal coverage and allows chargedparticle track reconstruction in the pseudorapidity range |η| < 0.8, as well as particle identification via the specific ionization energy loss dE/dx associated with each track. In addition to the TPC, the Time-Of-Flight (TOF) detector was used for identification of particles with transverse momentum p T > 0.5 GeV/c. The TPC has 18 sectors covering full azimuth with 159 pad rows radially placed in each sector. Tracks with at least 80 space points in the TPC were used in this analysis. Tracks compatible with a decay in flight (kink topology) were rejected. The track quality was determined by the χ 2 of the Kalman filter fit to the reconstructed TPC clusters [21]. The χ 2 per degree of freedom was required to be less than 4. For primary track selection, only trajectories passing within 3.2 cm from the primary vertex in the longitudinal direction and 2.4 cm in the transverse direction were used. Based on the specific ionization energy loss in the TPC gas compared with the corresponding Bethe-Bloch curve, and the time of flight in TOF, a probability for each track to be a pion, kaon, proton, or electron was determined. Particles for which the pion probability was the largest were used in this analysis. This resulted in an overall purity above 95%, with small contamination from electrons in the region where the dE/dx for the two particle types overlap. Pions were selected in the pseudorapidity range |η| < 0.8 and 0.15 < p T < 1.5 GeV/c.
The correlation function C (q) was calculated as where q = p 1 − p 2 is the relative momentum of two pions, A(q) is the distribution of particle pairs from the same event, and B(q) is the background distribution of uncorrelated particle pairs. The background distribution is built by using the mixed-event technique [22] in which pairs are made out of particles from three different events with similar centrality (less than 2% difference), event-plane angle (less than 6 • difference), and event vertex position along the beam direction (less than 4 cm difference). Both the A(q) and B(q) distributions were measured differentially with respect to the third harmonic event-plane angle EP,3 . Note, that measurements relative to EP,3 will smear any contribution from elliptic flow as the elliptic and triangular event planes are uncorrelated [23]. The third harmonic event-plane angle EP,3 was determined using TPC tracks. To avoid auto-correlation each event was split into two subevents (−0.8 < η < 0 and 0 < η < 0.8). Pairs were chosen from one subevent and the third harmonic eventplane angle EP,3 was estimated using the particles from the other subevent, and vice-versa, with the event plane resolution determined from the correlations between the event planes determined in different subevents [4]. Requiring a minimum value in the twotrack separation parameters ϕ * = |ϕ * 1 − ϕ * 2 | and η = |η 1 − η 2 | reduces two-track reconstruction effects such as track splitting or track merging. The quantity ϕ * is defined in this analysis as the azimuthal angle of the track in the laboratory frame at the radial position of 1.6 m inside the TPC. Splitting is the effect when one track is reconstructed as two tracks, and merging is the effect of two tracks being reconstructed as one. Also, to reduce the splitting effect, pairs that share more than 5% of the TPC clusters were removed from the analysis. It is observed that at large relative momentum the correlation function is a constant, and the background pair distribution is normalized such that this constant is equal to unity. The analysis was performed for different collision centralities in several ranges of k T , the magnitude of the pion-pair transverse where ϕ pair is the pair azimuthal angle. The Bertsch-Pratt [3,24] out-side-long coordinate system was used with the long direction pointing along the beam axis, out along the transverse pair momentum, and side being perpendicular to the other two. The three-dimensional correlation function was analyzed in the Longitudinally Co-Moving System (LCMS) [25], in which the total longitudinal momentum of the pair is zero, p 1,L = −p 2,L .
To isolate the Bose-Einstein contribution in the correlation function, effects due to final-state Coulomb repulsion must be taken into account. For that, the Bowler-Sinyukov fitting procedure [26,27] was used in which the Coulomb weight is only applied to the fraction of pairs (λ) that participate in the Bose-Einstein correlation. In this approach, the correlation function is fitted by (2) where N is the normalization factor. The function G(q, ϕ) describes the Bose-Einstein correlations and K (q) is the Coulomb part of the two-pion wave function integrated over a source function corresponding to G(q). In this analysis the Gaussian form of G(q, ϕ) [28] was used where the parameters R out , R side , and R long are traditionally called HBT radii in the out, side, and long directions. The cross-terms R 2 os , R 2 sl , and R 2 ol describe the correlation in the out-side, side-long, and out-long directions, respectively.
The systematic uncertainties on the extracted radii, discussed below, vary in k T and centrality. They include uncertainties related to the tracking efficiency and track quality, momentum resolution, different values for pair cuts ( ϕ * and η), and correlation function fit ranges [29]. Similarly to the azimuthally inclusive analysis [29], different pair cuts were used, with the default values chosen based on a Monte Carlo study. The difference in the results from using different pair cuts rather than the default pair cuts were included in the systematic uncertainties (1-4%). For different k T and centrality ranges, different fitting ranges of correlation function were used as the width of the correlation function depends on k T and centrality range. The difference in the results from using different fit ranges are due to the contamination of electrons in the particle identification and the non-perfect Gaussian source (1-3%).
We also studied the difference in the results by using positive and negative pion pairs separately as well as data obtained with two opposite magnetic field polarities of the ALICE L3 magnet. They have been analyzed separately and a small difference in the results (less than 3%) has been also accounted for in the systematic uncertainty. The total systematic uncertainties were obtained by adding in quadrature the contributions from all various sources mentioned above. The systematic uncertainty associated with the event plane determination is negligible compared to other systematic uncertainties; the procedure for the reaction plane resolution correction of the results is described in the next section.

Results
Fig . 1 presents the dependence of R 2 out , R 2 side , and R 2 long on the pion emission angle relative to the third harmonic event plane for centrality 20-30% and different k T ranges. Note that R 2 out and R 2 side exhibit in-phase oscillations (for a quantitative analysis, see below). Within the uncertainties of the measurement, R 2 long oscillations, if any, are insignificant. Oscillations of R 2 ol and R 2 sl radii (not shown) are found to be consistent with zero, as expected due to the source symmetry in longitudinal direction, and are not further investigated. The curves represent the fits to the data using the functions [12]: Fitting the radii's azimuthal dependence with the functional forms of Eq. (4) allows us to extract the average radii and the amplitudes of oscillations. The χ 2 per number of degree of freedom is 0.3-1.8 depending on k T and centrality range. The results for the average radii R 2 out,0 , R 2 side,0 , and R 2 os,0 were found to be consistent with those reported previously in [11] in azimuthally inclusive analysis. The extracted amplitudes of oscillations have to be corrected for the finite event plane resolution. There exist several methods for such a correction [30], which produce consistent results [31] well within uncertainties of this analysis. The results shown below have been obtained with the simplest method first used by the E895 Collaboration [9], in which the amplitude of oscillation is divided by the event plane resolution. In this analysis the event plane resolution correction factor is about 0.6-0.7, depending on centrality. Fig. 2 shows the oscillation parameters R 2 out,3 , R 2 side,3 , R 2 long,3 , and R 2 os,3 for different centrality and k T ranges. All radii oscillations exhibit weak centrality dependence, likely reflecting the weak centrality dependence of the triangular flow itself. The k T dependence is different for different radii oscillations: while the magnitudes of R 2 out,3 and R 2 os,3 are smallest for the smallest k T range, it is opposite for R 2 side,3 (and, possibly for R 2 long,3 ), where the oscillations become stronger. The parameter R 2 long,3 is consistent with zero within the systematic uncertainties while R 2 os,3 is positive for all centralities and k T ranges except for the lowest k T range 0.2-0.3 GeV/c. Note that R 2 out,3 and R 2 side,3 are negative for all centralities and k T ranges. In the toy model simulations [13] such phases of radii oscillations were observed only in the socalled "flow anisotropy dominated case" (a circular source with the radial expansion velocity including the third harmonic modulation) and not for "geometry dominated" case (triangular shape source with radial expansion velocity proportional to radial distance from the center, with corners having largest expansion velocity). Fig. 3 shows the relative amplitudes of radius oscillations R 2 out,3 /R 2 side,0 , R 2 side,3 /R 2 side,0 , and R 2 os,3 /R 2 side,0 . Similar to the previous analyses and theoretical calculations [14] we report all the radii oscillations relative to the side radius the least affected by the emission time duration. There exist no obvious centrality dependence. As the average radii decrease with increasing k T , the k T dependence of relative oscillation amplitudes appear much stronger for "out" and "out-side" radii, while "side" radius relative amplitude exhibits no k T dependence with the uncertainties. The shaded bands in Fig. 3 indicate the results of 3+1D hydrodynamical calculations [14]. These calculations assume constant shear viscosity to entropy ratio η/s = 0.08 and bulk viscosity that is nonzero in the hadronic phase ζ /s = 0.04, and the initial density from a Glauber Monte Carlo model. The parameters of the model, were tuned to reproduce the measured charged particle spectra, the elliptic and triangular flow. We find that the relative amplitudes R 2 side,3 /R 2 side,0 agree with these results rather well, while the relative amplitudes R 2 out,3 /R 2 side,0 and R 2 os,2 /R 2 side,0 agree only qualitatively. According to the 3+1D hydrodynamical calculations, the negative signs of R 2 side,3 and R 2 out,3 parameters are an indication that the initial triangularity has been washed-out or even reversed at freeze-out due to triangular flow [14].
To investigate further on the final source shape, we compare our results to the Blast-Wave model calculations [16]. In that model, the spatial geometry of the pion source at freeze-out is parameterized by R(φ) = R 0 1 − ∞ n=2 a n cos(n(φ − n )) , (5) where n 's denote the orientations of the n-th order symmetry planes. The amplitudes a n and the phases n are model parameters. The magnitude of the transverse expansion velocity is parameterized as v t = tanh ρ, where the transverse rapidity ρ [13,16] is ρ(r, φ b ) = ρ 0r 1 + ∞ n=2 2ρ n cos (n(φ b − n )) . (6) Here r = r/R(φ), and φ b (φ) is the transverse boost direction assumed to be perpendicular to the surface of constant r. The results of this model presented below were obtained assuming a kinetic freeze-out temperature of 120 MeV, and maximum expansion rapidity ρ 0 = 0.8, tuned to describe single particle spectra. Fig. 4 shows the relative amplitudes of the radius oscillations   the fit of the model to the data. It is observed that the final source anisotropy is close to zero, significantly smaller than the initial triangular eccentricities that are typically of the order of 0.2-0.3. The   Fig. 5. Blast-Wave model [16] source parameters, final spatial (a 3 ) and transverse flow (ρ 3 ) anisotropies, for different centrality ranges, as obtained from the fit to AL-ICE radii oscillation parameters. The contours represent the one sigma uncertainty. negative values of the final source anisotropy would be interpreted as that the triangular orientation at the initial-state is reversed at freeze out.

Summary
We have reported a measurement of two-pion azimuthallydifferential femtoscopy relative to the third harmonic event plane in Pb-Pb collisions at √ s NN = 2.76 TeV. The observed oscillations of the HBT radii unambiguously indicate a collective expansion of the system and anisotropy in collective velocity fields at freeze-out. Clear in-phase oscillations of R out and R side , with both R 2 out,3 and R 2 side,3 parameters (as defined in Eq. (4)) being negative, have been observed for all centralities and k T ranges. According to model calculations [13] the observed R out and R side in-phase oscillations are characteristics of the source with strong triangular flow and close to zero spatial anisotropy. This conclusion is further confirmed by a detailed comparison of our results with the Blast-Wave model calculations [16], from which the parameters of the source, the spatial anisotropy and modulations in the radial expansion velocity, have been derived, with spatial triangular anisotropy being more than an order of magnitude smaller than the typical initial anisotropy values. The oscillation amplitudes exhibit weak centrality dependence, and in general decrease with decreasing k T except for R 2 side,3 which on opposite is the largest in the smallest k T bin. The results of the 3+1D hydrodynamic calculations [14] are in a good qualitative agreement with our measurements but, quantitatively, the model predicts a stronger dependence of R 2 out,3 oscillations on k T than observed in the data.

Acknowledgements
We thank J. Cimerman and B. Tomasik for providing us with the results of the Blast-Model calculations [16].
The ALICE Collaboration would like to thank all its engineers and technicians for their invaluable contributions to the construction of the experiment and the CERN accelerator teams for the outstanding performance of the LHC complex. The ALICE Collaboration gratefully acknowledges the resources and support provided by all Grid centres and the Worldwide LHC Computing Grid (WLCG) collaboration. The ALICE Collaboration acknowledges the following funding agencies for their support in building and running the ALICE detector: A.I. Alikhanyan National Science Laboratory