Multi-partonic medium induced cascades in expanding media

Going beyond the simplified gluonic cascades, we introduce both gluon and quark degrees of freedom for partonic cascades inside the medium. We then solve the set of coupled evolution equations numerically with splitting kernels calculated for static, exponential, and Bjorken expanding media to arrive at medium-modified parton spectra for quark and gluon initiated jets. Using these, we calculate the inclusive jet RAA\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_\mathrm {AA}$$\end{document} where the phenomenologically driven combinations of quark and gluon jet fractions are included. Then, the rapidity dependence of the jet RAA\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_\mathrm {AA}$$\end{document} is examined. We also study the path-length dependence of jet quenching for different types of expanding media by calculating the jet v2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v_2$$\end{document}. Additionally, we study the sensitivity of observables on effects from nuclear modification of parton distribution functions, vacuum-like emissions in the plasma, and the time of the onset of the quenching. All calculations are compared with recently measured data.

Not all observables have the same sensitivity to the degree of complexity implemented in the theoretical description and therefore not the same sensitivity to the underlying physics which is modeled. Consequently, it is important to study the impact of individual aspects of complex modeling on individual observables. In this paper, we focus mainly on the canonical jet spectrum in lead-lead collisions at the LHC and the related nuclear modification factor, R AA . These observables are mainly sensitive to the energy lost by the leading particles in the jet, due to the bias from the steeply falling hard cross-section of jet production [31]. To quantify in better detail the dependence on path length of jets, we also focus on the azimuthal asymmetry of jet production quantified in the second harmonic coefficient of the jet multiplicity distribution, v 2 . All these observables have been analyzed in various frameworks before. However, in this paper we present a novel implementation of the effects of expanding medium which allows to quantify the sensitivity to the plasma evolution (for related work, see also [32]), and we present the calculation of the inclusive jet suppression evaluated differentially in the full set of basic kinematic quantities: transverse momentum, rapidity, and the azimuth.
In a recent work [33], we have presented a study of the distribution of medium-induced gluons and the jet calculated using the evaluation of in-medium evolution with splitting kernels derived from the gluon emission spectra for expanding profiles. Scaling behavior of splitting kernels was derived for low-x and high-x regimes in the asymptote of large times and its impact on the resulting jet R AA was discussed. 1 For the full phase space of the radiation, the scaling of jet R AA with an effective quenching parameter was introduced.
In this paper, the main aim is to study the impact of the medium expansion on observed jet suppression within more realistic description of the jet production. The description is improved by: (1) by including both gluon and quark degrees of freedom for partonic cascades into the evolution equations; (2) by improving on the description of the initial parton spectra. This should provide more realistic calculation of nuclear modification factor and more realistic values of quenching parameter,q, to be extracted from the data. Within this description we evaluate the rapidity dependence of R AA and path-length dependence of the energy loss (via jet v 2 ) and study if measurements of inclusive jet suppression differential in rapidity and azimuthal angle are sensitive to the way how the medium expands. Furthermore, we study the sensitivity of observables on effects from nuclear modification of parton distribution functions, effects from vacuum-like emissions in the plasma, and sensitivity on the time of the onset of the quenching.
The paper is organized as follows. In Sect. 2 we summarize the single parton emission spectra and the in-medium splitting rates for possible flavor combinations at the leading order. Next, we introduce the coupled evolution equations to arrive at the medium evolved parton spectra with the in-medium splitting functions as the kernels of the evolution. In Sect. 3, we first calculate jet R AA from moments of medium evolved partonic spectra for quark and gluon initiated jets. Then, in the next four subsections we analyze the transverse momentum dependence, initial-time dependence, rapidity dependence, and path-length dependence of the jet suppression. We summarize and conclude in Sect. 4.

Coupled parton cascade evolution equations
Medium-induced radiation was traditionally considered in the context of soft gluon emissions [34][35][36][37]. Naturally, the turbulent nature of the cascade was first discussed in the context of gluon fragmentation [38,39]. Nevertheless, because of the important role played by quasi-democratic splittings, it is natural to consider the interplay of quark and gluon degrees of freedom in the cascade [40,41]. From our point of view, this is also a natural demand when considering quenching effects for a realistic initial hard spectrum of initial partons, which is dominated by gluons at low-and intermediatep T and by quarks at highp T . For consistency, one should also consider the feedback from gluon to quark-antiquark splittings.
The main quantity of interest is the single-inclusive inmedium parton distributions, defined as where x = ω/E denotes the energy fraction carried by a parton of frequency ω with respect to the initial energy E of the original parton 2 and i refers to the parton flavor. This is measured at the dimensionless "time" τ ≡ q 0 /E L, where L is the in-medium path-length. 3 Here,q 0 is the value of the jet quenching coefficientq at an initial reference time. Hence, the distribution D(x, τ ) also intrinsically depends on this time scale. For a medium with constant or exponentially decaying density the natural initial time is simply t = 0 and this dependence is trivial. However, this is not the case for power-like decaying medium profiles, which are associated with a finite starting time. This will be discussed in more detail below. Concretely, we focus on the gluon D g (x) and quark singlet distributions , where the sum runs over all active flavors, which describes the energy distribution of quarks inside a jet. The evolution equations for the inclusive in-medium parton distribution can be written as [40] Each equation involves a positive gain term and a negative loss term which describe the production of a parton with energy fraction x from the splitting of a parent parton with energy fraction x/z and the decay of a parton with energy fraction x into softer fragments with fractions zx and (1−z)x, respectively. The splitting kernels K i j quantify the rate of splittings per unit evolution time τ . They can be expressed as where dI i j /dz is the medium-induced emission of parton i from parton j, carrying a longitudinal momentum fraction z.
The main purpose of this work is to explore the evolution of the cascades for expanding media. Below, we focus on three main scenarios: Static medium: In a static medium, where L is the medium length. The splitting rate reads For the case of soft gluon emissions, i.e. 1 − z, z 1, in a static medium, the splitting function is given as, Exponentially decaying medium: For exponentially decaying media the profile of the jet quenching parameter is given byq(t) =q 0 e −t/L . The splitting rate reads, Power-law decaying medium: For power-law decaying media, the jet quenching parameter is defined aŝ with α = 1 corresponding to the Bjorken expansion [42].
The splitting rate reads in this case, where with τ 0 = q 0 /Et 0 . Due to this additional dependence on τ 0 in the splitting function, the resulting distribution naturally also depends on it, i.e. D(x, τ ) → D(x, τ, τ 0 ).
The splitting rates have been written in terms of the colorkinematical factors κ i j (z) and (unregularised) Altarelli-Parisi splitting functions P i j (z) and can be found explicitly in Appendix A.

Numerical results for medium-induced parton spectra
We use the coupled evolution Eqs. (2) and (3) to numerically calculate the evolution of the quark and gluon spectra for both quark-initiated and gluon-initiated jets as a function of energy fraction x for different evolution time τ . The rates and solution for the pure gluon cascade was discussed in [33]. In the current work, we solve the coupled evolution equations using finite difference methods. For future reference, the resulting distributions are evaluated for two sets of initial conditions, as follows: gluon-initiated jets: x). Figure 1 shows the medium evolved spectra for the gluon (D g ) and the quark (D S ) distributions separately for gluoninitiated ( Fig. 1a) and quark-initiated ( Fig. 1b) jets as a function of x evolved up to τ = 1.0. The evolution has been performed in media with different density profiles, i.e. performing the evolution with the different kernels in Eqs. (6), (7) and (9), see the figure legends for details. For the Bjorken case, we have chosen τ 0 = 0.5. One can see that gluon spectra follow the 1/ √ x turbulent behavior for both quark and gluon initiated jets and for both the static medium and the expanding medium (the static case was already discussed in [40]). However, for the quark spectra this scaling is only approximate for all the medium profiles including the static profile. This is mainly because the stationary behavior, i.e. D g (x) ∼ c 1 / √ x and D S ∼ c 2 / √ x (where c 1 and c 2 are approximately constants), has not yet had sufficient time to establish itself [40]. Note, that the speed of the evolution is slower when considering the splitting rates for finite media, as done above [33].
Outside the region x ∼ 1, the exponential spectra have a greater amplitude than the full static spectra and similar (a) (b) Fig. 1 Evolution of the quark singlet spectra (D S ) and the gluon spectra (D g ) for the "gluon initiated" (left panels) and the "quark initiated" jets (right panels) for a particular evolution time τ = 1.0 for different medium profiles amplitude as the static soft spectra. Spectra for the Bjorken expanding media have the least amplitude due to a smaller rate of splittings compared to the other profiles. We have shown previously [33] that the "singular" parts of gluon spectra (z → 0, z → 1) for different medium profiles scale with the static one when using an effective quenching parameter, defined as 4q 0 and 4q 0 t 0 /L for the exponential and Bjorken profiles, respectively. This scaling between different medium profiles holds for individual contributions discussed here as well. However, such a scaling does not work when employing the full evolution kernels as obtained in Ref. [33] for gluon-only cascades.
In general, the D S and D g dominate the high-x part of the spectra for quark-initiated and gluon-initiated jets, respectively. For gluon-initiated jets, the probability of finding a gluon in the parton cascade at high as well as small x is significantly larger than a quark for all the cases of expanding media. For quark-initiated jets, the probability to find a gluon is smaller than a quark at high x and larger at small x. This is a generalization of similar observations done in Ref. [40] for static media.

Inclusive jet suppression
We intend to advance understanding the inclusive jet suppression in several directions in this work. Baseline of this work is to quantify the impact of medium expansion on p T (Sect. 3.1), rapidity (Sect. 3.3), and azimuthal dependence (Sect. 3.4) of the inclusive jet suppression in a realistic treatment of input parton spectra. Then we quantify the sensitivity of observables on the time of the onset of the quenching (Sect. 3.2).
Beside this, we also quantify the impact of including nuclear parton distribution functions (nPDF), impact of including effects from early vacuum-like emissions in the plasma, and sensitivity on input parton spectra connected with the choice of free parameters in underlying MC modeling. Although we consider scenarios with expanding media, we assume a fixed in-medium path length, L = 5 fm, for all results except those presented in Sect. 3.4.

Transverse momentum dependence of jet suppression
To quantify the suppression factor of jets in heavy-ion collisions one usually assumes a separation between the shortdistance production of hard partons and their subsequent propagation through the medium. We define medium evolved spectra of quark and gluon cascades for quark-initiated jets (denoted by "{q}") and gluon-initiated jets (denoted by "{g}") as One of the key observables quantifying inclusive jet suppression is the jet nuclear modification factor, measured by the LHC experiments [11,12,[43][44][45]. In our approach, the yield for the inclusive jet suppression can be obtained as a convolution of the D i (x, τ ) distribution with the initial parton spectra [31,33,46,47]. For a given parton species i that is produced in the hard cross section dσ 0 i /(d p T dy), we write for i = q, g and where the evolution time τ = q 0 / p T L is written in terms of the final, measured transverse momentum p T . Here, dσ 0 /d p T denotes the hard spectrum for producing a parton i in proton-proton collisions, or in "vacuum". For the Bjorken case, where the distribution has an explicit dependence on the initial time τ 0 , we approximate since τ 0 is already very small and the integral (14) is dominated by x ∼ 1. The quenching factor for parton species i, Q i ( p T ) , defined as the ratio of partonic medium and vacuum spectra, is then measured at the same final p T and the same rapidity y (or rapidity range). If we approximate the initial parton spectrum by a power law dσ 0 /d p T ∝ p −n T , with fixed power index n =const., the quenching factor is then a (n − 1) st moment of the in-medium parton distribution, i.e.
In this study, we go beyond this simple approximation and extend the power law into a more general form where which provides a good characterization of the initial parton spectra [48]. Here β i , γ i , and δ i are additional free parameters to be fitted separately for quark-and gluon-initiated jets. The rapidity dependence is parameterized via the coefficients n i,0 , β i , γ i and δ i that are extracted for each rapidity selection. In this case, the quenching factor for a given parton species i, is written as This factor describes the effect that energy loss has on the spectrum of a single parton. However, a jet is a multi-parton state due to the large phase space for radiation, and at highp T many of these emissions will take place with momentum scales that are much higher than the typical medium scales attainable via multiple scattering in the medium [49]. Every emission of this type corresponds to a resolved color charge that contributes toward the total energy loss of the full jet.
The resummation of such emissions in an expanding medium is derived in Appendix B, see also [32] for a related discussion. Implementation of these vacuum-like emissions (VLE) can be turn on and off allowing to assess their impact on measurable quantities as discussed further in this section. The full nuclear modification factor contains the contributions from both quarks and gluons to the final jet spectra, and reads Finally, the combined nuclear modification factor R AA reads where Above, we have explicitly included the dependence of jet cone-size R on the generated spectra and the quenching factors. Here σ 0,med i is the vacuum cross-section for the hard production which takes into account nPDF effects. If these are to be neglected, this factor is replaced by σ 0 i .  Parameters of the extended power law description of input parton spectra were obtained from fits of jet spectra from PYTHIA8 [50]. In general, three versions of spectra were used to study the sensitivity of observables on modeling the input parton spectra: (1) spectra based on PYTHIA8.185 with AU2 tune [51] and CT10 PDFs [52] ( configuration labeled "vacuum 1"), (2) spectra based on default version of PYTHIA 8.306 (configuration labeled "vacuum 2"), and (3) spectra based on default version of PYTHIA 8.306 with nPDFs from EPS09LO [53] (configuration labeled "nPDF"). The power-law parameters for the three configurations are summarized in Tables 1 and 4. Jet spectra were obtained by running FastJet [54] using the anti-k t algorithm [55] with resolution parameter R = 0.4 for the center of mass energy of √ s NN = 5.02 TeV. Jet spectra were reconstructed separately for several rapidity selections specified in Table 1 and used later in the study. Distributions calculated inclusively in rapidity were obtained for rapidity selection of |y| < 2.8. The jets were matched to one of two outgoing partons from the leading order hard-scattering process by choosing the parton with the smallest angular distance, R = φ 2 + η 2 , and the flavor of the jet was assigned to be that of the matched parton. Furthermore, for all the numerical results presented below we set α s = 0.1. The expressions for quenching factor, jet R A A , and flavor fraction in the configuration with no nPDF effects and no VLE as well as power-law parameters for simple "vacuum 1" configuration are provided in Appendix C.
Summary of calculated p T dependence of nuclear modification factor R AA is plotted in Fig. 2 along with experimental data on measured anti-k T jets with R = 0.4 [44]. For each configuration, we use value ofq 0 obtained from the χ 2 minimization of the difference in the nuclear modification factor between the data and the theory. We fixed the parameter L to 5 fm for all the medium profiles.
Left panel of Fig. 2 compares results for static medium for four configurations: two versions of input parton spectra ("vacuum 1", "vacuum 2"), input parton spectra with nPDF effects included ("nPDF"), and input parton spectra with nPDF effects and vacuum-like emissions included ("nPDF+VLE"). One can see that shape of the R AA can be fully described only after implementing nPDF effects. Comparison of "vacuum 1" and "vacuum 2" configurations shows minimal sensitivity of the shape of R AA ( p T ) on PYTHIA version of in-vacuum input parton spectra. Comparison of "nPDF" and "nPDF+VLE" configuration implies that adding vacuum-like emissions to the calculation has an important impact on both the shape of R AA and its overall normalization which has an impact on the extracted values ofq 0 .
Right panel of Fig. 2 shows R AA for four types of medium expansion in the configuration "nPDF+VLE". One can see that the difference in the shape and magnitude of R AA among different medium profiles is rather small, less than 10% for p T > 100 GeV. The obtainedq 0 values for realistic quark and gluon initial spectra in the configuration of nPDF+VLE are tabulated together with the values obtained for a gluon-only calculation with a simple power-law spectrum and no nPDF and VLE effects as published in Ref. [33] in Table 2. The extracted value of the jet quenching parameter differs significantly between the simple configuration and the more realistic treatment. The driving factor of this difference is the use of both the quark and gluon degrees of freedom and improved description of input parton spectra. This illustrates and quantifies a simple fact that the value of jet quenching parameter has an important dependence on the level of completeness of the jet energy loss description and approximations used. Since the R AA distributions do not scale when using the above quotedq 0 values, we confirm the breakdown of the "effective" scaling laws for R AA , which was discussed in Ref. [33], also for more realistic treatment presented in this paper. Fig. 2 Left: jet R AA for static medium in four configurations: two version of vacuum input parton spectra ("vacuum 1,2"), spectra with nPDF effects ("nPDF") and spectra with nPDF effects and vacuum-like emissions included ("nPDF+VLE"). Right: jet R AA for "nPDF+VLE" configuration for four medium profiles: static soft (black), static (green), exponential (blue), and Bjorken with t 0 = 0.1 fm (red) Theq 0 values were optimized for each profile to reproduce the data (in gray) [44] Table 2 The jet quenching parameter for different media profiles obtained from fitting the calculation to the data from Ref. [44]. Full description of the input parton spectra (first line) can be compared with results using gluon-only cascades and power-law description of input parton spectra (second line), published in Ref. [ The results in Fig. 2 shows that the single-parton description of the jet R AA used here describes the shape of the experimental data in the p T region below ∼ 250 GeV. On the contrary, at p T 250 GeV the calculated R AA is steeper compared to the data. This was already observed in Ref. [33] for the gluon-only calculation and the more complete description of jet energy loss used here does not lead to an improvement in the description of the data. Comparison of "nPDF+VLE" configuration with the data suggests that the flatting of R AA at highp T is due to the combination of nuclear PDFs, that produce a suppression at highp T due to the EMC effect (see also [56,57]), and the onset of the quenching of multiple resolved partons in a jet [58]. These appear due to hard vacuum-like splittings inside the medium [49], see also [59]. A resummation of such fluctuations leads to stronger suppression at highp T due to enhanced energy loss as discussed in Refs. [60][61][62][63].
The inability to describe the full p T region of measured R AA in the configuration without the nPDF+VLE effects is in contrast with parametric modeling presented in Refs. [44,48], where the use of improved description of input parton spectra and use of appropriate quark and gluon mixing was shown to lead to a good description of the R AA over a large p T range under the assumption of p T -dependence of the energy loss which was close to √ p T -dependence. Within the BDMPS, the √ p T -dependence of the lost energy holds at lowp T , but does not hold at highp T [31], where the singlescattering contribution to the medium-induced energy loss dominates in calculations [64], see also [62,63], and where nPDF effects play a role. The power of 0.55 used in parametric modeling of Refs. [48] may effectively account for these phenomena, however, without providing ability to distinguish them.

Initial-time dependence of jet suppression
The value of t 0 in the Bjorken medium quantifies the time during which no interaction with the medium occurs. This reflects a situation with a finite time needed for the onset of the quenching interactions in the medium. In this section, we discuss the impact of the choice of t 0 on the observable jet suppression by comparing two different starting times of the quenching in the Bjorken medium, namely t 0 = 0.1 fm and t 0 = 1.0 fm, in the most simple configuration with no nPDF effects and vacuum-like emissions included. Left panel of Fig. 3 shows a comparison of the q for the Bjorken media with these two choices of t 0 evaluated as a function of L calculated as for Bjorken medium with L = 5 fm and withq 0 values optimized to reproduce the experimental data [44]. Note that we have not included nPDF effects in this particular figure whereq(t) =q 0 t 0 /t for the Bjorken scenario. One can see that delaying the onset of quenching leads to a larger average value ofq at L > t 0 compared to an earlier time of the quenching. In general, the later the quenching sets in, the larger the q is, which makes the medium more opaque at the later stages of the evolution. We use again the χ 2 minimization procedure to find thê q 0 value which minimizes the difference in the R AA between the data and the theory prediction for the two t 0 values with medium length fixed to L = 5 fm. The results are shown in the right panel of Fig. 3 (note that we have left out nPDF effects in this particular figure). For t 0 = 0.1 and 1.0 fm, thê q 0 is found to be approximately 1.8 and 1.0 GeV 3 , respectively. That is for ten times larger initial time, theq 0 has to be roughly two times smaller to obtain the same R AA . While the choice of t 0 can be somehow absorbed into an adjusted value ofq 0 , we will see below that it affects other observables. In the next two sections, we use optimizedq 0 values, and their respective t 0 parameters, to calculate more differential observables derived from R AA which characterize the inclusive jet suppression and quantify how these settings affect them.

Rapidity dependence of jet suppression
Following the measurement done in Ref. [44], we study the rapidity dependence of the jet suppression. The optimized q 0 values from Sects. 3.1 and 3.2 are used to calculate the rapidity dependent R AA . The results are presented in Fig. 4 in terms of the ratio of the R AA in a given rapidity region to the R AA in central rapidity region (|y| < 0.3). Figure 4a shows the ratio of R AA in 2.1 < |y| < 2.8 and |y| < 0.3 as a function of p T for the static medium for four configurations: two different versions of input parton spectra ("vacuum 1", "vacuum 2"), input parton spectra with nPDF effects included ("nPDF"), and both nPDF and VLE effects included ("nPDF+VLE"). One can see that the resulting ratio differs significantly among different configurations.
Left panel of Fig. 4b detail the rapidity dependence of R AA ratio for a fixed p T , namely 316 GeV < p T < 562 GeV. One can see that inclusion of nPDF and VLE effects leads to qualitative agreement between trends seen in the data and calculations. At the same time a trend of decreasing R AA ratio is partially included in the second version of PYTHIA used in this study as well. Here, the differences between the different PYTHIA spectra should be taken as an indication of an overall theoretical uncertainty of the rapidity behavior.
Right panel of Fig. 4b shows the same quantity as the middle panel evaluated for default version of input parton spectra configuration for five medium profiles. Results for static, exponential, and two Bjorken profiles discussed in Sect. 3.2 are found to differ by less then 10% for all the rapidity bins (and also for all the p T values which is not shown in the figure).
The fact that the rapidity dependence of the inclusive jet suppression is very similar for different medium profiles implies that the origin of the rapidity dependence is somewhat universal. To check this, we have calculated the ratio of R AA for gluon-only evolution. When evaluated as a function of |y| we see the same trends as in the default case that includes both quark and gluon fragmentation. This implies that the rapidity dependence is in part a consequence of the change in the steepness of input parton spectra. However, both the nPDF and VLE effects have clearly larger impact on this observable quantity. Black, blue, green, and red color encodes static medium for "vacuum 1", "vacuum 2", static medium with nPDF effects, and with nPDF+VLE effects, respectively, in the left and middle panels. Black, green, blue, red and pink color encodes static soft, static, exponential, and two Bjorken medium profiles, respectively, in the right panel. The ATLAS data taken from Ref. [44] 3.4 Path-length dependence of jet suppression and jet v 2 The path-length dependence of jet suppression can be studied using the jet v 2 evaluated in a simple approximation [65] as where L in and L out are path lengths in the direction of the event plane and in the direction perpendicular to the event plane, respectively. Values of L in and L out were obtained from a hard sphere model where vertices of hard processes and directions of back-to-back initial partons in the azimuthal plane were randomly generated in the overlapping region of two spheres. The mapping between the impact parameter and N coll values is done using a Glauber model [66]. It should be stressed that this approach represents a tool for basic investigation of the relation between the medium expansion and the jet v 2 , rather than a modeling with a full description of the problem. The values L in and L out for different centralities are listed in Table 3. As the impact of nPDFs and vacuumlike emissions has no intrinsic azimuthal dependence and it cancels in the ratio (22), we show the results only for the default input parton configuration. The upper left panel of Fig. 5 shows the jet v 2 evaluated as a function of p T for central collisions (N coll = 1580) along with the result of measurements from Refs. [67,68]. One can see that the magnitude and the trends of the jet v 2 distribution obtained from the calculation match those seen in the data. The upper right panel of Fig. 5 shows the jet v 2 evaluated as a function of N coll . Both p T and N coll dependence of jet v 2 show rather small differences among static, static soft, and exponential medium. The results for Bjorken medium with t 0 = 0.1 fm differ significantly from results for other medium profiles. Even larger difference can be seen for the Bjorken medium with two different choices of t 0 shown in lower panels of Fig. 5. These observations have similar implications as those discussed in the case of the rapidity dependence of R AA : the impact of the medium expansion can be largely scaled out by a suitable choice ofq 0 , which is however not the case for the starting time of the quenching; the jet v 2 remains sensitive to choice of t 0 . This is in agreement with findings of the sensitivity of v 2 on t 0 published in Ref. [69] which was done in more complex modeling of the collision geometry, but less complex modeling of the medium induced showering than in this paper.

Summary and conclusions
Both quark and gluon degrees of freedom were used to calculate the evolution of partonic cascades in expanding media (static, exponential, and Bjorken). This together with an improved description of input parton spectra makes the modeling of medium induced radiation in the limit of mul- tiple soft scattering more realistic compared to gluon-only cascades and power-law modeling of input parton spectra. Using numerical estimates of medium-induced parton spectra, the nuclear modification factor R AA of jets at a fixed cone-size R = 0.4 was calculated and jet quenching parameter,q, was extracted by comparing the calculations with the data. The differences between extractedq in the improved description and in the simple, gluon-only scenario [33], were found to be significant (e.g. for Bjorken medium with t 0 = 0.1 fm theq 0 value is by 30% smaller in the improved description compared the simple gluon-only case). At the same time, the shape of the R AA is not modified significantly in the improved description compared to the simple scenario. Only when including effects from nuclear parton distribution functions and vacuum-like emissions, the shape of the calculated R AA starts matching the measured p T dependence observed in the data.
The rapidity dependence of the jet suppression is studied and it is found that the decrease of the R AA seen in the forward region by ATLAS [44] is consequence of multiple effects, namely the change in the steepness of input parton spectra and impact of effects from nuclear parton distribution functions and vacuum-like emissions. When comparing results for different medium profiles we conclude that the rapidity dependence of the inclusive jet suppression is not very sensitive to the way how the medium expands. Further, the jet v 2 is studied in an approximation of hard-sphere geometry. It is found that the differences in the jet v 2 between full calculation of static medium, static soft, and exponential medium are within 15%. On the contrary, jet v 2 differs in maximum roughly by a factor of two between Bjorken scenario with t 0 = 0.1 fm and t 0 = 1.0 fm. These findings indicate that the simultaneous study of jet v 2 and rapidity dependence of the jet production may help disentangling various effects entering the observed jet quenching.
The presented study, which quantifies the inclusive jet suppression differentially in the full set of basic kinematic quantities-transverse momentum, rapidity, and the azimuth, represents one more step towards realistic modeling of the parton energy loss. It should contribute to the effort of disentangling the impact of different effects on measured quantities characterizing inclusive jet suppression.
Acknowledgements We would like to thank the referee of our paper, whose comments and questions helped improve the presentation of our results. KT

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The experimental data used in this publication has already been published by experimental collaborations cited in the publication. Data related to calculations presented in the publication can be reproduced by following the procedures detailed in the text and references therein.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

Appendix A: Kinematical factors in splitting rates
The kinematical factors κ i j (z), appearing in the splitting rates, read Finally, the un-regularized Altarelli-Parisi splitting functions at leading order are given by

Appendix B: Multi-parton quenching in expanding medium
Here we provide details about the calculation of the resummed jet quenching factor in an expanding medium. First, we have to calculate the phase space for vacuum-like emissions (VLEs) inside the quark-gluon plasma. The relevant time scale for decoherence can be read off the decoherence parameter, or the correlator of a dipole at a fixed opening angle, 4 which corresponds to the probability for the dipole to remain in a color-correlated state. For our discussion, the dipole cross section v(r) is given simply by the harmonic oscillator approximation, i.e.
v(r, s) = 1 4q (s)r 2 , (B.10) and the dipole follows the trajectory r(s) θ 0 s, where θ 0 is the opening angle. We can define a decoherence time by considering the argument of the exponential in (B.9), namely where the index α depends on the density profile of the medium and follows from the resulting length dependence. As we will see below, α = 3 for static and exponential profiles, while it is α = 3 − γ for media whose density decays as t −γ . Processes at time scales shorter than t d are not resolved by the medium. Therefore the phase space for VLEs is given by the set of conditions [49,59] t f t d L (B.12) at DLA accuracy. At this accuracy, the available phase space for VLEs can be written as where the formation time is t f = 2/(ωθ 2 ), C i is the color charge of the parent parton and we treat the running coupling as a constant in this approximation. The quenching of a full jet is described by the collimator function [49,62]. In the approximation of small energy losses, we use the linearized approximation to obtain where i = q, g. We now proceed to work out in detail the scales of the phase space that allows for VLE's inside the medium for different medium profiles.

Static profile & exponential profile :
For these cases there is only a difference in a numerical constant.
The decoherence time is is the decoherence time [32]. From the condition t d = L, we also deduce the characteristic decoherence angle The well known results for static medium, γ = 0, are straightforwardly reproduced. Furthermore, γ = 1 corresponds to the Bjorken expansion scenario. The set of conditions for VLEs [49,59], t f t d L, demands that In the Bjorken case, the first condition leads to a restriction on the transverse momentum of splittings, namely combined with a restriction on the angle, given by θ > θ c = 8/(q 0 t 0 L 2 ). Emissions that fail to meet these conditions are resolved in the medium, and should be treated as mediuminduced processes. Finally, we also have to demand that the splittings are perturbative, i.e. k ⊥ > Q 0 ≈ 1 GeV.
The available phase space for VLEs at DLA accuracy can be written as where C i is the color charge of the parent parton and we treat the running coupling as a constant in this approximation. Note, that for lowp T jets, i.e. when p T ≤ 1 4q (L)L 2 , the requirement on the minimal angle becomes θ > θ d , where where Q s (L) 2 =q(L)L/2. Changing variables in (B.26), we finally obtain in = 2 where R min = max[θ c , θ d ] and k ⊥,min = max[Q s (L), Q 0 ]. This phase space counts the number of VLE, created inside the medium, that will contribute to the energy loss of a jet.

Appendix C: Nuclear modification factor with no nPDF effects
In this section we summarize analytic formulae for nuclear modification factor for the case of not including effects from nuclear modifications of parton distribution functions. In that