Parton shower evolution in medium and nuclear modification of photon-tagged jets in Pb+Pb collisions at the LHC

We study the medium modification of jets correlated with large transverse momentum photons at the LHC via a transport and perturbative QCD hybrid model which incorporates the contributions from both elastic collisions and radiative energy loss experienced by the parton showers. Calculations are performed for the modification of the photon-tagged jet yield, the photon-jet energy imbalance, and the azimuthal distribution of away-side jets. The modifications of photon-tagged jets with different $x_T = p_{T, J}/p_{T, \gamma}$ values are studied and they exhibit different centrality and jet cone size dependence due to traversing different medium lengths and density profiles. We further investigate the influence of transverse and longitudinal jet transport coefficients on the nuclear modification of photon-tagged jet production and jet shape observables.


Introduction
The energetic probes produced in early hard scattering processes provide valuable tools for studying the properties of the quark-gluon plasma (QGP) created in relativistic heavy-ion collisions. Large transverse momentum (p T ) jets are particularly interesting as they directly interact with the constituents of dense QCD matter, and are modified in the process [1]. Jet quenching has been confirmed by multiple experimental observations, such as the suppression of single inclusive hadron production at high p T in nucleus-nucleus collisions compared to binary collision scaled proton-proton collisions at the same energies [2,3,4].
A lot of theoretical effort has been devoted to understand the modification of jets in hot and dense stronginteraction QCD matter. The medium modification of jets originates from a combination of elastic collisions with the medium constituents and induced gluon emission while jets propagate through such highly excited nuclear medium [5,6,7,8]. There are currently a few perturbative QCD based schemes for calculating radiative and collisional parton energy loss [9,10,11,12,13,5,6,14,15,16]. A systematic comparison of these different jet quenching schemes has been carried out in Ref. [17] in the context of the "brick" problem. Various phenomenological studies have been performed for the suppression of single inclusive high p T hadron production [18,19,16,20,21], nuclear modification of dihadrons [22,23], photon-hadron [24,25,26] correlations in high energy nucleus-nucleus collisions, etc.
In recent years, sophisticated experimental techniques have been utilized to reconstruct full jets emitted in relativistic heavy-ion collisions, providing both challenges and opportunities for our understanding of jet modification in dense nuclear matter in terms of both leading and subleading fragments of the jet showers [27]. A few Monte-Carlo generators have also been developed to simulate the modification of full jet propagation through dense nuclear medium [28,29,30,31,32,33,34].
Meanwhile, the large kinematics available at the LHC allows us to investigate medium effects on jets with transverse energies over a hundred GeV. The strong modification of the momentum imbalance distribution between the correlated jet pairs first measured in Pb+Pb collisions at the LHC [35,36] indicated that the awayside subleading jets experienced significant energy loss during their propagation through the produced hot and dense medium. These results have triggered wide interest in studying full jets and their in-medium evolution and modification as opposed to single particle observables. Many calculations have suggested that significant jet cone energy loss and the observed dijet asymmetry may come from the loss of low energy gluons from the jet cone to the medium and are then not recovered by the jet algorithm. These gluons escape either by deflection or by energy dissipation and thermalization due to the interaction with the medium constituents [37,38,39,40,41,42,43,44].
Another exciting result is the measurement of full jets correlated with high p T photons in Pb+Pb collisions at the LHC [45]. Jets tagged by photons and electro-weak bosons have been regarded as the "golden" channel for studying jet quenching due to the fact that the boson triggers, once produced, will escape the medium without further interaction. Thus they are expected to put a stronger constraint on the momentum distributions of the away-side tagged jets than hadron or jet triggers [46]. Strong modification of the photon-jet momentum imbalance distribution has been obtained from model calculations as well [47,48].
The goal of various jet quenching studies is to acquire a better understanding of jet-medium interaction and obtain the transport properties of the hot and dense matter created in relativistic nuclear collisions. Much of current attention has been focused on the quantitative extraction of jet transport coefficients, such aŝ [8,49,50]. These coefficients not only control how energy and momentum are lost from jets and dumped into the medium [51,52,53], but are necessary ingredients for the investigation of medium response to jet propagation. The knowledge of these coefficients and their (relative) sizes may provide much insight into the internal structures and the properties of the hot and dense matter probed by the high energy jets.
With this motivation, we develop a Monte-Carlo transport and perturbative QCD hybrid model to simulate jet shower evolution and jet modification in dense QCD medium and apply it to study the medium modification of jets correlated with high p T photons at the LHC energies. The effect of elastic collisions on the parton shower is encoded by the above mentioned transport coefficients, namely the longitudinal drag and diffusion as well as transverse momentum broadening. The same transport coefficients also control the radiative energy loss which is simulated by employing the single gluon emission spectrum obtained from the perturbative QCD-based higher twist jet energy loss calculation [12,13].
Our paper is organized as follows. In section II, we provide some details about how we simulate the evolution and energy loss of jet showers in medium. The contributions of elastic collisions and gluon radiation to jet energy loss are compared. The numerical calculations for photon-jet correlations in Pb+Pb collisions at the LHC energies are presented in section III, including the nuclear modification of tagged jet yield, the photon-jet momentum imbalance and the azimuthal distribution of the away-side jets. To study the medium effects on jets traversing different medium lengths and density profiles, we compare the nuclear modification of photontagged jets with different values of x T = p T,J /p T,γ . We further investigate the influence of longitudinal and transverse jet transport coefficients on the nuclear modification of photon-tagged jets. The last section contains our summary.

Parton shower evolution in medium
In this section, we describe our Monte-Carlo method for simulating the transport and modification of a partonic jet shower in dense QCD matter. We assume each parton of the jet shower follows classical trajectories between each jet-medium interaction that induces momentum change of the parton (either by an elastic collision or the radiation of a gluon). We keep track of all the propagating partons and the radiative gluons until they exit the dense medium.
The probability for a parton at each time step ∆t to have an momentum exchange by an elastic collision with the medium constituents is determined by, where t mf is an "effective" mean free time, thus 1/t mf is the rate of the elastic collisions (the exchange of energy and momentum with medium). If there is a jet-medium interaction through an elastic collision, the energy and momentum exchange is sampled from a probability distribution of P (∆E, ∆p ⊥ , t, ∆t). Generally, the above mean free time t mf and the distribution of the momentum and energy exchange can be determined if the details about the medium structure and jet-medium interaction are known. For example, the mean free time can be obtained from the density ρ of the medium and the cross section σ, i.e., t mf = 1/(ρσ). In this work, we take the multiple soft scattering limit in which the momentum exchange through each elastic collision is small. In this limit, the probability distribution P (∆E, ∆p ⊥ , t, ∆t) is a Gaussian for both longitudinal and transverse directions. Thus we may set the "effective" mean free time in our simulation equal to the time step (t mf = ∆t), i.e., in each time step, the parton will experience a momentum exchange due to a soft elastic scattering. The mean and variance of the distribution are then determined by the longitudinal dragêt mf , longitudinal diffusionê 2 t mf (and transverse broadeningqt mf ). We have checked that the variation of t mf values does not affect the results if a Gaussian form is used for the momentum exchange in a sufficiently long medium (∆t ≤ t mf ≪ L). We note that within our framework one may go beyond the multiple scattering limit, such as a few hard scatterings, with the use of more realistic values of t mf and non-Gaussian form of momentum exchange [54,55]. We leave this to a future study.
We simulate the radiation of gluons by employing the single gluon emission spectrum from higher twist calculations [12,13], where x and l ⊥ are the fractional energy and transverse momentum carried by the radiated gluon, P (x) is the vacuum parton splitting function,q A is the gluon transport coefficient, and t form = 2Ex(1 − x)/l 2 ⊥ is the gluon formation time. In the above equation, the gluon radiation is purely stimulated by the transverse momentum broadening of the jet. Generally, the longitudinal momentum exchange (drag and diffusion) may give rise to medium-induced gluon radiation as well. Such contribution is straightforward to include in the current Monte-Carlo simulation when it becomes available [56].
The probability for a parton in each time step ∆t to radiate a gluon is determined from medium-induced gluon radiation spectrum: We choose the value of ∆t to ensure that the average radiated gluon number is smaller than 1 in a time step ∆t (0.01 fm in our simulation). For N g (t, ∆t) ≪ 1, our treatment is equivalent to the Poisson distribution for the gluon emission (1 − e − Ng ≈ N g ). If there is a gluon emitted in a given time step, the energy and the momentum of the radiated gluon are sampled from the following probability distribution, When simulating jet shower propagating through the hot and dense medium, we allow the radiated gluons to interact with the medium after their formation times t form are reached. We impose a minimum energy cutoff πT for the medium-induced radiation spectrum in order to take into account the balance between gluon radiation and absorption from the hydrodynamic medium.
To simulate multiple gluon emission, the initial time t i in the above equation is reset to zero when a gluon is radiated so that the probability of radiating another gluon starts to accumulate again with time. We note that the framework described here may be applied to other energy loss formalisms as well, as long as the energy/momentum distribution of the radiated gluons is provided.
In our simulation, the evolution of parton showers in dense QGP medium are controlled by three transport coefficientsê,ê 2 andq. These coefficients can be determined if detailed information about medium structure and jet-medium interaction are known [57]. In this work, the longitudinal drag and diffusion are related by applying the fluctuation-dissipation relationê 2 ≈ 2Tê. Unless otherwise stated, we assume the diffusion in transverse and longitudinal directions to be equal, i.e,q ≈ 2ê 2 , We will investigate the influence of varying relative sizes of transverse and longitudinal transport coefficients on the jet modification in the last section.
In Fig. 1, we show the energy loss experienced by a quark jet propagating through a static quark gluon plasma with a temperature of 400 MeV. A quark with initial energy of 100 GeV enters into the medium at t = 0, and a shower of partons will develop due to the medium-induced gluon radiation. In this plot, the jet transport coefficient is set asq = 2ê 2 = 4Tê = 2 GeV/fm 2 . We show the energy loss experienced by the leading parton and the full jet shower (with cone size R = 0.3), with and without the inclusion of the contribution from elastic collisions.
The energy loss of the leading quark is found to be dominated by radiative energy loss; the inclusion of collisional energy loss does not have much effect on leading parton energy loss. However, the energy loss of the full jet is increased by almost a factor of 2 when the contribution from elastic collisions is taken into account. This is due to the fact that for full jet shower energy loss, one needs to take into account not only the collisional energy loss from the leading parton, but the 4 effect of elastic collisions on radiated partons as well. Such an effect accumulates with time as a result of more and more partons in the shower being developed by the medium-induced radiation. This leads to stronger length dependence for full jet energy loss than leading parton energy loss. Similar effect has been found in Ref. [51,52].

Nuclear modification of photon-tagged jets
In this section, we present the numerical results for the medium modification of jets correlated with high p T photons in Pb+Pb collisions at √ s N N = 2.76TeV at the LHC. To simulate the medium modification of jet shower in relativistic heavy-ion collisions, we need to provide the initial conditions before the jet-medium interaction. In this application, the production of the photon-jet pairs from early hard scattering processes are generated from PYTHIA8.1 [58] for p+p collisions at √ s N N = 2.76TeV. We use anti-k T algorithm to reconstruct the full jets via the FASTJET package [59], from which we also obtain the detailed content of the full jets, such as the energy and momentum for each parton in a jet. The reconstructed vacuum jets are provided as inputs to our Monte-Carlo code which simulates in-medium jet evolution and medium modification. We note that during the propagation of jets through medium, energy can flow in and out of jet cone. To take into account such effect, when studying the medium effect for jets with given size (e.g., R = 0.3), the vacuum jets with a larger jet cone size (e.g., R = 1) are used as initial conditions for simulating medium modification.
The initial photon-jet production points are sampled according to the distribution of the binary collisions obtained from a Glauber model simulation. The space-time evolution profiles (energy/entropy density, temperature and flow velocities) of the bulk QGP medium that jets interact with are provided by relativistic hydrodynamics simulation. In this work, we employ a (2+1)D viscous hydrodynamics model (VISH2+1) developed by The Ohio State University group [60,61,62], with twocomponent Glauber model for hydrodynamics initial conditions. The code version and parameter tunings for Pb+Pb collisions at LHC energies are taken as in Ref. [62]. When the local temperature of the medium drops below the transition temperature of 160MeV, jets are decoupled from the medium. After jets exit the hydrodynamic medium, detailed information of the modified jets is again provided to the FASTJET package to perform the full jet reconstruction (to be compared with these jets without medium modification). In this work, we do not perform the hadronization process for the parton showers which we postpone to a future effort. To study the medium modification on the photontagged jets, we may define two distribution functions which differ only by normalization, where x T is the tagged jet momentum fraction x T = p T,J /p T,γ . Given the kinematic cuts for the trigger photons and associated jets, the former distribution is normalized to the fraction of photons having the associated jet pair R Jγ and the latter gives unity when integrating out the momentum fraction x T . The nuclear modification factor I AA for the tagged jet distribution is usually defined as In Fig. 2 we show the x T -integrated nuclear modification factor I AA = R AA Jγ /R pp Jγ (left) and the average fractional momentum loss ∆x T = x T | pp − x T | AA (right) as a function of centrality (the participant number N part ). The data points shown here are obtained from CMS measurements of R Jγ and x T for Pb+Pb collisions and PYTHIA+HYDJET references [45]. For all the presented results in this work, we take the CMS kinematic cuts, i.e., the photon momentum p T,γ > 60GeV, jet momentum p T,J > 30GeV, and the azimuthal angle between trigger photons and associated jets ∆φ = |φ γ − φ J | < 7π/8.
Here we scale the transport coefficients to the medium local temperature according to their dimensionsq = 2ê 2 = 4Tê ∝ T 3 , with the overall normalization to be fitted to one data point. In next section, we will investigate the influence of relative sizes and different parametrization forms of jet transport coefficients on medium modification of photon-tagged jets. By fitting to the value of I AA ≈ 0.7 in most central (0 − 10%)   Pb+Pb collisions; we obtain the gluon transport coefficientq A = 6.5GeV 2 /fm at the hydrodynamics initial time τ 0 = 0.6fm/c at the LHC (which corresponds to gluonq A = 2.9GeV 2 /fm at the same initial time at RHIC). Comparable values of jet transport coefficients have been obtained from other model calculations [18,63]. Due to stronger medium effect from peripheral to central collisions, we observe the decrease of the nuclear modification factor I AA and the increase of the average fractional momentum loss ∆x T for photon-triggered jets.
The distribution P (x T ) of the momentum imbalance variable x J between triggered photons and the away-side jets is shown in Fig. 3. One observes a shift of the distribution to smaller x J values for the distribution in Pb+Pb collisions compared to p+p collisions. This indicates that the away-side tagged jets lose some The azimuthal angle between trigger photons and associated jets ∆φ = |φ γ − φ J | is another interesting quantity. In Fig. 4, the normalized distribution P (∆φ) is shown for both p+p collisions and most central (0 − 10%) Pb+Pb collisions. While there is strong medium effect for the momentum imbalance x J distribution as shown in Fig. 3, no strong modification is observed for photon-jet azimuthal angle ∆φ distribution. This is partly due to the kinematic cuts applied on the selection of photon-jet events. Strongly deflected jets are usually those having experienced larger energy loss in the medium, and many of them may not be selected in the final presented results after applying the kinematic cuts. Such effect may be seen from the dashed curve which represents the results for events with no p T cut applied to the associated jet momentum after the inmedium evolution. One also observes some broadening for photon-jet azimuthal angle ∆φ distribution for this case.
We have presented above the medium modification of the photon-triggered jet yield, the photon-jet momentum imbalance, and their azimuthal angle distri- bution for Pb+Pb collisions at the LHC. We may further study photon-tagged jets by classifying the events according to their final x T values. In such way, one may use photon-tagged jets to probe different regions of the hot and dense QGP created in the collisions. This can be clearly seen in Fig. 5, where we show the density distribution of initial photon-jet production points (x ini , y ini ) in the transverse plane. Here we take those events with triggered photons that propagate along outof-plane directions (|φ γ − π/2| < π/12). We compare the initial photon-jet production point distribution for different x T values: x T = [0.5, 0.6] for (a) and (c); x T = [0.9, 1] for (b) and (d). Results for two different centralities are also shown for comparison: 0 − 10% for (a) and (b); 20 − 30% for (c) and (d). If one takes photon-jet events for all values of x T and propagating directions, the density distribution of initial photon-jet production points would follow the binary collision distribution (if no kinematic cuts are applied after medium evolution). From the plots we may see that jets with larger values of final x T typically have traversed shorter medium length than smaller x T jets. The above effect can be made more quantitative. In Fig. 6, we show the average for the photon-jet initial production points y ini if one looks at the triggered photons that propagate along out-of-plane directions (|φ γ − π/2| < π/12). One may see that on average jets with smaller x T values may have traversed about a few fm longer distance than larger x T jets. Such medium length difference traversed by different x T jets is found larger in more central collisions (the left panel) and for larger jet cone sizes (the right panel).
Since jets with different x T values have traversed different medium lengths and density profiles (and lost different amount of energy), they are expected to ex-  hibit different medium modification patterns. In Fig. 7, the photon-triggered nuclear modification factor I AA is plotted as a function of the momentum fraction x T for Pb+Pb collisions at the LHC. One may observe that the yield for the associated jets with large x T values is suppressed in Pb+Pb collisions due to the interaction with the QGP medium, while the yield for lower x T jets is enhanced due to the fact that the momentum fraction variable x T distribution shifts from higher x T to lower x T values. Such a medium modification effect is found to increase from peripheral collisions to central collisions.
The centrality dependence can be more clearly seen in Fig. 8   momentum spectra, we observe larger suppression and stronger centrality dependence for larger values of x T , and smaller suppression (or enhancement) for smaller x T jets.
The medium modification of photon-tagged jets depends on the cone sizes as well. In Fig. 9, the photontagged nuclear modification factor I AA is plotted as a function of x T for most central (0 − 10%) Pb+Pb collisions at the LHC. Results for different jet sizes (R = 0.2, 0.3, 0.4, 0.5) are shown for comparison. We observe that the suppression of the associated jet yield at high x T is similar for different jet sizes, whereas the enhancement of small x T jet yield is larger when jet cone size is increased. This indicates that the medium modification of smaller x T jets are more sensitive to the change of jet cone size due to traversing longer medium lengths.
We may see such cone size dependence due to medium path lengths in Fig. 10, where we show the x T -integrated nuclear modification factor I AA as a function of jet cone size for most central (0-10%) Pb+Pb collisions at the LHC. Results for different values of x T bins are shown for comparison. The nuclear modification of larger x T jets show weak jet cone size dependence, while a stronger dependence is observed for smaller x T jets. It is also of great interest to study the nuclear modification of jet shape observables. The jet shape ρ(r) in the transverse direction is defined as where p T,J is the transverse momentum of the reconstructed jet. In the above equation, p T,i represents the transverse momentum of the individual partons within the jet, r i = (φ i − φ J ) 2 + (η i − η J ) 2 , and ∆r is the bin size. The jet shapes for single inclusive jets have been recently measured at the LHC by CMS Collaboration [64]. A redistribution of the jet energy inside the cone has been observed in Pb+Pb collisions as compared to p+p collisions: an excess of energy fraction at large r, a depletion at intermediate r and little modification at very small r.
Here we study the jet shape observables for jets tagged by direct photons. In Fig. 11, we show the nuclear modification of jet shapes ρ(r) of photon-tagged jets for different values of x T . While some broadening effect (i.e., an excess at large r) is observed for small x T jets, for the large x T jets, jet shape ρ(r) with large r values might be suppressed. This may be understood from the different lengths traversed by different x T jets. As jets propagate through the medium, the gluon radiation (vacuum and medium-induced) at larger angles (larger r) interact with medium longer due to smaller formation times, while at later time, jet broadening effect takes over, leading to the enhancement at larger r values for smaller x T jets. These results indicate that jet shapes are sensitive to many details of full jets and their interaction with medium constituents. It is interesting to note that Refs. [65,66,43] have argued that the jet is mostly a coherent core, and thus loses energy in medium like a single parton without modifying jet inner structure. We note that the hydrodynamic response of the medium to jet transport and the redistribution of lost energy may influence the final observed jet structure. The inclusion of this contribution would be essential before performing quantitative comparison with the jet structure measurements. This will be pursued in future effort.

Dependence on jet transport coefficients
Jet transport coefficients encode information about the nature and properties of hot and dense matter probed by high energy jets. In the previous section, we related jet transport coefficients to the medium temperature according toq ∝ T 3 and we assumed the relative sizes of the longitudinal and transverse transport coefficients aŝ q = 2ê 2 = 4Tê. In general, the relative sizes of longitudinal and transverse transport coefficients may depend on detailed structure of the medium, such as the mass of medium constituents [67]. Jet transport coefficients may depend on the properties of the probes as well, such as the energy/momentum of the propagating jet showers. In this section, we investigate how the medium modification of photon-tagged jets depends on the relative sizes and different parameterizations of transverse and longitudinal jet transport coefficients.
We first look at the dependence of photon-tagged jet nuclear modification of the relative sizes between transverse and longitudinal jet transport coefficients. Here the transport coefficients are still related to the medium temperature according toq ∝ T 3 . In Fig. 12, we show the x T -integrated photon-tagged jet I AA as a function of the jet cone size for most central (0-10%) Pb+Pb collisions at the LHC. The results for different values of α ≡ 4Tê/q: α = 2, 1, 0.5 are shown for comparison. To see the pure effect on jet cone size dependence more clearly, we have refitted photon-tagged jet I AA to the same values for R = 0.2 (the smallest jet size that experiments use) when changing the relative sizes ofê andq. With the increase of the relative size ofq, a stronger jet cone size dependence is observed for the nuclear modification I AA of the photon-tagged jet yield. This is due to the fact that the relative contribution from medium-induced radiation and jet broadening effect become larger when increasing the relative size of q (the decrease of α = 4Tê/q).
In Fig. 13, we show the nuclear modification of the jet shape observable ρ(r) for photon-triggered jets as a function of jet cone size for most central (0-10%) Pb+Pb collisions at the LHC. A larger broadening effect is seen as we increase the relative size ofq. We also observe some suppression at larger values of r whenq becomes small (ê becomes large). For such case, the energy loss of jet showers is more dominated by the collisional energy loss; the radiated gluons at larger angles interact with medium earlier due to their smaller formation times.
We now investigate the dependence of photon-tagged jet modification on the different parameterizations of jet transport coefficients. For such purpose, we fix the relative sizes between longitudinal and transverse transport coefficients asq = 4Tê. We compare three different parameterizations:q ∝ T 3 , T 3 ln(E/T ) and T 5/2 E 1/2 to study the influence of energy dependence of jet trans- port coefficients. The energy dependence of the form T 3 ln(E/T ) originates from the leading-log Hard-Thermal-Loop results [16], and the use of power dependence on energy is motivated from a lattice calculation of jet quenching parameter [57]. In Fig. 14, we show the x T -integrated photon-triggered jet I AA as a function of jet cone size for most central (0-10%) Pb+Pb collisions at the LHC. Again, we have refitted the value of I AA to the same values for R = 0.2 when different parameterizations are used. One can see that the power law (E 1/2 ) parametrization gives the strongest jet cone size dependence for photon-tagged jet I AA . This may be understood from the fact that the radiated gluons with smaller energies are more difficult to interact with the medium constituents than much higher energy leading parton, thus the energy loss of jet shower is more dominated by the medium-induced radiation as compared to the other two scenarios. In Fig. 15, we show the nuclear modification of jet shape observable ρ(r) for photon-triggered jets as a function of jet cone size for most central (0-10%) Pb+Pb collisions at the LHC. Again a larger jet broadening effect is observed as we increase the energy dependence for jet transport coefficients.

Summary
We have studied the medium modification of jets correlated with high p T photons in Pb+Pb collisions at the LHC. A transport and perturbative QCD hybrid model is developed for simulating jet shower evolution and modification in hot and dense QCD medium with the inclusion of the contributions from both induced gluon radiation and elastic collisions. The effect of elastic collisions experienced by jet shower partons is encoded by a few transport coefficients, while medium-induced gluon radiation is simulated with the use of the single gluon emission spectrum obtained from higher twist jet energy loss calculations.
With our event-by-event parton shower simulation, we have presented the numerical results for the nuclear modification of the photon-tagged jet production yield, the photon-jet momentum imbalance and their azimuthal angle distribution. We have studied photontagged jets with different values of the momentum fraction x T , and they exhibit different medium modification patterns due to traversing different medium lengths and density profiles. This makes photon-triggered jets very useful for tomographic study of jet modification in dense matter created in relativistic nuclear collisions.
We have further investigated the influence of transverse and longitudinal jet transport coefficients on the modification of photon-tagged jets and found that both the relative sizes and different parametrization forms of e andq may affect the nuclear modification patterns of photon-tagged jets, especially the jet cone size dependence and jet shape observables. Such results indicate that the study of cone size dependence and jet shape observables may provide strong constraints on the quantitative extraction of transverse and longitudinal jet transport coefficients. The determination of these coefficients will provide insights into the internal structure and various transport properties of the hot and dense medium produced in relativistic heavy-ion collisions.
The above analysis may be applied for studying nuclear modification of high p T dijets and single inclusive jets. The systematic comparison to various experimental observables such as those done in Ref. [68] will provide more physics on jet modification. The inclusion of a hadronization process for final partonic jet showers will enable us to directly compare to the experimental measurements of the longitudinal momentum (fraction) distribution of hadron fragments of jets. The energy and momentum deposition profiles obtained from our simulation of in-medium jet shower evolution may be applied in hydrodynamics simulation to investigate the response of the medium to jet transport. Incorporating both jet energy loss and medium response components is essential to perform more quantitative comparison with final observed jet structure and correlations. These studies will be very helpful for more complete understanding of jet-medium interaction and are left for future effort.