$Z^0$ boson associated b-jet production in high-energy nuclear collisions

The production of vector boson tagged heavy quark jets potentially provides new tools to probe the jet quenching effect. In this paper, we present the first theoretical study on the angular correlations ($\Delta\phi_{bZ}$), transverse momentum imbalance ($x_{bZ}$), and nuclear modification factor ($I_{AA}$) of $Z^0$ boson tagged b-jets in heavy-ion collisions, which was performed using a Monte Carlo transport model. We find that the medium modification of the $\Delta\phi_{bZ}$ for $Z^0\,+\,$b-jet has a weaker dependence on $\Delta\phi_{bZ}$ than that for $Z^0\,+\,$jet, and the modification patterns are sensitive to the initial jet $p_T$ distribution. Additionally, with the high purity of the quark jet in $Z^0\,+\,$(b-)jet production, we calculate the momentum imbalance $x_{bZ}$ and the nuclear modification factor $I_{AA}$ of $Z^0\,+\,$b-jet in Pb+Pb collisions. We observe a smaller $\Delta \left\langle x_{jZ} \right\rangle$ and larger $I_{AA}$ of $Z^0\,+\,$b-jet in Pb+Pb collisions relative to those of $Z^0\,+\,$jet, which may be an indication of the mass effect of jet quenching and can be tested in future measurements.


I. INTRODUCTION
The formation of quark-gluon plasma (QGP), which was produced in the early stages of the high-energy nucleus-nucleus collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC), offers a new possibility to test quantum chromodynamics (QCD) under an extremely hot and dense deconfined state of nuclear matter. The high-p T partons (quarks and gluons) produced in the initial hard scattering strongly interact with the QGP and dissipate their energy to the medium, which is referred to as the jet quenching effect [1][2][3][4][5]. Consequently, the "quenched jet" observables are used to quantify the properties [6] of the hot and dense QCD matter by investigating their medium modifications in heavy-ion collisions (HIC) relative to their p+p baselines.
In particular, new measurements of the associated production of the Z 0 boson and b-jet (denoted as Z 0 + b-jet) in p+p collisions at the LHC have been performed by ATLAS and CMS [39][40][41][42][43][44][45], since the final state b-jet associated with the Z 0 boson is the dominant background * bwzhang@mail.ccnu.edu.cn of the associated production of Higgs and Z 0 bosons (Z 0 + H → Z + bb) within the standard model (SM) [46] and can test many physics scenarios beyond the SM that predict new generation mechanism of b quarks and Z 0 bosons [41]. It is noted that in HIC, the Z 0 boson tagged b-jet (as the initial energy of the b quark is well gauged by the vector boson and thus its energy loss can be directly obtained) is particularly suitable for exploring the quenching of the heavy flavor jet [47]. The "dead-cone" effect [48] of heavy quarks in QGP may lead to a smaller energy loss compared to light flavors, which is known as the mass effect of jet quenching and has attracted intense investigations [49][50][51][52][53][54]. At the particle level, the latest measurements indicate that the yield of the B meson appears to be less suppressed than that of the D meson in nucleus-nucleus collisions [55][56][57]. However, no clear evidence was found at the full-jet level in the previous experimental measurements, such as the comparison of the R AA between the inclusive jet and b-jet [58,59] and the comparison of the p T imbalance (x J ) between inclusive dijets and bb dijets [60], except for some preliminary indication in the recent measurement implemented by the ATLAS collaboration [61]. The possible reasons for this can be manifold; e.g., the large contribution of gluoninitiated b-parton processes may thwart the attempts at solving the problem. Previous studies [21,47] have indicated that the dominant contribution of the Z 0 tagged jet is a quark-initiated jet, and the study of Z 0 + b-jet in HIC, especially its different medium modifications compared with that of Z 0 + jet, will provide a very useful tool to directly address the mass effect between the lightquark jet and massive bottom jet. Nevertheless, thus far, studies on the associated production of the b-jet and Z 0 boson in nucleus-nucleus collisions are lacking.
With this in mind, in this work, we present a Monte Carlo transport simulation including elastic (col-lisional) [62] and inelastic (radiative) [63][64][65][66] interactions of the energetic parton in the hot/dense QCD medium, while taking the next-to-leading order (NLO) plus parton shower (PS) generated initial hard parton spectrum [67] as input, to study the in-medium modification of the vector boson Z 0 tagged b-jets. This framework was employed to describe the heavy-flavor jet production of highenergy nuclear collisions in our previous studies [68][69][70][71][72][73]. We first present our numerical results for Z 0 + jet and compare them with the available experimental data to test the applicability of our model. Then, we calculate the angular correlations of Z 0 + b-jet in A+A collisions and demonstrate that the modifications of these correlations are sensitive to the initial b-jet p T distribution instead of the azimuthal angle. In contrast to the case of Z 0 + jet, the requirement of b-tagging excludes the contribution from multiple jets so that the azimuthal angular correlations of Z 0 + b-jet show distinct pattern modifications. With the high purity of the light-quark jet in Z 0 + jet events [21,47], we expect to address the mass dependence of the jet quenching effects between Z 0 + jet and Z 0 + b-jet.
The remainder of this paper is organized as follows. In Sec. II we present the productions of Z 0 + b-jet in p+p collisions calculated via the Monte Carlo event generator and comparisons with experimental data. In Sec. III, we discuss our treatments of the jet in-medium evolution in A+A collisions. In Sec. IV, we present the simulated results and discussions of the azimuthal angular correlation, transverse momentum imbalance, and nuclear modification factor of Z 0 + b-jet in HIC. Finally, Sec. V summarizes the study. Before we move into the study on Z 0 + b-jet production in HIC, we should address its production in p+p collisions. Fig. 1 shows a few processes [39] contributing to the associated production of a Z 0 boson and b-jet. In Figs. 1(a) and (b), an initial bottom quark from the parton distribution function (PDF) derived from the gluon distribution of one beam particle suffers the hard scattering and then turns into a b-jet and a Z 0 boson in the final state. In the bottom two diagrams, the bb pairs originate from the hard scattering and then turn into two b-jets associated with an emitted Z 0 boson in the final state.
In this work, we use the MC@NLO event generator SHERPA [67] to obtain the initial Z 0 + b-jet production in a p+p collision. The tree-level matrix elements are calculated using the internal modules Amegic [74] and Comix [75], and the one-loop virtual correction is calculated using the external program BlackHat [76]. The parton shower based on the Catani-Seymour [77] subtraction method is matched with NLO QCD matrix elements via the MC@NLO method [78]. The NLO PDF from NNPDF3.0 [79] with a 5-flavor scheme has been chosen in the calculations. FastJet [80] with an anti-k algorithm is used in the final-state jet reconstruction.
To compare our calculation results based on SHERPA with experimental data for p+p collisions, the same configurations implemented by the CMS collaboration [45] were used in our simulations. The Z 0 boson is reconstructed according to its decay channels Z 0 → e + e − and Z 0 → µ + µ − . The transverse momentum of the electron and muon candidates is required to be larger than 20 GeV. To exclude the barrel-endcap transition region, the electrons are selected within the pseudorapidity region |η| < 1.44 or 1.57 < |η| < 2.4, while muons are selected within |η| < 2.4. According to the requirement of the experiment, the events are considered only when the invariant mass of the electron or muon pairs lies in the region 70 < M ll < 111 GeV. The jets associated with the Z 0 boson are reconstructed by FastJet using the anti-k T algorithm with a cone size of R = 0.5. To reduce the contribution from the underlying event, the reconstructed jets must be in the pseudorapidity region |η jet | < 2.4 and have p T,jet > 30 GeV. The contribution from the underlying event is less than 5%, as estimated by CMS [45] because the production of softer jets is significantly suppressed by this requirement.
For events with at least one b-jet, we show the differential cross sections calculated via SHERPA as a function of the leading b-jet p T and Z 0 boson p T in Fig. 2. In the upper panel of Fig. 3, the azimuthal angular correlation between the leading b-jet and Z 0 boson (∆φ bZ = |φ b − φ Z |) is compared between the calculation and the CMS data. In the lower panel of Fig. 3, for events with two bjets (hereinafter Z 0 + 2 b-jets), we plot the differential cross sections as a function of the azimuthal angle between the two b-jets (∆φ bb = |φ b1 − φ b2 |). To quantify the deviation of our calculations from the experimental data points, we estimate the χ 2 /d.o.f of these

III. IN-MEDIUM JET EVOLUTION
In high-energy nuclear collisions, a droplet of an exotic state of nuclear matter, i.e., the QGP, is expected to be formed. The high p T partons produced in the hard scattering propagating in the QGP suffer both collisional and radiative energy loss as a result of the in-medium interaction. Numerous theoretical approaches [81][82][83][84][85][86][87][88] and Monte Claro models [89][90][91][92][93][94][95][96] have been developed in the last two decades to describe the heavy flavor meson production in HIC both at the RHIC and LHC. Among them, the Langevin transport equations have been employed effectively to describe the heavy quark evolution in the expanding QCD medium [62,68,69,[97][98][99][100]. Because complete treatment of the heavy quark jets propagating in the QGP usually needs a simultaneous description of space-time evolution for both light and heavy partons [52,87,101,102], we use the modified Langevin equations [62,68,69] to describe the propaga-tion of heavy quarks in QGP and take into account the in-medium energy loss of light partons [63-66, 103, 104]. Additionally, owing the lack of a unified theoretical approach that covers the full phase space of in-medium jet evolution in HIC, there are usually two methods to combine the vacuum parton shower with the medium-induced radiation [105,106]. The first one is to introduce medium modifications on the vacuum parton splitting at higher virtuality scale t ≥qτ f (whereq is the jet transport parameter and τ f is the formation length of the radiated gluon), as in the Q-PYTHIA [107], MATTER [108] and JEWEL [109] models. In this work, we employ the alternative treatment implemented in LBT [110] and MAR-TINI [111] to simulate the jet energy loss at low virtuality t ≤qτ f and high energy E q ⊥ (where q ⊥ represents the momentum exchange between the hard parton and the hot medium). According to this strategy, we take the p+p events produced by SHERPA with full vacuum parton shower as the input, sample their initial spatial positions using the MC-Glauber model [112], and then simulate the subsequent in-medium evolution.

A. Collisional energy loss
The movement of a heavy quark with large mass (M T ) propagating in the hot/dense nuclear matter and suffering a large number of random kicks from the medium can be modeled as a Brownian motion [97]. Hence a discrete Langevin equation can be utilized to describe the propagation of heavy quarks in the QCD medium [62,68,69,90]: where ∆t represents the timestep in the Monte Carlo simulation, and Γ(p) is the drag coefficient representing the dissipation effect and controlling the strength of quasielastic scattering. ξ(t) is the stochastic term that obeys a Gaussian probability distribution, and leads to The diffusion coefficient κ is related to the drag coefficient Γ by the fluctuation-dissipation relation [113]: where D s is the spacial diffusion coefficient. The last term − p g is the recoil momentum due to the mediuminduced gluon radiation, which is discussed in the following section. At each timestep, we boost partons to the local rest frame of the expanding medium to update the four-momentum and then boost them back to the laboratory frame to update the spatial position. Note that the procedure is done for T > T c , where T c = 165 MeV is the QCD transition temperature [114,115]. The spacetime evolution profile of the bulk medium in Pb+Pb collision is provided by the smoothed VISHNU [116] code. Even though the event-by-event fluctuation effects on the jet energy loss are small, the initial geometry fluctuation may be non-negligible for other observables such as particle collective flow (v n ) [117][118][119]. Meanwhile, the calculation for leading logarithmic accuracy at Hard-Thermal-Loop approximation [103,104] is employed in our framework to take into account the collisional energy loss of light quarks and gluons: where L represent the transport path of the partons along the propagating direction, α s is the strong coupling constant, C s is the quadratic Casimir in color representation, and µ D represents the Debye screening mass in the QCD medium. Note that Eq. (7) is only employed for light partons to consider their collisional energy loss because they cannot be treated as massive particles to evolve in the medium with the Langevin equations. During each timestep, the amount of collisional energy loss of a light parton can be calculated by integrating Eq. (7). Because medium-induced gluon radiation is the dominant energy loss mechanism for light partons, this treatment can be regarded as an effective approximation.

B. Medium induced gluon radiation
The inelastic scattering also plays an important role in the in-medium energy loss of energetic partons [120,121]. In our work, the Higher-Twist (HT) radiated gluon spectra [63][64][65][66] is implemented to simulate the mediuminduced gluon radiation when a parton propagates in the dense and hot QCD matter: where x and k ⊥ represent the energy fraction and the transverse momentum of the radiated gluon, respectively. α s is the strong coupling constant, which is fixed at α s = 0.3 in our calculations, C s is the quadratic Casimir in color representation, and P (x) is splitting function [122] for the splitting processes q → q + g and g → g + g (g → q +q process is negligible owing to its low probability [100]).
τ f is the radiated gluon formation time defined as τ f = 2Ex(1−x)/(k 2 ⊥ +x 2 m 2 ), and t−t i is the time interval between two instead of inelastic scattering. In the rigorous Higher-Twist calculations, the splitting function P (x) of heavy quarks should be mass dependent [65]. Nevertheless, as a simplified treatment in many transport models, such as LBT [93,110], QLBT [123], and BAMPS [124,125], it's convenient to uniformly write P (x) as a massless form for heavy and light quarks, while the dominant mass effect of the gluon radiation spectra of heavy quarks can be approximately presented as an overall fac- Additionally,q is the jet transport coefficient [126]: whereq 0 denotes the value ofq at the center of the bulk medium at the initial time τ 0 = 0.6 fm/c, and ρ QGP (τ, r) is the parton number density where the parton is probed.
To take into account the radial flow effect [127], the fourmomentum of parton p µ and the four flow velocity of the medium in the collision frame u µ act as a modification forq in an expanding nuclear medium. We still hold that q/T 3 = const. in the current framework as an effective approximation. More studies focusing on the puzzle of the temperature dependence ofq/T 3 can be found in the Refs. [6,92,[128][129][130]. The last term of Eq. (8) represents the dead-cone effect [51,65] which suppresses the gluon radiation of heavy quarks at a small angle (θ ∼ M/E) owing to their large mass. We also noticed the up-to-data development of the HT approach in Ref. [131] achieved by including both transverse and longitudinal momentum exchanges between hard partons and the QCD medium and will consider their results in our framework in the future. An imposed cut-off of the radiated gluon energy fraction x min = µ D E is taken to avoid the divergence near x → 0, where µ 2 D = 4πα s (1 + n f 6 )T 2 is the Debye screening mass induced by the QGP medium. Therefore, following the method introduced in [93], one can estimate the mean number of radiated gluons N (t, ∆t) during a timestep ∆t by integrating the phase space of x, k ⊥ , and t in Eq. (8): By assuming that the multiple gluon radiation is a Poisson process, we obtain the probability distribution of the radiation number P (n, t, ∆t) during a time step, as well as the total inelastic scattering probability P rad (t, ∆t): In our Monte Carlo simulation, during every time step, the P rad (t, ∆t) is first evaluated to determine whether the radiation occurs. If accepted, the Possion distribution function P (n, t, ∆t) is used for the sampling of the radiated gluon number. Finally, the four-momentum of the radiated gluon can be sampled according to the spectrum dN/dxdk 2 ⊥ expressed in Eq. (8). In Fig. 4, for a consistent comparison between our Monte Carlo simulation and the analytical calculation, we estimate the collisional and radiative energy losses of the gluons, light quarks, charm, and bottom in a static medium (T = 400 MeV).
Here we fix the parton energy (50 GeV) at each evolution time step and also restore the initial time t i in Eq. (8) to be 0 (same as the treatment in Ref. [91]) because their variations during the Monte Carlo simulation are not automatically included in the analytical calculation. We find that the MC results agree well with the analytical calculations. For the collisional energy loss, the ∆E coll of the gluon is 9/4 times that of the light quark owing to the large color factor, and in our framework, the ∆E coll of the heavy quark is comparable to that of the light quark. For the radiative energy loss, a clear mass hierarchy for different parton species can be found: ∆E g rad > ∆E q rad > ∆E c rad > ∆E b rad . For a long propagation time t = 4 fm, we find that the radiative energy loss dominates the total parton energy loss because of its quadratic dependence on the path length.
In general, there are two parameters in our framework that need to be determined; the jet transport coefficientq and the diffusion coefficient D s . We treatq and κ as two independent parameters to be constrained by experimental data. First, the value ofq is determined via a global extraction of the single hadron production in Pb+Pb collisions [132], in which q 0 = 1.2 GeV 2 /fm is obtained at the LHC energy. Afterq is fixed, we extract the best value D s (2πT ) ∼ 4 via a χ 2 fitting to the D meson R AA data [133,134], which is consistent with the results of D s (2πT ) = 3.7 ∼ 7 reported by the Lattice QCD [135].

IV. NUMERICAL RESULTS AND DISCUSSIONS
In this section, to estimate the medium modification of jet observables in nucleus-nucleus collisions, we use the p+p events provided by SHERPA as the input of our simulation within the hydrodynamic background to study the in-medium jet evolution. Before proceeding to the Z 0 tagged b-jet, we calculate the azimuthal angular correlation (∆φ jZ = |φ jet −φ Z |) and transverse momentum imbalance (x jZ = p jet T /p Z T ) of Z 0 + jet, as well as the nuclear modification factor R AA of the inclusive b-jet, and compare our theoretical results with the available experimental data. Then, we calculate the Z 0 + b-jet observables, including the azimuthal angular correlation between the Z 0 boson and b-jet (∆φ bZ = |φ b-jet − φ Z |), angle separation between the two Z-tagged b-jets (∆φ bb = |φ b1 −φ b2 |), transverse momentum (x bZ = p b-jet T /p Z T ), and nuclear modification factor I AA [21,35] defined as Here, N bin denotes the average number of binary nucleon-nucleon collisions among A+A collisions calculated via the Glauber model [112]. In Fig. 5, we show our calculated ∆φ jZ distributions both in p+p and 0 − 30% Pb+Pb collisions at √ s N N =5.02 TeV compared with the CMS experimental data [31]. Additionally, χ 2 /d.o.f of our model fitting to the CMS data both in p+p and Pb+Pb are presented. The same configurations in the jet reconstruction used by the CMS are employed. All final-state jets are reconstructed by FastJet using the anti-k T algorithm with R = 0.3 and require p jet T > 30 GeV. The selected Z 0 boson are reconstructed by the electron or muon pairs based on their decay channels (Z 0 → e + e − and Z 0 → µ + µ − ) and require p Z T > 60 GeV. Note that these distributions are normalized by the number of Z 0 events, and the transverse momentum imbalances are subjected to Gaussian smearing [31] to take into account the detector resolution effects. This reveals that the distribution of azimuthal angular correlation in Pb+Pb collisions suffers a suppression in a small ∆φ jZ region relative to the p+p baseline, which is consistent with the CMS measurement. However, in a large angle region (∆φ jZ ∼ π, where the Z 0 boson and jet are almost back-to-back), this suppression is not very apparent. The reason for this behavior has been discussed in detail [26,34]; i.e., the small ∆φ jZ region is dominated by the multiple jets processes and the large ∆φ jZ region is dominated by soft/collinear radiation. Usually, the jet energy of multiple-jet processes is relatively low and easier to be shifted below the jet selection threshold (30 GeV) because of parton energy loss [27]. We also notice two significant differences between the p+p and Pb+Pb CMS data near ∆φ jZ = π/2 and ∆φ jZ = 7π/8, and our results cannot fit these points well. However, similar abnormal behaviors of the CMS data are not found in the measurements of the γ−jet [136]; hence, we guess that they may be caused by the statistical fluctuations in the experiment. In Fig. 6, we compute x jZ distribution for Z 0 + jet in both p+p and 0 − 30% Pb+Pb collisions at √ s N N =5.02 TeV compared with the CMS data. The χ 2 /d.o.f are also presented in the figures indicating that our calculations are consistent with the experimental data, but the p+p baseline needs to be improved. Note here that selected Z 0 + jet pairs are required to be almost backto-back (∆φ jZ > 7π/8). Relative to the p+p baseline, in Pb+Pb collisions, we find that the x jZ distribution is shifted toward smaller values, exhibiting an enhancement at 0 < x jZ < 0.7 and suppression at 0.7 < x jZ < 2. As x jZ represents the transverse momentum imbalance of Z 0 and jet, that for each Z 0 + jet pair, the values of x jZ decrease owing to the jet energy loss and thus are shifted to a smaller x jZ observed in the final state. In Fig. 7, we investigate the nuclear modification factor R AA of the inclusive b-jet in Pb+Pb collisions at √ s N N = 2.76 TeV in comparison with the CMS measurements [58] to test our model calculations. The values of χ 2 /d.o.f are shown in each panel of Fig. 7 to quantify the deviation of our calculations from the CMS data for different centrality bins. We find that our calculations are essentially consistent with the experimental data, but we may pay attention to some discrepancies. First, our theoretical results show weak p T dependence but the experimental R AA appears to increase visibly with the jet   p T in the 10 − 30% and 30 − 50% centrality bins. Second, some data points (p T = 100 GeV at 0 − 10% and 10 − 30%, p T = 210 GeV at 30 − 50% and 50 − 100%) cannot be well described by our calculations. The good agreement between our model calculations and the data of the Z 0 + jet and inclusive b-jet makes it possible to study the medium modification of Z 0 + bjet in nuclear-nuclear collisions. In Fig. 8, we calculate the azimuthal angular correlation of the Z 0 boson and the b-jet in p+p and 0 − 10% Pb+Pb collisions at √ s N N =5.02 TeV. The b-jets associated with the Z 0 boson are reconstructed by FastJet using the anti-k T algorithm with a cone size of ∆R = 0.5, |η jet | < 2.4 and p jet T > 30 GeV for both p+p and Pb+Pb collisions. Note that these distributions are normalized by the initial Z 0 + b-jet event number (in p+p collision) to address the medium modification. We observe an overall suppression in Pb+Pb collisions relative to the p+p baseline. We show their ratio PbPb/pp in the middle panel of Fig. 8, and find that the suppression for the Z 0 + b-jet has a far weaker dependence on ∆φ bZ than that for the Z 0 + jet, which exhibits stronger suppression in a small ∆φ jZ region where multiple jets dominate. All the selected jets must first be b quark tagged, and this requirement significantly reduces the contribution from multiple-jet pro- cesses when we consider the azimuthal angular (∆φ bZ ) distribution.
To address the key factor that leads to the flat suppression on ∆φ bZ distribution, we estimate the averaged b-jet transverse momentum p T as a function of ∆φ bZ , which can be calculated by: p T (∆φ) = dσ dp T d∆φ p T dp T dσ dp T d∆φ dp T .
The decrease of the selected event number in A+A collisions results from the in-medium energy loss, which shifts lower p T jet below the kinematic selection cut. The initial p T distribution actually reflects the ∆φ bZ dependence of this shift. It turns out that the distribution of p T versus ∆φ bZ is nearly a constant value at 55 GeV as shown in the bottom panel of Fig. 8, which leads to the rather flat suppression on ∆φ bZ distribution. Here the band of p T distribution represents the statistical standard errors in the simulations.
As mentioned in Sec. II, the azimuthal angular separation ∆φ bb of the two b-jets tagged by Z 0 boson is also a useful observable to distinguish the contribution from subprocesses where Z 0 boson is emitted from one of the final state b quark or gluon splitting (g → bb) [41], as shown in diagrams (c) and (d) of Fig. 1. Note that these two categories of contributions corresponding to the cases that the two b-jets are almost back-to-back or collinear. What interests us is how the ∆φ bb distribution of these two categories of Z 0 + 2 b-jet would be modified in the QGP. As shown in the top panel of Fig. 9, we plot the ∆φ bb distributions both in p+p and 0 − 10% Pb+Pb collisions at √ s N N =5.02 TeV and also plot the ratio PbPb/pp in the middle panel. In the upper panel, we can find a kink at ∆φ bb ∼ π/5 and a peak at ∆φ bb ∼ π in Fig. 9, which present the two contributions of Z 0 + 2 bjets production: the two b-jets are almost back-to-back or collinear. The two-peaks distribution is similar to what we have observed in the angular correlations of bb dijets [68], and the peaks at smaller and larger ∆φ bb regions corresponding to the gluon splitting (GSP) processes and flavor creation (FCR) processes. In the middle panel, we observe an upward trend of the ratio from 0.5 to 0.7 as ∆φ bb increases. To figure out the ∆φ bb dependence of the ratio PbPb/pp, we estimate the initial jet p T of the leading and sub-leading b-jet in p+p collisions, as shown in the bottom panel of Fig. 9. We find that, for both the leading and subleading one, p T is increasing with ∆φ bb , which shows a similar trend with that of the ratio PbPb/pp versus ∆φ bb . It may indicate that the medium modification of ∆φ bb in Pb+Pb has a close connection with the initial b-jet p T distribution versus ∆φ bb in p+p collisions. The associated production of Z 0 + jet may shed new light on the mass dependence of the jet quenching effect in nuclear matter, owing to the high purity of lightquark-initiated jets. In Ref. [47], the contributions from light-quark-jets and gluon-jets in the Z 0 + jet production are approximately 70% and 30%, respectively. To verify this point, in Fig. 10, we estimate the gluon-jet fraction in four categories of jets in p+p collisions at 5.02 TeV: inclusive jet, inclusive b-jet, Z 0 tagged jet, and Z 0 tagged b-jet. The gluon-jet in the events can be identified by requiring that the gluon is the leading parton in the jets, which can be easily implemented in the FastJet program. We find that at p jet T ∼ 50 GeV, the gluon-jet fraction is approximately 50% for the inclusive jet and 30% for the Z 0 tagged jet. The Z 0 -tagging requirement considerably decreases the gluon jet contribution (by 40%), especially at a low p T . More importantly, for the inclusive b-jet and Z 0 tagged b-jet, the contributions from gluon-jet are significantly suppressed owing to the requirement of b-quark tagging and show almost equal values.
We plot the scaled x jZ distributions of Z 0 + jet in both p+p and 0 − 10% Pb+Pb collisions at √ s N N =5.02 TeV in Fig. 11(a) and with the same configuration, the x bZ distributions of Z 0 + b-jet in Fig. 11(b). All the selected jets (b-jets) must satisfy ∆φ jZ > 7π/8 (∆φ bZ > 7π/8) to guarantee that they are almost back-to-back with the Z 0 boson. We observe that the distributions of x J shift toward smaller values for both Z 0 + jet and Z 0 + b-jet in Pb+Pb collisions relative to their p+p baselines, owing to the energy loss of the tagged jets. To perform a more intuitive comparison between Z 0 + jet and Z 0 + b-jet, we show the difference of the x J distribution in Pb+Pb from   and negative values at 0.8 < x J < 1.6. Then we find that the absolute value of the difference (PbPb-pp) of Z 0 + jet is larger than that of Z 0 + b-jet. Furthermore, we estimate the shifting of the mean value of momentum imbalance ∆ x jZ = x jZ pp − x jZ P bP b for Z 0 + jet, and ∆ x bZ = x bZ pp − x bZ P bP b for Z 0 + b-jet, as shown in Tab. I.
Here, J denotes different processes. Note that the standard errors of x J in our simulations are also presented in Tab. I. We find that within the statistical uncertainties, ∆ x jZ (∼ 0.136 ± 0.0108) for Z 0 + jet is visibly larger than ∆ x bZ (∼ 0.092 ± 0.012) for Z 0 + b-jet, indicating stronger modifications to the light-quark jet compared with the b-jet. Note that x jZ (x bZ ) is defined by the ratio of jet (b-jet) p T to Z 0 boson p T . The shifting of x jZ (x bZ ) toward smaller values is directly related to the amount of jet (b-jet) energy loss; hence, ∆ x jZ > ∆ x bZ indicates that light quark jets lose more energy than b-jets.
The nuclear modification factor I AA is practically another good observable to address the mass hierarchy and flavor dependence of the jet quenching effect. Comparisons of I AA between Z 0 + jet and Z 0 + b-jet would provide more reliable evidence of the mass effect of jet quenching. For this purpose, we present the calculations of I AA of Z 0 + jet and Z 0 + b-jet in 0 − 10% Pb+Pb collisions at √ s N N =5.02 TeV as a function of jet p T within three p Z T windows in Fig. 12. First, we find different shapes of I AA in these three panels. In the left panel (40 < p Z T < 60 GeV) the two curves of I AA for Z 0 + jet and Z 0 + b-jet are flat, but in the right panel (80 < p Z T < 120 GeV) the curves show enhancement at p jet T < 80 GeV. This is because if we constrain 80 < p Z T < 120 GeV in the event selection, the cross section at p jet T < 80 GeV is steeper, falling with jet p T . Then, the jets shifting from higher p T to a lower p T because of the in-medium energy loss naturally lead to the relatively large I AA values in the lower p jet T region, even larger than one. Additionally, we observe smaller values of I AA for Z + jet relative to those for Z + b-jet in the three panels of Fig. 12, which indicates that the yield of Z 0 tagged light-quark jets suffers stronger suppression after traversing the QCD matter than that of Z 0 tagged b-jets. Note that I AA is directly related to the p T shifting of the tagged jets due to in-medium energy loss. The smaller value of I AA for Z 0 + jet indicates a larger jet p T shift compared with that of Z 0 + b-jet, suggesting that the energy loss of Z 0 tagged b-jets is smaller than that of Z 0 tagged light-quark jets. Our conclusions are consistent with the recent ATLAS measurements [61], which indicate that the R AA of the b-jet is higher than that of the inclusive jet. We hope that our complementary predictions can be tested by future measurements at the LHC, which may be helpful for solving the puzzle of the mass hierarchy of jet quenching.

V. SUMMARY
Vector boson-tagged heavy quark jets are promising new tools for studying the jet quenching effect. In this work, we present a Monte Carlo transport simulation, which takes into account the elastic and inelastic jet interactions within a hydrodynamic background, to study the in-medium modification of Z 0 tagged b-jets. The NLO+PS event generator SHERPA was been used to provide the p+p baseline of Z 0 + b-jet production, which agreed well with the CMS measurements. This framework has been proven to give good descriptions of medium modifications of ∆φ jZ and x jZ of Z 0 + jet, as well as the R AA of the inclusive b-jet, measured in Pb+Pb collisions at the LHC.
The angular correlation between the vector boson and heavy quark-tagged jets may be a new promising observable for studying the in-medium jet interaction. We present the first calculation of the azimuthal angular correlation ∆φ bZ of Z 0 + b-jet in both p+p and 0 − 10% Pb+Pb collisions at √ s N N = 5.02 TeV. We observe a flat suppression factor versus ∆φ bZ , in contrast to the case of Z 0 + jet, because the requirement of b-tagging excludes the contribution from multiple-jet processes. Additionally, we calculate the medium modification of the azimuthal angular correlation ∆φ bb in central Pb+Pb collisions at √ s N N = 5.02 TeV and observe stronger suppression in a smaller ∆φ bb region of the distribution relative to that at ∆φ bb ∼ π. By analyzing the p T dis-tribution of the tagged b-jets, we find that the medium modification pattern on ∆φ bZ (∆φ bb ) in Pb+Pb has a close connection with the initial p T distribution versus ∆φ bZ (∆φ bb ) in p+p collisions. These investigations may help us to understand the experimental measurements of jet angular correlations at the LHC in recent years, e.g., for the γ+jet and Z 0 +jet. Finally, we predict that the mass effect of jet quenching can be addressed by comparing the medium modifications of Z 0 + jet and Z 0 + b-jet. With the high purity of the quark jet in Z 0 + (b-)jet events, we estimate the medium modification of the transverse momentum imbalance x jZ (x bZ ) and the nuclear modification factor I AA for both Z 0 + jet and Z 0 + b-jet in Pb+Pb collisions. We find a larger shift of x jZ and stronger suppression for I AA of Z 0 + jet than for Z 0 + b-jet indicating that b-jets lose less energy than light quark jets. These predictions can be tested via future measurements at the LHC and may provide a key to solving the puzzle of the mass hierarchy of jet quenching.