One-dimensional charged kaon femtoscopy in p-Pb collisions at √sNN = 5.02 TeV

This is a self-archived version of an original article. This version may differ from the original in pagination and typographic details. Author(s): Title: Year: Version: Copyright: Rights: Rights url: Please cite the original version: CC BY 4.0 https://creativecommons.org/licenses/by/4.0/ One-dimensional charged kaon femtoscopy in p-Pb collisions at √sNN = 5.02 TeV © 2019 CERN, for the ALICE Collaboration Published version ALICE Collaboration ALICE Collaboration. (2019). One-dimensional charged kaon femtoscopy in p-Pb collisions at √sNN = 5.02 TeV. Physical Review C, 100(2), Article 024002. https://doi.org/10.1103/PhysRevC.100.024002 2019 PHYSICAL REVIEW C 100, 024002 (2019) One-dimensional charged kaon femtoscopy in p-Pb collisions at √sNN =5.02 TeV S. Acharya et al.∗ (ALICE Collaboration) (Received 5 April 2019; published 22 August 2019) The correlations of identical charged kaons were measured in p-Pb collisions at √sNN =5.02 TeV by the ALICE experiment at the LHC. The femtoscopic invariant radii and correlation strengths were extracted from one-dimensional kaon correlation functions and were compared with those obtained in pp and Pb-Pb collisions at √s=7 TeV and √sNN =2.76 TeV, respectively. The presented results also complement the identical-pion femtoscopic data published by the ALICE collaboration. The extracted radii increase with increasing chargedparticle multiplicity and decrease with increasing pair transverse momentum. At comparable multiplicities, the radii measured in p-Pb collisions are found to be close to those observed in pp collisions. The obtained femtoscopic parameters are reproduced by the EPOS 3 hadronic interaction model and disfavor models with large initial size or strong collective expansion at low multiplicities. DOI: 10.1103/PhysRevC.100.024002 I. INTRODUCTION The Bose-Einstein enhancement of the production of two identical pions at low relative momenta (or, in other words, quantum statistics correlations) was first observed in the pp annihilation more than 50 years ago [1]. These correlations encode information about the space-time structure of the interaction region of particles created in collisions at kinetic freeze-out (particle-emitting source) [2–4]. Since that time the correlation method has been developed [5,6] and it is now known as correlation femtoscopy. Femtoscopy measures the apparent width of the distribution of the relative separation of emission points, which is conventionally called the radius parameter. The method was successfully applied to the measurement of the space-time characteristics of particle production processes at high energies in particle [7,8] and heavy-ion collisions (see, e.g., Refs. [4,9] and references therein). Identical boson correlations, especially of identical charged pions, have been used extensively over the years to experimentally study properties of the emitting source created in various collision systems [10]. Identical charged kaon femtoscopy studies were also carried out, for example, in Au-Au collisions at √sNN =200 GeV by the STAR [11] and PHENIX [12] collaborations and in ppcollisions at √s= 7 TeV and Pb-Pb collisions at √sNN =2.76 TeV by the ALICE collaboration [13,14]. The study of femtoscopic correlations in asymmetric collision systems is particularly interesting because it provides a bridge between small (pp) and large (A-A) collision systems, and may lead to additional constraints on model scenarios, ∗Full author list given at the end of the article. Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. which were successfully used to describe pp and A-Acollisions. The A-Afemtoscopy results are interpreted within the hydrodynamic framework as a signature of collective radial flow [10,15–17]. Attempts to describe the pp data in the same framework have not been successful so far and it is speculated that additional effects related to the uncertainty principle may play a role in such small systems [18]. The results obtained in asymmetric collisions are difficult to interpret unambiguously. For instance, the femtoscopic study of the data obtained at RHIC for d-Au collisions [19,20] suggest that a hydrodynamic evolution may be present in such a system, while at the LHC the ALICE three-pion [21] and three-dimensional two-pion [22] analyses in p-Pb collisions at √sNN =5.02 TeV demonstrate the more important role of the initial state shape and size of the created system. The excellent particle identification capabilities of the ALICE detector [23] and the data sample collected in p-Pb collisions at √sNN =5.02 TeV in 2013 allow one to perform the K±K±femtoscopic analysis. Kaons are a convenient tool to study Bose-Einstein correlations because they are less influenced by resonance decays than pions and therefore more effectively probe femtoscopic correlations of directly produced particles. The comparison of kaon and pion correlation radii [14,21,22] as a function of pair transverse momentum kT=|pT,1+pT,2|/2 or transverse mass mT=√k2 T+m2, where pT,1(pT,2) is the transverse momentum of the first (second) particle and mis the kaon or pion mass, allows one to understand the collective dynamics (collective flow) of the source created in high-energy collisions. In particular, in the system created by colliding heavy ions, the decrease of the correlation radii with increasing kT(mT) is usually considered as a manifestation of the strong collective expansion of the matter created in such collisions. If the dependence of the interferometry radii on pair momentum in p-Pb collisions followed the trends seen in heavy-ion collisions, it would be an indication of collectivity or the creation of a hot and dense system expanding hydrodynamically [24–26]. In 2469-9985/2019/100(2)/024002(15) 024002-1 ©2019 CERN, for the ALICE Collaboration S. ACHARYA et al. PHYSICAL REVIEW C 100, 024002 (2019) addition, comparing kaon femtoscopic results in pp,p-A, and A-Acollision systems can provide experimental constraints on the validity of hydrodynamic [24,26] and/or color glass condensate [27,28] approaches proposed for the interpretation of the p-Pb data. In this work, the kaon femtoscopic radii in p-Pb collisions at √sNN =5.02 TeV are shown as a function of kTand multiplicity, and are compared with those in pp and Pb-Pb collisions at √s=7TeV[13] and √sNN =2.76 TeV [14], respectively. The presented data are also compared with the EPOS 3.111 model [26], an event generator based on a (3 +1)-dimensional [(3 +1)D] viscous hydrodynamical evolution starting from flux tube initial conditions, which are generated in the Gribov-Regge multiple scattering framework. The approach contains a full viscous hydrodynamical simulation and a more sophisticated treatment of nonlinear effects in the parton evolution by considering individual (per Pomeron) saturation scales than in previous EPOS versions [25,29]. There are also changes in the core-corona procedure [30] crucial in proton-nucleus collisions, so that the initial energy of the flux tubes is separated into a part, which constitutes the initial conditions for hydrodynamic expansion (core) and the particles, which leave the matter (corona). This model reasonably reproduces multiplicity distributions, transverse momentum spectra, and flow results, and it gives the best description of kaon spectra [25,29,31–33]. The paper is organized as follows: Sec. II shortly describes the ALICE experimental setup and charged kaon selection criteria used in the presented work. In Sec. III,thefemtoscopic correlation function analysis is described in detail and the sources of systematic uncertainties are discussed. The extracted radii and correlation strengths are shown and compared with model predictions in Sec. IV. The obtained results are summarized in Sec. V. II. DESCRIPTION OF THE EXPERIMENT AND DATA SELECTION A. Experiment The ALICE detector and its performance in the LHC Run 1 (2009–2013) are described in Refs. [23,34], respectively. About 55 ×106p-Pb collision events collected in 2013 at a center-of-mass energy per nucleon-nucleon pair of √sNN = 5.02 TeV were analyzed in this work. Given the energies of the colliding pand Pb beams, the nucleon-nucleon center-ofmass system is shifted with respect to the ALICE laboratory system by 0.465 units of rapidity in the direction of the proton beam. Throughout this paper ηrepresents the pseudorapidity measured in the laboratory frame. The analyzed events were classified according to their multiplicity [35,36] using the measured energy deposition in the V0 detectors [37], which consist of two arrays of scintillators located along the beam line installed on each side of the interaction point and covering 2.8<η<5.1 (V0A, located on the Pb-remnant side) and −3.7<η<−1.7(V0C)[38]. Charged kaons were reconstructed with the central barrel detectors placed inside a solenoidal magnet providing a 0.5 T field parallel to the beam direction, namely the Time Projection Chamber (TPC) [39] and the Inner Tracking System (ITS) [23]. The primary vertex was obtained from the ITS. Its position along the beam direction (the zposition) was required to be within ±10 cm of the center of the ALICE detector to ensure uniform tracking performance. The TPC was used to reconstruct tracks and their momenta. The TPC is divided by the central electrode into two halves, each of which is composed of 18 sectors (covering the full azimuthal angle) with 159 pad rows placed radially in each sector. A track signal in the TPC consists of space points (clusters), each of which is reconstructed in one of the pad rows. The TPC covers an acceptance of |η|<0.8 for tracks, which reach the outer radius of the detector and |η|<1.5 for shorter tracks. The parameters of the track were determined by performing a Kalman fit to a set of clusters with an additional constraint that the track passes through the primary vertex. The quality of the fit is required to have χ2/NDF less than 2. The transverse momentum of each track was determined from its curvature in the uniform magnetic field. The track selection criteria based on the quality of the track reconstruction fit and the number of detected tracking points in the TPC [34,40] were used to ensure that only well-reconstructed tracks were considered in the analysis. Particle identification (PID) for reconstructed tracks was carried out using the TPC together with the time-of-flight (TOF) [40] detector. The TOF is a cylindrical detector consisting of 18 azimuthal sectors divided into five modules along the beam axis at a radius r≃380 cm. The active elements are multigap resistive plate chambers. 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, charge, and momentum. For PID based on TOF information, the particle mass was used to calculate the expected time-of-flight as a function of track length and momentum. For each PID method, the signal for each reconstructed particle is compared with the one expected for a kaon taking into account the detector resolution. The allowed deviations (nσ) depend on the momentum of the particle [14,34,40]. B. Charged kaon selection Track reconstruction for the charged kaon analysis was performed using the signals in the TPC. To ensure a good momentum resolution, each track was required to be composed of at least 70 out of the 159 TPC clusters. Tracks were selected based on their distance of closest approach (DCA) to the primary vertex, which was required to be less than 2.4 cm in the transverse plane and less than 3.0 cm in the longitudinal direction. The kinematic range for kaons selected in this analysis is 0.14 <pT<1.5GeV/cand |η|<0.8. Charged tracks with momentum p<0.5GeV/cwere identified as kaons if they satisfied the requirement nσ,TPC <2 in the TPC. Tracks with p>0.5GeV/cwere required to match to a signal in the TOF, and satisfy nσ,TPC <3aswellasthefollowing momentum-dependent nσselection: nσ,TOF <2for0.5< p<0.8GeV/c,nσ,TOF <1.5for0.8<p<1.0GeV/cand nσ,TOF <1for1.0<p<1.5GeV/c. All selection criteria are listed in Table I. The estimation of purity for p<0.5GeV/cwas performed by parametrizing the TPC dE/dx distribution of the 024002-2 ONE-DIMENSIONAL CHARGED KAON FEMTOSCOPY … PHYSICAL REVIEW C 100, 024002 (2019) TABLE I. Charged kaon selection criteria. pT0.14 <pT<1.5 GeV/c |η|<0.8 DCAtransverse to primary vertex <2.4cm DCAlongitudinal to primary vertex <3.0cm nσ,TPC (for p<0.5 GeV/c)<2 nσ,TPC (for p>0.5 GeV/c)<3 nσ,TOF (for 0.5<p<0.8 GeV/c)<2 nσ,TOF (for 0.8<p<1.0 GeV/c)<1.5 nσ,TOF (for 1.0<p<1.5 GeV/c)<1.0 Number of track points in TPC ⩾70 experimental data in momentum slices and computing the fraction of particle species that could mistakenly contribute to the kaon signal [34]. The use of momentum-dependent values for nσ,TPC and nσ,TOF was the result of studies to obtain the best kaon purity, defined as the fraction of accepted kaon tracks that correspond to true kaons, while retaining a decent efficiency of the PID. The dominant contamination for charged kaons comes from e±in the momentum range 0.4<p<0.5GeV/c. The parameters of the function, which fits the TPC distribution in momentum slices depend on the fit interval and can be a source of the systematic uncertainty associated with the single purity. The purity for p>0.5GeV/c, where the TOF information was employed, was studied with DPMJET [41] simulations using GEANT [42] to model particle transport through the detector. Based on the results of this study, the nσ,TOF values were chosen to provide a charged kaon purity greater than 99%. The momentum dependence of the single-kaon purity in the region of maximal contamination is shown in Fig. 1(a). The pair purity is calculated as the product of two single-particle purities for pairs with qinv <0.25 GeV/c, where the momenta are taken from the experimentally determined distribution. The obtained K±pair purity is shown in Fig. 1(b) as a function of kT. It can be seen from the figure that, despite the lower purity for single kaons in the range 0.45 <p<0.5GeV/c, the pair purity remains high in the wide kTbins used in the analysis due to the effects of averaging over low-purity 0.45 <p<0.5GeV/cand high-purity p<0.45 GeV/cor p>0.5GeV/cbins in the full single-kaon momentum range. The systematic uncertainties of the single purity values lead, in turn, to systematic uncertainties of the obtained pair purity. The analysis was performed in three event multiplicity classes [35,36,43]: 0–20%, 20–40%, and 40–90% and two pair transverse momentum kTbins: (0.2–0.5) and (0.5–1.0) GeV/c. The multiplicity was determined based on the sum of the signal amplitudes of V0A and V0C detectors, commonly referred to as V0M. Table II shows the corresponding mean charged-particle multiplicity densities dNch/dηaveraged over |η|<0.5 using the method presented in Ref. [35]. The dNch/dηvalues were not corrected for trigger and vertexreconstruction inefficiency, which is about 4% for nonsingle diffractive events [44]. At least one particle in the event had to be reconstructed and identified as a charged kaon. The correlation signal was constructed from events having at least two identical charged kaons. Events with a single kaon were included in the event mixing procedure to determine the reference distribution. The femtoscopic correlation functions (CFs) of identical particles are sensitive to two-track reconstruction effects because the considered particles are close in momentum and have close trajectories. Two kinds of two-track effects were investigated. Track splitting occurs when one track is mistakenly reconstructed as two. Track merging is the effect when two different tracks are reconstructed as one. To remove these effects, pairs with relative pseudorapidity |η|<0.02 and relative azimuthal angle |ϕ∗|<0.045 were rejected. The modified azimuthal angle ϕ∗takes into account the bending of the tracks inside the magnetic field and was calculated at a radial distance of 1.2 m [45]. )c (GeV/p 0.3 0.35 0.4 0.45 0.5 0.55 purity ± single K 0.5 0.6 0.7 0.8 0.9 1 pairs ± K ± K (a) = 5.02 TeV NN s Pb -ALICE p 20%-0 40%-20 90%-40 )c (GeV/ T k 0.2 0.4 0.6 0.8 pair purity ± K 0.8 0.85 0.9 0.95 1 = 5.02 TeV NN sPb -ALICE p pairs ± K ± K (b) 20%-0 40%-20 90%-40 FIG. 1. (a) Single and (b) pair K±purities for different event multiplicity classes. The systematic uncertainties associated with the purity correction are shown as boxes. Statistical uncertainties are negligible. The momentum p(kT) values for lower multiplicity classes (blue and green symbols) are slightly offset for clarity. 024002-3 S. ACHARYA et al. PHYSICAL REVIEW C 100, 024002 (2019) TABLE II. V0M event classes and their corresponding dNch/dη[35]. The given uncertainties are systematic only since the statistical uncertainties are negligible. Event class dNch/dη,|η|<0.5 0–20% 47.3 ±0.7 20–40% 24.3 ±0.7 40–90% 17.3 ±1.5 III. ANALYSIS TECHNIQUE The correlation function of two particles with momenta p1 and p2is defined as a ratio C(p1,p2)=A(p1,p2) B(p1,p2)(1) of the two-particle distribution in the given event A(p1,p2) to the reference distribution B(p1,p2)[46]. The reference distribution is formed by mixing events containing at least one charged kaon, where each event is mixed with five other events, which have similar zposition of the primary vertex and similar multiplicity [10]. The mixed particles come from events for which the vertex positions in the beam direction agree within 2 cm and the multiplicities do not differ by more than 1/4 of the width of the given multiplicity class. The correlation function is measured as a function of the invariant pair relative momentum qinv =√|q|2−q2 0, where q0=E1−E2and q=p1−p2are determined by the energy components E1,E2and momenta p1,p2of the particles, respectively. The correlation function is normalized to unity such that C→1 in the absence of a correlation signal. The obtained correlation function Craw was also corrected for purity before the fit [14,47] according to Ccorrected =(Craw −1+P)/P,(2) where the pair purity Pis taken from Fig. 1(b). A. Correlation function parametrization The CFs can be parametrized by various formulas depending on the origin of correlations between the considered particles. The pairwise interactions between K±K±that form the basis for femtoscopy are quantum statistics and the Coulomb interaction. Strong final-state interactions between kaons are negligible [48]. Assuming a Gaussian distribution of a particle source in the pair rest frame, the fit of the kaon CF is performed using the Bowler-Sinyukov formula [49,50] C(qinv)=N1−λ+λK(r,qinv) ×1+exp −R2 invq2 inv D(qinv).(3) The factor K(r,qinv) describes the Coulomb interaction with a radius r,D(qinv) parametrizes the baseline including all nonfemtoscopic effects, for instance resonance decays, and Nis a normalization coefficient. The Coulomb interaction is determined as K(r,qinv)=CQS+Coulomb CQS ,(4) 0.5 1 1.5 = 5.02 TeVsPb -ALICE p pairs KK c<0.5 GeV/ T k0.2< 20%– 0 0.5 1 1.5 c<0.5 GeV/ T k0.2< 40%– 20 ) fit inv q(C ) fit inv q(D 0 0.5 1 0.5 1 1.5 c<0.5 GeV/ T k0.2< 90% – 40 experiment EPOS c<1.0 GeV/ T k0.5< 20% – 0 c<1.0 GeV/ T k0.5< 40% – 20 0 0.5 1 c<1.0 GeV/ T k0.5< 90% – 40 )c (GeV/ inv q ) inv qC( FIG. 2. K±K±experimental correlation functions corrected for purity according to Eq. (2) (red points) and EPOS 3 model baselines [26] (black points) versus pair relative invariant momentum qinv. The CFs are presented in three event multiplicity classes: 0–20%, 20–40%, and 40–90% and two pair transverse momentum kTbins: (0.2–0.5) and (0.5–1.0) GeV/c. The black line shows the fit of EPOS 3 by a first-order polynomial for 0 <qinv <1.0 GeV/c.Thered line shows the subsequent fit of the CF up to qinv <0.5 GeV/cby Eq. (3). The CFs are normalized to unity in the range 0.5<qinv < 1.0 GeV/c. Statistical (lines) and systematic uncertainties (boxes) are shown. where CQS is a theoretical CF calculated with pure quantum statistical (QS) weights (wave function squared) and CQS+Coulomb corresponds to quantum statistical plus Coulomb weights [49,51]. The parameters Rinv and λdescribe the size of the source and the correlation strength, respectively. B. Fitting procedure The parameters Rinv and λcan be extracted using Eq. (3) with various assumptions to handle the nonfemtoscopic baseline Dfrom background effects outside the femtoscopic peak region. There are various methods to deal with the baseline. The simplest way is to assume that it is flat, D(qinv)=1, which can be reasonable in cases where nonfemtoscopic effects are negligible. As can be seen in Fig. 2, the baseline of the obtained experimental functions is not flat. Therefore, it would be more reasonable to describe it, for instance, by a first-order polynomial function D(qinv)=N(1 + aqinv), which reproduces this baseline slope. It interpolates the baseline behavior at high qinv taking into account all nonfemtoscopic effects, which make it nonflat. Then being extrapolated to low qinv, it is supposed to imitate the existing nonfemtoscopic effects. The most natural way to describe the 024002-4 ONE-DIMENSIONAL CHARGED KAON FEMTOSCOPY … PHYSICAL REVIEW C 100, 024002 (2019) baseline is to use Monte Carlo (MC) models where events are generated from physical considerations and contain all but QS and Coulomb effects. A suitable MC model has to reasonably describe the baseline at high qinv, where nonfemtoscopic effects are significant, and also contain nonfemtoscopic effects at low qinv.TheEPOS3[26] model without QS and Coulomb interaction effects included was used to describe the baseline D(qinv). As seen from Fig. 2, EPOS 3 describes the experimental CF outside the correlation peak. The extracted parameter values depend on the fit range, which should be chosen taking into account the characteristic width of the femtoscopic effect observed. In this analysis, the EPOS 3 baseline was fit with a first-order polynomial in 0 <qinv <1.0GeV/c (to flatten statistical uncertainties) and then the experimental CF was fit with Eq. (3) in the range 0 <qinv <0.5GeV/c. The Coulomb interaction radius was set to r=1.5 fm, which is on average close to the extracted radii values. C. Systematic uncertainties The effects of various sources of systematic uncertainty on the extracted parameters were studied as functions of multiplicity and kT. The systematic uncertainties were estimated by varying the selection criteria used for the events, particles, and pairs (with variation limits up to ±20%). The influence of the fit range was estimated by variation of the qinv upper limit by ±40%. Another source of systematic uncertainty is the misidentification of particles and the associated purity correction. A ±10% variation of the parameters (Sec. II B) used for the purity correction estimation was performed. To reduce the electron contamination, the PID criteria were tightened, in particular by extending the momentum range where the TOF signal was used and the energy-loss measurement was required to be consistent with the kaon hypothesis within nσ,TPC <1. There is also an uncertainty associated with the choice of the radius of the Coulomb interaction. It was set to 1.5 fm as a result of averaging of the three radii values that were extracted from the respective multiplicity bins and varied by ±0.5 fm. The relative difference was taken as a systematic TABLE III. Minimum and maximum uncertainty values for various sources of systematic uncertainty (in percent), the punctuation “–” means that the contribution from the given source is negligible. Note that each value is the minimum-maximum uncertainty from a specific source, but can pertain to different multiplicity or kTbins. Thus, the maximum total uncertainties are smaller than (or equal to) the sum of the maximum individual uncertainties shown in this table. Systematic uncertainties whose statistical significance level exceeds 50% were included in the total systematic uncertainty value. Rinv (%) λ(%) Single particle selection 0–1.5 0–3.2 PID and purity – 0–0.6 Pair selection 0–3 0–6 Baseline 1–4.8 0.2–4.1 Fit range 0.6–7 0.5–5.7 Coulomb function 0–2.3 1.8–3.8 Momentum resolution 0–1 0–1 uncertainty. Uncertainties associated with momentum resolution were estimated using a MC simulation with the DPMJET 3.05 [41] model. The effect is limited to low pair relative momentum, where it smears the correlation function and is especially pronounced for narrow femtoscopic peaks. In p-Pb collisions the qinv region of the femtoscopic effect is one order of magnitude wider than the region affected by this inefficiency and, consequently, the corresponding uncertainty is minor. As was explained in Sec. III B, the fitting procedure requires knowledge of the nonfemtoscopic background shape and magnitude. In this analysis, the EPOS 3 model was used for this purpose. The systematic uncertainty associated with the baseline was estimated using an alternative MC model, DPMJET, as well as the two methods based on the use of polynomials described in Sec. III B. Table III presents the uncertainty range for all considered sources of systematic uncertainty, where the minimum (maximum) was chosen from all available values in all multiplicity and kTbins. For each source and each multiplicity and )c (GeV/ T k 0.4 0.6 0.8 (fm) inv R 1 1.5 2ALICE EPOS w/ casc EPOS w/o casc 20%- 0 40%-20 90%-40 = 5.02 TeV NN sPb -p pairs ± K ± K (a) )c (GeV/ T k 0.4 0.6 0.8 λ 0.5 1 ALICE EPOS w/ casc EPOS w/o casc 20%- 0 40%-20 90%-40 = 5.02 TeV NN sPb -p pairs ± K ± K (b) FIG. 3. (a) Experimental K±K±invariant radii Rinv and (b) correlation strengths λshown versus pair transverse momentum kTfor three multiplicity classes and compared with the EPOS 3 model predictions with and without the hadronic cascade phase. Statistical (lines) and systematic uncertainties (boxes) are shown. The points for lower multiplicity classes (blue and green symbols) are slightly offset in the x direction for clarity. 024002-5 S. ACHARYA et al. PHYSICAL REVIEW C 100, 024002 (2019) )c (GeV/ T k 0.2 0.4 0.6 0.8 1 (fm) inv R 1 1.2 1.4 1.6 1.8 ALICE pairs ± K ± K = 17.2〉 ch N〈 = 7 TeV, s pp = 17.3〉 ch N〈 = 5.02 TeV, NN s Pb -p FIG. 4. Comparison of femtoscopic radii, as a function of pair transverse momentum kT, obtained in pp [13]andp-Pb collisions. Statistical (lines) and systematic uncertainties (boxes) are shown. kTbin, the maximum deviation from the parameters obtained with the optimal data selection criteria and fitting methods was taken and applied symmetrically as the uncertainty. The limited data sample for p-Pb collisions leads to quite high statistical uncertainty values and most of the systematic uncertainty contributions were found to be much smaller than the quadratic difference of the statistical uncertainties. Therefore, the systematic uncertainty values were added in quadrature, considering only those whose statistical significance level exceeded 50% [52]. As can be seen from Table III,themain sources of systematic uncertainty on the extracted parameters are the pair selection criteria, the influence of the fit range, the radius of the Coulomb interaction, and the baseline description. All of them contribute to the uncertainty associated with the radii. The extracted correlation strengths have higher statistical uncertainties than the radii and, consequently, for them the pair selection criteria is the only source of systematic uncertainty, which exceeds the statistical significance level chosen in this analysis. IV. RESULTS AND DISCUSSION The extracted Rinv and λparameters are depicted in Figs. 3(a) and 3(b), respectively. Statistical and systematic uncertainties as described in Sec. III C are shown for all results. Figure 3also shows comparisons with the EPOS 3 model (with femtoscopic effects included [29,53,54]) for the same collision system and energy in the same multiplicity and kTbins as the experimental data. Two cases are considered, one with and another one without the hadronic cascade (UrQMD) phase [55]. The EPOS 3 calculations for the radii without the cascade exhibit practically no kTdependence and do not describe the data, while the data are well reproduced by the full EPOS 3 model calculations thereby showing the importance of the hadronic cascade phase at LHC energies. This observation agrees with the conclusion from the threedimensional K±femtoscopic analysis in Pb-Pb collisions at √sNN =2.76 TeV [47] where the hydrokinetic model HKM [56] with the hadronic rescattering phase described the charged kaon femtoscopic radii well. The extracted experimental λvalues are about 0.45, whereas the EPOS 3 ones are about 0.65, i.e., apparently larger than the experimental λvalues. The value of the λparameter may be influenced by non-Gaussian features of the correlation function [57], by a finite coherent component of kaon emission [51,58] and also the contribution of kaons from K∗decays (≈ 50 MeV, where is the decay width) and from other longlived resonances [59]. The reason for the difference between the experimental correlation strengths and those obtained with EPOS 3 could be that the model does not accurately account for all contributions of kaons from various resonance decays [60]. Another explanation could be a partial coherence of the real emitting source [50,58,61,62], which is not taken into account in the EPOS 3 model. In Fig. 4, the radii from pp collisions at √s=7TeV[13] and p-Pb collisions at √sNN =5.02 TeV at similar multiplicity are compared as a function of pair transverse momentum kT. The corresponding radii in Pb-Pb collisions at √sNN = 2.76 TeV [14] are not shown because they were obtained for multiplicities, which are not available in this study. The figure shows that at the same multiplicity, the radii in p-Pb collisions are consistent with those in pp collisions within uncertainties. The statistical significance of this observation (4–15%) does not allow this result to be precisely compared with the results of the one-dimensional three-pion cumulants [21] and three-dimensional two-pion [22] analyses where the 1/3 〉 ch N〈 0 5 10 15 (fm) inv R 0 2 4 6 8 c = 0.42 GeV/〉 T k〈 = 7 TeV, spp c = 0.36 GeV/〉 T k〈 = 5.02 TeV, NN sPb -p c = 0.35 GeV/〉 T k〈 = 2.76 TeV, NN sPb -Pb ALICE pairs ± K ± K (a) 1/3 〉 ch N〈 0 5 10 15 (fm) inv R 0 2 4 6 8 c = 0.8 GeV/〉 T k〈 = 7 TeV, spp c = 0.7 GeV/〉 T k〈 = 5.02 TeV, NN sPb -p c = 0.75 GeV/〉 T k〈 = 2.76 TeV, NN sPb -Pb ALICE pairs ± K ± K (b) FIG. 5. Comparison of femtoscopic radii, as a function of the measured charged-particle multiplicity density Nch1/3,at(a)lowand (b) high kTobtained in pp [13], p-Pb, and Pb-Pb [14] collisions. Statistical (lines) and systematic uncertainties (boxes) are shown. 024002-6 ONE-DIMENSIONAL CHARGED KAON FEMTOSCOPY … PHYSICAL REVIEW C 100, 024002 (2019) 1/3 〉 ch N〈 0 5 10 15 λ 0 0.5 1 ALICE pairs ± K ± K (a) c = 0.42 GeV/〉 T k〈 = 7 TeV, spp c = 0.36 GeV/〉 T k〈 = 5.02 TeV, NN sPb -p c = 0.35 GeV/〉 T k〈 = 2.76 TeV, NN sPb -Pb 1/3 〉 ch N〈 051015 λ 0 0.5 1 ALICE pairs ± K ± K (b) c = 0.8 GeV/〉 T k〈 = 7 TeV, spp c = 0.7 GeV/〉 T k〈 = 5.02 TeV, NN sPb -p c = 0.75 GeV/〉 T k〈 = 2.76 TeV, NN sPb -Pb FIG. 6. Comparison of correlation strengths λ, as a function of the measured charged-particle multiplicity density Nch1/3,at(a)lowand (b) high kTobtained in pp [13], p-Pb, and Pb-Pb [14] collisions. Statistical (lines) and systematic uncertainties (boxes) are shown. radii in pp collisions were obtained to be 5–15% and 10–20% smaller than those in p-Pb collisions, respectively. Figure 5compares femtoscopic radii as a function of the measured charged-particle multiplicity density Nch1/3,at low [Fig. 5(a)] and high [Fig. 5(b)]kTin pp [13], p-Pb, and Pb-Pb [14] collisions. The obtained radii increase with Nch and follow the multiplicity trend observed in pp collisions. The radii are equal in p-Pb and pp collisions at similar multiplicity within uncertainties. This result could indicate that the dynamics of the source in p-Pb collisions at low multiplicities is similar to that in ppcollisions. In particular, if there is a collective expansion of the sources created in ppand p-Pb collisions, these results indicate that the expansion is not significantly stronger in p-Pb than in pp collisions [24]. As seen from the figure, the radii in p-Pb and Pb-Pb collisions were obtained in very different ranges of multiplicity and cannot be compared at the same Nch. In order to make a stronger conclusion between different collision systems, as )c (GeV/ T k 0 0.5 1 λ 0 0.5 1 1.5 ALICE pairs ± K ± K pp = 7 TeVs Pb–p = 5.02 TeV NN s Pb–Pb = 2.76 TeV NN s >22 ch N20%–010%–0 22– 12 ch N40%–20 30%–10 11– 1 ch N90%–40 50%–30 FIG. 7. The K±K±correlation strengths λin pp [13], p-Pb, and Pb-Pb [14] collisions versus pair transverse momentum kTin all multiplicity and kTbins. Statistical (lines) and systematic uncertainties (boxes) are shown. The data points for lower multiplicity classes (blue and green symbols) are slightly offset in kTwith respect to the highest multiplicity classes (red symbols) for better visibility. was done in the pion correlation analyses [21,22], a larger experimental data set should be considered. Figures 6(a) and 6(b) show the correlation strengths λin pp [13], p-Pb, and Pb-Pb [14] collisions at low and high kT, respectively. All λvalues are less than unity probably due to the influence of long-lived resonances and a non-Gaussian shape of the kaon CF peak. It can be noticed from the figure that the correlation strength parameters in Pb-Pb collisions tend to be higher than those in pp and p-Pb collisions. That could point to a more Gaussian source created in Pb-Pb collisions. Figure 7compares correlation strengths λin pp[13], p-Pb, and Pb-Pb collisions as a function of kTfor all available multiplicity bins. As seen from the figure, the correlation strengths in all multiplicity and all kTbins do not show any noticeable kTor multiplicity dependence. The systematic uncertainty values obtained for the compared collision systems are visibly different since even the same source of uncertainty gives a rather different contribution to the total uncertainty value in every collision system. V. SUMMARY In this work, one-dimensional identical charged kaon correlations were obtained and analyzed for the first time in proton-nucleus collisions, that is in p-Pb at √sNN =5.02 TeV. The source size Rinv and correlation strength λwere extracted from a correlation function parametrized in terms of the invariant pair relative momentum qinv. The obtained radii Rinv decrease with increasing pair transverse momentum kT and with decreasing event multiplicity. This is similar to the behavior of pion radii in the three-dimensional two-pion correlation analysis in p-Pb collisions at √sNN =5.02 TeV and one-dimensional three-pion cumulant results in pp collisions at √s=7TeV,p-Pb collisions at √sNN =5.02 TeV and Pb-Pb collisions at √sNN =2.76 TeV. The obtained radii Rinv are reproduced well by the EPOS 3 model (including femtoscopic effects) calculations with the hadronic rescattering phase, whose importance was also demonstrated in the three-dimensional femtoscopic analysis of K±pair correlations in Pb-Pb collisions. The values of the 024002-7 S. ACHARYA et al. PHYSICAL REVIEW C 100, 024002 (2019) correlation strength parameters λin EPOS 3 are apparently larger than the experimental λvalues, which could be due to coherent sources not incorporated in EPOS 3 and long-lived resonances not taken into account accurately enough in this model. The kaon Rinv values in p-Pb and pp collisions show the same trend with multiplicity. However, it is difficult to say whether the same is true for the Pb-Pb points because of a large gap in multiplicities available in p-Pb and Pb-Pb collisions. The results disfavor models, which incorporate substantially stronger collective expansion in p-Pb collisions compared to pp collisions at similar multiplicity. The correlation strength λdoes not show any trends with multiplicity or kT. The fact that the correlation strength in Pb-Pb collisions tends to be higher than in pp and p-Pb collisions could be an indication of a more Gaussian source created in Pb-Pb collisions. However, a stronger conclusion is prevented due to large statistical and systematic uncertainties, especially for the Pb-Pb data. ACKNOWLEDGMENTS 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 (Yerevan Physics Institute) Foundation (ANSL), State Committee of Science and World Federation of Scientists (WFS), Armenia; Austrian Academy of Sciences, Austrian Science Fund (FWF): [M 2467-N36] and Nationalstiftung für Forschung, Technologie und Entwicklung, Austria; Ministry of Communications and High Technologies, National Nuclear Research Center, Azerbaijan; Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Universidade Federal do Rio Grande do Sul (UFRGS), Financiadora de Estudos e Projetos (Finep) and Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Brazil; Ministry of Science & Technology of China (MSTC), National Natural Science Foundation of China (NSFC) and Ministry of Education of China (MOEC), China; Croatian Science Foundation and Ministry of Science and Education, Croatia; Centro de Aplicaciones Tecnológicas y Desarrollo Nuclear (CEADEN), Cubaenergía, Cuba; Ministry of Education, Youth and Sports of the Czech Republic, Czech Republic; The Danish Council for Independent Research - Natural Sciences, the Carlsberg Foundation and Danish National Research Foundation (DNRF), Denmark; Helsinki Institute of Physics (HIP), Finland; Commissariat à l’Energie Atomique (CEA), Institut National de Physique Nucléaire et de Physique des Particules (IN2P3) and Centre National de la Recherche Scientifique (CNRS) and Rlégion des Pays de la Loire, France; Bundesministerium für Bildung, Wissenschaft, Forschung und Technologie (BMBF) and GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany; General Secretariat for Research and Technology, Ministry of Education, Research and Religions, Greece; National Research, Development and Innovation Office, Hungary; Department of Atomic Energy Government of India (DAE), Department of Science and Technology, Government of India (DST), University Grants Commission, Government of India (UGC) and Council of Scientific and Industrial Research (CSIR), India; Indonesian Institute of Science, Indonesia; Centro Fermi - Museo Storico della Fisica e Centro Studi e Ricerche Enrico Fermi and Istituto Nazionale di Fisica Nucleare (INFN), Italy; Institute for Innovative Science and Technology, Nagasaki Institute of Applied Science (IIST), Japan Society for the Promotion of Science (JSPS) KAKENHI and Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan; Consejo Nacional de Ciencia (CONACYT) y Tecnología, through Fondo de Cooperación Internacional en Ciencia y Tecnología (FONCICYT) and Dirección General de Asuntos del Personal Academico (DGAPA), Mexico; Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO), Netherlands; The Research Council of Norway, Norway; Commission on Science and Technology for Sustainable Development in the South (COMSATS), Pakistan; Pontificia Universidad Católica del Perú, Peru; Ministry of Science and Higher Education and National Science Centre, Poland; Korea Institute of Science and Technology Information and National Research Foundation of Korea (NRF), Republic of Korea; Ministry of Education and Scientific Research, Institute of Atomic Physics and Ministry of Research and Innovation and Institute of Atomic Physics, Romania; Joint Institute for Nuclear Research (JINR), Ministry of Education and Science of the Russian Federation, National Research Centre Kurchatov Institute, Russian Science Foundation and Russian Foundation for Basic Research, Russia; Ministry of Education, Science, Research and Sport of the Slovak Republic, Slovakia; National Research Foundation of South Africa, South Africa; Swedish Research Council (VR) and Knut & Alice Wallenberg Foundation (KAW), Sweden; European Organization for Nuclear Research, Switzerland; National Science and Technology Development Agency (NSDTA), Suranaree University of Technology (SUT) and Office of the Higher Education Commission under NRU project of Thailand, Thailand; Turkish Atomic Energy Agency (TAEK), Turkey; National Academy of Sciences of Ukraine, Ukraine; Science and Technology Facilities Council (STFC), United Kingdom; National Science Foundation of the United States of America (NSF) and United States Department of Energy, Office of Nuclear Physics (DOE NP), United States of America. [1] G. Goldhaber, S. Goldhaber, W.-Y. Lee, and A. Pais, Influence of Bose-Einstein statistics on the anti-proton proton annihilation process, Phys. Rev. 120,300 (1960). [2] G. I. Kopylov and M. I. Podgoretsky, Correlations of identical particles emitted by highly excited nuclei, Sov. J. Nucl. Phys. 15, 219 (1972) [Yad. Fiz. 15, 392 (1972)]. 024002-8 ONE-DIMENSIONAL CHARGED KAON FEMTOSCOPY … PHYSICAL REVIEW C 100, 024002 (2019) 116Technical University of Košice, Košice, Slovakia 117Technische Universität München, Excellence Cluster ‘Universe’, Munich, Germany 118The Henryk Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, Cracow, Poland 119The University of Texas at Austin, Austin, Texas, United States 120Universidad Autónoma de Sinaloa, Culiacán, Mexico 121Universidade de São Paulo (USP), São Paulo, Brazil 122Universidade Estadual de Campinas (UNICAMP), Campinas, Brazil 123Universidade Federal do ABC, Santo Andre, Brazil 124University College of Southeast Norway, Tonsberg, Norway 125University of Cape Town, Cape Town, South Africa 126University of Houston, Houston, Texas, United States 127University of Jyväskylä, Jyväskylä, Finland 128University of Liverpool, Liverpool, United Kingdom 129University of Science and Techonology of China, Hefei, China 130University of Tennessee, Knoxville, Tennessee, United States 131University of the Witwatersrand, Johannesburg, South Africa 132University of Tokyo, Tokyo, Japan 133University of Tsukuba, Tsukuba, Japan 134Université Clermont Auvergne, CNRS/IN2P3, LPC, Clermont-Ferrand, France 135Université de Lyon, Université Lyon 1, CNRS/IN2P3, IPN-Lyon, Villeurbanne, Lyon, France 136Université de Strasbourg, CNRS, IPHC UMR 7178, F-67000 Strasbourg, France, Strasbourg, France 137Université Paris-Saclay Centre d’Etudes de Saclay (CEA), IRFU, Départment de Physique Nucléaire (DPhN), Saclay, France 138Università degli Studi di Foggia, Foggia, Italy 139Università degli Studi di Pavia, Pavia, Italy 140Università di Brescia, Brescia, Italy 141Variable Energy Cyclotron Centre, Homi Bhabha National Institute, Kolkata, India 142Warsaw University of Technology, Warsaw, Poland 143Wayne State University, Detroit, Michigan, United States 144Westfälische Wilhelms-Universität Münster, Institut für Kernphysik, Münster, Germany 145Wigner Research Centre for Physics, Hungarian Academy of Sciences, Budapest, Hungary 146Yale University, New Haven, Connecticut, United States 147Yonsei University, Seoul, Republic of Korea aDeceased. bDipartimento DET del Politecnico di Torino, Turin, Italy. cM.V. Lomonosov Moscow State University, D.V. Skobeltsyn Institute of Nuclear, Physics, Moscow, Russia. dDepartment of Applied Physics, Aligarh Muslim University, Aligarh, India. eInstitute of Theoretical Physics, University of Wroclaw, Poland. 024002-15