Published for SISSA by Springer
Received: October 14, 2016
Revised: January 18, 2017
Accepted: March 12, 2017 Published: March 27, 2017
Jorge Casalderrey-Solana,a,b Doga Can Gulhan,c Jos e Guilherme Milhano,d,e,f Daniel Pablosb and Krishna Rajagopalg
aRudolf Peierls Centre for Theoretical Physics, University of Oxford, 1 Keble Road, Oxford OX1 3NP, United Kingdom
bDepartament de F sica Qu antica i Astrof sica & Institut de Ci encies del Cosmos (ICC), Universitat de Barcelona, Mart i Franqu es 1, 08028 Barcelona, Spain
cCERN, EP Department,
CH-1211 Geneva 23, Switzerland
dCENTRA, Instituto Superior T ecnico, Universidade de Lisboa,
Av. Rovisco Pais, P-1049-001 Lisboa, Portugal
eLaborat orio de Instrumentaco e F sica Experimental de Part culas (LIP), Av. Elias Garcia 14-1, P-1000-149 Lisboa, Portugal
f Theoretical Physics Department, CERN,
Geneva, Switzerland
gCenter for Theoretical Physics, Massachusetts Institute of Technology,
Cambridge, MA 02139 U.S.A.
E-mail: mailto:[email protected]
Web End [email protected] , mailto:[email protected]
Web End [email protected] , mailto:[email protected]
Web End [email protected] , mailto:[email protected]
Web End [email protected] , [email protected]
Abstract: Within the context of a hybrid strong/weak coupling model of jet quenching, we study the modi cation of the angular distribution of the energy within jets in heavy ion collisions, as partons within jet showers lose energy and get kicked as they traverse the strongly coupled plasma produced in the collision. To describe the dynamics transverse to the jet axis, we add the e ects of transverse momentum broadening into our hybrid construction, introducing a parameter K ^
q=T 3 that governs its magnitude. We show that, because of the quenching of the energy of partons within a jet, even when K [negationslash]= 0
the jets that survive with some speci ed energy in the nal state are narrower than jets with that energy in proton-proton collisions. For this reason, many standard observables are rather insensitive to K. We propose a new di erential jet shape ratio observable in which the e ects of transverse momentum broadening are apparent. We also analyze the response of the medium to the passage of the jet through it, noting that the momentum
Open Access, c
[circlecopyrt] The Authors.
Article funded by SCOAP3. doi:http://dx.doi.org/10.1007/JHEP03(2017)135
Web End =10.1007/JHEP03(2017)135
Angular structure of jet quenching within a hybrid strong/weak coupling model
JHEP03(2017)135
lost by the jet appears as the momentum of a wake in the medium. After freezeout this wake becomes soft particles with a broad angular distribution but with net momentum in the jet direction, meaning that the wake contributes to what is reconstructed as a jet. This e ect must therefore be included in any description of the angular structure of the soft component of a jet. We show that the particles coming from the response of the medium to the momentum and energy deposited in it leads to a correlation between the momentum of soft particles well separated from the jet in angle with the direction of the jet momentum, and nd qualitative but not quantitative agreement with experimental data on observables designed to extract such a correlation. More generally, by confronting the results that we obtain upon introducing transverse momentum broadening and the response of the medium to the jet with available jet data, we highlight the importance of these processes for understanding the internal, soft, angular structure of high energy jets.
Keywords: Heavy Ion Phenomenology, Jets
ArXiv ePrint: 1609.05842
JHEP03(2017)135
Contents
1 Introduction 1
2 Brief summary of the hybrid model 5
3 Transverse kicks and jet broadening 73.1 Introduction of broadening into the hybrid model 93.2 Insensitivity of jet observables to broadening 103.3 An observable that is sensitive to broadening 17
4 Medium backreaction: observable consequences of a wake in the plasma 214.1 The spectrum of hadrons from a medium perturbed by the passage of a jet 244.2 Implementation of backreaction, background model, and jet hadronization 264.3 The e ect of backreaction on jet observables 294.4 Recovering the \lost" jet energy and \missing" pT deposited by a dijet pair 37
5 Discussion and outlook 46
A Transverse kicks kinematics 51
B Background subtraction procedure 52
1 Introduction
High energy heavy ion collisions provide a unique opportunity to explore the properties of hot, decon ned, strongly interacting matter, called quark-gluon plasma (QGP). The study of these collisions at the Relativistic Heavy Ion Collider (RHIC) and at the Large Hadron Collider (LHC) has demonstrated that matter at temperatures above the crossover between hot hadronic matter and hotter QGP exhibits strong collective phenomena [1{7] which can be described successfully by hydrodynamic simulations of the rapid expansion and cooling of the initially lumpy droplets of matter produced in the collisions [8{18]. Such strong collectivity has also recently been observed in smaller colliding systems, including p-Pb, p-p or 3He-Au [19{29], for which hydrodynamic simulations also seem to be successful [30{ 37]. The applicability of hydrodynamics from early times in the evolution and for small systems suggests that the matter formed in these ultrarelativistic collisions is a strongly coupled liquid. Support for this picture comes from analyses of collisions in strongly coupled gauge theories with a dual holographic description which show that collisions of objects with transverse size R produce a droplet of strongly coupled liquid that can be described hydrodynamically as long as the collisions are energetic enough that the temperature of the liquid at the time that it hydrodynamizes, Thydro, satis es RThydro [greaterorsimilar] 1 [38, 39] and that in
{ 1 {
JHEP03(2017)135
these collisions and in the collisions of objects that are in nite in transverse extent [40{42] hydrodynamization occurs at a time of order 1=Thydro after the collision.
The discovery that the QGP that lled the microseconds-old universe and that is recreated in collisions at RHIC and the LHC is a strongly coupled liquid challenges us to understand how such a liquid emerges from an asymptotically free gauge theory. When probed at very short length scales, the strongly coupled QGP of QCD (unlike the strongly coupled plasmas in theories with holographic descriptions) must consist of weakly coupled quarks and gluons. This makes constraining the microscopic nature of QGP via studying its interaction with energetic probes an important and interesting long-term goal. Some of the most salient examples of such probes are QCD jets. As a partonic jet shower propagates through the strongly coupled plasma created in a heavy ion collision, it su ers a strong process of energy loss as a result of its interactions with the plasma. The partons in the jet also pick up momentum transverse to their direction of motion as they are jostled during their passage through the medium. These interactions lead to a reduction in the jet energy (or quenching) and to modi cations of the properties of jets produced in heavy ion collisions relative to those of their counterparts produced in proton-proton collisions, that propagate in vacuum. These interactions also result in the transfer of energy and momentum to the plasma: the jets create a wake as they lose energy. By pursuing a large suite of jet measurements, the di erent LHC collaborations have observed strong modi cation of di erent jet observables in heavy ion collisions [43{61], making jets promising probes for medium diagnostics. The rst experimental constraints on jet quenching came from hadronic measurements at RHIC [62{64]. Analyses of jets themselves and their modi cation are also being performed at RHIC [65{68] and are one of the principal scienti c goals of the planned sPHENIX detector [69].
To fully exploit their potential as tomographic probes, a detailed understanding of the interactions of jet showers with hot QCD matter is needed. A complete theoretical description of these processes is a challenging task, due to the multi-scale nature of jet probes. On the one hand, the production of jets and the processes via which an initial hard parton fragments into a shower are governed by short distance physics that is weakly coupled. On the other hand, the interaction of jets with the medium, as well as the dynamics of softer components within jets, are sensitive to the strongly coupled dynamics of the plasma at scales of order its temperature.
One class of approaches toward making progress on this di cult theoretical problem that has been pursued intensively starts from the premise that the dynamics of the plasma itself are weakly coupled, as if the temperature of the plasma were asymptotically large, and from the premise that the interactions of energetic partons and jets with the plasma are also governed entirely by weakly coupled physics. (See refs. [70{75] for reviews.) Based on these approaches, Monte Carlo tools for analyzing jet observables are being developed [76{ 83] and many phenomenological studies of jets in medium have been confronted with LHC measurements of a variety of jet observables [81{109], including intrajet observables like those that we shall focus on [85, 89, 99, 101, 104, 106{109].
However, the observation that QGP is a strongly coupled liquid tells us that physics at scales of order its temperature is governed by strong coupling dynamics. This realiza-
{ 2 {
JHEP03(2017)135
tion has lead to many fruitful connections between the physics of the QCD plasma and the gauge/gravity duality [110]. This technique allows us rigorous and quantitative access to nonperturbative, strongly coupled, physics | including thermodynamics, transport coe cients, hydrodynamics, thermalization, response to hard probes and other real time dynamics far from equilibrium | in a large family of non-abelian gauge theory plasmas that have a dual holographic description in terms of a black hole spacetime in a gravitational theory with one higher dimension. Although the current formulation of the duality has not been shown to apply to QCD, the study of the plasmas in gauge theories that do have a holographic description has led to many insights into the dynamics of hot decon ned matter in QCD. (See refs. [111{113] for reviews.) Within this context, there have been many interesting studies that address varied aspects of the interaction between high energy probes and strongly coupled plasma [114{147]. However, in all the examples that are currently directly accessible via gauge/gravity duality, the gauge theory remains strongly coupled in the ultraviolet, which limits the direct applicability of these results for phenomenological applications.
To address the multifaceted dynamics of QCD jets in strongly coupled plasma more fully, in refs. [148, 149] we introduced a phenomenological hybrid strong/weak coupling approach to analyzing jet quenching. In this approach, we treat di erent physics processes of relevance for the interaction of developing jet showers with the quark gluon plasma di erently. In our model, the production and evolution of the jet shower is treated perturbatively, because the physics governing these processes is expected to be weakly coupled. And, we model the interaction between each of the partons formed in the shower with the medium using the result for the rate of energy loss of an energetic quark in strongly coupled plasma obtained via holographic calculations in refs. [143, 145]. The Monte Carlo implementation of this hybrid model has been successful in describing LHC measurements of a large suite of single jet, dijet and photon+jet observables [148, 149] and has been used to make predictions for more such observables and for Z+jet observables [149]. To date, the model has contained only a single free parameter, with all its successes and predictions having been obtained after tting this parameter to a single measured data point.
In this paper, after reviewing the construction of our hybrid model brie y in section 2 we will extend the model in order to be able to use it to address the angular distribution of the energy within a jet and its modi cation via its passage through the plasma, as well as the angular distribution of the energy lost by jets during their passage through the plasma. To do so, we will supplement our model with two important physics processes which were absent in our previous implementations. First, in section 3 we add \transverse momentum broadening", namely the de ection to the direction of propagation of partons as a result of the exchange of momentum transverse to their direction of motion between the parton and the medium. We assume a Gaussian distribution of the magnitude of the momentum transfer and introduce one further model parameter to specify the width of the Gaussian. Second, in section 4 we add a simpli ed model for the collective response of the medium to the passage of the jet, namely the wake in the plasma that carries the energy and momentum lost by the jet and deposited in the plasma. We shall frequently refer to this as the \backreaction of the medium". We shall not actually focus on the form of the wake itself,
{ 3 {
JHEP03(2017)135
focusing instead on the soft hadrons in the nal state that result from the hadronization of the plasma including the wake and that carry the net momentum and energy lost by the jet. We make simplifying assumptions that allow us to add a representation of the e ects of the wake on the nal state hadrons that respects energy and momentum conservation, without introducing any new parameters into our hybrid model. Our principal assumption is that the energy lost by the jet thermalizes quickly, in the sense that it quickly becomes a hydrodynamic wake in the plasma that carries the \lost" energy and momentum, which in turn after freezeout becomes soft particles spread over a calculable and wide range of angles relative to the jet. This is an immediate and natural consequence of strong coupling dynamics [123]. Something similar can happen at weak coupling even though the energy is initially lost by gluon radiation because these radiated gluons can experience a cascade of reinteractions that converts the energy into soft particles at large angles [150{156]. The e ects of both transverse momentum broadening [85, 86, 107{109] and the backreaction of the medium [95, 157{162] on jet observables have also been studied within the context of perturbative energy loss mechanisms.
After adding broadening and backreaction to our hybrid model, over the course of sections 3 and 4 we will assess their e ects on measurements of di erent jet and dijet observables that are sensitive to the angular structure of jets, including the dependence of jet suppression on the reconstruction parameter R which controls the angular size of the jets that are found and reconstructed in a sample of events, several observables based upon the jet shape including a new di erential jet shape ratio that we introduce, the dijet acoplanarity, and the balance of momenta along the dijet axis carried by particles in the event with a given momentum. We also return to the jet fragmentation function, as we wish to see how it is modi ed by the soft particles coming from the backreaction of the medium.
One of our signi cant ndings is that jets with a given energy that were produced in a heavy ion collision and so have had to traverse a droplet of QGP have a narrower hard jet core than jets with the same energy that were produced in vacuum, even if a substantial degree of transverse momentum broadening is turned on. This is a consequence of the fact that wider jets typically contain more, and less hard, partonic fragments and lose more energy than narrower jets with the same energy. The same phenomenon has been observed in calculations of jet quenching that are entirely done at strong coupling [145, 147] and in Monte Carlo calculations of radiative energy loss that are entirely done at weak coupling [105]. We nd that this observation leads to two unexpected consequences. First, jets with a larger reconstruction parameter R are more suppressed (the suppression factor RAA of jets is pushed downward) at least for moderate R. Second, the intrajet angular distribution of energy in the quenched jets is remarkably independent of the amount of transverse momentum broadening. Note that we are only able to see that energy loss causes jets with a given energy to be narrower in heavy ion collisions than in proton-proton collisions because we have incorporated event-by-event (maybe better phrased jetby-jet) variations in the fragmentation pattern of jets; this physical e ect, and its striking consequences, are absent in approaches based on analyzing how an average jet is modi ed by passage through the plasma, as for example in refs. [107, 109].
{ 4 {
JHEP03(2017)135
Although we nd that both the standard jet shape and the dijet acoplanarity are remarkably independent of the amount of transverse momentum broadening, we are able to construct a new observable | essentially a jet shape ratio that is partially di erential in pT | that is sensitive to the amount of transverse momentum broadening that we introduce. This points out a path toward detecting experimental evidence for the e ects of transverse momentum broadening within a jet, which would be a very important rst step toward using jets to resolve the microscopic structure of quark-gluon plasma.
Our last signi cant nding is that when we implement the collective response of the medium to the passage of the jet, the energy lost from the jet ends up in the form of soft particles separated from the jet axis by very large angles and is in qualitative agreement with experimental measurements of observables that are referred to as \missing-pT" distributions which have recently been reported by CMS [60]. Careful comparison between our calculations of these and several other observables and experimental data in section 4 indicates that our simpli ed treatment of the wake produces slightly too many very soft particles (pT < 2 GeV) at large angles and not enough particles with momenta in the 2 4 GeV range. This is not entirely unexpected because the approximations via which we
treat the particles coming from the wake are reliable only for particles with momenta that are not a lot larger than the freezeout temperature, but it may also be an indication that the wake in the plasma does not actually thermalize as fully as we assume in our simpli ed treatment.
We close in section 5 by discussing the various results of our analyses, identifying further improvements of our implementation of the in-medium dynamics of jets in strongly coupled plasma for the future, and looking ahead at the path toward using jets to resolve the microscopic properties of quark-gluon plasma.
2 Brief summary of the hybrid model
In this section, we provide a brief description of the hybrid model which we will employ to describe the modi cation of jets produced in heavy ion collisions that propagate through a droplet of hot matter relative to those produced in proton-proton collisions that propagate in vacuum. A more detailed account of the model may be found in refs. [148, 149].
The main motivation for introducing this model is to separate the strongly coupled dynamics of quark-gluon plasma itself and of interactions between it and partons plowing through it from the weakly coupled dynamics governing the production, showering and relaxation of virtuality of high energy QCD jets. Since for any parton that showers and forms a jet the initial virtuality of the parton is much larger than any scale associated with the medium, of order its temperature T , the rst assumption of the model is that the evolution of the jet proceeds as in vacuum, with the branching of the parton shower unmodi ed by the presence of the strongly interacting plasma. Upon making this assumption, the modi cation of jet showers is only due to the interaction of each of the partons in the jet with the strongly coupled medium. After associating each parton in the jet with a life-time determined via a formation time argument [148, 149, 163], we compute the energy lost by each parton as it propagates in the strongly coupled plasma. Since the interactions of each
{ 5 {
JHEP03(2017)135
of these partons with the medium is sensitive to the medium scale, the rate of energy loss is controlled by strongly coupled dynamics. The second assumption of the model is that the rate at which a parton loses energy can be modeled by the rate of energy loss of light quark jets in the strongly coupled plasma of N = 4 supersymmetric Yang-Mills (SYM)
theory which has been computed via holography and is given by [143, 145]
dE dx
[vextendsingle][vextendsingle][vextendsingle][vextendsingle]
qx2stop x2; (2.1)
with xstop the distance over which the light quark jet would lose all of its energy if it propagated through plasma at a constant temperature T . In N = 4 SYM theory, jets with
a given initial energy Ein can have a wide range of initial opening angles, with the narrower jets having larger values of xstop [145]. There is a minimum possible initial opening angle,
corresponding to the maximum possible xstop for jets with a given initial energy [145] that
was computed holographically in refs. [131, 137, 141] and is given by
xstop = 1 2 sc
1
=
4 Ein
x2 x2stop
strongly coupled
JHEP03(2017)135
E1/3in
T 4/3 ; (2.2)
with sc = 1:05 1/6. In the hybrid model, we apply (2.1) parton-by-parton to each parton in a QCD DGLAP shower as described by Pythia [164], rather than attempting to use the
N = 4 SYM jets that lose energy at the rate (2.1) themselves as models for QCD jets in
plasma, as in ref. [147]. Because we apply (2.1) to individual partons, we use the form (2.2) for xstop appropriate for the skinniest possible N = 4 SYM jets. We shall further assume
that the most salient di erences between the strongly coupled limit of N = 4 SYM theory
and QCD can be incorporated via varying the value of sc, which becomes the only tting parameter of the hybrid model formulated in refs. [148, 149].
We have implemented this hybrid model into a Monte Carlo simulation in which hard jets, showering as described by Pythia [164], are embedded within a droplet of hot matter produced in a heavy ion collision, expanding and cooling as described by relativistic viscous hydrodynamics. To generate the hard QCD jet shower, we employ Pythia8.183 [164], allowing the DGLAP shower to evolve down to a minimum transverse momentum of pTmin = 1 GeV. We distribute these hard events in the transverse plane of the heavy ion collisions according the number of binary nucleon-nucleon collisions. The trajectory of each jet is tracked from the generation point as it propagates through the expanding cooling plasma until the jet reaches a region where the temperature has dropped below a temperature Tc, below which we assume no further energy loss occurs. This Tc is not sharply de ned but it should presumably be near the crossover between quark-gluon plasma and hadronic matter and we therefore vary its value over the range 145 < Tc < 170 MeV.
Seeing how our results vary as we vary Tc over this range serves as a gauge of some of the uncertainties in our model. The energy loss rate is computed via eq. (2.1), with xstop
evaluated according to eq. (2.2) at the local temperature as given by the hydrodynamic simulations of refs. [16, 17]. Flow e ects are taken into account by evaluating the rate of energy loss in the local uid rest frame, as explained in ref. [149]. In refs. [148, 149],
{ 6 {
we have tted the value of the parameter sc to a single measurement, the suppression of the number of jets with one transverse momentum in LHC heavy ion collisions with one centrality, and have then successfully confronted this hybrid model with measurements of many single jet, dijet and -jet observables as functions of jet transverse momentum and collision centrality and made predictions for many further measurements of these types.
In the next sections, we shall extend our implementation of the hybrid model to include two new physics processes, transverse momentum broadening in section 3 and the response of the medium to the jet in section 4, and will evaluate their consequences for intra-jet observables, in particular those related to the angular structure of jets.
3 Transverse kicks and jet broadening
Previously in refs. [148, 149], all the e ects of the strongly coupled medium on the properties of the jets arise as a consequence of the energy lost by the partons in the jet shower as they plow through the medium. In this section, we augment our hybrid model by adding a second physical process, and hence a second free parameter, namely the kicks transverse to their direction of motion that the partons in the jet receive as a consequence of plowing through the medium. This process has long been referred to as \transverse momentum broadening" based upon the expectation that the consequence of the kicks in random transverse directions received by the many partons in a jet shower will be broadening of the jet. In perturbative calculations, transverse momentum broadening arises from the multiple soft exchanges of momentum that a parton su ers as it traverses a medium leading to a random change in its momentum and in particular providing the parton with some additional momentum perpendicular to its original direction of propagation.
If the energetic parton su ers multiple soft exchanges as it traverses the medium, the distribution of the momentum transferred via this stochastic process is well approximated as Gaussian. As a consequence, the transverse momentum distribution of partons that have traversed a medium of length L is approximately Gaussian with a width that scales with the medium length, Q2? = ^
qL. The quantity ^
q that arises here is called the momentum broadening parameter; this property of the medium codi es the typical squared momentum that the medium transfers to the probe per unit length. It has dimension 3 and in a plasma in thermal equilibrium with temperature T it is proportional to T 3, up to a possible logarithmic dependence on the ratio of the parton energy to T . In perturbative calculations of energy loss via gluon radiation, the medium parameter ^
q also determines the intensity of the gluon radiation induced by the medium, and hence is related to energy loss as well as to momentum broadening. In our model, we introduce ^
q only as a way of parametrizing
momentum broadening.
In the strongly coupled plasma of N = 4 SYM theory, transverse momentum broaden
ing has been calculated holographically for both heavy quarks [116, 121, 122] and massless quarks [115, 132]. Although there is no notion of scattering centers and no notion of multiple discrete transfers of momentum in the strongly coupled limit, both heavy and massless quarks pick up transverse momentum as they propagate through the hot strongly coupled liquid and the resulting transverse momentum distribution is Gaussian with a width
{ 7 {
JHEP03(2017)135
q / p T
3, with the t Hooft coupling. However, unlike at weak coupling there is no strong correlation between the dynamics responsible for transverse momentum broadening and that responsible for parton energy loss.
We shall introduce broadening into our hybrid model by assuming a Gaussian distribution for the transverse momentum picked up by each parton in the shower, with the momentum squared picked up per distance travelled given by
^
q = KT 3 (3.1)
with T the local temperature of the medium at the location of a given parton at a given time and with K a theory-dependent constant that we shall treat as a free parameter that should ultimately be determined via tting to data. Since the medium is dynamical, with longitudinal and transverse expansion, T and hence ^
q varies with position and time. We shall obtain this dependence from the hydrodynamic simulation of the expanding cooling droplet of plasma in our hybrid model.
For massless or very energetic particles, in N = 4 SYM theory in the large number
of colors (Nc) limit the value of K has been calculated holographically [115, 132] and is given by K = KN=4 [similarequal] 24 for = 10, a value of the t Hooft coupling = g2Nc that
corresponds to g2=4 [similarequal] 0:27 for Nc = 3. KN=4 is proportional to p as noted above. In
a large class of conformal theories KCFT=KN=4 =
psCFT=sN=4 [120], s being the entropy density. This suggests that KQCD is likely to be less than KN=4 since, at least at weak
coupling, sQCD=sN=4 [similarequal] 0:40. An alternative approach to gaining an expectation for the
likely value of KQCD is to start from a perturbative analysis of parton energy loss, in which the value of K controls energy loss via gluon radiation and, via this relation, can be related to experimental observables that are sensitive to parton energy loss like for example the suppression in the number of high-momentum hadrons in heavy ion collisions as compared to proton-proton collisions. The JET Collaboration has pursued this approach [165]; the value of ^
q that they have found corresponds to a value of K given by Kpert [similarequal] 5. To date, nobody has extracted a value of K via comparison to data on experimental observables that are directly sensitive to transverse momentum kicks and jet broadening. (A recent pioneering attempt [109] yields values ranging from 0 to several times larger than the value of Kpert obtained in ref. [165].) Our goal in this section is to introduce transverse momentum broadening into our hybrid model, treating K as a parameter that should ultimately be determined via comparison between calculations of the (K-dependent) observable consequences of momentum broadening to experimental measurements of observables that are directly sensitive to this physical process.
In this section, we will use our hybrid model to analyze the consequences of transverse momentum broadening for various di erent jet observables. Based upon the discussion above, we expect that K actually lies somewhere around 5 to 20. But, in order to better understand the consequences of broadening, we shall investigate the e ects on observables of varying K over the wide range 0 K 100. Our principal conclusion, after investigating a
suite of jet observables, will be that most observables, even those tailored to measuring the angular structure of jets, are remarkably insensitive to broadening, showing little sensitivity to K over the full range that we explore. This conclusion becomes less surprising once we
{ 8 {
Q2? = ^
qL with ^
JHEP03(2017)135
recall that even if every individual jet is broadened by its passage through the medium, jets with a given energy can end up being narrower in heavy ion collisions than in proton-proton collisions. This happens if wider jets lose more energy than narrower jets (as is the case in perturbative calculations in QCD [166] and in holographic calculations of jets in strongly coupled N = 4 SYM theory [145]) and if the probability distribution for the jet production
absent any medium e ects is a steeply falling function of jet energy, as is the case. Recent weakly coupled Monte Carlo calculations of an ensemble of QCD jets in heavy ion collisions as compared with those in proton-proton collisions [105] and recent holographic analyses of how propagation through plasma modi es the energy and opening angle distribution of an ensemble of N = 4 SYM jets with initial distributions as in perturbative QCD [147] both
provide clear illustrations of how jets with a given energy can end up narrower in heavy ion collisions even though every individual jet broadens as it propagates through plasma. We shall nd the same in our hybrid model. The resulting insensitivity of jet observables to the value of K will make it quite challenging to extract the value of this medium parameter from data.
In section 3.3 we shall introduce a new observable which does exhibit considerable sensitivity to in-medium transverse momentum broadening, proposing this observable as a possible route to using future experimental measurements to constrain the value of K.
3.1 Introduction of broadening into the hybrid model
We shall assume that each parton in a jet shower picks up some transverse momentum as it propagates through the plasma for a time dt, with the transverse direction chosen randomly and with the magnitude of the momentum chosen from a Gaussian distribution with a width ^
q dt with ^
q speci ed in terms of the local temperature of the plasma T and the parameter K according to (3.1). Many previous computations of in-medium energy loss and broadening of jets have been performed in the limit of a static uid at rest. However, we shall study the interactions of jets with an expanding, cooling, droplet of plasma as described by hydrodynamics. As in our previous analysis of jet observables related to parton energy loss [149], we shall apply results appropriate to a static medium in the local uid rest frame, meaning in the present case that we are neglecting any e ects of gradients in the uid on transverse momentum broadening. This prescription implies that all modi cations to the momentum of the parton (loss of longitudinal momentum according to (2.1) as previously and the transverse momentum kick that we are introducing here) are computed in the rest frame of the uid at the location of the parton at a given time. Note that the transverse momentum kicks are transverse to the direction of motion of the parton in the local uid rest frame, meaning that they need not be transverse to its momentum in the collision center-of-mass frame. We relegate the details of the transformation between the uid frame and the collision frame to appendix A.
Following our Monte Carlo approach, we perform a full simulation of an ensemble of jets. For the analysis performed in this section, we studied 500,000 jet events, generated and evolved by Pythia 8.183 [164]. The point of origin in the transverse plane of the hard processes is distributed by sampling the binary collison probability density while initial transverse direction and rapidity are retained from Pythia.
{ 9 {
JHEP03(2017)135
We embed these jets in the hydrodynamic background of refs. [16, 17]. We follow the trajectories of all the partons in the Pythia jet shower as they traverse the hot matter created in the collision, as described in the hydrodynamic simulation. Between the time when each parton is created at a branching event and the time when each parton itself branches, we discretize its trajectory and at each point we add to its momentum a random transverse momentum chosen according to a Gaussian distribution of width ^
q dtF with dtF
the length of the discretized time interval in the uid frame. We have checked that our results do not depend on our choice of dtF .
As described in detail in appendix A, we implement broadening in the local uid rest frame by assuming that in this frame the energy and virtuality of a parton do not change when the parton is kicked; the only thing that changes due to the tranverse kicks it receives is the direction of the partons momentum vector. Again as described in detail in appendix A, upon boosting back from the local uid rest frame to the collision center of mass frame there will in general be a change to the collision frame energy of the parton. It turns out that more often than not in the collision frame the partons lose a small amount of energy as a consequence of the transverse kicks they receive. In addition, the transverse kicks they receive may push some partons outside the jet. (Since the jet shape falls with distance from the jet axis, although transverse kicks may also push partons into the jet this is less likely.) These e ects together mean that, summed over the whole jet evolution, the transverse momentum kicks that we are adding will result in a slight increase to the overall jet energy loss at a xed value of sc. That is, the dynamics of broadening leads to a small increase in the quenching of jets. This means that for each nonzero value of K we need to re t the value of the parameter sc that, through (2.1) and (2.2), determines the amount of energy loss in our hybrid model. We do so in gure 1, nding that the e ect is small. For values of K around 5 to 20, the e ect of broadening implies a reduction of less than 5% in the value of sc relative to that reported in our previous work. This e ect is much smaller than the uncertainties represented by the width of the band in gure 1. Even for the extreme value K = 100, sc is only reduced by about 10%.
3.2 Insensitivity of jet observables to broadening
Having xed the dependence of the quenching parameter sc on the broadening parameter K, we can now begin our exploration of the e ect of broadening on di erent observables. With an ensemble of events in hand, generated by Pythia and modi ed to include energy loss as described in section 2 and broadening as described above, the rst step in the calculation of any jet observable is the nding and reconstruction of jets in each of the events in the ensemble. We do so using the anti-kt algorithm [167] as implemented in the
FastJet package [168].
The rst observable that we consider is the suppression factor of jets (namely the ratio RjetAA of the number of jets with speci ed kinematics in heavy ion collisions to the number of jets with the same kinematics in proton-proton collisions) as a function of the jet reconstruction parameter R that arises in the anti-kt algorithm [167]. The anti-kt algorithm groups particles within R of each other in the ( ; ) plane into what it de nes as a single jet whereas if R were smaller it may reconstruct the same particles as several smaller
{ 10 {
JHEP03(2017)135
0.6
0.5
0.4
sc
0.3
0.2
0.1
0 20 40 60 80 100
K
JHEP03(2017)135
Figure 1. When we introduce transverse momentum broadening via a nonzero value of the broadening parameter K, this introduces a small increase in energy loss. This means that for each nonzero value of K we need to re t the value of the energy loss parameter sc to the measured value [49]
of RjetAA for jets with 100 GeV< pT < 110 GeV and 2 < < 2, as in refs. [148, 149]. The resulting
dependence of sc on the broadening parameter K is mild. As in refs. [148, 149], the width of the band of values of sc in this gure (and the consequent widths of the bands depicting our hybrid model predictions in subsequent gures) comes both from the error bar on the experimentally measured data point used to x sc and from varying the crossover temperature Tc as described in section 2 in order to get some sense of the systematic uncertainties in the hybrid model.
1.2
R=0.2 R=0.3 R=0.4 R=0.5
Data
0-10% Centrality K=0
0 100 120 140 160 180 200 220 240 260 280 300
1.2
R=0.2 R=0.3 R=0.4 R=0.5
Data
0-10% Centrality K=40
0 100 120 140 160 180 200 220 240 260 280 300
1
1
0.8
0.8
JetR AA
0.6
JetR AA
0.6
0.4
0.4
0.2
0.2
PT (GeV)
PT (GeV)
Figure 2. Dependence of RjetAA on the anti-kt jet reconstruction parameter R for K = 0 (no broadening, left panel) and K = 40 (right panel). For comparison, we show the suppression of jets reconstructed using the anti-kt algorithm with R = 0:3 as measured by CMS in the same interval of rapidity [49].
narrower jets. Hence, choosing a larger R translates into reconstructing an ensemble of jets that tend to be wider in angle. Because of this, the dependence of the suppression factor RjetAA on R is often considered a proxy for modi cation of the angular structure of jets as a consequence of their interaction with the medium. A naive expectation, then, would be that turning on transverse momentum broadening should make jets broader and that this in turn should leave some imprint in the R-dependence of RjetAA, increasing the suppression for smaller R. This is not at all what we observe.
In gure 2 we show the R dependence of the suppression factor RjetAA for K = 0,i.e. no broadening, and for K = 40 for jets with [notdef] [notdef] < 2 as a function of pT . The data
points come from measurements of jets reconstructed with R = 0:3; the colored bands
{ 11 {
are the results from our model for jets reconstructed with four di erent values of R. The left-most data point in each panel is the one that we have used to constrain the value of sc, as we described in section 3.1. The calculations shown in the two panels therefore have di erent values of sc, although the di erence is small. In both panels, although the suppression factor shows only a very weak dependence on the reconstruction radius R, wide (larger R) jets tend to be somewhat more suppressed than narrow (smaller R) jets. This is a trend that we foreshadowed in the introduction to section 3. When comparing jets at the same energy, wide jets contain more fragments (in our hybrid model, simply more partons) than narrow jets. Since the larger the number of partons traversing the medium the larger the lost energy, wide jets are naturally more suppressed. The same qualitative dependence of the suppression factor on the jet size has recently been observed in holographic computations [145, 147]. Measurements of the R-dependence of RAA for
these high energy jets by the CMS collaboration show very little sensitivity to the value of R [49], as in our simulations, although the uncertainties in our calculations and in the measurements preclude a quantitative comparison at present. We shall return to the R-dependence of RAA in section 4.
In spite of the extreme transverse momentum broadening introduced by choosing K = 40, the R-dependence of the observable plotted in gure 2 is almost identical in the two panels. The origin of this lack of sensitivity to K lies in the strong quenching of jets by the plasma, and in particular in the fact that, as we have already noted, wider jets are more strongly quenched than narrower jets. In particular, the softer partons within a wider jet that could serve to further broaden the jet as they are kicked in transverse momentum instead lose almost all of their energy. This means that the jet sample that ends up dominating the inclusive jet spectrum ratio RjetAA is biased such that most jets in the sample contain only a few or even just one hard parton. For such jets, transverse momentum kicks, even with an extreme value of K, serve only to change the direction of the jet axis, not to broaden the jet.
We turn now to the second of the three observables that we shall analyze in this section, one that we will use to look for exactly the change in the direction of the jet axis due to transverse momentum broadening, namely due to the vector-summed e ect of the transverse momentum kicks felt by each of the partons in the jet. If we consider a dijet pair, the change in the direction of propagation of the two jets in the pair will in general di er, since broadening is a stochastic process and also since the temperature as a function of time along the trajectory of each of the jets will in general be di erent. If all dijet pairs were produced back-to-back in azimuthal angle, with = , deviations from = due to di erent broadening-induced kicks to the two jets in the pair could be used as a direct measurement of broadening. Reality is not this simple. The hard scattering processes that produce dijets often include radiative production of more than two partons, and in many events where two jets are reconstructed there may in fact have been a third or even fourth jet present also. This means that even in proton-proton collisions there is a nontrivial distribution of , centered around but with considerable width. We shall start from this distribution, and then look at the e ects on it due to propagation through the plasma originating from energy loss and broadening.
{ 12 {
JHEP03(2017)135
5
3
K=0 K=20 K=40 K=100 Vacuum
P L
t > 120 GeV, P S
t > 30 GeV, | | < 2
0 10%
4.5
K=0 K=20 K=40 K=100 Vacuum
P L
t > 35 GeV, P S
t > 10 GeV, | | < 2
0 10%
4
2.5
3.5
2
EventFraction
2.5
EventFraction
3
1.5
2
1.5
1
1
0.5
0.5
0 0 0.5 1 1.5 2 2.5 3
0 0 0.5 1 1.5 2 2.5 3
JHEP03(2017)135
Figure 3. Dijet acoplanarity distribution for high-energy (left) and low-energy (right) dijets in LHC heavy ion collisions with ps = 2:76 ATeV for two di erent values of the broadening parameter K. For comparison, the black dots show the acoplanarity in proton-proton collisions as simulated by Pythia.
Deviations of away from are termed acoplanarity because the two jet axes and the beam axis are not coplanar. The black dots in gure 3 shows the acoplanarity distribution for dijets in proton-proton collisions (propagating in vacuum) for two di erent cuts on the transverse momentum of the dijet pair: pTleading > 120 GeV and pTsubleading > 30 GeV (left)
and pTleading > 35 GeV and pTsubleading > 10 GeV (right) as simulated by Pythia. The
-distribution is di erent in the two momentum regions displayed, with the high-energy acoplanarity slightly narrower (closer to = ) than the low energy one. This occurs because the fragmentation of higher energy jets leads to a narrower angular distribution of fragments than for lower energy jets, and in a case where only two jets are produced the acoplanarity will be less if the jets are narrower.
As already pointed out in refs. [89, 108, 109], in the case of the high energy dijet pairs the e ects of broadening on their acoplanarity is much smaller than the width of the vacuum acoplanarity distribution. Indeed, in the left panel of gure 3 we see that our results with K = 0 (no-broadening) and K = 100 (extreme broadening) are both indistinguishable from the vacuum distribution.
In the right panel of gure 3, we look at the acoplanarity distribution for dijets with much lower energy, choosing dijets with pTleading > 35 GeV and pTsubleading > 10 GeV. It
is challenging, perhaps prohibitively so, to measure jets with energies as low as this in LHC heavy ion collisions, but even as an academic study the results are interesting. First, our K = 0 calculation apparently yields a narrower distribution than in vacuum. One contribution to the origin of this apparent e ect lies in the increase in the number of dijet pairs close to = 0; since what is plotted is a normalized probability distribution, this tends to lower the curve near = . The enhancement of almost collinear jet pairs is a consequence of energy loss and arises when a jet that had been propagating somewhere near = loses so much energy that its energy drops below that of a third jet produced close to = 0, meaning that what is reconstructed is a dijet with 0. We have checked
that if we restrict the dijet distribution to pairs of jets moving in opposite hemispheres, the K = 0 and vacuum acoplanarity distributions are much more similar. The second reason
{ 13 {
8
7
6
5
4
3
2
1
K=0 K=5 K=20 K=50 Vacuum
P L
t > 35 GeV, P S
t > 10 GeV, | | < 0.5
0 10%
Figure 4. Dijet acoplanarity distribution for dijets in heavy ion collisions at RHIC for di erent values of the broadening parameter K. For comparison, the black dots show the acoplanarity distribution in proton-proton collisions as simulated by Pythia.
why the distribution of around is slightly narrower in heavy ion collisions than in vacuum is that a quenched jet seen in a heavy ion collision with a given energy began with a larger energy, and the acoplanarity distribution for higher energy dijets as produced in vacuum is narrower.
We can now look at the acoplanarity distribution with K = 100 for the low energy dijets in the right panel of gure 3. We see that introducing this extreme degree of broadening does make the acoplanarity distribution very slightly wider than the K = 0 distribution, perhaps by coincidence bringing it back into agreement with the vacuum acoplanarity distribution. The principal conclusion, though, is that in LHC heavy ion collisions, even for dijets with very low energies and even with substantially more broadening than is expected, the e ects of broadening on the acoplanarity distribution are very small.
As noted in refs. [108, 109], the e ects of broadening on acoplanarity distributions are larger for low energy jets at RHIC, where the smaller soft background makes these measurements more feasible. Furthermore, the vacuum acoplanarity distribution is much narrower in RHIC heavy ion collisions than at the LHC for dijets with a comparable range in transverse momentum. As we now explain, this can be attributed to the fact that the jet spectrum at a given jet energy is more steeply falling in lower energy RHIC collisions than it is at the LHC. One way in which dijets with a given energy pick up acoplanarity is if the subleading jet in the pair started out at a higher energy and radiated a relatively hard gluon, which became a third jet that balances the acoplanarity in the momenta of the reconstructed dijet pair. If at a given energy the spectrum is a more steeply falling function of energy this means that there are fewer initially higher energy jets that could radiate and end up at the given energy. Thus, a steeper spectrum as in RHIC collisions yields a narrower acoplanarity distribution. This is illustrated by the black dots in gure 4 that show the acoplanarity distribution in RHIC heavy ion collisions for dijet pairs with pTleading > 35 GeV and pTsubleading > 10 GeV, the same dijet energies as in the right panel
of gure 3 although here we have chosen dijets with [notdef] [notdef] < 0:5. Note that here the vacuum
distribution shows a small accumulation of events at [lessorsimilar] 0:3. This is a result of the narrow rapidity coverage [notdef] [notdef] < 0:5: in a small fraction of events, the jet that balances
{ 14 {
EventFraction
0 0 0.5 1 1.5 2 2.5 3
JHEP03(2017)135
most of the transverse momentum of the leading jet falls outside the accepted range and what is reconstructed as a dijet is the leading jet and a subleading jet pointing in a similar direction.
We see in gure 4 that the RHIC acoplanarity distribution becomes visibly narrower in heavy ion collisions if we neglect broadening, setting K = 0. This e ect is perhaps hinted at in the right panel of gure 3 but becomes clearly visible at RHIC (energy loss, and therefore this e ect, are not included in the calculations of refs. [108, 109]). It is a consequence of strong energy loss and a steeply falling jet spectrum in sum. Wider jets with a given energy loses more energy than narrower jets, meaning that the jets that start out at in a given energy bin and stay in that bin are the narrow jets. Because of the steeply falling spectrum, there are not many jets that originate with higher energies, lose energy, and end up in the given energy bin. The result is a narrowing of jets that remain in a given energy bin, something that has been seen previously in both perturbative [105] and holographic [145, 147] analyses. And, we see from the K = 0 results in gure 4 that narrower dijets are less acoplanar. When we now turn on broadening, in addition to energy loss, the jets get broader and the acoplanarity increases. The e ect of broadening is small for realistic values of K in the 5-20 range; it takes unrealistically large values of K to broaden the jets su ciently that the acoplanarity distribution becomes as wide as in vacuum.
The lesson here is that acoplanarity is to some degree sensitive to both energy loss and broadening, more so for lower energy jets in lower energy collisions where the spectrum is more steeply falling. But, even in the best case, this observable exhibits little sensitivity to broadening, with the narrowing of the acoplanarity distribution due to energy loss being greater than the broadening of the acoplanarity distribution due to momentum broadening with realistic values of K. The competition between these two opposing e ects would complicate any attempt to use acoplanarity measurements in isolation to see e ects of broadening, leave apart to extract the value of the medium broadening parameter K.
The third observable that we shall analyze is called the jet shape and is an intrajet observable that is a measure of the angular distribution of the energy within a jet. The jet shape is de ned as the fraction of the jet energy in jets reconstructed with a given anti-kt parameter R that is contained within an annulus of radius r and width r (in space)
centered on the jet axis. Following the analysis in ref. [51], we de ne the di erential jet shape as
(r)
1 Njets
JHEP03(2017)135
1 r
Xjets
Pi2 r[notdef] r/2 pi,trackt
pjett
(3.2)
for r < R, where the tracks in the sum dont necessarily have to belong to the jet constituents de ned through the anti-kt clustering. For this reason the nal jet shape distribution is multiplied by the event averaged factor
Dpjett=
PNbinsi=0 ptrackt(ri)
[angbracketrightBig]
so that it is
normalized to one. We show the result of this analysis in gure 5, where we compare the ratio of the jet shape for the quenched jets in PbPb collisions to that for the unquenched jets in proton-proton collisions. In this analysis, the sum over jets in (3.2) includes all
{ 15 {
100 < P jetT < 300 GeV 0.3 < | | < 2, r < 0.3
P partonT > 1 GeV
K=100
K=40 K=20
K=0 CMS Data
0-10% Centrality
K=100
K=40 K=20
K=0 CMS Data
10-30% Centrality
1.8
1.6
1.8
1.6
1.4
1.4
PbPb/pp
PbPb/pp
1.2
1.2
1
1
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0 0.05 0.1 0.15 0.2 0.25 0.3
0 0.05 0.1 0.15 0.2 0.25 0.3
r
r
Figure 5. Left: ratio of the jet shape in PbPb collisions to that in proton-proton collisions for di erent values of the broadening parameter K as compared to CMS data from heavy ion collisions with ps = 2:76 ATeV [51] and 0-10% centrality (left) or 10-30% centrality (right).
jets with pjett > 100 GeV and 0:3 < [notdef] [notdef] < 2. For reference, we also show the experimental results for this ratio, as measured by the CMS collaboration [51].
The hybrid model calculations with K = 0 shown in gure 5 provide clear con rmation that, as we have discussed, energy loss serves to narrow the angular size of jets with a given energy in heavy ion collisions relative to that of jets with the same energy in proton-proton collisions. Again, this arises because wider jets lose more energy than narrower jets, leaving behind a sample that is dominated by narrower jets. Adding transverse momentum broadening by turning on a nonzero value of K serves to broaden the jets in the sample, slightly. The e ect of broadening is very small even for the unrealistically large choice K = 100.1 (Note that transverse momentum picked up by the one or few hardest partons in the jet serves to de ect the angle of the reconstructed jet, a ecting the acoplanarity. This has little e ect on the jet shape since r is measured relative to the center of the reconstructed jet, not relative to whatever the original direction of its parent parton was.) It is apparent in gure 5 that our analysis in this section does not do a good job of describing the jet shape ratio measured in experimental data [51], in particular at larger values of r. As we shall discuss in section 4, at larger and larger r the partons in the reconstructed jet are softer and softer. Our hybrid model fails to describe a medium-induced enhancement in the production of soft particles at large angles relative to the jet direction seen in heavy ion collisions relative to proton-proton collisions. This enhancement, which does not contribute much to the overall jet energy, points to the existence of soft modes moving in the same direction as the jet, which is to say a moving wake in the plasma that the jet itself excites as it loses energy and momentum to the plasma. In our implementation of the hybrid
1In the results of ref. [107], introducing broadening results in larger e ects on the jet shapes. There are at least two reasons for this. First, the analysis of ref. [107] focuses on the modi cation of average jets. Leaving out jet-by-jet uctuations means that this analysis cannot see that wider jets lose more energy than narrower ones, and so does not include the consequent narrowing of jets with a given energy in PbPb collisions relative to jets with the same energy in proton-proton collisions. And, second, in ref. [107] all the partons present in a Pythia shower just before hadronization are assumed to have passed through the medium, and felt the e ects of broadening, whereas in our calculation of the development of the jet shower in spacetime we see that many partons are produced at splittings that occur after the jet has already departed from the medium, and therefore do not feel any transverse momentum kicks.
{ 16 {
JHEP03(2017)135
model to this point, no e ects of such a wake are included. We are in e ect making the assumption that all the energy lost by the jet is so fully thermalized that it ends up as a little bit more plasma or a little bit hotter plasma, becoming a distribution of soft thermal particles in the nal state that is isotropic, uncorrelated with the jet direction. This cannot be the whole story since, after all, momentum is conserved and the jet loses momentum to the plasma as well as energy. We shall return to this in section 4.
3.3 An observable that is sensitive to broadening
Before returning to the question of where the momentum and energy lost by the jet ends up and how this a ects angular observables like the jet shape, we close this section by identifying a more di erential observable that is directly sensitive to transverse momentum broadening, meaning to the value of the broadening parameter K.
Let us recapitulate why the observables that we looked at above are not sensitive to broadening. Kicks received by the few highest momentum partons in a jet contribute to the acoplanarity, as we have discussed, but because these partons have such a high momentum the e ects on the acoplanarity are quantitatively very small. The conventional jet shape observable at larger values of r is dominated by very soft particles. In our hybrid model as formulated in this section, these softest particles tend to originate from the last fragmentation events in the shower, as harder partons fragment after they have already left the medium. Soft partons that are produced earlier, in the medium, rapidly lose all their energy. Soft partons at large r that are produced outside the medium cannot be a ected by broadening. For di erent reasons, therefore, both the acoplanarity (dominated by the few hardest partons) and the jet shape at large r (dominated by soft partons which survive because they were produced late, after the shower exits the plasma) are insensitive to transverse momentum broadening.
The key is to focus on the angular distribution of the energy in the jet that is carried by particles in an intermediate interval of transverse momentum. We must focus on semi-hard partons that are su ciently soft that they can be de ected signi cantly by the transverse momentum kicks that they receive from the medium but that are su ciently hard that they survive propagation through the plasma and emerge from it. The new observable that we de ne is, in essence, a more di erential version of the jet shape. Instead of determining the r distribution of all the jet energy, we focus on the distribution of jet energy carried by particles in a certain interval of transverse momentum. And, we construct this observable using only the subleading jets in dijet pairs, speci cally subleading jets with pTsubleading >
30 GeV in a dijet whose leading jet has pTleading > 120 GeV and where the two jets are separated by > 5=6. We make this choice because jets in a sample of subleading jets have on average lost more energy than in a sample of leading or inclusive jets. (We have checked that if we use a sample of leading or inclusive jets, we do see the e ects that we shall describe below but they are smaller in magnitude than in the sample of subleading jets.) We denote our new observable by subleadingpT(r), with
subleadingpT(r)
{ 17 {
JHEP03(2017)135
1
Nsubleading jets
1 r
Xsubleading jets
Pi2 r[notdef] r/2; pTi,track2 range pTi,track
pTjet ; (3.3)
where we take the particles that we include to be only those that have a transverse momentum in a speci ed range. We take this range to be 10 GeV < pT < 20 GeV, a choice whose motivation we describe below. (Note that we rst nd and reconstruct jets using all particles, in the standard fashion, and it is from the fully reconstructed jets we obtain pTjet for the leading and subleading jets. After we have a sample of reconstructed sub-leading jets, we evaluate subleadingpT(r) using only the particles within the speci ed range10 GeV < pT < 20 GeV.) We shall look at the observable subleadingpT(r) for angles up to r < 1. Recent preliminary results regarding somewhat similar jet shape observables, shown in gure 12 of ref. [169], indicate that the measurement of subleadingpT(r) will require background subtraction. We have chosen to follow the procedure used in ref. [169]. We discretize the event in ( ; ) space, where and are distances from the jet axis in rapidity and azimuthal angle, with [notdef] [notdef] < 2:5, in bins of width 0:025 [notdef] 0:025, building a
two-dimensional energy density distribution using only tracks with pT in the desired range. We take reconstruction and background subtraction e ects into account by smearing the jet energy, computing the -dependent average energy density far away from the jet axis in rapidity by summing over all bins with 1:5 < [notdef] [notdef] < 2:5, for each value of , in each event,
and then subtracting this average long-range energy density from the energy density in each of the bins with [notdef] [notdef] < 1:5, for each value of , in each event. After this background
subtraction, we then construct the observable subleadingpT(r) by summing the energy lying in the annulus at a distance r =
p 2 + 2 from the jet. Finally, to maintain consistency with the de nition of the standard jet shape we multiply subleadingpT(r) by the factor
DpTjet=
PNbinsi=0 pTtrack(ri)
JHEP03(2017)135
[angbracketrightBig]
in order to ensure that subleadingpT(r) is normalized to one.
In gure 6 we show the prediction of our hybrid model for the ratio of the new observable subleadingpT(r) in PbPb collisions to that in proton-proton collisions for several values of the broadening parameter K. We show results for subleadingpT(r) with pT in the interval10 GeV < pT < 20 GeV. The selection of these cuts follows from several requirements. We want the partons to be hard with respect to the medium temperature to ensure that no possible collective e ects (including those we will focus on in section 4) can a ect the measurement; this is safely achieved by the lower pT cut of 10 GeV. We also want the partons to be energetic enough to survive propagation through enough plasma that they pick up some transverse momentum kicks, which is to say so that they can exhibit sensitivity to the broadening parameter K. This motivates trying even somewhat larger values of the lower pT cut. From the other side, we want to choose an upper pT cut so that the partons used in the de nition of subleadingpT(r) do have their direction of propagation signi cantly altered by the transverse momentum that they pick up from the plasma. We have found that the range 10 GeV < pT < 20 GeV serves our purposes well. However, we have also investigated subleadingpT(r) with 20 GeV < pT < 40 GeV and this works almost as well. Unlike for the less di erential observables described previously, we see from gure 6 that subleadingpT(r) shows signi cant sensitivity to the value of K. Indeed, as K grows, the partons in this momentum range are more and more kicked out to angles further away from the center of the jet, populating the large r region and depleting the region nearest to the jet axis (r = 0).
{ 18 {
0-10% Centrality
10 < P partonT < 20 GeV
K=100
K=40
K=20
K=0
3
2.5
PbPb/pp
2
1.5
1
JHEP03(2017)135
0.5
0 0.2 0.4 0.6 0.8 1
r
Figure 6. Di erential jet shape ratio constructed from the subleading jets in dijet pairs satisfying pTleading > 120 GeV, pTsubleading > 30 GeV, and > 5=6. The analysis only includes partons
whose pT lies within the intermediate range 10 < pT < 20 GeV. The e ect of broadening on this observable is apparent.
The main features of the results plotted in gure 6 may be understood as coming from the dynamics of broadening and parton energy loss. Let us begin by looking at the K = 0 curve. In the absence of any transverse momentum kicks, the dynamics of energy loss in the hybrid model depletes the number of partons at large angles r. This happens because these partons are produced early and are on average softer, and for both reasons they lose more energy in the plasma. This pushes the PbPb jet shape down at large r. Because the jet shape is normalized, it is pushed up at small r. Since there is no parton energy loss for the jets in pp collisions, the modi cation of the numerator in the jet shape ratio shows up directly in the K = 0 jet shape ratio itself, plotted in gure 6. Now let us ask what happens when we turn on K [negationslash]= 0, adding transverse kicks felt by all the partons in the
jet as they propagate through the medium, namely broadening. As mentioned above, the depletion of partons at small angles r from the jet axis due to broadening results in the reduction of the jet shape at small r in PbPb collisions relative to that in pp collisions with increasing K. Since the shape of jets in vacuum is a rapidly falling function of r, kicking some partons from smaller r to larger r serves to enhance the jet shape in PbPb collisions, and hence pushes the jet shape ratio plotted in gure 6 upwards at larger r, again as a function of increasing K. Furthermore, as K increases partons that are kicked from smaller r can end up at larger and larger values of r, meaning that at larger and larger values of K the increase in the jet shape ratio seen in gure 6 extends further and further to the right. These K-dependent e ects are apparent in gure 6, and they are of course the reason why we have selected and are highlighting this observable.
In experiment, of course, the range in pT employed in the de nition of subleadingpT(r) will need to be a range of momenta of nal state hadrons, not a range of momenta of
{ 19 {
1.6
0-10% Centrality
5 < P hadronT < 10 GeV
K=100
K=40
K=20
K=0
1.4
1.2
PbPb/pp
1
0.8
JHEP03(2017)135
0.6
0.4
0.2
0 0.2 0.4 0.6 0.8 1
r
Figure 7. Di erential jet shape ratio subleadingpT(r) constructed from charged hadrons with 5 GeV< pT <10 GeV, for subleading jets that satisfy the same cuts as in gure 6.
partons. Further investigation of the e ects of hadronization, which are not under good theoretical control, is clearly important. Nevertheless, the observed sensitivity to transverse momentum broadening that we see in gure 6, and that we have not seen in any of the other observables we have studied, encourages us to explore hadronic versions of this observable that are sensitive to the angular distribution of the partons in jets with 10 < pT < 20 GeV
that contribute in gure 6, with the goal of a direct extraction of the broadening parameter K, a key characteristic of the medium. In gure 7 we present the prediction from our hybrid model for a hadronized version of gure 6, namely the observable subleadingpT(r) computed for hadrons rather than partons, with the speci ed range for the pT of the hadrons in the analysis taken to be 5 < pT < 10 GeV. (We use the hadronization prescription described at the end of section 4.2.) We have selected a lower momentum range for the hadrons entering the analysis in gure 7 than for the partons in the analysis in gure 6 simply because hadronization turns partons into softer hadrons. Other choices of momentum range can be investigated.
The main features observed in gure 6 are also observed in gure 7. As in the partonic case, broadening means that as K increases we see a decreasing jet shape ratio at small r and an enhancement in the intermediate r region, with the region where enhancement is seen extending to larger r as K increases further. However, in the hadronic case, the separation between the curves with di erent values of K values is less pronounced than in the partonic case. This is due, in part, to the fact that hadrons within any given range in momentum originate from partons with a wider range of momenta. The behavior of the pT-dependent partonic jet shape subleadingpT(r) is di erent in di erent pT regions, and in particular it becomes less sensitive to K for partons with very large momenta. Since
{ 20 {
some of the hadrons with momenta in the 5 GeV< pT <10 GeV range that we have used to construct the hadronic di erential jet shape ratio in gure 7 come from the hadronization of partons with much higher momentum, this hadronic observable plotted in gure 7 shows less sensitivity to broadening than that seen in gure 6. Because hadronization can turn a single hard parton at very small r into several softer partons that are more spread out in r, it tends to spread the largest changes seen at very small r in gure 6 over a wider range of r in gure 7, reducing their amplitude in the plot. Despite these e ects which serve to dilute the K-dependence seen in the partonic observable, the hadronic observable in gure 7 displays sensitivity to transverse momentum broadening dynamics.
Although it is of course necessary to do further investigations of hadronization, it seems promising that by selecting hadrons from subleading jets in an intermediate pT-
range, 5 GeV< pT <10 GeV in gure 7, the e ects of the transverse momentum kicks on the partons in the jet on the observable are enhanced, and are considerable, relative to what we have seen for less di erential observables such as the usual jet shapes in gure 5. The further investigation of observables that, like the pT-di erential jet shape ratio we have constructed from subleading jets, focus on the transverse dynamics of partons in the jet that are neither very soft (as in that regime other e ects that we investigate in the next section dominate) nor very hard (as in that regime the e ects of transverse kicks received from the medium are hard to see) holds considerable promise. Such observables represent the only path that we can see at present toward an experimental determination of the value of K directly from its de nition in terms of momentum broadening rather than indirectly via inference from measurements of energy loss. And, K = ^
q=T 3 is a key parameter in our model or in any model as it quanti es a central property of the medium that is related to how strongly coupled the uid is. We may not yet have hit upon the precise de nition of the optimal observable, as doing so requires balancing choices of jet selection, pT cuts, jet shape measure, and background subtraction in the face of potentially competing goals: maximizing sensitivity to K while at the same time optimizing the statistics, signal-to-noise, and utility of the analysis of any speci c data set. We look forward to seeing this done in consultation between experimentalists and theorists. Our results in gure 7 provide strong motivation for doing so.
4 Medium backreaction: observable consequences of a wake in the plasma
In this section we will take the rst steps toward incorporating another important feature of jet quenching dynamics in plasma into our hybrid model: the conservation of the momentum and energy lost by the jet. In our implementation of the model up to now, we have assumed that the energy lost by the energetic partons in the jet thermalizes to such a complete extent that, after hadronization, it becomes particles moving in random directions that are completely uncorrelated with the jet direction. The assumption of rapid thermalization is well motivated by the many lines of evidence indicating that quark-gluon plasma is a strongly coupled liquid. It is also motivated by the holographic analysis of the stress tensor that describes the response of the N = 4 SYM uid to the passage of
an energetic parton through it [123, 125] which showed that after a short time of order
{ 21 {
JHEP03(2017)135
1=T all the energy dumped into the medium was rapidly converted into a hydrodynamic excitations of the system. The possible hydrodynamic collective response of the plasma to the passage of an energetic parton through it was characterized earlier in refs. [171{ 174]; the holographic computation in refs. [123, 125] provided a concrete realization in a full quantum eld theoretical calculation, con rming that the anticipated hydrodynamic response could indeed be excited via the passage through the uid of an energetic probe whose size is much smaller than 1=T . This notwithstanding, it cannot be correct to assume that the energy lost from the jet, after thermalization, has no memory of the direction of motion of the jet. After all, momentum is conserved and the jet loses momentum to the plasma, as well as energy. And indeed, the hydrodynamic response of the plasma includes both a wake and sound waves (a Mach cone) that carry momentum in the direction of the jet, as well as energy. In this section, we will provide an initial account of this collective response by providing a simpli ed description of the backreaction of the medium to the passage of the jet that respects energy and momentum conservation without introducing any additional parameters into our hybrid model. We leave a full treatment (which would involve sourcing, propagating, and hadronizing a hydrodynamic wake in the expanding cooling hydrodynamic uid) to the future.
The mechanism we will describe is not exclusive to strong coupling, although this scenario provides a completely natural realization for this collective response. Even in perturbative QCD analyses in which a hard parton loses energy by radiating gluons the radiated gluons themselves interact with the medium and radiate further softer gluons, with the result being a rapid degradation of the emitted gluon momenta. This degradation may be viewed as a rapid transfer of energy from hard (jet) modes to soft (medium) modes and it has been described recently in detail in refs. [150{156]. In an in nite medium, this process leads to the eventual (local) thermalization of the lost energy after which the dynamics may be best understood in terms of hydrodynamics. In a nite medium, the local thermalization process may not be complete and it is possible that the medium modes excited by the jet may remember more than just the energy and momentum that the jet lost and they gained, but keeping track of their momentum and energy as we shall do is certainly a good rst step. One circumstance in which our analysis in this section would not be relevant is for a medium that is su ciently thin and su ciently weakly coupled that gluons radiated by the hard parton do not interact after they are radiated, but this is an extreme scenario.
Following this discussion, in section 4.1 we will describe a simple implementation of the backreaction of the plasma to the passage of a jet of energetic partons through it. We will refer to the backreaction of the plasma in generic terms as the wake of the jet, as we shall not need to focus on the distinction between a di usive wake (moving and perhaps heated uid) and propagating sound waves since upon making the approximations that we describe below the perturbations to the spectrum of hadrons in the nal state that we compute cannot distinguish between di usive and sound modes. We will characterize the medium response to the transfer of momentum to it from the jet, that is the wake, by analyzing the induced velocity and temperature variations of the hydrodynamic behavior of the quark-gluon plasma uid, upon making several simplifying assumptions. Since the
{ 22 {
JHEP03(2017)135
total energy deposited in the medium by a typical jet in heavy ion collisions (O 20 GeV)
is small when compared to the total energy per unit rapidity in the event (O 1 TeV),
as in refs. [171{174] we shall treat the additional momentum acquired by the medium as a consequence of the passage of the jet as a small perturbation. (For studies of the nonlinear response of the plasma, see refs. [158, 161, 175, 176]). We shall also assume that the small perturbation of the velocity and temperature of the medium translates into small perturbations to the resulting distribution of particles at all momenta. These physical assumptions can be described somewhat loosely by saying that we are assuming that the energy deposited into the medium as a wake thermalizes to the maximum degree allowed by conservation of momentum and energy, turning into a perturbation to the spectra of the hadrons in the nal state that remembers the energy and momentum deposited by the jet into the medium but nothing else about its origin. As we have discussed, this is natural at strong coupling, or at weak coupling if the gluons radiated by the hard partons in the jet themselves radiate many times further, but would not be valid if the energy lost by the jet is carried by only a few particles.
For simplicity, we will also assume that the unperturbed uid is well described by a boost invariant ow. And, although in our simulation of the amount of energy that the jet loses we do take the transverse ow of the medium into account and we do use a background obtained via solving viscous hydrodynamics, in this initial study of the wake we will for simplicity neglect the transverse ow of the unperturbed uid and employ ideal hydrodynamics. All that said, our aim will actually not be a description of the perturbations of the hydrodynamics per se: we wish to focus instead directly on the modi cation to the spectrum of hadrons produced after the perturbed hydrodynamic uid freezes out. The last simplifying assumption that we make is that the plasma freezes out along a constant proper time hypersurface.
After making these various approximations, in section 4.1 we will derive a simple expression for the modi cation to the spectrum of hadrons formed as the hydrodynamic uid perturbed by the passage of a jet through it freezes out, an expression that is determined solely by the amount of energy and momentum lost by the jet. In subsequent subsections, we will describe how we implement this expression in our hybrid model, and look at its consequences for a number of jet observables.
It is inevitable that when a jet reconstruction algorithm is used to nd and reconstruct jets in heavy ion collisions, some of the energy and momentum that is counted as part of a jet in fact comes from hadrons formed from the plasma as it freezes out, given that the plasma includes a moving heated wake that, by momentum conservation, is owing in the same direction as the jet [95, 158, 160{162]. Furthermore, since any background subtraction procedure involves comparing events with a jet or jets to events that do not contain jets, and since events that do not contain jets also do not contain wakes, the particles from the hadronization of the wake that end up reconstructed as part of a jet will not be removed by background subtraction. Since there is no way for experimentalists to remove them from the jets they reconstruct, theorists must add them to the jets in their calculations. This is our goal in this section.
{ 23 {
JHEP03(2017)135
4.1 The spectrum of hadrons from a medium perturbed by the passage of a jet
In a boost invariant uid with no transverse velocity, the wake associated with the passage of a jet may be characterized by a perturbation to the velocity eld with the form
u = 0; ui; u
[parenrightbig]
; (4.1)
where ui with i = 1; 2 and u are the variations of the velocity eld in the transverse plane and in the space-time rapidity direction. The disturbance will also be associated with a change in the temperature of the plasma T . These perturbations are functions of time and of all the space coordinates, and need not be boost invariant. (Note that because u + u is normalized, the u component in (4.1) is xed by, and quadratic in, ui and u and we have therefore neglected it.) At a xed proper time , the total momentum stored in the perturbation of the plasma is
P i? = w [integraldisplay]
d2x? d ui? ; P = w [integraldisplay]
d2x? d u ; (4.2)
where we have used the fact that, to leading order in the perturbation, the variation of the stress tensor of ideal hydrodynamics takes the form T a = w ua for a = 1; 2; with w = "+P is the enthalpy of the unperturbed uid, which is related to the entropy density of the unperturbed uid through w = T s and which is constant on a xed- surface in a uid which is boost invariant and has no transverse expansion. Because we are simplifying this analysis by assuming ideal hydrodynamics, we are neglecting any production of entropy during the hydrodynamic evolution, for example as the sound waves excited in the plasma are damped by viscosity. The only entropy dumped into the plasma is that dumped into the wake initially by the jet itself, and this entropy is then equal to the entropy associated with the perturbation computed at late time by computing the ux of the entropy current across a xed hyper-surface
S = s
c2s
[integraldisplay]
T
T ; (4.3)
where cs is the speed of sound of the unperturbed uid which, within the boost invariant assumption, only depends on the proper time .
Both eqs. (4.2) and (4.3) are valid on any xed proper time hypersurface including, in particular, on the freeze-out hypersurface. These expressions describe all the momentum and entropy that was dumped into the plasma by the disturbance, as they are at the freeze-out time. Following our assumption that all the momentum lost by the jet is incorporated into the plasma as a wake, P i? and P are given by the transverse momentum and
momentum rapidity lost by the jet as it traverses the plasma, respectively. Since in our implementation of energy loss we have assumed that the rapidity of the jet remains constant (or almost constant if we turn on transverse momentum broadening) we set P = 0 and
use our hybrid model calculation to give us the P i? injected into the hydrodynamic uid by
each of the jets we analyze. We shall see below that the entropy production S associated with the jet passage through the medium may also be constrained by the total energy carried by the particles produced by the perturbation.
{ 24 {
JHEP03(2017)135
d d2x?
Our goal is to determine the e ect of the perturbations ua and T that describe the wake in the uid that the passage of the jet creates on the spectrum of particles produced when the uid freezes out. To do so, we will employ the standard Cooper-Frye prescription [177]
E dN
d3p =
1 (2)3
[integraldisplay]
d p f(up) ; (4.4)
where the integral is over the freeze-out hypersurface and where for simplicity we shall assume a Boltzmann distribution f(E) = exp (E=T ). Expanding to leading order in the
perturbation, we obtain the expression
E d N
d3p =
JHEP03(2017)135
(2)3
[integraldisplay]
mT
T cosh(y )
d2x? d mT cosh(y ) exp [bracketleftBig]
[bracketrightBig]
(4.5)
[notdef]
(pi? ui?
T + 2p u
T +mTT TT cosh(y )
);
qm2 + p2T is the transverse mass of the emitted thermal particle. Note that this perturbative expression is only valid for particles whose momenta p are comparable to the temperature T on the freeze-out surface, where the small perturbations in the hydrodynamic quantities results in a small perturbation on the resulting distribution of particles. For momenta far above T , where the thermal distributions are exponentially small, the perturbations in the hydrodynamic quantities can have large relative e ects; in this regime, the expression (4.5) is not valid and, because it is based on a linear approximation to an exponential, it in fact underestimates the particle production from the wake. However, this regime of the spectrum contributes little to overall yields.
The expression in eq. (4.5) is general and independent of the space-time dependence of the perturbed hydrodynamic elds. To proceed further we will assume that during the space-time evolution of the perturbation over the boost invariant background, the space-time rapidity of the disturbance remains approximately constant. Since high-energy jets propagate at a xed space-time rapidity j equal to their momentum rapidity yj, j = yj, this assumption implies that the perturbation is narrow around the momentum rapidity of the jet, which allows us to perform the integration in eq. (4.5) by replacing ! yj. We perform this integration, use eqs. (4.2) and (4.3) to relate the three terms in eq. (4.5) to P i?, P (which vanishes) and S, and then impose that the energy and momentum of
the emitted particles equals the energy and momentum lost by the jet
E = [integraldisplay]
where mT
d3p d Nd3p p?,i ; (4.6)
where the integration of the rst of these equations leads to the relation S = E=(T cosh yj). After these manipulations, we are able to express the spectrum of particles eq. (4.5) emitted from the boosted, heated up, wake in the uid as
E d N
d3p =
1 32
d3p d Nd3p E ; P?,i =
[integraldisplay]
mT
T 5 cosh(y yj) exp
hmTT cosh(y yj) [bracketrightBig]
[notdef]
(p? P? cos( j) +13mT MT cosh(y yj)
):(4.7)
{ 25 {
where pT , mT , and y are the transverse momentum, transverse mass, azimuthal angle and rapidity of the emitted thermal particles whose distribution we have obtained, and where PT and MT = E= cosh yj are the transverse momentum and transverse mass transferred from the jet (whose azimuthal angle and rapidity are j and yj) to the wake in the uid. Note that the distribution (4.7) is a small correction that must be added to the one-body distribution of particles that the unperturbed hot plasma would have emitted. In particular, the distribution (4.7) may be negative. For example, this occurs for particles emitted in the direction opposite to the direction of the jet. Negative values simply mean that the perturbed thermal uid emits less particles in the direction opposite to the direction in which the jet was propagating than the unperturbed uid would have. This is a direct consequence of the fact that the jet loses momentum to the uid, exciting a wake of uid moving with net momentum along the direction of the jet.
The closed form expression (4.7), which only depends on the momentum lost by the jets in the plasma and on the kinematics of the jet, will be the basis for our analysis of the observable consequences of the wake in the plasma, which is to say of the backreaction of the medium. Subject to the assumptions employed in its derivation, the spectrum (4.7) will allow us to provide estimates of the observable e ects of the collective response of the plasma to the passage of the jet through it without having to model the complicated processes (pre-hydrodynamic and hydrodynamic) via which the energy lost by the jet relaxes. As stated above, even within the assumptions employed in its derivation the expression (4.7) is only valid for particles emitted with a momentum comparable to the mean thermal momentum in a uid cell at freeze-out; it need not be valid for semi-hard particles produced in the plasma. Perhaps the most important assumption in its derivation is the assumption that the wake in the plasma is fully thermalized, with the only aspects of its origins that it remembers being its total energy and momentum. This too need not be valid for semi-hard particles, some of which will originate from the energy lost by the jet either near the edge of the plasma or shortly before freezeout that do not thermalize. Nevertheless, we will use the spectrum (4.7) at all momenta to obtain rst estimates of the observable consequences of the presence of a wake in the plasma in events with jets. In the next subsection, we shall explain how we have implemented this spectrum in our hybrid model analysis of jets propagating within hot QCD plasma.
4.2 Implementation of backreaction, background model, and jet hadronization
The implementation of the simpli ed expression (4.7) for the spectrum of particles resulting from the wake that is the backreaction of the medium to the presence of the jet demands further modelling for a proper description in heavy ion collisions. Three aspects that we will need to incorporate into our description are: (i) the e ect of the radial ow and chemical composition of the unperturbed reball on the particles resulting from the backreaction perturbation; (ii) adding a background of particles coming from the freeze-out of the unperturbed reball to our hybrid model in order to be able to properly account for the negative contribution from the perturbation (4.7); and (iii) the generation of particles consistent with the one-body distribution eq. (4.7). Furthermore, since the particles produced after decoupling are hadrons, we will also need to consider the hadronization of our
{ 26 {
JHEP03(2017)135
quenched jets, which we have not needed to include in our previous implementation of the hybrid model. We will discuss these four aspects sequentially in this subsection.
As we have stated, our derivation of eq. (4.7) neglects the e ect of transverse ow in the unperturbed uid, meaning that we have neglected radial ow, as well as elliptic ow and higher azimuthal harmonics of the transverse ow. The e ects of elliptic ow and higher harmonics are small in the most central collisions, but radial ow cannot be neglected as it has important consequences for the spectrum of particles produced by the reball. The radial boost in the spectrum of particles due to the radial ow in the uid from which they are formed yields a blue-shifted spectrum which is harder than in the absence of the boost. Since heavier particles pick up more momentum than lighter particles when all are boosted by the same velocity, another consequence of radial ow is di ering spectra for particles with di erent masses, with heavier hadrons getting harder spectra than lighter hadrons. The result (4.7) was derived as a perturbation on a background that does not include either e ect. As a crude step towards including both e ects, we will employ the spectrum eq. (4.7) but instead of using the temperature T at the time of freezeout we will use species-dependent, momentum-dependent, empirical expressions for T that provide a good description of the measured particle spectra in ref. [178] upon tting the measured spectra to \thermal" spectra without radial ow. Speci cally, we assume a proton to pion ratio of 0:05, neglect hadrons other than protons and pions, and use the following momentum-dependent \temperatures" in eq. (4.7):
T(pT) =
0:15 GeV if pT < 0:07 GeV
0:33 pT
GeV
These empirical expressions provide a good description of both the proton and pion spectra in central heavy ion collisions with ps = 2:76 TeV at the LHC for pT < 3 GeV. (The largest deviations, around 10%, occur for pT [lessorsimilar] 0:5 GeV, where the measured pion yields contain large contributions from resonance decays.) We shall use them in eq. (4.7) for pT < 5 GeV.
This approach is admittedly crude, but it is of value in this exploratory investigation as it allows us to use the closed form expression (4.7) rather than attempting a full hydrodynamic calculation of the wake induced by each of the jets in our ensemble of events.
As we have already discussed, a characteristic feature of eq. (4.7) is that the spectrum of particles coming from the perturbation, from the wake in the uid, can become negative at large azimuthal angles with respect to the jet direction, particulary for particles with small mT or for jets with small MT . This re ects the fact that the wake is made up of uid moving with a net momentum in the jet direction, with the negative contribution to the perturbation representing a depletion of the momentum in directions opposite to that of the jet. This means that in order to implement eq. (4.7) in our hybrid model, we need to embed our jet sample in a background whose purpose is to provide su cient thermal particles such that where the perturbation (4.7) is negative we have some thermal
{ 27 {
JHEP03(2017)135
(0:19 GeV if pT < 0:7 GeV 0:21 pT
GeV
0.28 GeV if pT > 0:7 GeV
(4.8)
Tp(pT) =
8
>
>
<
>
>
:
0.3 GeV if 0:07 GeV < pT < 1:9 GeV 0:4 GeV if pT > 1:9 GeV
(4.9)
particles that we can remove. In previous implementations of our hybrid model in which we neglected the back reaction of the medium in response to the presence of a jet we ignored the background on the basis that in experimental analyses of jet data it would be subtracted. Contributions from the wake cannot be subtracted, which is why we are now adding them. Since these perturbations can be negative, we now need a baseline background too. The background that we use is oversimpli ed. It is constructed by generating an ensemble of pions and protons which is at in and that reproduces the measured particle yields and spectra [178]. The addition of this contribution will force us to introduce a background subtraction procedure in our analysis of in-medium jets. We will describe this procedure in detail in subsection 4.3 and appendix B
For each event in our sample, we determine the momentum lost by each of the partons in the jet shower as well as in initial state radiation following our hybrid strong/weak coupling model as described in section 2. (For the backreaction analysis, we include all partons in the jet shower and in initial state radiation, whether or not they end up being reconstructed in a jet.) Since each of the propagating partons loses energy into the medium, each of them generates its own wake induced by its lost momentum. In the linearized approximation that we have employed, the multiple wakes do not alter each other, and the nal spectrum is the superposition of the spectra generated by the wake of each propagating parton. At hadronization, each of the induced wakes generates an ensemble of particles with the one body distribution (4.7) computed from the kinematics of each parton and its lost momentum.
In order to incorporate the e ect of uctuations in the reconstructed jets, we generate the ensemble of particles coming from the wake in the medium via a simple Metropolis algorithm designed to satisfy the conservation of the lost jet energy while drawing particles from the distribution (4.7). First we generate an independent list of particles, including protons with 5% probability and pions with 95% probability, from the one body distribution (4.7), until the sum of their energies reaches the lost energy. This ensemble generally contains particles in the region of azimuthal angle in which eq. (4.7) becomes negative, which we will call negative particles. Whenever a negative particle is produced, its contribution to the net energy and momentum of the ensemble is negative. We later neutralize these negative particles by removing a particle from the background which is su ciently close in ( ; ) space and in transverse energy.2 From this initial ensemble, whose four-momentum in general will not coincide with the momentum lost by the jet, we randomly select a particle which we replace by a new particle drawn from the distribution (4.7). If
2Negative particle neutralization proceeds as follows. We identify the background particles within an angle of r = 0.3 in ( , ) space of the negative particle we wish to neutralize. Among the candidate background particles, we choose the best candidate in terms of ET and angular position in plane,
which means we minimize ET and r, the di erence in transverse energy and angular position with respect to that of the negative particle, respectively. We do the minimization by starting with any one of the candidate background particles and then checking each one of the other candidates to see whether choosing it instead reduces ET while not increasing r by more than 0.05 or reduces r while not increasing ET by more than 0.1 GeV. We have checked that after using this procedure we end up with 90% of the negative particles neutralized to better than 0.1 GeV in transverse energy via subtraction of a background particle that was within r = 0.13 of it.
{ 28 {
JHEP03(2017)135
the change improves four-momentum conservation, it is accepted. Otherwise, the change may be accepted with a probability distribution
e(pnew ensemble P )2
e(pensemble P
)2 ; (4.10)
where pnew ensemble is the four momentum of the candidate ensemble including the newly drawn particle, pensemble the four-momentum of the previous ensemble, and P is the four-momentum lost by the jet. We repeat the procedure until each of the four components of pensemble is within 0:4 GeV of the momentum lost by the jet. (We have explicitly
checked that changing this threshold does not signi cantly change our results.) The ensemble generated after this procedure conserves energy and momentum (within the tolerance above) and possesses a one-body distribution identical to eq. (4.7), something that we have checked explicitly.
Since the medium response manifests itself in the form of modi ed distributions for hadrons, namely pions and protons in our approach, we are forced to consider the hadronization of jet showers in order to properly incorporate the particles from the jet itself and those from its wake in each event before then reconstructing jets and calculating observables. Hadronization leads, generically, to a softening of the typical jet fragments. Because of its nonperturbative nature, hadronization even in vacuum remains a fundamental problem and presents serious challenges to phenomenological modelling. Furthermore, it is not presently understood how the phenomenological hadronization models that have been applied successfully to QCD processes in vacuum should be modi ed due to the presence of a heavy ion environment. As an example, changes in how color ows in the jet shower resulting from soft exchanges between partons in the shower and the medium lead to signi cant modi cations to subsequent hadronization in certain regions of phase space [179, 180], but the overall description of these e ects remains to be determined. Because of all these uncertainties, in this work we will adopt a simpli ed model for the in-medium hadronization, as in much of the literature. We simply assume that hadronization of high energy jets occurs in heavy ion collisions as in vacuum, in particular keeping the same color correlations between partons in the shower even though we know that in reality these must change as the partons interact with the medium. Although several di erent prescriptions for hadronization in vacuum exist, in this work we will employ the Lund string model as implemented in Pythia, feeding the showers to this hadronization model after they have been modi ed according to our hybrid approach to energy loss and broadening. We defer comparisons between di erent models for hadronization to future work.
4.3 The e ect of backreaction on jet observables
After the implementation of the energy loss su ered by the partons in jet showers produced in hard processes as they propagate through the hydrodynamic medium via the hybrid model described in section 2, the incorporation of the e ects of transverse momentum broadening described in section 3 if we choose K [negationslash]= 0, the incorporation of a thermal-like
background of particles and the perturbation to that background corresponding to the
{ 29 {
W (pnew ensemble) =
JHEP03(2017)135
e ects of the wake in the hydrodynamic medium as it responds to the passage of the jet as described in subsection 4.2, we now have a full event simulation from which to extract calculations of medium-modi ed jet observables. This will be our goal in this and the next subsections. In contrast to our previous publications, the inclusion of a background forces us to implement a background subtraction procedure, making our analyses of quenched events more similar to the actual analysis of jet data at the LHC. The fact that the medium responds to the passage of the jet, and in particular the fact that the medium picks up momentum in the jet direction, makes this complicated procedure absolutely necessary.
For hard jets produced together with a soft background that is completely uncorrelated with the jet direction, there are a number of established techniques that allow for systematic removal of the e ects of background particles from jet observables. (See, for example, refs. [181{183].) These procedures, generically referred to as background subtraction, are routinely applied to jet measurements at the LHC and, at least in proton-proton collisions, they e ciently remove the e ects of soft (non-perturbative) backgrounds that may be large but that are uncorrelated with the jet, allowing the measurement of theoretically controlled hard processes. However, in heavy ion collisions the fact that the medium includes a wake that carries momentum in the jet direction means, in e ect, that a component of the background is correlated with the jet direction. This makes it impossible for a background subtraction procedure to separate the jet (which has been modi ed, via energy loss and broadening) from the medium (which has been modi ed, via the wake). In order to compare to experimental measurements, therefore, we have added a background and a wake and must now perform a background subtraction as if the background were uncorrelated with the jet direction, followed by jet reconstruction, just as in an experimental analysis. This procedure is not necessary for jet observables that are dominated by the harder components of a jet. This procedure is important for the softer components, since the softer components of what is reconstructed as a jet will include contributions from the jet itself and from the wake in the medium. In particular, this procedure is critical to gauging the e ects of the wake on observables. We have implemented a full background subtraction procedure to analyze the events produced within our framework. In particular, we have implemented a version of the so called noise/pedestal background subtraction procedure [48, 184] and then done a jet energy scale correction; the details of our implementation can be found in appendix B.
As in section 3, we rst consider the R-dependence of jet suppression. In gure 8 we show the jet RAA as a function of pT for central events for a wide range of jet momenta and for di erent anti-kt reconstruction parameters R. We set K = 0, neglecting transverse momentum broadening since, as we saw in section 3, it has little e ect on these observables. In the right panel, we show our results for 70 < pT < 300 GeV in comparison with CMS measurements of RAA for R = 0:3 jets in the range 100 < pT < 300 GeV. Note that, as when we included broadening in section 3, when we incorporate the e ect of the backre-action of the medium to the jet this alters the jet suppression, meaning that we had to retune the energy loss parameter sc in our hybrid model. In this case, we only needed to modify the value of sc at the percent level, which is very much smaller than the theoretical
{ 30 {
JHEP03(2017)135
1.2
1.2
0-10% Centrality
R=0.2 R=0.3 R=0.4 R=0.5 ATLAS Data R=0.4
R=0.2 R=0.3 R=0.4 R=0.5 CMS Data R=0.3
0-10% Centrality
1
1
0.8
0.8
JetR AA
0.6
JetR AA
0.6
0.4
0.4
0.2
0.2
0 20 40 60 80 100 120
0 60 80 100 120 140 160 180 200 220 240 260 280 300
PT (GeV)
PT (GeV)
JHEP03(2017)135
Figure 8. Dependence of jet RAA on the anti-kt parameter R used in reconstructing the jet and on the transverse momentum of the jet. We have set K = 0, turning o transverse momentum broadening which we saw in section 3 has only very small e ects on this observable. We have included the e ects of the wake in the hydrodynamic uid, which is to say the backreaction of the medium to the presence of the jet, and have added a background and implemented two di erent background subtraction procedures as described in the text and in appendix B. The left plot shows this observable for jets with 40 < pTjet < 120 GeV as measured by the ATLAS collaboration [56], with their background subtraction procedure, while the right one is for jets with 70 < pTjet <
300 GeV as measured by the CMS collaboration [49], with their background subtraction procedure. Both panels are for the 0-10% most central collisions with ps = 2:76 ATeV. We have extended the pT-range in the right panel down to 70 GeV even though present CMS measurements are for jets in the range 100 < pTjet < 300 GeV to make it possible to compare the results from our model in the pT-range where the two panels overlap in order to see the e ect on RAA of choosing between the two di erent background subtraction procedures. The di erence between the two panels is small, but visible.
uncertainty corresponding to the widths of the bands in all our plots in this paper and in our previous publications. As we found when we included broadening in section 3, for the high energy jets in the right panel of gure 8 the suppression factor RAA shows only a very small dependence on R, consistent with LHC data [49].
As we saw in section 3, the suppression factor RAA shows a small decrease (i.e. increase in suppression) with increasing R, corresponding to the fact that with increasing R the angular size of the jets that are reconstructed increases together with the fact that wider jets lose more energy. However, this e ect is milder here than it was in gure 2 because the e ect of the backreaction of the medium that we are incorporating here is that some particles coming from the wake in the plasma, which is to say some of the energy that the jet lost, ends up reconstructed as if it were still part of the jet. Nevertheless, the wide angular distribution of the spectrum of particles from the wake given by eq. (4.7) implies that, even for the relatively large value of R = 0:5 explored in gure 8, the fraction of the energy lost by the jet that is recovered by the reconstruction procedure is small. For even larger values of the reconstruction parameter R, the full jet energy would in principle be recovered and RAA would approach unity. However, for such large values of R the uctuations in the background make the reconstruction algorithm we have implemented unreliable. (We have tested the recovery of jet energy by embedding our simulations in a homogeneous background with no uctuations.)
{ 31 {
1.6
Strong Coupling R=0.3 Strong Coupling R=0.4 Strong Coupling R=0.5
PYTHIA R=0.3 PYTHIA R=0.4 PYTHIA R=0.5
1.4
dN
dP Td (R)
0-10% Centrality
| | < 2.1
1.2
1
dN
dP Td (R=0.2)/
0.8
0.6
0.4
JHEP03(2017)135
0.2
0
40 60 80 100
Jet PT (GeV)
Figure 9. Ratio of the spectra of jets reconstructed with the anti-kt parameter R = 0:2 to the spectrum of jets reconstructed with various larger values of R as a function of the reconstructed jet pt. The colored bands correspond to the prediction of our hybrid model, with no broadening (K = 0) and including the e ects of the backreaction of the medium. The dots are the predictions for jets in vacuum, in proton-proton collisions as described by Pythia. As explained in the text, the fact that the three colored bands are closer together than the three sets of dots are is due to the fact that wider jets lose more energy than narrower jets do.
In the left panel of gure 8, we show the jet RAA for smaller jet pT for the di erent values of R. The calculated suppression factor shows a mild increase with pT. It is consistent with all but the lowest pT ATLAS data point within the theoretical and experimental uncertainties, although the data themselves dont show evidence for such an increase. In spite of the lower jet energy, the e ect of including the particles coming from the backre-action of the medium on RAA is small for all the values of R we have explored, meaning that the R-dependence of our results for jet RAA is similarly small to that seen in the right panel of the gure. This is in apparent contrast with RCP measurements reported by ATLAS in ref. [48], which seem to indicate that the suppression of jets decreases with R in this range of momenta, albeit with sizeable systematic uncertainties. Since ATLAS reports the R-dependence of RCP (which can be thought of as the ratio of RAA for central events to that for peripheral events) rather than of RAA, and since our simulations do not reproduce RAA for peripheral events well (since the small amount of energy lost in the hadronic phase, which we do not include in our hybrid model, is proportionally more relevant in peripheral events where the total energy lost is smaller [148, 149]) we will not attempt to make any quantitative comparisons with these data. We look forward to anticipated further measurements of the R-dependence of RAA itself, where a direct comparison will be possible.
To further study the R-dependence of jet suppression in our model, and motivated by the ALICE analysis reported in ref. [58], in gure 9 we show the ratio of spectra of jets in PbPb collisions reconstructed with di erent values of the anti-kt parameter R to that
{ 32 {
for R = 0:2 jets. For comparison, we show the same ratio in proton-proton collisions as simulated by Pythia. An important advantage of this observable is that it is constructed with PbPb data only. This leads to a signi cant reduction in the theoretical uncertainties in our model, since the errors in the spectra with di erent R are correlated. This is the reason why the widths of the colored bands depicting our theoretical predictions displayed in gure 9 are signi cantly narrower than those for our predictions of RAA. Because of this reduction in the theoretical uncertainty, as we shall now explain this observable is su ciently discriminating to show that, within the range of R studied, the wider jets that are reconstructed with larger values of R lose more energy than narrower jets. To see this, let us rst understand the behavior of this observable for jets produced in proton-proton collisions, which evolve in vacuum. Predictions from Pythia for these jets are shown as the colored dots in gure 9. In vacuum, the number of jets with a given pT increases as R increases, since reconstructing wider jets incorporates a larger fraction of the initial partonic energy into the jet. Therefore, the ratio of the spectra of jets reconstructed with R = 0:2 to that of jets reconstructed with a given R decreases with increasing R. In the medium produced in PbPb collisions, this general trend is also observed. However, this ratio decreases more slowly with increasing R, meaning that the number of jets with a given pT now increases more slowly with increasing R than was the case in vacuum. This means that the wider jets reconstructed with larger R have lost more energy than the narrower jets have.
Another interesting feature of this observable is that the deviation between the vacuum and medium ratios increases as the momentum of the jet increases. This, too, is a consequence of the fact that wider jets which contain more in-medium partons lose more energy than narrower jets. Low pT jets contain, on average, a smaller number of partons propagating simultaneously in medium. The extreme case is a jet that consists only of a single parton while it is in the medium, although it may branch later. In this case, the energy loss is independent of the reconstruction parameter R: if this jet is reconstructed, its energy loss is the energy lost by that single parton no matter what the value of R. The fact that the three colored bands in gure 9 come closer together as pT increases indicates that the lowest pT jets are dominated by jets that contain only very few partons while they are in the medium. As the pT increases, the population of jets includes some that are wider, with more partons, that lose more energy. And, at any given pT more of these are reconstructed at larger R. Although a direct comparison with the semi-inclusive jet measurements performed by ALICE [58] that motivated us to make gure 9 is not possible, our analysis highlights the potential that precise measurements of this type will have to discern the mechanisms of jet quenching in future.
Paralleling the discussion in section 3, to explore the e ects of medium-modi cation of jets on their angular structure we now turn to the jet shapes ratio, shown in gure 10 for two di erent centralities. We set K = 0 since, as we saw in section 3, transverse momentum broadening has little e ect on this jet shapes ratio. We describe the way in which we subtract the thermal background in appendix B. Since the uctuations in our simpli ed background do not coincide with those of an actual heavy ion collision, we need to correct for the di erence in jet energy resolution in order to do a fair comparison with CMS data [51].
{ 33 {
JHEP03(2017)135
1.8
1.8
0-10%
Backreaction No Backreaction
CMS Data
1.6
1.6
100 < P jetT < 300 GeV 0.3 < | | < 2, r < 0.3
P hadron
T > 1 GeV
10-30%
1.4
1.4
PbPb/pp
PbPb/pp
1.2
1.2
1
1
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0 0.05 0.1 0.15 0.2 0.25 0.3
0 0.05 0.1 0.15 0.2 0.25 0.3
r
r
Figure 10. Ratio of the jet shape in PbPb collisions with ps = 2:76 ATeV with 0-10% centrality (left) and 10-30% centrality (right) to the jet shape in proton-proton collisions. The two colored bands show the results of our hybrid model calculation with no broadening, with both jets and background hadronized, and with our background subtraction procedure for high-pT jets applied.
In the calculation shown as the red band we include the e ects of backreaction, namely the particles coming from a wake in the medium. We compare our calculation with and without backreaction to data from CMS [51].
This amounts to smearing the jet energies with a Gaussian whose width corresponds to the di erence between the jet energy resolution in the presence of our background and the jet energy resolution measured by CMS; we describe the procedure in appendix B. Last, we subtract background tracks in the jet cone following a simple procedure from ref. [51] in which we subtract the -re ection of each event from that event. This procedure does not work for jets near = 0; this is why [notdef] [notdef] < 0:3 is excluded from both our analysis and
the measurement reported in [51].
To gauge the e ects of adding our simpli ed background, performing the background subtraction procedure, and hadronization on one hand, and the e ects due to the backreaction of the medium, namely the particles coming from the wake in the plasma, on the other in both panels we show the jet shape ratio computed at the hadronic level with and without backreaction. As we saw in section 4, energy loss serves to narrow the angular size of jets in a given window of energies in heavy ion collisions relative to that of jets with the same energies in proton-proton collisions. As a consequence, without backreaction the e ect of energy loss is to increase the importance of narrow jets in the quenched jet sample, leading to a depletion of the jet shape at large angles r. Note that the only di erences between the simulations without backreaction in gure 10 and the K = 0 simulations displayed in gure 5 are: adding the simpli ed but uctuating background that we are employing, performing our background subtraction and jet reconstruction, and adding hadronization. The partonic distributions whose ratio is plotted in gure 5 give rise to narrower distributions that the hadronic ones that go into gure 10, a natural consequence of the non-trivial angular distribution of the Lund strings connecting the hard partons within the jet which means that hadronization broadens the jet somewhat. (See for example ref. [185].)
Despite the hadronic uncertainties, the jet shape ratio shows a clear increase at larger values of the angular variable r when we include backreaction, con rming the expectation that some of the particles from the wake in the plasma do end up reconstructed as part
{ 34 {
JHEP03(2017)135
of the jet, and con rming the expectation that they are less tightly focused in angle than the jet itself was. That said, it is also clear from gure 10 that the data features a much stronger increase in the jet shape ratio at large angular distance r from the jet axis than we obtain from our calculation. (Although the e ect is large when plotted in this way, it is important to keep in mind that the enhancement in the ratio of jet shapes at large r does not mean that in-medium jets are signi cantly wider than vacuum jets: in both populations, most of the jet energy is concentrated at small r; at large r what is plotted is the ratio of two quantities that are both small. The data indicate that in PbPb jets the fraction of the jet energy at r > 0:2 is larger than in proton-proton jets, but both are small.) Comparing the measured jet shape ratio to our calculations in gure 10 tells us that our treatment of backreaction substantially underestimates the amount of energy that ends up correlated with the jet axis but separated from it by a relatively large value of r. The hadronization of the wake in the hydrodynamic medium accomplishes this, but at least in our treatment it does not do enough. It could be that we are underestimating the amount of energy deposited in the medium at these angles, or it could be that we are missing energy at these angles that corresponds to energy lost by the jet that hydrodynamizes only partially or not at all, or it could be that our treatment of background uctuations, background subtraction, and jet reconstruction is subtracting away more of the particles originating from the backreaction of the medium than happens in analyses of experimental data. And, as we will elaborate in section 5, the jet shape analysis is sensitive to semi-hard particles in the region of pT > 2 GeV, in which the small-perturbation assumption behind eq. (4.7) starts to break down. We shall continue this discussion in section 5.
Next, we now turn from the angular jet shape variable to another intrinsically hadronic observable that assesses the longitudinal structure of jets: jet fragmentation functions. These show the distribution of the z of the tracks in a jet, where z ptrack cos =pjet is the ratio of the longitudinal momentum of a single charged hadron in the jet (a single track) to the momentum of the whole jet. ( is the angle between the track and the jet axis.) Unlike in our previous publications [148, 149], where we only analyzed partonic fragmentation functions, the inclusion of hadronization allows us to do a direct comparison with experimental measurements. However, we must keep in mind that the hadronization process is not under good theoretical control and that we have only implemented a simpli- ed prescription which ignores changes in the color ow of jets that shower in a medium. We defer the study of e ects of di erent prescriptions for hadronization in vacuum and in medium within our model to the future.
Results from our hybrid model calculations of jet fragmentation functions are shown in gure 11 for LHC heavy ion collisions with 0-10% and 10-30% centrality. Both panels display the ratio between hadronic fragmentation functions in PbPb collisions to those in proton-proton collisions compared to experimental measurements of this quantity by CMS [51]. We include two bands, one for the full calculation including the e ects of the backreaction of the medium, and the other without it. The overall background subtraction and hadronization are the same in both calculations, as is the jet energy resolution correction and the subtraction of background tracks in the jet cone by -re ection.
{ 35 {
JHEP03(2017)135
2
2
1.8
Backreaction No Backreaction
CMS Data
0 10%
1.8
1.6
1.6
10 30%
100 < P jetT < 300 GeV 0.3 < | | < 2, r < 0.3
P hadron
T > 1 GeV
1.4
1.4
PbPb/pp
PbPb/pp
1.2
1.2
1
1
0.8
0.8
0.6
0.6
0.4
0.4
0 1 2 3 4 5
0 1 2 3 4 5
ln(1/z)
ln(1/z)
Figure 11. Ratio of the jet fragmentation function in PbPb collisions at ps = 2:76 ATeV to that in proton-proton collisions. As in gure 10, we compare the results of our hybrid model calculation with and without the inclusion of the particles coming from the backreaction of the medium in 0-10% centrality (left) and 10-30% centrality (right) collisions to data measured by CMS [51].
By comparing the results of our simulations with and without backreaction, the jet fragmentation functions clearly show where the particles resulting from the wake in the plasma that get reconstructed as part of the jet end up. The hard part of the jet is practically una ected by the backreaction of the medium, with an almost identical distribution of hard fragments in the two calculations. (The small di erences arise from the the small change in reconstructed jet momentum associated with the addition of soft particles from the backreaction of the medium to the jet.) Both of the simulations show an enhancement of hard fragments at the largest values of z (smallest log 1=z) in PbPb collisions. As we have seen in other ways, wide jets (with more softer fragments) lose more energy, so at any given energy the jets that remain tend to be narrower, and tend to contain fewer, and therefore more energetic fragments, than in proton-proton collisions. Such an enhancement therefore seems generic to any mechanism of energy loss which signi cantly reduces the soft, large angle, components of jets. (See ref. [105] for a similar e ect in a perturbative-based jet quenching Monte Carlo [76, 77, 81, 82].) The small-z region of the fragmentation function is sensitive to backreaction e ects. The emission of soft particles by the jet-induced wake compensates the suppression of soft fragments due to energy loss and leads to an overall enhancement of soft tracks in the PbPb jets relative to proton-proton jets. The comparison between our calculations with and without the particles coming from the backreaction of the medium also shows the range of momenta at which back reaction contributes signi cantly to the particles reconstructed as part of a jet, namely pT up to pT 2:5 GeV. At this scale our approximate approach to the wake distribution eq. (4.7)
underestimates particle production, as we discussed after eq. (4.5). At softer momenta, by neglecting the e ects of viscosity as the sound waves produced by the jet damp out and heat the plasma we are underestimating the particle production also.
Comparing the fragmentation function ratios that we have obtained in our calculations including the particles from the backreaction of the medium to the ratios measured by CMS [51], we see qualitative similarities, but not quantitative agreement. At large z (small ln(1=z)) our calculated ratio is above 1 while there is no evidence for this in the data.
{ 36 {
JHEP03(2017)135
However, in this regime the uncertainties in both data and theory are signi cant. At small z (large ln(1=z)) we see that the particles from the wake in the plasma that are reconstructed as part of the jet turn the fragmentation function ratio upwards, but they do not turn it as far upward as in the data, and they turn it upward at a larger value of ln(1=z) than in the data. This suggests that our treatment of the back reaction of the medium is missing an increase in the production of few GeV particles. We have already suggested several possible interpretations of this above. It is quite reasonable to expect that whatever physics it is that our approximations have missed that would push the large ln(1=z) fragmentation function ratios up would also push the jet shapes up at larger values of r in gure 10. We are missing some soft to few GeV particles at these values of r.
Certainly a more complete description of the backreaction dynamics, from the wake in the plasma to the resulting particles reconstructed as part of a jet, would be helpful. It is also worth mentioning, though, that we have found a disagreement between the fragmentation functions from Pythia simulations of jets in proton-proton collisions and the fragmentation functions measured in small R jets in these collisions, for example with the fragmentation functions from Pythia high by as much as 30% at small z, at least with the Pythia tunes we have explored (Monash, 4C and 4C with a modi ed Q0 parameter).
A similar disagreement has also been found between state-of-the-art parametrizations of fragmentation functions and proton-proton data [186]. These uncertainties are comparable in size to the size of medium modi cations themselves. This means that understanding fragmentation function ratios more quantitatively will require a better understanding of hadronization of jets, both those in vacuum and those that have been modi ed by passage through a medium. Regardless, it is already an inescapable conclusion that backreaction dynamics has a signi cant impact on the jet fragmentation function at small z.
In the next section we will explore a suite of observables that are even more sensitive to the soft particles from the wake in the plasma that end up reconstructed as a part of the jet than the jet shape or the jet fragmentation function.
4.4 Recovering the \lost" jet energy and \missing" pT deposited by a dijet pair
We have seen in the previous subsection that the backreaction of the medium in response to the passage of a jet, namely the wake in the plasma, contributes to single jet observables like the jet RAA, jet shape, and jet fragmentation function because some of the particles that result from the hadronization of the wake in the plasma must of necessity end up counted as part of the jet after the jet is reconstructed and the background is subtracted. In these observables, the e ects of the particles coming from hadronizing the wake in the plasma constitute small corrections overall, although they can play a signi cant role in certain kinematic regimes (like in the soft region of the fragmentation function ratio and the large r region of the jet shape ratio). In this subsection, we focus on a suite of observables which are dominated by e ects originating from the backreaction of the medium in response to the passage of two jets, a dijet, through it.
The principles behind our simpli ed calculation of the spectrum of particles produced as the medium (including the wake therein) hadronizes are that the medium acquires
{ 37 {
JHEP03(2017)135
the energy and momentum lost by the jets passing through it, which thermalizes subject to energy and momentum conservation. This implies collective motion of the medium, a wake in the quark-gluon plasma. After hadronization, this energy and momentum is recovered in the form of soft particles with a wide angular distribution given by eq. (4.7). This makes it natural for us to calculate the so-called \missing-pT observables" introduced recently by CMS [60], and compare our calculations to measurements made by CMS. These observables are de ned to extract the spectrum and angular distribution of the particles that correspond to the net momentum lost by the two jets in a dijet pair. Each jet in the pair loses energy and momentum to the plasma, but in general one will lose more than the other. By constructing observables designed to recover the net momentum lost by the pair of jets we obtain observables that are dominated by the response of the medium to the jets. These observables are constructed from all the tracks in an event, not just those reconstructed as part of a jet. We shall set K = 0 throughout this subsection, neglecting transverse momentum broadening. We have obtained results with K [negationslash]= 0, but
the K-dependent e ects that we see in these missing-pT observables are all much smaller in magnitude than the e ects of particles coming from the hadronization of the wake in the plasma. As these observable consequences of the backreaction of the medium are our focus in this subsection, there is no reason to introduce K [negationslash]= 0.
The observables that we look at in this subsection are intrinsically hadronic and, therefore, are hard to bring under full theoretical control without signi cant modelling of the di erences between hadronization in heavy ion collisions and proton-proton collisions. Nevertheless, we will observe several qualitative features of our results which are similar to what is seen in the experimental data. This indicates that our treatment, although simpli ed, captures some of the main aspects of the collective reaction of the plasma to the deposited jet energy.
Following ref. [60], the missing pT analysis consists in studying the conservation of momentum in heavy ion events which include hard jets reconstructed with a speci ed anti-kt reconstruction parameter R. The rst step in the analysis is to use the anti-kt algorithm to nd a sample of events containing at least two jets and to determine the pT and the jet axis of the leading and subleading jet in each event. We then de ne the dijet angle dijet
as the bisection between leading and subleading, where (sub)leading is the azimuthal angle of the (sub)leading jet. The only use of the jet reconstruction algorithm is the selection of the event sample and the determination of the jet pTs and axes and hence dijet. In the remainder of the analysis, all charged tracks in the events in the sample are used.
The next step in the analysis is to project the momenta of each of the charged tracks in an event along the dijet axis. This projection is de ned by
=
p[bardbl]T pT cos (dijet ) ; (4.11)
where pT and are the transverse momentum and azimuthal angle of the track. With this convention, =
p[bardbl]T is positive for tracks in the subleading jet hemisphere and negative for tracks in the leading jet hemisphere. By momentum conservation, the sum of all the =
p[bardbl]T for
all the tracks in an event must be zero. However, in the experimental analysis only charged
{ 38 {
JHEP03(2017)135
p[bardbl]T need not cancel identically.
Nevertheless, by studying the approximate cancellation of this momentum component as a function of which tracks we include in the analysis, out to what angular distance from the dijet axis, we can extract valuable information about the distribution of the lost jet energy.
We start by computing the =
p[bardbl]T distribution sliced in bins, where is the distance in the ( ; ) plane between the track in question and either the leading jet axis or the subleading jet axis, whichever the track is closest to. We consider dijet pairs with leading and subleading transverse momenta pTleading > 120 GeV and pTsubleading > 50 GeV respec
tively, and with both jets within [notdef] [notdef] < 2. We also limit our sample to events in which
the two jets are back-to-back, our criterion being > 5=6. Finally, after identifying the dijet pair, only those dijets in which both jets have [notdef] [notdef] < 0:6 are considered. (We
start by nding jets within a much larger [notdef] [notdef] window in order to be as con dent as we
can be that when we have identi ed a dijet in which both jets have [notdef] [notdef] < 0:6 these two
jets really are the leading and subleading jets in the event.) All these speci cations of our event sample are the same as those in the experimental analysis of ref. [60]. The result of this analysis is shown in gure 12 for our full calculations including the medium response (left) and, for comparison, without medium backreaction (right). In the right panel, we con rm that momentum is not conserved in our hybrid model by its construction, when we do not include the particles coming from the wake in the plasma. The left-most bin in shows the momentum imbalance between the high-momentum tracks in the two dijets themselves, with more momentum in the leading jet direction and less in the subleading jet direction. In our hybrid model the two jets have lost di erent amounts of momentum, and in the right panel we have neglected the fact that the lost momentum is deposited as two wakes in the medium. The left panel of gure 12 clearly shows how the inclusion of the response of the medium, modeled via eq. (4.7), results in a population of soft particles spread broadly in with a net momentum in the subleading jet hemisphere (positive =
p[bardbl]T).
The subleading jet has lost more momentum, and hence the wake in the plasma moving in the subleading jet direction contains more momentum than the wake moving in the leading jet direction.
Another way to see the e ect of including the backreaction of the medium is to compare the solid black curves in the two panels of gure 12. While the accumulated =
p[bardbl]T tends to
zero as increases in the calculations in the left panel, in the right panel the accumulated =
p[bardbl]T is 20 GeV after all bins have been summed. This value, which corresponds to the
average imbalance in the energy lost by the two jets in the dijet pair, compares well with simple estimates [85]. The medium response transforms the \lost" energy from the two jets into two wakes, and hence into softer particles in the range of 0:5 2 GeV whose net
momentum is preferentially in the hemisphere of the more quenched, subleading, jet. This is shown by the enhancement of soft =
p[bardbl]T tracks observed in the left panel of the gure. The dashed curve, which is the same in both panels, shows the accumulated =
p[bardbl]T for
proton-proton collision events as described by our Pythia tune. Recalling that our analysis employs all the tracks in an event, not just those reconstructed as part of a jet, it is perhaps not too surprising that the =
p[bardbl]T distribution generated by Pythia that is described by the
{ 39 {
tracks with [notdef] [notdef] < 3 are considered, and therefore, the net =
JHEP03(2017)135
10
0
10
20
30
40
Quenching + Medium Response
1.0-2.0
2.0-4.0
4.0-8.0
8.0-300.0
Figure 12. Results from our calculations of the missing =
p[bardbl]T for tracks that are away from the dijet axis in ( ; )-space, in PbPb collisions containing a pair of R = 0:3 jets satisfying kinematic criteria described in the text. The rst nine bins in are each 0.2 wide, together covering 0 < < 1:8. The tenth bin, plotted at = 1:9, contains the missing =
p[bardbl]T for all the tracks with > 1:8. Note that we limit our analysis to tracks with [notdef] [notdef] < 3. In the left panel, we show the results of our full
calculation, with energy loss calculated within our hybrid model with no momentum broadening, and with particles coming from the backreaction of the medium included. The right panel shows our results if we leave out the particles coming from the backreaction of the medium. In both plots, the solid black points show the net =
p[bardbl]T of all charged tracks with [notdef] [notdef] < 3 in each bin. The
contribution to this net momentum of tracks in di erent momentum bins are codi ed by the colored histograms. In both plots, the solid line shows the cumulative sum of all the =
p[bardbl]T in a bin and all the bins to its left. The black points, the colored histograms, and the solid line are all calculated for 0-30% central PbPb collisions with ps = 2:76 ATeV. For comparison, we have repeated the calculation for proton-proton collisions as described by our Pythia tune, and the dashed line in each plot shows the same cumulative sum of =
p[bardbl]T for proton-proton collisions that the solid line shows
for PbPb collisions.
dashed curve is not in full agreement with measurements made in proton-proton collisions reported in ref. [60]. This small discrepancy illustrates the di culty in obtaining full theoretical control over this inherently hadronic observable. What we shall do in subsequent gures is to subtract the =
p[bardbl]T distribution in our Pythia calculation of proton-proton calculations from the =
p[bardbl]T distribution in our PbPb calculation including black-reaction. That is, we shall subtract the dashed curve from the solid curve in the left panel of gure 12. In this way, we shall focus explicitly on e ects that are due to medium-induced energy loss, and the wake in the medium, both of which our model is intended to describe.
By comparing the solid and dashed curves in the left panel of gure 12, we see that the net momentum lost by the pair of jets is distributed over a wide angular region, over [lessorsimilar] 1. We can also see from the similarity between the red histograms in the two panels in the gure that the backreaction of the medium does not a ect the hard tracks in the event, namely the hard components of the jet. We also see that almost the entire imbalance in the
{ 40 {
10
0
Quenching Only
p[bardbl] T [angbracketright] [GeV]
[angbracketleft]/
p[bardbl] T [angbracketright] [GeV]
10
20
[angbracketleft]/
30
JHEP03(2017)135
40
0 0.5 1 1.5
0 0.5 1 1.5
high-momentum tracks is an imbalance at small with the net momentum in the leading jet hemisphere, whereas almost the entire imbalance in soft tracks is in the hemisphere of the event in which that corresponds to particles from the (larger) wake of the subleading jet. The experimentally measured distributions [60] are quite similar to what we see in the left panel of gure 12, but because of the discrepancy between our calculations and data for proton-proton collisions we will postpone a quantitative comparison until we have results for the di erence between the missing momentum in PbPb and proton-proton collisions, below.
Following ref. [60], we sharpen the analysis by studying how the =
p[bardbl]T distribution depends on the dijet asymmetry AJ (pTleading pTsubleading)=(pTleading + pTsubleading). In events in which the dijet pair features a larger AJ we should expect a correspondingly larger momentum asymmetry in the soft particles coming from the hadronization of the wakes of the two jets. As discussed above and as in ref. [60], we shall also focus on the di erence between the =
p[bardbl]T distribution in PbPb collisions and that in proton-proton collisions containing dijets satisfying the same kinematic criteria. Jet quenching makes the
AJ-distribution wider in PbPb collisions, but it is already rather wide in proton-proton collisions because dijets are often produced in association with a third jet or other radiation, making it important for us to focus on the di erence between PbPb and proton-proton collisions.
In the upper and middle panels of gure 13 we show the value of =
p[bardbl]T obtained from our calculations upon integrating over all angular separations in PbPb collisions minus the same obtained in proton-proton collisions, in both cases for dijet events reconstructed using the anti-kt algorithm with R = 0:3 in bins of the dijet asymmetry AJ, in both cases for collisions with two di erent centralities. In the upper panels, we include the modi cations to the jets because of energy loss but leave out the contributions from the backreaction of the medium. The middle panels show the results from our full calculation, including both the modi cations to the jets and the wakes in the plasma. The black points in the upper panels show the momentum imbalance in dijet events when the momentum carried by the wakes are neglected. The black points in the middle panels con rm that total momentum is conserved regardless of the asymmetry AJ | once we include the particles coming from the wakes in the plasma. We see from the red, green and orange histograms in the upper panels of gure 13, though, that events with a larger and larger dijet asymmetry AJ feature more and more missing =
p[bardbl]T, in particular for tracks with pT > 2 GeV, with a sign corresponding to there being more tracks with these (hard) momenta in the leading jet direction. That is, jet quenching | here described in our hybrid model | suppresses the hard particle contribution in the subleading jet hemisphere, as the subleading jet has lost more energy. The energy and momenta lost by both jets have been thermalized as wakes, which in turn are transformed into softer hadrons when the medium including the wakes freezes out. We see the contributions of the momentum carried by hadrons coming from the wakes in the plasma by comparing the upper and middle panels of gure 13. In the blue and yellow histograms in the upper panels of the gure, we see the missing =
p[bardbl]T
for tracks with pT < 2 GeV with a sign corresponding to there being more tracks with these (soft) momenta in the subleading jet direction. Since the subleading jet lost more momentum, its wake contains more momentum and this shows up in the missing =
p[bardbl]T of the
{ 41 {
JHEP03(2017)135
PbPb-pp R=0.3, 0-10%
No Backreaction
40
h/p[bardbl]T [angbracketright]
0.5-1.0
1.0-2.0
2.0-4.0
4.0-8.0
8.0-300.0
PbPb-pp R=0.3, 10-30%
20
p[bardbl] T [angbracketright] [GeV]
0
[angbracketleft]/
-20
-40
JHEP03(2017)135
0 0.1 0.2 0.3 0.4
0 0.1 0.2 0.3 0.4
AJ AJ
40
PbPb-pp R=0.3, 0-10%
PbPb-pp R=0.3, 10-30%
h/p[bardbl]T i
0.5-1.0
1.0-2.0
2.0-4.0
4.0-8.0
8.0-300.0
20
p[bardbl] T i [GeV]
0
h/
20
40 0 0.1 0.2 0.3 0.4
0 0.1 0.2 0.3 0.4
AJ AJ
40
40
PbPb - pp R = 0.3, 0-10%
20
0
-20
-40
PbPb - pp R = 0.3, 10-30%
0.5-1.01.0-2.02.0-4.04.0-8.08.0-300.0
p
|| T
D
20
[GeV]
0
||
T
p
-20
-40
0.1 0.2 0.3 0.4
0.1 0.2 0.3 0.4
A
A
J
J
Figure 13. The middle two panels show the results of our full calculation, while the upper two panels show the results from our hybrid model for energy loss with the contribution from the backreaction of the medium left out. Black points with error bars show the di erence between the missing =
p[bardbl]T in PbPb and pp collisions, integrated over all , summed over all track momenta, binned in the dijet asymmetry AJ, for events containing R = 0:3 dijets in collisions with 0-10% (left) and 10-30% (right) centrality. The colored histograms show the missing =
p[bardbl]T for tracks whose momenta is in the indicated range. The lower two panels, obtained from ref. [60], show the results from the same analysis performed on experimental data by the CMS collaboration.
{ 42 {
soft tracks in the event. All these e ects are larger in the central collisions on the left than in the more peripheral collisions on the right, as expected since due to the smaller size of the reball the dijets are on average less quenched in more peripheral collisions.
Almost all the features seen in the middle panels of gure 13 that we have described above are present in the missing =
p[bardbl]T measurements reported by the CMS collaboration [60] that we show as the lower panels in gure 13, at the least at a qualitative level and in most cases to a degree that constitutes agreement within the uncertainties.
The biggest di erence between the results of our calculations and the experimental results is found for semi-hard particles in the 2 GeV< pT < 4 GeV momentum range. By comparing the orange histograms in the upper and middle panels of gure 13, we see that in our calculation the wake contributes little in this momentum range: the \orange contribution" to the missing =
p[bardbl]T imbalance due to jet quenching that is seen in the upper panels remains almost unmodi ed in the middle panels. Instead, in the data there is almost no missing =
p[bardbl]T imbalance in this momentum range. Furthermore, by comparing the upper and middle panels we see that in our calculation the contribution to the missing =
p[bardbl]T
from soft particles with pT < 2 GeV is dominated by the e ects of the backreaction of the medium, and then by comparing the middle and lower panels we see that our calculations yield a greater contribution from these soft particles than in the data. So, although our implementation of medium backreaction does restore momentum conservation | the =
p[bardbl]T
excess in hard particles in the direction of the leading jet due to jet quenching is balanced by the =
p[bardbl]T excess in soft particles in the direction of the subleading jet coming from the wake in the plasma | our calculation of the particles coming from the wake seems to yield somewhat more particles with pT < 2 GeV than in the data and substantially fewer particles with 2 GeV < pT < 4 GeV. We will defer further discussion of the origin of this discrepancy to section 5. For the present it su ces to remember from our discussion earlier in this section that several of the approximations that we used in setting up our crude treatment of the particles from the wake break down for particles with pT substantially greater than the temperature of the medium at freezeout, and they break down in a way such that we anticipated underpredicting the production of semi-hard particles from the wake.
Finally, we turn to the dependence of the missing =
p[bardbl]T distribution on the value of the parameter R used in nding and reconstructing the jets, recalling that using a large R results in nding a sample of jets that have a larger average angular size. In gure 14 we show the =
p[bardbl]T distribution as a function of the angular separation for di erent values of the jet reconstruction parameter. For the purposes of focusing on the e ects of jet quenching and in order to make a comparison with experimental data, we show the di erence between PbPb and proton-proton collisions. As in the previous plots, the colored histograms show the contributions of tracks in di erent pT-ranges, and the black points show =
p[bardbl]T for
tracks with all values of pT in a given -bin. (As always when we show results from our calculations, the error bars include the uncertainties coming from the error bar on the experimental data point that we used to x the sc parameter in our hybrid model and from our estimate of the systematic uncertainty in our model that we make by varying the temperature Tc below which we turn o parton energy loss.) As in gure 12, in the middle
{ 43 {
JHEP03(2017)135
4
0
p[bardbl] T i [GeV]
4
h/
1.0-2.0
2.0-4.0
4.0-8.0
8.0-300.0
8
PbPb-pp, 0-30%
No Backreaction
12 0 0.5 1 1.5
0 0.5 1 1.5
0 0.5 1 1.5 0 0.5 1 1.5
4
JHEP03(2017)135
0
R=0.2
PbPb-pp, 0-30%
R=0.3
R=0.4
R=0.5
h/p[bardbl]T i
0.5-1.0
1.0-2.0
2.0-4.0
4.0-8.0
8.0-300.0
p[bardbl] T i [GeV]
4
h/
8
12 0 0.5 1 1.5
0 0.5 1 1.5
0 0.5 1 1.5
0 0.5 1 1.5
4
4
-8
PbPb-pp, 0-30%
-12
4
-8
-12
4
-8
-12
0
[GeV]
-4
T
||p
-8
4.0-8.08.0-300.0
-12
0.5 1 1.5
0.5 1 1.5
0.5 1 1.5
0.5 1 1.5
D
D
D
D
Figure 14. Evolution in R of the missing =
p[bardbl]T observable, AJ inclusive and sliced in , for the di erence between PbPb and pp collisions. Top four panels are the results of our calculations with no contributions from the backreaction of the medium. Middle four panels are the results of our full calculations, in which the =
p[bardbl]T imbalance in hard particles due to jet quenching is balanced by that in soft particles coming from the wakes of the two jets. The R = 0:3 panels are the di erence between the PbPb results shown in gure 12 and the same for proton-proton collisions. Bottom four panels come from the analysis of experimental data by the CMS collaboration reported in ref. [60]. There is some hint of R-dependence, despite the big error bars, especially in the rst bin.
panels of gure 14 momentum conservation is (approximately) recovered after summing over all bins, although not precisely because we only use tracks with pT > 0:5 GeV and
| [notdef] < 3 in the analysis.
As already mentioned, in this analysis the anti-kt reconstruction algorithm is only employed to determine the sample of events over which the =
p[bardbl]T distribution is computed, as well as the jet and dijet axes in these events. All tracks in the events are then used in the analysis, whether or not the reconstruction algorithm counted them as being within a jet. By increasing the anti-kt parameter R while keeping the momentum cuts on the dijet
{ 44 {
selection xed, we include jets that are on average wider in angle and that include more, and hence softer, fragments. Since wider jets containing a larger number of fragments lose more energy, the missing =
p[bardbl]T in the hardest pT tracks in the smallest bin grows in magnitude with increasing R. Despite the large error bars in our calculations and in data, this enhancement is present in both, as seen in gure 14. This is further evidence for the fact that wider jets lose more energy than narrower jets do.
Beyond the rst bin, we see in gure 14 that the distribution of =
p[bardbl]T is dominated
by softer particles with net =
p[bardbl]T in the subleading jet direction. The -dependence of the softest components of the =
p[bardbl]T distribution shows very little sensitivity to R. This may be understood as a consequence of the wide angular region in which energy is recovered, controlled by eq. (4.7). If we look carefully at the =
p[bardbl]T distribution for the hardest particles in the 0:2 < < 0:8 range of angles, we do see an R-dependence in the results from our calculations and in data: the =
p[bardbl]T imbalance in this kinematic regime decreases with increasing R. In our calculations and in the data, this imbalance comes from the proton-proton collisions: there is no signi cant imbalance in the hard particles at these angles in the PbPb collisions for any R. This indicates that the e ect is a sign of the presence of third jets at these angles. When R is small, a structure that would have been reconstructed as a single jet if R were larger can be reconstructed instead as two nearby jets, which end up being counted as the subleading jet in a dijet and as a third jet. The third jet results in an imbalance in the =
p[bardbl]T distribution for hard particles away from the dijet axis in the proton-proton collisions, but this weak third jet is greatly quenched in the PbPb collisions. The result in the di erence between PbPb and proton-proton collisions is as seen in gure 14, and the e ect is reduced when jets are reconstructed with larger R.
Remarkably, for all the features in gure 14 that we have described to this point there is good qualitative agreement between the results of our calculations in the middle panels and the experimental distributions reported in [60] that we show in the lower panels. We do wish to note two qualitative disagreements, however. First, our results for the missing =
p[bardbl]T distribution di er from the experimental results in the last bin in which, we recall, corresponds to the =
p[bardbl]T imbalance for particles with all > 1:8. We have checked that this disagreement is entirely due to a di erence between our results and experimental results for proton-proton collisions; we do not know how to interpret this particular di erence between Pythia and experiment. Second, as in gure 13 we see that our calculation of the =
p[bardbl]T distribution of the particles coming from the backreaction of the plasma to the jets, (which we can discern by comparing the middle panels to the top panels) is lacking in \orange particles" with 2 GeV < pT < 4 GeV relative to what we see in the data in the lower panels. Correspondingly, our calculation restores momentum conservation with a greater imbalance in the softest particles with pT < 2 GeV than in the data. Our calculation of the particles coming from the wakes in the plasma has signi cantly fewer semi-hard particles and somewhat more soft particles. We will discuss this disagreement further in the next section.
Before we move on, it is important to look back at the way in which we de ned the new observable that we introduced at the end of section 3 that, as we argued in our discussion
{ 45 {
JHEP03(2017)135
around gure 6 and gure 7, is particularly sensitive to the e ects of transverse momentum broadening. In constructing the hadronic observable in gure 7, we only used particles with 5 GeV < pT < 10 GeV. We can now see clearly why we chose not to use particles with pT < 5 GeV in this observable. We want this observable to only be sensitive to (modi cations of) the jets themselves. We want it to be una ected by the wake in the plasma and, more generally, by the energy lost from the jet wherever that lost energy ends up. In our calculation, that lost energy dominates the missing =
p[bardbl]T distribution for pT < 2 GeV; in the data, it is clear that it is important for pT < 4 GeV. In order to avoid all the physics that is of interest to us here in section 4, in gure 7 in section 3 we ignored particles with pT < 5 GeV.
5 Discussion and outlook
In this paper we have described the origin and observable consequences of two important e ects upon implementing and describing them in our hybrid model of jet quenching: transverse momentum broadening, as the partons in the jet are kicked transverse to their direction of motion during their passage through the medium, and the backreaction of the medium to the momentum and energy deposited in it by the passage of the jet through it.
Our discussion of transverse momentum broadening is based on a Gaussian distribution of kicks in transverse momentum experienced by each parton in the jet. This approximation captures the bulk of the distribution of transverse momentum-transfer but neglects rare scatterings by larger angles that impart more transverse momentum. Quite surprisingly, we have found that adding Gaussian momentum broadening has little impact on typical jet observables such as the jet spectrum, jet shapes and dijet acoplanarity, even when we choose unrealistically large values of K ^
q=T 3, the new second parameter that we add to our model to quantify the strength of the transverse momentum broadening. The reason why it is hard to see e ects of broadening in these observables can be traced back to the parton energy loss and consequent jet quenching, described successfully by our hybrid model. We nd that in our hybrid model parton energy loss results in jets with a given energy observed in PbPb collisions being narrower than jets with the same energy in proton-proton collisions. This happens because jets that are wider in opening angle lose more energy than those that are narrower and because the jet energy spectrum is steeply falling, and has also been seen recently in entirely weakly coupled analyses [105] and in entirely strongly coupled analyses [147]. This narrowing of jets with a given energy because of energy loss turns out to be a larger e ect than the broadening of jets due to kicks to the transverse momentum of partons in the jets. Also, for narrower jets the principal e ect of transverse kicks is a change in the overall direction of the jet, rather than changes to its shape. (As an extreme case, think of a jet containing only one parton.) For a high energy jet the resulting acoplanarity is small, but for narrow jets the change to their shape is even smaller. (In the extreme case of a single parton, the shape of the resulting jet does not change at all if the direction of that parton changes.) Although initially surprising, once understood these considerations make it clear that it will be exceedingly di cult to extract the value of K, or the jet quenching parameter ^
q, directly from the e ect of transverse momentum
{ 46 {
JHEP03(2017)135
broadening on traditional jet observables. (In weakly coupled analyses but not in general, ^
q is related to parton energy loss and so can be constrained indirectly via observables that are sensitive to energy loss.)
We have, however, found a new kind of jet shape ratio that is partially di erential in the pT of the tracks entering its de nition that does exhibit sensitivity to the value of K. In gure 6 we restrict the pT-range of the partons entering the analysis to be 10 < pTparton < 20 GeV, as among the options we have investigated it is for partons in this momentum range that we have found the most sensitivity to the value of K. We then need a hadronic observable that is sensitive to the angular distribution of partons in this momentum range, and in gure 7 we construct a partially di erential jet shape ratio observable using hadrons in the pT-range 5 < pThadron < 10 GeV. The results we obtain indicate that experimental measurements of this observable, or similar observables, could be used to constrain the value of the broadening parameter K, a key property of the medium. Restricting attention to pT above a minimum value (like the 5 GeV we have chosen) is important because it makes the observable insensitive to the other physical e ect that we describe in this paper, namely the reaction of the medium to the passage of the jet. More on that below. Restricting attention to pT below a maximum value (like the10 GeV we have chosen) ensures that we are looking at jets containing many tracks and are looking at tracks that are not too energetic and thus can be a ected visibly by transverse momentum broadening.
Although the successes of our model to date [148, 149] motivate the investments in improving it that we describe in this paper, in the long run the most instructive use of a simple few-parameter model like our hybrid model is to discover and understand the ways in which it fails to describe data. Let us look ahead by imagining one particularly interesting way in which this could happen in future. Let us imagine a day when the measurement of observables like the partially di erential jet shape ratio subleadingpT(r) plotted in gure 7 have been used to discover experimental evidence for medium-induced transverse momentum broadening and to constrain the value of K. What will the next step be? At such a point in time, the key next question will be to ask whether the event-by-event and parton-by-parton distribution of the transverse kicks contributing to an ensemble-averaged observable like subleadingpT(r) is indeed Gaussian. This distribution is Gaussian by hypothesis in our hybrid model, meaning that our model will ultimately serve as a baseline against which to look for evidence of a non-Gaussian distribution. If the quark-gluon plasma in QCD were strongly coupled on all length scales, including short length scales, the distribution of transverse kicks would be Gaussian for all values of the momentum transfer, including large momentum transfers. This is not what is expected in an asymptotically free theory like QCD. If the strongly coupled liquid QCD plasma can be probed at high enough momentum transfer, the weakly coupled quarks and gluons present within it at short length scales should manifest themselves. Once the overall strength of the broadening K has been constrained, the way to look for point-like scatterers within the strongly coupled plasma will be to look for rare events in which a parton within a jet is scattered by an angle that would be improbably large if the distribution of transverse momentum kicks were a Gaussian with overall strength speci ed by K [151, 170]. The hallmark of
{ 47 {
JHEP03(2017)135
point-like scatterers is a distribution of momentum transfer that has a power-law tail at large momentum transfers. If in future our hybrid model, with its Gaussian distribution of transverse kicks with strength K, fails to describe data because it neglects rare (but only power-law-rare) large angle scatterings, their discovery may potentially be used to determine at which small length scale the uctuations within strongly coupled quark-gluon plasma behave as weakly coupled quark and gluon quasiparticles.
We nd that standard intrajet observables, like the conventional jet shape and jet fragmentation function, are sensitive to the presence of the wake in the plasma left behind by the passing jet. Since the wake carries the momentum lost by the jet, the particles in the nal state that come from the hadronization of the medium including the wake must have net momentum in the direction of the jet, meaning that some component of what experimentalists reconstruct as a jet, even after background subtraction, must in fact be composed of soft particles coming from the wake in the plasma. Motivated by expectations coming from explicit calculations of the wake in strongly coupled plasma [123, 125] and by recent weakly coupled calculations in which the energy lost by a jet cascades into multiple soft particles [150{156] we analyze the wake upon making the simplifying assumptions that all the energy and momentum lost by the jet is incorporated into hydrodynamic motion and that the resulting perturbations to the temperature and velocity of the hydrodynamic uid are small. That is, we assume that all of the lost momentum and energy go into a hydrodynamic wake that is as thermalized as it can be by the time of freezeout, yielding only small corrections to the spectra of the hadrons in the nal state, a nal state that remembers nothing about the origin of this perturbation other than the momentum and energy dumped into the plasma. (In a weakly coupled calculation as in [156] in which a part of the lost energy hydrodynamizes, these assumptions are relevant for this component of the lost energy, denoted Ethermal in [156].) In our calculation of the corrections to the hadronic
spectra coming from the wake in the plasma, we have made further simplifying assumptions (for example assuming boost invariant longitudinal expansion and neglecting transverse expansion) such that we have been able to determine the corrections to the hadron spectra in the nal state coming from the hadronization of the wake in the plasma analytically, and without needing to incorporate any additional parameters into our model. By comparing our computations with and without the e ects of the backreaction of the medium, we have established that, indeed, these collective e ects have important consequences for jet shapes, fragmentation function and missing-pT observables as discussed in sections 4.3 and section 4.4.
Including the e ects of the backreaction of the medium results in enhanced jet shapes at large angles and enhanced fragmentation functions at soft momenta. These e ects go in the directions favored by experimental data, but in our treatment they are not as large in magnitude as the data seem to require. It should also be noted that these observables are not dominated by these e ects.
The distribution of the missing-pT observable =
p[bardbl]T constructed from dijet events is particularly interesting because at all but the smallest angles it is dominated by e ects coming directly from the backreaction of the medium. Even with our simpli ed approach to the backreaction of the medium and to hadronization, we obtain a rather good qualitative de-
{ 48 {
JHEP03(2017)135
scription of the shape, pT-dependence, angular dependence, and AJ- and R-dependence of
the =
p[bardbl]T distributions, see gure 13 and gure 14. This indicates that our simpli ed analysis of the backreaction of the medium is on the right track, and con rms the absolute necessity of including this physics if one aims to understand the soft and/or large-angle components of jets as reconstructed in heavy ion collisions.
The main discrepancy that we nd between the results of our calculations and the experimentally measured =
p[bardbl]T distributions is that our calculation does not yield enough semi-hard particles coming from the wake in the plasma in the 2 < pT < 4 GeV regime at angles less than 0:5 radians from the jet axis. Given that it has momentum conservation
built into it, our calculation must overproduce particles in some other regime, and indeed we nd that our wake results in somewhat more particles with pT < 2 GeV than indicated in the data, over an even broader range of angles. Although we cannot be sure of this, we expect that the remaining discrepancies between our results and experimental data on jet shapes and fragmentation functions can also be traced to our calculation not producing enough particles with 2 < pT < 4 GeV.
We are aware of three possible reasons for the discrepancy that we have found in the semi-hard momentum regime, for momenta well above T but not so high that particles coming from the wake can be neglected. First, we know that our perturbative expression (4.5) for the correction to the spectrum of particles coming from the plasma due to the wake is not valid for pT well above T , and furthermore we know that because it is based on a linear approximation to an exponential the expression (4.5) necessarily underestimates particle production from the wake in this regime. And, even if the modi cations of the ow eld induced by the wake are small their e ect on the spectrum of particles produced at freezeout can be signi cant on the tail of the thermal distribution. Going beyond our simpli ed expressions (4.5) and (4.7) in order to improve upon our description of the semi-hard momentum regime requires taking account of the explicit space-time dependence of the hydrodynamic wake in the plasma [174]. It would be interesting to investigate the effect of the full ow pro le as predicted by holographic computations of the hydrodynamic wake [123, 125] on the nal particle production spectrum to see how this modi es our results in the semi-hard momentum regime.
The second possible reason for the discrepancy that we have found in the semi-hard momentum regime could be that the energy lost by the jet does not fully hydrodynamize. That is, a part of the discrepancy at semi-hard momenta between our results for the =
p[bardbl]T
distribution and the experimental results in gure 13 and gure 14 could be a manifestation of non-equilbrium, and non-hydrodynamic, physics in quark-gluon plasma that has been disturbed by the passage of the jet. This possibility can be thought of in two ways that sound di erent but are operationally equivalent in the present context. One way of describing this possibility is to say that the energy and momentum deposited by the jet in the plasma has not all had time to thermalize as fully as we assume in our simpli ed analysis. This may be particularly relevant to energy deposited not long before the jet exits the droplet of plasma. Another way of describing the same possibility is to say that our hybrid model neglects medium-induced gluon radiation: in our model, the partons in the jet lose energy to a wake in the plasma and are kicked in transverse momentum but the
{ 49 {
JHEP03(2017)135
branching structure of the jet shower is not modi ed in the way that happens in entirely weakly coupled treatments of jet quenching. In a weakly coupled analysis in which energy is lost via medium-induced gluon radiation, a fraction of the energy lost in this way (called Ethermal in ref. [156]) is carried by radiated gluons that themselves further radiate and
cascade down to the medium scale and form a hydrodynamized wake in the plasma [156]. With regards to this component of the energy lost via medium-induced gluon radiation, neglecting the medium-modi cation of the branching structure of the jet shower as we do in our model makes absolutely no di erence: our model includes the loss of energy from the partons in the jet and includes the hydrodynamized wake in the plasma. Where these details will matter is for those radiated gluons that do not thermalize as fully as we assume, as may be particularly relevant to gluons radiated not long before the jet exits the droplet of plasma. The rst description uses words that arise in strongly coupled analyses of jet quenching. The second description uses words that arises in weakly coupled analyses. The two are operationally equivalent here.
Another physical process which has not been incorporated into our approach is the e ect of imperfect resolution. In our hybrid model as currently implemented, as soon as one parton in the jet shower splits into two, the two daughter partons are treated as separate losers of energy. That is, the medium can resolve the separation of one parton into two as soon as it happens, when the two daughters are still very close to each other. This cannot be correct. In reality, the medium can only resolve two partons as two separate losers of energy after they have separated by some resolution length that can reasonably be expected to be of order the Debye length of the plasma. The main e ect that the inclusion of the e ect of imperfect resolution would have in our model is that semi-hard partons within jets, in particular within narrow jets, would start losing their energy rapidly on their own later than in our present model, and so would end up less quenched. Hence, our neglect of these e ects of imperfect resolution is a third possible reason for the discrepancy that we have found in the semi-hard momentum regime.
It would be very interesting to assess how much of the discrepancy that we nd in the semi-hard momentum regime is hydrodynamic (and due to our simpli cations), or arises because some of the energy lost by the jet does not fully hydrodynamize, or is related to our neglect of resolution e ects. Each explanation is quite interesting in its own right and deserves investigation. We leave this to future work, but note here that this example serves very well to make our point that the most instructive uses of a few-parameter model like the one we have constructed are to discover, and subsequently understand, the ways in which it fails to describe data. This example also serves to highlight the full power of the suite of experimental data on intrajet structure now available, including jet shapes, fragmentation functions, and the missing-pT observables. We encourage theorists pursuing other approaches to calculate all these observables and compare to data. Consider for example the weakly coupled calculations in which energy is lost to radiated gluons which themselves are quickly degraded down to soft gluons with momenta of order the temperature [150{156], which leads to rapid hydrodynamization of the emitted energy [156]. These approaches yield qualitative agreement with the =
p[bardbl]T-distribution that we have shown in gure 12 [152], but fail to reproduce the soft enhancement of the fragmentation function
{ 50 {
JHEP03(2017)135
unless further soft radiation attributed to decoherence is added [187]. It would be very interesting to see how well they predict the AJ-dependence and the R-dependence of the =
p[bardbl]T-
distribution as in gure 13 and gure 14, in particular in the soft and semi-hard momentum regimes. Weakly coupled approaches in which it is assumed that the emitted gluons do not further interact at all are able to reproduce the enhancement of jet shapes at large angles [104]; it would be interesting to analyze the =
p[bardbl]T-distributions in these approaches.
Although the physics of the wake in the plasma is of interest in its own right, if in the longer term our goal is to quantify broadening, determine the value of the key medium property K, and ultimately to see rare but not-too-rare large angle scattering of partons within jets that may allow us to see the length scale at which strongly coupled quark-gluon plasma emerges from the uctuations at even smaller length scales that behave as weakly coupled quark and gluon quasiparticles, then we will want to focus on observables sensitive to the angular distribution of 10-20 GeV partons within jets as in gure 7, where the wake makes no contribution.
A Transverse kicks kinematics
As a parton propagates though the plasma, it receives a random transverse kick in the local uid rest frame that we shall denote by q. In this appendix, we provide a precise speci cation of what we mean by this. Denoting the momenta of the parton in the uid rest frame before and after a transverse kick by P = EF (1; wF ) and P [prime] = EF (1; w[prime]F ), and assuming that the transverse kick serves only to change the direction of wF , meaning w[prime]2F = w2F , the momentum after the kick is given by
w[prime]F =
s1 q2E2F w2FwF + qEF e? ; (A.1)
where e? is a vector perpendicular to the parton velocity in the rest frame of the uid.
Note that the virtuality of the parton has not changed. And, in the uid rest frame, neither has its energy.
Since our computation is performed in the collision center-of-mass frame, we need to express the acquired momentum in that frame. Let us de ne the four-vector W = P=E, with E the energy of the parton in the collision frame. Denoting the uid velocity in the collision center-of-mass frame by u = F (1; v), the energies in the uid and collision frames are related by EF = E F (1 wv), with w the velocity of the parton in the collision frame.
We construct a four-vector transverse to the uid velocity
WT = 1 W 0F
JHEP03(2017)135
(W (W [notdef] u)u) ; (A.2)
which, in the uid frame, has components WT = (0; wF ). We use this vector to express the change in four-momentum associated with the kick as
P [prime] = P EF W T + qe? ; 1 [radicalBigg]
1
q2 E2F w2F
; (A.3)
{ 51 {
where E2F w2F = E2F E2(1 w2) and the four-vector e? satis es the conditionsu [notdef] e? = 0 ; W [notdef] e? = 0 ; e2? = 1 : (A.4)
It is then possible to show that e? can be written as a linear combination of two orthogonal
vectors satisfying the conditions
e1 = (0; w [notdef] v |w [notdef] v[notdef]
) ; e2 = 1
pN l2 + W ?
[parenrightbig]
; (A.5)
with
JHEP03(2017)135
l2 = (0; w
|w[notdef] [notdef]
w [notdef] v |w [notdef] v[notdef]
) ; W? = W
W 2 u [notdef] W
u ; (A.6)
=
(l2 [notdef] u) (u [notdef] W ) (u [notdef] W )2 W 2
; N = (u [notdef] W )2 W 2(1 + (l2 [notdef] u)2)
(u [notdef] W )2 W 2
: (A.7)
These expressions allow us to determine the change in momentum in terms of collision frame quantities alone, once q and the angle in the (e1; e2) plane have been chosen. After each time interval dt in the collision frame, we select a random value of q according to a Gaussian probability distribution with width Q2? = ^
q dtF , where the relation between dtF and dt is dtF = dt F (1 wv). The angle is chosen randomly with a uniform distribution.
This procedure is repeated for each parton in the shower as long as it is in a location where the temperature of the expanding cooling hydrodynamic uid still satis es T > Tc.
The kick that the parton receives is transverse only in the local uid rest frame. In the collision center-of-mass frame, the kick has a transverse component but it also has a longitudinal component which results in energy loss given by
E = E 2F (vw v2) q(2) [notdef]
w[notdef] [notdef]w [notdef] v[notdef]
w2 vw
; (A.8)
where q(2) is the component of the momentum transfer along the e2 direction. (Note that this formula is badly behaved when the velocities of the particle and the uid are identical, since in that case it is not possible to de ne the transverse direction in the uid frame. We have tested that in our simulation this does not occur over the whole propagation of the parton in the plasma.)
B Background subtraction procedure
As we have described in section 4, in order to analyze the observable consequences of the wake in the hydrodynamic medium, we need to add a background of particles that reproduces the measured particle yields and spectra so that we can then incorporate the e ects of the wake as a perturbation on these spectra. So, our Pythia simulations of jet showers, modi ed as in our hybrid model, are now embedded in a (perturbed) background, meaning that we must run a background subtraction algorithm. We have used an iterative noise/pedestal subtraction procedure similar to the one implemented by the experimental collaborations. Since the algorithms employed by ATLAS and CMS are di erent, we have
{ 52 {
employed di erent version of the procedure for the analysis of observables for jets with pT > 100 GeV, where we follow CMS [184], and for the observables that cover jets with pT < 100 GeV, where we follow ATLAS [48]. See the gures in section 4 for the results we have obtained after following these procedures, procedures which we detail below.
For high pT jet observables, our procedure follows the following steps:
(i) Discretize the ( ; ) space in cells of size 0:091 [notdef] 0:087. Sum the transverse energy
ET of all particles falling into the same cell.
(ii) Compute the average transverse energy and the standard deviation in the transverse energy for all the cells in a strip with a given rapidity, i.e. [angbracketleft]ET ( )[angbracketright] and ( )
[radicalBig][angbracketleftbig]
E2T ( )
(iv) Run the anti-kt clustering algorithm using all cells whose ^
EiT is di erent from zero.
Each such cell is given to the anti-kt algorithm as a null four-vector with transverse momentum ^
EiT located at an and corresponding to the geometric center of the cell.
(v) Repeat step (ii) excluding all cells that the anti-kt algorithm has already identi ed as belonging to a jet with transverse energy above EcutT. Jets with energy above this cut are considered signal jets. We now make a second pass, re-evaluating the background outside these jets.
(vi) Repeat the background subtraction (iii) and the jet nding step(iv) using the values obtained in previous step. The jets found by the anti-kt algorithm on this second pass, which may have transverse energies above or below EcutT, are added to the collection of signal jets for this event.
This procedure involves choosing two parameters, namely and EcutT. The best choice for these parameters will in general depend on the anti-kt reconstruction parameter R, the jet pT range under study, and the magnitude of the cell-to-cell uctuations of the background. The factor controls the e ect of background uctuations. For a homogeneous background with no uctuations, one should choose = 0. For a cell-to-cell uctuating background, one needs to increase accordingly, at the risk of potentially removing some of the signal of interest. The value of EcutT determines whether a group of cells corresponds to signal and should therefore be excluded from the background estimation. The choice of these two parameters should be guided by the criterion that the pT of a jet that was arti cially embedded into the background and then reconstructed upon performing the background subtraction procedure above is, on average, as close as possible to the original pT of the
{ 53 {
[angbracketleft]ET ( )[angbracketright]2.
(iii) For each cell, subtract the average for the strip in which the cell is found and subtract a contribution proportional to the standard deviation for that strip, with the proportionality constant a parameter that we shall return to below. If the result is negative, set it to zero instead. That is, for the ith cell in the strip compute:
^
EiT max
EiT [angbracketleft]ET ( )[angbracketright] ( ); 0
[bracketrightbig]
(B.1)
[angbracketrightbig]
JHEP03(2017)135
arti cial jet, within 2% for jets with pT 100 GeV. This is referred to as ensuring that
one has a good jet energy scale (JES). We have done this test using Pythia jets. We choose EcutT = 30 GeV throughout, and have then picked a value of for each centrality and each value of the anti-kt parameter R so as to opimize the JES. For events with 0-10% centrality, we choose = 0:48; 0:72; 1:00; 1:17 for R = 0:2; 0:3; 0:4; 0:5 while for events with 10-30% centrality the corresponding values that we choose are = 0:42; 0:71; 0:94; 1:13.
For jets in the lower pT region, with pT < 100 GeV, the procedure that we adopt begins with the same discretization of ( ; ) space and then continues as follows:
(i) Reconstruct jets using the anti-kt algorithm for many di erent values of the reconstruction parameter R, in each case using the uncorrected the ET in each cell in ( ; ), introducing null four-vectors with transverse momentum ET with a and corresponding to the geometric center of the cell. Because in this procedure we start by reconstructing jets before subtracting background, after subtracting background we must then also remove combinatorial jets.
(ii) Select a set of seed jets with R = 0:2 which have at least one constituent cell with ET > 3 GeV and whose cell with the maximum transverse energy EmaxT satis es EmaxT > 4 [angbracketleft]ET [angbracketright], where [angbracketleft]ET [angbracketright] is the average transverse energy of the cells within the
seed jet.
(iii) Compute the average transverse energy of each rapidity strip, i.e. [angbracketleft]ET ( )[angbracketright], but excluding all those cells that belong to a seed jet.
(iv) Subtract from each cell the average transverse energy for that rapidity strip, ^
EiT =
EiT [angbracketleft]ET ( )[angbracketright]. This subtraction is applied to all cells in a strip, including those in
a seed jet, meaning that it modi es the pT of the seed jets. This is the rst of two subtractions.
(v) The second subtraction, below, will employ a subset of the seed jets with R = 0:2 from above whose transverse momenta after the rst subtraction above satisfy pjetT >
25 GeV. It will also employ a set of seed \track jets" with pjetT > 10 GeV. Track jets are built using only charged tracks with ptrackT > 4 GeV and are reconstructed using the anti-kt algorithm with R = 0:4.
(vi) Recompute the average transverse energy [angbracketleft]E[prime]T ( )[angbracketright], now excluding all those cells that
lie within r = 0:4 from the axis of any of the seed jets or seed track jets, where R is the angular distance in ( ; ) space.
(vii) Subtract the new average energy [angbracketleft]E[prime]T ( )[angbracketright] from all cells and update the kinematics of
all jets reconstructed with all values of R. Keep only those jets with ET > 20 GeV.
(viii) In order to suppress the contribution of combinatorial jets, we impose that the reconstructed jets have to lie within r = 0:2 of a seed track jet (de ned above) with pT > 7 GeV, making it probable that the signal jets include one or a few hard particles in them.
{ 54 {
JHEP03(2017)135
Depending on which observables we want to compute and compare to data, we will need to apply one or more further corrections to the jets in PbPb events that we have extracted via the background subtraction and jet reconstruction procedures detailed above.
For all observables, we will follow the experimental analyses and apply a JES correction which takes into account the remaining average disagreement between the energy of a proton-proton jet from Pythia and the energy of that jet after it has been embedded into one of our PbPb events and reconstructed as above. For the rst of the two procedures above, the one that we employ for jets with pT > 100 GeV, the JES correction is less than 2%, as we noted above. For the procedure that we employ for jets with pT < 100 GeV, the JES correction can be as large as 12% or 20% for R = 0:4 or R = 0:5 jets with pT = 40 GeV.
If we wish to compare to experimental data that has not been unfolded, as we shall do in the case of several high-pT observables measured by CMS, we must include an additional jet energy resolution (JER) correction. We start by doing a Monte Carlo study in which we insert jets from Pythia into an event, t the distribution of the energy reconstruction e ciency (reconstructed jet energy over generator level jet energy as a function of generator level jet energy) for each jet pT bin with a Gaussian and extract the corresponding standard deviation , which tells us how much the reconstructed jet energy varies from event to event and jet to jet. (This is unrelated to the in the background subtraction procedure above; both are conventionally referred to as .) When we do this calculation using the events in our model including its simpli ed background, we denote what we obtain by model. CMS
has done this analysis on data from heavy ion collisions, obtaining LHC which they have tabulated as a function of jet pT in ref. [46]. We are now ready to correct for the fact that the JER in the real background measured by CMS is di erent from that in our simpli ed model background. We do so by smearing the jet energies in our model calculation with a
Gaussian whose width corresponds to extra
q2LHC 2model.
If we wish to compare to data that has been unfolded, we will perform the simplest version of the so called bin-by-bin unfolding. This a ects the jet spectrum measurements. The correction applied consists of multiplying the spectrum of the medium-modi ed jets that we have reconstructed via one of the two procedures above by the ratio of two other spectra: the jet spectrum obtained directly from Pythia divided by the jet spectrum obtained after embedding Pythia jets into heavy ion collision events and reconstructing them.
Acknowledgments
We thank J. Brewer, P. Chesler, Z. Hulcher, P. Jacobs, R. Kunnawalkam Elayavalli, Y.-J. Lee, Y. Mehtar-Tani, A. H. Mueller, G. Qin, A. Sadofyev, K. Tywoniuk, W. van der Schee and K. Zapp for helpful discussions. We also thank J. Balewski, M. Goncharov and C. Paus for their assistance with computing. JCS is a Royal Society University Research Fellow. The work of JCS was also supported by a Ram on y Cajal fellowship. The work of JGM was partly supported by Fundaco para a Ci^encia e a Tecnologia (Portugal) under project CERN/FIS-NUC/0049/2015 and contract \Investigador FCT - Development Grant". DP acknowledges the hospitality of the Center for Theoretical Physics, MIT, where part of this work was completed. JCS and DP acknowledge funding from grants FPA2013-46570 and
{ 55 {
JHEP03(2017)135
MDM-2014-0369 of ICCUB (Unidad de Excelencia Mar a de Maeztu) from the Spanish MINECO, from grant 2014-SGR-104 from the Generalitat de Catalunya and from the Consolider CPAN project. KR acknowledges the hospitality of the CERN Theory division during the completion of this work. The work of KR was supported by the U.S. Department of Energy under grant Contract Number DE-SC0011090.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/
Web End =CC-BY 4.0 ), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.
References
[1] PHENIX collaboration, K. Adcox et al., Formation of dense partonic matter in relativistic nucleus-nucleus collisions at RHIC: Experimental evaluation by the PHENIX collaboration, http://dx.doi.org/10.1016/j.nuclphysa.2005.03.086
Web End =Nucl. Phys. A 757 (2005) 184 [https://arxiv.org/abs/nucl-ex/0410003
Web End =nucl-ex/0410003 ] [http://inspirehep.net/search?p=find+EPRINT+nucl-ex/0410003
Web End =INSPIRE ].
[2] BRAHMS collaboration, I. Arsene et al., Quark gluon plasma and color glass condensate at RHIC? The Perspective from the BRAHMS experiment, http://dx.doi.org/10.1016/j.nuclphysa.2005.02.130
Web End =Nucl. Phys. A 757 (2005) 1 [https://arxiv.org/abs/nucl-ex/0410020
Web End =nucl-ex/0410020 ] [http://inspirehep.net/search?p=find+EPRINT+nucl-ex/0410020
Web End =INSPIRE ].
[3] B.B. Back et al., The PHOBOS perspective on discoveries at RHIC, http://dx.doi.org/10.1016/j.nuclphysa.2005.03.084
Web End =Nucl. Phys. A 757 http://dx.doi.org/10.1016/j.nuclphysa.2005.03.084
Web End =(2005) 28 [https://arxiv.org/abs/nucl-ex/0410022
Web End =nucl-ex/0410022 ] [http://inspirehep.net/search?p=find+EPRINT+nucl-ex/0410022
Web End =INSPIRE ].
[4] STAR collaboration, J. Adams et al., Experimental and theoretical challenges in the search for the quark gluon plasma: The STAR Collaborations critical assessment of the evidence from RHIC collisions, http://dx.doi.org/10.1016/j.nuclphysa.2005.03.085
Web End =Nucl. Phys. A 757 (2005) 102 [https://arxiv.org/abs/nucl-ex/0501009
Web End =nucl-ex/0501009 ] [http://inspirehep.net/search?p=find+EPRINT+nucl-ex/0501009
Web End =INSPIRE ].
[5] ALICE collaboration, Elliptic ow of charged particles in Pb-Pb collisions at 2.76 TeV, http://dx.doi.org/10.1103/PhysRevLett.105.252302
Web End =Phys. Rev. Lett. 105 (2010) 252302 [https://arxiv.org/abs/1011.3914
Web End =arXiv:1011.3914 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1011.3914
Web End =INSPIRE ].
[6] ATLAS collaboration, Measurement of the pseudorapidity and transverse momentum dependence of the elliptic ow of charged particles in lead-lead collisions at psNN = 2:76
TeV with the ATLAS detector, http://dx.doi.org/10.1016/j.physletb.2011.12.056
Web End =Phys. Lett. B 707 (2012) 330 [https://arxiv.org/abs/1108.6018
Web End =arXiv:1108.6018 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1108.6018
Web End =INSPIRE ].
[7] CMS collaboration, Measurement of the elliptic anisotropy of charged particles produced in PbPb collisions at psNN =2.76 TeV, http://dx.doi.org/10.1103/PhysRevC.87.014902
Web End =Phys. Rev. C 87 (2013) 014902 [https://arxiv.org/abs/1204.1409
Web End =arXiv:1204.1409 ]
[http://inspirehep.net/search?p=find+EPRINT+arXiv:1204.1409
Web End =INSPIRE ].
[8] P. Huovinen, P.F. Kolb, U.W. Heinz, P.V. Ruuskanen and S.A. Voloshin, Radial and elliptic ow at RHIC: Further predictions, http://dx.doi.org/10.1016/S0370-2693(01)00219-2
Web End =Phys. Lett. B 503 (2001) 58 [https://arxiv.org/abs/hep-ph/0101136
Web End =hep-ph/0101136 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0101136
Web End =INSPIRE ].
[9] D. Teaney, J. Lauret and E.V. Shuryak, A Hydrodynamic description of heavy ion collisions at the SPS and RHIC, https://arxiv.org/abs/nucl-th/0110037
Web End =nucl-th/0110037 [http://inspirehep.net/search?p=find+EPRINT+nucl-th/0110037
Web End =INSPIRE ].
[10] T. Hirano, U.W. Heinz, D. Kharzeev, R. Lacey and Y. Nara, Hadronic dissipative e ects on elliptic ow in ultrarelativistic heavy-ion collisions, http://dx.doi.org/10.1016/j.physletb.2006.03.060
Web End =Phys. Lett. B 636 (2006) 299 [https://arxiv.org/abs/nucl-th/0511046
Web End =nucl-th/0511046 ] [http://inspirehep.net/search?p=find+EPRINT+nucl-th/0511046
Web End =INSPIRE ].
[11] P. Romatschke and U. Romatschke, Viscosity Information from Relativistic Nuclear Collisions: How Perfect is the Fluid Observed at RHIC?, http://dx.doi.org/10.1103/PhysRevLett.99.172301
Web End =Phys. Rev. Lett. 99 (2007) 172301 [https://arxiv.org/abs/0706.1522
Web End =arXiv:0706.1522 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0706.1522
Web End =INSPIRE ].
{ 56 {
JHEP03(2017)135
[12] M. Luzum and P. Romatschke, Conformal Relativistic Viscous Hydrodynamics: Applications to RHIC results at psNN = 200 GeV, http://dx.doi.org/10.1103/PhysRevC.78.034915
Web End =Phys. Rev. C 78 (2008) 034915 [Erratum ibid. C 79 (2009) 039903] [https://arxiv.org/abs/0804.4015
Web End =arXiv:0804.4015 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0804.4015
Web End =INSPIRE ].
[13] B. Schenke, S. Jeon and C. Gale, Elliptic and triangular ow in event-by-event (3+1)D viscous hydrodynamics, http://dx.doi.org/10.1103/PhysRevLett.106.042301
Web End =Phys. Rev. Lett. 106 (2011) 042301 [https://arxiv.org/abs/1009.3244
Web End =arXiv:1009.3244 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1009.3244
Web End =INSPIRE ].
[14] T. Hirano, P. Huovinen and Y. Nara, Elliptic ow in Pb+Pb collisions at psNN = 2:76 TeV: hybrid model assessment of the rst data, http://dx.doi.org/10.1103/PhysRevC.84.011901
Web End =Phys. Rev. C 84 (2011) http://dx.doi.org/10.1103/PhysRevC.84.011901
Web End =011901 [https://arxiv.org/abs/1012.3955
Web End =arXiv:1012.3955 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1012.3955
Web End =INSPIRE ].
[15] C. Gale, S. Jeon, B. Schenke, P. Tribedy and R. Venugopalan, Event-by-event anisotropic ow in heavy-ion collisions from combined Yang-Mills and viscous uid dynamics, http://dx.doi.org/10.1103/PhysRevLett.110.012302
Web End =Phys. http://dx.doi.org/10.1103/PhysRevLett.110.012302
Web End =Rev. Lett. 110 (2013) 012302 [https://arxiv.org/abs/1209.6330
Web End =arXiv:1209.6330 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1209.6330
Web End =INSPIRE ].
[16] C. Shen, Z. Qiu, H. Song, J. Bernhard, S. Bass and U. Heinz, The iEBE-VISHNU code package for relativistic heavy-ion collisions, http://dx.doi.org/10.1016/j.cpc.2015.08.039
Web End =Comput. Phys. Commun. 199 (2016) 61 [https://arxiv.org/abs/1409.8164
Web End =arXiv:1409.8164 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1409.8164
Web End =INSPIRE ].
[17] C. Shen, J.-F. Paquet, U. Heinz and C. Gale, Photon Emission from a Momentum Anisotropic quark-gluon Plasma, http://dx.doi.org/10.1103/PhysRevC.91.014908
Web End =Phys. Rev. C 91 (2015) 014908 [https://arxiv.org/abs/1410.3404
Web End =arXiv:1410.3404 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1410.3404
Web End =INSPIRE ].
[18] J.E. Bernhard, J.S. Moreland, S.A. Bass, J. Liu and U. Heinz, Applying Bayesian parameter estimation to relativistic heavy-ion collisions: simultaneous characterization of the initial state and quark-gluon plasma medium, http://dx.doi.org/10.1103/PhysRevC.94.024907
Web End =Phys. Rev. C 94 (2016) 024907 [https://arxiv.org/abs/1605.03954
Web End =arXiv:1605.03954 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1605.03954
Web End =INSPIRE ].
[19] CMS collaboration, Observation of Long-Range Near-Side Angular Correlations in Proton-Proton Collisions at the LHC, http://dx.doi.org/10.1007/JHEP09(2010)091
Web End =JHEP 09 (2010) 091 [https://arxiv.org/abs/1009.4122
Web End =arXiv:1009.4122 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1009.4122
Web End =INSPIRE ].
[20] CMS collaboration, Observation of long-range near-side angular correlations in proton-lead collisions at the LHC, http://dx.doi.org/10.1016/j.physletb.2012.11.025
Web End =Phys. Lett. B 718 (2013) 795 [https://arxiv.org/abs/1210.5482
Web End =arXiv:1210.5482 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1210.5482
Web End =INSPIRE ].
[21] ALICE collaboration, Long-range angular correlations on the near and away side in p-Pb collisions at psNN = 5:02 TeV, http://dx.doi.org/10.1016/j.physletb.2013.01.012
Web End =Phys. Lett. B 719 (2013) 29 [https://arxiv.org/abs/1212.2001
Web End =arXiv:1212.2001 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1212.2001
Web End =INSPIRE ].
[22] ATLAS collaboration, Observation of Associated Near-Side and Away-Side Long-Range Correlations in psNN = 5:02 TeV Proton-Lead Collisions with the ATLAS Detector, http://dx.doi.org/10.1103/PhysRevLett.110.182302
Web End =Phys.
http://dx.doi.org/10.1103/PhysRevLett.110.182302
Web End =Rev. Lett. 110 (2013) 182302 [https://arxiv.org/abs/1212.5198
Web End =arXiv:1212.5198 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1212.5198
Web End =INSPIRE ].
[23] PHENIX collaboration, A. Adare et al., Quadrupole Anisotropy in Dihadron Azimuthal Correlations in Central d+Au Collisions at psNN =200 GeV, http://dx.doi.org/10.1103/PhysRevLett.111.212301
Web End =Phys. Rev. Lett. 111 (2013)
http://dx.doi.org/10.1103/PhysRevLett.111.212301
Web End =212301 [https://arxiv.org/abs/1303.1794
Web End =arXiv:1303.1794 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1303.1794
Web End =INSPIRE ].
[24] CMS collaboration, Multiplicity and transverse momentum dependence of two- and four-particle correlations in pPb and PbPb collisions, http://dx.doi.org/10.1016/j.physletb.2013.06.028
Web End =Phys. Lett. B 724 (2013) 213 [https://arxiv.org/abs/1305.0609
Web End =arXiv:1305.0609 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1305.0609
Web End =INSPIRE ].
[25] CMS collaboration, Evidence for Collective Multiparticle Correlations in p-Pb Collisions, http://dx.doi.org/10.1103/PhysRevLett.115.012301
Web End =Phys. Rev. Lett. 115 (2015) 012301 [https://arxiv.org/abs/1502.05382
Web End =arXiv:1502.05382 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1502.05382
Web End =INSPIRE ].
[26] PHENIX collaboration, A. Adare et al., Measurements of elliptic and triangular ow in high-multiplicity 3He+Au collisions at psNN = 200 GeV, http://dx.doi.org/10.1103/PhysRevLett.115.142301
Web End =Phys. Rev. Lett. 115 (2015)
http://dx.doi.org/10.1103/PhysRevLett.115.142301
Web End =142301 [https://arxiv.org/abs/1507.06273
Web End =arXiv:1507.06273 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1507.06273
Web End =INSPIRE ].
{ 57 {
JHEP03(2017)135
[27] ATLAS collaboration, Observation of Long-Range Elliptic Azimuthal Anisotropies inps = 13 and 2.76 TeV pp Collisions with the ATLAS Detector, http://dx.doi.org/10.1103/PhysRevLett.116.172301
Web End =Phys. Rev. Lett. 116 (2016) http://dx.doi.org/10.1103/PhysRevLett.116.172301
Web End =172301 [https://arxiv.org/abs/1509.04776
Web End =arXiv:1509.04776 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1509.04776
Web End =INSPIRE ].
[28] CMS collaboration, Pseudorapidity dependence of long-range two-particle correlations in pPb collisions at psNN = 5:02 TeV, Submitted to: Phys. Rev. C (2016)
[https://arxiv.org/abs/1604.05347
Web End =arXiv:1604.05347 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1604.05347
Web End =INSPIRE ].
[29] CMS collaboration, Evidence for collectivity in pp collisions at the LHC, http://dx.doi.org/10.1016/j.physletb.2016.12.009
Web End =Phys. Lett. B 765 http://dx.doi.org/10.1016/j.physletb.2016.12.009
Web End =(2017) 193 [https://arxiv.org/abs/1606.06198
Web End =arXiv:1606.06198 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1606.06198
Web End =INSPIRE ].
[30] P. Bozek, Collective ow in p-Pb and d-Pd collisions at TeV energies, http://dx.doi.org/10.1103/PhysRevC.85.014911
Web End =Phys. Rev. C 85 http://dx.doi.org/10.1103/PhysRevC.85.014911
Web End =(2012) 014911 [https://arxiv.org/abs/1112.0915
Web End =arXiv:1112.0915 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1112.0915
Web End =INSPIRE ].
[31] J.L. Nagle et al., Exploiting Intrinsic Triangular Geometry in Relativistic He3+Au Collisions to Disentangle Medium Properties, http://dx.doi.org/10.1103/PhysRevLett.113.112301
Web End =Phys. Rev. Lett. 113 (2014) 112301 [https://arxiv.org/abs/1312.4565
Web End =arXiv:1312.4565 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1312.4565
Web End =INSPIRE ].
[32] B. Schenke and R. Venugopalan, Eccentric protons? Sensitivity of ow to system size and shape in p+p, p+Pb and Pb+Pb collisions, http://dx.doi.org/10.1103/PhysRevLett.113.102301
Web End =Phys. Rev. Lett. 113 (2014) 102301 [https://arxiv.org/abs/1405.3605
Web End =arXiv:1405.3605 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1405.3605
Web End =INSPIRE ].
[33] I. Kozlov, M. Luzum, G. Denicol, S. Jeon and C. Gale, Transverse momentum structure of pair correlations as a signature of collective behavior in small collision systems, https://arxiv.org/abs/1405.3976
Web End =arXiv:1405.3976 [http://inspirehep.net/search?p=find+EPRINT+arXiv:1405.3976
Web End =INSPIRE ].
[34] P. Romatschke, Light-Heavy Ion Collisions: A window into pre-equilibrium QCD dynamics?, http://dx.doi.org/10.1140/epjc/s10052-015-3509-3
Web End =Eur. Phys. J. C 75 (2015) 305 [https://arxiv.org/abs/1502.04745
Web End =arXiv:1502.04745 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1502.04745
Web End =INSPIRE ].
[35] P. Bozek and W. Broniowski, Hydrodynamic modeling of 3He-Au collisions at psNN = 200 GeV, http://dx.doi.org/10.1016/j.physletb.2015.05.068
Web End =Phys. Lett. B 747 (2015) 135 [https://arxiv.org/abs/1503.00468
Web End =arXiv:1503.00468 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1503.00468
Web End =INSPIRE ].
[36] M. Habich, G.A. Miller, P. Romatschke and W. Xiang, Testing hydrodynamic descriptions of p+p collisions at ps = 7 TeV, http://dx.doi.org/10.1140/epjc/s10052-016-4237-z
Web End =Eur. Phys. J. C 76 (2016) 408 [https://arxiv.org/abs/1512.05354
Web End =arXiv:1512.05354 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1512.05354
Web End =INSPIRE ].
[37] J.D. Orjuela Koop, R. Belmont, P. Yin and J.L. Nagle, Exploring the Beam Energy Dependence of Flow-Like Signatures in Small System d+Au Collisions, http://dx.doi.org/10.1103/PhysRevC.93.044910
Web End =Phys. Rev. C 93 http://dx.doi.org/10.1103/PhysRevC.93.044910
Web End =(2016) 044910 [https://arxiv.org/abs/1512.06949
Web End =arXiv:1512.06949 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1512.06949
Web End =INSPIRE ].
[38] P.M. Chesler, Colliding shock waves and hydrodynamics in small systems, http://dx.doi.org/10.1103/PhysRevLett.115.241602
Web End =Phys. Rev. Lett. http://dx.doi.org/10.1103/PhysRevLett.115.241602
Web End =115 (2015) 241602 [https://arxiv.org/abs/1506.02209
Web End =arXiv:1506.02209 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1506.02209
Web End =INSPIRE ].
[39] P.M. Chesler, How big are the smallest drops of quark-gluon plasma?, http://dx.doi.org/10.1007/JHEP03(2016)146
Web End =JHEP 03 (2016) 146 [https://arxiv.org/abs/1601.01583
Web End =arXiv:1601.01583 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1601.01583
Web End =INSPIRE ].
[40] P.M. Chesler and L.G. Ya e, Holography and colliding gravitational shock waves in asymptotically AdS5 spacetime, http://dx.doi.org/10.1103/PhysRevLett.106.021601
Web End =Phys. Rev. Lett. 106 (2011) 021601 [https://arxiv.org/abs/1011.3562
Web End =arXiv:1011.3562 ]
[http://inspirehep.net/search?p=find+EPRINT+arXiv:1011.3562
Web End =INSPIRE ].
[41] J. Casalderrey-Solana, M.P. Heller, D. Mateos and W. van der Schee, From full stopping to transparency in a holographic model of heavy ion collisions, http://dx.doi.org/10.1103/PhysRevLett.111.181601
Web End =Phys. Rev. Lett. 111 (2013) http://dx.doi.org/10.1103/PhysRevLett.111.181601
Web End =181601 [https://arxiv.org/abs/1305.4919
Web End =arXiv:1305.4919 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1305.4919
Web End =INSPIRE ].
[42] P.M. Chesler and L.G. Ya e, Numerical solution of gravitational dynamics in asymptotically anti-de Sitter spacetimes, http://dx.doi.org/10.1007/JHEP07(2014)086
Web End =JHEP 07 (2014) 086 [https://arxiv.org/abs/1309.1439
Web End =arXiv:1309.1439 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1309.1439
Web End =INSPIRE ].
{ 58 {
JHEP03(2017)135
[43] ATLAS collaboration, Observation of a Centrality-Dependent Dijet Asymmetry in Lead-Lead Collisions at psNN = 2:77 TeV with the ATLAS Detector at the LHC, http://dx.doi.org/10.1103/PhysRevLett.105.252303
Web End =Phys.
http://dx.doi.org/10.1103/PhysRevLett.105.252303
Web End =Rev. Lett. 105 (2010) 252303 [https://arxiv.org/abs/1011.6182
Web End =arXiv:1011.6182 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1011.6182
Web End =INSPIRE ].
[44] CMS collaboration, Observation and studies of jet quenching in PbPb collisions at nucleon-nucleon center-of-mass energy = 2.76 TeV, http://dx.doi.org/10.1103/PhysRevC.84.024906
Web End =Phys. Rev. C 84 (2011) 024906 [https://arxiv.org/abs/1102.1957
Web End =arXiv:1102.1957 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1102.1957
Web End =INSPIRE ].
[45] CMS collaboration, Jet momentum dependence of jet quenching in PbPb collisions at psNN = 2:76 TeV, http://dx.doi.org/10.1016/j.physletb.2012.04.058
Web End =Phys. Lett. B 712 (2012) 176 [https://arxiv.org/abs/1202.5022
Web End =arXiv:1202.5022 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1202.5022
Web End =INSPIRE ].
[46] CMS collaboration, Studies of jet quenching using isolated-photon+jet correlations in PbPb and pp collisions at psNN = 2:76 TeV, http://dx.doi.org/10.1016/j.physletb.2012.11.003
Web End =Phys. Lett. B 718 (2013) 773 [https://arxiv.org/abs/1205.0206
Web End =arXiv:1205.0206 ]
[http://inspirehep.net/search?p=find+EPRINT+arXiv:1205.0206
Web End =INSPIRE ].
[47] CMS collaboration, Measurement of jet fragmentation into charged particles in pp and PbPb collisions at psNN = 2:76 TeV, http://dx.doi.org/10.1007/JHEP10(2012)087
Web End =JHEP 10 (2012) 087 [https://arxiv.org/abs/1205.5872
Web End =arXiv:1205.5872 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1205.5872
Web End =INSPIRE ].
[48] ATLAS collaboration, Measurement of the jet radius and transverse momentum dependence of inclusive jet suppression in lead-lead collisions at psNN = 2.76 TeV with the
ATLAS detector, http://dx.doi.org/10.1016/j.physletb.2013.01.024
Web End =Phys. Lett. B 719 (2013) 220 [https://arxiv.org/abs/1208.1967
Web End =arXiv:1208.1967 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1208.1967
Web End =INSPIRE ].
[49] CMS collaboration, Nuclear modi cation factor of high transverse momentum jets in PbPb collisions at sqrt(sNN) = 2.76 TeV, http://cds.cern.ch/record/1472722
Web End =CMS-PAS-HIN-12-004 .
[50] ATLAS collaboration, Measurement of the Azimuthal Angle Dependence of Inclusive Jet Yields in Pb+Pb Collisions at psNN = 2:76 TeV with the ATLAS detector, http://dx.doi.org/10.1103/PhysRevLett.111.152301
Web End =Phys. Rev.
http://dx.doi.org/10.1103/PhysRevLett.111.152301
Web End =Lett. 111 (2013) 152301 [https://arxiv.org/abs/1306.6469
Web End =arXiv:1306.6469 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1306.6469
Web End =INSPIRE ].
[51] CMS collaboration, Modi cation of jet shapes in PbPb collisions at psNN = 2:76 TeV,
http://dx.doi.org/10.1016/j.physletb.2014.01.042
Web End =Phys. Lett. B 730 (2014) 243 [https://arxiv.org/abs/1310.0878
Web End =arXiv:1310.0878 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1310.0878
Web End =INSPIRE ].
[52] ALICE collaboration, Measurement of charged jet suppression in Pb-Pb collisions at psNN = 2:76 TeV, http://dx.doi.org/10.1007/JHEP03(2014)013
Web End =JHEP 03 (2014) 013 [https://arxiv.org/abs/1311.0633
Web End =arXiv:1311.0633 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1311.0633
Web End =INSPIRE ].
[53] CMS collaboration, Evidence of b-Jet Quenching in PbPb Collisions at psNN = 2:76 TeV,
http://dx.doi.org/10.1103/PhysRevLett.115.029903
Web End =Phys. Rev. Lett. 113 (2014) 132301 [https://arxiv.org/abs/1312.4198
Web End =arXiv:1312.4198 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1312.4198
Web End =INSPIRE ].
[54] CMS collaboration, Measurement of jet fragmentation in PbPb and pp collisions at psNN = 2:76 TeV, http://dx.doi.org/10.1103/PhysRevC.90.024908
Web End =Phys. Rev. C 90 (2014) 024908 [https://arxiv.org/abs/1406.0932
Web End =arXiv:1406.0932 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1406.0932
Web End =INSPIRE ].
[55] ATLAS collaboration, Measurement of inclusive jet charged-particle fragmentation functions in Pb+Pb collisions at psNN = 2:76 TeV with the ATLAS detector, http://dx.doi.org/10.1016/j.physletb.2014.10.065
Web End =Phys. Lett. B http://dx.doi.org/10.1016/j.physletb.2014.10.065
Web End =739 (2014) 320 [https://arxiv.org/abs/1406.2979
Web End =arXiv:1406.2979 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1406.2979
Web End =INSPIRE ].
[56] ATLAS collaboration, Measurements of the Nuclear Modi cation Factor for Jets in Pb+Pb Collisions at psNN = 2:76 TeV with the ATLAS Detector, http://dx.doi.org/10.1103/PhysRevLett.114.072302
Web End =Phys. Rev. Lett. 114 (2015)
http://dx.doi.org/10.1103/PhysRevLett.114.072302
Web End =072302 [https://arxiv.org/abs/1411.2357
Web End =arXiv:1411.2357 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1411.2357
Web End =INSPIRE ].
[57] ALICE collaboration, Measurement of jet suppression in central Pb-Pb collisions at psNN = 2:76 TeV, http://dx.doi.org/10.1016/j.physletb.2015.04.039
Web End =Phys. Lett. B 746 (2015) 1 [https://arxiv.org/abs/1502.01689
Web End =arXiv:1502.01689 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1502.01689
Web End =INSPIRE ].
[58] ALICE collaboration, Measurement of jet quenching with semi-inclusive hadron-jet distributions in central Pb-Pb collisions at psNN = 2:76 TeV, http://dx.doi.org/10.1007/JHEP09(2015)170
Web End =JHEP 09 (2015) 170
[https://arxiv.org/abs/1506.03984
Web End =arXiv:1506.03984 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1506.03984
Web End =INSPIRE ].
{ 59 {
JHEP03(2017)135
[59] ATLAS collaboration, Measurement of the production of neighbouring jets in lead-lead collisions at psNN = 2:76 TeV with the ATLAS detector, http://dx.doi.org/10.1016/j.physletb.2015.10.059
Web End =Phys. Lett. B 751 (2015) 376 [https://arxiv.org/abs/1506.08656
Web End =arXiv:1506.08656 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1506.08656
Web End =INSPIRE ].
[60] CMS collaboration, Measurement of transverse momentum relative to dijet systems in PbPb and pp collisions at psNN = 2:76 TeV, http://dx.doi.org/10.1007/JHEP01(2016)006
Web End =JHEP 01 (2016) 006 [https://arxiv.org/abs/1509.09029
Web End =arXiv:1509.09029 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1509.09029
Web End =INSPIRE ].
[61] CMS collaboration, Correlations between jets and charged particles in PbPb and pp collisions at psNN = 2:76 TeV, http://dx.doi.org/10.1007/JHEP02(2016)156
Web End =JHEP 02 (2016) 156 [https://arxiv.org/abs/1601.00079
Web End =arXiv:1601.00079 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1601.00079
Web End =INSPIRE ].
[62] PHENIX collaboration, K. Adcox et al., Suppression of hadrons with large transverse momentum in central Au+Au collisions at psNN = 130-GeV, http://dx.doi.org/10.1103/PhysRevLett.88.022301
Web End =Phys. Rev. Lett. 88 (2002)
http://dx.doi.org/10.1103/PhysRevLett.88.022301
Web End =022301 [https://arxiv.org/abs/nucl-ex/0109003
Web End =nucl-ex/0109003 ] [http://inspirehep.net/search?p=find+EPRINT+nucl-ex/0109003
Web End =INSPIRE ].
[63] STAR collaboration, C. Adler et al., Centrality dependence of high pT hadron suppression in Au+Au collisions at psNN = 130-GeV, http://dx.doi.org/10.1103/PhysRevLett.89.202301
Web End =Phys. Rev. Lett. 89 (2002) 202301 [https://arxiv.org/abs/nucl-ex/0206011
Web End =nucl-ex/0206011 ] [http://inspirehep.net/search?p=find+EPRINT+nucl-ex/0206011
Web End =INSPIRE ].
[64] STAR collaboration, C. Adler et al., Disappearance of back-to-back high pT hadron correlations in central Au+Au collisions at psNN = 200 GeV, http://dx.doi.org/10.1103/PhysRevLett.90.082302
Web End =Phys. Rev. Lett. 90 (2003) http://dx.doi.org/10.1103/PhysRevLett.90.082302
Web End =082302 [https://arxiv.org/abs/nucl-ex/0210033
Web End =nucl-ex/0210033 ] [http://inspirehep.net/search?p=find+EPRINT+nucl-ex/0210033
Web End =INSPIRE ].
[65] STAR collaboration, M. Ploskon, Inclusive cross section and correlations of fully reconstructed jets in psNN = 200 GeV Au+Au and p+p collisions, http://dx.doi.org/10.1016/j.nuclphysa.2009.10.095
Web End =Nucl. Phys. A 830 http://dx.doi.org/10.1016/j.nuclphysa.2009.10.095
Web End =(2009) 255C{258C [https://arxiv.org/abs/0908.1799
Web End =arXiv:0908.1799 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0908.1799
Web End =INSPIRE ].
[66] PHENIX collaboration, D.V. Perepelitsa, Reconstructed Jet Results in p + p, d + Au and Cu + Cu collisions at 200 GeV from PHENIX, http://dx.doi.org/10.1016/j.nuclphysa.2012.12.073
Web End =Nucl. Phys. A 910-911 (2013) 425 .
[67] STAR collaboration, L. Adamczyk et al., Jet-Hadron Correlations in psNN = 200 GeV p + p and Central Au + Au Collisions, http://dx.doi.org/10.1103/PhysRevLett.112.122301
Web End =Phys. Rev. Lett. 112 (2014) 122301 [https://arxiv.org/abs/1302.6184
Web End =arXiv:1302.6184 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1302.6184
Web End =INSPIRE ].
[68] STAR collaboration, P.M. Jacobs and A. Schmah, Measurements of jet quenching with semi-inclusive charged jet distributions in Au + Au collisions at psNN = 200 GeV, http://dx.doi.org/10.1016/j.nuclphysa.2016.02.050
Web End =Nucl.
http://dx.doi.org/10.1016/j.nuclphysa.2016.02.050
Web End =Phys. A 956 (2016) 641 [https://arxiv.org/abs/1512.08784
Web End =arXiv:1512.08784 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1512.08784
Web End =INSPIRE ].
[69] A. Adare et al., An Upgrade Proposal from the PHENIX Collaboration, https://arxiv.org/abs/1501.06197
Web End =arXiv:1501.06197 [ http://inspirehep.net/search?p=find+EPRINT+arXiv:1501.06197
Web End =INSPIRE ].
[70] P. Jacobs and X.-N. Wang, Matter in extremis: Ultrarelativistic nuclear collisions at RHIC, http://dx.doi.org/10.1016/j.ppnp.2004.09.001
Web End =Prog. Part. Nucl. Phys. 54 (2005) 443 [https://arxiv.org/abs/hep-ph/0405125
Web End =hep-ph/0405125 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0405125
Web End =INSPIRE ].
[71] J. Casalderrey-Solana and C.A. Salgado, Introductory lectures on jet quenching in heavy ion collisions, Acta Phys. Polon. B 38 (2007) 3731 [https://arxiv.org/abs/0712.3443
Web End =arXiv:0712.3443 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0712.3443
Web End =INSPIRE ].
[72] A. Majumder and M. Van Leeuwen, The Theory and Phenomenology of Perturbative QCD Based Jet Quenching, http://dx.doi.org/10.1016/j.ppnp.2010.09.001
Web End =Prog. Part. Nucl. Phys. 66 (2011) 41 [https://arxiv.org/abs/1002.2206
Web End =arXiv:1002.2206 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1002.2206
Web End =INSPIRE ].
[73] J. Ghiglieri and D. Teaney, Parton energy loss and momentum broadening at NLO in high temperature QCD plasmas, http://dx.doi.org/10.1142/S0218301315300131
Web End =Int. J. Mod. Phys. E 24 (2015) 1530013 [https://arxiv.org/abs/1502.03730
Web End =arXiv:1502.03730 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1502.03730
Web End =INSPIRE ].
[74] J.-P. Blaizot and Y. Mehtar-Tani, Jet Structure in Heavy Ion Collisions, http://dx.doi.org/10.1142/S021830131530012X
Web End =Int. J. Mod. Phys. http://dx.doi.org/10.1142/S021830131530012X
Web End =E 24 (2015) 1530012 [https://arxiv.org/abs/1503.05958
Web End =arXiv:1503.05958 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1503.05958
Web End =INSPIRE ].
[75] G.-Y. Qin and X.-N. Wang, Jet quenching in high-energy heavy-ion collisions, http://dx.doi.org/10.1142/S0218301315300143
Web End =Int. J. Mod. http://dx.doi.org/10.1142/S0218301315300143
Web End =Phys. E 24 (2015) 1530014 [https://arxiv.org/abs/1511.00790
Web End =arXiv:1511.00790 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1511.00790
Web End =INSPIRE ].
{ 60 {
JHEP03(2017)135
[76] K. Zapp, J. Stachel and U.A. Wiedemann, A Local Monte Carlo implementation of the non-abelian Landau-Pomerantschuk-Migdal e ect, http://dx.doi.org/10.1103/PhysRevLett.103.152302
Web End =Phys. Rev. Lett. 103 (2009) 152302 [https://arxiv.org/abs/0812.3888
Web End =arXiv:0812.3888 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0812.3888
Web End =INSPIRE ].
[77] K. Zapp, G. Ingelman, J. Rathsman, J. Stachel and U.A. Wiedemann, A Monte Carlo Model for Jet Quenching, http://dx.doi.org/10.1140/epjc/s10052-009-0941-2
Web End =Eur. Phys. J. C 60 (2009) 617 [https://arxiv.org/abs/0804.3568
Web End =arXiv:0804.3568 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0804.3568
Web End =INSPIRE ].
[78] N. Armesto, L. Cunqueiro and C.A. Salgado, Q-PYTHIA: A Medium-modi ed implementation of nal state radiation, http://dx.doi.org/10.1140/epjc/s10052-009-1133-9
Web End =Eur. Phys. J. C 63 (2009) 679 [https://arxiv.org/abs/0907.1014
Web End =arXiv:0907.1014 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0907.1014
Web End =INSPIRE ].
[79] B. Schenke, C. Gale and S. Jeon, MARTINI: An Event generator for relativistic heavy-ion collisions, http://dx.doi.org/10.1103/PhysRevC.80.054913
Web End =Phys. Rev. C 80 (2009) 054913 [https://arxiv.org/abs/0909.2037
Web End =arXiv:0909.2037 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0909.2037
Web End =INSPIRE ].
[80] I.P. Lokhtin, A.V. Belyaev and A.M. Snigirev, Jet quenching pattern at LHC in PYQUEN model, http://dx.doi.org/10.1140/epjc/s10052-011-1650-1
Web End =Eur. Phys. J. C 71 (2011) 1650 [https://arxiv.org/abs/1103.1853
Web End =arXiv:1103.1853 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1103.1853
Web End =INSPIRE ].
[81] K.C. Zapp, F. Krauss and U.A. Wiedemann, A perturbative framework for jet quenching, http://dx.doi.org/10.1007/JHEP03(2013)080
Web End =JHEP 03 (2013) 080 [https://arxiv.org/abs/1212.1599
Web End =arXiv:1212.1599 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1212.1599
Web End =INSPIRE ].
[82] K.C. Zapp, JEWEL 2.0.0: directions for use, http://dx.doi.org/10.1140/epjc/s10052-014-2762-1
Web End =Eur. Phys. J. C 74 (2014) 2762 [https://arxiv.org/abs/1311.0048
Web End =arXiv:1311.0048 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1311.0048
Web End =INSPIRE ].
[83] K.C. Zapp, Geometrical aspects of jet quenching in JEWEL, http://dx.doi.org/10.1016/j.physletb.2014.06.020
Web End =Phys. Lett. B 735 (2014) 157 [https://arxiv.org/abs/1312.5536
Web End =arXiv:1312.5536 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1312.5536
Web End =INSPIRE ].
[84] I. Vitev and B.-W. Zhang, Jet tomography of high-energy nucleus-nucleus collisions at next-to-leading order, http://dx.doi.org/10.1103/PhysRevLett.104.132001
Web End =Phys. Rev. Lett. 104 (2010) 132001 [https://arxiv.org/abs/0910.1090
Web End =arXiv:0910.1090 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0910.1090
Web End =INSPIRE ].
[85] J. Casalderrey-Solana, J.G. Milhano and U.A. Wiedemann, Jet Quenching via Jet Collimation, http://dx.doi.org/10.1088/0954-3899/38/3/035006
Web End =J. Phys. G 38 (2011) 035006 [https://arxiv.org/abs/1012.0745
Web End =arXiv:1012.0745 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1012.0745
Web End =INSPIRE ].
[86] G.-Y. Qin and B. Muller, Explanation of Di-jet asymmetry in Pb+Pb collisions at the Large Hadron Collider, https://doi.org/10.1103/PhysRevLett.106.162302
Web End =Phys. Rev. Lett. 106 (2011) 162302 [Erratum ibid. 108 (2012) 189904] [https://arxiv.org/abs/1012.5280
Web End =arXiv:1012.5280 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1012.5280
Web End =INSPIRE ].
[87] C. Young, B. Schenke, S. Jeon and C. Gale, Dijet asymmetry at the energies available at the CERN Large Hadron Collider, http://dx.doi.org/10.1103/PhysRevC.84.024907
Web End =Phys. Rev. C 84 (2011) 024907 [https://arxiv.org/abs/1103.5769
Web End =arXiv:1103.5769 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1103.5769
Web End =INSPIRE ].
[88] Y. He, I. Vitev and B.-W. Zhang, O( 3s) Analysis of Inclusive Jet and di-Jet Production in
Heavy Ion Reactions at the Large Hadron Collider, http://dx.doi.org/10.1016/j.physletb.2012.05.054
Web End =Phys. Lett. B 713 (2012) 224 [https://arxiv.org/abs/1105.2566
Web End =arXiv:1105.2566 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1105.2566
Web End =INSPIRE ].
[89] J. Casalderrey-Solana, J.G. Milhano and U. Wiedemann, Jet quenching via jet collimation,http://dx.doi.org/10.1088/0954-3899/38/12/124086
Web End =J. Phys. G 38 (2011) 124086 [https://arxiv.org/abs/1107.1964
Web End =arXiv:1107.1964 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1107.1964
Web End =INSPIRE ].
[90] T. Renk, On the sensitivity of the dijet asymmetry to the physics of jet quenching, http://dx.doi.org/10.1103/PhysRevC.85.064908
Web End =Phys. http://dx.doi.org/10.1103/PhysRevC.85.064908
Web End =Rev. C 85 (2012) 064908 [https://arxiv.org/abs/1202.4579
Web End =arXiv:1202.4579 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1202.4579
Web End =INSPIRE ].
[91] R.B. Neufeld and I. Vitev, The Z0-tagged jet event asymmetry in heavy-ion collisions at the CERN Large Hadron Collider, http://dx.doi.org/10.1103/PhysRevLett.108.242001
Web End =Phys. Rev. Lett. 108 (2012) 242001 [https://arxiv.org/abs/1202.5556
Web End =arXiv:1202.5556 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1202.5556
Web End =INSPIRE ].
[92] T. Renk, Energy dependence of the dijet imbalance in Pb-Pb collisions at 2.76 ATeV, http://dx.doi.org/10.1103/PhysRevC.86.061901
Web End =Phys. http://dx.doi.org/10.1103/PhysRevC.86.061901
Web End =Rev. C 86 (2012) 061901 [https://arxiv.org/abs/1204.5572
Web End =arXiv:1204.5572 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1204.5572
Web End =INSPIRE ].
{ 61 {
JHEP03(2017)135
[93] W. Dai, I. Vitev and B.-W. Zhang, Momentum imbalance of isolated photon-tagged jet production at RHIC and LHC, http://dx.doi.org/10.1103/PhysRevLett.110.142001
Web End =Phys. Rev. Lett. 110 (2013) 142001 [https://arxiv.org/abs/1207.5177
Web End =arXiv:1207.5177 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1207.5177
Web End =INSPIRE ].
[94] L. Apolinario, N. Armesto and L. Cunqueiro, An analysis of the in uence of background subtraction and quenching on jet observables in heavy-ion collisions, http://dx.doi.org/10.1007/JHEP02(2013)022
Web End =JHEP 02 (2013) 022 [https://arxiv.org/abs/1211.1161
Web End =arXiv:1211.1161 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1211.1161
Web End =INSPIRE ].
[95] X.-N. Wang and Y. Zhu, Medium Modi cation of -jets in High-energy Heavy-ion Collisions, http://dx.doi.org/10.1103/PhysRevLett.111.062301
Web End =Phys. Rev. Lett. 111 (2013) 062301 [https://arxiv.org/abs/1302.5874
Web End =arXiv:1302.5874 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1302.5874
Web End =INSPIRE ].
[96] G.-L. Ma, Dijet asymmetry in Pb+Pb collisions at psNN = 2:76 TeV within a multiphase transport model, http://dx.doi.org/10.1103/PhysRevC.87.064901
Web End =Phys. Rev. C 87 (2013) 064901 [https://arxiv.org/abs/1304.2841
Web End =arXiv:1304.2841 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1304.2841
Web End =INSPIRE ].
[97] J. Huang, Z.-B. Kang and I. Vitev, Inclusive b-jet production in heavy ion collisions at the LHC, http://dx.doi.org/10.1016/j.physletb.2013.08.009
Web End =Phys. Lett. B 726 (2013) 251 [https://arxiv.org/abs/1306.0909
Web End =arXiv:1306.0909 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1306.0909
Web End =INSPIRE ].
[98] F. Senzel, O. Fochler, J. Upho , Z. Xu and C. Greiner, In uence of multiple in-medium scattering processes on the momentum imbalance of reconstructed di-jets, http://dx.doi.org/10.1088/0954-3899/42/11/115104
Web End =J. Phys. G 42 http://dx.doi.org/10.1088/0954-3899/42/11/115104
Web End =(2015) 115104 [https://arxiv.org/abs/1309.1657
Web End =arXiv:1309.1657 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1309.1657
Web End =INSPIRE ].
[99] R. Perez-Ramos and T. Renk, In-medium jet shape from energy collimation in parton showers: Comparison with CMS PbPb data at 2.76 TeV, http://dx.doi.org/10.1103/PhysRevD.90.014018
Web End =Phys. Rev. D 90 (2014) 014018 [https://arxiv.org/abs/1401.5283
Web End =arXiv:1401.5283 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1401.5283
Web End =INSPIRE ].
[100] T. Renk, A study of the constraining power of high pT observables in heavy-ion collisions, https://arxiv.org/abs/1408.6684
Web End =arXiv:1408.6684 [http://inspirehep.net/search?p=find+EPRINT+arXiv:1408.6684
Web End =INSPIRE ].
[101] R. P erez-Ramos and T. Renk, A Monte Carlo study of jet fragmentation functions in PbPb and pp collisions at ps = 2:76 TeV, https://arxiv.org/abs/1411.1983
Web End =arXiv:1411.1983 [http://inspirehep.net/search?p=find+EPRINT+arXiv:1411.1983
Web End =INSPIRE ].
[102] Y.-T. Chien, A. Emerman, Z.-B. Kang, G. Ovanesyan and I. Vitev, Jet Quenching fromQCD Evolution, http://dx.doi.org/10.1103/PhysRevD.93.074030
Web End =Phys. Rev. D 93 (2016) 074030 [https://arxiv.org/abs/1509.02936
Web End =arXiv:1509.02936 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1509.02936
Web End =INSPIRE ].
[103] J. Huang, Z.-B. Kang, I. Vitev and H. Xing, Photon-tagged and B-meson-tagged b-jet production at the LHC, http://dx.doi.org/10.1016/j.physletb.2015.09.029
Web End =Phys. Lett. B 750 (2015) 287 [https://arxiv.org/abs/1505.03517
Web End =arXiv:1505.03517 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1505.03517
Web End =INSPIRE ].
[104] Y.-T. Chien and I. Vitev, Towards the understanding of jet shapes and cross sections in heavy ion collisions using soft-collinear e ective theory, http://dx.doi.org/10.1007/JHEP05(2016)023
Web End =JHEP 05 (2016) 023 [https://arxiv.org/abs/1509.07257
Web End =arXiv:1509.07257 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1509.07257
Web End =INSPIRE ].
[105] J.G. Milhano and K.C. Zapp, Origins of the di-jet asymmetry in heavy ion collisions, http://dx.doi.org/10.1140/epjc/s10052-016-4130-9
Web End =Eur. http://dx.doi.org/10.1140/epjc/s10052-016-4130-9
Web End =Phys. J. C 76 (2016) 288 [https://arxiv.org/abs/1512.08107
Web End =arXiv:1512.08107 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1512.08107
Web End =INSPIRE ].
[106] X. Zhang, L. Apolin ario, J.G. Milhano and M. P losko n, Sub-jet structure as a discriminating quenching probe, http://dx.doi.org/10.1016/j.nuclphysa.2016.02.028
Web End =Nucl. Phys. A 956 (2016) 597 [https://arxiv.org/abs/1512.09255
Web End =arXiv:1512.09255 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1512.09255
Web End =INSPIRE ].
[107] N.-B. Chang and G.-Y. Qin, Full jet evolution in quark-gluon plasma and nuclear modi cation of jet production and jet shape in Pb+Pb collisions at 2.76ATeV at the CERN Large Hadron Collider, http://dx.doi.org/10.1103/PhysRevC.94.024902
Web End =Phys. Rev. C 94 (2016) 024902 [https://arxiv.org/abs/1603.01920
Web End =arXiv:1603.01920 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1603.01920
Web End =INSPIRE ].
[108] A.H. Mueller, B. Wu, B.-W. Xiao and F. Yuan, Probing Transverse Momentum Broadening in Heavy Ion Collisions, http://dx.doi.org/10.1016/j.physletb.2016.10.037
Web End =Phys. Lett. B 763 (2016) 208 [https://arxiv.org/abs/1604.04250
Web End =arXiv:1604.04250 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1604.04250
Web End =INSPIRE ].
[109] L. Chen, G.-Y. Qin, S.-Y. Wei, B.-W. Xiao and H.-Z. Zhang, Probing TransverseMomentum Broadening via Dihadron and Hadron-jet Angular Correlations in Relativistic Heavy-ion Collisions, https://arxiv.org/abs/1607.01932
Web End =arXiv:1607.01932 [http://inspirehep.net/search?p=find+EPRINT+arXiv:1607.01932
Web End =INSPIRE ].
{ 62 {
JHEP03(2017)135
[110] J.M. Maldacena, The Large-N limit of superconformal eld theories and supergravity, http://dx.doi.org/10.1023/A:1026654312961
Web End =Int. http://dx.doi.org/10.1023/A:1026654312961
Web End =J. Theor. Phys. 38 (1999) 1113 [https://arxiv.org/abs/hep-th/9711200
Web End =hep-th/9711200 ] [http://inspirehep.net/search?p=find+EPRINT+hep-th/9711200
Web End =INSPIRE ].
[111] J. Casalderrey-Solana, H. Liu, D. Mateos, K. Rajagopal and U.A. Wiedemann,Gauge/String Duality, Hot QCD and Heavy Ion Collisions, https://arxiv.org/abs/1101.0618
Web End =arXiv:1101.0618 [http://inspirehep.net/search?p=find+EPRINT+arXiv:1101.0618
Web End =INSPIRE ].
[112] O. DeWolfe, S.S. Gubser, C. Rosen and D. Teaney, Heavy ions and string theory, http://dx.doi.org/10.1016/j.ppnp.2013.11.001
Web End =Prog. http://dx.doi.org/10.1016/j.ppnp.2013.11.001
Web End =Part. Nucl. Phys. 75 (2014) 86 [https://arxiv.org/abs/1304.7794
Web End =arXiv:1304.7794 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1304.7794
Web End =INSPIRE ].
[113] P.M. Chesler and W. van der Schee, Early thermalization, hydrodynamics and energy loss in AdS/CFT, http://dx.doi.org/10.1142/S0218301315300118
Web End =Int. J. Mod. Phys. E 24 (2015) 1530011 [https://arxiv.org/abs/1501.04952
Web End =arXiv:1501.04952 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1501.04952
Web End =INSPIRE ].
[114] C.P. Herzog, A. Karch, P. Kovtun, C. Kozcaz and L.G. Ya e, Energy loss of a heavy quark moving through N = 4 supersymmetric Yang-Mills plasma, http://dx.doi.org/10.1088/1126-6708/2006/07/013
Web End =JHEP 07 (2006) 013
[https://arxiv.org/abs/hep-th/0605158
Web End =hep-th/0605158 ] [http://inspirehep.net/search?p=find+EPRINT+hep-th/0605158
Web End =INSPIRE ].
[115] H. Liu, K. Rajagopal and U.A. Wiedemann, Calculating the jet quenching parameter fromAdS/CFT, http://dx.doi.org/10.1103/PhysRevLett.97.182301
Web End =Phys. Rev. Lett. 97 (2006) 182301 [https://arxiv.org/abs/hep-ph/0605178
Web End =hep-ph/0605178 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0605178
Web End =INSPIRE ].
[116] J. Casalderrey-Solana and D. Teaney, Heavy quark di usion in strongly coupled N = 4
Yang-Mills, http://dx.doi.org/10.1103/PhysRevD.74.085012
Web End =Phys. Rev. D 74 (2006) 085012 [https://arxiv.org/abs/hep-ph/0605199
Web End =hep-ph/0605199 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0605199
Web End =INSPIRE ].
[117] S.S. Gubser, Drag force in AdS/CFT, http://dx.doi.org/10.1103/PhysRevD.74.126005
Web End =Phys. Rev. D 74 (2006) 126005 [https://arxiv.org/abs/hep-th/0605182
Web End =hep-th/0605182 ][http://inspirehep.net/search?p=find+EPRINT+hep-th/0605182
Web End =INSPIRE ].
[118] H. Liu, K. Rajagopal and U.A. Wiedemann, An AdS/CFT Calculation of Screening in aHot Wind, http://dx.doi.org/10.1103/PhysRevLett.98.182301
Web End =Phys. Rev. Lett. 98 (2007) 182301 [https://arxiv.org/abs/hep-ph/0607062
Web End =hep-ph/0607062 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0607062
Web End =INSPIRE ].
[119] M. Chernico , J.A. Garcia and A. Guijosa, The Energy of a Moving Quark-Antiquark Pair in an N = 4 SYM Plasma, http://dx.doi.org/10.1088/1126-6708/2006/09/068
Web End =JHEP 09 (2006) 068 [https://arxiv.org/abs/hep-th/0607089
Web End =hep-th/0607089 ] [
http://inspirehep.net/search?p=find+EPRINT+hep-th/0607089
Web End =INSPIRE ].
[120] H. Liu, K. Rajagopal and U.A. Wiedemann, Wilson loops in heavy ion collisions and their calculation in AdS/CFT, http://dx.doi.org/10.1088/1126-6708/2007/03/066
Web End =JHEP 03 (2007) 066 [https://arxiv.org/abs/hep-ph/0612168
Web End =hep-ph/0612168 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0612168
Web End =INSPIRE ].
[121] S.S. Gubser, Momentum uctuations of heavy quarks in the gauge-string duality, http://dx.doi.org/10.1016/j.nuclphysb.2007.09.017
Web End =Nucl.
http://dx.doi.org/10.1016/j.nuclphysb.2007.09.017
Web End =Phys. B 790 (2008) 175 [https://arxiv.org/abs/hep-th/0612143
Web End =hep-th/0612143 ] [http://inspirehep.net/search?p=find+EPRINT+hep-th/0612143
Web End =INSPIRE ].
[122] J. Casalderrey-Solana and D. Teaney, Transverse Momentum Broadening of a Fast Quark in a N = 4 Yang-Mills Plasma, http://dx.doi.org/10.1088/1126-6708/2007/04/039
Web End =JHEP 04 (2007) 039 [https://arxiv.org/abs/hep-th/0701123
Web End =hep-th/0701123 ] [
http://inspirehep.net/search?p=find+EPRINT+hep-th/0701123
Web End =INSPIRE ].
[123] P.M. Chesler and L.G. Ya e, The Wake of a quark moving through a strongly-coupled plasma, http://dx.doi.org/10.1103/PhysRevLett.99.152001
Web End =Phys. Rev. Lett. 99 (2007) 152001 [https://arxiv.org/abs/0706.0368
Web End =arXiv:0706.0368 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0706.0368
Web End =INSPIRE ].
[124] S.S. Gubser, S.S. Pufu and A. Yarom, Sonic booms and di usion wakes generated by a heavy quark in thermal AdS/CFT, http://dx.doi.org/10.1103/PhysRevLett.100.012301
Web End =Phys. Rev. Lett. 100 (2008) 012301 [https://arxiv.org/abs/0706.4307
Web End =arXiv:0706.4307 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0706.4307
Web End =INSPIRE ].
[125] P.M. Chesler and L.G. Ya e, The Stress-energy tensor of a quark moving through a strongly-coupled N = 4 supersymmetric Yang-Mills plasma: Comparing hydrodynamics and AdS/CFT, http://dx.doi.org/10.1103/PhysRevD.78.045013
Web End =Phys. Rev. D 78 (2008) 045013 [https://arxiv.org/abs/0712.0050
Web End =arXiv:0712.0050 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0712.0050
Web End =INSPIRE ].
[126] D.M. Hofman and J. Maldacena, Conformal collider physics: Energy and charge correlations, http://dx.doi.org/10.1088/1126-6708/2008/05/012
Web End =JHEP 05 (2008) 012 [https://arxiv.org/abs/0803.1467
Web End =arXiv:0803.1467 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0803.1467
Web End =INSPIRE ].
[127] S.S. Gubser, D.R. Gulotta, S.S. Pufu and F.D. Rocha, Gluon energy loss in the gauge-string duality, http://dx.doi.org/10.1088/1126-6708/2008/10/052
Web End =JHEP 10 (2008) 052 [https://arxiv.org/abs/0803.1470
Web End =arXiv:0803.1470 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0803.1470
Web End =INSPIRE ].
[128] Y. Hatta, E. Iancu and A.H. Mueller, Jet evolution in the N = 4 SYM plasma at strong coupling, http://dx.doi.org/10.1088/1126-6708/2008/05/037
Web End =JHEP 05 (2008) 037 [https://arxiv.org/abs/0803.2481
Web End =arXiv:0803.2481 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0803.2481
Web End =INSPIRE ].
{ 63 {
JHEP03(2017)135
[129] F. Dominguez, C. Marquet, A.H. Mueller, B. Wu and B.-W. Xiao, Comparing energy loss and p-perpendicular | broadening in perturbative QCD with strong coupling N = 4 SYM theory, http://dx.doi.org/10.1016/j.nuclphysa.2008.07.004
Web End =Nucl. Phys. A 811 (2008) 197 [https://arxiv.org/abs/0803.3234
Web End =arXiv:0803.3234 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0803.3234
Web End =INSPIRE ].
[130] P.M. Chesler, K. Jensen and A. Karch, Jets in strongly-coupled N = 4 super Yang-Mills
theory, http://dx.doi.org/10.1103/PhysRevD.79.025021
Web End =Phys. Rev. D 79 (2009) 025021 [https://arxiv.org/abs/0804.3110
Web End =arXiv:0804.3110 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0804.3110
Web End =INSPIRE ].
[131] P.M. Chesler, K. Jensen, A. Karch and L.G. Ya e, Light quark energy loss in strongly-coupled N = 4 supersymmetric Yang-Mills plasma, http://dx.doi.org/10.1103/PhysRevD.79.125015
Web End =Phys. Rev. D 79 (2009) 125015
[https://arxiv.org/abs/0810.1985
Web End =arXiv:0810.1985 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0810.1985
Web End =INSPIRE ].
[132] F. DEramo, H. Liu and K. Rajagopal, Transverse Momentum Broadening and the JetQuenching Parameter, Redux, http://dx.doi.org/10.1103/PhysRevD.84.065015
Web End =Phys. Rev. D 84 (2011) 065015 [https://arxiv.org/abs/1006.1367
Web End =arXiv:1006.1367 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1006.1367
Web End =INSPIRE ].
[133] P. Arnold and D. Vaman, Jet quenching in hot strongly coupled gauge theories revisited: 3-point correlators with gauge-gravity duality, http://dx.doi.org/10.1007/JHEP10(2010)099
Web End =JHEP 10 (2010) 099 [https://arxiv.org/abs/1008.4023
Web End =arXiv:1008.4023 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1008.4023
Web End =INSPIRE ].
[134] P. Arnold and D. Vaman, Jet quenching in hot strongly coupled gauge theories simpli ed,http://dx.doi.org/10.1007/JHEP04(2011)027
Web End =JHEP 04 (2011) 027 [https://arxiv.org/abs/1101.2689
Web End =arXiv:1101.2689 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1101.2689
Web End =INSPIRE ].
[135] P. Arnold and D. Vaman, Some new results for jet stopping in AdS/CFT: long version, http://dx.doi.org/10.1088/0954-3899/38/12/124175
Web End =J. http://dx.doi.org/10.1088/0954-3899/38/12/124175
Web End =Phys. G 38 (2011) 124175 [https://arxiv.org/abs/1106.1680
Web End =arXiv:1106.1680 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1106.1680
Web End =INSPIRE ].
[136] M. Chernico , J.A. Garcia, A. Guijosa and J.F. Pedraza, Holographic Lessons for QuarkDynamics, http://dx.doi.org/10.1088/0954-3899/39/5/054002
Web End =J. Phys. G 39 (2012) 054002 [https://arxiv.org/abs/1111.0872
Web End =arXiv:1111.0872 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1111.0872
Web End =INSPIRE ].
[137] P.M. Chesler, Y.-Y. Ho and K. Rajagopal, Shining a Gluon Beam Through quark-gluonPlasma, http://dx.doi.org/10.1103/PhysRevD.85.126006
Web End =Phys. Rev. D 85 (2012) 126006 [https://arxiv.org/abs/1111.1691
Web End =arXiv:1111.1691 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1111.1691
Web End =INSPIRE ].
[138] P. Arnold, P. Szepietowski and D. Vaman, Coupling dependence of jet quenching in hot strongly-coupled gauge theories, http://dx.doi.org/10.1007/JHEP07(2012)024
Web End =JHEP 07 (2012) 024 [https://arxiv.org/abs/1203.6658
Web End =arXiv:1203.6658 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1203.6658
Web End =INSPIRE ].
[139] P. Arnold, P. Szepietowski, D. Vaman and G. Wong, Tidal stretching of gravitons into classical strings: application to jet quenching with AdS/CFT, http://dx.doi.org/10.1007/JHEP02(2013)130
Web End =JHEP 02 (2013) 130 [https://arxiv.org/abs/1212.3321
Web End =arXiv:1212.3321 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1212.3321
Web End =INSPIRE ].
[140] P.M. Chesler, M. Lekaveckas and K. Rajagopal, Heavy quark energy loss far from equilibrium in a strongly coupled collision, http://dx.doi.org/10.1007/JHEP10(2013)013
Web End =JHEP 10 (2013) 013 [https://arxiv.org/abs/1306.0564
Web End =arXiv:1306.0564 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1306.0564
Web End =INSPIRE ].
[141] A. Ficnar and S.S. Gubser, Finite momentum at string endpoints, http://dx.doi.org/10.1103/PhysRevD.89.026002
Web End =Phys. Rev. D 89 (2014) http://dx.doi.org/10.1103/PhysRevD.89.026002
Web End =026002 [https://arxiv.org/abs/1306.6648
Web End =arXiv:1306.6648 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1306.6648
Web End =INSPIRE ].
[142] A. Ficnar, S.S. Gubser and M. Gyulassy, Shooting String Holography of Jet Quenching atRHIC and LHC, http://dx.doi.org/10.1016/j.physletb.2014.10.016
Web End =Phys. Lett. B 738 (2014) 464 [https://arxiv.org/abs/1311.6160
Web End =arXiv:1311.6160 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1311.6160
Web End =INSPIRE ].
[143] P.M. Chesler and K. Rajagopal, Jet quenching in strongly coupled plasma, http://dx.doi.org/10.1103/PhysRevD.90.025033
Web End =Phys. Rev. D 90 http://dx.doi.org/10.1103/PhysRevD.90.025033
Web End =(2014) 025033 [https://arxiv.org/abs/1402.6756
Web End =arXiv:1402.6756 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1402.6756
Web End =INSPIRE ].
[144] R. Rougemont, A. Ficnar, S. Finazzo and J. Noronha, Energy loss, equilibration and thermodynamics of a baryon rich strongly coupled quark-gluon plasma, http://dx.doi.org/10.1007/JHEP04(2016)102
Web End =JHEP 04 (2016) 102 [https://arxiv.org/abs/1507.06556
Web End =arXiv:1507.06556 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1507.06556
Web End =INSPIRE ].
[145] P.M. Chesler and K. Rajagopal, On the Evolution of Jet Energy and Opening Angle inStrongly Coupled Plasma, http://dx.doi.org/10.1007/JHEP05(2016)098
Web End =JHEP 05 (2016) 098 [https://arxiv.org/abs/1511.07567
Web End =arXiv:1511.07567 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1511.07567
Web End =INSPIRE ].
{ 64 {
JHEP03(2017)135
[146] J. Casalderrey-Solana and A. Ficnar, Holographic Three-Jet Events in Strongly CoupledN = 4 Yang-Mills Plasma, https://arxiv.org/abs/1512.00371
Web End =arXiv:1512.00371 [http://inspirehep.net/search?p=find+EPRINT+arXiv:1512.00371
Web End =INSPIRE ].
[147] K. Rajagopal, A.V. Sadofyev and W. van der Schee, Evolution of the jet opening angle distribution in holographic plasma, http://dx.doi.org/10.1103/PhysRevLett.116.211603
Web End =Phys. Rev. Lett. 116 (2016) 211603 [https://arxiv.org/abs/1602.04187
Web End =arXiv:1602.04187 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1602.04187
Web End =INSPIRE ].
[148] J. Casalderrey-Solana, D.C. Gulhan, J.G. Milhano, D. Pablos and K. Rajagopal, A HybridStrong/Weak Coupling Approach to Jet Quenching, http://dx.doi.org/10.1007/JHEP10(2014)019
Web End =JHEP 10 (2014) 019 [Erratum ibid. 1509 (2015) 175] [https://arxiv.org/abs/1405.3864
Web End =arXiv:1405.3864 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1405.3864
Web End =INSPIRE ].
[149] J. Casalderrey-Solana, D.C. Gulhan, J.G. Milhano, D. Pablos and K. Rajagopal, Predictions for Boson-Jet Observables and Fragmentation Function Ratios from a Hybrid Strong/Weak Coupling Model for Jet Quenching, http://dx.doi.org/10.1007/JHEP03(2016)053
Web End =JHEP 03 (2016) 053 [https://arxiv.org/abs/1508.00815
Web End =arXiv:1508.00815 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1508.00815
Web End =INSPIRE ].
[150] J.-P. Blaizot, E. Iancu and Y. Mehtar-Tani, Medium-induced QCD cascade: democratic branching and wave turbulence, http://dx.doi.org/10.1103/PhysRevLett.111.052001
Web End =Phys. Rev. Lett. 111 (2013) 052001 [https://arxiv.org/abs/1301.6102
Web End =arXiv:1301.6102 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1301.6102
Web End =INSPIRE ].
[151] A. Kurkela and U.A. Wiedemann, Picturing perturbative parton cascades in QCD matter,http://dx.doi.org/10.1016/j.physletb.2014.11.054
Web End =Phys. Lett. B 740 (2015) 172 [https://arxiv.org/abs/1407.0293
Web End =arXiv:1407.0293 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1407.0293
Web End =INSPIRE ].
[152] J.-P. Blaizot, Y. Mehtar-Tani and M.A.C. Torres, Angular structure of the in-medium QCD cascade, http://dx.doi.org/10.1103/PhysRevLett.114.222002
Web End =Phys. Rev. Lett. 114 (2015) 222002 [https://arxiv.org/abs/1407.0326
Web End =arXiv:1407.0326 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1407.0326
Web End =INSPIRE ].
[153] L. Fister and E. Iancu, Medium-induced jet evolution: wave turbulence and energy loss,http://dx.doi.org/10.1007/JHEP03(2015)082
Web End =JHEP 03 (2015) 082 [https://arxiv.org/abs/1409.2010
Web End =arXiv:1409.2010 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1409.2010
Web End =INSPIRE ].
[154] J.-P. Blaizot, L. Fister and Y. Mehtar-Tani, Angular distribution of medium-induced QCD cascades, http://dx.doi.org/10.1016/j.nuclphysa.2015.03.014
Web End =Nucl. Phys. A 940 (2015) 67 [https://arxiv.org/abs/1409.6202
Web End =arXiv:1409.6202 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1409.6202
Web End =INSPIRE ].
[155] J.-P. Blaizot and Y. Mehtar-Tani, Energy ow along the medium-induced parton cascade,http://dx.doi.org/10.1016/j.aop.2016.01.002
Web End =Annals Phys. 368 (2016) 148 [https://arxiv.org/abs/1501.03443
Web End =arXiv:1501.03443 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1501.03443
Web End =INSPIRE ].
[156] E. Iancu and B. Wu, Thermalization of mini-jets in a quark-gluon plasma, http://dx.doi.org/10.1007/JHEP10(2015)155
Web End =JHEP 10 (2015) http://dx.doi.org/10.1007/JHEP10(2015)155
Web End =155 [https://arxiv.org/abs/1506.07871
Web End =arXiv:1506.07871 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1506.07871
Web End =INSPIRE ].
[157] R.B. Neufeld and I. Vitev, Parton showers as sources of energy-momentum deposition in the QGP and their implication for shockwave formation at RHIC and at the LHC, http://dx.doi.org/10.1103/PhysRevC.86.024905
Web End =Phys. http://dx.doi.org/10.1103/PhysRevC.86.024905
Web End =Rev. C 86 (2012) 024905 [https://arxiv.org/abs/1105.2067
Web End =arXiv:1105.2067 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1105.2067
Web End =INSPIRE ].
[158] Y. Tachibana and T. Hirano, Momentum transport away from a jet in an expanding nuclear medium, http://dx.doi.org/10.1103/PhysRevC.90.021902
Web End =Phys. Rev. C 90 (2014) 021902 [https://arxiv.org/abs/1402.6469
Web End =arXiv:1402.6469 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1402.6469
Web End =INSPIRE ].
[159] S. Floerchinger and K.C. Zapp, Hydrodynamics and Jets in Dialogue, http://dx.doi.org/10.1140/epjc/s10052-014-3189-4
Web End =Eur. Phys. J. C 74 http://dx.doi.org/10.1140/epjc/s10052-014-3189-4
Web End =(2014) 3189 [https://arxiv.org/abs/1407.1782
Web End =arXiv:1407.1782 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1407.1782
Web End =INSPIRE ].
[160] Y. He, T. Luo, X.-N. Wang and Y. Zhu, Linear Boltzmann Transport for Jet Propagation in the quark-gluon Plasma: Elastic Processes and Medium Recoil, http://dx.doi.org/10.1103/PhysRevC.91.054908
Web End =Phys. Rev. C 91 (2015) http://dx.doi.org/10.1103/PhysRevC.91.054908
Web End =054908 [https://arxiv.org/abs/1503.03313
Web End =arXiv:1503.03313 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1503.03313
Web End =INSPIRE ].
[161] Y. Tachibana and T. Hirano, Interplay between Mach cone and radial expansion and its signal in -jet events, http://dx.doi.org/10.1103/PhysRevC.93.054907
Web End =Phys. Rev. C 93 (2016) 054907 [https://arxiv.org/abs/1510.06966
Web End =arXiv:1510.06966 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1510.06966
Web End =INSPIRE ].
[162] S. Cao, T. Luo, G.-Y. Qin and X.-N. Wang, Linearized Boltzmann transport model for jet propagation in the quark-gluon plasma: Heavy quark evolution, http://dx.doi.org/10.1103/PhysRevC.94.014909
Web End =Phys. Rev. C 94 (2016) http://dx.doi.org/10.1103/PhysRevC.94.014909
Web End =014909 [https://arxiv.org/abs/1605.06447
Web End =arXiv:1605.06447 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1605.06447
Web End =INSPIRE ].
{ 65 {
JHEP03(2017)135
[163] J. Casalderrey-Solana, J.G. Milhano and P. Quiroga-Arias, Out of Medium Fragmentation from Long-Lived Jet Showers, http://dx.doi.org/10.1016/j.physletb.2012.02.066
Web End =Phys. Lett. B 710 (2012) 175 [https://arxiv.org/abs/1111.0310
Web End =arXiv:1111.0310 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1111.0310
Web End =INSPIRE ].
[164] T. Sjostrand, S. Mrenna and P.Z. Skands, A Brief Introduction to PYTHIA 8.1, http://dx.doi.org/10.1016/j.cpc.2008.01.036
Web End =Comput. http://dx.doi.org/10.1016/j.cpc.2008.01.036
Web End =Phys. Commun. 178 (2008) 852 [https://arxiv.org/abs/0710.3820
Web End =arXiv:0710.3820 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0710.3820
Web End =INSPIRE ].
[165] JET collaboration, K.M. Burke et al., Extracting the jet transport coe cient from jet quenching in high-energy heavy-ion collisions, http://dx.doi.org/10.1103/PhysRevC.90.014909
Web End =Phys. Rev. C 90 (2014) 014909 [https://arxiv.org/abs/1312.5003
Web End =arXiv:1312.5003 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1312.5003
Web End =INSPIRE ].
[166] J. Casalderrey-Solana, Y. Mehtar-Tani, C.A. Salgado and K. Tywoniuk, New picture of jet quenching dictated by color coherence, http://dx.doi.org/10.1016/j.physletb.2013.07.046
Web End =Phys. Lett. B 725 (2013) 357 [https://arxiv.org/abs/1210.7765
Web End =arXiv:1210.7765 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1210.7765
Web End =INSPIRE ].
[167] M. Cacciari, G.P. Salam and G. Soyez, The Anti-kt jet clustering algorithm, http://dx.doi.org/10.1088/1126-6708/2008/04/063
Web End =JHEP 04 http://dx.doi.org/10.1088/1126-6708/2008/04/063
Web End =(2008) 063 [https://arxiv.org/abs/0802.1189
Web End =arXiv:0802.1189 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0802.1189
Web End =INSPIRE ].
[168] M. Cacciari, G.P. Salam and G. Soyez, FastJet User Manual, http://dx.doi.org/10.1140/epjc/s10052-012-1896-2
Web End =Eur. Phys. J. C 72 (2012) http://dx.doi.org/10.1140/epjc/s10052-012-1896-2
Web End =1896 [https://arxiv.org/abs/1111.6097
Web End =arXiv:1111.6097 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1111.6097
Web End =INSPIRE ].
[169] CMS collaboration, Decomposing energy balance contributions for quenched jets in PbPb versus pp collisions at psNN = 2:76 TeV, http://cds.cern.ch/record/2055633
Web End =CMS-PAS-HIN-15-011 .
[170] F. DEramo, M. Lekaveckas, H. Liu and K. Rajagopal, Momentum Broadening in Weakly Coupled quark-gluon Plasma (with a view to nding the quasiparticles within liquid quark-gluon plasma), http://dx.doi.org/10.1007/JHEP05(2013)031
Web End =JHEP 05 (2013) 031 [https://arxiv.org/abs/1211.1922
Web End =arXiv:1211.1922 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1211.1922
Web End =INSPIRE ].
[171] H. Stoecker, Collective ow signals the quark gluon plasma, http://dx.doi.org/10.1016/j.nuclphysa.2004.12.074
Web End =Nucl. Phys. A 750 (2005) 121 [https://arxiv.org/abs/nucl-th/0406018
Web End =nucl-th/0406018 ] [http://inspirehep.net/search?p=find+EPRINT+nucl-th/0406018
Web End =INSPIRE ].
[172] J. Casalderrey-Solana, E.V. Shuryak and D. Teaney, Conical ow induced by quenchedQCD jets, http://dx.doi.org/10.1016/j.nuclphysa.2006.06.091
Web End =J. Phys. Conf. Ser. 27 (2005) 22 [https://arxiv.org/abs/hep-ph/0411315
Web End =hep-ph/0411315 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0411315
Web End =INSPIRE ].
[173] L.M. Satarov, H. Stoecker and I.N. Mishustin, Mach shocks induced by partonic jets in expanding quark-gluon plasma, http://dx.doi.org/10.1016/j.physletb.2005.08.102
Web End =Phys. Lett. B 627 (2005) 64 [https://arxiv.org/abs/hep-ph/0505245
Web End =hep-ph/0505245 ] [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0505245
Web End =INSPIRE ].
[174] J. Casalderrey-Solana, E.V. Shuryak and D. Teaney, Hydrodynamic ow from fast particles, https://arxiv.org/abs/hep-ph/0602183
Web End =hep-ph/0602183 [http://inspirehep.net/search?p=find+EPRINT+hep-ph/0602183
Web End =INSPIRE ].
[175] B. Betz, M. Gyulassy, D.H. Rischke, H. Stocker and G. Torrieri, Jet Propagation and MachCones in (3+1)d Ideal Hydrodynamics, http://dx.doi.org/10.1088/0954-3899/35/10/104106
Web End =J. Phys. G 35 (2008) 104106 [https://arxiv.org/abs/0804.4408
Web End =arXiv:0804.4408 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0804.4408
Web End =INSPIRE ].
[176] B. Betz, J. Noronha, G. Torrieri, M. Gyulassy, I. Mishustin and D.H. Rischke, Universality of the Di usion Wake from Stopped and Punch-Through Jets in Heavy-Ion Collisions, http://dx.doi.org/10.1103/PhysRevC.79.034902
Web End =Phys. http://dx.doi.org/10.1103/PhysRevC.79.034902
Web End =Rev. C 79 (2009) 034902 [https://arxiv.org/abs/0812.4401
Web End =arXiv:0812.4401 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0812.4401
Web End =INSPIRE ].
[177] F. Cooper and G. Frye, Comment on the Single Particle Distribution in the Hydrodynamic and Statistical Thermodynamic Models of Multiparticle Production, http://dx.doi.org/10.1103/PhysRevD.10.186
Web End =Phys. Rev. D 10 (1974) http://dx.doi.org/10.1103/PhysRevD.10.186
Web End =186 [http://inspirehep.net/search?p=find+J+%22Phys.Rev.,D10,186%22
Web End =INSPIRE ].
[178] ALICE collaboration, Centrality dependence of , K, p production in Pb-Pb collisions at psNN = 2:76 TeV, http://dx.doi.org/10.1103/PhysRevC.88.044910
Web End =Phys. Rev. C 88 (2013) 044910 [https://arxiv.org/abs/1303.0737
Web End =arXiv:1303.0737 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1303.0737
Web End =INSPIRE ].
[179] A. Beraudo, J.G. Milhano and U.A. Wiedemann, Medium-induced color ow softens hadronization, http://dx.doi.org/10.1103/PhysRevC.85.031901
Web End =Phys. Rev. C 85 (2012) 031901 [https://arxiv.org/abs/1109.5025
Web End =arXiv:1109.5025 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1109.5025
Web End =INSPIRE ].
[180] A. Beraudo, J.G. Milhano and U.A. Wiedemann, The Contribution of Medium-Modi ed Color Flow to Jet Quenching, http://dx.doi.org/10.1007/JHEP07(2012)144
Web End =JHEP 07 (2012) 144 [https://arxiv.org/abs/1204.4342
Web End =arXiv:1204.4342 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1204.4342
Web End =INSPIRE ].
{ 66 {
JHEP03(2017)135
[181] M. Cacciari and G.P. Salam, Pileup subtraction using jet areas, http://dx.doi.org/10.1016/j.physletb.2007.09.077
Web End =Phys. Lett. B 659 (2008) http://dx.doi.org/10.1016/j.physletb.2007.09.077
Web End =119 [https://arxiv.org/abs/0707.1378
Web End =arXiv:0707.1378 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0707.1378
Web End =INSPIRE ].
[182] M. Cacciari, G.P. Salam and G. Soyez, The Catchment Area of Jets, http://dx.doi.org/10.1088/1126-6708/2008/04/005
Web End =JHEP 04 (2008) 005 [https://arxiv.org/abs/0802.1188
Web End =arXiv:0802.1188 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0802.1188
Web End =INSPIRE ].
[183] M. Cacciari, J. Rojo, G.P. Salam and G. Soyez, Jet Reconstruction in Heavy Ion Collisions,http://dx.doi.org/10.1140/epjc/s10052-011-1539-z
Web End =Eur. Phys. J. C 71 (2011) 1539 [https://arxiv.org/abs/1010.1759
Web End =arXiv:1010.1759 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1010.1759
Web End =INSPIRE ].
[184] O. Kodolova, I. Vardanian, A. Nikitenko and A. Oulianov, The performance of the jet identi cation and reconstruction in heavy ions collisions with CMS detector, http://dx.doi.org/10.1140/epjc/s10052-007-0223-9
Web End =Eur. Phys. J. http://dx.doi.org/10.1140/epjc/s10052-007-0223-9
Web End =C 50 (2007) 117 [http://inspirehep.net/search?p=find+J+%22Eur.Phys.J.,C50,117%22
Web End =INSPIRE ].
[185] S.D. Ellis, J. Huston, K. Hatakeyama, P. Loch and M. Tonnesmann, Jets in hadron-hadron collisions, http://dx.doi.org/10.1016/j.ppnp.2007.12.002
Web End =Prog. Part. Nucl. Phys. 60 (2008) 484 [https://arxiv.org/abs/0712.2447
Web End =arXiv:0712.2447 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:0712.2447
Web End =INSPIRE ].
[186] D. dEnterria, K.J. Eskola, I. Helenius and H. Paukkunen, Confronting current NLO parton fragmentation functions with inclusive charged-particle spectra at hadron colliders, http://dx.doi.org/10.1016/j.nuclphysb.2014.04.006
Web End =Nucl. http://dx.doi.org/10.1016/j.nuclphysb.2014.04.006
Web End =Phys. B 883 (2014) 615 [https://arxiv.org/abs/1311.1415
Web End =arXiv:1311.1415 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1311.1415
Web End =INSPIRE ].
[187] Y. Mehtar-Tani and K. Tywoniuk, Jet (de)coherence in Pb-Pb collisions at the LHC, http://dx.doi.org/10.1016/j.physletb.2015.03.041
Web End =Phys. http://dx.doi.org/10.1016/j.physletb.2015.03.041
Web End =Lett. B 744 (2015) 284 [https://arxiv.org/abs/1401.8293
Web End =arXiv:1401.8293 ] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1401.8293
Web End =INSPIRE ].
JHEP03(2017)135
{ 67 {
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Journal of High Energy Physics is a copyright of Springer, 2017.
Abstract
(ProQuest: ... denotes formulae and/or non-USASCII text omitted; see image)
Abstract
Within the context of a hybrid strong/weak coupling model of jet quenching, we study the modification of the angular distribution of the energy within jets in heavy ion collisions, as partons within jet showers lose energy and get kicked as they traverse the strongly coupled plasma produced in the collision. To describe the dynamics transverse to the jet axis, we add the effects of transverse momentum broadening into our hybrid construction, introducing a parameter ...... that governs its magnitude. We show that, because of the quenching of the energy of partons within a jet, even when K≠0 the jets that survive with some specified energy in the final state are narrower than jets with that energy in proton-proton collisions. For this reason, many standard observables are rather insensitive to K. We propose a new differential jet shape ratio observable in which the effects of transverse momentum broadening are apparent. We also analyze the response of the medium to the passage of the jet through it, noting that the momentum lost by the jet appears as the momentum of a wake in the medium. After freezeout this wake becomes soft particles with a broad angular distribution but with net momentum in the jet direction, meaning that the wake contributes to what is reconstructed as a jet. This effect must therefore be included in any description of the angular structure of the soft component of a jet. We show that the particles coming from the response of the medium to the momentum and energy deposited in it leads to a correlation between the momentum of soft particles well separated from the jet in angle with the direction of the jet momentum, and find qualitative but not quantitative agreement with experimental data on observables designed to extract such a correlation. More generally, by confronting the results that we obtain upon introducing transverse momentum broadening and the response of the medium to the jet with available jet data, we highlight the importance of these processes for understanding the internal, soft, angular structure of high energy jets.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer