Jet shape modification at LHC energies by JEWEL

Jet shape measurements are employed to explore the microscopic evolution mechanisms of parton-medium interaction in ultra-relativistic heavy-ion collisions. In this study, jet shape modifications are quantified in terms of the fragmentation function , relative momentum , density of charged particles , jet angularity , jet momentum dispersion , and for proton-proton (pp) collisions at 0.9, 2.76, 5.02, 7, and 13 TeV, as well as for lead-lead collisions at 2.76 TeV and 5.02 TeV by JEWEL. A differential jet shape parameter is proposed and studied at a smaller jet radius . The results indicate that the medium has the dominant effect on jet shape modification, which also has a weak dependence on the center-of-mass energy. Jet fragmentation is enhanced significantly at very low , and fragmented jet constituents are linearly spread to larger jet-radii for . The waveform attenuation phenomena is observed in , , and distributions. The results obtained for from to , where the wave-like distribution in collision is ahead of collisions at small jet-radii, indicates a strong medium effect.


Introduction
One of the objectives of jet physics is to explore the microscopic properties of hot-dense quantum chromodynamic (QCD) matter, quark-gluon plasma (QGP), generated in ultra-relativistic heavy-ion collisions [1][2][3]. Jets produced by partons from colliding nuclei in the early stage travel through the QCD matter carrying multi-scale physics evolution information. The interplay between elementary scattering, subsequent branching process, and strongly coupled parton-medium interactions can lead to the modification of jet shapes with respect to jet fragmentation in the absence of the medium. The jet structure modifications resulting from proton-proton (pp) collisions to heavy-ion collisions can be investigated in terms of a set of jet shape quantifications, which will shed light on the jet energy loss mechanism in the medium, color coherence, and fundamental medium properties.
Latest results from the LHC experiments conducted on charged particle nuclear modification factor measurement in a wide range of momenta exhibit sectionalized behavior that gives rise to an apparent soft and hard pQCD region [4]. The observations of more highly unbalanced di-jet with increasing event centrality [5], the suppression of inclusive jet yield by about a factor of two in central heavy-ion collisions relative to peripheral collisions, as well as the correlation of the jet suppression with missing transverse momentum [6][7][8] indicate that the jet energy redistributes within the jet in medium relative to the jet in vacuum, and that there is a non-zero fraction of jet energy loss especially at large angles relative to the jet axis. The fragmentation yield is found to decrease at z z r < 0.2 intermediate and increase at small in central collisions relative to peripheral collisions [9,10]. The azimuthal dihadron correlation excess yield is found to be more pronounced in the sub-leading jet, predominantly from several to 20 GeV arising from soft particles [11]. In [12], the jet shapes modifications are studied at small jet cone sizes ( ) and found an inconsistent with a fully coherent energy loss picture. These results provide abundant inputs for theoretical calculations and phenomena modeling with regard to jet-medium interactions.

Pb + Pb
In theory, most studies of the phenomena concerned with jet energy loss are based on pQCD, which can efficiently explain data at the intermediate momentum range, but is incapable of covering the entire kinematic range in the non-pQCD region [13,14]. One important reason is the treatment of gluon radiation by assuming collinear radiation and model dependent on momentum exchange among the medium and parton dynamical evolution for leading order or next leading order contribution, especially for central collisions [15]. Moreover, several physical mechanisms involve co-existence and competition of several effects when partons traverse the medium, such as recoil and non-recoil, coherence and de-coherence, depending on QCD local equilibrium temperature, medium density, path length, gluon formation time, and momentum exchange etc., which redistribute the parton energy at a relative larger angle with respect to its initial directions [16].
Jet-medium interactions cover the entire kinematic range and include perturbative and non-perturbative effects at the LHC energies. It is essential to study the jet shape modification in various collision systems with a large set of jet shape quantities. In this study, the measurement of the jet shape observables, such as the fragmentation function , relative momentum , density of charged particles , jet angularity , jet momentum dispersion , and the difference between leading and sub-leading constituent transverse momenta for proton-proton collisions at 0.9, 2.76, 5.02, 7, and 13 TeV, as well as for lead-lead collisions at 2.76 TeV and 5.02 TeV, are systematically studied by JEWEL [17]. One additional differential jet shape parameter of is proposed and studied at to investigate the jet shape evolution.

Jet shape measurements
Jet shape measurements have been suggested for the exploration of the microscopic evolution mechanism of parton-medium interaction in ultra-relativistic heavy-ion collisions. For different physical motivations, several jet shape observations have been proposed using jet constituents, overall jet-by-jet quantities, and jet cluster historical information. For example, the fragmentation function, p T Pb + Pb pp missing method and jet-track angular correlations measurements using intra or inter-jet distributions are dedicated to investigating the longitudinal share of the energy within the jet and large angle radiations. The overall jet shape observations built on the jet-by-jet function of jet constituent four-momenta, such as jet mass, jet granularity, and jet momentum dispersion etc., could probe the jet energy loss dependence on large angle soft particle emission. Jet grooming measurements based on historical information of the jet cluster could locate the splitting phase space where medium-induced effects are expected [18,19]. This study presents a jet shape investigation in and collisions by inclusive constituents of the jet and jet-by-jet quantity to describe jet evolution and jet energy loss within QCD matter.
In [20], ATLAS presented the fragmentation function and transverse profiles in = 7 TeV collisions in a wide jet momentum range of 25 GeV/c 500 GeV/c and found discrepancies between various models and data. Fragmentation functions show a reduction yield at and an enhancement for in 2.76 TeV collisions by ATLAS [9]. This measurement was also conducted in ALICE [21] from dihadron correlations in = 7 TeV and = 5.02 TeV collisions, where no significant cold nuclear matter effect was observed. In [12,22], the jet shape modifications were studied in small-radius jets to discriminate the relative quark and gluon jet fractions, which also suggested that the medium was able to resolve the jet structure at angular scales smaller than . However, a detailed and systematic study is highly required to quantify the gluon induced fragmentation modifications, better constrain energy loss models, and isolate the quark and gluon jet fraction at the intermediate momentum range for multiple collision systems. In this study, the below observations are used.
The jet fragmentation function is defined as the probability that a hadron carries a longitudinal momentum fraction of in the experiment. is defined: which connects color partons and colorless hadrons to constrain QCD-motivated models from experiments. The quantity of is measured as: Here is the number of charged particles in the jet, and is the number of selected jets to be used for normalization. Two additional quantities and the density of charged particle are also studied. is the momentum of charged particles in a jet transverse to that of the jet's axis defined in Eq. (3), and its distribution is given in Eq. (4): The density of the charged particle in space is measured as a function of angular distance of charged particles from the jet axis, given by:

LeS ub
As a complementary study, the jet-by-jet quantities, namely the first radial moment or angularity (or girth), , the momentum dispersion, , and the difference between the leading and sub-leading track transverse momentum, are also studied. and are related to the moments of the so-called generalized angularities [23].
is not an IRC-safe observable but indicates robustness against soft background particles contributions. The angularity is defined as is the transverse momentum of the -th constituent and is the distance in ( , ) space between constituent and the jet axis. This shape is sensitive to the radial energy profile of the jet.
The momentum dispersion is defined as: LeS ub p T p T is defined as the difference of the leading track and sub-leading track , D girth girth r r + δr For jets with small radii r, it contributes to isolated pure energy loss effect from other medium effects and correlated background. Hence, the differential jet-by-jet quantity is proposed, which is defined as jet's distribution of inside an annulus of inner radius and outer radius around the jet axis, which is directly related to the jet energy topological structure and its evolution.

Simulation method
In this study, JEWEL [17] is used to simulate jet production, QCD scale evolution, and re-scattering of jets in heavy-ion collisions based on pQCD. JEWEL describes the jet evolution and jet-medium interactions simultaneously and dynamically, based on leading-order matrix elements and the parton shower method. All partons belonging to partons showers initiated by hard scattered partons undergo collisions with thermal partons in the medium, leading to both elastic and radiative energy loss. Soft gluon radiation, recoil effect, and scattering processes are governed by formation time, which has been shown to agree with analytical calculations within the appropriate limits. The Landau-Migdal-Pomeranchuk (LPM) effect is included by generalizing the probabilistic formulation in the Eikonal limit to general kinematics. For the heavy-ion environment, the medium density profile is based on a longitudinally expanding Glauber overlap; the local temperature is sampled to determine the density of scattering centers and their momentum distribution. The initial time = 0.6 fm and temperature = 0.4 GeV, the critical temperature = 170 MeV are fixed for central lead-lead collisions. Jets are reconstructed by using all final detectable charged and neutral particles with the antikt algorithm provided by the FastJet [24] package within and for the fragmentation distribution and jet-by-jet quantities measurements in separately. The analysis takes into account the recoil effect for heavy-ion collision, and the constituent subtraction method was used to remove the contribution of medium fragments from the parton energy. distributions indicate a strong enhancement at and a reduction at compared from to collisions with the recoil effect. For collisions and collisions without the recoil effect, the enhancement or reduction of depends on the collisional energy and participating nuclei.

Results and discussion
indicates the same performance that jet constituents spread to large jetradius with the increase of the center-of-mass energy of collisions, and the medium effect is significant for the modification at relatively large angles. The distribution decreases exponentially with the jet-radius. The distribution is independent of the collisional energy, however it is strongly modified by QCD medium and exhibits an obvious increase in the approximate range of .
To better illustrate the differences between the simulation and data, Figs. 2, 3, and 4 show the ratio plots with the reference data of = 7 TeV collisions for , , and separately at four jet momentum intervals It's obvious that in Fig. 2 is suppressed from to collisions by the so-called nuclear modification factor of about 15% at [25,40] GeV/c and of about 10% at [80,110] GeV/c when . When , the nuclear modification factor increases by a factor of 2 to 3 with a decrease in . For 60 GeV/c, the suppression increases at very low and then decreases when . For collisions, the ratios are all enhanced at low and then kept almost at 1 for large and 40 GeV/c, while for high momentum jets with , the ratio is suppressed by up to 50%. The fragmentation becomes harder for high and weak dependence on beam energies and medium effects.
The nuclear modification factor of charged particle density indicates a limited suppression in the small jet radius and a linear increase with the increase of the jet radius from to collisions. The slope is about 2, which is the same level of increase in . Whereas for the collisions, the ratio is almost constant and slightly differs with the center-of-mass energy. The weak dependence on the jet momentum and the center-of-mass energies of collisions, as well as the strong dependence on the medium are observed in the ratio plot of in Fig. 4. The nuclear modification factor of shows an enhancement up to 50% at low followed by a sharp decrease at larger . At very low in collisions, a fluctuation in the curve may indicate a flow effect in collisions, since and are directly related to the non-perturbative hadronization process controlled by perturbative QCD radiation. With the increase of , soft gluon radiation is important, as it will contribute to jet broadening and cause the As demonstrated in [12], in the present study, these quantities are also dedicated to the study of smaller jet-radii ( ) to explore the medium effect and jet evolution. The at small jet-radii with different configurations is presented in Fig. 5. The trend obtained by simulation is similar to the data, although it could not efficiently reproduce the data. The angularity is the first radial moment of the jets, which is very sensitive to its radial energy profile. The girth distributions indicate that with the increase of jet-radius, the peak value of shifts to larger values.
LeS ub LeS ub Interestingly, the distribution decreases exponentially with at different jet-radii. For small jetradii, asymmetric parton splitting is more probable, while symmetric branching can be found at large jet-radii, which is directly related to the gluon and quark jet fraction. The momentum dispersion distribution shows that is a local maximum at and then decreases with the increase of the jet radius at large , implying a large energy fraction within smaller jet-radius. Figure 6 shows the differential distribution iñ 0% -10% of central collisions at 2.76 TeV in [40,60] GeV/c. It is found that the differential girth evolution from the small jet-radius interval to the larger behaves as a waveform, which attenuates and becomes flatter at larger values. Since the small radius jet shows larger jet energy loss and smaller granularity.
is sensitive to soft particles or the correlated background arising from hot QCD medium. This wave-like pattern and its evolution may potentially be used to investigate the medium properties, such as the coherence effect, medium density, and flow effect, etc. with data [12].