Novel subjet observables for jet quenching in heavy-ion collisions

# Novel subjet observables for jet quenching in heavy-ion collisions

Liliana Apolinário LIP, Av. Prof. Gama Pinto, 2, 1649-003 Lisbon, Portugal IST University of Lisbon, Av. Rovisco Pais 1, 1049-001, Lisbon, Portugal liliana@lip.pt    José Guilherme Milhano LIP, Av. Prof. Gama Pinto, 2, 1649-003 Lisbon, Portugal IST University of Lisbon, Av. Rovisco Pais 1, 1049-001, Lisbon, Portugal Theoretical Physics Department, CERN, CH-1211, Geneva 23, Switzerland guilherme.milhano@ist.utl.pt    Mateusz Ploskon Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA mploskon@lbl.gov    Xiaoming Zhang Institute of Particle Physics, Central China Normal University, Wuhan 430079, Chinaxiaoming.zhang@cern.ch
Received: date / Revised version: date
###### Abstract

Using a novel observable that relies on the momentum difference of the two most energetic subjets within a jet  we study the internal structure of high-energy jets simulated by several Monte Carlo event generators that implement the partonic energy-loss in a dense partonic medium. Based on inclusive jet and dijet production we demonstrate that  is an effective tool to discriminate between different models of jet modifications over a broad kinematic range. The new quantity, while preserving the collinear and infrared safety of modern jet algorithms, it is experimentally attractive because of its inherent resilience against backgrounds of heavy-ion collisions.

###### pacs:
12.38.MhQuark-gluon plasma and 13.87.-aJets in large-Q2 scattering and 24.10.LxMonte Carlo simulations and 25.75.-qRelativistic heavy-ion collisions

## 1 Introduction

Interactions of high-energy partons with a strongly coupled hot partonic medium - a quark-gluon plasma (QGP) Satz:2000bn (); Bass:1998vz (); Shuryak:1984nq (); Cleymans:1985wb () - created in heavy-ion collisions, leading to modifications of the internal jet structure (jet quenching), was first proposed in Bjorken:1982tu () and is studied as a sensitive probe of the medium properties Wiedemann:2009sh (); Burke:2013yra (); Mehtar-Tani:2013pia (). Experiments at RHIC and the LHC observed a strong suppression of high transverse momentum particle yields Adams:2005dq (); Adcox:2004mh (); Arsene:2004fa (); Back:2004je (); Aamodt:2010jd (); Aamodt:2011vg (); CMS:2012aa (), suppression of inclusive and semi-inclusive yields of fully reconstructed jets Aad:2010bu (); Chatrchyan:2011sx (); Adam:2015ewa (); Adamczyk:2017yhe (); Adam:2015doa (), and, more recently, the internal structure of the jets Aad:2014wha (); Chatrchyan:2013kwa (); Chatrchyan:2012gw (); Sirunyan:2017bsd (); Cunqueiro:2015dmx () for detailed studies of jet quenching. However, in all these measurements the treatment of the background originating from the copiously produced particles not associated to hard scatterings poses an experimental challenge for precise and unbiased measurements. Previous works that addressed the effect of filtering on subjet analysis Rubin:2010fc () and the recent analytic calculations of the momentum distributions of subjets Kang:2017mda (); Larkoski:2015lea (); Larkoski:2017bvj () and groomed jet mass distributions Frye:2016aiz (); Marzani:2017mva (); Marzani:2017kqd () that were recently measured in proton-proton collisions CMS:2017tdn (); Aaboud:2017qwh (), provide a strong motivation for novel studies. In this writeup, following previous works Zhang:2015trf (), we propose observables that are sensitive to the internal jet structure but significantly alleviate the difficulties associated to the effects of the background. Our approach is attractive from the experimental point of view and obeys the theoretical requirements of the infrared and collinear safety.

While jet substructure techniques are extensively used in the high-energy collisions Chatrchyan:2013vbb (); ATLAS:2012am (); Tripathee:2017ybi () the first studies in the heavy-ion context Sirunyan:2017bsd (); Kauder:2017cvz () are recent (see Zhang:2015trf () for a first attempt of using subjets as a phenomenological tool for jet quenching studies). In this manuscript we propose an observable that uses only the highest and next-to-highest energetic fully reconstructed subjets within a jet. This choice aims to minimize the impact of the heavy-ion background on the extracted jet properties allowing for better experimental control.

## 2 Observable definition and setup

We introduce a new jet substructure observable defined as the difference between the fractions of transverse momentum of a jet carried by its leading (hardest) and subleading (second hardest) subjets. That is

 ΔS12=z1−z2, (1)

where

 zi=pT,i/pT,jet. (2)

The subjets used to evaluate eq. (1) are obtained as follows:

1. For each event, reconstruct jets with the anti- algorithm Cacciari:2008gp () provided by the FastJet package Cacciari:2011ma () with radius and within pseudo-rapidity ;

2. Within each jet, find subjets by reclustering the jet components with a smaller radius parameter . Retain the two hardest (highest-) subjets.

The subjet samples used in this study were obtained with and . In general, the reclustering of the jet components into subjets in step (ii) above can be carried out with a different jet algorithm from that chosen to reconstruct the jet to which they belong. We chose to use anti- after assessing the discriminating power of for subjets reconstructed with different algorithms and checking its sensitivity to hadronization effects (see Subsections 4.2 and 4.4) The subjet radius parameter was set to except when assessing, in Subsection 4.2, the dependence of on  (where the range was considered), and when comparing, in Subsection 4.5, with the analysis Sirunyan:2017bsd () (where we used ).

The bulk of soft particles produced in high-energy collisions is not a priori distinguishable from the particles produced from the hadronisation of an energetic parton shower. The presence of these background particles is the main experimental confounding factor when establishing the jet energy scale and jet energy resolution in jet quenching studies (see for example Abelev:2012ej ()). Moreover, unlike in measurements of proton-proton collisions with high event pile-up probability within the detectors, the background in heavy-ion collisions is complex. It consists of region to region fluctuations, modified particle production as compared to pp collisions, and particle correlations caused by the collective expansion of the QGP. In consequence, experimental observables at relatively low jet energies at the LHC ( GeV) are prone to systematic uncertainties related to complicated multi-dimensional unfolding procedures that are susceptible to large correction factors.

The substructure observable  defined in eq. (1) has been constructed to minimize correlated background contributions. Take Cacciari:2007fd ()

 (3)

where is the true subjet momentum, is the area of a subjet, is the level of noise corresponding to the amount of transverse momentum added to each subjet per unit area by the background, and the is the experimentally reconstructed subjet momentum containing the background contribution . For subjets reconstructed with the same radius parameter  (in our case ) with the anti- algorithm, the corresponding active areas are necessarily very similar . In an ideal case, where , the background term in the numerator of  vanishes. For real events, where subjets sit close by, and can only differ by very localized fluctuations and thus should be on average still very similar. Thus, the background effect in the numerator of eq. (1) is small.

A variety of observables similar to  can be defined. In particular,  in eq. (2) could be redefined by replacing the denominator by the sum of the momenta of the leading and subleading subjets, such that . Although such a definition could have some welcome consequences in reducing the influence of background effects in the reconstructed jet transverse momentum (denominator of eq. 2), all information on the overall hardness of the jet fragmentation, that is the fraction of jet momentum carried by the two hardest subjets, would be neglected.

## 3 Models

To assess the potential of the proposed observable we consider a set of Monte Carlo event generators which rely on different implementations of jet quenching. This allows both for a comparison between theoretical calculations that is not limited by systematic uncertainties of the putative experimental measurement and to assess the potential of the observable as a discriminant of different modelling scenarios.

Below we provide a short description of each event generator considered in this study — Q-Pythia v1.0.2 Armesto:2009fj (), Jewel v2.0.0 Zapp:2013vla (), and PyQuen v1.5.1 Lokhtin:2005px () — emphasising only the main characteristics and details of the setup we adopted (for further details please see the corresponding references). All samples used in this work were generated for central (0-10% most central) PbPb collisions at TeV.

Q-Pythia is a modification of Pythia 6.4 Sjostrand:2006za () where the splitting probability in the final state parton shower is enhanced by an additional term that follows the BDMPS-Z radiation spectrum Baier:1996sk (). The medium is modelled by a single parameter, a local in space and time transport coefficient that translates the averaged transverse momentum squared exchanged between a parton and the medium per mean free path in that medium, such that . The time and spatial variation of is modelled following the PQM prescription Dainese:2004te (). We considered two different average values ( GeV fm and GeV fm) known to capture the main jet quenching features observed in dijets Apolinario:2012cg ().

Jewel implements a description of jet evolution that takes into account both elastic and inelastic energy losses as all scatterings with the medium are described by infra-red continued leading order matrix elements for processes. Additional medium-induced radiation is also taken into account during the jet development and can be induced by several coherent scatterings, as predicted by the LPM effect Landau:1953um (); Migdal:1956tc (). We kept all default settings and used the medium implementation with Bjorken expansion described in detail in Zapp:2013zya () validated on a large set of jet quenching observables Zapp:2013vla ().

PyQuen is a modification (afterburner) of standard Pythia 6.4 jet events in which both radiative and collisional accumulated energy losses are applied during the parton shower development. The former is calculated for an expanding medium within the BDMPS framework, where the angular distribution follows three simple parameterisations (small, wide and collinear angular distributions) that are used for comparison purposes. The latter is calculated in the high-momentum transfer approximation. Additional in-medium gluon radiation is added at the end of the parton shower, before hadronization. We chose the internal parameters that characterise the QGP formation expected for central PbPb collisions at the LHC.

While a typical Monte Carlo reference for jet production in pp collisions is constructed with Pythia Sjostrand:2006za () each of the models provides their own implementation and/or modifications of Pythia original routines and consequently their own pp reference. Therefore, when comparing the medium-modified jets with jets showering in vacuum we take the model-provided proton-proton collision equivalent.

## 4 Results

We provide examples of how  can be used to discriminate between the different implementations of jet quenching (Sec. 4.1) and evaluate its sensitivity to the choice of algorithm for subjet clustering and subjet radius  (Sec. 4.2). Section 4.3 illustrates how  combined with a dijet analysis can be used to study jet quenching more differentially as compared to the inclusive measurements. In section 4.4 we show the robustness of the results against hadronization effects. Finally, Sec. 4.5 provides an overview of the relation between  and the recently explored observable in vacuum.

### 4.1 ΔS12 as a model discriminant

The distribution of the difference  between the fractions of the jet total transverse momentum carried by the leading and subleading subjets is shown in Figure 1.

In vacuum – Q-Pythia (vac) and Jewel (vac) (top panel), and Pythia 6 (bottom panel) – the distribution displays a pronounced maximum for and a tail towards lower  values. Medium effects in Q-Pythia and Jewel (top panel) modify the  distribution in incompatible directions. Jewel enhances significantly the maximum of the distribution and mildly depletes its tail. Q-Pythia softens the peak at high  and produces a flat tail towards values of   with the effects more pronounced for increasing . These observations are consistent with a collimation of jets in Jewel and broadening in Q-Pythia as compared to their vacuum references. PyQuen (Coll) (bottom panel) modifies the  similarly to Jewel, PyQuen (Small) gives a distribution with features resembling those found for Q-Pythia, and PyQuen (Wide) displays an intermediate behaviour.

While the above features directly reflect the behaviour of the  and  distributions (see appendix A), we emphasise that, from the experimental point of view, studies of  are more attractive as compared to the individual distributions since the difference removes, by construction, a large fraction of the correlated background. Although the effect of the uncorrelated background is enhanced in  ratio with respect to  on a jet-by-jet basis, we found that this effect is small when taking the integrated/inclusive distributions, and subsequently their moments. Moreover,  is more robust for low momentum jets for which becomes gradually (with decreasing ) dominated by background particles. u

To further expose the differences among models we now turn our focus to the quartiles of the  distribution, considering med/vac ratios , where med refers to calculations including jet quenching effects and vac to the corresponding model specific no-quenching baseline

 RQi[ΔS12]=Qi[ΔS12]medQi[ΔS12]vac. (4)

While the med/vac ratio of the medians of the  distributions as a function of  display a clear evolution and discrimination power among the models (data not shown), we find that to characterise the modifications to the subjet structure for models that show jet collimation – Jewel and PyQuen (Coll) – the ratios of the first quartile () of  distributions is preferable. The criteria of selecting the best discriminant was made by calculating the relative standard deviation (RSD) given by the models in each bin, i.e, the ratio of the standard deviation over the mean. The relative spread among the different models for each observable is thus quantified (a larger spread translates into a larger RSD) and it can be used as a guiding parameter to select the observable that maximizes the differences among jet quenching models.

The ratios are shown in Figure 2, in the upper panel, while the corresponding RSD in the bottom panel of the same figure.

Here we find a clear evolution with the jet momentum for jets with     for all models. Up to this jet all models show a suppression of reflecting more balanced momentum sharing between the two leading subjet structures than in the vacuum references. However, at high jet , this observable remains fairly constant and shows a strong sensitivity to models that produce jets with a more symmetric structure, such as Q-Pythia and PyQuen(Small) () separating them well apart from PyQuen (Wide and Coll) and Jewel ().

Further, we find that the interquartile range , that characterises the width of the  distribution, gives additional information. Figure 3 shows as a function of jet transverse momentum for the different quenching models with the corresponding RSD calculated in each bin. Here again, models that result in jet collimation, characterised by a similar or narrower  distribution than its vacuum reference ( ), are clearly separated from those that broaden the jet, where  is typically broader with respect to the vacuum reference ().

Moreover, for Q-Pythia, Jewel and PyQuen(Coll), converges quickly to a constant value with increasing jet . Importantly, it also allows to better discriminate between the two models that destroy the vacuum subjet asymmetry: while Q-Pythia is well separated from its vacuum reference for all   , PyQuen(Small) evolves slowly towards more asymmetric jets with increasing . Thus, provides relevant complementary information to identify the main characteristics of jet quenching within specific models, in particular for jets with  , where the first quartile of the  is suppressed with respect to the vacuum reference in all models. Comparing the RSD of and we find that for low jets (jets with  ), the has a higher discrimination power while is preferable for higher jets. Nonetheless, it should be noted that is also more sensitive to hadronization effects (see section 4.4)

We also investigated the evolution of  with the relative distance in space between the leading and subleading subjets, in particular of the med/vac of its median value and interquartile range . The corresponding RSD are calculated in each bin.

The median ratio (Figure 4) shows a clear separation between models – Q-Pythia and PyQuen(Small) – that broaden the jet structure. In Q-Pythia the two leading subjets become more symmetric with increasing (the ratio is below one and decreases). The same behaviour is seen for PyQuen(Small) up to , but interestingly, vacuum-like behaviour is recovered for larger separations. In contrast, the median ratio in Jewel and PyQuen(Coll) show a dependence on similar to their vacuum references with the ratio nearly independent of the distance between the two leading subjets. The large reduction of the interquartile range () for all observed (Figure 5) in these models provides another clear signature of the jet collimation effect.

From the RSD values the interquartile ratio allows to have a wider spread between the models, although the transition from PyQuen(Small) to vacuum-like behaviour is more noticeable through the median ratio.

### 4.2 Dependence on the choice of subjet clustering algorithm and subjet radius

To investigate the dependence of the proposed observable and its sensitivity to the effects of jet quenching we varied the subjet reconstruction algorithm as well as the subjet radius parameter .

We find no significant differences in  subjet distributions when changing the clustering algorithm from  Cacciari:2008gp () to  Ellis:1993tq () or Cambridge-Achen (C/A) Dokshitzer:1997in (). However, depends, by construction, on the reconstruction algorithm. Figure 6 shows the med/vac ratio of the medians of distributions for different models for   jets with the subjet radius set, as before, to , and the corresponding RSD for completeness. The integer values -1, 0, 1 on the x-axis correspond, respectively, to , C/A, and . Despite the finite differences between clustering algorithms, we find that the power of discrimination between the different models is largely independent of the choice of the algorithm. This observation, together with the results obtained in section 4.4, where we study the effect of different hadronization models on the reconstructed subjets, allow us to conclude that the anti- algorithm provides the most promising option when optimising for jet quenching effects. We therefore adopt this clustering algorithm as the standard setting for the remainder of this work.

Figure 7 shows the dependence on of the med/vac ratio of the medians of  distributions for subjets reconstructed with the  algorithm. The RSD is also shown in the bottom panel, now calculated for each algorithm separately. Here, we find that an increased discrimination between models resulting in jet collimation and the models preferring jet broadening is achieved with . Further, we find a clear difference in the energy distribution inside the jet that results from the different models. On the one hand, Q-Pythia and PyQuen(Small) increase the leading subjet when the subjet radius is increased, indicating broadening of the jet structure. On the other hand, in the models that produce collimated jets by medium effects, the energy in the leading subjet is nearly independent of the chosen subjet radius as it is highly concentrated close to the jet core.

### 4.3 Subjets in dijet pairs

In a back-to-back dijet pair propagating through the QGP, the sub-leading jet has typically lost more energy than its leading partner Milhano:2015mng (). This quenching asymmetry can be combined with  to experimentally further constrain the nature of jet quenching. We have performed an analysis of dijet pairs with  jets within with where the leading jet was required to have   and the recoil jet  . The jets in the pair were required to be separated in azimuth by at least . The med/vac ratios of medians of the  distribution, , as a function of are shown in Figure 8. The upper figures show results for leading jets and the bottom figures those for recoil jets. Again, the upper panels of each figure show the evolution of all models and the bottom panels the corresponding spread quantified through the RSD calculated in each asymmetry bin.

Models that collimate jets towards their core – Jewel and PyQuen (Coll) – display leading jets with a (slightly) enhanced asymmetric subjet momentum balance as compared to their vacuum references, while the recoil jets have a more balanced subjet momentum distribution than in vacuum. In contrast, in models that broaden the jet structure – Q-Pythia, PyQuen (Small), and to a more limited extent PyQuen (Wide) – both leading and recoil jets have a more balanced subjet momentum distribution than in vacuum with the modification significantly stronger for recoil jets. In all cases, the leading jet is modified independently of the pair asymmetry, while the momentum sharing between subjets becomes increasingly balanced (with respect to vacuum) with increasing dijet asymmetry (decreasing ) for recoil jets.

Similar findings are also present in the ratio of the interquartile range of the distributions, for both leading and recoil jets, shown in Fig. 9. However, differences between leading jets are more noticeable through the asymmetry of the  distribution while recoil jets show a larger spread among models through

All these observations are consistent with the findings from section 4.1, figure 2, where jets below always have a  that is smaller than its vacuum reference. This is the preferred kinematic region for the recoil jet in unbalanced dijet systems. Moreover, is fairly constant for jets above , where the leading jet (and recoil for balanced dijet systems) typically comes.

### 4.4 Hadronization effects on the reconstructed subjets

Small radii jets are known to be more sensitive to hadronization effectsDasgupta:2007wa (). For this reason, we investigate the role of different hadronization models in the distributions that were presented so far by using both PYTHIA 8 and HERWIG 7Bahr:2008pv (); Bellm:2015jjp (). The former is based solely on the Lund string fragmentation framework Andersson:1983ia () while the later applies a cluster modelKupco:1998fx () to hadronize the resulting partonic final state to produce hadrons. Although such study is not ideal to accurately assess the uncertainties induced by hadronization effects, including in-medium hadronization modifications Beraudo:2011bh (), it can provide an estimate of the robustness of the proposed observable, .

We have found that is almost insensitive to the hadronization model, with relative differences (taking PYTHIA 8 as reference) smaller than for any choice of subjet radius or clustering algorithm and jets with a transverse momentum  . For low momentum jets ( ), this difference goes up to for and for .

As for the first quartile of the distribution, , the relative change of HERWIG 7 with respect to PYTHIA 8 is for anti- subjets with . Any other choice of clustering algorithm or subjet radius provide a relative difference of independently of the jet transverse momentum.

Finally, the interquartile range, that is able to provide, in general, a larger dispersion between the jet quenching models, is also able to discriminate more among hadronization models. The relative change in low momentum jets ( ) between the interquartile range provided by the two Monte Carlo event generators is around for any clustering algorithm and subjets reconstructed with . For , this change increases to , independently of the clustering algorithm. In high momentum jets ( ) the relative difference is around . The lower bracketing is constantly observed for anti- and small radius subjets while the upper bracketing occurs for  and large radius subjets. For the chosen parameters of this manuscript (anti- subjets with ), the relative change is for any jet with .

The general large sensitivity of the interquartile range to the choice of the hadronization model comes from the fact that this observable is designed to promote the tails of the distributions. While it is the preferable region to tag energy loss modifications imprinted on the jet, it is also the region dominated by a fragmentation pattern that promotes the existence of one () or two () subjets mainly composed by very soft particles. Any modification on the hadronization mechanism would imply a stronger deviation on both (as observed from the increase of the relative differences with respect to ) and .

These observations validate our choice of using reclustered anti- subjets with as to maximize jet quenching phenomena with respect to hadronization effects.

### 4.5 Sub-jet momentum fraction zg and ΔS12

Recent studies of the momentum fraction in jets Sirunyan:2017bsd () prompt for a comparison of with . We have performed an analysis of (vacuum) Pythia jets with   using settings of the Soft Drop algorithm Dasgupta:2013ihk (); Larkoski:2014wba () as in Sirunyan:2017bsd (). Figure 10 shows the as a function of  for two cases: one, where all jets where used; and a second, where jets with between the subjets used to calculate are discarded.  was calculated with in both cases. We find a strong correlation between  and the calculation of when using the cut as in Sirunyan:2017bsd (). Without the cut on the distributions have two dominating structures. One is the diagonal, but the other is largely independent of the  at .

A comment on the differences of behaviour between these two observables in the presence of a QGP is in order. One of the proposed explanations Milhano:2017nzm () for the observed modification of the distribution in nucleus-nucleus collisions relies on the ability of subjets to collect contributions from the QGP. As argued in that work, this QGP backreaction process implies a distinctive increase in size (measured girth) of the subjets. Since in  we explicitly impose, by specifying a radius parameter for the subjet reconstruction, a size for the subjets, contributions from the QGP to both subjets will be of the same order, and thus, will cancel in . This makes  and , well correlated in vacuum, complementary observables in the presence of a QGP that can be used to disentangle the role of QGP backreaction from other dynamical processes that conceivably modify the jet substructure.

## 5 Discussion

Jewel. Results from Jewel are consistent with a jet collimation effect, i.e., most of the radiation is transported outside of the cone leaving the energy core of the jet almost un-modified with respect to the vacuum reference but in a narrower region of phase space. We note, that in this analysis we have used Jewel in its ”recoil-off’” mode which discards the medium partons that interacted with the jet. In this way, the results are independent of the medium-response, whose impact was recently analysed in Milhano:2017nzm (). As a consequence, the is closer to unity for the medium modified jets and it does not change with the chosen subjet subjet radius. Moreover, the comparison of properties of the leading and subleading jets from a dijet event shows that  increases for the leading jets and decreases for the (more strongly medium-modified/low momentum) subleading jet. This is an exclusive characteristic of the jet collimation phenomena and/or collisional energy loss as the same kind of behaviour is observed for PyQuen(Coll). In this model, since , all the energy that is lost outside of the cone is due to elastic energy loss.

Q-Pythia. In Q-Pythia, which is as an implementation of the BDMPS-Z spectrum (without account for destructive interferences), the emission rate is enhanced according to the quenching parameter  leading to a large modifications of the jet inner core. As a consequence, the distribution with a maximum for in vacuum shows a large tail to lower values due to in-medium interactions due to softening of the subjet spectrum (including the leading subjet). Such effect is visible for both leading and subleading jets. Moreover, the medium-induced gluon radiation is evenly distributed in phase space up to very large distances as  mean value is constantly below the vacuum reference without a significant change for .

PyQuen. PyQuen considerations are centred around three angular distributions for the in-medium radiation spectrum. For PyQuen(Small) the finite angle of the radiation () enhances the substructure and the impact on  is qualitatively similar to Q-Pythia. Nonetheless, a striking difference from this model with respect to Q-Pythia is the increasing asymmetry of the subjet structure when biasing the jet sample with . This could be due to the fact that since the radiation is displaced at a finite angle from the leading parton, the more the second hardest subjet is reconstructed away from the jet core, the less probable is to recover the energy. As such, the energy-momentum distribution inside of the jet is located at intermediate distances from the jet core in contrast to what happens in Q-Pythia. A similar, but much milder modification of  is observed for PyQuen(Wide). Since gluon radiation goes as , the radiation is essentially kept near the core with few particles going to very large angles. This angular distribution is similar to the jet vacuum development, which makes this model undistinguishable from the vacuum reference for jets with a large transverse momentum and/or leading jets in dijet systems. On the other hand, the PyQuen(Coll) mode, constrained to elastic energy loss only, affects mainly the softest jet constituents by such elastic collisions, whose energy is ”absorbed” by the medium. As such, jets become more collimated, as it happens in Jewel.

Considerations of the RSD distributions show that the carries the largest discrimination power for high- jets, although there is an associated uncertainty of due to hadronization effects. For low- jets () and/or recoil-jets in dijet systems the use of the and/or the median may prove more advantageous, with hadronization uncertainties that are smaller than for the chosen subjet parameters.

Finally, we reiterate that inclusive or semi-inclusive measurements of nuclear modification factor(s) for jets that fall within a range of  (and ) can provide a rather straightforward insight into the properties of jet quenching (see Zhang:2015trf () for example).

## 6 Conclusions

We have presented observables of subjet structure that by minimizing the impact of the particle backgrounds in heavy-ion collisions are advantageous from the experimental point of view. At the same time, the introduced  quantity preserves the collinear and infrared safety of modern jet algorithms. Using a number of Monte Carlo jet quenching models we have demonstrated that  distribution and can be used as a sensitive tool to discriminate between different quenching mechanisms. We have shown that it is possible to use the quartiles of those distributions, together with the widths and/or use of dijets to make an accurate assessment of the main jet quenching characteristics, in particular, to determine the angular structure of the medium-induced gluon radiation and to investigate further the role of collisional energy loss in the in-medium shower development.

## Acknowledgements

The authors would like to thank to N. Armesto, M. Verweij and K. Zapp for useful discussions.

This work was supported in part by the U.S. Department of Energy, Office of Science, Office of Nuclear Physics, under contract DE-AC02-05CH11231 (MP, XZ) and by the Fundação para a Ciência e Tecnologia (Portugal) under contracts CERN/FIS-NUC/0049/2015, Investigador FCT - Development Grant IF/00563/2012 (JGM) and SFRH/BPD/103196/2014 (LA).

The gross features of the differences between the models have been discussed in terms of  in Sec. 4.1. In this appendix we present the individual  distributions only for completeness and with a limited analysis. The distribution of the fraction of the jet total transverse momentum carried by the leading subjet in jets with   is shown in Fig. 11 for (top panel) Q-Pythia and Jewel, and (bottom panel) PyQuen with its three radiation pattern variants. The vacuum references for each model — Q-Pythia (vac), Jewel (vac), and Pythia (for PyQuen) — are also shown. Clearly Q-Pythia and Jewel modify the  distribution in incompatible directions. As noted for  these observations are consistent with a collimation of jets within Jewel and broadening in Q-Pythia as compared to their vacuum references. For PyQuen (botton panel) we find a clear separation of its different parametrisations of the angular distribution of medium induced radiation. The  distribution (the  fraction carried by the subleading subjet) shown in Fig. 12 is, by definition, limited to the interval. The differences among the  distributions obtained from the different models mirror those observed for . Globally, the  fraction  carried by the subleading subjet reflects the strongly peaked  distribution at large- which necessarily places the average to be below .

## References

• (1) H. Satz, Rept. Prog. Phys. 63, 1511 (2000).
• (2) S.A. Bass, M. Gyulassy, H. Stoecker, W. Greiner, J. Phys. G25, R1 (1999).
• (3) E.V. Shuryak, Phys.Rept. 115, 151 (1984)
• (4) J. Cleymans, R. Gavai, E. Suhonen, Phys.Rept. 130, 217 (1986).
• (5) J.D. Bjorken, (1982)
• (6) U.A. Wiedemann, pp. 521–562 (2010).
• (7) K.M. Burke, et al., Phys. Rev. C90(1), 014909 (2014).
• (8) Y. Mehtar-Tani, J.G. Milhano, K. Tywoniuk, Int. J. Mod. Phys. A28, 1340013 (2013).
• (9) J. Adams, et al., Nucl.Phys. A757, 102 (2005).
• (10) K. Adcox, et al., Nucl.Phys. A757, 184 (2005).
• (11) I. Arsene, et al., Nucl.Phys. A757, 1 (2005).
• (12) B. Back, et al., Nucl.Phys. A757, 28 (2005).
• (13) K. Aamodt, et al., Phys. Lett. B696, 30 (2011).
• (14) K. Aamodt, et al., Phys. Rev. Lett. 108, 092301 (2012).
• (15) S. Chatrchyan, et al., Eur.Phys.J. C72, 1945 (2012).
• (16) G. Aad, et al., Phys.Rev.Lett. 105, 252303 (2010).
• (17) S. Chatrchyan, et al., Phys.Rev. C84, 024906 (2011).
• (18) J. Adam, et al., Phys.Lett. B746, 1 (2015).
• (19) L. Adamczyk, et al., Phys. Rev. C96(2), 024905 (2017).
• (20) J. Adam, et al., JHEP 09, 170 (2015).
• (21) G. Aad, et al., Phys.Lett. B739, 320 (2014).
• (22) S. Chatrchyan, et al., Phys.Lett. B730, 243 (2014).
• (23) S. Chatrchyan, et al., JHEP 1210, 087 (2012).
• (24) A.M. Sirunyan, et al., (2017)
• (25) L. Cunqueiro, Nucl. Phys. A956, 593 (2016).
• (26) M. Rubin, JHEP 05, 005 (2010).
• (27) Z.B. Kang, F. Ringer, W.J. Waalewijn, JHEP 07, 064 (2017).
• (28) A.J. Larkoski, S. Marzani, J. Thaler, Phys. Rev. D91(11), 111501 (2015).
• (29) A. Larkoski, S. Marzani, J. Thaler, A. Tripathee, W. Xue, Phys. Rev. Lett. 119(13), 132003 (2017).
• (30) C. Frye, A.J. Larkoski, M.D. Schwartz, K. Yan, JHEP 07, 064 (2016).
• (31) S. Marzani, L. Schunk, G. Soyez, JHEP 07, 132 (2017).
• (32) S. Marzani, L. Schunk, G. Soyez, Eur. Phys. J. C78(2), 96 (2018).
• (33) C. Collaboration, (2017)
• (34) M. Aaboud, et al., (2017)
• (35) X. Zhang, L. Apolinário, J.G. Milhano, M. Ploskon, Nucl. Phys. A956, 597 (2016).
• (36) S. Chatrchyan, et al., JHEP 05, 090 (2013).
• (37) G. Aad, et al., JHEP 05, 128 (2012).
• (38) A. Tripathee, W. Xue, A. Larkoski, S. Marzani, J. Thaler, Phys. Rev. D96(7), 074003 (2017).
• (39) K. Kauder, Nucl. Part. Phys. Proc. 289-290, 137 (2017).
• (40) M. Cacciari, G.P. Salam, G. Soyez, JHEP 04, 063 (2008).
• (41) M. Cacciari, G.P. Salam, G. Soyez, Eur. Phys. J. C72, 1896 (2012).
• (42) B. Abelev, et al., JHEP 1203, 053 (2012).
• (43) M. Cacciari, G.P. Salam, Phys. Lett. B659, 119 (2008).
• (44) N. Armesto, L. Cunqueiro, C.A. Salgado, Eur. Phys. J. C63, 679 (2009).
• (45) K.C. Zapp, Eur. Phys. J. C74(2), 2762 (2014).
• (46) I.P. Lokhtin, A.M. Snigirev, Eur. Phys. J. C45, 211 (2006).
• (47) T. Sjostrand, S. Mrenna, P.Z. Skands, JHEP 05, 026 (2006).
• (48) R. Baier, Y.L. Dokshitzer, A.H. Mueller, S. Peigne, D. Schiff, Nucl.Phys. B484, 265 (1997).
• (49) A. Dainese, C. Loizides, G. Paic, Eur. Phys. J. C38, 461 (2005).
• (50) L. Apolinario, N. Armesto, L. Cunqueiro, JHEP 02, 022 (2013).
• (51) L.D. Landau, I. Pomeranchuk, Dokl. Akad. Nauk Ser. Fiz. 92, 535 (1953)
• (52) A.B. Migdal, Phys. Rev. 103, 1811 (1956).
• (53) K.C. Zapp, Phys. Lett. B735, 157 (2014).
• (54) S.D. Ellis, D.E. Soper, Phys. Rev. D48, 3160 (1993).
• (55) Y.L. Dokshitzer, G.D. Leder, S. Moretti, B.R. Webber, JHEP 08, 001 (1997).
• (56) J.G. Milhano, K.C. Zapp, Eur. Phys. J. C76(5), 288 (2016).
• (57) M. Dasgupta, L. Magnea, G.P. Salam, JHEP 02, 055 (2008).
• (58) M. Bahr, et al., Eur. Phys. J. C58, 639 (2008).
• (59) J. Bellm, et al., Eur. Phys. J. C76(4), 196 (2016).
• (60) B. Andersson, G. Gustafson, G. Ingelman, T. Sjostrand, Phys. Rept. 97, 31 (1983).
• (61) A. Kupco, pp. 292–300 (1998)
• (62) A. Beraudo, J.G. Milhano, U.A. Wiedemann, Phys. Rev. C85, 031901 (2012).
• (63) M. Dasgupta, A. Fregoso, S. Marzani, G.P. Salam, JHEP 09, 029 (2013).
• (64) A.J. Larkoski, S. Marzani, G. Soyez, J. Thaler, JHEP 05, 146 (2014).
• (65) J.G. Milhano, U.A. Wiedemann, K.C. Zapp, (2017).
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters