One-dimensional pion, kaon, and proton femtoscopy in Pb-Pb collisions at $\sqrt{s_{\rm {NN}}}$ =2.76 TeV

The size of the particle emission region in high-energy collisions can be deduced using the femtoscopic correlations of particle pairs at low relative momentum. Such correlations arise due to quantum statistics and Coulomb and strong final state interactions. In this paper, results are presented from femtoscopic analyses of $\pi^{\pm}\pi^{\pm}$, ${\rm K}^{\pm}{\rm K}^{\pm}$, ${\rm K}^{0}_S{\rm K}^{0}_S$, ${\rm pp}$, and ${\rm \overline{p}}{\rm \overline{p}}$ correlations from Pb-Pb collisions at $\sqrt{s_{\mathrm {NN}}}=2.76$ TeV by the ALICE experiment at the LHC. One-dimensional radii of the system are extracted from correlation functions in terms of the invariant momentum difference of the pair. The comparison of the measured radii with the predictions from a hydrokinetic model is discussed. The pion and kaon source radii display a monotonic decrease with increasing average pair transverse mass $m_{\rm T}$ which is consistent with hydrodynamic model predictions for central collisions. The kaon and proton source sizes can be reasonably described by approximate $m_{\rm T}$-scaling.


Introduction
Two-particle correlations at low relative momenta (commonly referred to as femtoscopy), which are sensitive to quantum statistics (in the case of identical particles) as well as strong and Coulomb finalstate interactions (FSI), are used to extract the space-time characteristics of the particle-emitting sources created in heavy-ion collisions [1][2][3]. The source radii extracted from these correlations describe the system at kinetic freeze-out, i.e. the last stage of particle interactions. Pion femtoscopy, which is the most common femtoscopic analysis, has shown signatures of hydrodynamic flow in heavy-ion collisions, manifesting as a decrease in the source radii with increasing transverse mass m T = k 2 T + m 2 [4,5], where k T = |p T,1 + p T,2 |/2 is the average transverse momentum of the pair. This behavior can be interpreted as one of the signatures of the formation of deconfined quark matter in these collisions [6]. However, a necessary condition for collective behavior is for all particles created in the collision, not just pions, to experience hydrodynamic flow. Thus, femtoscopic studies with particles other than pions are also needed. It was shown that the hydrodynamic picture of nuclear collisions for the particular case of small transverse flow leads to the same m T behavior of the longitudinal radii (R long ) for pions and kaons [7]. This common m T -scaling for π and K is an indication that the thermal freeze-out occurs simultaneously for π and K and that these two particle species are subject to the same Lorentz boost. Previous kaon femtoscopy studies carried out in Pb-Pb collisions at the SPS by the NA44, NA49, and CERES Collaborations [8][9][10] reported the decrease of R long with m T as ∼ m −0.5 T as a consequence of the boost-invariant longitudinal flow. Subsequent studies carried out in Au-Au collisions at RHIC [11][12][13] have shown the same level in the m T -dependencies for π and K radii, consistent with a common freeze-out hyper-surface. Like the SPS analysis, no exact universal m T -scaling for the 3D radii was observed at RHIC. In the case of the one-dimensional correlation radius R inv , only approximate scaling with m T is expected as an additional confirmation of hydrodynamic expansion [4]. In fact, R inv source sizes as a function of m T for different particle types (π, K, p...) follow the common curve with an accuracy of ∼ 10%.
The motivation for comparing femtoscopic analyses with different particle species is not limited to studying m T dependence. The kaon analyses also offer a cleaner signal compared to pions, as they are less affected by resonance decays, while the proton analysis provides a possibility for checking if baryons are included in the collective motion. Studying charged and neutral kaon correlations together provides a convenient experimental consistency check, since they require different detection techniques (charged tracks vs. decay vertex reconstruction) and call for different final-state interaction fitting parametrizations (Coulomb-dominated vs. strong interaction-dominated), yet they are predicted to exhibit the same femtoscopic parameters [14]. In addition to the charged kaon analyses at the SPS and RHIC, neutral kaon correlations were studied in Au-Au collisions at RHIC [15], and ALICE has performed analyses on both charged and neutral kaons in pp collisions [16,17]. Recent pion femtoscopic results were obtained at RHIC [18] and the LHC [5,[19][20][21], and proton femtoscopy has also been previously studied at RHIC [22]. This paper presents the results of femtoscopic studies of π ± π ± , K ± K ± , K 0 S K 0 S , pp, and pp correlations from Pb-Pb collisions at √ s NN = 2.76 TeV by the ALICE experiment at the LHC. The femtoscopic radii and λ parameters (the latter describe the decrease of the femtoscopic correlations due to e.g. long-lived resonances; see Sec. 3.1 and 4) are extracted from one-dimensional correlation functions in terms of the invariant momentum difference for a range of collision centralities and m T values. A hydrokinetic model [14] is used to compare the kaon experimental results with hydrodynamic predictions.
The organization of the paper is as follows. In Sec. 2, we describe the data selection criteria. In Sec. 3, the details of the correlation functions and the fitting process are discussed. The results of the analysis are shown in Sec. 4, and a summary is provided in Sec. 5.

Data analysis
The dataset analyzed in this paper is from Pb-Pb collisions at √ s NN = 2.76 TeV at the LHC measured by the ALICE detector [23]. About 8 million events from 2010 and about 40 million events from 2011 were used (2010 data were analyzed in the pion and K 0 S analyses only). Events were classified according to their centrality determined using the measured amplitudes in the V0 detectors [24]. Charged particle tracking is generally performed using the Time Projection Chamber (TPC) [25] and the Inner Tracking System (ITS) [23]. The ITS allows for high spatial resolution in determining the primary (collision) vertex. In the pion, charged kaon, and proton analyses, the determination of the momenta of the tracks was performed using tracks reconstructed with the TPC only and constrained to the primary vertex. Primary tracks were selected based on the distance of closest approach (DCA) to the primary vertex. Additional track selections based on the quality of the track reconstruction fit and the number of detected "hit" points in the TPC were used. Also, all primary pairs sharing more than 5% of TPC clusters were rejected. In the neutral kaon analysis, the secondary daughter tracks used global (TPC and ITS) track reconstruction and did not use any cuts based on track reconstruction quality or number of used or shared TPC clusters. The secondary vertex finder used to locate the neutral kaon decays employed the "onthe-fly" reconstruction method [26], which recalculates the daughter track momenta during the original tracking process under the assumption that the tracks came from a decay vertex instead of the primary vertex.
Particle identification (PID) for reconstructed tracks was carried out using both the TPC and the Timeof-Flight (TOF) detector [27] in the pseudorapidity range |η| < 0.8. For TPC PID, a parametrized Bethe-Bloch formula was used to calculate the specific energy loss (dE/dx) in the detector expected for a particle with a given mass and momentum. For PID with TOF, the particle mass was used to calculate the expected time-of-flight as a function of track length and momentum. For each PID method, a value N σ was assigned to each track denoting the number of standard deviations between the measured track information and the calculations mentioned above. Different cut values of N σ were chosen based on detector performance for the various particle types and track momentum (see Table 1 for specific values used in each analysis) [28].
The analysis details specific to each particle species used in this study are discussed separately below.

Pion selection
The main single-particle selection criteria used in the pion analysis are summarized in Table 1. Pion identification was performed using the TPC only. An overall purity of the pion candidate sample was estimated using TPC dE/dx distributions of the data and was found to be above 95%. The main source of contamination comes from e ± in the region where the dE/dx curves for pions and electrons intersect.
Femtoscopic correlation functions of identical particles are sensitive to the two-track reconstruction efficiency because the correlated particle pairs (i.e. those with small relative momentum) generally have close trajectories. The main two-track issues are splitting (two tracks reconstructed from one particle) and merging (one track reconstructed from two particles), which are generally avoided using a track separation cut. For pions, pairs were required to have a separation of |∆η| > 0.016 or ∆η 2 + ∆ϕ * 2 > 0.045 measured at the radial distance 1.2 m. Here, η is the pseudorapidity, and ϕ * is the azimuthal coordinate taking into account track bending due to the magnetic field.

Charged kaon selection
The main single-particle selection criteria used in the charged kaon analysis are listed in Table 1. K ± identification was performed using the TPC (for all momenta) and TOF (for p > 0.5 GeV/c) detectors. Figure 1 (a) shows the momentum dependence of the single kaon purity, defined as the fraction of accepted kaon tracks that correspond to true kaon particles. The    dE/dx distributions of the data and by studying HIJING [29] simulations using GEANT3 [30] to model particle transport through the detector. Like the pions, the dominant contamination for charged kaons in the momentum region 0.4 < p < 0.5 GeV/c comes from e ± . The pair purity is calculated as a product of two single-particle purities, where the momenta are taken from the experimental distribution. The K ± pair purity as a function of k T for three different centralities is shown in Fig.1 Regarding two-track selection criteria, charged kaon pairs were required to have a separation of |∆η| > 0.02 or |∆ϕ * | > 0.017 measured at the radial distance 1.6 m.   π + π − pairs showing the K 0 S peak for two centrality intervals. The 45-50% centrality is scaled so that both distributions have the same integral in the range 0.480 < m π + π − < 0.515 GeV/c 2 .

Neutral kaon selection
The decay channel K 0 S → π + π − was used for the identification of neutral kaons. The single-particle cuts for parents (K 0 S ) and daughters (π ± ) used in the decay-vertex reconstruction are shown in Table 1. PID for the pion daughters was performed using both TPC (for all momenta) and TOF (for p > 0.8 GeV/c). Figure 2 shows an example of the π + π − invariant mass distribution where the K 0 S peak is seen. The cuts used in this analysis were chosen to balance statistics and signal purity. The neutral kaon purity (defined as Sig./[Sig.+Bkg.] for 0.480 < m π + π − < 0.515 GeV/c 2 ) was found to be greater than 0.95. Two main two-particle cuts were used in the neutral kaon analysis. To resolve two-track inefficiencies associated with the daughter tracks, such as the splitting or merging of tracks discussed above, a separation cut was employed in the following way. For each kaon pair, the spatial separation between the same-sign pion daughters was tabulated at several points throughout the TPC (every 20 cm radially from 85 cm to 245 cm) and averaged. If the average separation of either pair of tracks was below 5 cm, the kaon pair was not used. Another cut was used to prevent two reconstructed kaons from using the same daughter track. If two kaons shared a daughter track, one of them was cut using a procedure which compared the two K 0 S candidates and kept the candidate whose reconstructed parameters best matched those expected of a true K 0 S particle in two of three categories (smaller K 0 S DCA to primary vertex, smaller daughter-daughter DCA, and K 0 S mass closer to the PDG value [31]). This procedure was shown, using HIJING+GEANT3 simulations, to have a success rate of about 95% in selecting a true K 0 S particle over a fake one. More details about K 0 S K 0 S analysis can be found in Ref. [16,32].

Proton selection
The single-particle cuts used in the proton analysis are summarized in Table 1. The proton analysis used tracks with 0.7 < p T < 4.0 GeV/c. The lower p T cut is used to suppress protons coming from weak decays and interactions with the detector material. Particle identification for p and p was performed using both TPC (for all momenta) and TOF (for p > 0.8 GeV/c). The proton purity was estimated using HIJING+GEANT3 simulations and was found to be greater than 95%. The used DCA criteria do not fully discriminate between primary protons and protons from weak decays. This may lead to a significant contamination from protons from lambda particles. The effect of this contamination is discussed in Sec. 3.4.
Regarding two-track selection criteria, pairs were required to have a separation of |∆η| > 0.01 or |∆ϕ * | > 0.045 measured at the radial distance 1.2 m.

Construction of the correlation functions and fitting procedures
The experimental two-particle correlation function is defined as C(q) = A(q)/B(q), where A(q) is the measured distribution of same-event pair momentum difference, q = p 1 − p 2 , and B(q) is the reference distribution of pairs from mixed events. The pairs in the denominator distribution B(q) are constructed by taking a particle from one event and pairing it with a particle from another event with a similar centrality and primary vertex position along the beam direction. Each event is mixed with five (ten) others for the K 0 S (π ± , K ± , p) analysis. The available statistics of proton pairs with low q (< 0.2 GeV/c) allowed us to perform the analyses only for the one-dimensional correlation function C(q), where q = |q| in the Pair Rest Frame (PRF). In the case of pions and kaons, the statistics were high enough for three-dimensional studies, but these are beyond the scope of this paper; here, only the one-dimensional analysis is presented in order to compare results with heavier particles. The numerator and denominator are normalized such that C(q) → 1 as q → ∞. Pair cuts have been applied in exactly the same way for the same-event (signal) and mixed-event (background) pairs. All correlation functions have been corrected for momentum resolution effects. The correction factors were determined using HIJING events to build simulated correlation functions using theoretical correlation functions as weights. The ratio of the correlation functions using HIJING generated momenta to those using HIJING+GEANT3 reconstructed momenta forms the correction factor.
For the analyses presented in this paper, the theoretical femtoscopic correlation function is defined as the square of the two-particle wavefunction averaged over the relative distance r * of the emitters in the PRF. This is performed using the Koonin-Pratt equation [33,34] For the one-dimensional analysis, we assume a spherically symmetric Gaussian distribution of the particle emitter spatial separation r * in the PRF with size R inv [35], The two-particle wavefunction is (anti-)symmetrized for identical bosons (fermions) and may include terms incorporating Coulomb or strong final-state interactions, depending on the type of particles being studied.
The methods used in constructing and fitting the various correlation functions are discussed separately below.

Pions
Pion correlation functions were fitted using the Bowler-Sinyukov formula [36,37]: where N is the normalization factor. The λ parameter (also used in the other analyses) can be affected by long-lived resonances, coherent sources [19,38,39], and non-Gaussian features of the particle-emission distribution. K(q) is a symmetrized K-factor calculated according to Ref. [19,37] as where C(QS) and C(QS + Coulomb) are the theoretical correlation functions calculated with THERMI-NATOR 2 [40] using the quantum statistics ("QS") and "QS+Coulomb" weights (i.e. squared wavefunction), respectively [41]. The effect of the strong interaction is neglected here, since for like-sign pions, the contribution is small for the expected source sizes [41]. Figure 3 shows an example π + π + correlation function with the corresponding line of best fit. More details about the pion analysis may be found in [42]. Figure 4 shows an example K ± K ± correlation function with the corresponding line of best fit. A purity correction was applied to the correlation function according to

Charged kaons
where the pair purity P is taken from Fig. 1. K ± K ± correlation functions were fitted using the Bowler-Sinyukov formula of Eq. 3; the procedure is essentially the same as for pions. There are no available experimental data for K ± K ± strong FSI. The influence of the strong interaction to the correlation function was estimated with the s-wave scattering length calculated within the fully-dynamical lattice QCD [43]. The systematic uncertainty assigned to this effect was determined to be 4%. Figure 5 shows an example K 0 S K 0 S correlation function with the corresponding line of best fit. K 0 S K 0 S correlation functions were fitted with a parametrization which includes Bose-Einstein statistics as well as strong final-state interactions (FSI) [15,35],

Neutral kaons
where ) c (GeV/ and f (q) is the s-wave scattering amplitude for the K 0 K 0 system; we neglect the scattering for K 0 K 0 and K 0 K 0 due to small scattering lengths ≈ 0.1 fm [15]. The factor of 1/2 in Eq. 8 is due to the fact that half of the K 0 S K 0 S pairs come from K 0 K 0 . The strong FSI have a significant effect on the K 0 K 0 contribution to the K 0 S K 0 S correlation function due to the near-threshold resonances, f 0 (980) and a 0 (980). For the scattering amplitude, only s-wave contributions were taken into account; the higher-order corrections were small and therefore neglected [44]. The scattering amplitude f (q) is calculated using a a two-channel parametrization which accounts for the elastic transition K 0 K 0 → K 0 K 0 and the inelastic transition K + K − → K 0 K 0 (See Ref. [15] for more detailed expressions describing the fit function). Equation 6 also includes an additional factor to account for non-femtoscopic background correlations at large q, with a and b being free parameters in the fit. Figure 6 shows an example pp correlation function with the corresponding line of best fit. The femtoscopic correlations of pp and pp pairs are due to a combination of Fermi-Dirac statistics, Coulomb and strong FSI. A distinct maximum is seen at q ≈ 40 MeV/c [35]; this enhancement is due to the strong interaction, as both quantum statistics and Coulomb interaction present a negative correlation. Due to the fact that feed-down from weak decays cannot be neglected in high-energy heavy-ion collisions, the effects of residual correlations related to the pΛ system are taken into account. The proton daughter of a Λ decay has similar momentum to the Λ itself and may survive the experimental selection for primary protons. Thus, it may contribute to the measured correlations by forming a pair with a primary proton. As can be seen in Fig. 6, attempting to fit the measured correlation functions with the theoretical pp (pp) functions alone was unsuccessful due to the additional positive correlation observed in the range 60 < q < 160 MeV/c. Thus, a method of simultaneous fitting of pp (pp) and pΛ (pΛ) correlations was applied. Contributions from heavier baryon-baryon pairs are not taken into account since the original correlation between the parent particles is not known due to unknown interaction parameters, for example for the ΛΛ pair. Moreover such residual correlations are more smeared compared with pΛ because of larger decay momentum. In addition, the fraction of baryons heavier than Λ decaying to protons is smaller than the fraction of Λ's. Finally, comparing with baryon-antibaryon pairs analysed in [45], the width of the correlation for baryon-baryon pairs is much smaller, and therefore the effect is much more smeared due to decay kinematics.

Protons
The experimental correlation function of pp and pp systems were fitted with [45] where λ pp is the fraction of correlated pp pairs where both particles are primary, and λ pΛ is the fraction of correlated pp pairs where one particle is primary and the other is a daughter of Λ decay. The theoretical proton-proton correlation function was calculated as This formulation takes into account the necessary (anti-)symmetrization of the wavefunction for a pp pair in the singlet (triplet) spin state with a corresponding weight of 1/4 (3/4). The pp pair wavefunction may be written as [44] where r * is the spatial separation of particle emitters at generally different emission moments in the PRF, δ c = arg Γ(1 + iη) is the Coulomb s-wave phase shift, A c (η) = 2πη(e 2πη − 1) −1 is the Gamow factor (also referred to as Coulomb penetration factor), η = ( 1 2 aq pp ) −1 , a = (µz 1 z 2 e 2 ) −1 is the twoparticle Bohr radius taking into account the sign of the interaction (a = 57.6 fm for pp pair), F is the confluent hypergeometric function, ξ = 1 2 q pp r * (1 + cos θ * ), θ * is the angle between q pp and r * ,G is the combination of the regular and singular s-wave Coulomb functions, and ρ = 1 2 q pp r * . The amplitude of the low-energy s-wave elastic scattering due to the short range interaction f c (q pp ) may be expressed as For the feed-down term, the theoretical pΛ correlation function for a given R pΛ transformed into the pp momentum space is obtained from the Lednicky-Lyuboshitz model [35] and calculated as where C pΛ (q pΛ ; R pΛ ) = 1 +C strongFSI (q pΛ ; R pΛ ), and T (q pp , q pΛ ) are the transformation factors related to Λ decay kinematics, calculated with THERMINATOR 2 [40]. Here, a spin-dependent version of Eq. 8 is used [35]: where f S (q) is the spin-dependent scattering amplitude, ρ s is the fraction of pairs in each total spin state S, and d S 0 is the effective radius of the interaction. It is assumed that the radii of pp and pΛ sources are equal. Therefore, there are three free fit parameters in Eq. 10: λ pp , λ pΛ , and R. Theoretical pp and pΛ correlation functions were calculated using several values of the free parameters, and the fit function (for the set of parameters given during each fit iteration) was formed by a quadratic interpolation of the calculated correlation functions.

Systematic uncertainties
The effects of various sources of systematic uncertainty on the extracted fit parameters were studied as functions of centrality and k T . Table 2 shows the minimum and maximum uncertainties from each source. The values of the total uncertainty are not necessarily equal to the sum of the individual uncertainties, as the latter can come from different centrality or k T bins. All four analyses studied the effects of changing the selection criteria for the events, particles and pairs used (variation of cut values up to ±25%) and varying the range of q values over which the fit is performed (variation of q limits up to ±25%). Uncertainties associated with momentum resolution corrections are included in the π, K ± , and p analyses; the K 0 S analysis also studied this and found the uncertainties to be negligible. The K ± , K 0 S , and p analyses encountered uncertainties associated with the non-flat background seen at large-q for high-k T pairs in peripheral collisions (estimated by using different parametrizations (linear or polynomial) to fit the large-q region). Strong FSI uncertainties affect both kaon analyses. For K 0 S , the strong FSI uncertainty /pair selection 2-13  6  3-5  3-5  1-4 2-10 2-27 12-58  Non-flat background  --0.2-5 0.2-5 0-5  0-4  2-3  1-9  Fit range  10  33  1-5  1-5  0-4  0-3  3-26  3-57  Strong FSI  --4  4 1-2 5-10 n/a n/a Coulomb function 3 3 2 4 n/a n/a n/a n/a PID and purity ---5- 28  --4-18 13-41  Momentum resolution  2  3  3-5  5-10  --3 1-8 Fixing λ pp n/a n/a n/a n/a n/a n/a 1-29 n/a R pp /R pΛ ratio n/a n/a n/a n/a n/a n/a 1-13 20-52 Total (quad. sum) 11-21 34 6-9 10-32 2-7 7-15 10-40 30-80 Table 2: Minimal and maximal uncertainty values for various sources of systematic uncertainty (in percent). The λ for the proton analysis refers to the sum of λ pp and λ pΛ . Please note that each value is the minimum (maximum) uncertainty from a specific source, but each can be from a different centrality or k T bin. Thus, the minimum (maximum) total uncertainties are greater (smaller) than (or equal to) the sum of the minimum (maximum) individual uncertainties. 'n/a' denotes that the given descriptor of the systematic uncertainty is not applicable for the specific pair type, and '-' means that the contribution from the given source is negligible.
comes from the fact that several sets of f 0 (980) and a 0 (980) parameters are available [46][47][48][49]; each set is used to fit the data, the results are averaged, and the maximum difference was taken as the systematic error. The π and K ± analyses have uncertainties associated with the choice of the Coulomb function used in the fitting procedure. The K ± analysis had additional uncertainties due to the misidentification of particles and the associated purity correction. The p analysis also had uncertainties associated with the uncertainty in the R pp /R pΛ ratio and attempts to fix the λ pp parameter using the single-particle purity. All of the analyses were performed separately for the two different signs of the ALICE dipole magnetic field, but the resulting systematic uncertainty was found to be negligible in all cases.
Systematic uncertainties on correlation functions (Figs. 4, 5 and 6) were derived from the variation of single-and two-particle cuts.

Results
Figures 7 and 8 present the extracted fit parameters from π ± π ± , K ± K ± , K 0 S K 0 S , and pp correlations for several intervals of centrality and transverse mass. Both statistical and systematic uncertainties are shown. The quality of the fits used to extract the shown parameters can be assessed using the χ 2 /NDF values, which are in the ranges of 1.2-5.0, 0.8-3.5, 0.6-1.5, and 0.8-3.2 for the pion, charged kaon, neutral kaon, and proton analyses, respectively. Figure 7 shows the extracted λ parameters vs. m T for several centralities. The proton λ is the sum of λ pp and λ pΛ from Eq. 10. The values for all species measured lie mostly in the range 0.3-0.7 and show no significant centrality dependence. The values of λ are less than unity due to long-lived resonances which dilute the correlation functions and also lead to non-Gaussian shapes of the correlation functions, especially in the one-dimensional case [20]. Results for kaons and protons are consistent with each other at similar m T . Values of λ for pions are lower than for kaons due to the stronger influence of resonances; an additional cause could be a partial coherence of pions [19]. Figure 8 shows the extracted R inv parameters vs. m T for several centralities. For overlapping m T , the radius parameters are mostly consistent with each other within uncertainties, though the pion radii are generally larger than the kaon radii. The K 0 S radii are slightly higher than K ± radii for central collisions, but the difference is less than the systematic uncertainties. The radius parameters show increasing size  0-10% 10-30% 30-50% Fig. 8: (Color online) R inv parameters vs. m T for the three centralities considered for π ± π ± , K ± K ± , K 0 S K 0 S , pp, and pp. Statistical (thin lines) and systematic (boxes) uncertainties are shown. with increasing centrality as would be expected from a simple geometric picture of the collisions. They also show a decreasing size with increasing m T as would be expected in the presence of collective radial flow [6]. Both of these dependences can be seen in previous π ± π ± femtoscopic measurements [4,5] and also reinforce the interpretation that collective flow is present in these collisions for pions, kaons (neutral and charged), and protons alike. Deviations from exact m T -scaling of R inv can be explained as a consequence of the increase of the Lorentz factor with decreasing particle mass. In a hydrodynamic model [50], scaling is observed for the three-dimensional radii measured in the Longitudinally Co-Moving System (LCMS). The transformation from LCMS to PRF involves a boost along the outward direction only, where the boost value is proportional to the transverse velocity of the pair and inversely proportional to the particle mass (for similar m T ). Thus, a smaller mass leads to an increase in the boosted R out and, subsequently, R inv in the PRF. Indeed, we observe such an effect in the data, as pion radii are systematically higher than kaon radii at the same m T .
A comparison of a hydrodynamic flow + kinetics model, HKM [14], with the measured R inv and λ parameters for 0-5% centrality is shown in Fig. 9. The HKM values in Fig. 9 are specifically from K ± K ± , but the predictions for K 0 S K 0 S and K ± K ± are consistent with each other. For R inv , the charged kaon data show very good agreement with the predictions. The experimental data for the neutral kaons are again slightly higher than for the charged kaons, but this difference is still within systematic uncertainties. For λ , both sets of kaon data match the decreasing trend with increasing k T exhibited by the HKM points, but the model slightly overpredicts the data. It is shown in [14] that the most important resonances for KK pairs, K*(890) and φ (1020), do not significantly influence the λ parameter (due to their low contribution), and the discrepancy between the model and experimental data can be explained by the lower experimental kaon purity and deviations of the experimental correlation function shape from a Gaussian distribution. For protons, the HKM prediction is compatible with the data. HKM calculations for one-dimensional pion radii are currently not available, but three-dimensional radii were reasonably reproduced by this model [51].

Summary
Results from femtoscopic studies of π ± π ± , K ± K ± , K 0 S K 0 S , pp, and pp correlations from Pb-Pb collisions at √ s NN = 2.76 TeV with ALICE at the LHC have been presented. The femtoscopic radii and λ parameters were extracted from one-dimensional correlation functions in terms of the invariant momentum difference. It was found that the emission source sizes of kaons and protons measured in these collisions exhibit transverse mass scaling within uncertainties, which is consistent with hydrodynamic model predictions assuming collective flow. The deviation from the scaling for the pions can be explained as a consequence of the increase of the Lorentz factor with decreasing particle mass during the transformation from LCMS to PRF systems [50]. The extracted λ parameters are found to be less than unity, as is expected due to long-lived resonances and non-Gaussian correlation functions. The predictions of the hydrokinetic model (HKM) for the one-dimensional femtoscopic radii for charged and neutral kaons and protons coincide well with the observations.

Acknowledgements
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