Two-particle Bose–Einstein correlations in pp collisions at \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} =}$$\end{document}s= 0.9 and 7 TeV measured with the ATLAS detector

The paper presents studies of Bose–Einstein Correlations (BEC) for pairs of like-sign charged particles measured in the kinematic range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_\mathrm{T}>$$\end{document}pT> 100 MeV and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|\eta |<$$\end{document}|η|< 2.5 in proton collisions at centre-of-mass energies of 0.9 and 7 TeV with the ATLAS detector at the CERN Large Hadron Collider. The integrated luminosities are approximately 7 \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\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}-1, 190 \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\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}-1 and 12.4 nb\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}-1 for 0.9 TeV, 7 TeV minimum-bias and 7 TeV high-multiplicity data samples, respectively. The multiplicity dependence of the BEC parameters characterizing the correlation strength and the correlation source size are investigated for charged-particle multiplicities of up to 240. A saturation effect in the multiplicity dependence of the correlation source size parameter is observed using the high-multiplicity 7 TeV data sample. The dependence of the BEC parameters on the average transverse momentum of the particle pair is also investigated.


Introduction
Particle correlations play an important role in the understanding of multiparticle production.Correlations between identical bosons, called Bose-Einstein correlations (BEC), are a well-known phenomenon in high-energy and nuclear physics (for reviews see [1][2][3][4][5][6][7][8][9]).The BEC are often considered to be the analogue of the Hanbury-Brown and Twiss effect [10] in astronomy, describing the interference of incoherently emitted identical bosons [11][12][13].They represent a sensitive probe of the space-time geometry of the hadronization region and allow the determination of the size and the shape of the source from which particles are emitted.
The production of identical bosons that are close together in phase space is enhanced by the presence of BEC.The first observation of BEC effects in identically charged pions produced in pp collisions was reported in Ref. [14].Since then, BEC have been studied for systems of two or more identical bosons produced in various types of collisions, from leptonic to hadronic and nuclear collisions (see Refs. [1][2][3][4][5][6][7] and references therein).
Studies of the dependence of BEC on particle multiplicity and transverse momentum are of special interest.They help to understand the multiparticle production mechanism.The size of the source emitting the correlated particles has been observed to increase with particle multiplicity.This can be understood as arising from the increase in the initial geometrical region of overlap of the colliding objects [15]: a large overlap implies a large multiplicity.While this dependence is natural in nucleusnucleus collisions, the increase of size with multiplicity has also been observed in hadronic and leptonic interactions.In the latter, it is understood as a result of superposition of many sources [6,[16][17][18] or related to the number of jets [19,20].High-multiplicity data in proton-proton interactions can serve as a reference for studies of nucleusnucleus collisions.The effect is reproduced in both the hydrodynamical/hydrokinetic [21][22][23] and Pomeron-based [24,25] approaches for hadronic interactions where highmultiplicities play a crucial role.The dependence on the transverse momentum of the emitter particle pair is another important feature of the BEC effect [26].In nucleusnucleus collisions the dependence of the particle emitter size on the transverse momentum is explained as a "collective flow", which generates a characteristic fall-off of the emitter size with increasing transverse momentum [27,28] while strong space-time momentum-energy correlations may offer an explanation in more "elementary" leptonic and hadronic systems [4,5,7,[21][22][23]26] where BEC measurements serve as a test of different models [21][22][23][29][30][31][32][33][34][35][36].
In the present analysis, studies of one-dimensional BEC effects in pp collisions at centre-of-mass energies of 0.9 and 7 TeV, using the ATLAS detector [37] at the Large Hadron Collider (LHC), are presented.At the LHC, BEC have been studied by the CMS [38,39] and ALICE [40,41] experiments.In the analysis reported here, the studies are extended to the region of high-multiplicities available thanks to the high-multiplicity track trigger.The results are compared to measurements at the same or lower energies.

Two-particle correlation function
Bose-Einstein correlations are measured in terms of a twoparticle correlation function, where p 1 and p 2 are the four-momenta of two identical bosons in the event, ρ is the two-particle density function, and ρ 0 is a two-particle density function (known as the reference function) specially constructed to exclude BEC effects.
In order to compare with data over the widest possible range of centre-of-mass energies and system sizes, the density function is parameterized in terms of the Lorentzinvariant four-momentum difference squared, Q 2 , of the two particles, The BEC effect is usually described by a function with two parameters: the effective radius R and the strength parameter λ [42], where the latter is also called the incoherence or chaoticity parameter.A typical functional form is For fully coherent emission of identical bosons, λ = 0, while for incoherent (chaotic) emission, λ = 1.The QR dependence comes from the Fourier transform of the distribution of the space-time points of boson emission.Several different functional forms have been proposed for Ω(λ, QR).Those used in this paper are described in Sect.2.4.The fitted parameter ε takes into account long-distance correlations not fully removed from ρ 0 .Finally, C 0 is a normalization constant, typically chosen such that C 2 (Q) is unity for large Q.In this paper, the density function ρ is calculated for like-sign charged-particle pairs, with both the ++ and −− combinations included, ρ(Q) ≡ ρ(++, −−).

Coulomb correction
The long-range Coulomb force causes a momentum shift between the like-sign and unlike-sign pairs of particles.The density distributions are corrected for this effect using the Gamow penetration factor [43][44][45], where the Gamow factor G(Q) is given by with the dimensionless parameter ζ defined as Here α is the electromagnetic fine-structure constant and m is the particle mass.In the following, the pion mass is considered.The sign of ζ is positive for like-sign pairs and negative for unlike-sign pairs.The size of this correction is found not to exceed 20% for Q < 0.03 GeV.Neither the Coulomb interaction nor the BEC effect are present in the generation of MC event samples which are used in the analysis.The Coulomb correction is thus not applied to MC events.

Reference sample
A good choice of the reference sample is important to allow the experimental detection of the BEC signal.Ideally, ρ 0 (Q) should include all relevant momentum correlations, from both the local and global conservation laws.Thus, several different choices have been studied to construct an appropriate reference sample.Most of the proposed approaches use random pairing of particles, such as mixing particles from different events (the "mixed event" technique [46]), or choosing them from the same event but from opposite hemispheres or by rotating the transverse momentum vector of one of the particles of the like-sign pair [7].Although these mixing techniques reproduce the topology and some properties of the event under consideration and destroy BEC, they violate energy-momentum conservation.Moreover, there are many possible ways to construct the pairs, such as mixing the particles randomly, or keeping some topological constraints such as the event multiplicity, the invariant mass of the pair or the rapidity of the pair.All of these introduce additional biases in the BEC observables.
A natural choice is to use the unlike-sign particle pairs from the same events that are used to form pairs of likesign particles, i.e., ρ 0 (Q) ≡ ρ(+−), called in the following the unlike-charge reference sample.This sample has the same topology and global properties as the like sign sample ρ(++, −−), but is naturally free of any BEC effect.However, this sample contains hadron-pairs from the decay of resonances such as ρ, η, η ′ , ω, φ, K * , which are not present in the like-sign combinations.These contribute to the low-Q region and can give a spurious BEC signature with a large effective radius of the source [47][48][49][50][51].
In this paper, the unlike-charge reference sample is used.To account for the effects of resonances, the twoparticle correlation function C 2 (Q) is corrected using Monte Carlo simulation without BEC effects via a double-ratio R 2 (Q) defined as

The parameterizations of BEC
Various parameterizations of the Ω(λ, QR) function can be found in the literature, each assuming a different shape for the particle-emitting source.In the studies presented here, the data are analysed using the following parameterizations: the Goldhaber parameterization [14] of a static Gaussian source in the plane-wave approach, which assumes a spherical shape with a radial Gaussian distribution of the emitter; -the exponential parameterization of a static source which assumes a radial Lorentzian distribution of the source.This parameterization provides a better description of the data at small Q values, as discussed in [7].
The first moment of the Ω(QR) distribution corresponds to 1/R for the exponential form and to 1/(R √ π) for the Gaussian form.To compare the values of the radii obtained from the two functions, the R value of the Gaussian should be compared to R/ √ π of the exponential form.
3 Experimental details

The ATLAS detector
The ATLAS detector [37] is a multi-purpose particle physics experiment operating at one of the beam interaction points of the LHC.The detector 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 detector (ID), which has full coverage in φ and covers the pseudorapidity range |η| < 2.5. 1 It consists of a silicon pixel detector (Pixel), a silicon microstrip detector (SCT) and a transition radiation tracker (TRT).These detectors are immersed in a 2 T solenoidal magnetic field.The Pixel, SCT, and TRT detectors have typical position resolutions of 10, 17 and 130 µm for the r-φ coordinate, respectively.In the case of the Pixel and SCT, the resolutions are 115 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 θ as η = − ln tan(θ/2).and 580 µm, respectively, for the second measured coordinate.A track from a charged particle traversing the full radial extent of the ID would typically have three Pixel hits, eight or more SCT hits and more than 30 TRT hits.
The ATLAS detector has a three-level trigger system: Level 1 (L1), Level 2 (L2) and Event Filter (EF).For this measurement, the trigger relies on the L1 signals from the Beam Pickup Timing devices (BPTX) and the Minimum-Bias Trigger Scintillators (MBTS).The BPTX are composed of electrostatic button pick-up detectors attached to the beam pipe and located 175 m from the centre of the ATLAS detector in both directions along the beam pipe.The MBTS are mounted at each end of the detector in front of the liquid-argon end-cap calorimeter cryostats at z = ±3.56 m.They are segmented into eight sectors in azimuth and two rings in pseudorapidity (2.09 < |η| < 2.82 and 2.82 < |η| < 3.84).The MBTS trigger used for this paper is configured to require one hit above threshold from either side of the detector, referred to as a single-arm trigger.The efficiency of this trigger was studied with events collected with a separate prescaled L1 BPTX trigger, filtered by ID requirements at L2 and at EF level in order to obtain inelastic interactions [52].
High-multiplicity track (HM) events were collected at 7 TeV using a dedicated high-multiplicity track trigger.At L1, the collisions were triggered using the summed transverse energy (ΣE T ) in all calorimeters, calibrated at the electromagnetic energy scale [53].The minimum-bias (MB) events were required to have ΣE T > 20 GeV.A high number of hits in the SCT was required at L2, while at the EF level at least 124 tracks with p T > 400 MeV were required to originate from a single vertex.

Data and Monte Carlo samples
The study is carried out using the pp-collision datasets at the centre-of-mass energies √ s = 0.9 and 7 TeV that were used in previously published ATLAS studies of minimumbias interactions [52,54].
The event and track selection criteria are the same as the ones used for the ATLAS minimum-bias multiplicity analysis [54] with the same minimum-bias trigger and quality criteria for the track reconstruction.All events in these datasets are required to have at least one vertex [64], formed from a minimum of two tracks with p T > 100 MeV and consistent with the average beam spot position within the ATLAS detector (primary vertex) [65].The tracks satisfying the above-mentioned selection criteria are used as the input to determine the corrected distributions, as described in Sect.3.3.The multiplicity of selected tracks with p T > 100 MeV and |η| < 2.5 within an event is denoted by n sel .
The contributions from beam-gas collision and from non-collision background (cosmic rays and detector noise) were investigated in Ref. [52] and found to be negligible.Events with more than one primary vertex (less than 0.3% of the sample) are rejected in order to prevent a bias from multiple proton-proton interactions (pile-up) in the colliding proton bunches.
The same event selection criteria are applied to highmultiplicity events, which are defined to be those with at least 120 selected tracks.To estimate the possible influence of multiple pp interactions in the 7 TeV highmultiplicity track trigger data, the distribution of the distances ∆z between the z coordinates of primary and pileup vertices are studied.The study shows that fewer than one pile-up track out of the selected tracks contributes to the HM trigger sample, which is negligible.
For the measurements at √ s = 0.9 TeV, about 3.6×10 5 events with a total of more than 4.5 × 10 6 tracks are after selection, and in the case of √ s = 7 TeV, about 10 7 events with about 2.1×10 8 tracks overall are after selection.This corresponds to integrated luminosities of ∼7 µb −1 and ∼190 µb −1 at 0.9 TeV and 7 TeV, respectively.For the measurements at 7 TeV with the high-multiplicity track trigger, about 1.8 × 10 4 events with more than 2.7 × 10 6 tracks overall were after selection.This corresponds to integrated luminosity of ∼12.4 nb −1 .Large Monte Carlo samples of minimum-bias events were generated using the PYTHIA 6.421 Monte Carlo event generator [55] with the ATLAS MC09 set of optimised parameters (tune) [56] with non-diffractive, singlediffractive and double-diffractive processes.As discussed in Sec.2.2, no simulation of the BEC effect is implemented in the generator.This is the baseline Monte Carlo generator which reproduces single-particle spectra [52].The generated events were passed through the ATLAS simulation and reconstruction chain; the detector simulation program [58] is based on GEANT4 [59].Dedicated sets of high-multiplicity events were also generated.
For the study of systematic effects, additional Monte Carlo samples were produced using the PHOJET 1.12.1.35generator [60], PYTHIA with the Perugia0 tune [61]; and the EPOS 1.99 v2965 generator [36] for the high-multiplicity analysis.The PHOJET program uses the Dual Parton Model [62] for low-p T physics and is interfaced to PYTHIA for the fragmentation of partons.The EPOS generator is based on an implementation of the QCDinspired Gribov-Regge field theory describing soft and hard scattering simultaneously, and relies on the same parton distribution functions as used in PYTHIA.The EPOS LHC tune is used with parameters optimised to describe the LHC minimum-bias data [63].

Data correction procedure
Following the procedure applied in the previous ATLAS minimum-bias measurements [52,54], each track is assigned a weight which corrects for the track reconstruction efficiency, for the fraction of secondary particles, for the fraction of the primary particles2 outside the kinematic range and for the fraction of fake tracks. 3In addition, the effect of events lost due to trigger and vertex reconstruction inefficiencies is corrected for using an event-by-event weight applied to pairs of particles in the Q distribution.The efficiency of the high-multiplicity track trigger has been studied in data as a function of the number of reconstructed tracks and is found to be 5% for 120 selected tracks and to reach a plateau at 100% once 150 tracks are selected.The measured trigger inefficiency is used to correct the experimental distributions and is found to have negligible impact on the extraction of the BEC parameters discussed in Sec. 5.
The multiplicity distributions are corrected to the particle level using an iterative method that follows the Bayesian approach [66] as it is described in Refs.[52,54].An unfolding matrix reflecting the probability of reconstructing n sel charged tracks in an event with generated chargedparticle multiplicity n ch is populated using Monte Carlo simulation and applied to the data.The unfolding matrix is built using the ATLAS MC09 PYTHIA tune [56].The unfolding procedure converges after the fifth iteration.It is found that the corrected multiplicity distribution agrees well with the published result [52,54].The unfolding procedure of the 7 TeV high-multiplicity data follows the same technique and unfolding matrix used in the previous analysis of minimum-bias data in Ref. [52], restricted to the region of high charged particle multiplicity specific to this analysis, and convolved with a normalised Gaussian distribution to account for the experimental resolution on the number of selected tracks.It is found that a number of 120 selected tracks at detector level, n sel , corresponds to about 150 charged particle, n ch , at particle level.Momentum distributions are unfolded in a similar way.
For all distributions, closure tests are carried out using Monte Carlo samples corrected according to the same procedure as used in the data.The difference obtained between the reweighted distributions and those at the particle level is due to tracking effects such as a smaller reconstruction efficiency for pairs of tracks with very small opening angle.These effects are small for correlation functions constructed using data, typically 1-3%, and are included in the systematic uncertainty.In the case of the unfolded Q distributions, the data are corrected for the bias from secondary tracks using Monte Carlo simulation and the corresponding systematic uncertainty is obtained by variation of the amount of material in the inner detector by ±10%.

Systematic uncertainties
The systematic uncertainties of the inclusive fit parameters, R and λ, of the exponential model are summarized in Table 1 .The following contributions to the systematic uncertainties on the fitted parameters are considered.
The systematic uncertainties resulting from the track reconstruction efficiency, which are parameterized in bins of p T and η, were determined in earlier analyses [52,54].These cause uncertainties in the track weights of particle pairs in the Q distributions entering the correlation functions.
The effects of track splitting and merging are sizeable only for very low Q values (smaller than 5 MeV), and are found to be negligible for the measurements with Q ≥ 20 MeV.
The leading source of systematic uncertainty is due to differences in the Monte Carlo generators used to calculate the R 2 correlation function from the C 2 correlation function.The corresponding contribution to the systematic uncertainty is estimated as the root-mean-squared (RMS) spread of the results obtained for the different Monte Carlo datasets.The statistical uncertainties arising from the Monte Carlo datasets are negligibly small.
The systematic uncertainty due to Coulomb corrections is estimated by varying the corrections by ±20%.
The influence of the fit range is estimated by changing the upper bound of the Q range from the nominal 2 GeV: decreasing it to 1.5 GeV and increasing it up to 2.5 GeV.The latter better estimates the uncertainty due the longrange correlations.This contribution is taken into account by the value of ε, the parameter in the linear term of Eq. ( 3) describing the long-range correlations.
Other effects contributing to the systematic uncertainties are the lowest value of Q for the fit, the bin size and exclusion of the interval 0.5 ≤ Q ≤ 0.9 GeV due to the overestimate of the ρ meson contribution in the Monte Carlo simulations, as discussed in the following Sec.5.1.These uncertainties are estimated by varying the lowest Q value in the fit by ±10 MeV, by changing the bin size by ±10 MeV, and by broadening the excluded interval by 100 MeV on both sides.
The background of photon conversions into e + e − pairs was studied and found to be negligible.
To test the effect of treating all charged particles as pions, the double-ratio correlation functions R 2 are also obtained using only identical particles in the Monte Carlo sample to compute the correction.The resulting BEC parameters fitted to the R 2 functions defined this way show negligible differences to the nominal result and no further systematic uncertainties are assigned.
Finally, the systematic uncertainties are combined by adding them in quadrature and the resulting values are given in the bottom row of Table 1.
The same sources of uncertainty are considered for the differential measurements in n ch and the average transverse momentum k T of a pair, and their impact on the fit parameters is found to be similar in size.

Two-particle correlations
In Fig. 1 the double-ratio R 2 (Q) distributions, measured for 0.9 and 7 TeV, are compared with Gaussian and exponential fitting functions, Eqs. ( 8) and (9).The fits are performed in the Q range 0.02 GeV to 2 GeV and with a bin width of 0.02 GeV.The upper Q limit is chosen to be far away from the low-Q region, which is sensitive to BEC effects and resonances.Around Q ∼ 0.7 GeV there is a visible bump which is due to an overestimate of ρ → π + π − decays in the Monte Carlo simulation.Therefore the region 0.5 ≤ Q ≤ 0.9 GeV is excluded from the fits.As seen in Fig. 1, the Gaussian function does not describe the low-Q region while the exponential function provides a good description of the data.
The resolution of the Q variable is better than 10 MeV for the region most sensitive to BEC effect, Q < 0.4 GeV.The Q resolution is included in the fit of R 2 by convolving the fitting function with a Gaussian detector resolution function.The change in the fit results from those with no convolution applied is found to be negligible.
In the process of fitting R 2 (Q) with the exponential function, large χ 2 values are observed, in particular for the 7 TeV sample where statistical uncertainties on the fitted data points are below 2-4%.These large χ 2 values can be traced back to a small number of individual points or small cluster of points.The removal of these points does not change the results of the fit while the χ 2 substantially improves.In the analysis of the 7 TeV data, for most of the considered cases, the expected statistical uncertainties are small compared to the systematic ones, therefore only total uncertainties on the fitted parameters are given.The latter include the statistical uncertainties rescaled by χ 2 /ndf [57].For consistency, the same treatment is applied to the 0.9 TeV analysis where the statistical uncertainties are of the same order of magnitude as the systematic ones.
The results of BEC parameters for fits of the twoparticle double-ratio correlation function R 2 (Q) for events with the unlike-charge reference sample are The values of the fitted parameters are close to the values obtained by the CMS [39] and ALICE [40] experiments.

Multiplicity dependence
The R 2 (Q) functions defined in Eq. ( 7), are shown for various multiplicity intervals in Fig. 2 for 0.9 TeV, 7 TeV and 7 TeV high-multiplicity data.The multiplicity intervals are chosen so as to be similarly populated and comparable to those used by other LHC experiments [38][39][40][41].Only the exponential fit is shown.As in the fit procedure for the inclusive case, the detector Q resolution is included in the fits.
Within the multiplicity studies, the BEC parameters are also measured by excluding the low-multiplicity events, n ch < 8, expected to be contaminated by diffractive physics [52].No noticeable changes in the strength parameters and radii for n ch ≥ 8 are observed compared to the full multiplicity range for n ch ≥ 2.
The multiplicity dependence of the λ and R parameters is shown in Fig. 3.
The λ parameter decreases with multiplicity, faster for 0.9 TeV than for 7 TeV interactions.The decrease of the λ parameter with n ch is found to be well fitted with the exponential function λ(n ch ) = p 0 e −p1n ch .The fit parameter values are presented in Table 2 for 0.9 TeV and for the combined nominal and high-multiplicity 7 TeV data.
The R parameter increases with multiplicity up to about n ch ≃ 50 independently of the center of mass energy.For higher multiplicities, the measured R parameter is observed to be independent of multiplicity.For n ch ≤ 82 at 0.9 TeV and n ch < 55 at 7 TeV the n ch dependence of R is fitted with the function R(n ch ) = p 0 3 √ n ch , similar to that used in heavy-ion studies [41,68].The results of the fit are presented in Table 2 and are close to the CMS results [39].The n ch dependence of R at 7 TeV is fitted with a constant R(n ch ) = p 0 for n ch > 55; the resulting value is given in Table 2.
In Fig. 3, the multiplicity dependence of the R is compared with the results of CMS at the same energies and with the results of UA1 in pp collisions at √ s = 200 GeV.The radii measured by CMS are obtained from the exponential fit to the double-ratio correlation functions with a reference sample of mixed events [39].The radius values obtained by UA1 are from the Gaussian fit to the single-ratio with a mixed-event reference sample of the same multiplicity range and are multiplied by √ π.The data from the three experiments are in good agreement up to n ch ≃ 50.
The observed change of the fitted parameters with multiplicity has been predicted in Refs.[7,[16][17][18], and is similar to the one also observed in e + e − interactions [19], however the saturation of R for very high multiplicity is observed for the first time.
The saturation of R at high multiplicities is expected in a Pomeron-based model [24,25] as the consequence of the overlap of colliding protons, with the value of the radius at n ch ≈ 70 close to the one obtained in the present studies.However, the same model predicts that above n ch ≈ 70, R will decrease with multiplicity, returning to its lowmultiplicity value which is not supported by the data.

Dependence on the transverse momentum of the particle pair
The average transverse momentum k T of a particle pair is defined as half of the magnitude of the vector sum of the two transverse momenta, k T = |p T,1 + p T,2 |/2.The study is performed in the k T intervals which are chosen in a way to be similarly populated and, as for the multiplicity bins, to be similar to the intervals used by other LHC experiments [38][39][40][41].
As an example, the R 2 (Q) distributions for the 500 ≤ k T ≤ 600 MeV interval for the 0.9 TeV, 7 TeV and highmultiplicity 7 TeV samples are shown in Fig. 4 together with the results of the corresponding exponential fit.For the R 2 (Q) correlation function measured at 7 TeV (see Fig. 4(b)), there is an indication that the Monte Carlo simulation overestimates the production and decay of the ω-meson in the Q region of 0.3-0.44GeV.This region is thus excluded from the fit range for k T > 500 MeV bin results.
In the region most important for the BEC parameters, the quality of the exponential fit is found to deteriorate as k T increases.This is due to the fact that at large k T values, the characteristic BEC peak becomes steeper than the exponential function can accommodate.Despite the deteriorating fit quality, the behaviour of the fitted parameters is presented for comparison with previous experiments.
The fit values of the λ and R parameters are shown in Fig. 5 as a function of k T .The values of both λ and R decrease with increasing k T .
The decrease of λ with k T is well described by an exponential function, λ(k T ) = p 0 e −p1kT .The k T dependence of the R parameter is also found to follow an exponential decrease, R(k T ) = p 0 e −p1kT .The shapes of the k T dependence are similar for the 7 TeV and the 7 TeV highmultiplicity data.The results of the fits are presented in Table 2.
In Fig. 5(b), the k T dependence of the R parameter is compared to the measurements performed by the E735 [69] and the STAR [70] experiments with mixed-event reference samples.These earlier results were obtained from Gaussian fits to the single-ratio correlation functions and therefore the values of the measured radii are multiplied by √ π as discussed in Sect.2.4.The values of the parameters are observed to be energy-independent within the uncertainties.
In Fig. 6, the k T dependence of λ and R, obtained for the 7 TeV data, is also studied in various multiplicity regions: 2 ≤ n ch ≤ 9; 10 ≤ n ch ≤ 24; 25 ≤ n ch ≤ 80; and 81 ≤ n ch ≤ 125.The decrease of λ with k T is nearly independent of multiplicity for n ch > 9 and the same as for the inclusive case.For n ch ≤ 9 no conclusions can be drawn due to the large uncertainties.The R-parameter decreases with k T and exhibits an increase with increasing multiplicity as was observed for the fully inclusive case.√ s =0.9 and 7 TeV, compared to the equivalent measurements of the CMS [38,39] and UA1 [67] experiments.The solid and dashed curves are the results of (a) the exponential and (b) 3  √ n ch for n ch < 55 fits.
The dotted line in (b) is a result of a constant fit to minimum-bias and high-multiplicity events data at 7 TeVfor n ch ≥ 55.The error bars represent the quadratic sum of the statistical and systematic uncertainties.

Summary and conclusions
The two-particle Bose-Einstein correlations of like-sign hadrons with p T > 100 MeV and |η| < 2.5 produced in pp collisions recorded by the ATLAS detector at 0.9 and 7 TeV at the CERN Large Hadron Collider are studied.In addition to minimum-bias data, high-multiplicity data recorded at 7 TeV using a dedicated trigger are investigated.The integrated luminosities are about 7 µb −1 , 190 µb −1 and 12.4 nb −1 for 0.9 TeV, 7 TeV minimum-bias and 7 TeV high-multiplicity data samples, respectively.The studies were performed using the double-ratio correlation function.In the double-ratio method, the singleratio correlation function obtained from the data is divided by a similar single-ratio calculated using Monte Carlo events, which do not have BEC effects.The reference sample for each of the two single-ratios is constructed from unlike-sign charged-particle pairs.
A clear signal of Bose-Einstein correlations is observed in the region of small four-momentum difference.To quantitatively characterize the BEC effect, Gaussian and exponential parametrizations are fit to the measured correlation functions.As observed in studies performed by other experiments, the Gaussian parameterization provides a poor description of the BEC-enhanced region and hence the exponential parameterization is used for the final results.
The BEC parameters are studied as a function of the charged-particle multiplicity and the transverse momentum of the particle pair.A decrease of the correlation strength λ along with an increase of the correlation source size R are found with increasing charged-particle multiplicity.On the other hand no dependence of R on the centre-of-mass energy of pp collisions is observed.For the first time a saturation of the source size parameter is observed for multiplicity n ch ≥ 55.The correlation strength λ and the source size R are found to decrease with increasing average transverse momentum of a pair.The study of BEC in (n ch , k T ) bins at 7 TeV shows a decrease of the R parameter with k T for different multiplicity ranges, while the R values increase with multiplicity.The λ parameter is found to decrease with k T independently of the multiplicity range.These resemble the dependences for the inclusive case at 7 TeV for minimum-bias and high-multiplicity data.
A comparison is made to the measurements by other experiments at the same and lower energies where possible.The measurements presented here complement the earlier measurements by extending the studies to higher multiplicities and transverse momenta.This has allowed a first observation of a saturation in the magnitude of the source radius at high charged-particle multiplicities, and confirms the exponential decrease, observed in previous measurements of the radii with increasing pair transverse momenta.

Fig. 1 .
Fig. 1.The two-particle double-ratio correlation function R2(Q) for charged particles in pp collisions at (a) √ s =0.9 TeV, (b) 7 TeV and (c) 7 TeV high-multiplicity events.The lines show the Gaussian and exponential fits as described in the legend.The region excluded from the fits is indicated.The error bars represent the statistical uncertainties.

Fig. 2 .
The two-particle double-ratio correlation function R2(Q) for charged particles in pp collisions for multiplicity intervals (a) 36 ≤ n ch < 45 at √ s =0.9 TeV, (b) 68 ≤ n ch < 79 at 7 TeV and (c) 183 ≤ n ch < 197 at 7 TeV high-multiplicity events.The lines show the results of the exponential fit.The region excluded from the fits is indicated.The error bars represent the statistical uncertainties.

λFig. 3 .
Fig.3.Multiplicity, n ch , dependence of the parameters (a) λ and (b) R obtained from the exponential fit to the two-particle double-ratio correlation functions R2(Q) at √ s =0.9 and 7 TeV, compared to the equivalent measurements of the CMS[38,39]

Fig. 4 .Fig. 5 .Fig. 6 .
Fig. 4. The two-particle double-ratio correlation function R2(Q) for charged particles in pp collisions for 500 ≤ kT < 600 MeV interval at (a) √ s =0.9 TeV, (b) 7 TeV and (c) 7 TeV high-multiplicity events.The average transverse momentum kT of the particle pairs is defined as kT = |pT,1 + pT,2|/2.The lines show the exponential fits.The region excluded from the fits is indicated.The error bars represent the statistical uncertainties.

Table 1 .
Systematic uncertainties on λ and R for the exponential fit of the two-particle double-ratio correlation function R2(Q) in the full kinematic region at √ s =0.9 and 7 TeV for minimum-bias and high-multiplicity (HM) events.

Table 2 .
Results of fitting the multiplicity, n ch , and the transverse momentum of the pair, kT, dependence of the BEC parameters R and λ with different functional forms and for different data samples.The error represent the quadratic sum of the statistical and systematic uncertainties.