Mass suppression effect in QCD radiation and hadron angular distribution in jet

The finite mass of the heavy quark suppresses the collimated radiations, which is generally referred to as the dead cone effect. In this paper, we study the distribution of hadron multiplicity over the hadron opening angle with respect to the jet axis in various flavors of jets. The corresponding measurement can be the most straightforward and simplest to explore the dynamical evolution of the radiations in the corresponding jet, which can expose the mass effect. We also propose the transverse energy-weighted angular distribution which sheds light on the interplay between perturbative and nonperturbative effects in the radiation. With Monte-Carlo simulations, our calculation shows that the dead cone effect can be clearly seen by taking the ratio between the b jet and the light-quark (inclusive) jet, promising to be measured at the LHC in the future.


I. INTRODUCTION
Measuring jet production offers unique opportunities to study the perturbative and non-perturbative behavior of QCD [1][2][3].Heavy flavor jets have been used to verify perturbative QCD and explore the non-perturbative effect, thereby enhancing our understanding of jet evolution in vacuum or QCD media [4,5].Currently at the LHC, identification techniques of heavy flavor jets make it possible to discriminate the jets originating from b or c quarks and those from light flavor quarks or gluons.Experimental collaborations at the LHC have measured heavy flavor jet production [6][7][8][9][10][11][12][13][14][15][16] and more results are expected in the near future, especially from LHC Run 3.
Many efforts on the theoretical side have been devoted to studying heavy flavor jets at various colliders.The inclusive p T spectrum of heavy flavor jets [17,18] can be predicted with the help of the semi-inclusive jet function [19,20].The heavy flavor jet with high transverse momentum can be used to understand the evolution of a massless quark when the mass is small enough to be ignored in comparison to the jet energy.In addition, the mass impact on the perturbative and non-perturbative nature during the evolution of heavy flavor quarks can be monitored by measuring the radiation pattern inside jets.In this regard, the most famous phenomenon is the dead cone effect [21][22][23], which is a direct consequence of the suppression of the collinear radiation due to the mass of the radiator, i.e., the heavy quark.In recent years, the mass effect using heavy flavor jets has drawn a lot of attention in both theoretical and experimental studies.This mass effect has been analyzed in gauge theory models such as QED and QCD [24][25][26][27][28][29][30].Besides, many studies have been devoted to quantifying the dead cone effect in heavy-ion and electron-ion collisions [18,[31][32][33][34][35][36][37].
In general, the non-zero mass of heavy quarks can control the infrared behavior of the radiation, leading to a specific perturbative radiation effect.In the collinear limit, the splitting of a massive quark can be described by effective field theory as [38] where k ⊥ is the transverse momentum of the emitted gluon and z is the energy fraction of the gluon relative to the parent massive quark.According to Eq. ( 1) if we compare with the massless quark splitting, the smallangle radiation is significantly suppressed by the mass term, i.e., the dead cone effect [22].Jet properties and substructures have been widely used to study this effect.In [28,39,40] the authors investigated z g and θ g distributions for the groomed b jet that the b-quark mass effect plays an important role in the collinear splitting pattern.Many other developments have been proposed on jet substructures for heavy flavor jet production [34,41,42].There are also studies that applied heavy flavor jet observables at the electron-ion collider [43]; see also Ref. [44] and the references therein.By inspecting the splitting inside the jet, the dead cone effect has already been measured using the substructure of the charm jet by AL-ICE [13][14][15].
Eq. ( 1) also shows that generally the mass effect is more significant at relatively low energy scales/small angles.To obtain a better energy/angular resolution for this effect, the most straightforward approach is to measure the angular distribution of hadron multiplicity in the jet.The dead cone effect can be explored by comparing the hadron angular distribution between heavy flavors and light quark jets.As a physical phenomenon that typically occurs at preconfinement scale [45,46], hadron multiplicity is not an infrared safe observable.However, we expect that the angular distribution of the hadrons preserves most of the perturbative effect from the heavy quark mass.We also propose a transverse energy weighted angular distribution that connects the perturbative and non-perturbative multiplicity distributions.The proposed observable can be measured at the LHC and provides a new approach for investigating the mass effect for QCD radiation and jet formation.Furthermore, this observable can be used to test or tune the hadronization models of Monte-Carlo event generators.
The rest of the paper is organized as follows.We provide a definition of the observable in Sec.II.Section III presents numerical results from Monte-Carlo event generators and discusses the mass effect on jet evolution.Section IV concludes the paper..

II. DEFINITION OF THE OBSERVABLE
In general, the non-zero mass of a radiator or the radiated particles can control the infrared behavior of the radiation.Given that heavy quark mass does not originate from confinement dynamics, it can lead to a specific perturbative radiation effect.The averaged charged particle multiplicities in e + e − collision were used to investigate the mass effect [23,[47][48][49].Currently, leveraging the LHC and high luminosity LHC, it is possible to take a closer look inside jets.In this paper, we propose the simplest approach to comprehensively analyze the dead cone effect for heavy flavor jets: an averaged multiplicity distribution defined as where θ is the opening angle between the jet axis and the moving direction of the charged hadron ch and dP ch /dθ is the probability distribution as a function of θ for a charged hadron ch.According to the dead cone effect the collinear radiation inside heavy flavor jets is suppressed at a small angle.A detailed study of d⟨N ch ⟩/dθ for jets with various transverse momenta can be used to quantitatively identify the energy scale of the dead cone effect.Although this observable is not infrared safe, the distribution is expected to reflect the pattern from perturbative radiations; in particular, the difference between heavy flavor and light quark jets.This observable may depend on the definition of the jet axis.In this study, we used the traditional energy combination scheme to retain the correlation between the directions of the jet axis and momentum of the parent parton.
Regarding the non-perturbative nature of d⟨N ch ⟩/dθ , it is interesting to investigate to what extent this angle distribution can be affected by nonperturbative QCD.Therefore, we introduce a variation to Eq. ( 2) where p T,ch and p T,j are the transverse momenta for the charged hadron and jets, respectively; κ is a free parameter.When κ = 0, it is reduced to the multiplicity distributions inside jets.For κ = 1, d⟨N ch ⟩(κ = 1)/dθ measures the θ dependence of the energy deposit inside the jet cone and r 0 dθ d⟨N ch ⟩(κ=1) dθ corresponds to the infrared safe observable jet shape, or jet transverse energy profile [50].For a variation in the range of 0 ≤ κ ≤ 1 a bridge between infrared unsafe and safe observables is established, which can be utilized to probe the nonpertubative effect.

III. NUMERICAL RESULTS
Numerical results for d⟨N ch ⟩(κ)/dθ are obtained from simulations with PYTHIA [51,52].The default setting in PYTHIA 8.306 was adopted for parton showers and hadronizaiton.The simulation was performed for dijet production at the 13 TeV LHC.To investigate the physics of jets with different flavors, the multiparton interaction was switched off.The jets were constructed with anti-k T algorithms [53] and jet radius R = 0.4 using only the charged tracks with p T > 1 GeV and |η| < 2.0.The recombination of jets was achieved using the package FastJet [54].We classify the jet as a b jet if there was at least one B hadron in its component.
Figure 1 shows the distributions of d⟨N ch ⟩/dθ for quark, gluon, and b jets for various intervals of jet transverse momentum.The width of the bands corresponds to the uncertainties from varying the scales in parton showers by a factor of two.The integration over the θ/R axis provides the average charged multiplicity of the jet.As expected, he radiation from light quarks and gluons is greater.In particular, there is more radiation for gluon jets at large angles, leading to a broader distribution.At a large angle, the b jet behaves like the light quark jet; this is more evident for the high p T jet.A clear suppression of the radiation can be observed for b jet at small angle radiation.Needless to say, the suppression stems from the dead cone effect.Roughly speaking, the distributions for the quark and b jet are expected to be different when the typical scale of the collinear splitting inside jets p T,j θ z(1 − z) ≤ p T,j θ/41 is close to the b quark mass m b .In the first and last panels of Fig. 1 the black lines show the multiplicity distribution for b jets if we do not count the b hadrons inside jets.We find that b hadron tends to stay close to the jet axis.As a result, there is a growing trend for very small angles, particularly when the jet p T,j is small.
One of the interesting features of the multiplicity distribution is the scaling behavior in the limit θ → 0. Figure 2 shows the multiplicity distributions for extremely small opening angles for jets with the transverse momentum 30 < p T,j < 40 GeV and 600 < p T,j < 800 GeV.When p T,j θ ≤ Λ QCD , the phase space is extremely limited and the distributions are supposed to be dominated by non-perturbative features of QCD.Surprisingly, for jet multiplicities, we find that the quark, gluon, and b jets exhibit a similar scaling behavior in the small θ region.This phenomenon can be explained by the existence   of uniformly distributed hadrons in the collinear limit of jets; as a result d⟨N ch ⟩/dθ ∝ θ, similar to energy-energycorrelators first reported by Ref. [55].This interesting feature might originate from kinematics, given that the dynamic evolu-tion is frozen below Λ QCD .Consistent with what we find in Fig. 1, the multiplicity distribution of the b jet is suppressed and then enhanced when de-creasing the angle relative to the jet axis.The enhancement in the small angle for the b jet mainly comes from that the B hadron tends to stay close to the jet axis.By varying κ in Eq. ( 2) we observe relatively small κ dependence in the small angle region for b jet.From Figs. p T jets.
To investigate the mass effect in detail, we present the ratio of the multiplicity distributions between b and quark/inclusive jets at the LHC in Fig. 3.The mass effect is supposed to become smaller with increasing θ because the typical QCD scale for the splitting is larger, while at an extremely small angle, there is an enhancement observed for heavy flavor jet given that the collinear radiation is suppressed and the parent particle tends to stay close to the jet axis.As a result of the interplay between these effects, there are dips in the distributions of the ratio.The dead cone effect can be captured by the position of the dip, which depends on the transverse momentum of the jet.The right plot of fig. 3 displays the ratios between the distributions of b and inclusive jets.We found that the ratios are similar to those between the b and quark jet; this is an approach to test the dead cone effects directly.A measurement of the multiplicity distributions of jets can be used to investigate the mass effect on the QCD dynamic evolution.
For κ < 1, d⟨N ch ⟩/dθ is not an infrared-safe observable, while for κ = 1 d⟨N ch ⟩/dθ corresponds to the differential jet shape and perturbatively calculable.The quantitative effect from non-perturbative physics can be explored by varying κ in Eq. ( 2).The distribution of d⟨N ch ⟩(κ)/dθ heavily depends on κ; for larger values of κ, the distribution is smaller.To further analyze the bquark mass and the nonperturbative effect, we present the dependence of the ratio between b-jet and lightquark/inclusive jet with 30 < p T,j < 40GeV in Fig. 4. Remarkably, we found that in the large θ region the nonperturbative effect seems to be canceled in the ratios, thereby setting a guideline for the nonperturbative corrections in heavy flavor jet substructures.In the small θ region the differences between distinct κ settings are large, as expected.The results for different values of κ demonstrate a good and smooth transition of nonperturbative and perturbative QCD.As shown in the right plot of fig. 4, the ratio between b jet and inclusive jets still keeps the mass effect as for the case of the ratio between b and quark jets.

IV. CONCLUSION
In this paper, we present the most straightforward and simplest approach to expose the dead cone effect of heavy flavor jets: a multiplicity distribution over the opening angle between the hadron inside the jet and the jet axis.Although this multiplicity is not infrared safe, we expect that this distribution reveals the mass effect on both perturbative and non-perturbative evolutions of the heavy flavor jets.To address the non-perturbative effect, we propose a transverse energy weighted multiplicity distribution that sheds light on the interplay between perturbative and non-perturbative effects.
We present simulations of the multiplicity distributions for b, light quark, and gluon originated jets with PYTHIA using charged particles in the events.In comparison to light quark and gluon jets, the radiation of b jets is suppressed for θ ≈ m b /p T j ; this is the dead-cone effect.Given the b jets resemble a hard B hadron dressed with some soft radiations, we observed an enhancement of the distribution and a scaling behavior indicating uniformly distributed hadrons at an extremely small angle.By calculating the ratio between the b jet and light quark or inclusive jets, the dead cone effect can be clearly observed in the simulations and is expected to be measured at the LHC in the future.We also examined the κ dependence of the ratio of the multiplicity distributions between the b jet and light quark or inclusive jets, and found that the non-perturbative effect cancels at large angles and is important in small-angle regions.One important fact is that the dead cone or mass suppression on small-angle radiation is not removed by increasing the quark energy, which is also clearly observed from the b fragmentation function (see recent study [56]).
Last but not least, we remark that the multiplicity distribution reveals the mass effect during the dynamical evolution of heavy flavor jets.The shape of the distributions indicates the energy scale of the dead cone effect.It would be interesting to apply this observable to heavy ion collisions, which can be used to reveal the scale of the interactions of colored partons with quark-gluon plasma.

FIG. 1 :
FIG. 1: Angle distribution of charged-particle multiplicity in different intervals of jet transverse momentum.The x-axis is the opening angle between the charged particle and the jet axis which is normalized by the jet radius R = 0.4.The red, blue, and green bands represent the b, quark, and gluon jets.The bands are Monte Carlo uncertainties.The black lines are the multiplicity distributions for b-jet without counting b hadrons.

1 and 2 FIG. 3 :FIG. 4 :
FIG. 3: Ratios of the multiplicity distributions of b-jet to quark jet (left) and b jet to inclusive jet (right) with various jet transverse momenta p T,j .