Two-particle Bose–Einstein correlations in pp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${ pp }$$\end{document} collisions at s=13\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf {\sqrt{s} = 13}$$\end{document} TeV measured with the ATLAS detector at the LHC

This paper presents studies of Bose–Einstein correlations (BEC) in proton–proton collisions at a centre-of-mass energy of 13 TeV, using data from the ATLAS detector at the CERN Large Hadron Collider. Data were collected in a special low-luminosity configuration with a minimum-bias trigger and a high-multiplicity track trigger, accumulating integrated luminosities of 151 μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu $$\end{document}b-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document} and 8.4 nb-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document}, respectively. The BEC are measured for pairs of like-sign charged particles, each with |η|<2.5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|\eta | < 2.5$$\end{document}, for two kinematic ranges: the first with particle pT>100\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_{\mathrm {T}} > 100$$\end{document} MeV and the second with particle pT>500\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_{\mathrm {T}} > 500$$\end{document} MeV. The BEC parameters, characterizing the source radius and particle correlation strength, are investigated as functions of charged-particle multiplicity (up to 300) and average transverse momentum of the pair (up to 1.5 GeV). The double-differential dependence on charged-particle multiplicity and average transverse momentum of the pair is also studied. The BEC radius is found to be independent of the charged-particle multiplicity for high charged-particle multiplicity (above 100), confirming a previous observation at lower energy. This saturation occurs independent of the transverse momentum of the pair.


Introduction
Particle correlations play an important role in the understanding of multiparticle production in hadron-hadron, hadron-nucleus and nucleus-nucleus collisions. Correlations between two identical bosons, called Bose-Einstein correlations, (BEC), are a well known phenomenon in highenergy and nuclear physics (for reviews see Refs. [1][2][3][4][5][6][7][8]). The BEC are often considered to be the analogue of the Hanbury Brown and Twiss effect [9][10][11] in astronomy, describing the interference of incoherently emitted identical bosons [12][13][14][15]. The BEC constitute a sensitive probe of the space-time geometry of the hadronization region, and allow the determination of the size and shape of the source from which particles are emitted. High-multiplicity data in proton-proton e-mail: atlas.publications@cern.ch interactions can serve as a reference for studies of nucleusnucleus collisions [16]. The effect of BEC is reproduced in both the hydrodynamical/hydrokinetic [17][18][19] and the Pomeron-based [20,21] approaches for modelling of highmultiplicity hadronic interactions. The BEC between a pair of particles depend on the particle multiplicity in the event and on the average transverse momentum of the pair. The dependence on these parameters is particularly sensitive to the space-time features of the hadronization region, as discussed by the authors of various models [22][23][24][25][26][27].
The BEC can be studied for two, three, or more identical bosons, assuming a hadronization region parametrized by one, two, or three size parameters. At the Large Hadron Collider (LHC), the BEC with one size parameter, the source radius, have been studied by the ATLAS Collaboration in pp collisions at the centre-of-mass energies √ s = 0.9 and 7 TeV [28], and in p+Pb collisions at √ s NN = 5.02 TeV per nucleon-nucleon pair [29]. The BEC source radius parameter, R, was observed by ATLAS to saturate, i.e. to reach a plateau level, at high charged-particle multiplicity, in √ s = 7 TeV pp collisions [28]. The BEC in pp collisions have also been studied in the ALICE [30][31][32][33][34], CMS [35][36][37][38] and LHCb [39] experiments. The CMS results at 13 TeV also show a plateau in R at high multiplicity [38].
Pomeron-based models predict a plateauing of R at intermediate multiplicities as a consequence of the overlap of the colliding protons, followed by the return to a lower value of R at high multiplicities [20,21]. The IP-Glasma model [40] predicts some tendency to R saturation at high multiplicities for pp interactions.
In a different approach, two-particle correlations have been analysed at √ s = 7 TeV by ATLAS in the framework of a phenomenological model based on ordered hadron chains. This approach provides an alternative interpretation of the correlations commonly attributed to Bose-Einstein interference [41].
In this paper, studies of the BEC in pp collisions at √ s = 13 TeV using the ATLAS detector [42] are presented. The same methodology is used in this analysis as was used in the previous ATLAS BEC studies [28], and the same minimum-bias (MB) trigger is used as was used in an earlier study of charged-particle distributions at 13 TeV [43,44]. In addition, the study is extended to the region of high-multiplicity collisions by using the high-multiplicity track (HMT) trigger previously used for studies of long-range elliptic azimuthal anisotropies [45].
This paper is organized as follows. The analysis method is described in Sect. 2. The ATLAS detector is described in Sect. 3. The Monte Carlo (MC) and data samples are presented in Sect. 4. The track selection criteria and data correction procedure are presented in Sect. 5 and the reference sample is discussed in Sect. 6. The systematic uncertainties are presented in Sect. 7. The BEC source radius and correlation strength dependence on, separately, charged-particle multiplicity and average transverse momentum of the pair is shown in Sect. 8. This section also presents the doubledifferential dependence of the BEC parameters on multiplicity and average transverse momentum. A comparison with previous, lower-energy ATLAS results is presented, as well as a comparison with CMS results at 13 TeV. The summary and conclusions are given in Sect. 9.

Analysis method
The BEC are observed at small values of the square of the four-momentum difference, Q 2 , between two identical bosons, where Q 2 = −( p 1 − p 2 ) 2 = M 2 12 − 4m 2 , p 1 and p 2 are the four-momenta of the particles, M 12 is the invariant mass of the particle pair and m is the particle mass. The BEC parameters are measured in terms of a two-particle correlation function, where ρ(Q) is the distribution of Q formed from all likesign charged-particle pairs, assuming in this analysis that all charged particles are pions, and ρ 0 (Q) is a reference distribution specially constructed to exclude the effect of BEC. In this analysis, ρ 0 (Q) is constructed from unlike-sign chargedparticle pairs, as discussed in Sect. 6. The distributions ρ(Q) and ρ 0 (Q) are normalized to unity, i.e. they are probability density functions.
To account for the effects of resonances (for example η, η , ω, ρ, K * , φ), the C 2 (Q) correlation function is corrected by dividing by the corresponding quantity constructed from MC simulations that contain resonances but not the BEC: where C data 2 (Q) and C MC 2 (Q) are the C 2 (Q) functions reconstructed from data and MC samples, respectively.
The BEC effect is usually described by a function, , with two parameters: the effective R parameter, and the strength parameter, λ, where the latter is also called the incoherence or chaoticity parameter [46]. If measurement can be made down to a very small value of Q ( 10 MeV), the parameter λ may also be sensitive to a "halo" contribution due to longlived resonances [5], but this is not possible in this analysis. A typical functional form is: In a simplified scheme for fully coherent emission of identical bosons, λ = 0, while for incoherent (chaotic) emission, λ = 1. Various parameterizations of the (λ, Q R) function in Eq. (3) can be found in the literature (see for example Refs. [3,47]), each assuming a different shape for the particleemitting source. The parameter c 0 is a constant depending on the normalizations of ρ(Q) distributions and is chosen such that R 2 (Q) is unity at large Q.
A Lévy-type parameterization for introduces a further free parameter, α: with 0 < α ≤ 2. The Lévy distribution is a simple exponential when α = 1 and a Gaussian when α = 2. Fits to R 2 (Q) with α free usually return a value between 1 and 2, and the fit quality, in terms of χ 2 /nd f , is generally better than the fit using an exponential or Gaussian. However, the interpretation of the R parameter in terms of the distribution of space-time points is then less clear. Also, in studies for the 7 TeV analysis it was found that α varied when R 2 (Q) was studied more differentially in terms of charged multiplicity or pair transverse momentum. So, in common with the other LHC experiments, this analysis constrains α to be 1 or 2. Investigation of the Gaussian parameterization is found to give a significantly poorer description of the data in the low-Q region. (See Sect. 8.1). So, in the study presented here, the BEC parameters are extracted using an exponential parameterization, of a static source assuming a radial Lorentzian distribution of the source. The fitted parameter ε in Eq. (3) takes into account possible long-distance correlations that have not been removed fully by ρ(Q).

The ATLAS detector
ATLAS is a multipurpose particle physics experiment operating at one of the beam interaction points at the LHC. The ATLAS detector 1 covers almost the whole solid angle around the collision point with layers of tracking detectors, calorimeters and muon chambers. It is designed to study a wide range of physics topics at LHC energies. For the measurements presented in this paper, the tracking devices and the trigger system are of particular importance. The innermost part of the ATLAS detector is the inner tracking detector (ID), which has full coverage in φ and covers the pseudorapidity range |η| < 2.5. The ID is immersed in the 2 T axial magnetic field of a superconducting solenoid and measures the trajectories of charged particles. It consists of a silicon pixel detector (Pixel), a silicon microstrip detector (SCT) and a straw-tube transition radiation tracker (TRT), each split into a barrel and two endcap components. The Pixel, SCT and TRT are located around the interaction point spanning radial distances of 33-150 mm, 299-560 mm and 563-1066 mm, respectively. The barrel (each endcap) consists of four (three) pixel layers, four (nine) double layers of silicon microstrips and 73 (160) layers of TRT straws. The Pixel, SCT and TRT have (r, φ)-position resolutions of 10 µm, 17 µm, and 130 µm, respectively.
For Run 2 a new innermost pixel layer, the insertable Blayer (IBL) [48], was installed around a new, narrower and thinner beam pipe. The IBL is composed of 14 lightweight staves arranged in a cylindrical geometry, each made of 12 silicon planar sensors in its central region and 2 × 4 threedimensional sensors at the ends. The IBL pixel dimensions are 50 µm in the φ-direction and 250 µm in the z-direction (compared with 50 µm by 400 µm for the other pixel layers). The intrinsic spatial resolution of the IBL readout is 10 µm in (r, φ)-position and 75 µm in z-position [49]. The smaller radius and the reduced pixel size result in improvements in both the transverse and longitudinal impact parameter resolutions [43,44]. The services for the existing pixel detector were upgraded, significantly reducing the amount of material in the region |η| > 1.5, in particular at the boundaries of the active tracking volume.
A track from a charged particle traversing the barrel detector typically has 12 silicon measurement points (hits), of which 4 are Pixel and 8 are SCT, and more than 30 TRT straw hits. Requirements on an IBL hit and on impact parameters strongly suppress the number of tracks from secondary particles.
The ATLAS detector has a two-level trigger system: the first-level (L1) trigger and the high-level trigger (HLT) [50]. Events used in this analysis were required to satisfy L1 trig-1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upward. Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the beam pipe. The pseudorapidity, η, is defined in terms of the polar angle θ with η = − ln tan(θ/2). gers using the minimum-bias trigger scintillators (MBTS). These are mounted at each end of the detector in front of the liquid-argon endcap-calorimeter cryostats at z = ±3.56 m, and are segmented into two rings in pseudorapidity (2.07 < |η| < 2.76 and 2.76 < |η| < 3.86) with 8 azimuthal sectors in the inner ring and 4 in the outer. The MB events were selected on the basis of the MBTS alone, while the HMT events satisfied additional HLT criteria. Further details are given in Sect. 4.1.
An extensive software suite [51] is used in the reconstruction and analysis of real and simulated data, in detector operations and in the trigger and data acquisition systems of the experiment.

Data samples and event selection
This study uses the pp-collision datasets at √ s = 13 TeV that were used in previous ATLAS studies of MB events [43,44,52] and HMT events [53]. The MB and HMT datasets were taken in a special configuration of the LHC in June 2015, with low beam currents and reduced beam focusing, producing a low mean number of interactions per bunch-crossing, μ : in the range 0.003-0.007 for MB events and 0.01-0.04 for HMT events. This configuration and the event selection used are documented in detail in the ATLAS √ s = 13 TeV MB analysis [43,44]. For the MB data sample, events were selected from colliding proton bunches using the MBTS trigger. This trigger selection required one counter above threshold from either side of the detector, and is referred to as a single-arm trigger. The HMT trigger is a combination of an L1 requirement on the number of hits in the MBTS and an HLT requirement on the multiplicity of HLT-reconstructed charged-particle tracks. This trigger required more than 900 SCT space-points and more than 60 reconstructed goodquality charged tracks with p T > 400 MeV associated with the reconstructed vertex with the highest multiplicity in the event. These selections correspond to integrated luminosities of 151 µb −1 and 8.4 nb −1 for the MB and HMT triggers, respectively.
In the offline selection, track candidates are reconstructed in the Pixel and SCT detectors and extended to include measurements in the TRT [54]. The selected tracks satisfy the following criteria: p T > 100 MeV and |η| < 2.5; at least one Pixel hit and an IBL hit if expected 2 ; at least two, four or six SCT hits for p T < 300 MeV, p T < 400 MeV or p T > 400 MeV, respectively, in order to account for the dependence of track length on p T . All events are required to have at least one vertex [55], formed from a minimum of two tracks. The vertex position is required to be consistent with the beam-spot position within the ATLAS detector [54]. If more than one vertex is reconstructed, the primary vertex (PV) is defined as the vertex with the largest number of associated reconstructed tracks, weighted by the square of their transverse momenta [54]. Selected tracks are then also required to satisfy |d BL 0 | < 1.5 mm, where the transverse impact parameter d BL 0 is calculated with respect to the measured beam line (BL) [56]; and |z BL 0 × sin θ | < 1.5 mm, where z BL 0 is the difference between the longitudinal coordinate of the track along the beam line at the point where d BL 0 is measured and the longitudinal coordinate of the PV, and θ is the polar angle of the track. High-momentum tracks with mismeasured p T are removed by requiring the trackfit χ 2 probability to be larger than 0.01 for tracks with p T > 10 GeV. A special configuration of the track reconstruction algorithms is used in this analysis to reconstruct low-momentum tracks [54]. The track reconstruction efficiency at p T ≈ 100 MeV is about 28% and increases rapidly to almost 75% at p T ≈ 300 MeV [44]. For the MB trigger, the lowp T selected dataset consists of 9.3 × 10 6 events with a total of 2.4 × 10 8 tracks with p T > 100 MeV. The highp T selection requires in addition at least two tracks with p T > 500 MeV. This yields 8.9 × 10 6 events with 1.1 × 10 8 tracks with p T > 500 MeV. For the HMT trigger, which requires at least 60 tracks with p T > 400 MeV at the trigger level, the lowp T selection requires at least 60 tracks with p T > 100 MeV and results in 9.1 × 10 6 events with 9.8 × 10 8 tracks with p T > 100 MeV. For p T > 500 MeV the requirement is at least 40 tracks with p T > 500 MeV; the sample consists of 9.1 × 10 6 events with 5.3 × 10 8 tracks with p T > 500 MeV.
To reduce contamination from multiple proton-proton interactions (pile-up) in the colliding bunches, events with a second vertex containing four or more tracks are removed. Such events represent less than 0.3% of the sample. Secondary vertices can be reconstructed with high efficiency provided the distance, z, between the z-coordinates of the primary and pile-up vertices satisfies | z| ≥ 4 mm. To estimate any possible residual influence of multiple pp interactions in MB-and HMT-triggered data, a dedicated study of the z distribution was performed, using a sample with no pile-up rejection. This distribution was extrapolated to estimate the level of pile-up that will survive the maximum trackmultiplicity limit on a second vertex. For the MB events, the fraction of events with pile-up vertices in the region | z| ≤ 4 mm around the PV is about 0.006%. For the HMT events, it is about 0.05%. 3 The average track multiplicities 3 Since μ is much smaller than 1 for both the MB and HMT samples, the ratio of the number of events with more than one vertex to the per pile-up vertex for the MB and HMT samples are 9.4 and 23, respectively. For the primary vertices these multiplicities are 26 and 108. The fractions of surviving pile-up tracks in the MB and HMT samples are approximately 0.002% and 0.01% per event, respectively, and these levels have a negligible influence on the BEC studies.
The contributions from beam-gas collisions and from non-collision background (cosmic rays and detector noise) were investigated in Ref. [57] and found to be negligible.

Monte Carlo samples
The Pythia 8 MC generator [58] was used to calculate the C MC 2 (Q) correlation functions (Eq. (2)). The generator was used with parameter values set to the A2 tune [59] and with the MSTW2008lo PDF set [60]. Pythia 8 models the effects of colour coherence, which is important in dense parton environments and effectively reduces the number of particles produced in multiple parton-parton interactions. In Pythia 8, the simulation is split into non-diffractive and diffractive processes. The non-diffractive process is dominated by t-channel gluon exchange, comprising approximately 80% of the selected events. The diffractive processes, which are further split into single-and double-diffractive processes, are described by a Pomeron-based approach [61]. The contributions from non-diffractive, single-diffractive and double-diffractive processes were included in proportion to the cross sections predicted by Pythia 8 with the A2 tune. Large MC samples of MB and HMT events were generated and passed through the ATLAS simulation program [62], which is based on Geant4 [63], and the reconstruction chain, which is exactly that used for collision dataset. In addition, all selection criteria applied to the data were also applied as appropriate to the MC samples.
In order to estimate the systematic uncertainty due to the choice of MC model, large MC samples, comparable to the size of the collision dataset, were also generated with three other programs. These were the Pythia 8 Monash tune [64], the EPOS MC generator [65] and Herwig++ [66] using the UE-EE-5 tune [67]. The Monash underlying-event tune is based on the NNPDF2.3lo PDF [68] and incorporates updated fragmentation parameters. EPOS implements a parton-based Gribov-Regge [69] theory, which is an effective field theory describing both hard and soft scattering. The EPOS generator was used with the LHC tune [70]. The Her-wig++ UE-EE-5 tune is based on the CTEQ6L1 PDF. All sets of generated events were passed through the ATLAS simulation program and reconstruction chain.
Footnote 3 continued number of events with at least one vertex is ≈ μ /2. The requirement | z| ≤ 4 mm reduces this fraction by a further factor of ≈ 100.
The purity of the analysis sample of particle pairs is defined as the fraction of identical particle pairs (π ± π ± , K ± K ± , pp andpp) in the sample relative to all like-charge particle (LCP) pairs. The purity of the analysis sample in terms of identical boson pairs is estimated from Pythia 8 with the A2 tune to be about 74% (where about 73% are π ± π ± and about 1% are K ± K ± ), using the selection criteria described above for MB p T > 100 MeV. In studies for the 7 TeV analysis [28] it was found that the impurity reduces the value of λ by a factor that is approximately equal to the impurity of the sample, where a 100% purity has no effect, but the parameter R is unaffected.

Data correction procedure
To correct for the inefficiencies in the trigger and in the reconstruction of tracks and vertices, event-by-event and track-bytrack weighting procedures are applied as appropriate and as used in the previous ATLAS measurements [43,44].
The efficiency of the MB trigger was studied with an independent control trigger [43,44]. The control trigger selects events randomly at L1, which are then filtered in the HLT by requiring at least one reconstructed track with p T > 200 MeV. The MB trigger efficiency increases from 97% for events with low track multiplicity, n sel = 2, to 100% for events with higher multiplicity, n sel ≥ 7, where n sel is the number of tracks reconstructed offline and satisfying the selection criteria given in Sect. 4.1. The efficiency of the HMT trigger is studied using MB-triggered events as the control sample. The efficiency is found to be 15% for events with n sel 85 and reaches a plateau at 100% for n sel 110. The measured MB and HMT trigger inefficiencies as a function of n sel are used to correct each event. These corrections are found to have a negligible impact on the extraction of the BEC parameters discussed in Sect. 8. The effect of the vertex-reconstruction efficiency on the BEC parameters is also negligible.
In addition, each track is assigned a weight that corrects for the track reconstruction efficiency, for the fraction of secondary particles, for the effect of detector resolution smearing primary particles 4 into and out of the fiducial region ( p T > 100 MeV, |η| < 2.5) and for the fraction of fake tracks formed by a random combination of hits [43,44]. The effect of the loss of reconstructed pairs at small opening angles is negligible in comparison with the track reconstruction corrections.
The Q resolution is found to be similar for all MC generators. For the BEC region Q < 0.2 GeV, where the Q distribution rises steeply, the Q resolution is better than 5 MeV and much smaller than the bin size. For Q > 0.2 GeV, where the Q distribution is essentially flat, the resolution increases linearly with Q, rising to 40 MeV for Q 2 GeV. Overall, no unfolding of Q is required.
The track reconstruction inefficiency is such that n sel differs significantly from the true charged-particle multiplicity, n ch . In order to present and study R 2 (Q) (Eq. (2)) as a function of the true charged-particle multiplicity, n ch , the correlation functions R 2 (Q, n sel ) must be unfolded. In the earlier ATLAS MB study at 13 TeV [43,44] an unfolding matrix that reflects the probability of reconstructing n sel charged tracks in an event with a generated charged-particle multiplicity n ch was constructed, using MC simulation with the Pythia 8 ATLAS A2 tune. For the present analysis this unfolding matrix was applied to both the MB and HMT data. It is found that, on average, an event with n sel = 100 reconstructed tracks corresponds to n ch = 127.6 ± 0.1 charged particles at particle level.
As a check of the entire analysis chain, the distribution of C 2 (Q) (Eq. (1)) using the generated particles, C MC,gen 2 (Q), was compared with the C MC,rec 2 (Q) distribution using the simulated detector response to reconstruct the tracks and vertices, which were then corrected for inefficiencies as described above, and unfolded to the particle level. This study was carried out for several n ch and k T intervals, where k T is the average transverse momentum of the pair: In the small-Q region the difference between the unfolded and particle-level distributions varies, depending on n ch and k T , but is a few percent on average. The effect is significant for the first two Q-bins used, 0.02 ≤ Q < 0.04 GeV and 0.04 ≤ Q < 0.06 GeV, and for these bins it is accounted for by adding it in quadrature to the statistical uncertainty in the R 2 (Q) distribution. It is therefore not included in the systematic uncertainties discussed in Sect. 7 and listed in Table 1. A large part of the difference in the first Q-bin, 0.02 ≤ Q < 0.04 GeV, is due to photon conversions in the detector material, which are included in the simulated detector response but are not part of the set of generated particles. Since C 2 is a ratio of like-sign pairs to unlike-sign pairs, the denominator in C MC,rec 2 (Q) contains the contribution of pairs from photon conversions, which show up at small Q, whereas C MC,gen 2 (Q) does not. For MB events the difference for the first bin is 4%, of which 3.5% is the contribution of photon conversions. For HMT events, the corresponding numbers are 3.25% and 2%. Conversion photons are less significant for the HMT events, because the number of track combinations increases quadratically with n ch while the number of conversion photons increases only linearly. The effect of photon conversions in the region Q > 0.04 GeV is negligible. The long-range Coulomb force causes a momentum shift between particles in the LCP and unlike-charge particle (UCP) pairs. The density distributions are corrected for this effect by applying the Gamow penetration factor per track pair with a weight 1/G(Q) such that ρ corr (Q) = ρ(Q)/G(Q) [71][72][73] (for a review see Ref. [74]). The Gamow factor is given by with the dimensionless parameter ζ defined as ζ = ±αm/Q, where α is the electromagnetic fine-structure constant and m is the pion mass. The corrected two-particle density function, ρ corr (Q), is normalized to unity. The sign of ζ is positive for LCP pairs and negative for UCP pairs. The Gamow factor particularly affects pairs with Q < 0.1 GeV, and the correction is 24% for 0.02 ≤ Q < 0.04 GeV, 14% for 0.04 ≤ Q < 0.06 GeV and 10% for 0.04 ≤ Q < 0.06 GeV. The Coulomb interaction is not included in the generated MC event samples, and so the correction is applied only to data events.

Reference sample
A good choice of reference sample is important in obtaining an unbiased estimate of the BEC signal. Ideally, ρ 0 (Q) should include all momentum correlations present in the signal dataset except for those arising from BEC. Several differ-ent methods of constructing an appropriate reference sample have been investigated by various collaborations and authors.
The main methods are: • take UCP pairs from the same event; • take LCP pairs from the same event, but invert the momentum of one of the selected particles: . This is the opposite hemisphere (OHP) technique [7]; • take LCP pairs from the same event, but invert the transverse momentum of one of the selected particles: . This is the z-rotation method [7]; • take LCP pairs from different events, randomly mixed; this is the 'mixed event' (MIX) technique [75].
In addition, most of the recently used methods compare the data with a reference sample having global event parameter values that are the same as, or similar to, those in data -such as the event multiplicity, the invariant mass of the pair, or the rapidity of the pair. For the previous ATLAS analysis [28] the UCP, OHP and MIX procedures were investigated in detail using dedicated MC studies, and it was observed that MIX and OHP behaved similarly. For this analysis, the UCP and OHP were studied. As the MC models do not implement BEC, the singleratio correlation functions C MC 2 (Q) are expected to be close to 1, and to be a rather flat function of Q, particularly in the low-Q, BEC-sensitive region. This is indeed the behaviour observed in the earlier analysis [28] over much of the param-eter space in n ch and k T for the MIX or OHP technique, and for the OHP technique in this analysis. However, in highmultiplicity and high-k T intervals, which are relevant for this analysis, C MC 2 (Q) constructed with the MIX or OHP technique exhibits an increase in the low-Q, BEC-sensitive region.
As an example, Fig. 1a shows the correlation functions using the OHP technique for the reference sample. The correlation functions shown are the single ratios C data 2 (Q) and C MC 2 (Q), and the double ratio R 2 (Q), all for HMT events, and for the average transverse momentum of the pair 1000 < k T ≤ 1500 MeV. The C MC 2 (Q) exhibits a rise similar to that seen in C data 2 (Q) in the BEC-sensitive region, although no such effect should be present in C MC 2 (Q). Such a rise would have a marked effect on the BEC parameters extracted from R 2 (Q). Figure 1a also shows that R 2 (Q) is less than 1 in the region 0.7 < Q < 1.2 GeV.
The use of UCP pairs for the reference sample was also investigated in detail. As in the other techniques, C MC 2 (Q) behaves satisfactorily over much of the parameter space in n ch and k T , although resonances distort the region Q ≈ 0.5 GeV, particularly at low multiplicity or at high k T . An example is shown in Fig. 1b, using the UCP reference sample, and for the same n ch and k T intervals as in Fig. 1a. Although there is some distortion in the resonance region, where C MC 2 (Q) and C data 2 (Q) do not match perfectly, the resultant R 2 (Q) matches well the BEC region of C data 2 (Q). Accordingly, the UCP method is used in this analysis, as it was for the previous ATLAS analysis [28], and the other methods are not suitable.
The UCP reference sample used in this analysis is taken from the same events used to form pairs of like-charge particles, i.e. ρ 0 (Q) ≡ ρ( p ± 1 , p ∓ 2 ). By construction, this UCP sample has the same topology and global parameters as the LCP sample ρ( p ± 1 , p ± 2 ), but is naturally free of any BEC effect. However, the UCP sample contains oppositely charged hadron pairs from the decay of resonances such as η, η , ω, ρ, K * , and φ, which are not present in the LCP combinations. Details concerning resonances are given in Appendix A. These resonances contribute to the low-Q region, and can give a spurious BEC signature with a large effective R [76][77][78][79][80][81][82]. The effect of the resonant states in the Q variable was checked in MC studies, and the most affected regions are excluded from the analysis of the BEC effect, as described in Sect. 8.1.

Systematic uncertainties
In this section the systematic uncertainties in values of the BEC parameters R and λ are summarized. These parameter values are obtained from a fit to the exponential model (4). This fit is applied to the 'inclusive' distributions, correspond-ing to 2 ≤ n ch ≤ 250 for MB and 100 < n ch ≤ 300 for HMT, and all k T . The fit is also applied to selected ranges of n ch or k T , and to the double-differential (n ch , k T )-distributions. The discussion in this section is mainly in terms of the 'inclusive' fits, but the systematic uncertainties have been evaluated for all selections applied to the data.
The systematic uncertainties resulting from the track reconstruction efficiency, which are parameterized as a function of p T and η, were determined in earlier analyses [43,44] and these are implemented as uncertainties in the track weights for particle pairs in the Q distributions that enter the correlation functions. The effects of track splitting and merging have been studied in detail. They are largest for Q values below 0.04 GeV, but even for the first used bin (0.02 ≤ Q < 0.04 GeV) the effect is only 0.4% and is negligible in comparison with other uncertainties.
The effects of uncertainties in the unfolding matrix on the BEC parameters λ and R are found to be negligible for n ch < 200. For the highest-multiplicity MB (201 ≤ n ch < 250) and HMT (231 ≤ n ch < 300) events this systematic uncertainty is very small and not more than 0.01% for λ and 0.03% for R.
The leading source of uncertainty in λ and R is the choice of MC generator used in the calculation of R 2 (Q). Figure 2 shows a selection of plots of the BEC parameters as a function of n ch or k T , for different choices of MC generator. The values of the parameters λ and R obtained using EPOS LHC, Pythia 8 Monash or Herwig++ UE-EE-5 are compared with the parameter values obtained using Pythia 8 A2 in terms of the ratios shown in the lower panel of each plot in Fig. 2.
As can be seen in this figure, there are some regions of n ch and k T where there are pronounced and significant differences between the results obtained using the different generators. In order to assign the systematic uncertainty due to the MC modelling, a conservative approach is adopted: for each generator the weighted deviation is evaluated with respect to Pythia 8 A2, where the weight is inversely proportional to the statistical uncertainty for that generator. The systematic uncertainty is taken as the largest (in magnitude) of these weighted deviations, treating positive and negative deviations separately, which leads to asymmetric uncertainties. This procedure is followed for MB and HMT events separately.
An important source of systematic uncertainty for the determination of the parameter λ arises from the Coulomb correction. This uncertainty is estimated by varying the Gamow factor G(Q) by ±0.15, which is the magnitude of the effect on UCP or LCP pairs at Q = 0.02 GeV. This is the standard variation for a conservative estimate of the systematic uncertainty in the Coulomb correction [28,35], and accounts also for effects that include the extended size of the emission source, as discussed in Ref. [83]. The resulting  Fig. 1 Comparison of single-ratio two-particle correlation functions, C data 2 (Q) and C MC 2 (Q), with the two-particle double-ratio correlation function, R 2 (Q), for the high-multiplicity track (HMT) events using a the opposite hemisphere (OHP) like-charge particles pairs refer-ence sample for k T -interval 1000 < k T ≤ 1500 MeV and b the unlike-charge particle (UCP) pairs reference sample for k T -interval 1000 < k T ≤ 1500 MeV systematic uncertainties for the inclusive fit of the parameter λ are the same for both triggers and are equal to 1.8%. For the parameter R the resulting systematic uncertainties are small compared to the cumulative systematic uncertainties, and equal to 0.1% and 0.3% for MB and HMT events, respectively.
The effect of the upper limit of the fit range in Q is estimated by changing it from the nominal 2 GeV by ±0.1 GeV, which corresponds to a total variation of about five times the nominal Q-resolution for Q ≈ 2 GeV. This gives an estimate of the residual uncertainty due to long-range correlations accounted for by ε, the parameter in the linear term of Eq. (3). Other effects contributing to the systematic uncertainties are the lowest value of Q for the fit, the bin size, and the exclusion of Q regions where resonances are not modelled perfectly, as discussed in Sect. 8.1. These uncertainties are estimated by varying the lowest Q value in the fit by ±0.01 GeV, by changing the bin size by ±0.01 GeV, and by broadening the excluded intervals by 0.02 GeV on both sides.
The systematic uncertainties for the inclusive study are summarized in Table 1. The systematic uncertainties are combined by adding them in quadrature, and the resulting values are given in the bottom row of the table. The cumulative systematic uncertainties are +3.9 −5.9 % and +5.6 −2.6 % for the inclusive determination of λ, and +5.0 −5.8 % and +5.7 −2.5 % for R, for the MB and HMT events, respectively.
The same sources of uncertainty are considered for the single-and double-differential measurements in n ch and k T . The results for the n ch distribution are also shown in Table 1.
The pattern of uncertainties is similar to the inclusive case, but shows a wider spread. Away from the BEC region, the MC description of C data 2 improves as the multiplicity increases, and this is reflected in Table 1 by smaller systematic uncertainties for the HMT results. This is also true for the k T distribution and the double-differential distribution.

Results
The BEC parameters, λ and R, are extracted from a fit to the R 2 (Q) distribution in intervals of n ch and k T . In addition, there are two possible selections on the p T threshold of the final-state particles. In order to give a reasonably complete picture of all of this information, the parameters describing the dependencies of λ and R on variables of interest are themselves fitted to simple functions (usually exponentials), and the parameters of these fits are presented.
In Sect. 8.1 the procedure used to fit the R 2 (Q) distribution is described. The dependence on multiplicity is presented in Sect. 8.2. The comparison with previous ATLAS results at 7 TeV is discussed in Sect. 8.3. The behaviour with k T is presented in Sect. 8.4. The 'double-differential' dependence of the BEC parameters on multiplicity and k T at 13 TeV is explored in Sect. 8 Figure 3 shows examples of the R 2 (Q) correlation functions measured at √ s = 13 TeV for two different multiplicity intervals, and the fits to the Gaussian and exponential forms. In the lower panels of Fig. 3 the differences between the measured R 2 (Q) and the fit to the exponential form, normalized to the experimental uncertainty σ (Q), , are shown for both multiplicity intervals. It is evident from Fig. 3 that the Gaussian function does not give a good description of the low-Q, BEC-sensitive region, Q ≤ 0.1 GeV, while the exponential function matches the data well in this region. Therefore, the quantity R 2 (Q)/σ (Q) is not shown for the Gaussian case, and the Gaussian fit is not considered further. Figure 3 shows structure in R 2 (Q), compared with the smooth fitting functions, in the regions Q ≈ 0.25 GeV and 0.4 ≤ Q ≤ 0.9 GeV. If the MC modelling of resonance production, particularly the η, ρ and ω, does not agree very well with data, then this Q region will be affected, as appears to be the case. Accordingly, the regions 0.2 ≤ Q ≤ 0.3 GeV and 0.4 ≤ Q ≤ 0.9 GeV are excluded from the fits to the R 2 (Q) correlation functions for the MB and HMT events. For the MB events, similar arguments lead to the additional exclusion of the region 1.0 ≤ Q ≤ 1.16 GeV, as it appears to be influenced by the reflection of the f 0 (980) and f 2 (1270) mesons in the case of 100 ≤ k T ≤ 200 MeV and in the multiplicity region n ch ≤ 40. These excluded regions lie outside the BEC-sensitive region and the main effect of the exclusions is to markedly reduce the χ 2 without significantly changing the fitted BEC parameters.
The results for the BEC parameters obtained using the exponential fit with the UCP reference sample for the MB and HMT samples, summed over all multiplicities for p T > 100 MeV, are In both cases the uncertainties are the statistical and systematic uncertainties summed in quadrature. The overall quality of the fits is not good, as the high number of events in the inclusive distributions exposes the imperfections in the models used. As noted in Sect. 2, the fits would be better using a Lévy distribution with an additional free parameter. For all fits throughout Sect. 8, the statistical uncertainty is multiplied by the factor χ 2 /nd f if this factor is significantly greater than 1 [84]. This applies in relatively few cases, see e.g. 4 of the 8 entries in Table 3.

Multiplicity dependence
The dependence of R and λ on n ch is studied in multiplicity intervals that are chosen to be well populated, and comparable to those used by other LHC experiments [30,31,35,36,38].
Using the exponential fit, the extracted BEC parameters λ(m ch ) and R(m ch ) at 13 TeV are shown as a function of scaled multiplicity in Fig. 4 for p T > 100 and 500 MeV. The MB and HMT data agree well where they overlap.
For p T > 100 MeV the λ parameter decreases with increasing m ch , and the dependence is well described by the exponential function: The R parameter increases with multiplicity up to about m ch 3.08. For m ch ≤ 1.9, the m ch dependence of R is well described by: which is also the function used in heavy-ion studies [8,31]. A linear increase of R with 3 √ m ch follows naturally if the hadronization volume is proportional to the number of produced particles.
For higher multiplicities, the measured R parameter is observed to be independent of multiplicity. This saturation of the BEC R parameter at high multiplicity is similar to that observed by ATLAS at √ s = 7 TeV [28] and by CMS at √ s = 13 TeV [38]. The dependence of R for the region m ch > 3.08 is then fitted simply by a constant: For the combined MB and HMT data at √ s = 13 TeV and p T > 100 MeV the parameters of the m ch dependence of R and λ are given in the upper part of Table 2, with statistical and asymmetric systematic uncertainties summed in quadrature.
For p T > 500 MeV the variation of the BEC parameters is fitted to the same functional forms as used for p T > 100 MeV (Eqs. (5)-(7)), and the results are shown in the lower part of Table 2.
From Fig. 4a it can be seen that λ(m ch ) shows only a weak dependence on multiplicity for p T > 500 MeV. It is constant within measurement uncertainties, indicating a relatively high level of incoherence. Figure 4b,

Comparison of C 2 (Q) functions at 13 and 7 TeV
The BEC parameters λ and R were measured previously by ATLAS at 0.9 and 7 TeV [28]. The data sample at 7 TeV was large enough to explore high-multiplicity events, and a clear plateau in R was observed for n ch ≥ 50. This plateau at high multiplicity is confirmed by the present analysis at 13 TeV, as shown in Fig. 4b, c. However, a direct comparison of the BEC parameters obtained at the two centre-of-mass energies is complicated by small but significant differences between the Monte Carlo generators used to construct C MC 2 (Q), and hence R 2 (Q), in the two analyses (see Eqs. (1) and (2)). The BEC parameters obtained at 7 TeV are systematically lower than those at 13 TeV by about 1.5 standard deviations of the estimated 7 TeV uncertainties. For the previous 7 TeV analysis only limited MC modelling was available. For this paper, several MC models and tunings are used to estimate the modelling uncertainty. This finding is supported by a comparison between the quantities C data 2 (Q) at the two centre-of-mass energies. Figure 5a, b show C data 2 (Q) for MB events at 7 and 13 TeV. Figure 5a shows a representative plot for a scaled-multiplicity interval (3.09 < m ch ≤ 3.86) and Fig. 5b shows a representative plot for a k T interval (400 < k T ≤ 500 MeV). In Fig. 5 the uncertainties shown are the statistical and systematic uncertainties combined in quadrature. The system- The plots shown in Fig. 5a, b are representative of all multiplicity and k T intervals. The data from the two centre-of-mass energies are in agreement everywhere within the uncertainties, and the ratio of 7 TeV to 13 TeV data is 1 to better than 1% except in the very first one or two Q bins.

Dependence on the pair average transverse momentum
The BEC are studied in k T intervals from 0.1 to 1.5 GeV that are chosen so that they are populated roughly equally.
The fitted dependences of λ(k T ) and R(k T ) for the two p T thresholds are shown in Fig. 6 for p T > 100 and 500 MeV.
In the small-Q, BEC-sensitive region, the minimum k T is highly correlated with the p T threshold, and so k T thresholds of 100 and 500 MeV are applied for the p T thresholds of 100 and 500 MeV, respectively. The values of λ(k T ) and R(k T ) decrease with increasing k T . The decrease of λ(k T ) with k T is described by the function:     The R(k T ) parameter shows similar behaviour: The fitted parameters in Eqs. (8) and (9) are given in Table 3. 8.5 Double-differential dependence on m ch and k T There is sufficient data at 13 TeV to study the doubledifferential dependence: λ(m ch , k T ) and R(m ch , k T ). Figure 7 shows the two-dimensional (2D) dependence of λ and R on m ch and k T , for p T > 100 MeV. In order to convey a clear picture of the 2D dependence, the values are presented without uncertainties, but they are generally less than 10%.
The following features can be seen in Fig. 7: • The parameter λ decreases monotonically with k T in all m ch intervals, similar to the behaviour seen in Fig. 6a. • The variation of λ with m ch is rather flat for k T 500 MeV, which may be compared with the clear decrease of λ as a function of m ch seen in Fig. 4a for p T > 100 MeV. This is discussed below. • The variation of λ with m ch for k T > 500 MeV indicates a slight rise with m ch similar to what is seen in Fig. 4a for p T > 500 MeV. • The parameter R also decreases monotonically with k T in all m ch intervals, with the decrease being more pronounced in lower m ch intervals. A similar decrease is seen in Fig. 6b. • R increases with m ch in all k T intervals, starting from a lower value at larger k T .  Fig. 7 The two-dimensional dependence on m ch and k T for p T > 100 MeV for a the correlation strength, λ, and b the source radius, R, obtained from the exponential fit to the R 2 (Q) correlation functions using the MB sample for m ch ≤ 3.08 and the HMT sample for m ch > 3.08 • This increase of R flattens off to a plateau level in all k T intervals. This is similar to the plateau seen in Fig. 4b, c.
Within the precision of the present data, the onset of the plateau occurs at m ch ∼ 3 in all k T intervals.
These features can be seen more clearly in selected 1D projections in Figs. 8 and 9 for λ and R, respectively. MB and HMT data agree where they overlap. Figure 8c, d show the essentially flat behaviour of λ as a function of m ch for low k T , and the slight rise at higher k T , which are described above. As also noted above, Fig. 4a shows a slight but clear decrease of λ with m ch for p T > 100 MeV when integrated over all k T . This arises because particle p T and hence k T of the pair increase with m ch , but λ decreases as k T increases. This leads to the overall decrease with m ch seen in Fig. 4a for p T > 100 MeV.
The 2D behaviour can be parameterized according to Eqs. (8) and (9): and The parameters μ, ν, ξ and κ are shown in Figs. 10 and 11 for both the p T > 100 MeV and p T > 500 MeV selections. The 2D dependence has also been studied for the p T > 500 MeV selection and is shown in Fig. 12. As explained above, the p T selection constrains the BEC parameters to the region k T > 500 MeV.
The main features are similar to those observed in Fig. 7. In particular, the value of R tends to a plateau at large m ch in all k T intervals, but the plateau level is systematically lower than for p T > 100 MeV, in line with what is observed in Fig. 4b, c. Selected 1D projections are shown in Figs. 13 and 14 for λ and R, respectively. As with the lower p T selection, MB and HMT data agree where they overlap.
The 2D behaviour is parameterized according to Eqs. (10) and (11). The fit parameters μ, ν, ξ and κ for p T > 500 MeV are shown in Figs. 10  In the CMS publication comparison is made to the ATLAS 7 TeV results. To do so, CMS has adjusted its results for R to match the ATLAS kinematic region, p T > 100 MeV, |η| < 2.5, by extrapolating the n ch values to the lower p T threshold and excluding the same Q regions as ATLAS in the fit to R 2 (Q). These adjusted CMS values are shown together with the ATLAS results in Fig. 15. The behaviour of R in the two experiments is qualitatively similar. However, there is a clear difference in the results, particularly at low n ch . Results from this analysis (Fig. 4b) and earlier studies for the lower-energy ATLAS analysis [28] indicate a direct dependence of R on the p T threshold. So the effect of the difference in kinematic coverage as well as the differing approaches to constructing the reference sample requires further investigation.

Summary and conclusions
Results are presented of measurements of two-particle BEC of like-charge hadrons with track p T -thresholds of 100 and 500 MeV and |η| < 2.5 produced in pp collisions at √ s = 13 TeV recorded by the ATLAS detector at the CERN LHC.
The study is carried out using data collected with the MB and HMT triggers. The integrated luminosities are 151 µb −1 and 8.4 nb −1 for the MB and HMT data samples, respectively. The studies are performed using the R 2 (Q) correlation function, which is the ratio of the C data 2 (Q) correlation function obtained from the data to the C MC 2 (Q) calculated using MC events with no BEC effect. The C 2 function is itself a ratio of like-charge particle pairs to a reference sample, and the reference sample is constructed using unlike-charge particle pairs. A clear indication of BEC is observed in the region of small four-momentum difference, Q 0.2 GeV. The BEC parameters, R and λ, are studied as functions of a scaled charged-particle multiplicity m ch , the pair average transverse momentum (k T ) and in (m ch , k T )-intervals.
The parameter R(m ch ) is found to increase as α 3 √ m ch for low multiplicity for both p T thresholds. For p T > 100 MeV, α = 2.54 +0.12 −0.22 (total) fm, for scaled multiplicity up to m ch ≈ 2. For m ch ≥ 3, the source radius R saturates at a value R = 3.35 +0.20 −0.09 (total) fm, confirming the previous observation of a high-multiplicity plateau by ATLAS at √ s = 7 TeV. For p T > 500 MeV, the behaviour of R is similar, but systematically lower.
The parameter λ(m ch ) decreases with multiplicity for the lower p T threshold and is lower for the higher p T threshold but increases slightly with multiplicity. The parameters R(k T ) and λ(k T ) both decrease with increasing pair average transverse momentum k T . Considering the combined (m ch , k T ) dependence, λ decreases with k T in all m ch inter-vals, for both p T thresholds. The decrease is more pronounced at lower multiplicity. The radius parameter R also decreases with k T for p T > 100 MeV, decreasing more strongly in the lower m ch intervals. For p T > 500 MeV the behaviour of R is rather flat with k T . As a function of m ch , R increases and reaches a plateau at large multiplicity in all m ch intervals and for both p T thresholds. Within the uncertainties of this analysis, the turn-on of the plateau occurs at the same value of m ch in all k T intervals and for both p T thresholds. A comparison of the R parameter as a function of n ch with CMS results at 13 TeV shows qualitatively similar behaviour and the differences at low n ch require further investigation.
The measurements presented here complement the ATLAS measurements at √ s = 0.9 and 7 TeV by extending the studies to higher multiplicities and pair average transverse momenta, and provide a more detailed study of the BEC parameters by exploring the double-differential (m ch , k T )dependence. The present results confirm the saturation of the source radius at higher multiplicities with n ch 100. and  [85].

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All ATLAS scientific output is published in journals, and preliminary results are made available in Conference Notes. All are openly available, without restriction on use by external parties beyond copyright law and the standard conditions agreed by CERN. Data associated with journal publications are also made available: tables and data from plots (e.g. cross section values, likelihood profiles, selection efficiencies, cross section limits, ...) are stored in appropriate repositories such as HEPDATA (http:// hepdata.cedar.ac.uk/). ATLAS also strives to make additional material related to the paper available that allows a reinterpretation of the data in the context of new theoretical models. For example, an extended encapsulation of the analysis is often provided for measurements in the framework of RIVET (http://rivet.hepforge.org/).] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.

Appendix A: Effect of resonances on R 2 ( Q) correlation function
The potential influence of resonances on the R 2 (Q) is evident by studying the Q distributions of the UCP pairs (h ± h ∓ ) where both particles come from the same species of resonance; this is shown in Fig. 16. The presented distributions are for track pairs with p T > 100 MeV. The reflections of the resonance peaks can be seen clearly, and arise mainly from  Fig. 16 The inclusive Q distribution of unlike-charge particle pairs for track pairs originating from the decays of various species of single resonances. The full range of multiplicities contribute to each distribution contributions from the η, ω, ρ and K * 0 resonances, with much smaller contributions from the η and φ mesons. The contributions of f 0 (980) → π + π − and f 2 (1270) → π + π − are insignificant.