Eur. Phys. J. Plus (2016) 131: 52 DOI 10.1140/epjp/i2016-16052-4

THE EUROPEAN PHYSICAL JOURNAL PLUS

Review

Heavy-ion collisions at the Large Hadron Collider: A review of the results from Run

Néstor Armesto1 and Enrico Scomparin2'3

1 Departamento de Física de Partículas and IGFAE, Universidade de Santiago de Compostela, 15706 Santiago de Compostela, Galicia-Spain

2 INFN, Sezione di Torino - 10125 Torino, Italy

Received: 9 November 2015 / Revised: 4 January 2016 Published online: 2 March 2016

© The Author(s) 2016. This article is published with open access at Springerlink.com

Abstract. We present an overview of the results obtained in pPb and PbPb collisions at the Large Hadron Collider during Run 1. We first discuss the results for global characteristics: cross sections, hadron multiplicities, azimuthal asymmetries, correlations at low transverse momentum, hadrochemistry, and femtoscopy. We then review hard and electroweak probes: particles with high transverse momentum, jets, heavy quarks, quarkonia, electroweak bosons and high transverse momentum photons, low transverse momentum photons and dileptons, and ultraperipheral collisions. We mainly focus on the experimental results, and present very briefly the main current theoretical explanations.

1 Introduction

Quantum Chromodynamics (QCD) [1] is the quantum field theory of strong interactions, see the recent books [2, 3]. While its perturbative domain [4, 5] has been widely explored and calculations are constantly being improved towards higher orders in the expansion in the strong-coupling constant, our knowledge on the non-perturbative regime is far from being satisfactory, see, e.g., the review [6]. Numerical lattice computations are currently able to describe the hadronic spectrum [7] and even light nuclei, but not to deal with dynamic quantities like transport coefficients or cross sections, or with nuclei composed of more than very few nucleons. Therefore, an analytic understanding of some basic features of QCD matter in its usual state, such as confinement and chiral symmetry breaking, is badly needed. One way of getting information about these aspects of QCD is by studying matter in phases where quarks and gluons are no more confined into hadrons and where chiral symmetry is restored [8-13].

Such a phase of matter, named Quark-Gluon Plasma (QGP, see the reviews [14-24]), where the degrees of freedom are quarks and gluons, can be created by colliding heavy ions at the Large Hadron Collider (LHC) at CERN. Its detailed characterization should provide insight into the mentioned, yet unexplained features of QCD that are crucial for understanding hadron and nuclear properties, both static (mass generation and spectra) and dynamic (cross sections). In addition, new high-density regimes of QCD are theoretically expected and signatures of their existence can be explored in both proton-nucleus and nucleus-nucleus collisions. Furthermore, such studies offer the possibility of investigating basic physics questions like the number of interactions between particles required to consider a collection of them as a true medium and the corresponding transition from a microscopic to a macroscopic description, in a relativistic quantum system.

Experiments performed since 2000 at the Relativistic Heavy Ion Collider (RHIC) at BNL, where heavy-ions are collided up to a center-of-mass energy per nucleon-nucleon interaction ^sNN = 200 GeV, have shown that a system is created [25, 26] in which particles are produced with a large degree of coherence, far from a mere superposition of nucleon-nucleon collisions. The obtained energy densities well exceed those above which lattice QCD calculations predict the formation of a QGP, e ~ 1GeV/fm3 [27]. The measured particle spectra are reproduced by viscous relativistic hydrodynamics, if the hydrodynamic regime is assumed to hold at short times after the collision and the shear viscosity over entropy density ratio is small. The yields and expected back-to-back correlations of high-energy partons are strongly suppressed, compared to expectations from extrapolations from pp collisions and from

* Contribution to the Focus Point "Taking stock: measurements, searches and discoveries at the LHC" edited by G. Dissertori.

a e-mail: scompar@to.infn.it

Time [weeks in physics]

Fig. 1. Delivered integrated luminosity per nucleon-nucleon pair at the LHC for collisions involving ions, during Run 1. Taken from [34].

perturbative QCD calculations that include nuclear effects not due to the presence of a hot deconfined medium. Such cold nuclear matter (CNM) effects include in primis the nuclear modification of parton densites (nPDF) [28, 29]. From these findings, it was argued [25, 26] that the formed deconfined partonic medium behaves like an almost ideal liquid very shortly after the collision, and is very opaque to the high-energy partons traversing it.

At the LHC, the much larger ^snn (up to 2.76 TeV for PbPb collisions) is expected to lead to the creation of a denser, hotter, longer-lived medium with a much more abundant rate of various rare particles that can be used to probe it. Previous to the first PbPb and pPb runs, predictions and proposals were compiled in [30-33]. In 2010 and 2011 two PbPb runs took place, and (after a pilot run in September 2012) a pPb run, at y^NN — 5.02 TeV in early 2013. The corresponding delivered integrated luminosities per nucleon-nucleon pair are shown in fig. 1, taken from [34]. In addition, two proton-proton runs at a/s — 2.76 TeV were taken in March 2011 and February 2013, to serve as a reference for the study of nuclear collisions. The delivered luminosities were ~ 250 nb_1 for the former, and 5.6 pb_1 for the latter.

In this paper we review the main experimental results from the full Run 1 of the LHC with ions (other reviews can be found in [35-37]), as well as some of the theoretical discussions triggered by the measurements. The manuscript is organised as follows. In the next section we discuss soft probes of the QGP: cross sections, hadron multiplicities and global event characteristics, azimuthal asymmetries, correlations at low transverse momentum, hadrochemistry, and intensity interferometry between identical particles (femtoscopy). In sect. 3, we review hard and electroweak probes at the LHC, namely: particle production at large transverse momentum, jets, heavy flavours, quarkonia, electroweak bosons and large transverse momentum photons, low transverse momentum photons and dileptons, and ultraperipheral collisions. Note that some of these processes, particularly when studied in pPb collisions, provide the required benchmarking for such hard probes. Needless to say, the organisation that we have chosen is not unique and the discussion of some observables will extend across different subsections. Due to the experimental findings that indicate that pPb collisions cannot be simply considered a reference system but show features in common with PbPb, we address pPb results inside each item related to an observable as an integral part of the ion programme, without separating a specific section for proton-nucleus. We conclude with a summary.

2 Soft probes

One possibility to characterise the medium created in high-energy heavy-ion collisions is through the use of external probes, which we will address in sect. 3. Nevertheless, properties of the medium can also be extracted from the behaviour of the bulk of particle production, i.e. global characteristics like multiplicities or cross sections, correlations,

chemical composition... Such observables, that will be addressed in this section, are called soft probes, as they are dominated by momentum scales of the order of the typical ones of the particles that compose it, as average transverse momentum1. We refer the reader to the reviews [14-17, 30, 31] for extensive information.

2.1 Cross sections

A key quantity to be measured is the total inelastic (or at least non-single diffractive, NSD) cross section that affects the normalisation of many different observables, the translation from yields to cross sections and the determination of the centrality of the collisions. Such a measurement is challenging already in pp collisions. ALICE has measured [38] the visible cross section (i.e. the cross section for events with activity in the kinematic region covered by the experiment) in pPb collisions via van de Meer scans for various trigger conditions. In particular, asking for a coincidence between particle production on both sides of the interaction point, in the pseudorapidity regions 2.8 < n < 5.1 and —3.7 < n < — 1.7, the value GpPb = 2.09 ± 0.07b has been obtained. This quantity was shown to correspond, within less than 1%, to the NSD cross section [39]. CMS [40] has performed an extrapolation to the total inelastic cross section, obtaining 2.061 ± 0.003(stat.) ± 0.034(syst.) ± 0.072(lumi.) b. This estimate was shown to be in agreement with calculations of the same quantity performed in the frame of the Glauber model, starting from the pp inelastic cross section app = 0.070 ± 0.0015 b [41].

Moving to PbPb collisions at ^sNN = 2.76 TeV, ALICE [42] has measured the cross section for electromagnetic dissociation (at large impact parameters beyond the range of the strong interaction) with neutron emission, a very important measurement for luminosity calibration and for understanding the lifetime of the ion beams. The obtained value, gsingle emd = 187.4±0.2(stat)+112(syst) b is much larger than the corresponding hadronic cross section, GPbPb = 7.7 ± 0.1(stat.)+0'| (syst.)b. The latter measurement has been used for centrality determination and validation of the Glauber model in [43].

2.2 Hadron multiplicities and global characteristics

The study of the observables connected with the global characteristics of the event is usually among the first measurements that are performed in the study of ultrarelativistic hevy-ion collisions when a new energy domain opens up. In particular, the measurements of the charged hadron multiplicity as a function of the centrality of the collision and of the pseudorapidity are sensitive to the relative contribution from hard scatterings, which is expected to be proportional to the number of nucleon-nucleon collisions, and soft processes, usually assumed to be proportional to the number of participant nucleons. They also represent an interesting test for theoretical calculations, as they are sensitive to the modelling of the initial state of the collision in terms of gluon saturation. Results from ALICE [44], ATLAS [45] and CMS [46] shown in fig. 2 (left, top) indicate a steep rise in the charged hadron multiplicity density per participant pair (dNch/dn)|n=o/((Npart}/2) when moving from top RHIC energy to the LHC, a feature possibly related to the significant increase of the contribution of hard processes. In particular, the approximately logarithmic increase observed from AGS up to RHIC does not hold up to LHC energy, while a power-law dependence (s0 15) approximately reproduces the data. As was the case at lower energy, dNch/dn|n=0/((Npart}/2) is systematically higher for heavy-ion collisions with respect to pp, showing that nuclear collisions cannot be described as a mere superposition of nucleon-nucleon scatterings. Another observable connected with the global properties of the event is the pseudorapidity distribution of the transverse energy ET [51], shown in fig. 2 (right). It reaches values up to ~ 2.1 TeV for central PbPb collisions and it has a stronger increase (a factor 3.07 ± 0.24) than the charged-particle multiplicity (2.17 ± 0.15) when moving from RHIC to LHC energy, implying that there is a significant increase of the mean transverse energy per particle. Estimates of the energy density according to the Bjorken formula [52] e = (dET/dn)|n=0/(A ■ t0), where A is the overlap area of the incident nuclei and t0 is the parton formation time, yield e =14 GeV/fm3 for a time t0 = 1 fm/c. Such a value is larger by more than one order of magnitude than the estimates for the critical energy density for the phase transition to a deconfined state [27].

The multiplicity per colliding nucleon pair was found, at RHIC energies, to have a mild increase with the collision centrality [50, 53-56]. This situation persists at LHC energy [45,46, 57], where only a moderate variation, by less than a factor 2 from peripheral (75-80%) to central (0-2%) events, is observed, as can be seen in fig. 2 (left, bottom). Remarkably, when the corresponding RHIC results are scaled in such a way to match the observed values for central collisions [57], a good agreement is found over all the centrality range. A similar behaviour was already found in the past when comparing results at RHIC and SPS energy [58]. Two main classes of calculations have been compared to the data, either two-component models based on a combination of perturbative QCD processes and soft interactions [59, 60], or saturation models with various parametrisations for the energy and centrality dependence of the saturation scale [6163]. As can be seen in fig. 3 (left), calculations from both classes of models were able to qualitatively predict the observed behaviour.

1 It is a standard practice to use temperature but it should be kept in mind that the question of the thermodynamical

equilibration of the medium is still an open issue.

• ........ ......... 1 ATLAS Pb+Pb 0-6% ....... ......... -

a ALICE Pb+Pb 0-5% (shifted)

- 0 CMS Pb+Pb 0-5% (shifted)

▼ PHOBOS Au+Au

— ■ SPS Pb+Pb

_ ▲ AGS Au+Au

V ALICE p+p (inel. & NSD) /

" 0 ATLAS p+p

_ A CMS p+p (NSD) / y

NSD p+p/p+p / &

" □ inel. p+p/p+p .,•••',--' .................... t -

— ...... Landau-Carruthers

- --- ALICE power-law JP^

- Log extrapolation g> 9 - atlas _

mil ......... , Pb+Pb \JiNN=2.76 TeV _ ......1 ........

10 102 103 104 Sis and\[sNN [GeV]

0 10 O 10

********

ATLAS Pb+Pb 2.76 TeV ALICE Pb+Pb 2.76 TeV CMS Pb+Pb 2.76 TeV RHIC Au+Au 200 GeV (x 2.15)

"6 4 2 0

FOPI, 0-1% AuAu E802, 0-5% AuAu - NA49, 0-7% PbPb WA98, 0-5% PbPb PHENIX, 0-5% AuAu CMS, 0-5% PbPb RHIC parametrization

■0.46 sNN, V/Snn * 8.7 GeV

50 100 150 200 250 300 350 400

x part'

10 102 \/Snn (GeV)

Fig. 2. Left, top: yS-dependence of the charged particle dWch/dn per colliding nucleon pair from a variety of measurements in pp and pp and central nucleus-nucleus collisions, including the ATLAS 0-6% centrality measurement for \n\ < 0.5 and the 0-5% centrality ALICE and CMS measurements (points shifted horizontally for clarity). The curves show different expectations for the yS-dependence in nucleus-nucleus collisions: results of a Landau hydrodynamics calculation [47] (dotted line), an s015 extrapolation of RHIC and SPS data (dashed line), a logarithmic extrapolation of RHIC and SPS data from [48,49] (solid line); left, bottom: (dNch/dn)\n=0/((Npart}/2) vs. Npart for 2% centrality intervals over 0-20% and 5% centrality intervals over 2080%. Error bars represent combined statistical and systematic uncertainties on the dNch/dn measurements, whereas the shaded band indicates the total systematic uncertainty including Npart uncertainties. The RHIC measurements have been multiplied by 2.15 to allow comparison with the y/snn = 2.76 TeV results. The inset shows the Npart < 60 region in more detail. Right: transverse energy density normalized by ((Npart}/2) for central collisions at n = 0 vs. sNN. The vertical bar on the CMS point is the full systematic uncertainty. The statistical uncertainty is negligible. The solid line is a parametrization of lower-energy data compiled by the PHENIX Collaboration [50]. The dashed line corresponds to a power law fit sNN for ,S

snn > 8.7 GeV.

Measurements of the centrality dependence of the pseudorapidity density distributions were performed for a wide —5.0 < n < 5.5 range [67,68], allowing the investigation of the fragmentation region, i.e., the region of phase space of particles where they have longitudinal momentum close to that of the original colliding nuclei. The main result is shown in fig. 3 (right), where a comparison with theoretical models [63-66] is also performed. The Color-Glass Condensate (CGC) model well describes the data in its (limited) domain of applicability, while no model is able to satisfactorily reproduce both shape and absolute values of the yields over the whole n-range. A double Gaussian fit, with a ratio of the widths independent of centrality, well reproduces the data. When compared to RHIC results [54, 69] as a function of the variable n — ybeam, a scaling is observed in the large-n region (longitudinal scaling). Finally, the total integrated charged multiplicity can be obtained from these data. It reaches Nch = 17165 ± 772 for the 0-5% most central collisions, and its evolution with centrality is identical to the one observed at RHIC [69] once a 2.67 scaling factor is introduced.

2.3 Azimuthal asymmetries

In heavy-ion collisions, the impact-parameter vector and the beam direction define the reaction plane. If the particle momenta do not only depend on the local conditions at their production point, but also on the overall event geometry, a correlation between the azimuthal angle of the produced particles and the reaction plane (anisotropic transverse flow) may be detected, and represents an unambiguous signature of a collective behaviour in the created medium. In particular, in non-central heavy-ion collisions, the initial shape of the created fireball has a geometrical anisotropy, which induces different pressure gradients in the collectively expanding medium. As a consequence, an anisotropy of the azimuthal distributions of the produced particles arises, if interactions in the medium are strong and consequently the mean free path is small enough for relativistic hydrodynamics to be applicable [20, 70, 71] very early —an issue still lacking a theoretical explanation [18, 23]. The coefficients of a Fourier expansion of such distributions [72] can be

С n a.

•v i

) [111 20) [11]

--Armesto et al. [12]

.....- Kharzeev et al. [13]

.....Albacete et al. [14]

200 300 400

4 part'

22002000180016001400-з/ 1200-^ 1000800600400 200

• ALICE

.......AMPT

-- UrQMD — CGC based model

0-5% .5-10% 10-20%

Л в?'*

,&¡¿y .........

__ 20-30%__

-5 -4 -3 -2

Fig. 3. Left: comparison of (dNch/dn)/((Npart)/2) with model calculations for PbPb at y/sNN = 2.76 TeV. Uncorrelated uncertainties are indicated by the error bars, while correlated uncertainties are shown as the grey band. Statistical errors are negligible. The HIJING 2.0 curve is shown for two values of the gluon shadowing (sg) parameter. Right: dNch/dn in different centrality classes compared to model predictions [63-66].

0.08 0.06 0.04 0.02 0

-0.02 -0.04 -0.06 -0.08

-1—I I I I I

Г1-1-1—I I I I I I I-1-1—I I I I I 11-1-1—I I I I I I

/ ALICE —

f STAR -

T i PHOBOS —

PHENIX -

NA49 -

T CERES

t E877 -

é EOS -

E895 -j

f FOPI -

—I_I I I I I 11_I_I_I I I I 111_I_I_I I I I I il_I_I_I I I I I I

0.3 0.25

0.25 0.2 0.15 0.1 0.05

Ж v2{2} 40-50% ▼ v2{4} 40-50% Qv2{4} (STAR)

• 10-20% 20-30% 30-40% 110-20% (STAR) I 20-30% (STAR) ^^ 30-40% (STAR)

V®nn (GeV).....p (GeV/c-

Fig. 4. Left: the ALICE result on pT-integrated elliptic flow at 2.76 TeV in Pb-Pb 20-30% centrality class [80] compared with results from lower energies taken at similar centralities [81,82]. Right: (a) v2{pT), measured by ALICE, for the centrality bin 40-50% from the 2- and 4-particle cumulant methods and for Au-Au collisions at ^JsNn = 200 GeV (STAR). (b) v2{4}(pT) for various centralities compared to STAR measurements.

used to characterize the observed asymmetry. Specifically, the second coefficient of the expansion (elliptic flow, v2) is strongly related to the pressure gradient difference, but also higher-order coefficients (v3, v4,... ) provide valuable information on the evolution of the medium and on the geometric fluctuations of the initial state. At RHIC [73], the observation of v2 values at low pT in good agreement with ideal hydrodynamic calculations (i.e. zero viscosity) in semi-central and central collisions were at the basis of the claim for the production of a strongly interacting QGP phase behaving as an ideal liquid [74-77]. In addition, a hierarchy of v2 values was observed for identified particles, with higher values, at a given pT, for lighter particles [78]. This characteristic ordering can be mainly understood in terms of a radial flow effect and models including a partonic phase, characterized by a harder equation of state (and consequently a higher speed of sound) with respect to the hadronic phase, showed quantitative agreement with the results [79].

At the LHC, anisotropic flow measurements have reached an unprecedent accuracy, due to the larger charged multiplicities and to the large acceptance of the experiments. Concerning v2, it was observed (see fig. 4 (left)) that the pT-integrated values increase by about 30% with respect to RHIC, the increase being mainly due to the larger (pT) (as v2 increases with pT at relatively low pT) [80]. At fixed pT, as shown in fig. 4 (right), the observed values are similar at RHIC and LHC [80,83,84]. It is interesting to note that such a similarity extends down to low RHIC energies [85].

CMS Lint = 150 nb-1 °-1°% PbPb (t

Pt (GeV/c)

Pt (GeV/c)

pT [GeV]

Fig. 5. Left: the single-particle azimuthal anisotropy, v2, as a function of the charged-particle transverse momentum from 1 to 60 GeV/c with \n\ < 1 for six centrality ranges in Pb-Pb collisions at ^JsNN = 2.76 TeV, measured by the CMS experiment (solid markers) [84]. Error bars denote the statistical uncertainties, while the grey bands correspond to the small systematic uncertainties. Comparison to results from the ATLAS (open squares) and CMS (open circles) experiments using data collected in 2010 is also shown. Right: v"!"/vl/2 vs. pT for several centrality intervals [86], measured by ATLAS. The shaded bands indicate the total systematic uncertainties.

The pT-reach of the elliptic flow measurements has been pushed to very-high pT (> 50 GeV/c) [84], as shown in fig. 5 (left). Here the interpretation in terms of collective effects does not hold any more and the observed non-zero v2 is mainly related to the different pathlength in the medium for particles emitted in-plane and out-of-plane, which induces a different energy loss (see sect. 3.1). Higher-order azimuthal asymmetries have also received considerable attention, with measurements extending up to v6 [86-88]. Their study is particularly important to help constraining the initial conditions of the collision, which play a key role in the determination of the viscosity. The measured flow coefficients show little dependence on centrality, consistent with an anisotropy primarily associated with fluctuations in the initial geometry. The approximate scaling v^n(pT) « V2/2(pT), expected in a hydrodynamic scenario, has been observed for all but the most central events (see fig. 5 (right) [86]).

At LHC energies, first systematic studies of event-by-event flow coefficients have been carried out [89]. They are particularly relevant for the assessment of fluctuations in the initial geometry, and were extensively compared to models. The results show that the vn distributions broaden from central to peripheral collisions, and are in good agreement with a fluctuation-only scenario, except for v2. Calculations based on eccentricity distributions predicted in the Glauber model [90] or in the MC-KLN model (which includes saturation effects) [91] do not describe the data, as can be seen in fig. 6 (top). Other approaches, such as EKRT [93] and IP-Glasma [94], were then shown to reach a better agreement with the ATLAS results.

Detailed studies of the elliptic flow of identified particles (n±, K±, KS0, p+p, A+A, 5- + 5 +, Q- + Q+) have also been carried out [92], see fig. 6 (bottom). At low momentum, the mass ordering already seen at RHIC, and related to an interplay of radial flow (i.e., azimuthally symmetric) and anisotropic flow, has been confirmed. At higher pT values (pT > 3 GeV/c) particles tend to group according to their type (mesons vs. baryons). However, systematic comparisons of v2/nq, where nq is the number of constituent quarks, as a function of pT/nq (or (mT — m0)/nq, where m0 is the mass of the particle), show that this scaling is only approximate at LHC energies. Comparisons with hydrodynamical calculations coupling viscous fluid dynamics for the QGP with a hadronic cascade model (VISHNU [95]) show a similar mass ordering, indicating that a rather small value of the viscosity to entropy density ratio is favoured (n/s = 0.16, about two times the 1/4n lower bound [96], is used at the beginning of the hydrodynamic evolution of the system).

: jwii^. ' ' ATLAS Pb+Pb :

V5NN=2.76 TeV '

Lint = 7 ^b-1 -

_ centrality: 0-1% ^^ -

; -»-data pT>0.5GeV,|ti|<2.5 :

. Glauber 0.36e2 W\ "

_-- MC-KLN 0.31e2 0.02 0.04 . . i . >\ 0.06

. ■^tigg. ATLAS Pb+Pb

^ &2.76 TeV :

Lint = 7 ^b-1 ■

I centrality: 5-10%

-».data p>0.5GeV,|n <2.5 #

: Glauber 0.46e2 i V \:

; - - MC-KLN 0.30e2

ATLAS Pb+Pb

(5NN=2.76 TeV = 7 ^b'1

: —^ ATLAS Pb+Pb ■

: y i5NN= 'S Lin 2.76 TeV • = 7 ^b-1

L< 1 o J centrality: 40-45% * J

E .^.data PT>0.5GeV,|n|<2.5 \

' Glauber 0.31e2

: - -MC-KLN 0.26E2 ■ .... i .... i .1. T

0.15 ( 10

r ATLAS Pb+Pb -

: /V A 7f' /•/ &2.76 TeV : L„ = 7 ^b-1 ; \»\ ^

■ centrality: 20-25% • ■ " • • V T -»- data p>0.5GeV,|n|<2.5 ^ : : Glauber 0.40e2 :

' - - MC-KLN 0.29e2 0.05 0.1 0.15 0.2 0.

: -«■""^' ■"■n^, atlas pb+pb :

f V » \T5NN=2.76 TeV ■ » Lint = 7 lib-1 "1 \ •

; centrality: 55-60% ' \ •

.»data pT>0.5GeV,n|<2.5 i • _

: Glauber 0.24e2 ' !

' - -MC-KLN 0.21e2 ■ il t -

"F 0.2 oo 0.1

0.30.2 0.1 0

ALICE Pb-PbSNi = 2.76 TeV |y| < 0.5

10-20%

20-30%

Particle species

-e-rc1

30-40%

40-50%

50-60%

p (GeV/c)

Fig. 6. Top: the event-by-event v2 distributions measured by ATLAS [89] and compared with the e2 eccentricity distributions from two initial geometry models: a Glauber model (solid lines) and the MC-KLN model (dashed lines). The e2 distributions have been rescaled to the same mean values. The scale factors are indicated in the legends. Bottom: the pT-differential v2 for different particle species grouped by centrality class in Pb-Pb collisions at sNN = 2.76 TeV, from ALICE [92].

Fig. 7. Top: the Fourier coefficient v2{2PC,sub} for hadrons (black squares), pions (red triangles), kaons (green stars) and protons (blue circles) as a function of pT from the correlation in the 0-20% multiplicity class after subtraction of the correlation from the 60-100% multiplicity class, from ALICE [97]. Error bars show statistical uncertainties while shaded areas denote systematic uncertainties. Bottom: the coefficients v2 (top row), v3 (middle row) and v4 (bottom row) measured by ATLAS as a function of pT compared between p-Pb collisions with 220 < N^h < 260 [98] and Pb-Pb collisions in 55-60% centrality from ref. [86]. The left column shows the original data with their statistical (error bars) and systematic uncertainties (shaded boxes). In the right column, the same Pb-Pb data are rescaled horizontally by a constant factor of 1.25, and the v2 and v4 are also down-scaled by an empirical factor of 0.66 to match the p-Pb data.

Finally, important results have been obtained in the study of identified particle v2 in pPb collisions, estimated from the study of two-particle correlations. Here, for high-multiplicity pPb collisions (0-20%) one observes (see fig. 7 (top) [97]) a mass ordering similar to the one seen in PbPb, with vp < v^ for 0.5 < pT < 1.5GeV/c and larger beyond pT ~ 2.5GeV/c. There is also a hint of vK < v2 below 1GeV/c. The striking similarity between PbPb and (central) pPb results, with the former being described by hydrodynamical calculations, suggests a collective behaviour also in pPb collisions at LHC energies. For unidentified charged particles, the coefficients from v2 to v5 have also been measured [98, 99]. Comparing central pPb results with the corresponding peripheral PbPb collisions at a similar hadronic multiplicity, and after applying a scale factor to take into account the difference in the (pT) between the two systems [100], a good agreement of the vn distributions is found (fig. 7 (bottom)), in particular in the low-pT region where the statistical uncertainties are small.

CMS Lint = 3.9 nb-1 PbPb \|sNN = 2.76 TeV 3.0 < p"'" < 3.5 GeV/c 1.0 < p^80" < 1.5 GeV/c

Pb-Pb 2.76 TeV, 0-2% central

1.015 1.01 1.005 1

1.002 1

- h / V ,» \ 2 < pT < 2.5 GeV/c 1.5 <pT < 2 GeV/c 0.8 < |An| < 1.8

- f \ A^. A

r \ / \»A

\ 4 V w * X2/ndf = 33.3 / 35

A^ [rad]

Fig. 8. Left: two-dimensional (2D) per-trigger particle associated yield of charged hadrons as a function of \A"q\ and for 3 < pTriS < 3.5 GeV/c and 1 < pTssoc < 1.5GeV/c, for twelve centrality ranges of PbPb collisions at 2.76 TeV/nucleon (the near-side peak is truncated in the two most peripheral distributions to better display the surrounding structure), from CMS. Right: correlation for particle pairs at \A^\ > 0.8, with the Fourier harmonics for n = 1 to 5 superimposed in colour, their sum shown as the dashed curve and the ratio of data to the n < 5 sum shown in the lower panel, from ALICE. Taken from [117] and [115].

2.4 Correlations at low transverse momentum

Correlations in azimuthal angle were treated in the previous subsection, and those in momentum for particles with similar rapidities will be treated in subsect. 2.6. Here we review event-by-event fluctuations and correlations, chargeevent plane correlations, and long-range rapidity correlations. The latter could be included in the previous subsection but here we focus on their long-range nature in rapidity, while usually measurements of azimuthal asymmetries are performed without such consideration.

The analysis of event-by-event fluctuations was proposed as a probe of the matter generated in high-energy heavy-ion collisions, see the review [101]. In the presence of either a second-order phase transition from the QGP to a hadron gas, or the existence of a critical point in the phase diagram of strongly interacting matter, strong fluctuations of thermodynamic quantities, such as temperature, or of global conserved charges, may occur. This could be reflected in dynamical event-by-event fluctuations of the mean transverse momentum of final-state charged particles or in netcharge fluctuations studied by ALICE in [102] and [103], respectively. The results indicate the existence of non-trivial collective effects in the collisions but are not yet conclusive about the signal of a well-defined, finite-order phase transition2.

On the other hand, charge separation with respect to the event plane has been suggested as a signal of the Chiral Magnetic Effect that would point to the existence of domains in the created medium in which discrete symmetries like parity are violated [104]. The creation of strong magnetic fields in heavy-ion collisions, together with the existence of an imbalance in the populations of different helicity states, would produce an electric current and thus charge separation. ALICE has measured such separation [105] with a magnitude similar to that found at RHIC. Whether this is a signal of the Chiral Magnetic Effect or not is still an open issue.

Besides, the charge balance functions (describing the conditional probability that a particle in one momentum bin will be accompanied by a particle of opposite charge in another momentum bin) have been suggested to be sensitive probes of the properties of the system, providing valuable insight into the charge creation mechanism and into fundamental questions concerning hadronisation in heavy-ion collisions [106]. These balance functions are sensitive to the time of creation of the pairs of opposite-charge partons or particles. The charge correlations measurements by ALICE [107,108] indicate that charge generation gets a sizeable contribution from pre-hadronic stages.

Now we turn to long-range correlations in rapidity. Simple causal arguments indicate that such correlations must have their origin in the initial stages of the collisions [109]. In AuAu collisions at RHIC [110,111], a large correlation in the difference in rapidity and azimuthal angle for pairs of particles, extended along several units of rapidity and collimated in azimuth at 0 (near side) and 180 (away side) degrees (thus named "ridge"), was observed for triggered

2 Lattice calculations at high temperatures [27] indicate that the transition from confined to deconfined matter is a crossover.

CMS pPb x/s^ = 5.02 TeV, N°™ne < 35

1 < pT < 3 GeV/c

CMS pPb \/sNN = 5.02 TeV, N°™ne > 110

1 < pT < 3 GeV/c

- 1 1 1 1 I '_1 1 1 I f1 1 ' , '

atlas p+Pb ^NN=5.02 TeV, J L = 1 цЬ-1 0.5<pT,b<4 GeV, 2<| n|<5A -9-EEpb>80 GeV bCYAM=14.3

□ □ □ 1

AtLaS p+Pb 1 ysNN=S.02 TeV, J L = 1 ¿-1 '

0.5<pT,b<4 GeV, 2<| n|<5a (d) -Щ- Near: |Дф|<п/3

Away: |Дф|>2п/3 я

-♦- Difference ■ '

♦ ♦ ♦

Д3ф|

<2EPb > [GeV]

Fig. 9. Top: 2D two-particle correlation functions for 5.02 TeV/nucleon pPb collisions for pairs of charged particles with 1 < pT < 3GeV/c; results are shown (a) for low-multiplicity events (№®lne < 35) and (b) for a high-multiplicity selection (N°fine > 110), and the sharp near-side peaks from jet correlations have been truncated to better illustrate the structure outside that region. Bottom: (c) the per-trigger yield distribution together with pedestal levels for peripheral (bpYAM) and central (&zYAM) events, and (d) integrated per-trigger yield as a function of the sum of the transverse energy on the Pb side for pairs in 2 < \ An\ < 5. Taken from CMS [118] and ATLAS [120].

-B- 2E7<20 GeV bPYAM=3.2

and untriggered analysis, with a peak at (0, 0) in (Ац,Аф). Later, a strikingly similar structure was found in high-multiplicity pp collisions at the LHC by CMS [112,113] and ATLAS [114], and then in PbPb collisions [115-117]. The correlation becomes stronger with increasing centrality. As shown by ALICE, its structure in azimuth can be characterised by a Fourier decomposition even for pairs with large rapidity separations, see fig. 8, which provided a new view on this type of correlations.

Both the rapidity and azimuth structures have also been observed in pPb collisions [98,118-121], with again a clear transition from the absence of a near-side ridge to a large signal with increasing centrality, see fig. 9, while the jet signal remains roughly constant with centrality [122]. Finally, the Fourier azimuthal coefficients of identified particles extracted in pPb [123] showed the mass ordering expected if such azimuthal correlations would come from a collective flow driven by relativistic hydrodynamics. In addition, correlations between strange particles and charged hadrons in pPb and PbPb were also measured [124]. A comparison between the yields in pPb and PbPb at the same multiplicity [99], see fig. 10, shows different features and magnitudes with increasing multiplicity and transverse momenta for correlations of short and long range in rapidity, but a very similar size for both colliding systems. All these observations, together with the measurement of up to 8-particle correlations in pPb collisions [125], have triggered an intense debate on the origin of this phenomenon and on the possibility of large degrees of collectivity in very high-multiplicity pp and pPb collisions at LHC energies.

Many theoretical explanations have been proposed for the ridge. While not the only ones, see others in [126-129], two main lines are currently under debate. One is the pure hydrodynamical explanation that considers that, provided an initial condition that contains long-range correlations in rapidity, the coupling to an expanding medium automatically produces the azimuthal correlation [130-134]. The other one, based on ideas of the CGC that indeed provides a long-range rapidity correlation, considers that the azimuthal asymmetry originates from the initial correlations in the wave function of the colliding hadrons [135-140], e.g., Bose enhancement for gluon and the existence of domains in the hadron with oriented chromoelectric fields or density gradients. While explanations based on CGC ideas have

(a) |ani>2 PbPb \TsNN = 2.76 TeV -# pPb \rSNN = 5.02 TeV, 2013-CMS DpPb (iNN = 5.02 TeV, 2012 " O pp lTi = 7 TeV

pPb CGC

"—Q0(proton)=0.336 GeV2 ■ ■ • Q°0(proton)=1.008 GeV2 -..Q0(proton)=1.680 GeV2

(b) |an|<1 minus |ar)|>2

1 < pT 1 < pT

g < 2 GeV/c ioc < 2 GeV/c

■S 0.2-

Ï' ^HlUÎ'miV

(a) |ani>2 220 < NorkMne < 260

-cms 1 < pTiioc < 2 GeV/c

PbPb ^iNN = 2.76 TeV • pPb IiZ = 5.02 TeV

(b) |an|<1 minus |an|>2

200 .offline

200 ..offline

5 10 0

pt^GeV/c)

pt^GeV/c)

Fig. 10. Left: associated event-normalized yields for the near-side correlation function as a function of multiplicity Nffline for 1 < p^18 < 2 GeV/c and 1 < p!ssoc < 2 GeV/c averaged over the (a) long-range {\A."q\ > 2) region and (b) short-range < 1) region, from which the event-normalized yield of the long-range region is subtracted; the results for 7TeV pp collisions (open circles) and 5.02TeV/nucleon pPb collisions from the 2012 run (open squares), as well as calculations from the Color-Glass Condensate (CGC) theory (curves), are also shown. Right: associated event-normalized yield for the near-side correlation function integrated over the region \A$\ < 1.2, averaged over the (a) long-range (\A^\ > 2) region and (b) short-range < 1) region, from which the event-normalized yield of the long-range region is subtracted; the results are shown

as a function of p!^8 at 1 < pTssoc < 2 GeV/c for events with 220 < Ntoffline < 260 for 5.02/nucleon TeV pPb collisions (filled circles) and 2.76TeV/nucleon PbPb collisions (filled squares). Taken from [99].

been successful in reproducing quantitatively the ridge in pp, flow is required to reproduce its magnitude in PbPb and the situation in pPb is still unclear. This intense activity is linked with the attempts to understand the origin of the fluctuations in the initial conditions for relativistic hydrodynamics, see, e.g., [141]. In any case, the apparent collectivity suggested by pp and pPb data is among the largest surprises from Run 1 at the LHC, although such possibility had already been put forward [142,143].

2.5 Hadrochemistry

The study of the chemical composition of the particles produced in heavy-ion collisions has been a key feature of the experimental programmes and theoretical discussions for decades [144]. The success of the description of hadrochemistry assuming an equilibrated system ruled by the grand canonical ensemble, the statistical model of hadronisation [145], has been taken as evidence of both the approximate equilibration and the partonic nature of the produced medium, as well as of the simple statistical nature of the process.

In PbPb collisions at the LHC, ALICE has made a wealth of pT-differential analyses of identified particle production. In [146] and [147], the relative yields of multistrange hadrons and of $ mesons have been measured. As shown in fig. 11, the $ mesons exhibit, with respect to pions, a relative enhancement with increasing pT that closely follows that of protons, suggesting that a mass effect is present, as resulting from hydrodynamic calculations where the radial flow is linked to the mass of the particles (but other options exist, see, e.g., [151]). The relative ratios of multistrange particles and A's with respect to pions increase gradually with centrality, and agree with the expectations of the statistical models for mid-central and central collisions. The same trends have been observed in the A/K0 ratios [152], or for the transverse momentum and centrality dependence of pions, kaons and protons [153,154] that show a behaviour compatible with relativistic hydrodynamics.

In pPb collisions, see fig. 12, both ALICE [157] and CMS [156] have measured a hardening of the transverse momentum spectra of different particle species, larger for particles with larger masses, and smoothly increasing with centrality. This hardening with increasing hadron multiplicity is present in both pp, pPb and PbPb and is more pronounced for smaller systems, a feature also observed for unidentified charged-particle production. Noticeably, the observed mean transverse momentum or the extracted effective temperature show mass ordering like in PbPb. Besides, ALICE has measured the spectrum of $ mesons [158], that shows a large enhancement (compared to scaled pp) for intermediate transverse momentum 3-4GeV/c) in the backward (Pb-going) region and suppression at smaller transverse momentum in the forward (p-going) region.

The spectra of at mid-rapidity has been measured for different centralities in PbPb collisions by ALICE [159]. Their nuclear modification factor (see the definition in sect. 3) shows a suppression slightly larger than that found at RHIC and compatible with the one of charged hadrons that we will discuss in more length in subsect. 3.1. Interestingly, see fig. 13, ALICE [160,161] has also studied the nuclear modification factor for pions, kaons and protons for different centralities and up to pT ~ 20GeV/c. For all centralities, the nuclear modification factors for all species, and the relative chemical composition in pp and PbPb, become equal above pT ~ 8GeV/c, indicating that hadronisation of such particles is not affected by the medium and that the origin of the suppression is partonic —as already observed

Fig. 11. Left: ratios of p and 0 yields to charged pions as a function of pT for central PbPb collisions at 2.76 TeV/nucleon; the p/n ratio is presented using two pT binning schemes: the ratio with its original measured bins is shown along with a recalculated version that uses the same bins as the 0 meson pT distribution for 0.5 < pT < 5GeV/c. Right: ratios of particle yields to charged pion yields for PbPb collisions at 2.76 TeV/nucleon, AuAu collisions at 200 GeV/nucleon, and pp collisions at 200 GeV, 900 GeV and 7TeV; the lines show ratios given by grand-canonical thermal models with temperatures of 170 MeV [148] (upper dashed lines), 164MeV [149] (solid lines), and 156MeV [150] (lower dashed lines); the total uncertainties (including centrality-uncorrelated and centrality-correlated components) are shown as bars. Some of the measurements at Npart = 2 have been shifted horizontally for visibility; the two most central 0/n values for AuAu collisions are for overlapping centrality intervals (0-5% and 0-10%). Taken from ALICE [147].

at RHIC but for smaller transverse momenta [162-164]. This poses strong constrains on ideas that consider non-perturbative [165] or perturbative modifications of hadronisation inside the medium [166].

Finally, ALICE has also presented results [167-169] on weakly decaying bound states of A, A and (anti)nucleons, and on hypertritons, tritons, deuterons and their corresponding antinuclei. Furthermore, studies of the mass-overcharge differences between light nuclei and antinuclei confirm CPT invariance with improved precision in the nuclear domain [170].

2.6 Femtoscopy

The size of the particle-emitting region at freeze-out can be extracted via two-particle correlations of identical particles, usually pions (referred to as Bose-Einstein, or Hanbury-Brown-Twiss "HBT", correlations). This technique, also known as femtoscopy, was first used in particle physics in the 1960s [171] and measures the apparent width of the distribution of relative separation of emission points, which is conventionally called the "radius parameter". The radius can be determined independently for three directions: Rlong along the beam axis, R0ut along the pair transverse momentum, and Rside, perpendicular to the other two. Femtoscopic analyses were extensively carried out already at AGS, SPS and RHIC energies, providing valuable information [172]. Among the main findings in nucleus-nucleus collisions, a linear increase of the radius parameters was measured as a function of (dNch/dn)1/3. In addition, a significant decrease of the radii with the momentum of the pair was also established. This is a characteristic feature of expanding particle sources since the HBT radii describe the homogeneity length [173] rather than the overall size of the particle-emitting system. The homogeneity length is defined as the size of the region that contributes to the pion spectrum at a particular three-momentum p. The decrease of the size with kT, defined as half of the modulus of the vector sum of the two transverse momenta, is observed in experimental data from heavy-ion collisions at all centralities, various collision energies and colliding system types, and is well described quantitatively in hydrodynamic models [174,175].

At the LHC, the ALICE experiment has performed an extensive set of femtoscopic measurements [176-181]. In fig. 14 (left) the results for 0-5% Pb-Pb collisions at ^JsNN = 2.76 TeV [176] are compared, separately for the three radius parameters, to lower energy experiment results. The LHC values are higher than the top RHIC energy values by 10-35%, and the scaling with (dNch/dn)1/3 is seen to hold up to LHC energies. The product of the radii is connected to the volume of the homogeneity region and is shown in fig. 14 (right) for kT = 0.3GeV/c [176]. Compared with previous results it shows a linear dependence on the charged-particle pseudorapidity density and is two times larger at the LHC than at RHIC. Estimates of the decoupling time for pions, based on Riong measurements, exceeds 10fm/c, and are 40% larger than at RHIC.

Fig. 12. Left: average transverse momentum of identified charged hadrons (pions, kaons, protons) in the range \y\ < 1 as a function of the corrected track multiplicity for \n\ < 2.4, for pp collisions (open symbols) at several energies [155], and for pPb collisions (filled symbols) at 5.02 TeV/nucleon; both (pT) and yield ratios were computed assuming a Tsallis-Pareto distribution in the unmeasured range; lines are drawn to guide the eye (gray solid pp 0.9 TeV, gray dotted pp 2.76 TeV, black dash-dotted pp 7TeV, coloured solid pPb 5.02TeV/nucleon); the ranges of (pT) values measured by ALICE in various centrality PbPb collisions at 2.76 TeV/nucleon [154] are indicated with horizontal bands. Right: inverse slope parameters T' from fits of pion, kaon, and proton spectra (both charges) with a form proportional to pT exp(-mT/T'); results for a selection of multiplicity classes, with different (Ntracks) as indicated, are plotted for pPb data; the curves are drawn to guide the eye. Taken from [156].

Fig. 13. Left: nuclear modification factor RAA as a function of pT for different particle species; results for 0-5% (left) and 60-80% (right) collision centralities are shown. Right: particle ratios as a function of pT measured in pp and the most central, 0-5%, PbPb collisions; the theoretical predictions refer to PbPb collisions. Taken from [160].

As a further study, two- and three-pion correlations of same and mixed charge were measured at low relative momentum to estimate the coherent fraction of charged pions in PbPb collisions [177]. The presence of coherent production reduces the maximum of Bose-Einstein correlations and the estimated fraction was found to be 23 ± 8% for low triplet momentum. Comparisons of HBT parameters were also carried out between pp, pPb and PbPb (see fig. 15 (left), where the radii obtained from a Gaussian fit of the two- and three-pions correlation functions are reported) [178]. At similar multiplicity, a modest increase of the parameters was found between pp and pPb (5-15%), while the difference between pPb and PbPb is larger (35-55%). The weak increase between pp and pPb favours models with a CGC initial phase [132] with respect to those assuming a hydrodynamic evolution of the system [191]. Finally, femtoscopic analyses have also been carried out in Pb-Pb collisions for identified particles, including , K ±, Kg, p and p [180]. Common mT scaling of the n and K parameters, with a decrease of the radii with increasing mT, were predicted to be a robust signature of hydrodynamic flow [192]. Furthermore, K results are considered as a cleaner signal due to a smaller decay contribution, and (anti)proton results provide a possibility for checking if baryons are included in the collective motion. In fig. 15 (right), results on Rinv, the invariant radius from two-particle correlations evaluated in the pair-rest frame are shown. For overlapping mT, the radius parameters are mostly consistent with each other within uncertainties, though the pion radii are generally larger than the kaon radii (the latter effect can be explained as a consequence of the increase of the Lorentz factor with decreasing particle mass [193]).

cT E 400 1 1 1 1 1 1 1 1 1

▲ E895 2.7, 3.3, 3.8, 4.3 GeV

M— 350 — A NA49 8.7, 12.5, 17.3 GeV

c a ; ■ CERES 17.3 GeV

OC^ 300 — •¿r STAR 62.4, 200 GeV

■a 'cn ; □ PHOBOS 62.4, 200 GeV

CCr □ 250 E- • ALICE 2760 GeV

0 DC 200

150 100 50

<dN /dn )

Fig. 14. Left: pion HBT radii at kT = 0.3GeV/c for the 5% most central PbPb at ^snn = 2.76 TeV (red filled dot) and the radii obtained for central gold and lead collisions at lower energies at the AGS [182], SPS [183-185], and RHIC [49,186-190]. Right: product of the three pion HBT radii at kT = 0.3 GeV/c. The ALICE result (red filled dot) is compared to those obtained for central gold and lead collisions at lower energies.

O .E g

pp (%=7 TeV p-Pb ^=5.02 TeV Pb-Pb &2.76TeV

Two-Pions

o □ o -

Three-Pions-

• ■

te iaf -

O ObOoOOo

n~ n~ \|/ K± K±

0-10% 10-30% 30-50%

'it* o *

1.6 1.8 2

<mT) (GeV/ c2)

Fig. 15. Two- and three-pion Gaussian fit parameters vs. {Nch} in pp, pPb and PbPb collision systems for low average values of the two- and three-pion average transverse momenta, kT and KT,3. The top panel shows the Gaussian radii Rjnv and Rinv>3. The bottom panel (not discussed in the text) shows the effective Gaussian intercept parameters and AG3 (see [178] for details). Right: Rinv parameters vs. mT for the three centralities considered for , K± K±, KS0K0 , pp, and pp [180]. Statistical (thin

lines) and systematic (boxes) uncertainties are shown.

ALICE Pb-Pb ysNN = 2.76 TeV

The radius parameters show an increase with centrality as would be expected from a simple geometric picture of the collisions. They also show a decreasing size with increasing mT, as would be expected in the presence of collective radial flow. The results are found to be in good agreement with hydrodynamical models [194].

3 Hard and electroweak probes

A powerful way to analyse the medium created in high-energy heavy-ion collisions comes from the study of probes3 created in hard scatterings. Such probes, due to their large scale (mass, energy, transverse momentum, ..., larger than any medium scale like average transverse momentum, temperature, saturation scale, Aqcd, ...), are: (a) clearly separated from the medium that they analyse; (b) accurately computable in perturbative QCD so that the separation between production of the hard probe and propagation through the medium is under theoretical control, and (c) produced very early in the collision. In this way, it is possible to establish their initial flux and to study the change in observed yields or cross sections due to the presence of the medium. Experimentally, this is done by comparing results on pp and pA/AA collisions. The quantity usually employed to quantify such variation is the nuclear modification factor

with (Ncoll) the number of collisions in a given centrality class, usually computed in the frame of the Glauber model. With this normalisation, the nuclear modification factor would be 1 if no nuclear effect exists, i.e., if an AA collision were a mere superposition of nucleon-nucleon ones. A similar quantity, called Rcp, can also be defined using a peripheral class for reference instead of pp. The numerator and denominator may be not yields but cross sections, in which case the normalisation factor is different.

These probes are generically known as hard probes and include high transverse momentum particles, jets, heavy flavours and quarkonium states. In addition, there are also other penetrating probes, so-called electroweak - those provided by particles that do not interact strongly with the medium: electroweak bosons, and real and virtual photons at high and low transverse momentum. Electroweak bosons, and real and virtual photons at high transverse momentum, constitute control probes that check our understanding of perturbative QCD in nuclear collisions, in the absence of any medium, particularly for determining the nuclear modification of parton densities nPDFs4. Low transverse momentum photons and dileptons are instead directly sensitive to the medium, because they can be thermally emitted by the medium itself. We refer the reader to the reviews [14-17,19, 22, 24, 30, 31,195-200] for extensive information.

3.1 Particle production at large transverse momentum

The suppression of particles produced at large transverse momentum and the disappearance of back-to-back correlations was one of the golden measurements at RHIC that established the dense partonic nature of the medium produced in nucleus-nucleus collisions, see [201-205]. At the LHC, analogous results have been found for the nuclear modification factor [159-161,206-209] and for the suppression of back-to-back correlations [210].

In fig. 16, results for the nuclear modification factor of charged particles vs. transverse momentum and number of participants are shown, taken from [209]. The measurements extend to pT > 100 GeV that corresponds to much higher jet transverse momentum. The agreement between different experiments when the results are plotted for similar centralities and rapidity windows is notable. For that, the use of a reference pp spectrum measured at the same centre-of-mass energy, 2.76 TeV, has been crucial (see the uncertainties due to interpolation of the reference in [206]). A steady decrease with increasing centrality is seen for all transverse momenta.

In fig. 17, experimental results of CMS [207] and ALICE [208] for the nuclear modification factor of charged particles are compared with different models that take into account energy loss of the parton traversing the medium whose hadronisation leads to the observed charged particle. This energy loss is considered to be dominantly radiative (medium-induced gluon radiation) at high energies of the parton, but elastic energy loss is included in some models, see [22, 24,199, 200] and references therein. In all these calculations, the medium is characterised by its length and dynamical behaviour, and by the strength of the interaction of the parton with the medium constituents, usually given by the transport coefficient q (the average transverse momentum squared transferred from the medium to the parton per mean free path). Such interaction leads to an enhancement of radiation off the parton that implies a longitudinal medium-induced energy degradation. Several models where the energy loss mechanism is embedded in a realistic

3 Due to the transient nature of the medium, they must be self-generated.

4 The study of hard probes in pA collisions is mandatory not only for the determination of nPDFs but also for the understanding of other CNM effects like absorption of produced bound states in normal nuclear matter and, in general, of factorisation in nuclear collisions.

Raa(V,Pt )

(Ncoll)

d2Npp ' dyd^T

p+p, Pb+Pb

Iii, ^nn = 2.76 TeV

b = 0.15 nb

0 - 5 %

• ATLAS, M<1.0 CMS, M<10 □ ALICE, in!<0.8

PT [GeV]

' 1 ' atlas a, IÍÜN = 2.76 TeV

— „ p+p, Pb+Pb Lpp Hnt = 4.2 pb-1 —

:î Î + î - LPnbPb = 0.15 nb -1 -

" * 5 i - - ï +

_ " * ■ g + + _

- S S 1 I i 1 » «

" in! < 2 i -

maximum (1.7<p <2.0 GeV) s

minimum (6.7<pT<7.7 GeV) »

■ intermediate (19.9<pt<22.8 GeV) f *

+ plateau (59.8<p <94.8 GeV) , , , i , T, , , i ,

300 /N s 400

* pair

Fig. 16. The nuclear modification factor for charged particles measured by ATLAS [209], vs. transverse momentum (left, compared to CMS [207] and ALICE [208] results) and vs. the number of participants (right). Taken from [209].

(GeV/c)

100 200

ALICE, Pb-Pb, ^sNN = 2.76 TeV

charged particles, |n| < 0.8 norm.

Mil."""'"'

• ALICE (0-5%) ■ CMS (0-5%)

■ ■ HT (Chen et al.) lower density

■ ■ HT (Chen et al.) higher density

— HT (A.M.) -.... ASW (T.R.)

— YaJEM-D (T.R.)

_ elastic (T.R.) large P elastic (t.R.) small P^"

■ ■ WHDG (W.H.) n0 upper limit

■ ■ WHDG (W.H.) n0 lower limit

PT (GeV/c)

Fig. 17. The nuclear modification factor for charged particles compared to different models, for CMS [207] (left) and ALICE [208] (right) results. Taken from [207] and [208].

dynamical model for the medium were compared to experimental data. In this way, values of the transport coefficient lying between those expected for weak and strong coupling pictures of the medium have been extracted [211,212]. However, good agreement with the data can be found by various models that mainly differ in the inclusion or not of elastic energy loss, and in the way in which energy and momentum conservation is imposed on top of the theoretical calculation [213], done in the high-energy limit but applied to finite energies. This is due to the fact that the nuclear modification factor is subject to several biases related to the steeply falling partonic spectrum, the different production points in the medium (usually called surface bias) and the hardness of the fragmentation functions. Besides, models assume that hadronisation takes place outside the medium. Energy loss is therefore a purely partonic process, an assumption that may not be fully true up to sizeable transverse momenta that at the LHC may be as large as 10 GeV/c, see subsect. 2.5. Therefore, more differential measurements like those of jets are required to disentangle the various mechanisms of energy loss including full jet-medium interactions, and the weak- or strong-coupling nature of the medium.

To conclude this subsection, we comment on pPb collisions. The nuclear modification factor for neutral pions measured at RHIC was used to provide information about the nuclear modification of parton densities [28, 29]. The charged particle RpPb has been also measured at the LHC. While both ALICE [214,215] and CMS [216] observe a nuclear modification factor clearly compatible with unity for transverse momenta from roughly 3 to 30 GeV/c, CMS data5 show an enhancement as large as 1.4 for larger transverse momenta, while ALICE data —that extend to smaller pT than those of CMS— seem to be compatible with unity. Discussions are ongoing on the possible origin of the

5 Preliminary ATLAS data [217] also indicate an enhancement.

Fig. 18. Distributions in the energy asymmetry (left) and in the azimuthal angle (right) of the two hardest jets (with pt,i > 120GeV/c and pT,2 > 50GeV/c) in pp collisions at 7TeV, and PbPb collisions at 2.76TeV/nucleon for different centralities, see the legends on the plots. Taken from [225].

apparent discrepancy, like the different estimates of the reference pp spectra (based on interpolation procedures) used by the experiments. In such a situation, a pp run at yfs = 5.02 TeV can be considered as mandatory.

3.2 Jets

Jets —sprays of collimated hadrons— are the closest objects to partons that in QCD can be properly defined, see [218]. The study of their modification inside a coloured medium was proposed long ago as a sensitive proof of the mechanisms of energy loss, of the interaction strength with the medium and of its density [219, 220].

While jets have been measured in heavy-ion collisions at RHIC [221,222], their measurements there suffer from limitations due to the accessible energy range and to the limited detector acceptances. On the other hand, at the LHC [223-239] the much higher centre-of-mass energy provides jets with large transverse energy that can be reconstructed above the fluctuations of energy contributions from the background event that is produced in heavy-ion collisions. Furthermore, the large calorimeters of the experiments allow full jet reconstruction and detailed differential studies. In this respect, ATLAS and CMS have traditional large acceptance electromagnetic and hadronic calorimeters, while ALICE reconstructs jets with electromagnetic energy and tracks.

In fig. 18 we show the energy asymmetry AJ = (pT,i — pt,2)/(pt,i + pT,2) and the azimuthal distribution for dijets first studied by ATLAS [223] and then by CMS [225,226]. A clear increase in the dijet asymmetry without sizeable azimuthal decorrelation can be observed. The trend of the data was subsequently verified for 7 +jet final states [227]. Note that the fraction of energy loss due to medium effects (pT,2/pT,1) is of the order of 10%. The missing transverse momentum is recovered in the form of relatively soft particles at large angles [225, 239].

In fig. 19 the nuclear modification factor for jets (both charged [231] and fully calorimetric [224] jets) shows a suppression similar to that found for charged particles (with due caution on the comparison of the different transverse momentum scales), and no sizeable difference for different jet reconstruction parameter in the algorithm. The suppression shows a significant dependence with the azimuthal angle with respect to the reaction plane [230, 238], as expected from energy loss processes.

On the other hand, the reduction in the number of events in which the subleading jet has disappeared (either by degradation below the cut or by migration outside the acceptance) compared to pp collisions is moderate, of the order of 10% for most central collisions, and the contamination from background jets is small, see fig. 20, left [226]. All features behave smoothly and tend to become less and less pronounced with increasing jet energy. The fragmentation function [228, 229, 233, 234] shows an enhancement of soft particles, a suppression of particles with intermediate momentum fractions, and little modification of hard ones (the small enhancement is usually interpreted in terms of the modification of the denominator —the jet energy— due to some energy loss), see fig. 20, right.

Note that jet measurements in heavy-ion collisions require a procedure for background subtraction (see studies for characterising the background in [240]) that, for correlation observables, must be carefully considered as it can mimic medium effects [241].

From a theoretical point of view, the observation of i) a sizeable increase, with respect to the proton-proton case, of the energy asymmetry in dijet systems without modification of their azimuthal correlation, ii) the recovery of missing energy at large angles away from the jet axis in the form of soft particles, and iii) the lack of strong modification of the hard part of the jet fragmentation function, poses serious challenges for the standard explanation of jet quenching

£L 2 c°1.8 1.6 1.4 1.2 1

0.8 0.6 0.4 0.2 0

Ch. particles

(0-10%)/(50-80%) CMS

(0-5%)/(50-90%)

ALICE Ch. Jets R=0.3 (0-10%)/(50-80%) ATLAS Calo Jets R=0.3-(0-10%)/(60-80%)

pTack, j (GeV/c)

2 1.8 1.6 1.4 1.2 1

0.8 0.6 0.4 0.2 0

ALICE Pb-Pb S Anti-fcT

=2.76 TeV

p^ack > 0.15 GeV/c Leading track pT > 5 GeV/c

-,—i—i—i—i—i Pb-Pb 0-10% ^ Pb-Pb 50-80% PYTHIA -

p (GeV/c)

T,ch jet x '

Fig. 19. Left: nuclear modification factor of central collisions compared to peripheral ones, for charged jets, fully reconstructed jets and charged hadrons, vs. the corresponding transverse momentum of the jet or hadron. Right: ratio of spectra of charged jets defined with reconstruction parameter R = 0.2 and R = 0.3, vs. charged jet transverse momentum. Taken from [231].

„ 0.95

Z" 0.85

# 0.8 Z

0.75 _ 0.08

• ' 1 1 1 cms ;

M B * J sD . • J :# a PbPb ysNNN = 2.76 TeV, J*Ldt = 150 ^b"1 : T 0 pp Vs = 2.76 TeV, jLdt = 231nb-1 : : PYTHIA+HYDJET - ÍB D n □ ^ - • u □ - : • : ~ • ~ : • -

Centrality 0-20% pT,2 > 30 GeV/c » Â*,2 > h _ : * + °Anï.8 .BIB ■ * - i.-Tii." Pt, > 120 GeV/c Pt2 > 30 GeV/c _ Â^i2 > 2" « - " . d " ^ ■ D ffl, B 1 1 IBi 1 1 ni , , , , 1 , , ,

150 200 250 300 350 0

PT1 (GeV/c)

Pb+Pb\isNN=2.76 TeV 0.14 nb-1

anti-kT R=0.4 " pjet>100 GeV

Pb+P^^/sNN=2.76 TeV 0.14 nb-1

Pb+Pb\fsNN=2.76 TeV 0.14 nb-1

20-30%/60-80%

Pb+Pb\fsNN=2.76 TeV 0.14 nb-1

50-60%/60-80%

* 4 4 » » » 41

Fig. 20. Left, top: ratio of events with a second jet over the total number of leading jets, for pp and PbPb collisions vs. the transverse momentum of the leading jet and vs. the number of participants. Left, bottom: the same as for the left, top panels for the number of background jets over the number of leading jets. Right: ratio of jet fragmentation functions in PbPb and pp vs. the momentum fraction of the energy of the jet carried by the measured hadron. Taken from [226] and [234].

in terms of medium-induced gluon radiation in which energy loss and broadening are intrinsically connected and the emitted radiation is semihard. While several phenomenological explanations, see [242-251], have been put forward (that go from the collimation by scattering-away of the soft components of the jet in the medium, to the important role of colour coherence between emitters and its medium disruption), none provides a complete quantitative description of the data nor proper theoretical justification for some of the assumptions.

In pPb collisions [252-255], the nuclear modification for dijet [252] and single-jet distributions [253,255] in minimum-bias events is very mild and in agreement with the expectations from the nuclear modification of par-ton densities (actually these data can be used to constrain them [256]). Negligible nuclear effects are seen on the dijet acoplanarity [254]. On the other hand, studies on the centrality and pseudorapidity dependence of these observables [252, 253] point to the existence of non-trivial effects coupling the hard process to the centrality estimator, see [257-260] for several interpretation attempts.

3.3 Heavy flavours

Heavy quarks are an important probe of the medium created in nuclear collisions. Due to their mass, much larger than the temperature of the system, they are not created or destroyed in the medium, but rather produced by hard scattering in the early stage of the collision process. Therefore, they are particularly sensitive to the QGP phase. The study of the

A (ALICE) 8<Pt <16 GeV/c, |y |<0.8 ■ D mesons (ALICE) 8<pt<16 GeV/c, |y |<0.5 (empty) filled boxes: (un)correlated syst. uncert.

Djordjevic et al. Phys.Lett.B 737 (2014) 298 ■ D mesons

shifted by +10 in (N )

k i ' i Part

0 50 100 150 200 250 300 350 400

x part'

Pb-Pb,

: 2.76 TeV

D mesons (ALICE) 8<p <16 GeV/c, |y|<0.5 Non-prompt Jfy (CMS) 6.5<p <30 GeV/c, |y|<2.4

JHEP 05 (2012) 063

(empty) filled boxes: (un)correlated syst. uncert.

Djordjevic et al. Phys.Lett.B 737 (2014) 298

--D mesons

--Non-prompt J/^

......Non-prompt J/^ with c quark energy loss-

50 100 150 200 250 300

350 400

Fig. 21. Left: comparison of D meson and charged pion RAA, measured by ALICE in PbPb collisions at ^JsNN = 2.76 TeV [271]. Right: same D-meson results, compared to non-prompt J/^ RAA measured by CMS [278]. The results are compared to the theoretical model of ref. [282].

Pb-Pb, VsNN = 2.76 TeV

modification of their yields and of their angular distribution can give information on the degree of thermalization of the medium and on its transport coefficient. In particular, heavy quarks can lose energy and participate in the collective motion of the created system, the mainly related observables being their nuclear modification factor and their elliptic flow. A definite prediction of the radiative energy loss models is a substantial dependence of the specific energy loss on the mass of the quark [261,262], with the heavy quarks losing less energy than the light ones (often called "the dead cone effect"). Systematic studies of heavy quark production were started at RHIC, where the measurements were essentially based on the detection of the non-photonic electrons (muons), a name related to the main background source (photon conversions) for leptons from heavy-flavour decays. A strong suppression of non-photonic electrons, of the same order of magnitude as the one seen for light hadrons, was observed by PHENIX [263] and STAR [264], although these measurements could not separate the contributions from bottom and charm quarks. More recently, direct measurements of the hadronic decay of D mesons by STAR [265], opened the way towards a more detailed understanding of heavy-quark energy loss. Finally, thanks to the recent experiment upgrades (Si vertex tracking systems), PHENIX is now able to separate charm/bottom contributions in the lepton decay channel [266].

At the LHC, thanks to the much increased heavy-quark production cross section, the indirect measurements have been extended towards higher pT (ATLAS [267], ALICE [268-270]), a detailed study of hadronic decays for various charmed mesons has taken place (ALICE [271-276]), and results on beauty production via either detection of J/-0 from the decay B ^ J/^ + X (ALICE [277], CMS [278]) or hadronic decays (CMS [279]) are now available. In addition, the study of jets originated from a b-quark has now been extended for the first time from pp to heavy-ion collisions by CMS [280,281].

In fig. 21 (left) we show the Raa for D mesons measured in PbPb collisions at ^JsNN = 2.76 TeV by ALICE [271], in the region 8 <pT < 16 GeV/c, as a function of the centrality. The result is obtained as the average of D0, D+ and D*'+ Raa values. It is compared with the corresponding quantity for charged pions and to a calculation implementing energy loss for gluons, light and heavy quarks, including both radiative and collisional processes [282]. In fig. 21 (right) the result for D mesons is compared to CMS data for J/-0 from b-decays [278], in the pT range 6.5 <pT < 30 GeV/c, and to the same theoretical calculation. It can be seen that the D meson suppression reaches a factor 5 for central (0-10%) PbPb events, and is compatible with the corresponding result for charged pions. The expected smaller suppression for particles resulting from the fragmentation of heavy quarks is not seen neither in the data nor in the model calculations, and is interpreted (mainly) as a consequence of the different fragmentation function and pT spectrum of gluons, light and c quarks [262,283]. On the contrary, when comparing b- and c-quark suppression [271], a clear ordering is seen, with the heavier b quark being less suppressed. It has to be noted that the effect is even more evident when the ALICE results of fig. 21 are compared with preliminary higher-statistics results from CMS [284].

Further insight into the interaction processes of heavy quarks in the medium can be obtained by studying the pT-dependence of Raa and v2, as the two quantities are sensitive to the interplay of the various energy loss mechanisms. In fig. 22 we show (left panel) ALICE results on the average D-meson v2 measured in PbPb collisions in the 30-50% centrality class, and (right panel) the Raa for the 0-20% centrality class [274]. The results are compared with various model calculations [285-291]. A significant (> 5a, for pT integrated results) non-zero v2 is observed, an effect that can

0.4 0.3 0.2

ALICE D ,D+, D + average Pb-Pb, \fsNN = 2.76 TeV : ] Syst. from data ] Syst. from B feed-down

Centrality 30-50%

/ÍRn.,

.1- ^ f n -: láff&SSStmagj-*.....«HrtBra.iri«iiKiiiii.¿

...... WHDG rad+coll

POWLANG Cao, Qin, Bass

MCasHQ+EPOS, Coll+Rad(LPM) ----TAMU elastic

BAMPS --UrQMD

10 12 14 16

PT (GeV/c)

0.8 0.6 0.4 0.2 0

¡MSW ~

Pb-Pb, |/sNN = 2.76 TeV Centrality 0-20% WHDG rad+coll POWLANG

-------- Cao, Qin, Bass

MCasHQ+EPOS, Coll+Rad(LPM) BAMPS TAMU elastic UrQMD

'^S^sw«*«- ~ «i i^t-:T ■.'.......».v.v.v,v

^'-^r^:^.--..,-------¿t-""-' ^'

V.yulHW'.

10 12 14 16 PT (GeV/c)

Fig. 22. Left: the average D-meson v2 in the 30-50% centrality class, measured by ALICE [274]. Right: the average D-meson Raa on the 0-20% centrality class, measured by ALICE [274]. Both results refer to PbPb collisions at \JSnn = 2.76 TeV, and are compared to theoretical models.

GC .ó

0.8 0.6 0.4 0.2 0

1 ,—r

0-100%

^ = 2.76 TeV

PbPb, 150 p.b-1 pp, 5.3 pb-1

pQCD: PLB 726 (2013) 251-256

gmed = 2.0 gmed = 2.2

b-jet p (GeV/c)

£ < cl a.

pPb 35 nb-1 (5.02 TeV)

b jet Rpythia, -2.5 < n < 1.5

J pA CM

Huang, Kang, Vitev (Ref. [29]) -

pPb luminosity uncert. pp reference uncert.

0 50 100 150 200 250 300 350 400 jet pT [GeV/c]

Fig. 23. Left: the centrality-integrated 6-jet RAA as a function of pT, from CMS [280]. The data are compared to pQCD-based calculations from [293]. Right: the 6-jet RpA as a function of jet pT [281]. Filled boxes represent the systematic uncertainties. The pp reference and integrated luminosity uncertainties are shown as red and green bands around unity, respectively. A pQCD prediction from Huang et al. [293] is also shown.

gmed = 1

be related, at large pT, to the difference in the in-medium path length between D mesons emitted in-plane or out-of-plane. At low pT in-medium interactions of the produced heavy-quarks can also lead to a momentum anisotropy and consequently to a non-zero v2. Concerning Raa, a strong suppression, reaching a factor 3-4, is seen for pT > 5 GeV/c, with a tendency towards a smaller suppression for decreasing pT. The comparison with the models shows that it is still challenging to describe quantitatively both observables over the full pT range. Finally, it has to be noted that cold nuclear matter effects, and in particular parton shadowing in the nucleus [28, 29], cannot account for the observed D-meson suppression, since the corresponding measurements for pPb collisions give nuclear modification factor values compatible with 1 [292] in the whole pT range explored (pT > 1 GeV/c).

The study of the behaviour of heavy quarks in the medium has been extended to very large pT by studying b-jet production. Figure 23 (left) shows the centrality-integrated Raa measured by CMS [280]. A significant suppression is observed, qualitatively consistent with that of inclusive jets [224], showing that at high pT a large mass/flavour dependence for parton energy loss is unlikely, as expected by theoretical models for large parton energies. Recent CMS results on b-jet production in pPb [281] (see fig. 23 (right)) show a nuclear modification factor compatible with unity, indicating that the suppression observed in PbPb is not due to cold-nuclear-matter effects.

Fig. 24. Left: the inclusive J/^ RAA as a function of (Npart), measured in PbPb collisions at forward-^ by ALICE [312] and compared to PHENIX Au-Au results. Right: the ^-dependence of the inclusive J/^ RAA [312], compared to PHENIX results and to theoretical models.

3.4 Quarkonium

The study of bound states of heavy quarks, charmonia (cc) and bottomonia (66), plays an important role in the understanding of the properties of the medium created in nuclear collisions. A suppression of quarkonium production, related to the screening of the quark-antiquark binding in a colour-deconfined medium, was predicted long ago as a signature of the formation of a QGP [294]. In addition, states with increasing binding energies were predicted to melt for increasing temperatures of the medium [295]. The suppression of the charmonium states J/-0 and ^(2S) was indeed observed at SPS [296-298] and RHIC [299, 300] energies, although the interpretation of the results proved to be not straightforward, due to competing mechanisms which also reduce the charmonium yield but are not related to deconfinement. These include the break-up of the charmonia by interaction with cold nuclear matter and with the hot hadronic gas produced when the medium cools down below the deconfinement temperature [301-303]. Also initial state effects as nuclear shadowing [28, 29] have been proven to play a sizeable role on the quarkonium production rates. At the LHC, new features come into play. On the charmonium side, as a consequence of the large multiplicity of cc pairs (> 100 for central PbPb collisions), recombination mechanisms in the deconfined phase and/or at hadronization were predicted to counterbalance the suppression effects [304,305]. In addition, bottomonium states are more copiously produced at the LHC with respect to lower-energy machines, and an accurate study of the strongly bound Y (1S) state, along with the corresponding 2S and 3S resonances, becomes possible.

Heavy-quarkonium states in pPb and PbPb collisions were mainly studied by ALICE [277, 306-316] and CMS [278, 317-320], with contributions on selected topics also from ATLAS [321,322] and LHCb [323,324]. The experiments measure quarkonia in the dilepton decay mode (mainly , but also e+e-), both at central rapidity (CMS, with

maximum coverage corresponding to \y\ < 2.4) and forward rapidity (ALICE, 2.5 < y < 4). ALICE covers charmonium production down to zero pT while CMS acceptance starts at pT ~ 6 GeV/ c (down to a minimum of 3 GeV/ c for selected rapidity intervals). For bottomonium, all the experiments have coverage down to pT = 0.

In fig. 24 (left) the dependence of the inclusive J/-0 Raa on (Npart), measured by ALICE at forward rapidity [312], is compared to the corresponding measurement by PHENIX (Au-Au, 1.2 < \y\ < 2.2) [299]. For (Npart) > 100 the suppression observed in ALICE data is clearly smaller than at RHIC energy. Such a behaviour, also seen at midrapidity [308], is compatible with the presence of a significant regeneration of J/^. This mechanism is expected to enhance charmonium production mainly at low transverse momentum. When comparing the pT-dependence of Raa for central collisions between PHENIX and ALICE, as in fig. 24 (right) [312], one can clearly see that the difference is concentrated at low pT. Theoretical calculations (including transport and statistical models) qualitatively reproduce the ALICE data [325-328]. For transport models, which implement charmonium dissociation and regeneration in a thermally expanding fireball, more than 50% of the measured J/-0 is produced via regeneration. As for many observables, it is expected that cold-nuclear-matter effects may influence the measured Raa. Estimates based on the results from pPb collisions at yjsNN = 5.02 TeV show that cold-nuclear-matter effects are negligible at large pT, while at low pT their size becomes of the same order of magnitude as the observed suppression. Once such effects are corrected for, the low-pT enhancement of the J/^ production becomes even stronger [311].

In fig. 25 (left) preliminary results on the nuclear modification factor for Y(1S) and Y(2S), measured by CMS in \y\ < 2.4, are shown as a function of Npart [329]. A strong suppression, increasing with centrality, is seen, in particular for the relatively less bound 2S state. Such a behaviour is consistent with the observation of the sequential suppression of the quarkonium states according to their binding energy [295]. One should note that a significant fraction (originally estimated to be ~ 50% by CDF [330], but likely to be ~ 30% following more recent LHCb results [331]) of

PbPb 166 lib"', pp 5.4 pb"1

l(SNN = 2.76 TeV

0.8 O.67 0.4^ 0.2^ 0,

|y| < 2.4

Preliminary

; t(IS)

• T(2S)

0 50 100 150 200 250 300 350 400

1.4 1.2 1

0.8 0.6 0.4 0.2 0

PbPb 166 lb"1, pp 5.4 pb"1

1 Cent. 0-100%, |y| < 2.4

Preliminary

■ t(1S) • T(2S)

0 2 4 6 8 10 12 14 16 18 20

1.4 1.2 1

0.8 0.6 0.4 0.2 0

PbPb 166 lb"1, pp 5.4 pb"1

(SnN = 2.76 TeV

Ill 1 Cent. 0"100% CMS Preliminary

: ; t(1S)

; T(2S) _

- * 0 t -

"1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1"

40"80%, |n| < 1.37

.......I-►

40"80%, 1.52 < |n| < 2.37

■ JETPHOX "err". —IETPHOX Pb+Pb/pp —IETPHOX EPS09/pp <T AA> uncertainty

20"40%, |T|| < 137

.......I-►

20"40%, 1.52 < |T|| < 2.37

10"20%, |T|| < 1.37

.......I-►

10"20%, 1.52 < |n| < 2.37

0"10%, |n| < 1.37

......I-'

11 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1

ATLAS j,____b—k

0 0.4 0.8 1.2 1.6 2 2.4

Fig. 25. The Y(1S) and Y(2S) Raa as a function of Npart (left), y (center), pt (right), measured by CMS [329].

I 40 ^35 30 25 20 15 10 5 0

= .= * t

• W± Data —W± POWHEG CT10

r W+ i W

Ldt = 0.14"0.15 nb" Pb+Pb|[iNN = 2.76 TeV . I .... I .... I .... I .... I.■.■I .... I .... I .

30 50 100

Photon p [GeV]

30 50 100

Photon p [GeV]

30 50 100

Photon p [GeV]

30 50 100

Photon p [GeV]

0 50 100 150 200 250 300 350 400

< N rt>

x part '

Fig. 26. Left: fully corrected normalized yields of prompt photons by ATLAS as a function of pT in \n\ < 1.37 and 1.52 < |n| < 2.37 using tight photon selection, isolation cone size Riso = 0.3 and isolation transverse energy of less than 6 GeV, divided by JETPHOX predictions for pp collisions, which implement the same isolation selection. The combined scale and PDF uncertainty on the JETPHOX calculation is shown by the yellow area. In addition two other JETPHOX calculations are shown, also divided by the pp results: PbPb collisions with no nuclear modification (solid black line), and PbPb collisions with EPS09 nuclear modifications (blue area). Statistical uncertainties are shown by the bars. Systematic uncertainties on the photon yields are combined and shown by the upper and lower braces, and the uncertainty of the nuclear overlap in grey. Right: W boson production yield per binary collision as a function of the mean number of participants Npart for W+, W- and W± bosons for combined muon and electron channels, in PbPb collisions at 2.76TeV/nucleon from ATLAS. Taken from [336] and [338].

measured Y(1S) states comes from the decay of higher-mass bottomonium resonances (Y(2S), Y(3S), Xb) so that a large fraction of the observed 1S suppression is connected to such feed-down effects. The contribution of regeneration effects is expected to be much smaller in the bottomonium sector, due to the lower bb yields compared to charm. In fig. 25 (center, right) the pT- and y-dependence of Raa are shown [329]. No significant dependence on the two variables can be seen, within uncertainties. It is remarkable that also ALICE results on Y(1S), extending the y-coverage up to y = 4, exhibit the same level of suppression [315]. Theoretical models do not still reproduce quantitatively the size of the suppression and/or its rapidity dependence [332-334].

3.5 Electroweak bosons and large transverse momentum photons

High transverse momentum isolated photons have been measured at the LHC by CMS [335] and ATLAS [336]. Electroweak bosons have been measured in PbPb by ATLAS [321, 337, 338] and CMS [339-341] and in pPb by LHCb [342], CMS [343] and ATLAS [344].

In PbPb collisions, the main use of such measurements is as a control probe of our understanding of the calibration of the initial flux of hard probes that goes through the medium. Such a flux is calculated in perturbative QCD without modification due to the presence of such dense partonic medium. Note that in perturbative QCD we assume that collinear factorisation works, with the scaling between pp and AA collisions given by the number of binary collisions (or the average nuclear thickness) computed in the Glauber model. This is illustrated in fig. 26 where it is shown how the nuclear modification factor of high pT photons is compatible with unity and with perturbative QCD calculations, and how W production at different centralities scales very well with the number of collisions.

p+Pb 2013, L = 29 nb-1

VsNN = 5-02 TeV

Glauber(œa=0)

0-10% Centrality 10-40% Centrality P-40-90% Centrality —CT10(NLO)

■ 0.4 - 0.3; : 0.2; ~ 0.1 0; -0.1; -0.2 -0.3 -0.4

pPb 34.6 nb-'

I I I I I

Pnn = 5.02 TeV

■T JpWJpi Iivy.il ▼

- ___VffA

J*»*?* -

W ^ l + V pT > 25 GeV/c

• Data CT10 AW CT10+EPS09

Fig. 27. Left, top panel: the rapidity differential Z boson yields by ATLAS, scaled by (Ncoii), for three centrality ranges compared with the CT10 model calculation. The bars indicate statistical uncertainty and the shaded boxes systematic uncertainty. The scale uncertainty stemming from the centrality calculation for each bin is included in the systematic uncertainty. The uncertainty associated with the model is not plotted. Left, middle panel: the ratios of the data to the model. The uncertainty of the model added in quadrature with the scale uncertainty due to uncertainty in the inclusive NN cross section is shown as a band around unity. Left, bottom panel: RCP. The 0-10% and 40-90% centrality points are offset for visual clarity. The arrows in the lower panels indicate values outside the plotted axes. Right: lepton charge asymmetry (N+ — N— )/(Nt+ + N—) for W ^ l, vs. lepton

pseudorapidity (with the proton going towards positive pseudorapidities), Taken from [344] and [343].

from CMS, compared to theoretical calculations.

In pPb collisions, the main use of these observables is for constraining the nuclear modification of parton densities [345,346]. In fig. 27 we show some LHC results compared to predictions of perturbative QCD containing the state-of-the-art nuclear modification of parton densities [28]. In the backward (Pb-going) region, some discrepancies between theory and data can be observed, pointing to deficiencies in the current nuclear parton densities, most probably related with the assumption of isospin independence for their nuclear modifications. These data, together with those from jets, will be used for improving the nPDFs.

3.6 Low transverse momentum photons and dileptons

The measurement of the electromagnetic probes (real and virtual photons) plays an important role since these particles do not interact strongly with the medium and carry undistorted information on the various stages of the collision history including the early ones. In particular, low pT thermal photons (or dileptons) from the equilibrated QGP phase can give information on the early temperature of the system. Dielectrons/dimuons, thanks to their coupling with low-mass resonances and in particular with the p-meson, can also be sensitive to the partial restoration of chiral symmetry which may induce variations in the spectral functions of the p and its chiral partner ai [347]. At SPS energy, the most accurate results come from the NA60 experiment which observed a strong broadening of the p peak [348], with no mass shift and an enhancement, in the mass region between the $ and the J/^, which could be ascribed to thermal dimuon production corresponding to a temperature T ~ 190 MeV, above the critical temperature for deconfinement [349]. At RHIC energies, by comparing the measured invariant mass dilepton spectrum with the expected yield from hadronic decays (commonly referred to as the "hadronic cocktail"), a clear excess has been singled out in the low-mass region below the p and w. Results from STAR [350] are well reproduced by models [351, 352] which were able to describe the NA60 data. Very recent results from PHENIX [353] are consistent with STAR measurements. By treating the PHENIX low-mass excess below me+e- = 0.03 GeV/c2 as photon internal conversions [354], the slope of their transverse momentum distribution is compatible with T ~ 220 MeV and can be reproduced by hydrodynamical models [355] with initial temperatures in the range 300-600 MeV.

PT (GeV/c)

Fig. 28. Direct photon spectrum measured by ALICE for PbPb collisions at ^Jsnn = 2.76 TeV [358]. Comparisons with NLO pQCD calculations are shown.

At the LHC, such measurements are made much more difficult by the large background levels, mainly connected with the large increase in the combinatorial decays of charged pions/kaons and heavy flavour for dilepton measurements, and from n0 decays for real photons. For the moment, analyses of the dilepton invariant mass spectra have been carried out by ALICE in pp and pPb collisions, and the results have been compared with hadronic cocktail calculations [356]. Within uncertainties, an agreement is observed. In pp collisions, an extraction of the direct photon yield starting from the dilepton spectrum in the low-mass region, gives results consistent with NLO pQCD calculations [356]. In PbPb collisions, no results have been provided up to now, and a significant improvement of the experimental conditions is expected for Run 3, due to the foreseen upgrades of the ALICE Inner Tracking System and Time Projection Chamber [357]. Direct photon studies have been carried out by ALICE by using photon conversions to e+e~ or a calorimetric measurement [358]. The direct photon yield 7direct was extracted subtracting the decay contribution 7decay from the inclusive spectrum 7inci as 7direct = 7inci x (1 — Ydecay/7inci). NLO pQCD calculations [359-362] well describe the spectrum above pT ~ 5GeV/c, while in the region 0.9 < pT < 2.1 GeV/c a 2.6a excess is found for the 0-20% centrality class (see fig. 28). It is compatible with a thermal slope with T = 297 ± 12(stat) ± 41(syst) MeV. Preliminary results [363] show a non-zero v2 for direct photons, possibly reflecting the development of collective expansion at early stages.

3.7 Ultraperipheral collisions

Ultraperipheral collisions [364], usually analysed for exclusive final states, offer unique possibilities for understanding the production mechanisms and for determining parton densities in a hadron collider that otherwise could only be studied in lepton-hadron/nucleus machines. In nucleus-nucleus collisions, one or both nuclei act as quasireal photon sources that collide with the other nucleus or photon. Therefore, 7A collisions can be studied, in which the nucleus remains intact (coherent production) or gets excited and dissolves but without filling the rapidity gap between the nucleus and the produced particle (incoherent production).

ALICE has measured exclusive coherent and incoherent charmonium production [365-367] and exclusive coherent p0 production [368], see fig. 29. Data show a large discriminatory power on models6 and, in the case of charmonium, have been taken as a first direct evidence of the existence of gluon shadowing in nuclei [369]. They have also been used to test dipole models that describe data in lepton-hadron collisions [370,371].

Finally, the photon flux from Pb has been used in pPb [372] to study the energy dependence of J/-0 production in Yp collisions with results competitive with those obtained through pp [373, 374].

6 Due to lack of space we cannot provide a description of the different models. Let us simply note that STARLIGHT, a simulator widely used by the experimental collaborations, is a model that contains no shadowing of nPDFs.

Fig. 29. Left: cross section for coherent exclusive Jfrf production vs. rapidity in ultraperipheral PbPb collisions at 2.76 TeV/ nucleon from ALICE, compared to theoretical calculations. Right: cross section for coherent exclusive p0 production versus rapidity in ultraperipheral PbPb collisions at 2.76 TeV/nucleon from ALICE, compared to theoretical calculations. Taken from [366] and [368].

4 Summary

In this manuscript we have overviewed the main results obtained in pPb and PbPb collisions at the Large Hadron Collider during Run 1. We have first discussed the results for soft probes of the medium created in the collisions and then turned to hard and electroweak probes. We have focused on a selection of experimental results and, due to space limitations, the main current theoretical explanations have just been briefly presented.

Overall, the physical picture of the medium created in PbPb collisions that arises from the results in Run 1 is qualitatively similar to that obtained at RHIC: a system with a large energy density and strong collectivity, behaving very early close to a nearly ideal liquid and very opaque to energetic coloured particles that traverse it. But the larger collision energy provided by the LHC has led to a longer-lived medium with a larger energy density and to a more abundant rate of hard probes. In this situation, and thanks to the sophisticated detector systems and new experimental techniques, many of the observations at RHIC have been made again at the higher energy, with substantially increased precision. This is particularly true for azimuthal asymmetries where the higher energy density leads to higher multiplicities which, in association with the much larger detector acceptances, provides dramatically improved measurements of anisotropies. In particular, spectacular progress has been carried out in determining them on an event-by-event basis. Also for hard probes, where particles with higher masses and transverse energies can be precisely measured, new density effects appear that affect the production and in-medium dynamics of quarkonia, and new observables can be analysed due to the excellent performance of the detectors.

This new era of precision is triggering numerous theoretical efforts. First, data on soft and electromagnetic probes are stimulating new developments on i) the description of the initial state [375], ii) the understanding of the weak or strong dynamics leading to a stage that is isotropised and thermalised enough for relativistic hydrodynamics to be applicable [18,23], and iii) the relation between transport and hydrodynamics and the extraction of the transport coefficients characterising the medium [20]. Second, measurements on hard probes are leading to new developments on our understanding of i) energy loss processes for both light and heavy partons in order to get a unified description of particle and jet inclusive and differential observables [19, 22, 24], and ii) the formation of quarkonium bound states and the different physical mechanisms through which the presence of a coloured medium may affect them [21, 376]. All these elements are required for a quantitative characterisation of the new phase of matter produced in the collisions with controlled uncertainties. They go in parallel with new developments in lattice and perturbative QCD at finite temperature, and with our understanding of the benchmark for soft and hard probes.

Besides, the pPb runs have brought the bonus of the apparent collective behaviour of such collisions, that may also be at play in high-multiplicity pp ones. These findings raise new questions about the nature, thermal or not, and onset of such collectivity, together with new opportunities for benchmarking for hard processes like the determination of nPDFs, and for understanding of QCD at large partonic densities. The latter point naturally links the studies at the LHC with the proposals of future electron-ion colliders [377, 378].

With the increasing statistics and larger energy that will be achieved in Run 2 and, further in the future, with the detector upgrades [379-382] that will become operative in Run 3 and afterwards, heavy-ion physics will continue to play a very important role in the LHC programme, leading to a complete and quantitative comprehension of the high-temperature QGP phase of QCD matter.

We thank John Jowett for useful discussions and for providing fig. 1, and Anton Andronic, Elena Ferreiro, Matthew Luzum, Andreas Morsch, Matthew Nguyen and Peter Steinberg for a critical reading of the manuscript and most useful suggestions. The work of NA was supported by the People Programme (Marie Curie Actions) of the European Union Seventh Framework Programme FP7/2007-2013/ under REA grant agreement No. 318921, the European Research Council grant HotLHC ERC-2011-StG-279579, Ministerio de Ciencia e Innovacion of Spain under project FPA2014-58293-C2-1-P, Xunta de Galicia (Consellerla de Educacion and Conselleria de Innovacion e Industria - Programa Incite), the Spanish Consolider-Ingenio 2010 Programme CPAN and FEDER.

Open Access This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

References

1. H. Fritzsch, M. Gell-Mann, H. Leutwyler, Phys. Lett. B 47, 365 (1973).

2. B.L. Ioffe, V.S. Fadin, L.N. Lipatov, Quantum Chrom,odynam,ics: Perturbative and Nonperturbative Aspects (Cambridge University Press, 2010).

3. J. Collins, Foundations of Perturbative QCD (Cambridge University Press, 2013).

4. D.J. Gross, F. Wilczek, Phys. Rev. Lett. 30, 1343 (1973).

5. H.D. Politzer, Phys. Rev. Lett. 30, 1346 (1973).

6. J. Greensite, Lect. Notes Phys. 821, 1 (2011).

7. Z. Fodor, C. Hoelbling, Rev. Mod. Phys. 84, 449 (2012) (arXiv:1203.4789 [hep-lat]).

8. T.D. Lee, G.C. Wick, Phys. Rev. D 9, 2291 (1974).

9. J.C. Collins, M.J. Perry, Phys. Rev. Lett. 34, 1353 (1975).

10. J. Hofmann, H. Stocker, W. Scheid, W. Greiner, Report of the Workshop on BeV/nucleon Collisions of Heavy Ions: How and Why, Bear Mountain, New York, 29 November-1 December 1974.

11. B.A. Freedman, L.D. McLerran, Phys. Rev. D 16, 1169 (1977).

12. N. Cabibbo, G. Parisi, Phys. Lett. B 59, 67 (1975).

13. E.V. Shuryak, Phys. Lett. B 78, 150 (1978) (Sov. J. Nucl. Phys. 28, 408 (1978)).

14. R. Hwa (Editor), Quark-Gluon Plasma, Vol. 1 (World Scientific, Singapore, 1990).

15. R. Hwa (Editor), Quark-Gluon Plasma, Vol. 2 (World Scientific, Singapore, 1995).

16. R. Hwa, X.N. Wang (Editors), Quark-Gluon Plasma, Vol. 3 (World Scientific, Singapore, 2004).

17. R. Hwa, X.N. Wang (Editors), Quark-Gluon Plasma, Vol. 4 (World Scientific, Singapore, 2010).

18. P.M. Chesler, W. van der Schee, arXiv:1501.04952 [nucl-th].

19. J. Ghiglieri, D. Teaney, arXiv:1502.03730 [hep-ph].

20. S. Jeon, U. Heinz, arXiv:1503.03931 [hep-ph].

21. H.T. Ding, F. Karsch, S. Mukherjee, arXiv:1504.05274 [hep-lat].

22. J.P. Blaizot, Y. Mehtar-Tani, arXiv:1503.05958 [hep-ph].

23. F. Gelis, arXiv:1508.07974 [hep-ph].

24. G.Y. Qin, X.N. Wang, arXiv:1511.00790 [hep-ph].

25. M. Gyulassy, L. McLerran, Nucl. Phys. A 750, 30 (2005) (nucl-th/0405013).

26. B. Muller, J.L. Nagle, Annu. Rev. Nucl. Part. Sci. 56, 93 (2006) (nucl-th/0602029).

27. F. Karsch, Lect. Notes Phys. 583, 209 (2002) (hep-lat/0106019).

28. K.J. Eskola, H. Paukkunen, C.A. Salgado, JHEP 04, 065 (2009) (arXiv:0902.4154 [hep-ph]).

29. D. de Florian, R. Sassot, P. Zurita, M. Stratmann, Phys. Rev. D 85, 074028 (2012) (arXiv:1112.6324 [hep-ph]).

30. N. Armesto, N. Borghini, S. Jeon, U.A. Wiedemann, S. Abreu, V. Akkelin, J. Alam, J.L. Albacete et al., J. Phys. G 35, 054001 (2008) (arXiv:0711.0974 [hep-ph]).

31. N. Armesto, arXiv:0903.1330 [hep-ph].

32. C.A. Salgado, J. Alvarez-Muniz, F. Arleo, N. Armesto, M. Botje, M. Cacciari, J. Campbell, C. Carli et al., J. Phys. G 39, 015010 (2012) (arXiv:1105.3919 [hep-ph]).

33. J.L. Albacete, N. Armesto, R. Baier, G.G. Barnafoldi, J. Barrette, S. De, W.-T. Deng, A. Dumitru et al., Int. J. Mod. Phys. E 22, 1330007 (2013) (arXiv:1301.3395 [hep-ph]).

34. W. Fischer, J.M. Jowett, Rev. Accel. Sci. Technol. 7, 49 (2014).

35. B. Muller, J. Schukraft, B. Wyslouch, Annu. Rev. Nucl. Part. Sci. 62, 361 (2012) (arXiv:1202.3233 [hep-ex]).

36. G. Roland, K. Safarik, P. Steinberg, Prog. Part. Nucl. Phys. 77, 70 (2014).

37. E. Norbeck, K. Safarik, P.A. Steinberg, Annu. Rev. Nucl. Part. Sci. 64, 383 (2014).

38. ALICE Collaboration (B.B. Abelev et al.), JINST 9, P11003 (2014) (arXiv:1405.1849 [nucl-ex]).

39. ALICE Collaboration (b. Abelev et al.), Phys. Rev. Lett. 110, 032301 (2013) (arXiv:1210.3615 [nucl-ex]).

40. CMS Collaboration (V. Khachatryan et al.), arXiv:1509.03893 [hep-ex].

41. TOTEM Collaboration (G. Antchev et al.), EPL 101, 21004 (2013).

42. ALICE Collaboration (B. Abelev et al.), Phys. Rev. Lett. 109, 252302 (2012) (arXiv:1203.2436 [nucl-ex]).

43. ALICE Collaboration (B. Abelev et al.), Phys. Rev. C 88, 044909 (2013) (arXiv:1301.4361 [nucl-ex]).

44. ALICE Collaboration (k. Aamodt et al.), Phys. Rev. Lett. 105, 252301 (2010) (arXiv:1011.3916 [nucl-ex]).

45. ATLAS Collaboration (G. Aad et al.), Phys. Lett. B 710, 363 (2012) (arXiv:1108.6027 [hep-ex]).

46. CMS Collaboration (S. Chatrchyan et al.), JHEP 08, 141 (2011) (arXiv:1107.4800 [nucl-ex]).

47. P. Carruthers, M. Doung-van, Phys. Rev. D 8, 859 (1973).

48. W. Busza, J. Phys. G 35, 044040 (2008) (arXiv:0710.2293 [nucl-ex]).

49. PHOBOS Collaboration (B.B. Back et al.), Phys. Rev. C 74, 021901 (2006) (nucl-ex/0509034).

50. PHENIX Collaboration (S.S. Adler et al.), Phys. Rev. C 71, 034908 (2005) 71, 049901(E) (2005) (nucl-ex/0409015).

51. CMS Collaboration (S. Chatrchyan et al.), Phys. Rev. Lett. 109, 152303 (2012) (arXiv:1205.2488 [nucl-ex]).

52. J.D. Bjorken, Phys. Rev. D 27, 140 (1983).

53. BRAHMS Collaboration (I.G. Bearden et al.), Phys. Lett. B 523, 227 (2001) (nucl-ex/0108016).

54. BRAHMS Collaboration (i.G. Bearden et al.), Phys. Rev. Lett. 88, 202301 (2002) (nucl-ex/0112001).

55. PHOBOS Collaboration (B.B. Back et al.), Phys. Rev. C 65, 061901 (2002) (nucl-ex/0201005).

56. STAR Collaboration (J. Adams et al.), Phys. Rev. C 70, 054907 (2004) (nucl-ex/0407003).

57. ALICE Collaboration (K. Aamodt et al.), Phys. Rev. Lett. 106, 032301 (2011) (arXiv:1012.1657 [nucl-ex]).

58. PHOBOS Collaboration (B.B. Back et al.), Phys. Rev. C 70, 021902 (2004) (nucl-ex/0405027).

59. F.W. Bopp, R. Engel, J. Ranft, S. Roesler, arXiv:0706.3875 [hep-ph].

60. W.T. Deng, X.N. Wang, R. Xu, Phys. Lett. B 701, 133 (2011) (arXiv:1011.5907 [nucl-th]).

61. N. Armesto, C.A. Salgado, U.A. Wiedemann, Phys. Rev. Lett. 94, 022002 (2005) (hep-ph/0407018).

62. D. Kharzeev, E. Levin, M. Nardi, Nucl. Phys. A 747, 609 (2005) (hep-ph/0408050).

63. J.L. Albacete, A. Dumitru, arXiv:1011.5161 [hep-ph].

64. M. Mitrovski, T. Schuster, G. Graf, H. Petersen, M. Bleicher, Phys. Rev. C 79, 044901 (2009) (arXiv:0812.2041 [hep-ph]).

65. Z.W. Lin, C.M. Ko, B.A. Li, B. Zhang, S. Pal, Phys. Rev. C 72, 064901 (2005) (nucl-th/0411110).

66. J. Xu, C.M. Ko, Phys. Rev. C 83, 034904 (2011) (arXiv:1101.2231 [nucl-th]).

67. ALICE Collaboration (E. Abbas et al.), Phys. Lett. B 726, 610 (2013) (arXiv:1304.0347 [nucl-ex]).

68. ALICE Collaboration (j. Adam et al.), arXiv:1509.07299 [nucl-ex].

69. PHOBOS Collaboration (B. Alver et al.), Phys. Rev. C 83, 024913 (2011) (arXiv:1011.1940 [nucl-ex]).

70. D.A. Teaney, arXiv:0905.2433 [nucl-th].

71. P. Romatschke, Int. J. Mod. Phys. E 19, 1 (2010) (arXiv:0902.3663 [hep-ph]).

72. S. Voloshin, Y. Zhang, Z. Phys. C 70, 665 (1996) (hep-ph/9407282).

73. STAR Collaboration (K.H. Ackermann et al.), Phys. Rev. Lett. 86, 402 (2001) (nucl-ex/0009011).

74. BRAHMS Collaboration (I. Arsene et al.), Nucl. Phys. A 757, 1 (2005) (nucl-ex/0410020).

75. B.B. Back et al., Nucl. Phys. A 757, 28 (2005) (nucl-ex/0410022).

76. STAR Collaboration (J. Adams et al.), Nucl. Phys. A 757, 102 (2005) (nucl-ex/0501009).

77. PHENIX Collaboration (K. Adcox et al.), Nucl. Phys. A 757, 184 (2005) (nucl-ex/0410003).

78. STAR Collaboration (C. Adler et al.), Phys. Rev. Lett. 87, 182301 (2001) (nucl-ex/0107003).

79. P. Huovinen, Nucl. Phys. A 761, 296 (2005) (nucl-th/0505036).

80. ALICE Collaboration (K. Aamodt et al.), Phys. Rev. Lett. 105, 252302 (2010) (arXiv:1011.3914 [nucl-ex]).

81. S.A. Voloshin, A.M. Poskanzer, R. Snellings, in Landolt-Boernstein, Relativistic Heavy Ion Physics, Vol. 1/23 (SpringerVerlag, 2010) pp. 5-54.

82. FOPI Collaboration (A. Andronic et al.), Phys. Lett. B 612, 173 (2005) (nucl-ex/0411024).

83. ATLAS Collaboration (G. Aad et al.), Phys. Lett. B 707, 330 (2012) (arXiv:1108.6018 [hep-ex]).

84. CMS Collaboration (S. Chatrchyan et al.), Phys. Rev. Lett. 109, 022301 (2012) (arXiv:1204.1850 [nucl-ex]).

85. STAR Collaboration (L. Adamczyk et al.), Phys. Rev. C 86, 054908 (2012) (arXiv:1206.5528 [nucl-ex]).

86. ATLAS Collaboration (G. Aad et al.), Phys. Rev. C 86, 014907 (2012) (arXiv:1203.3087 [hep-ex]).

87. ALICE Collaboration (K. Aamodt et al.), Phys. Rev. Lett. 107, 032301 (2011) (arXiv:1105.3865 [nucl-ex]).

88. CMS Collaboration (S. Chatrchyan et al.), Phys. Rev. C 89, 044906 (2014) (arXiv:1310.8651 [nucl-ex]).

89. ATLAS Collaboration (G. Aad et al.), JHEP 11, 183 (2013) (arXiv:1305.2942 [hep-ex]).

90. M.L. Miller, K. Reygers, S.J. Sanders, P. Steinberg, Annu. Rev. Nucl. Part. Sci. 57, 205 (2007) (nucl-ex/0701025).

91. H.J. Drescher, A. Dumitru, A. Hayashigaki, Y. Nara, Phys. Rev. C 74, 044905 (2006) (nucl-th/0605012).

92. ALICE Collaboration (B.B. Abelev et al.), JHEP 06, 190 (2015) (arXiv:1405.4632 [nucl-ex]).

93. H. Niemi, K.J. Eskola, R. Paatelainen, arXiv:1505.02677 [hep-ph].

94. B. Schenke, P. Tribedy, R. Venugopalan, Nucl. Phys. A 926, 102 (2014) (arXiv:1312.5588 [hep-ph]).

95. H. Song, S.A. Bass, U. Heinz, T. Hirano, C. Shen, Phys. Rev. Lett. 106, 192301 (2011) 109, 139904(E) (2012) (arXiv:1011.2783 [nucl-th]).

96. P. Kovtun, D.T. Son, A.O. Starinets, Phys. Rev. Lett. 94, 111601 (2005) (hep-th/0405231).

97. ALICE Collaboration (B.B. Abelev et al.), Phys. Rev. C 90, 054901 (2014) (arXiv:1406.2474 [nucl-ex]).

98. ATLAS Collaboration (G. Aad et al.), Phys. Rev. C 90, 044906 (2014) (arXiv:1409.1792 [hep-ex]).

99. CMS Collaboration (S. Chatrchyan et al.), Phys. Lett. B 724, 213 (2013) (arXiv:1305.0609 [nucl-ex]).

100. G. Bafar, D. Teaney, Phys. Rev. C 90, 054903 (2014) (arXiv:1312.6770 [nucl-th]).

101. S. Jeon, V. Koch, in Quark Gluon Plasma,, edited by R.C. Hwa et al., pp. 430-490 (hep-ph/0304012).

102. ALICE Collaboration (B.B. Abelev et al.), Eur. Phys. J. C 74, 3077 (2014) (arXiv:1407.5530 [nucl-ex]).

103. ALICE Collaboration (b. Abelev et al.), Phys. Rev. Lett. 110, 152301 (2013) (arXiv:1207.6068 [nucl-ex]).

104. K. Fukushima, D.E. Kharzeev, H.J. Warringa, Phys. Rev. D 78, 074033 (2008) (arXiv:0808.3382 [hep-ph]).

105. ALICE Collaboration (B. Abelev et al.), Phys. Rev. Lett. 110, 012301 (2013) (arXiv:1207.0900 [nucl-ex]).

106. S.A. Bass, P. Danielewicz, S. Pratt, Phys. Rev. Lett. 85, 2689 (2000) (nucl-th/0005044).

107. ALICE Collaboration (B. Abelev et al.), Phys. Lett. B 723, 267 (2013) (arXiv:1301.3756 [nucl-ex]).

108. ALICE Collaboration (j. Adam et al.), arXiv:1509.07255 [nucl-ex].

109. A. Dumitru, F. Gelis, L. McLerran, R. Venugopalan, Nucl. Phys. A 810, 91 (2008) (arXiv:0804.3858 [hep-ph]).

110. PHOBOS Collaboration (B. Alver et al.), Phys. Rev. Lett. 104, 062301 (2010) (arXiv:0903.2811 [nucl-ex]).

111. STAR Collaboration (B.I. Abelev et al.), Phys. Rev. C 80, 064912 (2009) (arXiv:0909.0191 [nucl-ex]).

112. CMS Collaboration (V. Khachatryan et al.), JHEP 09, 091 (2010) (arXiv:1009.4122 [hep-ex]).

113. CMS Collaboration (v. Khachatryan et al.), arXiv:1510.03068 [nucl-ex].

114. ATLAS Collaboration (G. Aad et al.), arXiv:1509.04776 [hep-ex].

115. ALICE Collaboration (K. Aamodt et al.), Phys. Lett. B 708, 249 (2012) (arXiv:1109.2501 [nucl-ex]).

116. CMS Collaboration (S. Chatrchyan et al.), JHEP 07, 076 (2011) (arXiv:1105.2438 [nucl-ex]).

117. CMS Collaboration (s. Chatrchyan et al.), Eur. Phys. J. C 72, 10052 (2012) (arXiv:1201.3158 [nucl-ex]).

118. CMS Collaboration (s. Chatrchyan et al.), Phys. Lett. B 718, 795 (2013) (arXiv:1210.5482 [nucl-ex]).

119. ALICE Collaboration (B. Abelev et al.), Phys. Lett. B 719, 29 (2013) (arXiv:1212.2001 [nucl-ex]).

120. ATLAS Collaboration (G. Aad et al.), Phys. Rev. Lett. 110, 182302 (2013) (arXiv:1212.5198 [hep-ex]).

121. ALICE Collaboration (J. Adam et al.), arXiv:1506.08032 [nucl-ex].

122. ALICE Collaboration (B.B. Abelev et al.), Phys. Lett. B 741, 38 (2015) (arXiv:1406.5463 [nucl-ex]).

123. ALICE Collaboration (B.B. Abelev et al.), Phys. Lett. B 726, 164 (2013) (arXiv:1307.3237 [nucl-ex]).

124. CMS Collaboration (V. Khachatryan et al.), Phys. Lett. B 742, 200 (2015) (arXiv:1409.3392 [nucl-ex]).

125. CMS Collaboration (v. Khachatryan et al.), arXiv:1502.05382 [nucl-ex].

126. C.B. Chiu, R.C. Hwa, C.B. Yang, Phys. Rev. C 78, 044903 (2008) (arXiv:0801.2183 [nucl-th]).

127. J.D. Bjorken, S.J. Brodsky, A. Scharff Goldhaber, Phys. Lett. B 726, 344 (2013) (arXiv:1308.1435 [hep-ph]).

128. E. Shuryak, I. Zahed, Phys. Rev. D 89, 094001 (2014)(arXiv:1311.0836 [hep-ph]).

129. C. Andres, A. Moscoso, C. Pajares, Phys. Rev. C 90, 054902 (2014) (arXiv:1405.3632 [hep-ph]).

130. P. Bozek, W. Broniowski, Phys. Lett. B 718, 1557 (2013) (arXiv:1211.0845 [nucl-th]).

131. E. Shuryak, I. Zahed, Phys. Rev. C 88, 044915 (2013) (arXiv:1301.4470 [hep-ph]).

132. A. Bzdak, B. Schenke, P. Tribedy, R. Venugopalan, Phys. Rev. C 87, 064906 (2013) (arXiv:1304.3403 [nucl-th]).

133. K. Werner, I. Karpenko, T. Pierog, Phys. Rev. Lett. 106, 122004 (2011) (arXiv:1011.0375 [hep-ph]).

134. S. Gavin, L. McLerran, G. Moschelli, Phys. Rev. C 79, 051902 (2009) (arXiv:0806.4718 [nucl-th]).

135. A. Kovner, M. Lublinsky, Phys. Rev. D 83, 034017 (2011) (arXiv:1012.3398 [hep-ph]).

136. E. Levin, A.H. Rezaeian, Phys. Rev. D 84, 034031 (2011) (arXiv:1105.3275 [hep-ph]).

137. A. Dumitru, K. Dusling, F. Gelis, J. Jalilian-Marian, T. Lappi, R. Venugopalan, Phys. Lett. B 697, 21 (2011) (arXiv:1009.5295 [hep-ph]).

138. K. Dusling, R. Venugopalan, Phys. Rev. Lett. 108, 262001 (2012) (arXiv:1201.2658 [hep-ph]).

139. K. Dusling, R. Venugopalan, Phys. Rev. D 87, 054014 (2013) (arXiv:1211.3701 [hep-ph]).

140. T. Altinoluk, N. Armesto, G. Beuf, A. Kovner, M. Lublinsky, arXiv:1503.07126 [hep-ph].

141. A. Dumitru, L. McLerran, V. Skokov, Phys. Lett. B 743, 134 (2015) (arXiv:1410.4844 [hep-ph]).

142. L. Cunqueiro, J. Dias de Deus, C. Pajares, Eur. Phys. J. C 65, 423 (2010) (arXiv:0806.0523 [hep-ph]).

143. J. Casalderrey-Solana, U.A. Wiedemann, Phys. Rev. Lett. 104, 102301 (2010) (arXiv:0911.4400 [hep-ph]).

144. P. Koch, B. Muller, J. Rafelski, Phys. Rep. 142, 167 (1986).

145. P. Braun-Munzinger, K. Redlich, J. Stachel, in Quark Gluon Plasma, edited by R.C. Hwa et al., pp. 491-599 (nucl-th/0304013).

146. ALICE Collaboration (B.B. Abelev et al.), Phys. Lett. B 728, 216 (2014) 734, 409(E) (2014) (arXiv:1307.5543 [nucl-ex]).

147. ALICE Collaboration (B.B. Abelev et al.), Phys. Rev. C 91, 024609 (2015) (arXiv:1404.0495 [nucl-ex]).

148. J. Cleymans, I. Kraus, H. Oeschler, K. Redlich, S. Wheaton, Phys. Rev. C 74, 034903 (2006) (hep-ph/0604237).

149. A. Andronic, P. Braun-Munzinger, J. Stachel, Phys. Lett. B 673, 142 (2009) 678, 516(E) (2009) (arXiv:0812.1186 [nucl-th]).

150. J. Stachel, A. Andronic, P. Braun-Munzinger, K. Redlich, J. Phys. Conf. Ser. 509, 012019 (2014) (arXiv:1311.4662 [nucl-th]).

151. A. Ortiz Velasquez, P. Christiansen, E. Cuautle Flores, I. Maldonado Cervantes, G. Paic, Phys. Rev. Lett. 111, 042001 (2013) (arXiv:1303.6326 [hep-ph]).

152. ALICE Collaboration (B.B. Abelev et al.), Phys. Rev. Lett. 111, 222301 (2013) (arXiv:1307.5530 [nucl-ex]).

153. ALICE Collaboration (b. Abelev et al.), Phys. Rev. Lett. 109, 252301 (2012) (arXiv:1208.1974 [hep-ex]).

154. ALICE Collaboration (b. Abelev et al.), Phys. Rev. C 88, 044910 (2013) (arXiv:1303.0737 [hep-ex]).

155. CMS Collaboration (S. Chatrchyan et al.), Eur. Phys. J. C 72, 2164 (2012) (arXiv:1207.4724 [hep-ex]).

156. CMS Collaboration (s. Chatrchyan et al.), Eur. Phys. J. C 74, 2847 (2014) (arXiv:1307.3442 [hep-ex]).

157. ALICE Collaboration (B.B. Abelev et al.), Phys. Lett. B 728, 25 (2014) (arXiv:1307.6796 [nucl-ex]).

158. ALICE Collaboration (j. Adam et al.), arXiv:1506.09206 [nucl-ex].

159. ALICE Collaboration (B.B. Abelev et al.), Eur. Phys. J. C 74, 3108 (2014) (arXiv:1405.3794 [nucl-ex]).

160. ALICE Collaboration (B.B. Abelev et al.), Phys. Lett. B 736, 196 (2014) (arXiv:1401.1250 [nucl-ex]).

161. ALICE Collaboration (j. Adam et al.), arXiv:1506.07287 [nucl-ex].

162. PHENIX Collaboration (A. Adare et al.), Phys. Rev. Lett. 101, 232301 (2008) (arXiv:0801.4020 [nucl-ex]).

163. PHENIX Collaboration (a. Adare et al.), Phys. Rev. C 82, 011902 (2010) (arXiv:1005.4916 [nucl-ex]).

164. STAR Collaboration (G. Agakishiev et al.), Phys. Rev. Lett. 108, 072302 (2012) (arXiv:1110.0579 [nucl-ex]).

165. R.J. Fries, B. Muller, C. Nonaka, S.A. Bass, Phys. Rev. C 68, 044902 (2003) (nucl-th/0306027).

166. S. Sapeta, U.A. Wiedemann, Eur. Phys. J. C 55, 293 (2008) (arXiv:0707.3494 [hep-ph]).

167. ALICE Collaboration (J. Adam et al.), arXiv:1506.07499 [nucl-ex].

168. ALICE Collaboration (j. Adam et al.), arXiv:1506.08453 [nucl-ex].

169. ALICE Collaboration (j. Adam et al.), arXiv:1506.08951 [nucl-ex].

170. ALICE Collaboration (j. Adam et al.), Nat. Phys. (2015) (arXiv:1508.03986 [nucl-ex]).

171. G. Goldhaber, S. Goldhaber, W.Y. Lee, A. Pais, Phys. Rev. 120, 300 (1960).

172. M.A. Lisa, S. Pratt, R. Soltz, U. Wiedemann, Annu. Rev. Nucl. Part. Sci. 55, 357 (2005) (nucl-ex/0505014).

173. S.V. Akkelin, Y.M. Sinyukov, Phys. Lett. B 356, 525 (1995).

174. W. Broniowski, M. Chojnacki, W. Florkowski, A. Kisiel, Phys. Rev. Lett. 101, 022301 (2008) (arXiv:0801.4361 [nucl-th]).

175. I.A. Karpenko, Y.M. Sinyukov, K. Werner, Phys. Rev. C 87, 024914 (2013) (arXiv:1204.5351 [nucl-th]).

176. ALICE Collaboration (K. Aamodt et al.), Phys. Lett. B 696, 328 (2011) (arXiv:1012.4035 [nucl-ex]).

177. ALICE Collaboration (B.B. Abelev et al.), Phys. Rev. C 89, 024911 (2014) (arXiv:1310.7808 [nucl-ex]).

178. ALICE Collaboration (B.B. Abelev et al.), Phys. Lett. B 739, 139 (2014) (arXiv:1404.1194 [nucl-ex]).

179. ALICE Collaboration (j. Adam et al.), Phys. Rev. C 91, 034906 (2015) (arXiv:1502.00559 [nucl-ex]).

180. ALICE Collaboration (j. Adam et al.), arXiv:1506.07884 [nucl-ex].

181. ALICE Collaboration (j. Adam et al.), arXiv:1507.06842 [nucl-ex].

182. E895 Collaboration (M.A. Lisa et al.), Phys. Rev. Lett. 84, 2798 (2000).

183. NA49 Collaboration (C. Alt et al.), Phys. Rev. C 77, 064908 (2008) (arXiv:0709.4507 [nucl-ex]).

184. NA49 Collaboration (s.V. Afanasiev et al.), Phys. Rev. C 66, 054902 (2002) (nucl-ex/0205002).

185. CERES Collaboration (D. Adamova et al.), Nucl. Phys. A 714, 124 (2003) (nucl-ex/0207005).

186. STAR Collaboration (B.I. Abelev et al.), Phys. Rev. C 80, 024905 (2009) (arXiv:0903.1296 [nucl-ex]).

187. PHOBOS Collaboration (B.B. Back et al.), Phys. Rev. C 73, 031901 (2006) (nucl-ex/0409001).

188. B.B. Back et al., Phys. Rev. Lett. 91, 052303 (2003) (nucl-ex/0210015).

189. STAR Collaboration (J. Adams et al.), Phys. Rev. C 71, 044906 (2005) (nucl-ex/0411036).

190. STAR Collaboration (b.I. Abelev et al.), Phys. Rev. C 79, 034909 (2009) (arXiv:0808.2041 [nucl-ex]).

191. P. Bozek, W. Broniowski, Phys. Lett. B 720, 250 (2013) (arXiv:1301.3314 [nucl-th]).

192. A.N. Makhlin, Y.M. Sinyukov, Z. Phys. C 39, 69 (1988).

193. A. Kisiel, M. Galazyn, P. Bozek, Phys. Rev. C 90, 064914 (2014) (arXiv:1409.4571 [nucl-th]).

194. V.M. Shapoval, P. Braun-Munzinger, I.A. Karpenko, Y.M. Sinyukov, Nucl. Phys. A 929, 1 (2014) (arXiv:1404.4501 [hep-ph]).

195. A. Accardi, N. Armesto, M. Botje, S.J. Brodsky, B. Cole, K.J. Eskola, G.I. Fai, L. Frankfurt et al., hep-ph/0308248.

196. A. Accardi, F. Arleo, N. Armesto, R. Baier, D.G. d'Enterria, R.J. Fries, O. Kodolova, I.P. Lokhtin et al., hep-ph/0310274.

197. M. Bedjidian, D. Blaschke, G.T. Bodwin, N. Carrer, B. Cole, P. Crochet, A. Dainese, A. Deandrea et al., hep-ph/0311048.

198. F. Arleo, P. Aurenche, F.W. Bopp, I. Dadic, G. David, H. Delagrange, D.G. d'Enterria, K.J. Eskola et al., hep-ph/0311131.

199. A. Majumder, M. Van Leeuwen, Prog. Part. Nucl. Phys. A 66, 41 (2011) (arXiv:1002.2206 [hep-ph]).

200. Y. Mehtar-Tani, J.G. Milhano, K. Tywoniuk, Int. J. Mod. Phys. A 28, 1340013 (2013) (arXiv:1302.2579 [hep-ph]).

201. PHENIX Collaboration (S.S. Adler et al.), Phys. Rev. C 69, 034910 (2004) (nucl-ex/0308006).

202. STAR Collaboration (J. Adams et al.), Phys. Rev. Lett. 91, 172302 (2003) (nucl-ex/0305015).

203. STAR Collaboration (j. Adams et al.), Phys. Rev. Lett. 91, 072304 (2003) (nucl-ex/0306024).

204. PHENIX Collaboration (S.S. Adler et al.), Phys. Rev. Lett. 97, 052301 (2006) (nucl-ex/0507004).

205. STAR Collaboration (J. Adams et al.), Phys. Rev. Lett. 95, 152301 (2005) (nucl-ex/0501016).

206. ALICE Collaboration (K. Aamodt et al.), Phys. Lett. B 696, 30 (2011) (arXiv:1012.1004 [nucl-ex]).

207. CMS Collaboration (S. Chatrchyan et al.), Eur. Phys. J. C 72, 1945 (2012) (arXiv:1202.2554 [nucl-ex]).

208. ALICE Collaboration (B. Abelev et al.), Phys. Lett. B 720, 52 (2013) (arXiv:1208.2711 [hep-ex]).

209. ATLAS Collaboration (G. Aad et al.), arXiv:1504.04337 [hep-ex].

210. ALICE Collaboration (K. Aamodt et al.), Phys. Rev. Lett. 108, 092301 (2012) (arXiv:1110.0121 [nucl-ex]).

211. JET Collaboration (K.M. Burke et al.), Phys. Rev. C 90, 014909 (2014) (arXiv:1312.5003 [nucl-th]).

212. Z.Q. Liu, H. Zhang, B.W. Zhang, E. Wang, arXiv:1506.02840 [nucl-th].

213. N. Armesto, B. Cole, C. Gale, W.A. Horowitz, P. Jacobs, S. Jeon, M. van Leeuwen, A. Majumder et al., Phys. Rev. C 86, 064904 (2012) (arXiv:1106.1106 [hep-ph]).

214. ALICE Collaboration (B. Abelev et al.), Phys. Rev. Lett. 110, 082302 (2013) (arXiv:1210.4520 [nucl-ex]).

215. ALICE Collaboration (B.B. Abelev et al.), Eur. Phys. J. C 74, 3054 (2014) (arXiv:1405.2737 [nucl-ex]).

216. CMS Collaboration (V. Khachatryan et al.), Eur. Phys. J. C 75, 237 (2015) (arXiv:1502.05387 [nucl-ex]).

217. ATLAS Collaboration, ATLAS-CONF-2014-029.

218. G.P. Salam, Eur. Phys. J. C 67, 637 (2010) (arXiv:0906.1833 [hep-ph]).

219. R. Baier, Y.L. Dokshitzer, A.H. Mueller, D. Schiff, Phys. Rev. C 60, 064902 (1999) (hep-ph/9907267).

220. C.A. Salgado, U.A. Wiedemann, Phys. Rev. Lett. 93, 042301 (2004) (hep-ph/0310079).

221. PHENIX Collaboration (D.V. Perepelitsa), J. Phys.: Conf. Ser. 389, 012006 (2012).

222. STAR Collaboration (L. Adamczyk et al.), Phys. Rev. Lett. 112, 122301 (2014) (arXiv:1302.6184 [nucl-ex]).

223. ATLAS Collaboration (G. Aad et al.), Phys. Rev. Lett. 105, 252303 (2010) (arXiv:1011.6182 [hep-ex]).

224. ATLAS Collaboration (g. Aad et al.), Phys. Lett. B 719, 220 (2013) (arXiv:1208.1967 [hep-ex]).

225. CMS Collaboration (S. Chatrchyan et al.), Phys. Rev. C 84, 024906 (2011) (arXiv:1102.1957 [nucl-ex]).

226. CMS Collaboration (s. Chatrchyan et al.), Phys. Lett. B 712, 176 (2012) (arXiv:1202.5022 [nucl-ex]).

227. CMS Collaboration (s. Chatrchyan et al.), Phys. Lett. B 718, 773 (2013) (arXiv:1205.0206 [nucl-ex]).

228. CMS Collaboration (s. Chatrchyan et al.), JHEP 10, 087 (2012) (arXiv:1205.5872 [nucl-ex]).

229. CMS Collaboration (s. Chatrchyan et al.), Phys. Lett. B 730, 243 (2014) (arXiv:1310.0878 [nucl-ex]).

230. ATLAS Collaboration (G. Aad et al.), Phys. Rev. Lett. 111, 152301 (2013) (arXiv:1306.6469 [hep-ex]).

231. ALICE Collaboration (B. Abelev et al.), JHEP 03, 013 (2014) (arXiv:1311.0633 [nucl-ex]).

232. ALICE Collaboration (j. Adam et al.), Phys. Lett. B 746, 1 (2015) (arXiv:1502.01689 [nucl-ex]).

233. CMS Collaboration (S. Chatrchyan et al.), Phys. Rev. C 90, 024908 (2014) (arXiv:1406.0932 [nucl-ex]).

234. ATLAS Collaboration (G. Aad et al.), Phys. Lett. B 739, 320 (2014) (arXiv:1406.2979 [hep-ex]).

235. ATLAS Collaboration (g. Aad et al.), Phys. Rev. Lett. 114, 072302 (2015) (arXiv:1411.2357 [hep-ex]).

236. ALICE Collaboration (J. Adam et al.), arXiv:1506.03984 [nucl-ex].

237. ATLAS Collaboration (G. Aad et al.), arXiv:1506.08656 [hep-ex].

238. ALICE Collaboration (J. Adam et al.), arXiv:1509.07334 [nucl-ex].

239. CMS Collaboration (V. Khachatryan et al.), arXiv:1509.09029 [nucl-ex].

240. ALICE Collaboration (B. Abelev et al.), JHEP 03, 053 (2012) (arXiv:1201.2423 [hep-ex]).

241. M. Cacciari, G.P. Salam, G. Soyez, Eur. Phys. J. C 71, 1692 (2011) (arXiv:1101.2878 [hep-ph]).

242. J. Casalderrey-Solana, J.G. Milhano, U.A. Wiedemann, J. Phys. G 38, 035006 (2011) (arXiv:1012.0745 [hep-ph]).

243. G.Y. Qin, B. Muller, Phys. Rev. Lett. 106, 162302 (2011) 108, 189904(E) (2012) (arXiv:1012.5280 [hep-ph]).

244. Y. He, I. Vitev, B.W. Zhang, Phys. Lett. B 713, 224 (2012) (arXiv:1105.2566 [hep-ph]).

245. C. Young, B. Schenke, S. Jeon, C. Gale, Phys. Rev. C 84, 024907 (2011) (arXiv:1103.5769 [nucl-th]).

246. I.P. Lokhtin, A.V. Belyaev, A.M. Snigirev, Eur. Phys. J. C 71, 1650 (2011) (arXiv:1103.1853 [hep-ph]).

247. T. Renk, Phys. Rev. C 85, 064908 (2012) (arXiv:1202.4579 [hep-ph]).

248. L. Apolinario, N. Armesto, L. Cunqueiro, JHEP 02, 022 (2013) (arXiv:1211.1161 [hep-ph]).

249. K.C. Zapp, F. Krauss, U.A. Wiedemann, JHEP 03, 080 (2013) (arXiv:1212.1599 [hep-ph]).

250. J. Casalderrey-Solana, Y. Mehtar-Tani, C.A. Salgado, K. Tywoniuk, Phys. Lett. B 725, 357 (2013) (arXiv:1210.7765 [hep-ph]).

251. Y. Mehtar-Tani, K. Tywoniuk, Phys. Lett. B 744, 284 (2015) (arXiv:1401.8293 [hep-ph]).

252. CMS Collaboration (S. Chatrchyan et al.), Eur. Phys. J. C 74, 2951 (2014) (arXiv:1401.4433 [nucl-ex]).

253. ATLAS Collaboration (G. Aad et al.), arXiv:1412.4092 [hep-ex].

254. ALICE Collaboration (J. Adam et al.), Phys. Lett. B 746, 385 (2015) (arXiv:1503.03050 [nucl-ex]).

255. ALICE Collaboration (j. Adam et al.), arXiv:1503.00681 [nucl-ex].

256. H. Paukkunen, K.J. Eskola, C. Salgado, Nucl. Phys. A 931, 331 (2014) (arXiv:1408.4563 [hep-ph]).

257. A. Bzdak, V. Skokov, S. Bathe, arXiv:1408.3156 [hep-ph].

258. M. Alvioli, B. Cole, L. Frankfurt, M. Strikman, arXiv:1409.7381 [hep-ph].

259. D.V. Perepelitsa, P.A. Steinberg, arXiv:1412.0976 [nucl-ex].

260. N. Armesto, D.C. Gulhan, J.G. Milhano, arXiv:1502.02986 [hep-ph].

261. Y.L. Dokshitzer, D.E. Kharzeev, Phys. Lett. B 519, 199 (2001) (hep-ph/0106202).

262. N. Armesto, A. Dainese, C.A. Salgado, U.A. Wiedemann, Phys. Rev. D 71, 054027 (2005) (hep-ph/0501225).

263. PHENIX Collaboration (A. Adare et al.), Phys. Rev. Lett. 98, 172301 (2007) (nucl-ex/0611018).

264. STAR Collaboration (B.I. Abelev et al.), Phys. Rev. Lett. 98, 192301 (2007) 106, 159902(E) (2011) (nucl-ex/0607012).

265. STAR Collaboration (l. Adamczyk et al.), Phys. Rev. Lett. 113, 142301 (2014) (arXiv:1404.6185 [nucl-ex]).

266. A. Adare et al., arXiv:1509.04662 [nucl-ex].

267. ATLAS Collaboration (D.V. Perepelitsa), Nucl. Phys. A 904-905, 669c (2013).

268. ALICE Collaboration (B. Abelev et al.), Phys. Rev. Lett. 109, 112301 (2012) (arXiv:1205.6443 [hep-ex]).

269. ALICE Collaboration (j. Adam et al.), arXiv:1507.03134 [nucl-ex].

270. ALICE Collaboration (j. Adam et al.), arXiv:1509.07491 [nucl-ex].

271. ALICE Collaboration (j. Adam et al.), arXiv:1506.06604 [nucl-ex].

272. ALICE Collaboration (B. Abelev et al.), JHEP 09, 112 (2012) (arXiv:1203.2160 [nucl-ex]).

273. ALICE Collaboration (b. Abelev et al.), Phys. Rev. Lett. 111, 102301 (2013) (arXiv:1305.2707 [nucl-ex]).

274. ALICE Collaboration (B.B. Abelev et al.), Phys. Rev. C 90, 034904 (2014) (arXiv:1405.2001 [nucl-ex]).

275. ALICE Collaboration (j. Adam et al.), arXiv:1509.06888 [nucl-ex].

276. ALICE Collaboration (J. Adam et al.), arXiv:1509.07287 [nucl-ex].

277. ALICE Collaboration (j. Adam et al.), arXiv:1504.07151 [nucl-ex].

278. CMS Collaboration (S. Chatrchyan et al.), JHEP 05, 063 (2012) (arXiv:1201.5069 [nucl-ex]).

279. CMS Collaboration (v. Khachatryan et al.), arXiv:1508.06678 [nucl-ex].

280. CMS Collaboration (s. Chatrchyan et al.), Phys. Rev. Lett. 113, 132301 (2014) 115, 029903(E) (2015) (arXiv:1312.4198 [nucl-ex]).

281. CMS Collaboration (V. Khachatryan et al.), arXiv:1510.03373 [nucl-ex].

282. M. Djordjevic, M. Djordjevic, B. Blagojevic, Phys. Lett. B 737, 298 (2014) (arXiv:1405.4250 [nucl-th]).

283. M. Djordjevic, Phys. Rev. Lett. 112, 042302 (2014) (arXiv:1307.4702 [nucl-th]).

284. CMS Collaboration [CMS Collaboration], CMS-PAS-HIN-12-014.

285. W.A. Horowitz, M. Gyulassy, Nucl. Phys. A 872, 265 (2011) (arXiv:1104.4958 [hep-ph]).

286. M. Nahrgang, J. Aichelin, P.B. Gossiaux, K. Werner, Phys. Rev. C 89, 014905 (2014) (arXiv:1305.6544 [hep-ph]).

287. M. He, R.J. Fries, R. Rapp, Phys. Lett. B 735, 445 (2014) (arXiv:1401.3817 [nucl-th]).

288. W.M. Alberico, A. Beraudo, A. De Pace, A. Molinari, M. Monteno, M. Nardi, F. Prino, Eur. Phys. J. C 71, 1666 (2011) (arXiv:1101.6008 [hep-ph]).

289. J. Uphoff, O. Fochler, Z. Xu, C. Greiner, Phys. Lett. B 717, 430 (2012) (arXiv:1205.4945 [hep-ph]).

290. T. Lang, H. van Hees, J. Steinheimer, M. Bleicher, arXiv:1211.6912 [hep-ph].

291. S. Cao, G.Y. Qin, S.A. Bass, Phys. Rev. C 88, 044907 (2013) (arXiv:1308.0617 [nucl-th]).

292. ALICE Collaboration (B.B. Abelev et al.), Phys. Rev. Lett. 113, 232301 (2014) (arXiv:1405.3452 [nucl-ex]).

293. J. Huang, Z.B. Kang, I. Vitev, Phys. Lett. B 726, 251 (2013) (arXiv:1306.0909 [hep-ph]).

294. T. Matsui, H. Satz, Phys. Lett. B 178, 416 (1986).

295. S. Digal, P. Petreczky, H. Satz, Phys. Rev. D 64, 094015 (2001) (hep-ph/0106017).

296. NA50 Collaboration (B. Alessandro et al.), Eur. Phys. J. C 39, 335 (2005) (hep-ex/0412036).

297. NA60 Collaboration (r. Arnaldi et al.), Phys. Rev. Lett. 99, 132302 (2007).

298. NA60 Collaboration (e. Scomparin), Nucl. Phys. A 830, 239C (2009) (arXiv:0907.3682 [nucl-ex]).

299. PHENIX Collaboration (A. Adare et al.), Phys. Rev. C 84, 054912 (2011) (arXiv:1103.6269 [nucl-ex]).

300. STAR Collaboration (B.I. Abelev et al.), Phys. Rev. C 80, 041902 (2009) (arXiv:0904.0439 [nucl-ex]).

301. R. Vogt, Nucl. Phys. A 700, 539 (2002) (hep-ph/0107045).

302. B.Z. Kopeliovich, B.G. Zakharov, Phys. Rev. D 44, 3466 (1991).

303. D.C. McGlinchey, A.D. Frawley, R. Vogt, Phys. Rev. C 87, 054910 (2013) (arXiv:1208.2667 [nucl-th]).

304. P. Braun-Munzinger, J. Stachel, Phys. Lett. B 490, 196 (2000) (nucl-th/0007059).

305. R.L. Thews, M. Schroedter, J. Rafelski, Phys. Rev. C 63, 054905 (2001) (hep-ph/0007323).

306. ALICE Collaboration (B. Abelev et al.), Phys. Rev. Lett. 109, 072301 (2012) (arXiv:1202.1383 [hep-ex]).

307. ALICE Collaboration (e. Abbas et al.), Phys. Rev. Lett. 111, 162301 (2013) (arXiv:1303.5880 [nucl-ex]).

308. ALICE Collaboration (B.B. Abelev et al.), Phys. Lett. B 734, 314 (2014) (arXiv:1311.0214 [nucl-ex]).

309. ALICE Collaboration (B.B. Abelev et al.), JHEP 02, 073 (2014) (arXiv:1308.6726 [nucl-ex]).

310. ALICE Collaboration (B.B. Abelev et al.), JHEP 12, 073 (2014) (arXiv:1405.3796 [nucl-ex]).

311. ALICE Collaboration (j. Adam et al.), JHEP 06, 055 (2015) (arXiv:1503.07179 [nucl-ex]).

312. ALICE Collaboration (J. Adam et al.), arXiv:1506.08804 [nucl-ex].

313. ALICE Collaboration (j. Adam et al.), arXiv:1506.08808 [nucl-ex].

314. ALICE Collaboration (j. Adam et al.), arXiv:1509.08802 [nucl-ex].

315. ALICE Collaboration (B.B. Abelev et al.), Phys. Lett. B 738, 361 (2014) (arXiv:1405.4493 [nucl-ex]).

316. ALICE Collaboration (B.B. Abelev et al.), Phys. Lett. B 740, 105 (2015) (arXiv:1410.2234 [nucl-ex]).

317. CMS Collaboration (S. Chatrchyan et al.), Phys. Rev. Lett. 107, 052302 (2011) (arXiv:1105.4894 [nucl-ex]).

318. CMS Collaboration (s. Chatrchyan et al.), Phys. Rev. Lett. 109, 222301 (2012) (arXiv:1208.2826 [nucl-ex]).

319. CMS Collaboration (v. Khachatryan et al.), Phys. Rev. Lett. 113, 262301 (2014) (arXiv:1410.1804 [nucl-ex]).

320. CMS Collaboration (s. Chatrchyan et al.), JHEP 04, 103 (2014) (arXiv:1312.6300 [nucl-ex]).

321. ATLAS Collaboration (G. Aad et al.), Phys. Lett. B 697, 294 (2011) (arXiv:1012.5419 [hep-ex]).

322. ATLAS Collaboration (g. Aad et al.), arXiv:1505.08141 [hep-ex].

323. LHCb Collaboration (R. Aaij et al.), JHEP 02, 072 (2014) (arXiv:1308.6729 [nucl-ex]).

324. LHCb Collaboration (r. Aaij et al.), JHEP 07, 094 (2014) (arXiv:1405.5152 [nucl-ex]).

325. X. Zhao, R. Rapp, Nucl. Phys. A 859, 114 (2011) (arXiv:1102.2194 [hep-ph]).

326. K. Zhou, N. Xu, Z. Xu, P. Zhuang, Phys. Rev. C 89, 054911 (2014) (arXiv:1401.5845 [nucl-th]).

327. A. Andronic, P. Braun-Munzinger, K. Redlich, J. Stachel, J. Phys. G 38, 124081 (2011) (arXiv:1106.6321 [nucl-th]).

328. E.G. Ferreiro, Phys. Lett. B 731, 57 (2014) (arXiv:1210.3209 [hep-ph]).

329. CMS Collaboration, CMS PAS HIN-15-001.

330. CDF Collaboration (T. Affolder et al.), Phys. Rev. Lett. 84, 2094 (2000) (hep-ex/9910025).

331. LHCb Collaboration (R. Aaij et al.), Eur. Phys. J. C 74, 3092 (2014) (arXiv:1407.7734 [hep-ex]).

332. L. Grandchamp, S. Lumpkins, D. Sun, H. van Hees, R. Rapp, Phys. Rev. C 73, 064906 (2006) (hep-ph/0507314).

333. A. Emerick, X. Zhao, R. Rapp, Eur. Phys. J. A 48, 72 (2012) (arXiv:1111.6537 [hep-ph]).

334. M. Strickland, AIP Conf. Proc. 1520, 179 (2013) (arXiv:1207.5327 [hep-ph]).

335. CMS Collaboration (S. Chatrchyan et al.), Phys. Lett. B 710, 256 (2012)(arXiv:1201.3093 [nucl-ex]).

336. ATLAS Collaboration (G. Aad et al.), arXiv:1506.08552 [hep-ex].

337. ATLAS Collaboration (g. Aad et al.), Phys. Rev. Lett. 110, 022301 (2013) (arXiv:1210.6486 [hep-ex]).

338. ATLAS Collaboration (G. Aad et al.), Eur. Phys. J. C 75, 23 (2015) (arXiv:1408.4674 [hep-ex]).

339. CMS Collaboration (S. Chatrchyan et al.), Phys. Rev. Lett. 106, 212301 (2011) (arXiv:1102.5435 [nucl-ex]).

340. CMS Collaboration (s. Chatrchyan et al.), Phys. Lett. B 715, 66 (2012) (arXiv:1205.6334 [nucl-ex]).

341. CMS Collaboration (s. Chatrchyan et al.), JHEP 03, 022 (2015) (arXiv:1410.4825 [nucl-ex]).

342. LHCb Collaboration (R. Aaij et al.), JHEP 09, 030 (2014) (arXiv:1406.2885 [hep-ex]).

343. CMS Collaboration (V. Khachatryan et al.), arXiv:1503.05825 [nucl-ex].

344. ATLAS Collaboration (G. Aad et al.), arXiv:1507.06232 [hep-ex].

345. H. Paukkunen, C.A. Salgado, JHEP 03, 071 (2011) (arXiv:1010.5392 [hep-ph]).

346. F. Arleo, K.J. Eskola, H. Paukkunen, C.A. Salgado, JHEP 04, 055 (2011) (arXiv:1103.1471 [hep-ph]).

347. V. Koch, Int. J. Mod. Phys. E 6, 203 (1997) (nucl-th/9706075).

348. NA60 Collaboration (R. Arnaldi et al.), Phys. Rev. Lett. 96, 162302 (2006) (nucl-ex/0605007).

349. NA60 Collaboration (r. Arnaldi et al.), Phys. Rev. Lett. 100, 022302 (2008) (arXiv:0711.1816 [nucl-ex]).

350. STAR Collaboration (L. Adamczyk et al.), Phys. Rev. Lett. 113, 022301 (2014) 113, 049903(e) (2014) (arXiv:1312.7397 [hep-ex]).

351. R. Rapp, Phys. Rev. C 63, 054907 (2001) (hep-ph/0010101).

352. O. Linnyk, W. Cassing, J. Manninen, E.L. Bratkovskaya, C.M. Ko, Phys. Rev. C 85, 024910 (2012) (arXiv:1111.2975 [nucl-th]).

353. A. Adare et al., arXiv:1509.04667 [nucl-ex].

354. PHENIX Collaboration (A. Adare et al.), Phys. Rev. Lett. 104, 132301 (2010) (arXiv:0804.4168 [nucl-ex]).

355. S. Turbide, R. Rapp, C. Gale, Phys. Rev. C 69, 014903 (2004) (hep-ph/0308085).

356. ALICE Collaboration (M.K. Koehler), Nucl. Phys. A 931, 665 (2014) (arXiv:1407.7809 [nucl-ex]).

357. ALICE Collaboration (p. Reichelt), J. Phys. Conf. Ser. 612, 012028 (2015) (arXiv:1412.0619 [hep-ex]).

358. ALICE Collaboration (j. Adam et al.), arXiv:1509.07324 [nucl-ex].

359. M. Klasen, C. Klein-Bsing, F. Knig, J.P. Wessels, JHEP 10, 119 (2013) (arXiv:1307.7034 [hep-ph]).

360. L.E. Gordon, W. Vogelsang, Phys. Rev. D 48, 3136 (1993).

361. W. Vogelsang, M.R. Whalley, J. Phys. G 23, A1 (1997).

362. J.F. Paquet, C. Shen, G.S. Denicol, M. Luzum, B. Schenke, S. Jeon, C. Gale, arXiv:1509.06738 [hep-ph].

363. ALICE Collaboration (A. Morreale), arXiv:1409.4456 [nucl-ex].

364. A.J. Baltz, G. Baur, D. d'Enterria, L. Frankfurt, F. Gelis, V. Guzey, K. Hencken, Y. Kharlov et al., Phys. Rep. 458, 1 (2008) (arXiv:0706.3356 [nucl-ex]).

365. ALICE Collaboration (B. Abelev et al.), Phys. Lett. B 718, 1273 (2013) (arXiv:1209.3715 [nucl-ex]).

366. ALICE Collaboration (e. Abbas et al.), Eur. Phys. J. C 73, 2617 (2013) (arXiv:1305.1467 [nucl-ex]).

367. ALICE Collaboration (j. Adam et al.), arXiv:1508.05076 [nucl-ex].

368. ALICE Collaboration (j. Adam et al.), arXiv:1503.09177 [nucl-ex].

369. V. Guzey, E. Kryshen, M. Strikman, M. Zhalov, Phys. Lett. B 726, 290 (2013) (arXiv:1305.1724 [hep-ph]).

370. M.B. Gay Ducati, M.T. Griep, M.V.T. Machado, AIP Conf. Proc. 1654, 100002 (2015).

371. T. Lappi, H. Mantysaari, Phys. Rev. C 87, 032201 (2013) (arXiv:1301.4095 [hep-ph]).

372. ALICE Collaboration (B.B. Abelev et al.), Phys. Rev. Lett. 113, 232504 (2014) (arXiv:1406.7819 [nucl-ex]).

373. LHCb Collaboration (R. Aaij et al.), J. Phys. G 40, 045001 (2013) (arXiv:1301.7084 [hep-ex]).

374. LHCb Collaboration (r. Aaij et al.), J. Phys. G 41, 055002 (2014) (arXiv:1401.3288 [hep-ex]).

375. J.L. Albacete, C. Marquet, Prog. Part. Nucl. Phys. 76, 1 (2014) (arXiv:1401.4866 [hep-ph]).

376. A. Andronic et al., arXiv:1506.03981 [nucl-ex].

377. A. Accardi et al., arXiv:1212.1701 [nucl-ex].

378. LHeC Study Group Collaboration (J.L. Abelleira Fernandez et al.), J. Phys. G 39, 075001 (2012) (arXiv:1206.2913 [physics.acc-ph]).

379. ALICE Collaboration (B. Abelev et al.), J. Phys. G 41, 087001 (2014).

380. ATLAS Collaboration, CERN-LHCC-2011-012; CERN-LHCC-2012-022.

381. CMS Collaboration (K. Hoepfner), EPJ Web of Conferences 95, 03016 (2015).

382. The LHCb Collaboration [LHCb Collaboration], CERN-LHCC-2011-001.