Inclusive, prompt and non-prompt J/ψ production at midrapidity in p-Pb collisions at sNN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sqrt{s_{\mathrm{NN}}} $$\end{document} = 5.02 TeV

A measurement of inclusive, prompt, and non-prompt J/ψ production in p-Pb collisions at a nucleon-nucleon centre-of-mass energy sNN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sqrt{s_{\mathrm{NN}}} $$\end{document} = 5.02 TeV is presented. The inclusive J/ψ mesons are reconstructed in the dielectron decay channel at midrapidity down to a transverse momentum pT = 0. The inclusive J/ψ nuclear modification factor RpPb is calculated by comparing the new results in p-Pb collisions to a recently measured proton-proton reference at the same centre-of-mass energy. Non-prompt J/ψ mesons, which originate from the decay of beauty hadrons, are separated from promptly produced J/ψ on a statistical basis for pT larger than 1.0 GeV/c. These results are based on the data sample collected by the ALICE detector during the 2016 LHC p-Pb run, corresponding to an integrated luminosity L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{L} $$\end{document}int = 292 ± 11 μb−1, which is six times larger than the previous publications. The total uncertainty on the pT-integrated inclusive J/ψ and non-prompt J/ψ cross section are reduced by a factor 1.7 and 2.2, respectively. The measured cross sections and RpPb are compared with theoretical models that include various combinations of cold nuclear matter effects. From the non-prompt J/ψ production cross section, the bb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathrm{b}\overline{\mathrm{b}} $$\end{document} production cross section at midrapidity, dσbb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathrm{d}\sigma}_{\mathrm{b}\overline{\mathrm{b}}} $$\end{document}/dy, and the total cross section extrapolated over full phase space, σbb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\sigma}_{\mathrm{b}\overline{\mathrm{b}}} $$\end{document}, are derived.


JHEP06(2022)011
5.02 TeV, using the data sample collected in 2016, which is six times larger than that of 2013. Moreover, the cross section of inclusive J/ψ production measured in pp collisions at √ s = 5.02 TeV [7] is used to derive the R pA results instead of the interpolation procedure adopted in the previous p-Pb publication [57]. Therefore, the new results, which are significantly more precise and are obtained differentially in p T and in finer p T intervals, supersede the measurements published in refs. [51,57].

Data analysis
A complete description of the ALICE apparatus and its performance is presented in refs. [69,70]. The central-barrel detectors employed for the analysis presented in this paper are the Inner Tracking System (ITS) and the Time Projection Chamber (TPC). The ITS [71] provides tracking and vertex reconstruction close to the interaction point (IP). It is made up of six concentric cylindrical layers of silicon detectors surrounding the beam pipe with radial positions between 3.9 cm and 43.0 cm. The two innermost layers consist of Silicon Pixel Detectors (SPD), the two central layers are made up of Silicon Drift Detectors (SDD), and the two outermost layers of Silicon Strip Detectors (SSD). The TPC [72] consists of a large cylindrical drift chamber surrounding the ITS and extending from 85 cm to 247 cm along the radial direction and from −250 cm to +250 cm along the beam direction (z) relative to the IP. It is the main ALICE tracking device and allows also charged particles to be identified through specific energy loss (dE/dx) measurements in the detector gas. Both the TPC and ITS are embedded in a solenoidal magnet that generates a 0.5 T magnetic field along the beam direction. They cover the pseudorapidity interval |η| < 0.9 and allow J/ψ mesons to be reconstructed through the e + e − decay channel in the central rapidity region down to zero p T .
The measurements presented in this paper are based on the set of minimum bias (MB) p−Pb collisions at √ s NN = 5.02 TeV collected in 2016 during the LHC Run 2 data taking period, corresponding to an integrated luminosity L int = 292 ± 11 µb −1 . The latter is determined from the number of MB events and the MB-trigger cross section, which was measured via a van der Meer scan, with negligible statistical uncertainty and a systematic uncertainty of 3.7% [73]. Collisions were realized by delivering proton and Pb beams with energies of 4 TeV and 1.58 TeV per nucleon, respectively. The proton and Pb beams circulated in the LHC anticlockwise and clockwise, respectively, during the period of data taking considered for this analysis. The MB trigger condition is provided by the V0 detector [74]: a system made up of two arrays of plastic scintillators placed on either side of the IP and covering the full azimuthal angle and the pseudorapidity intervals 2.8 < η < 5.1 and −3.7 < η < −1.7. The trigger condition required at least one hit in both the two arrays during the nominal bunch crossing time frame, allowing non-singlediffractive p−Pb collisions to be selected with an efficiency higher than 99% [75]. The timing information from the V0 detectors is also used, in combination with that from the SPD, to implement an offline rejection of beam-induced background interactions occurring outside the nominal colliding bunch crossings. Events with more than one interaction per bunch crossing are reduced down to a negligible amount by means of a dedicated algorithm -3 -

JHEP06(2022)011
employing reconstructed tracks to detect the presence of multiple collision vertices. Only collision events with a reconstructed primary vertex lying within ±10 cm from the nominal IP along the beam direction are considered in order to obtain a uniform coverage for the central-barrel detectors. An event sample of about 6 × 10 8 MB events is obtained after the application of the above described selection criteria.

Inclusive J/ψ
Electron candidates are selected following similar procedures as those described in ref. [57]. The tracks, reconstructed with the ITS and TPC detectors, are required to have a transverse momentum p e T > 1.0 GeV/c and pseudorapidity |η e | < 0.9, as well as at least 70 (out of a maximum of 159) attached TPC clusters and a track fit χ 2 /dof < 2 in order to ensure a uniform tracking efficiency in the TPC. Electron identification is performed by requiring the measured dE/dx to be compatible with the expected specific energy loss for electrons within 3σ, with σ denoting the specific energy-loss resolution of the TPC. Tracks compatible with the pion and proton energy loss expectations within 3σ are rejected. At least one hit in either of the two SPD layers is required to remove background electrons produced from the conversion of photons in the detector materials at large radii. Additional suppression of this background is realized by discarding electron (positron) candidates, which are compatible with a photon conversion when combined with a positron (electron) candidate of the same event, through the application of dedicated topological selections. These selections were verified, employing Monte Carlo (MC) simulations, to have a negligible impact on the J/ψ signal. Finally, in order to reduce the overall background at low transverse momentum, a set of slightly tighter SPD and particle identification (PID) requirements is applied to electrons and positrons forming candidate pairs with p T < 3 GeV/c. A hit in the first SPD layer and a 3.5σ pion and proton rejection condition is required instead of 3σ for higher p T values. The sample of J/ψ candidates is obtained by combining the selected opposite-sign tracks in the same event and requiring the J/ψ rapidity to be within |y lab | < 0.9 in the laboratory system. Due to the energy asymmetry of the proton and lead beams, such a requirement corresponds to a selection of J/ψ candidates within −1.37 < y < 0.43 in the nucleon−nucleon centre-ofmass system. The resulting dielectron invariant mass (m e + e − ) distributions are shown in figure 1 for eight selected transverse momentum intervals, from 0 to 14 GeV/c. The signal component is characterised by an asymmetric shape, with a long tail towards low invariant masses due to the J/ψ radiative decay channel (J/ψ → e + e − γ ) and the bremsstrahlunginduced energy loss of daughter electrons in the detector material. The background component is composed of both a combinatorial and a correlated part, with the latter mainly originating from the semileptonic decay of correlated open heavy-flavour hadrons.
The inclusive J/ψ yield is determined from the invariant mass distributions using the same technique as described in ref. [7]. At first, the combinatorial background shape is modelled by means of a mixed event (ME) technique and then scaled to the invariant mass distribution of like-sign track pairs. Then, the combinatorial background is subtracted from the opposite-sign dielectron invariant mass distribution, and the correlated background is evaluated by fitting the resulting distribution with a two-component function composed of a MC template for the J/ψ signal and of an empirical function for the correlated back-  Figure 1. Opposite-sign dielectron invariant mass distributions for the p T -intervals used for this analysis. The signal plus total background (blue), the combinatorial background (red), and the correlated background (green), evaluated as described in the text, are shown separately in each panel. The χ 2 /ndf values of the signal template plus the total background function are also reported along with the raw yields in the range 2.92 < m e + e − < 3.16 GeV/c 2 . ground. The latter is defined to be either an exponential or a combination of an exponential and a polynomial. The former is obtained by a detailed MC simulation of J/ψ decays in the ALICE detectors, based on GEANT3 [76] and the full reconstruction chain as for real events, which is then also used to correct the raw yield for the selection procedure and detector inefficiencies and described in details in the next paragraph. After subtracting from the opposite-sign dielectron invariant mass distribution also the correlated background, the raw J/ψ yield is obtained by counting the number of entries within the invariant mass interval 2.92 < m e + e − < 3.16 GeV/c 2 . In figure 1, the different components used in the procedure to describe the opposite-sign dielectron invariant mass distributions are shown superimposed. An alternative method was also considered, where the invariant mass distribution after the subtraction of the correlated background is fitted with a Crystal Ball (CB) function [77] for the signal plus either an exponential or a combination of an exponential -5 -JHEP06(2022)011 and a polynomial for the background, and the raw yield is obtained from the integral of the best-fit CB function. The alternative method yields results compatible with those from the standard approach.
In order to correct the raw yield for the chosen selection procedure as well as for detector inefficiencies, a MC simulation was implemented by injecting J/ψ signal events into MB p−Pb collision events simulated with the EPOS-LHC model [78]. The J/ψ component was generated starting with p T and y distributions that match well a next-to-leading order (NLO) Colour Evaporation Model (CEM) calculations [79,80] with the inclusion of nuclear effects based on the EPS09 parameterisation [81]. The J/ψ decay into dielectrons was simulated using the EvtGen package [82] in combination with the PHOTOS model [83] in order to provide a proper description of the radiative decay channel. In the simulation, GEANT3 [76] was used to reproduce the propagation of particles through the ALICE experimental setup, taking into account the response of the detectors. The same reconstruction procedure used for data was then applied to the simulated events in order to evaluate the product of acceptance times efficiency (A× ), which accounts for: the detector acceptance, the track quality requirements, the electron identification criteria, and the fraction of the signal counted within the invariant mass interval 2.92 < m e + e − < 3.16 GeV/c 2 . The A × retrieved from MC exhibits a smooth and mild variation with the J/ψ p T , ranging from ∼8.5% to ∼16% in the p T range from 0 to 14 GeV/c. Consequently, the resulting A × correction factor, which is the average of A × over p T in a finite-size p T interval, shows a weak dependence on the p T shape assumed in the simulation for the J/ψ component. In the end, the final correction factors were computed by re-weighting the original MC distribution to best-fit the inclusive J/ψ spectrum already measured in p-Pb collisions at the same energy [57]. The p T -differential cross section for inclusive J/ψ production is calculated as where ∆y = 1.8 corresponds to the width of the analysed rapidity interval, ∆p T is the width of the considered p T interval, N J/ψ is the raw J/ψ yield in the interval, and BR J/ψ → e + e − = (5.97 ± 0.03)% is the branching ratio for J/ψ decaying into dielectrons [84]. The inclusive J/ψ nuclear modification factor is obtained, according to eq. (1.1), by dividing the p T -differential cross section by the reference cross section measured up to p T = 10 GeV/c in pp collisions at √ s = 5.02 TeV [7]. The rapidity shift of ∆y = 0.465 between the p-Pb and pp samples is expected to introduce a 1% effect on the R pPb , which is negligible with respect to the other uncertainties. An interpolation procedure, which is described in ref. [51], is adopted for the computation of the reference cross section in the last p T interval, 10 < p T < 14 GeV/c, that was not measured in pp collisions at this energy. The interpolated value of d 2 σ pp /dydp T for this interval amounts to 10 ± 2 nb/(GeV/c), where the quoted uncertainty refers to the total systematic uncertainty arising from the interpolation procedure and is uncorrelated with the uncertainties of the measured d 2 σ pp /dydp T for p T < 10 GeV/c.  Total (R pPb ) 12.9 11.5 11.9 12.5 13.4 13.9 19.7 23.7 Table 1. Summary of the systematic uncertainties, in percentage, of the inclusive J/ψ cross section d 2 σ J/ψ /dydp T and nuclear modification factor R pPb in different p T intervals. All contributions to the d 2 σ J/ψ /dydp T uncertainty are considered to be highly correlated over the p T bins, except that for the background subtraction which is considered as fully uncorrelated. The reported values for the measured σ pp reference cross section, which was determined up to p T = 10 GeV/c [7], are the total uncertainties, both of statistical and systematic origins. The uncertainty of σ pp for the interval 10 < p T < 14 GeV/c is also the total uncertainty from the interpolation procedure as discussed in the text.
The estimated systematic uncertainties affecting the inclusive J/ψ measurements are listed in table 1. The dominant sources of uncertainty are related to the tracking and electron identification procedures. The remaining contributions are related to the signal extraction procedure, the J/ψ input kinematic distributions used in the MC simulation, the dielectron decay channel branching ratio, and the integrated luminosity determination.
The uncertainty of the tracking procedure dominates at low p T values and is related to both the ITS-TPC matching efficiency and to the adopted track quality requirements. The first component is estimated by evaluating the discrepancy in the matching probability of TPC tracks to ITS hits between data and Monte Carlo [85]. The observed discrepancy is used to re-scale the tracking efficiency of electrons in MC simulations in order to evaluate the difference in the resulting number of reconstructed J/ψ candidates. The second component is assessed by employing several variations to the adopted track selection criteria and by computing the RMS of the corrected J/ψ yield distribution resulting after these variations. The sum in quadrature of the uncertainties related to both these components is taken as systematic uncertainty on the tracking procedure. The uncertainty related to the electron identification is estimated by evaluating the TPC electron PID response for a clean sample of topologically identified electrons from conversion processes in data and computing the difference with the corresponding quantity from MC simulations. This pertrack uncertainty is then propagated to the reconstructed J/ψ candidates with the use of -7 -JHEP06(2022)011 MC simulations. The resulting uncertainty on the J/ψ cross section increases up to ∼ 6% towards high p T values, where it is the largest uncertainty contribution. The systematic uncertainty related to the signal extraction procedure is due to both the background subtraction and the assumptions on the signal shape. It is estimated as the RMS of the yield distributions corresponding to variations of the mass interval used for the signal counting, the alternative parameterisations employed to fit the correlated background, and the alternative method using the CB function to fit the signal. The uncertainty on the signal shape ranges between ∼2% at low p T and ∼2.5% at high p T , whereas the uncertainty due to the background subtraction varies between ∼1% and ∼2% and is largest for the lowest p T intervals. The systematic uncertainty on the J/ψ p T distribution used as input for the computation of the efficiency corrections is determined by randomly varying, within one standard deviation contour, the parameters of a function fitted to the measured p T distribution in p−Pb collisions [57], taking into account their correlations. The functional form used for this fit is discussed in ref. [86] and very well describes the measured J/ψ p T distribution. The uncertainty of the integrated luminosity amounts to 3.7% and is determined from the visible p−Pb cross sections measured in van der Meer scans as detailed in ref. [73]. Both the uncertainty on the integrated luminosity and that on the branching ratio BR J/ψ → e + e − = (5.97 ± 0.03)% [84] constitute global uncertainties for the inclusive J/ψ cross section, fully correlated between all p T intervals. All the other discussed sources of uncertainty are considered to be highly correlated 1 over p T , with the exception of the background uncertainty, which is considered as uncorrelated.
The total relative uncertainty for the reference pp cross section, σ pp , is also reported in table 1. In this case, the values up to p T = 10 GeV/c, are the total uncertainties, of both statistical and systematic origin, associated to the measurement performed up to p T = 10 GeV/c [7], while that for the 10 < p T < 14 GeV/c interval is the relative uncertainty of the interpolated cross section (10 ± 2 nb/(GeV/c)) quoted before. The uncertainty of σ pp propagates only to the R pPb observable and it varies between ∼10% and ∼20% and is largest for the highest p T intervals.

Determination of the non-prompt J/ψ fraction
The fraction f b of the J/ψ yield originating from b-hadron decays is measured for p T > 1 GeV/c by discriminating, on a statistical basis, the reconstructed J/ψ candidates according to the displacement between their production vertex and the primary p−Pb collision vertex. The discrimination is realized by means of an unbinned two-dimensional likelihood fit, following the same technique adopted in previous analyses for the pp [13], p−Pb [51], and Pb-Pb [87] systems. In particular, it is performed by maximising the following loglikelihood function in which N indicates the number of e + e − pairs within the 2.32 < m e + e − < 4.00 GeV/c 2 invariant mass interval. The pseudoproper decay length x is introduced to separate J/ψ 1 With high correlation, we mean a Pearson coefficient larger than 0.7.

JHEP06(2022)011
originating from the decay of b-hadrons from prompt J/ψ. It is defined as where c is the speed of light, L xy = L · p T /p T is the signed projection of the pair flight distance, L, onto its transverse momentum vector, p T , and m J/ψ is the J/ψ pole mass value [84]. The terms F Sig (F Bkg ) and M Sig (M Bkg ) in eq. (2.2) represent the probability density functions (PrDFs) describing the signal (background) pair distributions as a function of x and m e + e − , respectively, whereas f Sig denotes the ratio of signal to all candidates within the considered mass interval. The signal x PrDF is given by indicating the prompt and non-prompt J/ψ PrDFs, and f b being the uncorrected fraction of J/ψ coming from b-hadron decays. The evaluation of the different PrDFs used in eq. (2.2) is performed relying either on data or on MC simulations and following the same procedures described in previous analyses [13,87]. For the MC simulations, the prompt J/ψ component was generated with the p T and y distributions obtained with a procedure analogue to that previously discussed for the inclusive J/ψ analysis, while the non-prompt J/ψ component was obtained using PYTHIA 6.4 [88] with Perugia-0 tuning [89] to simulate the production of beauty hadrons. Also in this case, the J/ψ decay into dielectrons was simulated using the EvtGen package [82] in combination with the PHOTOS model [83] in order to provide a proper description of the radiative decay channel. The background x PrDF, F Bkg (x), is determined in three invariant mass ranges by fitting the x distributions of dielectron candidates in the lower (2.32 < m e + e − < 2.68 GeV/c 2 ) and upper (3.20 < m e + e − < 4.00 GeV/c 2 ) side bands of the invariant mass distributions and by interpolating the resulting fit functions to the region under the invariant mass signal peak (2.68 < m e + e − < 3.20 GeV/c 2 ). The experimental resolution function, R(x), which is the key ingredient in the F prompt (x), F b (x) and F Bkg (x) PrDFs, is evaluated from the x distributions of prompt J/ψ in MC simulations, reconstructed after applying the same selection criteria as in data. In order to improve the resolution of the secondary decay vertices, it is required that at least one of the two J/ψ candidate decay tracks has a hit in the innermost SPD layer. A tune-on-data procedure [51] is applied to the MC sample in order to reproduce the observed single-track impact parameter distributions. This minimises the discrepancy between data and simulation, reducing the systematic uncertainty related to the R(x) determination. The F b (x) PrDF is obtained as the convolution of the R(x) function and a template of the x distribution for the mixture of b-hadrons decaying into J/ψ. The latter is obtained with a MC simulation study of the kinematics of the b-hadron decays. In this simulation, the p T -distribution of the b-hadrons is obtained from pQCD calculations at fixed order with next-to leading-log re-summation (FONLL) [90]. The decay description is based on the EvtGen package [82], and the relative abundance of b-hadron species as a function of p T is based on the precise measurements reported by the LHCb collaboration in pp collisions [91], which are consistent with those measured in p-Pb collisions [67]. The x resolution estimated from the MC simulations is characterised by a pronounced dependence as a function of the J/ψ p T : for events with both J/ψ decay tracks yielding a hit in the first SPD layer, the RMS of the R(x) distribution ranges from ∼140 µm at p T = 1.5 GeV/c to ∼50 µm at p T > 7 GeV/c. This allows the fraction of non-prompt J/ψ to be determined for events with J/ψ p T greater than 1 GeV/c as well as in five transverse momentum intervals (1-3, 3-5, 5-7, 7-10 and 10-14 GeV/c). The projections of the maximised likelihood fit function superimposed over the m e + e − and x distributions of J/ψ candidates with p T > 1 GeV/c are shown as an example in figure 2.
The f b fraction is obtained after correcting the f b values (eq. (2.5) ) to account for slightly different A × factors of prompt and non-prompt J/ψ: (2.5) This small correction is computed relying on MC simulations assuming prompt J/ψ to be unpolarised. A small residual polarisation, resulting from the admixture of the different b-hadron decay channels, is assumed for non-prompt J/ψ as predicted by EvtGen [82]. Under these conditions, the correction mainly originates from the difference in the p T distribution between the two components and is found to be significant only for the p Tintegrated case. Small relative variations of the corrected f b values, in the order of ∼1-4%, are expected in the case of a null-polarisation assumption for non-prompt J/ψ [51]. The variations estimated in extreme polarisation scenarios for prompt J/ψ are discussed in ref.
[13]. Considering the null or very small degree of polarisation measured in pp collisions at the LHC [17, [92][93][94], these variations are not further propagated to the final results, and only the choice of the p T shapes used as input for the MC simulations is taken into account for the systematic uncertainty evaluation, as discussed below.
The evaluated systematic uncertainties affecting the measurements of f b in the five p T intervals as well as in the p T -integrated range (p T > 1 GeV/c) are listed in table 2.

JHEP06(2022)011
Most of the listed contributions are due to incomplete knowledge of the different PrDFs used as input for the likelihood fits. An additional contribution originates from the assumptions on the p T distributions employed for the computation of the correction factor of eq. (2.5). The uncertainties affecting the evaluation of the resolution function and of the background PrDF constitute the largest contributions to the total systematic uncertainty of the f b measurements. The former is estimated by propagating to the R(x) PrDF the residual discrepancy of the single-track impact parameter distributions between data and MC simulations after the application of the previously discussed tuning procedure. The latter is evaluated by repeating the likelihood fits after varying the procedure used for the determination of the F Bkg (x) PrDFs, following the same approach described in ref. [87]. Both uncertainties increase towards low transverse momenta and are largest for the lowest p T interval, where they amount to 10% and 8.5%, respectively. The uncertainty related to the invariant mass PrDF of the J/ψ signal is estimated by changing the width of the CB function used to parameterise the M Sig PrDF so as to vary the fraction of signal enclosed within the 2.92 < m e + e − < 3.16 GeV/c 2 interval by ±2.5%. The likelihood fits are repeated, and the variation of the resulting f b values is taken as systematic uncertainty. The uncertainty related to the invariant mass background PrDF is estimated as the RMS of the f b value distributions obtained after employing different parameterisations and alternative fitting approaches for the evaluation of the M Bkg PrDF. The estimate of the systematic uncertainty affecting the non-prompt J/ψ x PrDF is performed by repeating the likelihood fits after employing PYTHIA 6.4 for the description of the p T -distribution, decay kinematics, and relative abundance of the beauty hadrons in the MC simulations used to model the F b (x) PrDF. The relative variation of the resulting f b values, assumed as systematic uncertainty, increases up to ∼3% towards low transverse momenta. This uncertainty also contemplates any conservative assumption for a rapidity dependence of the relative abundances of beauty hadrons at the LHC. The uncertainty related to the acceptance times efficiency correction procedure is assessed by testing different hypothesis for the kinematic p T -spectra used to compute the A × values that enter into eq. (2.5). Among the tested variations, a tune-on-data parameterisation based on the Run 1 measurement [51] for prompt J/ψ, a p T -distribution based on FONLL calculations [90] for the non-prompt component, and the inclusion or exclusion of nuclear shadowing modifications according to the EPPS16 parameterisation [95] are considered. The resulting variations of the correction factors are largest for the p T -integrated measurement, where they amount to ∼3%, while they are smaller than ∼1% within the analysed p T -intervals. The overall systematic uncertainty of the f b measurements is found to increase up to 13.7% towards low transverse momenta, mostly as a consequence of both the increasing combinatorial background and the worsening of the x resolution.
The prompt and non-prompt J/ψ nuclear modification factors are computed by combining the measurements of f b with the previously discussed nuclear modification factors R incl. J/ψ pPb of inclusive J/ψ: The value of f b in pp collision at √ s = 5.02 TeV, indicated as f pp b in eq. (2.6), is determined by means of the same interpolation procedure adopted in previous analyses for the p−Pb [51] and Pb-Pb [87] systems. The procedure consists of fitting to existing midrapidity f b measurements at √ s = 1.96 TeV (from CDF [3]) and √ s = 7 TeV (from ALICE [13], ATLAS [96], and CMS [97]) the semi-phenomenological function discussed in ref. [87], which includes FONLL predictions [90] for the non-prompt J/ψ production cross section. An energy interpolation is then performed to derive the f pp b (p T ) at √ s = 5.02 TeV as a function of p T . The average value of f pp b in a given p T interval is obtained by weighting f pp b (p T ) over the inclusive J/ψ spectrum in pp collisions in that p T interval. Compared to our previous estimates [51], a tune-on-data spectrum based on the inclusive J/ψ yield measured by ALICE in pp collisions at √ s = 5.02 TeV [7] is now employed for this purpose. The values of f pp b at √ s = 5.02 TeV, computed in the considered momentum intervals, are reported in table 3. The quoted uncertainties take into account the uncertainties of both data and FONLL predictions, as well as an additional systematic uncertainty due to the choice of the functional form (either a linear, or an exponential, or a power law function) employed for the energy interpolation procedure.

Results
The inclusive J/ψ cross section is measured in −1.37 < y < 0.43 both for p T > 0 and differentially in p T considering seven p T intervals, with the first and last bins being [0−1] GeV/c and [10−14] GeV/c, respectively. The value of the p T -integrated inclusive J/ψ cross section per unit of rapidity is dσ/dy = 999 ± 33 (stat.) ± 56 (syst.) µb. The p T -differential cross section of inclusive J/ψ per unit of rapidity, d 2 σ incl. J/ψ /dydp T , is shown in figure 3 in comparison with the cross section measured in pp collisions at √ s = 5.02 TeV [7] multiplied by A = 208. The latter extends up to p T = 10 GeV/c. The highest p T point for pp collisions, which is shown in the figure with the empty symbol, was obtained using the interpolation procedure as described before.  The fraction of J/ψ from b-hadron decays in the kinematic range p T > 1 GeV/c and −1.37 < y < 0.43, which is referred to as "visible region" in the following, is found to be f b = 0.125 ± 0.017 (stat.) ± 0.011 (syst.), where the first quoted uncertainty is statistical and the second one is systematic. The f b measurements in the five analysed p T intervals are shown in figure 4 in comparison with our previous results [51] and with the results from the ATLAS collaboration [60], measured for p T > 8 GeV/c within a similar rapidity interval (−1.94 < y < 0). The measurements from the CDF [3], ATLAS [98], and CMS [97] experiments in pp and pp collisions at midrapidity are also shown for comparison. With respect to our previous results [51], the present measurements are performed over a wider   and ATLAS [60] collaborations in the same collision system (blue and green closed circles, and the blue arrow that shows an upper limit at 95% CL in the range 1.3 < p T < 3 GeV/c p T range, with a more granular binning, and show a significantly improved precision, with about half of the statistical uncertainty within similar p T intervals.
The prompt and non-prompt J/ψ production cross sections are obtained from the combination of the f b fractions with the measurements of the inclusive J/ψ cross section σ J/ψ : The non-prompt J/ψ cross section in the visible region, σ vis J/ψ from h b = 201±28 (stat.)± 21 (syst.) µb, is computed using the inclusive J/ψ cross section for p T > 1 GeV/c, which amounts to 1603 ± 55 (stat.) ± 89 (syst.) µb.
In order to derive the p T -integrated values of the prompt and non-prompt J/ψ cross section at midrapidity, σ vis J/ψ from h b is extrapolated down to p T = 0 following the approach described in our previous work [51]. The extrapolation is performed assuming the shape of the p T distribution of b-quarks obtained from FONLL [90] with the CTEQ6.6 PDFs [99] modified according to EPPS16 nPDF parameterisation [95]. The fragmentation of b-quarks into hadrons is then modelled using PYTHIA 6.4 [88] with the Perugia-0 tune [89]. The ratio of the extrapolated cross section for p T > 0 and −1.37 < y < 0.43 to that in the visible region (p T > 1 GeV/c and −1.37 < y < 0.43) equals 1.127 +0.014 −0.025 , where the quoted uncertainty takes into account the FONLL, CTEQ6.6 and EPPS16 uncertainties, as described in ref. [51], as well as an additional uncertainty, which is related to that on the -14 -JHEP06(2022)011 In the left plot, the systematic uncertainty of the ALICE data point includes also the contribution from the extrapolation procedure to go from the visible region (p T > 1 GeV/c) to p T > 0, as described in the text. For the measurements as a function of p T , the data symbols are placed within each bin at the mean of the p T distribution determined from MC simulations. The results of a model [100][101][102][103] including nuclear shadowing based on the EPPS16 [95] and nCTEQ15 [104] nPDFs are shown superimposed on both panels (see text for details). In the right panel, the computations refer only to the ALICE rapidity range.
relative abundance of beauty hadron species in the extrapolated range. The latter is estimated to be about 0.4% after changing the assumed fractions of beauty hadrons according to the recent LHCb measurements [91]. Also, in this case, the considered variation largely includes a possible dependence of these fractions on rapidity. Thus, the measured cross section corresponds to more than 85% of the p T -integrated cross section at midrapidity. Dividing by the rapidity range ∆y = 1.8, the following value is derived for the non-prompt J/ψ cross section per unit of rapidity (p T > 0 and −1.37 < y < 0.43): dσ J/ψ from h b dy = 125.6 ± 17.6 (stat.) ± 13.3 (syst.) +1.6 −2.8 (extr.) µb.
The corresponding value for the prompt component is obtained as the difference between the inclusive J/ψ cross section, which is measured for p T > 0, and that of J/ψ from bhadron decays, as determined with the extrapolation procedure described above. It is (p T > 0 and −1.37 < y < 0.43) dσ prompt J/ψ dy = 873.1 ± 33.6 (stat.) ± 50.4 (syst.) +1.6 −2.8 (extr.) µb.
In figure 5 (left panel) this result is shown as a function of rapidity together with the results from the LHCb experiment at positive ("forward") and negative ("backward") rapidity [64], corresponding respectively to the p-going and Pb-going direction. The p Tdifferential cross section of prompt J/ψ is shown, in comparison with ATLAS measurements [14] at high p T and for −2 < y < 1.5, in the right panel of figure 5. The ALICE The results are compared to FONLL computations [90] with EPPS16 [95] nPDFs, highlighting the total theoretical uncertainty (empty band) and the contribution from EPPS16 (coloured band). In the right panel, model computations are obtained in the same rapidity range of the ALICE results, namely −1.37 < y < 0. 43. results, covering the low p T region at midrapidity, are complementary to the measurements from both the LHCb and ATLAS collaborations. The data are reported in comparison with model calculations for prompt J/ψ (Lansberg et al. [100][101][102][103]) based on the EPPS16 [95] and the nCTEQ15 [104] sets of nuclear parton distribution functions (nPDFs). In both cases, the shaded bands represent the envelope of the computations for different assumptions of the values of the pQCD factorisation (µ F ) and renormalisation (µ R ) scales (varied within 0.5 < µ F /µ R < 2) computed at the 90% confidence level. The predictions show good agreement with data within the large model uncertainties, which are dominated by those on the pQCD scales. The results of a Bayesian reweighting approach from the same authors [101], employing LHCb measurements of J/ψ [66,105] as a constraint for the computations, are also shown in the left panel of figure 5. Both the size of uncertainties and the difference between the nPDF sets are largely reduced after the reweighting. In figure 6, the cross sections of non-prompt J/ψ, computed either for p T > 0 (left panel) or differentially in p T (right panel), are reported together with the corresponding results from the LHCb [64] and ATLAS [14] collaborations. The results are compared with theoretical predictions based on FONLL pQCD calculations [90] with the inclusion of nuclear shadowing effects according to the EPPS16 nPDFs [95]. In each panel, the coloured curves delimit the total theoretical uncertainty on the production cross section, which is dominated by that of the b-quark mass and the pQCD scales, while the shaded bands refer to the theoretical uncertainty of the EPPS16 nPDFs. The nuclear modification factor of inclusive J/ψ, measured for p T > 0 and −1.37 < y < 0.43, amounts to 0.851 ± 0.028 (stat.) ± 0.079 (syst.). This quantity is obtained using the measured inclusive J/ψ cross section in pp collisions at √ s = 5.02 TeV [7]. Similarly, the p T -differential R pPb of inclusive J/ψ is obtained on the basis of the measured pp reference, except for the highest p T interval (10-14 GeV/c), where the statistical sample is limited and the interpolation procedure [57] is still used. In figure 7, the R pPb of prompt J/ψ is reported either for p T > 0 in comparison with LHCb measurements [64] at backward and forward rapidity (left panel) or as a function of p T , computed according to eq. (2.6), together with that of inclusive J/ψ (right panel) in comparison with ATLAS results [14]. The p T -integrated R pPb of prompt J/ψ at midrapidity (p T > 0 and −1.37 < y < 0.43) is measured to be smaller than unity and amounts to 0.860 ± 0.033 (stat.) ± 0.081 (syst.). Given also the relatively small fraction of J/ψ from b-hadron decays for p T < 14 GeV/c, the R pPb of inclusive J/ψ is comparable with that of the prompt component. As shown in the right panel of figure 7, both trends indicate that the suppression observed at midrapidity is a low-p T effect, concentrated for p T 3 GeV/c. The measurements are compared with results from various model predictions which embed different CNM effects into prompt J/ψ production. In addition to the previously described computations by Lansberg et al., which include a reweighting of the EPPS16 and nCTE15 nPDFs [101], the central values of a computation based on EPS09 nPDF with or without interaction with a nuclear medium (Ferreiro et al. [106]) are shown. A calculation including the effects of coherent energy loss (Arleo et al. [44]), with or without the introduction of nuclear shadowing effects according to EPS09 nPDF, provides a fairly good description of the measurements either The nuclear modification factor for non-prompt J/ψ, determined according to eq. (2.6), is shown in figure 8. The measured value of the p T -integrated R pPb for p T > 0 and −1.37 < y < 0.43 is found to be 0.79±0.11 (stat.)±0.13 (syst.) +0.01 −0.02 (extr.), suggesting the presence of nuclear effects also for the non-prompt J/ψ component. Within uncertainties, the measurements are found to be compatible with those of the LHCb collaboration [64] at both forward and backward rapidity as well as with those of the ATLAS collaboration [14] for p T 9 GeV/c. All data are in fair agreement with the mild degree of suppression predicted by FONLL computations employing the EPPS16 nPDFs. The nuclear modification factor, as predicted from the Bayesian reweighting approach [101] of the EPPS16 nPDFs previously introduced for the prompt component, also provides a good description of the measurements. The central value of an alternative parameterisation of the nuclear PDF, nDSgLO [109], is reported in the left-hand panel. Despite the larger relative uncertainties, the comparison with the results shown for the prompt component (right panel of figure 7) suggests that a reduced suppression as well as a less pronounced p T -dependence affect the component of J/ψ from b-hadron decays.
Similarly as for our previous work [51], the low p T coverage of the measured nonprompt J/ψ cross section at midrapidity, now extending down to p T = 1 GeV/c, allows the p T -integrated bb cross section per unit of rapidity, dσ bb /dy, and the total bb production cross section, σ(pPb → bb + X), to be derived with small extrapolation uncertainties. By -18 -JHEP06(2022)011 using FONLL with CTEQ6.6 and EPPS16 nPDFs as input model for the computation of the extrapolation factor, the bb production cross section at midrapidity is derived as Assuming the average branching ratio of J/ψ from b-hadron decays measured at LEP [110][111][112], BR(h b → J/ψ + X) = (1.16 ± 0.10)%, for the computation of σ vis, model J/ψ from h b , the resulting midrapidity cross section per unit of rapidity is dσ bb dy = 5.52 ± 0.77 (stat.) ± 0.75 (syst.) +0.07 −0.12 (extr.) mb.
With a similar approach, the total bb production cross section is obtained by extrapolating the visible cross section to the full phase space: where the factor α 4π is defined as the ratio of the yield of non-prompt J/ψ produced in the full phase space to that in the visible region, and the factor 2 takes into account that hadrons with b or b valence quark can decay into J/ψ. The value of the α 4π factor obtained from FONLL pQCD calculations with EPPS16 nPDFs, with the b-quark fragmentation performed using PYTHIA 6.4 with the Perugia-0 tune, is α 4π = 4.10 +0.14 −0.12 . Using PYTHIA 8 instead of FONLL for the generation of bb quark pairs provides the value 4.02, which is 2% smaller than the used value for the extrapolation based on FONLL. The total cross section is 2 σ(p + Pb → bb + X) = 35.5 ± 5.0 (stat.) ± 4.8 (syst.) +1.2 −1.0 (extr.) mb.
The reported results for the extrapolated bb production cross sections are consistent with our previous derivations [51]. The total uncertainty is reduced by a factor of 2 thanks to the larger data sample, the smaller systematic uncertainty, and the slightly extended coverage of the visible region where the non-prompt J/ψ cross section is measured. In p-Pb collisions at √ s NN = 5.02 TeV the LHCb Collaboration measured the nonprompt J/ψ production cross section at forward and backward rapidities for p T < 14 GeV/c, reporting [64] σ J/ψ from h b (1.5 < y < 4.0) = 166.0 ± 4.1 ± 8.2 µb and σ J/ψ from h b (−5.0 < y < −2.5) = 118.2±6.8±11.7 µb, respectively. A more precise estimate of the total bb cross section can be obtained by repeating the same procedure with including also these results from the LHCb collaboration [64], to obtain a wider visible region: (−5 < y < −2.5, p T < 14 GeV/c)∪(−1.37 < y < 0.43, p T > 1.0 GeV/c)∪(1.5 < y < 4, p T < 14 GeV/c). The cross section in this wider visible region is obtained as the sum of the cross sections measured in 2 The extrapolation uncertainties for the dσ bb /dy and for the total bb cross section include the contributions related to the FONLL, CTEQ6.6 and EPPS16 uncertainties as discussed in ref. [51], and also a minor contribution of about 0.4% related to the uncertainty on the p T -dependence of the relative abundance of b-hadron species.

JHEP06(2022)011
this work at central rapidity and those from LHCb. All the uncertainties are uncorrelated except that of the branching ratio. In this case, the α 4π factor, which is calculated as the ratio of the yield from the model in full phase space to that in the wider region covered by the ALICE and LHCb experiments, is reduced to 1.60 ± 0.02, and the corresponding total cross section is σ(p + Pb → bb + X) = 33.8 ± 2.0 (stat.) ± 3.4 (syst.) +0.4 −0.5 (extr.) mb (ALICE and LHCb).

Summary
The production of J/ψ mesons in p−Pb collisions at √ s NN = 5.02 TeV is studied based on a data sample about six times larger than that of previously published results yielding smaller uncertainties and extending the p T coverage. The inclusive J/ψ production cross section at midrapidity is measured down to p T = 0 after reconstructing the J/ψ mesons in the dielectron decay channel. The fraction of the inclusive J/ψ yield originated from b-hadron decays is then determined on a statistical basis, allowing the prompt and non-prompt J/ψ production cross sections at midrapidity to be derived for p T > 1 GeV/c and as a function of p T in five momentum intervals. The results are scaled to reference measurements from pp collisions at the same centre-of-mass energy in order to investigate the presence of nuclear effects on J/ψ production. The nuclear modification factor of prompt J/ψ shows a significant suppression for p T 3 GeV/c, whereas there is a hint of a less pronounced suppression of the non-prompt component over the inspected p T range. The results can be described by theoretical calculations including various combinations of cold nuclear matter effects, although a precise discrimination among the different models is impaired by the uncertainties affecting the currently available predictions. Finally, the measurement of the non-prompt J/ψ production cross section is used to derive the extrapolated midrapidity dσ bb /dy and total cross section, σ bb , of beauty quark production. Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.