Direct observation of the dead-cone effect in QCD

In particle collider experiments, elementary particle interactions with large momentum transfer produce quarks and gluons (known as partons) whose evolution is governed by the strong force, as described by the theory of quantum chromodynamics (QCD). These partons subsequently emit further partons in a process that can be described as a parton shower which culminates in the formation of detectable hadrons. Studying the pattern of the parton shower is one of the key experimental tools for testing QCD. This pattern is expected to depend on the mass of the initiating parton, through a phenomenon known as the dead-cone effect, which predicts a suppression of the gluon spectrum emitted by a heavy quark of mass $m_{\rm{Q}}$ and energy $E$, within a cone of angular size $m_{\rm{Q}}$/$E$ around the emitter. Previously, a direct observation of the dead-cone effect in QCD had not been possible, owing to the challenge of reconstructing the cascading quarks and gluons from the experimentally accessible hadrons. We report the direct observation of the QCD dead cone by using new iterative declustering techniques to reconstruct the parton shower of charm quarks. This result confirms a fundamental feature of QCD. Furthermore, the measurement of a dead-cone angle constitutes a direct experimental observation of the non-zero mass of the charm quark, which is a fundamental constant in the standard model of particle physics.


Introduction
In particle colliders, quarks and gluons are produced in high-energy interactions through processes with large momentum transfer, which are calculable and well described by quantum chromodynamics (QCD). These partons undergo subsequent emissions, resulting in the production of more quarks and gluons. This evolution can be described in the collinear limit by a cascade process known as a parton shower, which transfers the original parton energy to multiple lower energy particles. This shower then evolves into a multi-particle final state, with the partons combining into a spray of experimentally detectable hadrons known as a jet [6]. The pattern of the parton shower is expected to depend on the mass of the emitting parton, through a phenomenon known as the dead-cone effect, whereby the radiation from an emitter of mass m and energy E is suppressed at angular scales smaller than m/E, relative to the direction of the emitter. The dead-cone effect is a fundamental feature of all gauge field theories (see Ref. [3] for the derivation of the dead cone in QCD).
The dead-cone effect is expected to have sizeable implications for charm and beauty quarks, which have masses of 1.28 ± 0.02 GeV/c 2 and 4.18 +0.03 −0.02 GeV/c 2 (Ref. [1]) in the minimal subtraction scheme, respectively, at energies on the GeV scale. The emission probability in the collinear region, which is the divergent limit of QCD at which the radiation is most intense, is suppressed with increasing mass of the quark. This leads to a decrease in the mean number of particles produced in the parton shower. The DELPHI Collaboration at the LEP e + e − collider measured the multiplicity difference between events containing jets initiated by heavy beauty quarks and those containing light quarks (up, down or strange). They found that the differences depend only on the quark mass [7], which was attributed to the suppression of collinear gluon radiation from the heavy quark because of the dead-cone effect. A measurement of the momentum density of jet constituents as a function of distance from the jet axis was also performed by the ATLAS collaboration at CERN [8], which pointed to a depletion of momentum close to the jet axis that was ascribed as a consequence of the dead-cone effect. The mass of the beauty quark was also estimated through a phenomenological fit to the measured data [9]. As hard (large transverse momentum) emissions are preferentially emitted at small angles, and are therefore suppressed for massive emitters, heavy quarks also retain a larger fraction of their original momentum compared to lighter quarks, leading to a phenomenon known as the leading-particle effect. This has been well established experimentally, with the fraction of the jet momentum carried by the leading (highest transverse momentum) hadron containing a charm or beauty quark (heavy-flavour hadron) in jets, peaking at 0.6-0.7 and 0.8-0.9, respectively, whereas the corresponding fraction carried by the leading hadron in light quark-initiated jets peaks at smaller values [10][11][12][13][14].
Until now, a direct experimental measurement of the dead-cone effect has been subject to two main challenges. First, the dead-cone angular region can receive contributions from hadronization effects or particles that do not originate from the gluon radiation from the heavy-flavour quark, such as the decay products of heavy-flavour hadrons. The second difficulty lies in the accurate determination of the dynamically evolving direction of the heavy-flavour quark, relative to which the radiation is suppressed, throughout the shower process. The development of new experimental declustering techniques [4] enables these aforementioned difficulties to be overcome by reconstructing the evolution of the jet shower, giving access to the kinematic properties of each individual emission. These techniques reorganise the particle constituents of an experimentally reconstructed jet, to access the building blocks of the shower and trace back the cascade process. Isolated elements of the reconstructed parton shower that are likely to be unmodified by hadronization processes provide a good proxy for real quark and gluon emissions (splittings). These reclustering techniques have been demonstrated in inclusive (without tagging the initiating parton flavour) jets to successfully reconstruct splittings that are connected to or that preserve the memory of the parton branchings. This is demonstrated by measurements such as the groomed momentum balance [15][16][17][18], which probes the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi splitting function [19], and the Lund plane [20], which exposes the running of the strong coupling with the scale of the splittings. An experimental method to expose the dead cone in boosted top-quark events was also proposed in Ref. [21].
Reclustering techniques are extended in this work to jets containing a charm quark based on the prescription given in Ref. [22]. These jets are tagged through the presence of a reconstructed D 0 meson amongst their constituents, which has a mass of 1.86 GeV/c 2 (Ref. [1]) and is composed of a heavy charm quark and a light anti-up quark. The measurement is performed in proton-proton collisions at a centre-of-mass energy of √ s = 13 TeV at the Large Hadron Collider (LHC), using the ALICE (A Large Ion Collider Experiment) detector. Further details of the detector apparatus and data measured can be found in the Methods. As the charm-quark flavour is conserved through the shower process, this provides an opportunity to isolate and trace back the emission history of the charm quark. In this way, by comparing the emission patterns of charm quarks to those of light quarks and gluons, the QCD dead cone can be directly revealed.
2 Selecting jets containing a D 0 meson To select jets initiated by a charm quark, through the presence of a D 0 meson in their list of constituents, the D 0 mesons and jets need to be reconstructed in the events. The D 0 -meson candidates (and their anti-particles) were reconstructed in the transverse-momentum interval 2 < p D 0 T < 36 GeV/c, through the D 0 → K − π + (and charged conjugate) hadronic decay channel, which has a branching ratio of 3.95 ± 0.03% (Ref. [1]). The D 0 -meson candidates were identified by topological selections based on the displacement of the D 0 -meson candidate decay vertex, in addition to applying particle identification on the D 0 -meson candidate decay particles. These selection criteria largely suppress the combinatorial background of K ∓ π ± pairs that do not originate from the decay of a D 0 meson. Further details on the selection criteria are provided in Ref. [23].
Tracks (reconstructed charged-particle trajectories) corresponding to the D 0 -meson candidate decay particles were replaced by the reconstructed D 0 -meson candidate in the event, with the D 0 -meson candidate four-momentum being the sum of the decay-particle four-momenta. One benefit of this procedure is to avoid the case in which the decay products of the D 0 -meson candidate fill the dead-cone region. A jet-finding algorithm was then used to cluster the particles (tracks and the D 0 -meson candidate) in the event, to reconstruct the parton shower by sequentially recombining the shower particles into a single object (the jet). The jet containing the D 0 -meson candidate was then selected. The four-momentum of the jet is a proxy for the four-momentum of the charm quark initiating the parton shower. The jet-finding algorithm used was the anti-k T algorithm [24] from the Fastjet package [25], which is a standard choice for jet reconstruction because of its high performance in reconstructing the original parton kinematics. More details on the jet finding procedure can be found in the Methods.

Reconstructing the jet shower
Once jets containing a D 0 -meson candidate amongst their constituents are selected, the internal cascade process is reconstructed. This is done by reorganizing (reclustering) the jet constituents according to the Cambridge-Aachen (C/A) algorithm [26], which clusters these constituents based solely on their angular distance from one another. A pictorial representation of this reclustering process, which starts by reconstructing the smallest angle splittings, is shown in the top panels of Fig. 1. As QCD emissions approximately follow an angular-ordered structure [27], the C/A algorithm was chosen as it also returns an angular-ordered splitting tree.
This splitting tree is then iteratively declustered by unwinding the reclustering history, to access the building blocks of the reconstructed jet shower. At each declustering step, two prongs corresponding to a splitting are returned. The angle between these splitting daughter prongs, θ , the relative transverse Figure 1: Reconstruction of the showering quark. A sketch detailing the reconstruction of the showering charm quark, using iterative declustering, is presented. The top panels show the initial reclustering procedure with the C/A algorithm, in which the particles separated by the smallest angles are brought together first. Once the reclustering is complete, the declustering procedure is carried out by unwinding the reclustering history. Each splitting node is numbered according to the declustering step in which it is reconstructed. With each splitting, the charm-quark energy, E Radiator,n , is reduced and the gluon is emitted at a smaller angle, θ n , with respect to previous emissions. The mass of the heavy quark, m Q , remains constant throughout the showering process. At each splitting, gluon emissions are suppressed in the dead-cone region (shown by a red cone for the last splitting), which increases in angle as the quark energy decreases throughout the shower. momentum of the splitting, k T , and the sum of the energy of the two prongs, E Radiator , are registered. As the charm flavour is conserved throughout the showering process, the full reconstruction of the D 0 -meson candidate enables the isolation of the emissions of the charm quark in the parton shower, by following the daughter prong containing the fully reconstructed D 0 -meson candidate at each declustering step. This can be seen in the bottom part of Fig. 1, which shows the evolution of the charm quark reconstructed from the measured final state particles. Moreover, the kinematic properties of the charm quark are updated along the splitting tree, enabling an accurate reconstruction of each emission angle against the dynamically evolving charm-quark direction. It was verified that in more than 99% of the cases the prong containing the D 0 -meson candidate at each splitting coincided with the leading prong. This means that following the D 0 -meson candidate or leading prong at each step is equivalent, and therefore a complementary measurement for an inclusive jet sample, when no flavour tagging is available, can be made by following the leading prong through the reclustering history. As the inclusive sample is dominated by massless gluon and nearly massless light quark-initiated jets, it acts as a reference to highlight the mass effects present in the charm-tagged sample.

Extracting the true charm splittings
The selected sample of splittings has contributions from jets tagged with combinatorial K ∓ π ± pairs, which are not rejected by the applied topological and particle identification selections. The measured invariant mass of real D 0 mesons, which corresponds to the rest mass, is distributed in a Gaussian (because of uncertainties in the measurement of the momenta of the K ∓ π ± pairs) with a peak at the true D 0 -meson mass. This enables the implementation of a statistical two-dimensional side-band subtraction procedure, which characterizes the background distribution of splittings by sampling the background-dominated re-gions of the D 0 -meson candidate invariant mass distributions, far away from the signal peak. In this way the combinatorial contribution can be accounted for and removed. Furthermore, the selections on the D 0meson candidates also select a fraction of D 0 mesons originating as a product of beauty-hadron decays. These were found to contribute 10-15% of the reconstructed splittings, with only a small influence on the results which will be discussed later. The studies were performed using Monte-Carlo (MC) PYTHIA 6.425 (Perugia 2011) [28,29] simulations (this generator includes mass effects in the parton shower [30] and was used for all MC-based corrections in this work), propagating the generated particles through a detailed description of the ALICE detector with GEANT3 (Ref. [31]). The finite efficiency of selecting real D 0 -meson tagged jets, through the chosen selection criteria on the D 0 -meson candidates, as well as kinematic selections on the jets, was studied and accounted for through MC simulations. This efficiency was found to be strongly p D 0 T dependent and different for D 0 mesons originating from the hadronization of charm quarks or from the decay of beauty hadrons. Further details on these analysis steps can be found in the Methods.
As the reconstructed jet shower is built from experimentally detectable hadrons, as opposed to partons, hadronization effects must be accounted for. As hadronization processes occur at low non-perturbative scales, they are expected to distort the parton shower by mainly adding low-k T splittings [32]. A selection of k T > 200 MeV/c ensures that only sufficiently hard splittings are accepted and is used to suppress such hadronization effects. Other choices of k T selection were also explored, with stronger k T selections further removing non-perturbative effects from the measurement, at the expense of statistical precision. Other non-perturbative effects, such as the underlying event, contribute with extra soft splittings primarily at large angles and do not affect the small-angle region under study.
Detector effects also distort the reconstructed parton shower through inefficiencies and irresolution in the tracking of charged particles. However, these have been tested and largely cancel in the final observable and any residual effects are quantified in a data-driven way and included in the systematic uncertainties.
It should be noted that in addition to direct heavy-flavour pair creation in the elementary hard scattering, charm quarks can also be produced in higher-order processes as a result of gluon splitting. Therefore, the shower history of D 0 mesons containing such charm quarks will also have contributions from splittings originating from gluons. Furthermore, in the case of high transverse momentum gluons in which the charm quarks are produced close in angle to each other, the dead-cone region of the charm quark hadronizing into the reconstructed D 0 meson can be populated by particles produced in the shower, hadronization and subsequent decays of the other (anti-)charm quark. The influence of such contaminations through gluon splittings was studied with MC simulations and found to be negligible.

The observable R(θ )
The observable used to reveal the dead cone is built by constructing the ratio of the splitting angle (θ ) distributions for D 0 -meson tagged jets and inclusive jets, in bins of E Radiator . This is given by, where the θ distributions were normalized to the number of jets that contain at least one splitting in the given E Radiator and k T selection, denoted by N D 0 jets and N inclusive jets for the D 0 -meson tagged and inclusive jet samples, respectively. Expressing equation (1) in terms of the logarithm of the inverse of the angle is natural, given that at leading order the QCD probability for a parton to split is proportional to ln(1/θ ) ln(k T ).
A selection on the transverse momentum of the leading track in the leading prong of each registered splitting in the inclusive jet sample, p ch,leadingtrack T,inclusivejets ≥ 2.8 GeV/c, was applied. This corresponds to the transverse mass (obtained through the quadrature sum of the rest mass and transverse momentum) of a 2 GeV/c D 0 -meson and accounts for the p D 0 T selection in the D 0 -meson tagged jet sample, enabling a fair comparison of the two samples.
In the absence of mass effects, the charm quark is expected to have the same radiating properties as a light quark. In this limit, equation (1) can be rewritten as where the superscript LQ refers to light quarks, and the inclusive sample contains both light-quark and gluon-initiated jets. This indicates that the R(θ ) no dead−cone limit ratio depends on the differences between light-quark and gluon radiation patterns, which originate from the fact that gluons carry two colour charges (the charge responsible for strong interactions) whereas quarks only carry one. These differences result in quarks fragmenting at a lower rate and more collinearly than gluons. Therefore, in the limit of having no dead-cone effect, the ratio of the θ distributions for D 0 -meson tagged jets and inclusive jets becomes R(θ ) no dead−cone limit > 1, at small angles. This was verified through SHERPA v.2.2.8 (Ref. [33]) and PYTHIA v.8.230 (Tune 4C) [34] MC generator calculations, with the specific R(θ ) no dead−cone limit value dependent on the quark and gluon fractions in the inclusive sample. SHERPA and PYTHIA are two MC generators commonly used in high-energy particle physics and they use different shower prescriptions and hadronization models. Both models implement the dead-cone effect.

Exposing the dead cone
The measurements of R(θ ), in the three radiator (charm-quark) energy intervals 5 < E Radiator < 10 GeV, 10 < E Radiator < 20 GeV and 20 < E Radiator < 35 GeV, are presented in Fig. 2. Detector effects largely cancel out in the ratio and results are compared to particle-level simulations. Residual detector effects are considered in the systematic uncertainty together with uncertainties associated with the reconstruction and signal extraction of D 0 -meson tagged jets, as well as detector inefficiencies in the reconstruction of charged tracks in both the D 0 -meson tagged and inclusive jet samples. More details on the study of systematic uncertainties can be found in the Methods.
A significant suppression in the rate of small-angle splittings is observed in D 0 -meson tagged jets relative to the inclusive jet population. In Fig. 2, the data are compared with particle-level SHERPA (green) and PYTHIA v.8.230 (blue) MC calculations, with SHERPA v.2.2.8 providing a better agreement with the data. The no dead-cone baseline, as described in equation (2), is also provided for each MC generator (dashed lines). The suppression of the measured data points relative to the no dead-cone limit directly reveals the dead cone within which the charm-quark emissions are suppressed. The coloured regions in the plots correspond to the dead-cone angles in each E Radiator interval, θ dc < m Q /E Radiator , where emissions are suppressed. For a charm-quark mass m Q = 1.275 GeV/c 2 (Ref. [1]), these angles correspond to ln(1/θ dc ) ≥ 1.37, 2 and 2.75 for the intervals 5 < E Radiator < 10 GeV, 10 < E Radiator < 20 GeV and 20 < E Radiator < 35 GeV, respectively. These values are in qualitative agreement with the angles at which the data starts to show suppression relative to the MC limits for no dead-cone effect. The magnitude of this suppression increases with decreasing radiator energy, as expected from the inverse dependence of the dead-cone angle on the energy of the radiator.
A lower limit for the significance of the small-angle suppression is estimated by comparing the measured data to R(θ ) = 1, which represents the limit of no dead-cone effect in the case in which the inclusive sample is entirely composed of light quark-initiated jets. To test the compatibility of the measured data with the R(θ ) = 1 limit, a statistical test was performed by generating pseudodata distributions consistent with the statistical and systematic uncertainties of the measured data. A chi-square test was then carried out against this hypothesis for each of the pseudodata distributions. The mean P values correspond to significances of 7.7σ , 3.5σ and 1.0σ , for the 5 < E Radiator < 10 GeV, 10 < E Radiator < 20 GeV and 20 < E Radiator < 35 GeV intervals, respectively. A σ value greater than 5 is considered the criteria for a definitive observation, whereas the value of 1.0 is consistent with the null hypothesis.
The MC distributions shown were generated separately for prompt (charm-quark initiated) and nonprompt (beauty-quark initiated) D 0 -meson tagged jet production and were then combined using the prompt and non-prompt fractions in data calculated with POWHEG [35]+PYTHIA 6.425 [34] simulations. The non-prompt fraction was found to be independent of the splitting angle and corresponds to approximately 10% of the splittings in the 5 < E Radiator < 10 GeV interval and approximately 15% of the splittings in both the 10 < E Radiator < 20 GeV and 20 < E Radiator < 35 GeV intervals. It was verified through the MC simulations that non-prompt D 0 -meson tagged jets should exhibit a smaller suppression at small angles in R(θ ) compared with inclusive jets than their prompt counterparts. This is due to the additional decay products accompanying non-prompt D 0 -meson tagged jets that are produced in the decay of the beauty hadron. These may populate the dead-cone region, leading to a smaller observed suppression in R(θ ), despite the larger dead-cone angle of the heavier beauty quark.

Conclusions
We have reported the direct measurement of the QCD dead cone, using iterative declustering of jets tagged with a fully reconstructed charmed hadron. The dead cone is a fundamental phenomenon in QCD, dictated by the non-zero quark masses, whose direct experimental observation has previously remained elusive. This measurement provides insight into the influence of mass effects on jet properties and provides constraints for MC models. These results pave the way for a study of the mass dependence of the dead-cone effect, by measuring the dead cone of beauty jets tagged with a reconstructed beauty hadron.
A future study of the dead-cone effect in heavy-ion collisions, in which partons interact strongly with the hot QCD medium that is formed and undergo energy loss through (dominantly) medium-induced radiation, is also envisaged. If a dead cone were observed for these medium-induced emissions, it would be a confirmation of the theoretical understanding of in-medium QCD radiation, which is a primary tool used to characterize the high-temperature phase of QCD matter [36][37][38].
The quark masses are fundamental constants of the standard model of particle physics and needed for all numerical calculations within its framework. Because of confinement, their values are commonly inferred through their influence on hadronic observables. An exception is the top quark, which decays before it can hadronize, as its mass can be constrained experimentally from the direct reconstruction of the decay final states [39] (see Ref. [40] for a review of top mass measurements at the Fermilab Tevatron and CERN LHC).
By accessing the kinematics of the showering charm quark, before hadronization, and directly uncovering the QCD dead-cone effect, our measurement provides direct sensitivity to the mass of quasi-free charm quarks, before they bind into hadrons.
Furthermore, future high-precision measurements using this technique on charm and beauty-tagged jets, potentially in conjunction with machine-learning tools to separate quark and gluon emissions, could experimentally constrain the magnitude of the quark masses.

Acknowledgements
The ALICE Collaboration is grateful to Urs Wiedemann for his valuable suggestions and fruitful discussions. We are also grateful to Davide Napoletano for providing the SHERPA configuration files and useful discussions.
The ALICE Collaboration would like to thank all its engineers and technicians for their invaluable contributions to the construction of the experiment and the CERN accelerator teams for the outstanding performance of the LHC complex. The ALICE Collaboration gratefully acknowledges the resources and support provided by all Grid centres and the Worldwide LHC Computing Grid (WLCG) collaboration. The ALICE Collaboration acknowledges the following funding agencies for their support in building and running the ALICE detector:

A.1 Detector setup and data set
The analysis was performed with the ALICE detector at the CERN LHC [41]. The ALICE Inner Tracking System [42] and Time Projection Chamber [43] were used for charged-particle reconstruction, whilst particle identification (PID) was obtained using the combined information from the Time Projection Chamber and the Time-Of-Flight detectors [44]. These detectors are located in the ALICE central barrel, which has full azimuthal coverage and a pseudorapidity range of |η| < 0.9. The data set used in this analysis was collected in 2016, 2017, and 2018 in proton-proton collisions at √ s = 13 TeV, with a minimum-bias trigger condition defined by the presence of at least one hit in each of the two V0 scintillators [45]. This trigger accepts all events of interest for this analysis and the collected data sample corresponds to an integrated luminosity of L int = 25 nb −1 .

A.2 Jet finding and tagging
Jet finding was performed using the anti-k T algorithm, with a jet resolution parameter of R = 0.4. The E-scheme recombination strategy was chosen to combine the tracks of the jet by adding their fourmomenta, with a geometric constraint on the pseudorapidity of |η| < 0.5 enforced on the jet axis, to ensure that the full jet cone was contained in the acceptance of the central barrel of the ALICE detector. The ALICE detector has excellent tracking efficiency down to low p T (≈ 80% at p T = 500 MeV/c) which is homogeneous as a function of pseudorapidity and azimuthal angle [41], within the acceptance. The effect of track density on the tracking efficiency is also negligible [46]. The angular resolution is about 20% down to splitting angles of 0.05 radians, which motivated a track-based jet measurement as opposed to a full jet measurement using calorimetric information. Recent measurements [15,20] have shown that track-based jet observables are successful at reconstructing the parton shower information through declustering techniques, despite missing the information from the neutral component of the jet.
Jets with a transverse momentum in the interval of 5 ≤ p ch T,jet < 50 GeV/c were selected for this analysis. To mitigate against the cases in which two D 0 -meson candidates share a common decay track, jet-finding passes were performed independently for each D 0 -meson candidate in the event, each time replacing only the decay tracks of that candidate with the corresponding D 0 -meson candidate. In each pass the jet containing the reconstructed D 0 -meson candidate of that pass was subsequently tagged as a charminitiated jet candidate.
A.3 Subtraction of the combinatorial background in the D 0 -meson candidate sample To extract the true D 0 -meson tagged jet R(θ ) distributions and remove the contribution from combinatorial K ∓ π ± pairs surviving the topological and PID selections, a side-band subtraction procedure was used. This involved dividing the sample in p D 0 T intervals and fitting the invariant-mass distributions of the D 0 candidates in each interval with a Gaussian function for the signal and an exponential function for the background. The width (σ ) and mean of the fitted Gaussian were used to define signal and side-band regions, with the two-dimensional distributions of θ and E Radiator for D 0 -meson tagged jet candidates, ρ(θ , E Radiator ) D 0 jet candidate , obtained in each region. The signal region was defined to be within 2σ on either side of the Gaussian mean and contained most of the real D 0 mesons, with some contamination present from the combinatorial background. The side-band regions were defined to be from 4σ to 9σ away from the peak in either direction and were composed entirely of background D 0 -meson candidates. The combined ρ(θ , E Radiator ) D 0 jet candidate distributions measured in the two side-band regions represent the structural form of the contribution of background candidates to the ρ(θ , E Radiator ) D 0 jet candidate distribution measured in the signal region. In this way, the background component of the total ρ(θ , E Radiator ) D 0 jet candidate measured in the signal region can be subtracted, using the following equation, where the subscripts S and B denote the signal and side-band regions of the invariant-mass distributions, respectively. The A S and A B variables are the areas under the background fit function in the signal and combined side-band regions, respectively, and were used to normalize the magnitude of the background in the side-band regions to that in the signal region. The D 0 -meson tagged jet selection efficiency (discussed in more detail in the next section) is denoted by ε, with the index i running over the p D 0 T bins. As a result of this side-band subtraction, the true D 0 -meson tagged jet ρ(θ , E Radiator ) D 0 jet distributions are obtained, in the different intervals of p D 0 T .

A.4 D 0 -meson tagged jet reconstruction efficiency correction
The topological and PID selections used to identify the D 0 mesons, in the chosen jet kinematic interval, have a limited efficiency, which exhibits a strong p T dependence. Therefore, before integrating the sideband subtracted ρ(θ , E Radiator ) D 0 jet distributions across the measured p D 0 T intervals, the ρ(θ , E Radiator ) D 0 jet distributions were corrected for this efficiency. The efficiency, ε, was estimated from PYTHIA v.6 MC studies and varies strongly with p D 0 T , from approximately 0.01 at p D 0 T = 2.5 GeV/c to approximately 0.3 at p D 0 T = 30 GeV/c for prompt D 0 -meson tagged jets and from approximately 0.01 at p D 0 T = 2.5 GeV/c to approximately 0.2 at p D 0 T = 30 GeV/c for non-prompt D 0 -meson tagged jets. As the prompt and nonprompt D 0 -meson tagged jet reconstruction efficiencies were different, the final efficiency was obtained by combining the prompt and non-prompt D 0 -tagged jet reconstruction efficiencies, evaluated separately. These were combined with weights derived from simulations, corresponding to the admixture of prompt and non-prompt D 0 -meson tagged jets in the reconstructed sample. The fractions of this admixture were obtained in bins of p D 0 T by calculating the prompt and non-prompt D 0 -meson tagged jet production cross sections with POWHEG combined with PYTHIA v.6 showering.

A.5 Evaluation of systematic uncertainties
Considered sources of systematic uncertainty in the measurement relate to the reconstruction and signal extraction of D 0 -meson candidates, with the former contributing as the leading source. These uncertainties were estimated by varying the topological and PID selections, as well as the fitting and side-band subtraction configurations applied to the D 0 -meson candidate invariant mass distributions. Variations were chosen that tested the influence of selected analysis parameters as much as possible, while maintaining a reasonable significance in the signal extraction. For each of these categories, the root mean square of all deviations was taken as the final systematic uncertainty. Theoretical uncertainties in the prompt and non-prompt D 0 -meson tagged jet production cross sections from POWHEG were also considered in the calculation of the reconstruction efficiency, with the largest variation taken as the uncertainty. For each category, the final systematic uncertainty was symmetrized before adding up the uncertainties in quadrature across all categories to obtain the total systematic uncertainty of the D 0 -meson tagged jet measurement.
For the inclusive-jet results, the minimum p T requirement on the track with the highest transverse momentum within the leading prong of each splitting was varied. The magnitude of the variation was taken to be the resolution of the transverse momentum of a D 0 -meson with p D 0 T = 2 GeV/c, which was found to be 0.06 GeV/c. Variations above and below the nominal selection value were made and the largest deviation was symmetrized. Systematic detector effects are dominated by the tracking efficiency and were shown in detector simulations to affect both the D 0 -meson tagged jet and inclusive-jet samples equally, and they largely cancelled in the R(θ ) ratio. Therefore, the systematic uncertainty of R(θ ) because of detector effects was estimated directly on the ratio by randomly removing 15% of the reconstructed tracks, as given by the tracking efficiency of the ALICE detector, in the track samples used for clustering both the D 0 -tagged jets and inclusive jets. The ratio of the resulting R(θ ) distribution to the case with no track removal was taken, to obtain the uncertainty, which was symmetrized.
The relative uncertainty of R(θ ) resulting from the separate D 0 -meson tagged jet and inclusive-jet uncertainties was calculated, with the resulting absolute uncertainty added in quadrature to the detector effects uncertainty to obtain the total systematic uncertainty of the R(θ ) measurement. The magnitude of each of these sources of systematic uncertainty is shown in Table A.1, for the smallest-angle splittings corresponding to the interval 2 ≤ ln(1/θ ) < 3, in which the uncertainties are largest.