arxiv: v1 [nucl-ex] 3 Jan 2019

Transcription

1 EPJ Web of Conferences will be set by publisher DOI: will be set by publisher c Owned by authors, published by EDP Sciences, 9 arxiv:9.7v [nucl-ex] Jan 9 Correlations,, ridge in pp p-pb collisions Alice Ohlson,a Collaboration Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany Introduction Abstract. Measurements made by Collaboration of single- two-particle in high-energy pp p Pb collisions are used to characterize interactions in small collision systems, tune models of particle production in QCD, serve as a baseline heavy-ion observables. measurements of charged-particle density, dn ch /dη, are shown in pp p Pb collisions, including data from top center-of-mass energy achieved at Large Hadron Collider (LHC), s = TeV. Two-particle angular correlations in p Pb collisions are studied in detail to investigate long-range correlations in which are reminiscent of structures previously thought unique to heavy-ion collisions. In high-energy hadronic collisions, studies of inclusive single-particle are used to investigate particle production in QCD. charged-particle density in pp p Pb collisions is measured by over a range of centre-of-mass energies, including top LHC energy of s = TeV. are also shown, all experimental data is compared Monte Carlo models. Since produced is dominated by soft (low-momentum) particle production, which is in non-perturbative regime of QCD, se measurements can be used to furr constrain tune models. Beyond single-particle inclusive measurements, two-particle correlation studies have yielded surprising results in small collision systems, showing presence of correlations between particles over large ranges in in high- pp p Pb collisions. se correlations are reminiscent of features observed in heavy-ion collisions where y are commonly attributed to anisotropic flow (v n ). transverse momentum (p T ), (η), particle species dependence of v in p Pb collisions has been measured in. In particular, in analysis of correlations between ward muons mid-rapidity charged hadrons it is possible to measure v large values of in both proton-going Pb-going directions. se observations will be used to deepen our understing of possible collective effects in small collision systems ir implications heavy-ion physics. a alice.ohlson@cern.ch

2 EPJ Web of Conferences detector main subsystems of detector [] used in analyses reported here are: Inner Tracking System (ITS), Time Projection Chamber (TPC), Forward Muon Spectrometer (FMS), V system. ITS, used tracking vertex reconstruction, consists of six layers of silicon detectors; innermost two layers comprise Silicon Pixel Detector (SPD). Short track segments ( tracklets ) can be reconstructed using only SPD, are used in density muon-hadron correlations analyses below. Inmation from ITS TPC can also be combined to fully reconstruct charged particle tracks. Muons are detected in FMS, which has a coverage of < η <.. composition of parent particles of detected muons depends on transverse momentum: at low p T muons predominantly come from weak decays of pions kaons, while at high p T muons are largely result of heavy flavor decays. V detectors, located at ward rapidity ( VA at.8 < η <. VC at.7 < η <.7), are used triggering also to classify overall event activity. Symmetric coverage can also be achieved by utilizing only two of four rings in each V detector, innermost two rings of VC (.7 < η <.7) outermost two rings of VA (.8 < η <.9), as was done in muon-hadron analysis below. Multiplicity density In, charged-particle density, dn ch /dη has been measured across a wide range in center-of-mass energy, at s =.9,.6,.76, 7, 8, TeV. is measured in different classes of events, including inelastic events ( ), inelastic events at least one charged particle produced in η < ( > ), non-single-diffractive events ( NSD ). Figure shows results > classes, which demonstrate power-law scaling s. Results from p Pb collisions are also shown in Fig. []. Additionally, charged particle density has been measured as a function of in NSD events at s =.9.6 TeV, in > events at TeV, as shown in Fig.. are also compared to multiple Monte Carlo (MC) models Figure. mid-rapidity charged-particle density, dn ch /dη, is shown as a function of center-ofmass energy (left) pp (right) pp, p Pb, Pb Pb collisions at LHC [, 7].

3 XLVI International Symposium on Multiparticle Dynamics 7 dn ch/dη 6 Table. Mean of measured by UA [9] in proton antiproton Cq -moments () collisions at s =.9 TeV, by at s =.9 TeV.6 TeV, NSD events in three different intervals. first error is statistical second systematic. η <. (x ) Probability P (N ch) Probability P (N ch) dn ch/dη η <. (x ) UA pp pp s =.9 TeV s =.6 TeV η <. s =.6 TeV Nch pp.6. pp.6 ±. ±..7 ±. ±. s =.9 TeV NSD±. ± C pp.9 ±..96 ±. ±.6. ±. ±. pp NSD pp Pythia NSD±. η <. (x ) η <. (x ) C. ±. ±.. ±.6 ±..76 ±.9 ± UA pp NSD UA pp Phojet pp NSD 9 ± ± 8. ±. ±.6.6 ±.6 ±. CMS pp NSD CMS pp NSD η <. s =.9 TeV s =.6 TeV - - Nch 7.8 ±.8 ± ±. ± ±.6 ± NSD - η <. (x ) NSD - C.7 ±. ±..77 ±. ±..8 ±. ±.6 Pseudorapidity η - Pseudorapidity η η <. (x ) C. ±. ±.. ±. ±..6 ±.6 ±. -. ±.9 ±Multiplicity..6 ±. ±.9. ±. ±. Multiplicityof N chdnch /dη at N chnsd Fig. 6. Left: measured dependence s =.9 TeV (full symbols) (open η <pp.data at symbols) collisions. measurements (squares) are compared to UA pp data [] (triangles) to CMS Fig.LHC 9. Exped views of measured low- region of corrected NSD events, left N 9.7 ±. ±.9.86 ±. ±. ch.6 [] (stars). Right: dependence of dn s = TeV (full symbols) ch /dη at ch TeV data. gray bs indicate systematic uncertainty. Distribution.9 TeV right.6 NSD events are not. ±..79 ±. ±.7 compared NSD (open symbols) collisions. measurement (squares) NSD collisionscis to CMS NSD.7 data± [] (stars) normalized to unity but scaled down in such a way that NSDCevents match at high multiplicities,.9 ±systematic. ±.. ±.6 ±. to model predictions, PYTHIA tune D6T [9] (solid line) PHOJET [] (dashed line). For data, which makes difference at low clearly visible. Left: data at s =.9 TeV. Right: data at s =.6 TeV. Note.9 ±. ±.6.8 ±.7 ±. uncertainties are shown as. shaded areas; statistical uncertainties are invisible (smaller than data marks). For CMS data error that η <. η < have been scaled clarity by factor indicated. Figure. hdn /dηi vs η in pp collisions is shown at s =.9 (left),.6 (center), TeV (right) [7, 8]. η <. P(z ) <N ch> Reduced moment C q Probability P (N ch) Probability P (N ch) bars show statistical systematic uncertainties added in quadrature. UA. TeV Ratio.9 /.6 TeV Ratio Data / MC Ratio Data / MC - on due to limerror inferred from χ minimization, D6T (9) D6T (9) except in lowua.9 TeV -detector acceptance efficiency ited hit statistics current alignment of region. In that region, bias.9 TeV ATLAS-CSCprecision (6) ATLAS-CSC (6)is about - %, but.6 TeV detector is estimated by this method to be. %. is negligible. Perugia- () () ree, we C statistical uncertaintyperugia- uncertainty in background corrections was estimated ac--added estimated value ofphojet bias to statistical - PHOJET C cording to description in Section. uncertainty in this region. correction procedure is in η <. total systematic uncertainty on pseudorapid- of sensitive to shape of distribution ity density measurement at.9 TeV is smaller than. % - events, which produce response matrix. - collisions is about. % NSD collisions. Table summarizes systematic uncertainties UA s =.9 s =.6 TeV distribution measurements. Note that- NSD At.6 TeV, TeV corresponding uncertainties are below C - is a function of which is re- η <. 6.7 %.7 % NSD collisions, respectively.-uncertainty C <. are dominated by uncertainties in flected ηby <. C For all cases, ηy ranges of values. Furr details about cross sections of diffractive processes analysis, corrections, of systematic 6 8 ir kinematics. evaluation Energy s (GeV). To.evaluate systematic error on.uncertainties are in [8]. distribution, a new response matrix was generated each Both density dis.. Fig... Energy dependence of Cq -moments () of change listed above used to unfold measured spectribution measurements have cross-checked measured by UAbeen [9] at by a sec. trum.. difference between se unfolded spectra. ond analysis employing Time-Projection both energies NSD events in two different Chamber unfolded spectrum produced unaltered re- (TPC) intervals. represent combined statistical []. It uses tracks vertices reconstructed in 6 8 error bars uncertainties. data at.9 TeV are sponse matrix determines systematic uncertainty. z = N ch / <N ch> Multiplicity N ch systematic Multiplicity N ch pseutpc in region η displaced <.8. visibility. Additional systematic uncertainties originate from horizontally dorapidity density is corrected using a method similar to Fig.. Comparison of measured events to models range η <..Fig.. Comparison of in KNO variunfolding method itself, consisting of two contributions. usedtunes: analysis. results [] of ables measured two by UA [8,9] in proton antiproton Predictions are shown based on PHOJET model [] (solid line) that PYTHIA D6TSPD [9] (dashed line), ATLAS-CSC collisions first one from fluctuations duebars to independent analyses are consistent. (dotted line), arises Perugia- []statistical (dash-dotted line). error data points represent shadedat s =. TeV.9 TeV, by at s =.9 TeV Figure. are shown atstatistical s =uncertainties,.9 (left).6 TeV (center) MC model number of events used to Left: produce.9response areasfinite represent systematic uncertainties. data at TeV. Right: data at.6 TeV. For both cases ratios between.6 TeV, NSD events in η <.. In lower part matrix well as model limited number of events in summarized measuredasvalues calculations are shown in lower part insame convention. represent Table. For η < shaded. areas η <. ourratio between measurements at.9 TeV.6 TeV comparisons. KNO scaling is also shown (right). [8] combined statistical systematic uncertainties. measurement. unfolding procedure was repeated are results compared to UA measurement pp colli- is shown. error bars represent combined statistical times while romizing input measurement sions at6results s =.9 TeV [9]. Note that mean multiplici- systematic uncertainties. response matrix according to ir respective statistical ties quoted in this table are those calculated from mul stronger energy seen in charged-particle high multiplicities.9 TeV sample, PHO- tiplicity are ree slightly different uncertainties. resulting uncertainty due to re- In this rise section, density density is, at least partly, due a larger fraction of highjet model agrees well data. PYTHIA tunes from values givenresults in Table. to value of two pseudosponse matrix fluctuations is negligible. uncertainty distribution are presented centre-of-mass events. D6T Perugia- underestimate data at high mul- rapidity in -density -distribution density obtained when averaging on measured due in event energies compared to results of or experiments tiplicities ATLAS-CSCdistribution tune is above todata analyses are different. distribution models. η <.For is consistent obfrom se we value have calcustatistics reproduces uncertainty obtained to model comparisons we have used this region. At.6 TeV, ATLAS-CSC tune of PYTHIA tained in -density analysis. This is an Our data are consistent UA proton antiproton latedqgsm mean tunes first reduced moments minimization procedure, as expected. [6], three different of PYTHIA, tune D6T [9],, to some extent, PHOJET are close to data. important consistency check, since correction methods measurements at 9 GeV (Fig. 8a Table ). ena second arises from are influence of tune ATLAS-CSC [] tune Perugia- [], PHOratios of data contribution over Monte Carlo calculations very simq q have been developed by regularization distribution. bias introduced tunes ilar in all threeon ranges suggests that by JET []. Cq PYTHIA Nch / N () ch, regularization was estimated using prescription de- three independent groups extensively comparing Monte scribed in [9] is significantly lower than statistical Carlo to underlying-event minimum-bias including PYTHIA 6 [], PYTHIA 8 [], PHOJET [], EPOS LHC [6]. It can be observed in Fig. that model in best agreement s = TeV data is PYTHIA 6. se results will be used furr tuning of Monte Carlo generators. Multiplicity charged-particle, P(Nch ), were measured at s =.9.6 TeV, as shown in Fig.. experimental data were compared to results from PHOJET three PYTHIA 6 tunes (Perugia-, ATLAS-CSC, D6T). best agreement data is achieved by PHOJET at s =.9 TeV ATLAS-CSC tune of PYTHIA 6 at s =.6 TeV. Furrmore, were scaled by mean to obtain distribution of z = Nch /hnch i, also shown in Fig.. hyposis that of hnch ip(z) are independent of center-of-mass energy is known as KNO scaling [9], se experimental results indicate that KNO scaling holds up to approximately z =. Two-particle correlations Two-particle angular correlations, which are in relative azimuthal angle ( ϕ = ϕtrig ϕassoc ) relative ( η = ηtrig ηassoc ) between trigger associated particles, are

4 EPJ Web of Conferences used to study many aspects of physics of heavy-ion collisions, in particular jet fragmentation collective effects. In elementary collisions small collision systems such as pp y show characteristic features attributed to jet production, while in heavy-ion collisions same jet features are observed in addition to structures around ϕ = (nearside) ϕ = π (awayside) extended in η. se long-range correlations, known as ridges, are often attributed to hydrodynamic flow behavior in quark-gluon plasma (QGP) are typically quantified by coefficients of a Fourier cosine series, v n. It was ree surprising when a nearside ridge was observed in high collisions of small systems, pp [] p Pb []. Furrmore, it was observed that in p Pb collisions at s NN =. TeV nearside peak yields are mostly independent of [], meaning that same trigger associated p T same jet population is selected regardless of. This served as justification to subtract correlations in low- events from high- correlation functions in order to remove correlations due to jet minijet fragmentation. This subtraction procedure (illustrated in Fig. ) showed nearside ridge more clearly also revealed a symmetric ridge on awayside [, ]. This double ridge structure was decomposed into Fourier coefficients in order to extract parameter v in p Pb collisions. analysis was repeated identified particles it was observed that v shows similar mass ordering as was observed in Pb Pb collisions. Figure shows v in p Pb collisions as a function of p T unidentified hadrons, pions, kaons, protons []. Results from CMS show similar behavior KS mesons Λ baryons [6]. v in p Pb collisions was also measured two- multi-particle cumulant methods [7 9]. It is important to note, however, that while v measured in p Pb collisions shows qualitatively similar features as v measured in heavy ion collisions, physical mechanism leading to a non-zero v is still under oretical debate presence of v does not necessarily imply existence of hydrodynamics or a QGP in small collision systems.. Muon-hadron correlations In order to gain more inmation about potential collective effects constrain oretical calculations, it is important to measure strength of ridge to larger η to measure dependence of v on. Both of se points are addressed in muon-hadron analysis permed in [], in which correlation functions between muons at ward rapidities charged hadrons < p T,trig < p T,assoc < GeV/c < GeV/c p Pb s NN =. TeV % < p T,trig < p T,assoc < GeV/c < GeV/c p Pb s NN =. TeV 6 % < p T,trig < p T,assoc < GeV/c < GeV/c p Pb s NN =. TeV ( %) (6 %) ) (rad d N assoc d ηd ϕ N trig... ALI PUB 68 η ϕ (rad) ) (rad d N assoc d ηd ϕ N trig.6.. ALI PUB 6 η ϕ (rad) ) (rad d N assoc d ηd ϕ N trig ALI PUB 66 η ϕ (rad) Figure. two-particle correlation functions in p Pb collisions show a nearside ridge in high- collisions (left) while no ridge is visible in low- collisions (center). subtracted distribution (right) reveals a double ridge structure [].

5 XLVI International Symposium on Multiparticle Dynamics v {PC, sub}..... p Pb s NN =. TeV ( %) (6 %) h K π p Figure. v of π, K, p as well as inclusive unidentified hadrons was measured as a function of p T in p Pb collisions at s NN =. TeV subtraction method []. ALI PUB p (GeV/c) T at mid-rapidity are constructed in order to investigate long-range behavior of double ridge structure < η <.. correlations between muons detected in FMS tracklets reconstructed in ITS were measured in high- ( top % of analyzed event sample) in low- (6- %) events. As in [], low- correlations are subtracted from high- correlation functions to remove structures associated jet fragmentation. After subtraction, correlation functions were projected onto ϕ, n fit a Fourier cosine series to extract v muons detected at ward rapidities. resulting v µ {PC,sub} values are shown in Fig. 6 muons heading in proton- Pb-going directions. data are compared an AMPT [] simulation in which muon decay products are scaled to account efficiency of absorber in FMS. In Fig. 6 it is seen that while AMPT qualitatively describes p T -dependence at low p T, re are significant quantitative differences in p T -dependence η-dependence between data model. At high p T (above p T GeV/c), where muon production is dominated by heavy flavor decays, AMPT does not describe data well. This could be because heavy flavor muons have a non-zero v, or parent particle composition or v values in data AMPT are different. ratio of v µ {PC,sub} in Pb-going p-going directions is also shown in Fig. 6 where it is observed to be independent of p T in statistical systematic uncertainties. A constant fit to data points shows that v is (6 ± 6)% higher in Pb-going than in p-going direction. se results are qualitatively in agreement model predictions. However, current oretical calculations cannot be directly compared experimental results, because effects of absorber are included in experimental data (unfolding such effects could not be done in a model-independent way). Future model calculations should use efficiencies provided in [] in order to compare directly to experimental results. 6 Conclusions Single-particle inclusive two-particle correlation measurements are used to characterize pp p Pb collision systems. charged-particle density has been measured across a range of energies including top LHC energy of s = TeV. dependence of dn ch /dη has been compared MC generators in order to furr tune models. were also compared models demonstrate KNO scaling up to z. In two-particle measurements, long-range correlations in are observed up to η η. presence of se correlations is reminiscent of Pb Pb collisions where structures are frequently attributed to hydrodynamic flow, similar mass ordering being observed in both small large systems. v of ward muons in Pb-going direction is observed to be higher

6 EPJ Web of Conferences {PC,sub} µ v ALI PUB 99 p Pb s NN =. TeV VS: ( %) (6 %) Data, Pb going Data, p going AMPT, Pb going AMPT, p going.... p (GeV/c) T (p going) v (Pb going) / v ALI PUB 996 p Pb s NN =. TeV VS: ( %) (6 %) Data AMPT.... p (GeV/c) T Figure 6. (left) v µ {PC,sub} measured in muon-hadron correlation analysis is shown in p- Pbgoing directions compared results from AMPT. (right) ratio of v µ {PC,sub} results in Pb p-going directions is compared AMPT. than in p-going direction. While se features are similar to correlations observed in heavy-ion collisions, furr oretical phenomenological investigations are needed bee any inferences about collectivity in small systems can be drawn. References [] K. Aamodt et al. (), JINST, S8 (8) [] B. Abelev et al. (), Phys. Rev. Lett., (),.6 [] T. Sjostr, S. Mrenna, P.Z. Sks, JHEP, 6 (6), hep-ph/67 [] T. Sjöstr, S. Ask, J.R. Christiansen, R. Corke, N. Desai, P. Ilten, S. Mrenna, S. Prestel, C.O. Rasmussen, P.Z. Sks, Comput. Phys. Commun. 9, 9 (),. [] R. Engel, J. Ranft, S. Roesler, Phys. Rev. D, 9 (99) [6] T. Pierog, I. Karpenko, J.M. Katzy, E. Yatsenko, K. Werner, Phys. Rev. C9, 96 (), 6. [7] J. Adam et al. (), Phys. Lett. B7, 9 (6), 9.87 [8] K. Aamodt et al. (), Eur. Phys. J. C68, 89 (),. [9] Z. Koba, H. Nielsen, P. Olesen, Nuclear Physics B, 7 (97) [] V. Khachatryan et al. (CMS), JHEP 9, 9 (), 9. [] S. Chatrchyan et al. (CMS), Phys. Lett. B78, 79 (),.8 [] B.B. Abelev et al. (), Phys. Lett. B7, 8 (), 6.6 [] B. Abelev et al. (), Phys. Lett. B79, 9 (),. [] G. Aad et al. (ATLAS), Phys. Rev. Lett., 8 (),.98 [] B.B. Abelev et al. (), Phys. Lett. B76, 6 (), 7.7 [6] V. Khachatryan et al. (CMS), Phys. Lett. B7, (), 9.9 [7] B. Abelev et al. (), Phys. Rev. C9, 9 () [8] G. Aad et al. (ATLAS), Phys. Lett. B7, 6 (),.8 [9] V. Khachatryan et al. (CMS), Phys. Rev. Lett., (),.8 [] J. Adam et al. (), Phys. Lett. B7, 6 (6), 6.8 [] Z.W. Lin, C.M. Ko, B.A. Li, B. Zhang, S. Pal, Phys. Rev. C7, 69 (), nucl-th/