Abstract
The total charm-quark production cross section per unit of rapidity \(\textrm{d}\sigma ({{\textrm{c}}\overline{\textrm{c}}})/\textrm{d}y\), and the fragmentation fractions of charm quarks to different charm-hadron species \(f(\textrm{c}\rightarrow {\textrm{h}}_{\textrm{c}})\), are measured for the first time in p–Pb collisions at \(\sqrt{s_\textrm{NN}} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \) at midrapidity (\(-0.96<y<0.04\) in the centre-of-mass frame) using data collected by ALICE at the CERN LHC. The results are obtained based on all the available measurements of prompt production of ground-state charm-hadron species: \(\textrm{D}^{0}\), \(\textrm{D}^{+}\), \(\textrm{D}_\textrm{s}^{+}\), and \(\mathrm {J/\psi }\) mesons, and \(\Lambda _\textrm{c}^{+}\) and \(\Xi _\textrm{c}^{0}\) baryons. The resulting cross section is \( \textrm{d}\sigma ({{\textrm{c}}\overline{\textrm{c}}})/\textrm{d}y =219.6 \pm 6.3\;(\mathrm {stat.}) {\;}_{-11.8}^{+10.5}\;(\mathrm {syst.}) {\;}_{-2.9}^{+8.3}\;(\mathrm {extr.})\pm 5.4\;(\textrm{BR})\pm 4.6\;(\mathrm {lumi.}) \pm 19.5\;(\text {rapidity shape})+15.0\;(\Omega _\textrm{c}^{0})\;\textrm{mb} \), which is consistent with a binary scaling of pQCD calculations from pp collisions. The measured fragmentation fractions are compatible with those measured in pp collisions at \(\sqrt{s} = 5.02\) and 13 TeV, showing an increase in the relative production rates of charm baryons with respect to charm mesons in pp and p–Pb collisions compared with \(\mathrm {e^{+}e^{-}}\) and \(\mathrm {e^{-}p}\) collisions. The \(p_\textrm{T}\)-integrated nuclear modification factor of charm quarks, \(R_\textrm{pPb}({\textrm{c}}\overline{\textrm{c}})= 0.91 \pm 0.04\;\mathrm{(stat.)} ^{+0.08}_{-0.09}\;\mathrm{(syst.)} ^{+0.05}_{-0.03}\;\mathrm{(extr.)} \pm 0.03\;\mathrm{(lumi.)}\), is found to be consistent with unity and with theoretical predictions including nuclear modifications of the parton distribution functions.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
The processes governing the production of heavy-flavour hadrons (those containing at least one charm or beauty quark) in hadronic collisions have recently become a focal point of study at the CERN LHC. Due to their large masses (\(m_\textrm{c,b} \gg \Lambda _\textrm{QCD}\)), heavy-flavour hadron production is calculable using perturbative quantum chromodynamics (pQCD) under the assumption of factorisation [1, 2]. In this approach, the production cross section of a given hadron species is considered as a convolution of: (i) the parton distribution functions (PDFs) of the colliding hadrons; (ii) the partonic hard-scattering cross sections for a heavy quark to be produced; and (iii) the fragmentation functions, which describe the fraction of the heavy-quark momentum carried by the heavy-flavour hadron. The fragmentation term in the calculation also governs the relative abundances of each hadron species, which are often referred to as ‘fragmentation fractions’. Typically, in the factorisation approach, the hadronisation process is taken to be independent of the collision system. Therefore, the fragmentation functions measured in \(\textrm{e}^{+}\textrm{e}^{-}\) collisions are applied equally to calculations in proton–proton (pp) collisions. The assumption of fragmentation as a universal hadronisation mechanism would also result in identical fragmentation fractions between \(\textrm{e}^{+}\textrm{e}^{-}\) and pp collisions. This has been tested extensively by measuring the production cross section ratios of different charm-hadron species.
Measurements of the \(\textrm{D}^+/\textrm{D}^\textrm{0} \) and \(\textrm{D}_{\textrm{s}}^+/\textrm{D}^\textrm{0} \) meson-to-meson production ratios in pp collisions at different energies by the ALICE Collaboration [3, 4] have shown consistent production rates of charm mesons between hadronic collisions and earlier measurements from \(\textrm{e}^{+}\textrm{e}^{-}\) and \(\textrm{e}^-\textrm{p}\) collisions [5]. However, more recent measurements of the \(\Lambda _{\textrm{c}}^{+}/\textrm{D}^\textrm{0} \) and \({\Xi _{\textrm{c}}^{0,+}}/\textrm{D}^\textrm{0} \) production ratios in pp and heavy-ion collisions [4, 6,7,8,9,10,11] have shown a marked increase in charm-baryon production with respect to charm-meson production over the \(\textrm{e}^{+}\textrm{e}^{-}\) baseline in hadronic collisions at LHC energies. This was further quantified by determining the charm-hadron fragmentation fractions in pp collisions at \(\sqrt{s} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \) [3, 4] and \(\sqrt{s} = 13~\text {Te}\hspace{-1.00006pt}\textrm{V} \) [4]. In both cases, a significant relative enhancement of charm-baryon production was observed, along with a corresponding depletion in the charm-meson fragmentation fractions. The results are also consistent within uncertainties between \(\sqrt{s} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \) and \(\sqrt{s} = 13~\text {Te}\hspace{-1.00006pt}\textrm{V} \), implying that the collision energy does not play a significant role in the relative abundances of charm-hadron species in hadronic collisions at the LHC. The increase in charm-baryon production is also reflected in an increase of the total charm production cross section in pp collisions over earlier determinations using only the measured D-meson cross sections with the fragmentation fractions from leptonic colliders [12,13,14]. More recently, efforts have also been made to extrapolate and combine the available measurements of charm hadrons in pp collisions from different LHC experiments in complementary regions of phase space in order to determine the total charm production cross section [15].
Measurements of the charm fragmentation fractions in p–Pb collisions are also of particular interest, as they reveal whether the observed changes in relative hadronisation rates for \(\textrm{e}^{+}\textrm{e}^{-}\) and pp collisions are mainly due to the increase of the particle multiplicity between these collision systems. The mean charged-particle multiplicity per unit of pseudorapidity \(\left. \left\langle \textrm{d}N_\textrm{ch}/\textrm{d}\eta \right\rangle \right| _{|\eta | < 1.0}\) in collisions at the LHC increases both as a function of energy and system size – for minimum-bias pp collisions at \(\sqrt{s} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \), \(\left. \left\langle \textrm{d}N_\textrm{ch}/\textrm{d}\eta \right\rangle \right| _{|\eta | < 1.0} = 5.48\pm 0.05\;(\text {uncorr.\,syst.})\pm 0.05\;(\text {corr.\,syst.})\) [16], and for minimum-bias p–Pb collisions at \(\sqrt{s_\textrm{NN}} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \), \(\left. \left\langle \textrm{d}N_\textrm{ch}/\textrm{d}\eta \right\rangle \right| _{|\eta | < 1.0}=16.81\pm 0.71\;(\mathrm{syst.})\) [17]. The ALICE Collaboration measured the \(\Lambda _{\textrm{c}}^{+}/\textrm{D}^\textrm{0} \) yield ratio in intervals of charged-particle multiplicity in pp collisions at \(\sqrt{s} = 13~\text {Te}\hspace{-1.00006pt}\textrm{V} \) [18], showing that in a hadronic collision, even at relatively small multiplicities, charm-quark hadronisation proceeds differently from \(\textrm{e}^{+}\textrm{e}^{-}\) collisions. In the beauty sector, the measurements of \({\Lambda _{\text {b}}^{0}}\)-baryon production relative to that of B mesons at forward rapidity in pp and p–Pb collisions by the LHCb Collaboration show a modification of the fragmentation fractions among collision systems similar to that observed for charm quarks at midrapidity [19,20,21,22]. The ALICE Collaboration measured production cross sections of \(\textrm{D}^\textrm{0}\) and \(\Lambda _{\textrm{c}}^{+}\) hadrons originating from beauty-hadron decays at midrapidity in proton–proton collisions at a centre-of-mass energy \(\sqrt{s}\) = 13 TeV [23], and a baryon-to-meson enhancement similar to the one measured by LHCb at forward rapidity was observed.
The study of charm production in p–Pb collisions serves as an additional baseline to examine so-called cold nuclear matter (CNM) effects at the LHC. Due to different effects related to the presence of a nucleus in the colliding system, there may also be corresponding changes to the yield and transverse momentum (\(p_\textrm{T}\)) distribution of charm hadrons. These effects include modifications of the PDFs in nuclei (nPDFs) with respect to the proton PDFs [24, 25] or \(k_\textrm{T}\)-broadening due to multiple soft collisions before the heavy-quark pair is produced [26, 27]. The resulting modifications of the yields are typically examined using the nuclear modification factor \(R_\textrm{pPb}\) [28]. The ALICE measurements of \(\Lambda _{\textrm{c}}^{+}\) production in pp and p–Pb collisions at \(\sqrt{s_\textrm{NN}} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \) [6] revealed a significant modification of the \(p_\textrm{T}\) distribution of the charm baryon-to-meson ratio between p–Pb and pp collisions, which was not observed for the charm meson-to-meson ratios within the current measurement uncertainties. In addition, an increase of the mean transverse momentum of \(\Lambda _{\textrm{c}}^{+}\) baryons was observed in p–Pb collisions with respect to pp collisions, with no significant difference seen for \(\textrm{D}^\textrm{0}\) mesons. This hardening of the \(p_\textrm{T}\) spectrum of \(\Lambda _{\textrm{c}}^{+}\) baryons in p–Pb collisions indicates a possible influence of radial flow [29, 30] and/or hadronisation via coalescence [31,32,33]. These effects are consistent with observations from the light-flavour sector, such as measurements of the \(\Lambda /\textrm{K}^0_\textrm{S} \) baryon-to-meson ratio in p–Pb collisions [34]. However, the \(p_\textrm{T}\)-integrated \(R_\textrm{pPb}\) of \(\Lambda _{\textrm{c}}^{+}\) baryons was found to be consistent with that of \(\textrm{D}^\textrm{0}\) mesons and with unity. With the recent addition of \({\Xi _{\textrm{c}}^{0}}\) baryons to the available measurements [7], it has become possible to study potential nuclear modifications of the total \({\textrm{c}}\overline{\textrm{c}}\) production cross section to conclude whether nPDFs affect the overall creation of charm hadrons. A recent review with a more comprehensive overview of heavy-quark hadronisation can be found in Ref. [35].
In this article, the first measurements of the total charm production cross section per unit of rapidity and the charm-hadron fragmentation fractions at midrapidity in p–Pb collisions at the LHC are presented. Section 2 outlines the previous measurements by the ALICE Collaboration that were used to compute the total charm-quark production cross section and the fragmentation fractions of charm quarks to different charm-hadron species. It also details the theoretical models that were used to extrapolate the measured charm-hadron production cross sections down to \(p_\textrm{T} =0\), where necessary. Section 3 shows the results for the charm-hadron fragmentation fractions and \({\textrm{c}}\overline{\textrm{c}}\) production cross section in p–Pb collisions at \(\sqrt{s_\textrm{NN}} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \), along with the \(p_\textrm{T}\)-integrated nuclear modification factor \(R_\textrm{pPb}\). A summary is given in Sect. 4.
2 Dataset and analysis
The ALICE detector system is described in detail in Ref. [36]. The measured prompt charm-hadron production cross sections used in this article were obtained from a sample of minimum-bias p–Pb collisions collected in 2016 during Run 2 of the LHC with a total integrated luminosity of \(\mathcal {L}_\textrm{int} =287\pm 11\;\mu \textrm{b} ^{-1}\) [37]. The measurements were performed at midrapidity (defined as \(-0.96<y<0.04\) in the collision centre-of-mass frame). Prompt charm hadrons are those produced directly in the hadronisation of a charm quark or via the decay of a directly produced excited charm-hadron or charmonium state, while charm hadrons produced by the decays of beauty hadrons are denoted as ‘non-prompt’. D mesons were reconstructed in the decay channels \(\textrm{D}^\textrm{0} \rightarrow \textrm{K}^-\pi ^+\), \(\textrm{D}^+ \rightarrow \textrm{K}^-\pi ^+\pi ^+\), \(\textrm{D}_{\textrm{s}}^+ \rightarrow \phi \pi ^+\), and \(\textrm{D}^{*+} \rightarrow \textrm{D}^\textrm{0} \pi ^+\) (and respective charge conjugates) [28]. Charm baryons were reconstructed in the channels \(\Lambda _{\textrm{c}}^{+}\rightarrow \textrm{p} K^{0}_{\textrm{S}}\) [6, 38], \(\Lambda _\textrm{c}^{+}\rightarrow \textrm{p} K^-\pi ^+\) [38], and \({\Xi _{\textrm{c}}^{0}} \rightarrow \Xi ^-\pi ^+\) [7] (and respective charge conjugates). The \({\textrm{J}}/\psi \) mesons were reconstructed in the channel \({\textrm{J}}/\psi \rightarrow \textrm{e}^{+}\textrm{e}^{-} \) [39]. The details of the hadron reconstruction methods can be found in the corresponding citations.
In order to derive the total charm production cross section and the individual charm-hadron fragmentation fractions at midrapidity, the \(p_\textrm{T}\)-integrated rapidity-differential cross sections must be considered. For \(\textrm{D}^\textrm{0}\) mesons [28] and \(\Lambda _{\textrm{c}}^{+}\) baryons [6], the \(p_\textrm{T}\)-differential yields were measured down to \(p_\textrm{T} =0\) and the calculation of the \(p_\textrm{T}\)-integrated cross section did not require further extrapolation. For the remaining hadron species, where no measurement down to \(p_\textrm{T} =0\) is available, an extrapolation was performed based on model calculations.
For \(\textrm{D}^+\), \(\textrm{D}_{\textrm{s}}^+\), and \(\textrm{D}^{*+}\) mesons, POWHEG+PYTHIA 6 calculations were used to extrapolate the \(p_\textrm{T}\) spectra, which were measured down to \(p_\textrm{T} =1\) \(\text {Ge}\hspace{-1.00006pt}\textrm{V}/c\) for \(\textrm{D}^+\) and \(\textrm{D}^{*+}\), and \(p_\textrm{T} =2\) \(\text {Ge}\hspace{-1.00006pt}\textrm{V}/c\) for \(\textrm{D}_{\textrm{s}}^+\). POWHEG is an event generator with next-to-leading order (NLO) accuracy [40]. It is matched with PYTHIA 6 [41] to generate the parton shower and fragmentation. The CT14NLO parton distribution functions [42] with the EPPS16 parameterisation of the nPDFs [43] were used in the calculations. The POWHEG+PYTHIA 6 simulations describe the charm-meson production cross sections within uncertainties. The extrapolation factors were determined from the fraction of the theoretically calculated cross section in the visible region with respect to the \(p_\textrm{T}\)-integrated one, relying on the \(p_\textrm{T}\) shape of the production spectra without directly using the magnitude from theory. This was performed by taking the integral of the calculated production cross section for \(0<p_\textrm{T} <50\;\text {Ge}\hspace{-1.00006pt}\textrm{V}/c \) and dividing it by the corresponding integral in the measured \(p_\textrm{T}\) region for that species. The contribution to the total cross section for \(p_\textrm{T} >50\;\text {Ge}\hspace{-1.00006pt}\textrm{V}/c \) is negligible with respect to that for \(p_\textrm{T} <50\;\text {Ge}\hspace{-1.00006pt}\textrm{V}/c \), so any extrapolation to higher \(p_\textrm{T}\) values than this would have no impact on the final result. An uncertainty on this extrapolation factor was then determined by separately performing the same procedure using the following standard variations on the POWHEG+PYTHIA 6 model: (i) varying the factorisation and renormalisation scales in the POWHEG simulation (\(\mu _\textrm{r,f}\)) within the range \(0.5<\mu _{\textrm{r,f}}/\mu _0<2.0\) (where \(\mu _0 = \sqrt{m^2+p_\textrm{T} ^2}\)), with the additional constraint \(0.5<\mu _\textrm{r}/\mu _\textrm{f}<2\); and (ii) considering differences due to the uncertainty sets on the EPPS16 nPDF parametrisations. The maximum relative upward and downward differences from the central extrapolation factor provided by (i) and (ii) were added in quadrature to give the overall extrapolation uncertainty on the cross sections, which amounted to \( ^{+7.7}_{-2.2}\%\) for \(\textrm{D}^+\), \( ^{+18.8}_{-\phantom {0}7.3}\%\) for \(\textrm{D}_{\textrm{s}}^+\), and \( ^{+27.0}_{-\phantom {0}2.3}\%\) for \(\textrm{D}^{*+}\) mesons. In all cases, the contribution due to the values of \(\mu _\textrm{r,f}\) was much larger than that from the nPDF uncertainty sets. The extrapolation factors were \(1.27^{+0.10}_{-0.04}\) for \(\textrm{D}^+\) mesons, \(1.25^{+0.08}_{-0.04}\) for \(\textrm{D}^{*+}\) mesons, and \(2.08^{+0.40}_{-0.20}\) for \(\textrm{D}_{\textrm{s}}^+\) mesons.
For \({\Xi _{\textrm{c}}^{0}}\) baryons, the extrapolation was performed and reported in Ref. [7]. The POWHEG+PYTHIA 6 calculations were found to poorly describe the \(p_\textrm{T}\) shape of the measured differential cross sections of charm baryons in p–Pb collisions. Instead, for the case of \({\Xi _{\textrm{c}}^{0}}\), the quark (re)combination model (QCM) [44, 45] was chosen to perform the extrapolation, as it is tuned to reproduce previous measurements of \(\Lambda _{\textrm{c}}^{+}\)-baryon production from ALICE and provides a reasonable description of the shape of the measured \(\textrm{d}\sigma /\textrm{d}p_\textrm{T} \) of \({\Xi _{\textrm{c}}^{0}}\) baryons [7]. As the QCM model does not provide any theoretical uncertainties, the extrapolation uncertainty was defined by adding in quadrature (i) the maximum difference between the central extrapolation factor determined using QCM and those calculated using PYTHIA 8 with the Angantyr generator [46] and POWHEG+PYTHIA 6 (\( ^{+20.6}_{-\phantom {0}0.0}\%\)), and (ii) the envelope of the POWHEG+PYTHIA 6 variations for \({\Xi _{\textrm{c}}^{0}}\) as described above for the D mesons (\( ^{+23.0}_{-\phantom {0}8.6}\%\)). The resulting extrapolation uncertainty is \( ^{+30.9}_{-\phantom {0}8.6}\%\) of the total \({\Xi _{\textrm{c}}^{0}}\) production cross section at midrapidity. The final extrapolation factor for \({\Xi _{\textrm{c}}^{0}}\) baryons is \(1.74^{+0.54}_{-0.15}\). As no measurement of the \({\Xi _{\textrm{c}}^{+}}\) production cross section is available in p–Pb collisions at midrapidity, the contribution of \({\Xi _{\textrm{c}}^{+}}\) to the total charm cross section was considered as equal to that of \({\Xi _{\textrm{c}}^{0}}\) baryons, as also done in Refs. [3, 4]. This assumption was made based on isospin symmetry, and is further supported by the measurements of \({\Xi _{\textrm{c}}^{0,+}}\) in pp collisions at \(\sqrt{s} = 13~\text {Te}\hspace{-1.00006pt}\textrm{V} \) [4], which showed that the production cross sections are consistent between the two charge states.
For \({\textrm{J}}/\psi \) mesons, the extrapolation was performed and reported in Ref. [39]. While the inclusive \({\textrm{J}}/\psi \)-meson cross section (defined as the sum of prompt and non-prompt \({\textrm{J}}/\psi \) mesons) was measured for \({p_\textrm{T} >0}\), an extrapolation was required for the non-prompt \({\textrm{J}}/\psi \) cross section to account for the \(p_\textrm{T}\) region \({0<p_\textrm{T} <1\;\text {Ge}\hspace{-1.00006pt}\textrm{V}/c}\) in order to extract the prompt cross section. The extrapolation was performed using Fixed Order plus Next-to-Leading Logarithm (FONLL) pQCD calculations [47, 48], with the standard CTEQ6.6 PDFs [49] modified according to the EPPS16 nPDF parametrisation [43] to describe p–Pb collisions. The corresponding value for the prompt component is obtained as the difference between the inclusive \({\textrm{J}}/\psi \)-meson cross section and the non-prompt one, as determined with the extrapolation procedure described above. The \(\textrm{d}\sigma /\textrm{d}y\) value for prompt \({\textrm{J}}/\psi \)-meson production from Ref. [39] was used directly in this work.
In the calculation of the total charm-hadron yield, the contribution from \(\Omega _\textrm{c}^{0}\) baryons must also be considered. Currently, no measurements of the \(\Omega _\textrm{c}^{0}\) are available in p–Pb collisions. In addition, as there is no measurement of the absolute branching ratios for any of the decay modes of the \(\Omega _\textrm{c}^{0}\) baryon, it is not yet possible to directly measure its production cross section and so it is not included in the sum of the charm-hadron cross sections. The potential contribution from \(\Omega _\textrm{c}^{0}\) baryons was instead taken into account by assuming an extreme upper limit of \(\sigma (\Omega _\textrm{c}^{0}) = \sigma ({\Xi _{\textrm{c}}^{0}})\), and assigning this as a fully asymmetric uncertainty bound on the total charm-hadron yield, as was previously done in Refs. [3, 4]. The uncertainty related to the estimated \(\Omega _\textrm{c}^{0}\)-baryon production cross section is added in quadrature to the upper extrapolation uncertainty on the total cross section. There are further contributions to the total charm cross section from doubly- and triply-charmed baryons such as \(\Xi _\textrm{cc}^{+}\) and \(\Omega _\textrm{ccc}^{++}\) as well as charmonium states that do not decay to \({\textrm{J}}/\psi \) mesons. However, these are considered to be negligible compared to the precision of the current measurement of the total cross section.
When calculating the total \({\textrm{c}}\overline{\textrm{c}}\) production cross section, additional correction factors were applied to account for the differing rapidity distributions between charm hadrons, individual charm quarks, and \({\textrm{c}}\overline{\textrm{c}}\) pairs. These were computed using the same procedure as in Refs. [3, 4]. The scaling factor between charm hadrons and single charm quarks was determined from FONLL calculations of pp collisions at \(\sqrt{s} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \) [50] and found to be unity, with a 2% uncertainty assigned based on the differences between FONLL and PYTHIA 8 calculations. The scaling between single charm quarks and charm-anticharm pairs was determined using POWHEG+PYTHIA 6 calculations. In the range \(-0.96<y<0.04\), this correction factor was found to have a value of 1.03. An uncertainty was assigned by adding in quadrature the envelopes defined by (i) variations of the \(\mu _\textrm{r,f}\) scale parameters in POWHEG and (ii) the uncertainty sets of the EPPS16 nPDFs, with a final value of \(8.7\%\). Combining these two factors, the uncertainty on the total charm cross section due to the rapidity shape is 8.9%.
In calculating the fragmentation fractions and total \({\textrm{c}}\overline{\textrm{c}}\) production cross section at midrapidity, the systematic uncertainties related to tracking, feed-down from beauty-hadron decays, \(p_\textrm{T}\) extrapolation, and luminosity were propagated as fully correlated among different charm-hadron species; all other uncertainties were propagated as uncorrelated.
3 Results
The \(p_\textrm{T}\)-integrated \(\left. \textrm{d}\sigma (\textrm{h}_\textrm{c})/\textrm{d}y\right| _{-0.96<y<0.04}\) values for all of the considered charm-hadron species in p–Pb collisions at \(\sqrt{s_\textrm{NN}} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \) are listed in Table 1. For all species, the systematic uncertainties related to the extrapolation (extr.), luminosity determination (lumi.), and branching ratio (BR) are reported separately from the systematic uncertainties related to the data analysis, except where these values were combined or not applicable in their original publications.
The fragmentation fractions \(f(\textrm{c}\rightarrow {\textrm{h}}_{\textrm{c}})\) were computed from the integrated production cross sections of the charm hadrons according to
In the denominator, only ground-state charm hadrons were considered. The production cross section of \({\Xi _{\textrm{c}}^{0}}\) was considered twice to represent both of the \({\Xi _{\textrm{c}}} \) charge states.The fragmentation fractions of all measured charm-hadron species in p–Pb collisions are reported in Fig. 1 (left panel) and summarised in Table 2. The results are compared with those measured previously in pp collisions at \(\sqrt{s} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \) by the ALICE Collaboration. The original measurements in pp collisions published in Ref. [3] were updated in Ref. [4] by taking into account more recent measurements of the \(\Lambda _{\textrm{c}}^{+}\)-baryon cross section down to \(p_\textrm{T} =0\) and of the prompt \({\textrm{J}}/\psi \) cross section. No significant modification of the hadronisation between the two colliding systems is observed despite the larger system size and higher charged-particle multiplicity density in p–Pb collisions. The values are consistent with those measured in pp collisions at \(\sqrt{s} = 13~\text {Te}\hspace{-1.00006pt}\textrm{V} \) [4]. The measurements are also compared with those from \(\textrm{e}^{+}\textrm{e}^{-}\) and \(\textrm{e}^-\textrm{p}\) collisions [5]. As observed for pp collisions, the fragmentation fractions in p–Pb collisions also exhibit an overall enhancement in the relative production of charm baryons and a corresponding deficit in charm-meson production with respect to \(\textrm{e}^{+}\textrm{e}^{-}\) and \(\textrm{e}^-\textrm{p}\) collisions.
The total \({\textrm{c}}\overline{\textrm{c}}\) production cross section at midrapidity was computed as the sum of the measured \(p_\textrm{T}\)-integrated cross sections of the ground-state charm hadrons: \(\textrm{D}^\textrm{0}\), \(\textrm{D}^+\), \(\textrm{D}_{\textrm{s}}^+\), \(\Lambda _{\textrm{c}}^{+}\), \({\Xi _{\textrm{c}}^{0}}\), and \({\textrm{J}}/\psi \). The \({\Xi _{\textrm{c}}^{0}}\) contribution was considered twice to account for the \({\Xi _{\textrm{c}}^{+}}\) cross section, as discussed above. The rapidity shift between the measured charm hadrons and the originating charm-quark pairs was accounted for by multiplying the sum of the hadron species by the aforementioned correction factor of 1.03. The correlations of the systematic uncertainties between charm-hadron species were treated in the same way as for the fragmentation fractions discussed above. The resulting charm cross section is
The measured cross section is compared to results in pp collisions at other centre-of-mass energies from ALICE [3, 4], PHENIX [51], and STAR [52] in Fig. 1 (right). In order to be directly comparable with the pp results, the cross section measured in p–Pb collisions is scaled down by a factor 1/A with \(A=208\), under the assumption that the charm production cross section in proton–nucleus collisions scales with the mass number of the larger nucleus due to the binary scaling of hard processes with the number of nucleon–nucleon collisions. The charm production cross section per nucleon in p–Pb collisions at \(\sqrt{s_\textrm{NN}} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \) is consistent with the one measured in pp collisions at \(\sqrt{s} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \), and is also in agreement with the general trend of results across different collision energies. The results are compared with pQCD calculations in the FONLL [47, 48] and Next-to-Next-to-Leading Order (NNLO) [53,54,55] approaches. As observed in previous works relating to the charm production cross section in pp collisions, while the experimental results tend to lie at the upper edge of the theoretical uncertainties, the pQCD approaches describe the evolution of the \({\textrm{c}}\overline{\textrm{c}}\) production cross section as a function of \(\sqrt{s_{\textsc {NN}}}\).
The \(p_\textrm{T}\)-integrated nuclear modification factor \(R_\textrm{pPb}\) for \({\textrm{c}}\overline{\textrm{c}}\) production in p–Pb collisions at \(\sqrt{s_\textrm{NN}} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \) was computed from the total charm cross section discussed above. The \(R_\textrm{pPb}\) is calculated from the measured cross sections in p–Pb and pp collisions at the same centre-of-mass energy as
with \(A=208\) as the Pb-nucleus mass number. A value of \(R_\textrm{pPb} =1\) would imply that there is no modification of the total yield at midrapidity in p–Pb collisions with respect to pp collisions. The correction factor \(\alpha _y=1.01\) accounts for the different y coverage of the measurements in pp and p–Pb collisions due to the centre-of-mass rapidity shift between the two collision systems. It was estimated from FONLL calculations of the \(\textrm{d}\sigma /\textrm{d}y\) distribution of charm quarks. The uncertainties on the measured charm cross sections in pp and p–Pb collisions due to the branching ratios and the contribution from \(\Omega _\textrm{c}^{0}\) baryons were treated as fully correlated in the ratio, and the rapidity shape uncertainties as partially correlated. All other sources of uncertainty were propagated as fully uncorrelated.
The \(R_\textrm{pPb}\) for \({\textrm{c}}\overline{\textrm{c}}\) pair production at midrapidity at \(\sqrt{s_\textrm{NN}} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \) is shown in Fig. 2 and is
which is consistent with unity within 1\(\sigma \). This implies that the overall charm production rate in p–Pb collisions is consistent with a binary scaling of pp collisions without any significant modification from nuclear effects.
The result is also compared with those for \(\textrm{D}^\textrm{0}\) mesons [28], \(\Lambda _{\textrm{c}}^{+}\) baryons [6], and \({\textrm{J}}/\psi \) mesons [39]. The \(p_\textrm{T}\)-integrated \(R_\textrm{pPb}\) values are consistent within uncertainties between all species, implying that the individual particle species’ production rates at midrapidity are similarly unaffected by nuclear effects. The results are compared with calculations of the \(p_\textrm{T}\)-integrated \(\textrm{D}^\textrm{0}\) \(R_\textrm{pPb}\) using both the EPPS16 [43] and nCTEQ15\(_{\textrm{rwHF}}\) [56, 57] nPDF sets. In the latter case, the nCTEQ15 predictions are modified by a Bayesian reweighting, constrained to heavy-flavour hadron measurements from the LHC [57]. The predictions of the \(p_\textrm{T}\)-integrated \(R_\textrm{pPb}\) are considered to be equivalent for all charm-hadron species, as only the parton distribution functions are varied between the pp and p–Pb calculations, and the \(p_\textrm{T}\)-integrated fragmentation functions cancel in the ratios. The calculation using EPPS16 is consistent with unity within large uncertainties, while the nCTEQ15 nPDFs predict a slight suppression of charm production. The experimental results are described by both of the calculations within uncertainties. In the case of calculations with nCTEQ15, the experimental results are on the upper edge of the model uncertainty.
4 Summary
The total charm production cross section per unit of rapidity and charm-hadron fragmentation fractions were measured for the first time at midrapidity in p–Pb collisions at \(\sqrt{s_\textrm{NN}} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} \). The fragmentation fractions for all charm-hadron species were found to be consistent for pp and p–Pb collisions at the same collision energy, indicating that there is no significant modification of the charm-quark hadronisation process due to the different hadronic collisions’ system sizes. The total charm production cross section per nucleon in p–Pb collisions was measured to be consistent with pp collisions at the same centre-of-mass energy within the uncertainties, meaning that the assumption of binary scaling holds for charm-quark production. This was further confirmed by examining the nuclear modification factor \(R_\textrm{pPb}\), which has a value consistent with unity at midrapidity and compatible with both the individual \(p_\textrm{T}\)-integrated charm-hadron \(R_\textrm{pPb}\) results and with model calculations including nuclear modification of the PDFs. Future measurements during Run 3 of the LHC will allow for a reduction of the uncertainties possibly reducing the reliance on model-dependent extrapolations. A particular focus in this field will be the determination of charm and beauty fragmentation fractions in Pb–Pb collisions, in order to precisely quantify modification of hadronisation driven by the formation and presence of a quark–gluon plasma, like modification of the abundance of different charm-hadron species [58,59,60,61] and quarkonium states [62, 63]. Future measurements in p–Pb collisions with higher precision will further allow significant constraints to be drawn on nuclear modifications of the PDFs.
Data Availability Statement
This manuscript has associated data in a data repository. [Author’s comment: This manuscript has associated data in a HEPData repository at https://www.hepdata.net/record/ins2789567.]
Code Availability Statement
This manuscript has no associated code/software. [Author’s comment: Code/software sharing is not applicable to this paper.]
References
J.C. Collins, D.E. Soper, G.F. Sterman, Factorization of hard processes in QCD. Adv. Ser. Direct. High Energy Phys. Ser. 5, 1–91 (1989). https://doi.org/10.1142/9789814503266_0001. arXiv:hep-ph/0409313
S. Catani, M. Ciafaloni, F. Hautmann, High-energy factorization and small x heavy flavor production. Nucl. Phys. B Ser. 366, 135–188 (1991). https://doi.org/10.1016/0550-3213(91)90055-3
ALICE Collaboration, S. Acharya et al., Charm-quark fragmentation fractions and production cross section at midrapidity in pp collisions at the LHC. Phys. Rev. D 105, L011103 (2022). https://doi.org/10.1103/PhysRevD.105.L011103. arXiv:2105.06335 [nucl-ex]
ALICE Collaboration, S. Acharya et al., Charm production and fragmentation fractions at midrapidity in pp collisions at \( \sqrt{\rm s} \) = 13 TeV. JHEP 12, 086 (2023). https://doi.org/10.1007/JHEP12(2023)086. arXiv:2308.04877 [hep-ex]
M. Lisovyi, A. Verbytskyi, O. Zenaiev, Combined analysis of charm-quark fragmentation-fraction measurements. Eur. Phys. J. C Ser. 76, 397 (2016). https://doi.org/10.1140/epjc/s10052-016-4246-y. arXiv:1509.01061 [hep-ex]
ALICE Collaboration, S. Acharya et al., First measurement of \(\rm \Lambda _{c}^{+}\) production down to \({p_{\rm T}}=0\) in pp and \({\text{p--Pb}}\) collisions at \({\sqrt{s_{\rm NN}} = 5.02~ \text{ TeV }} \). Phys. Rev. C107, 064901 (2023). https://doi.org/10.1103/PhysRevC.107.064901. arXiv:2211.14032 [nucl-ex]
ALICE Collaboration, S. Acharya et al., Measurement of the production cross section of prompt \(\Xi ^0_{\rm c}\) baryons in p–Pb collisions at \(\sqrt{s_{\rm NN}}=5.02\) TeV. arXiv:2405.14538 [nucl-ex]
LHCb Collaboration, R. Aaij et al., Prompt \(\Lambda _{\text{ c }}^{+}\) production in \({\text{ p--Pb }}\) collisions at \({\sqrt{s_{{\rm NN}}} = 5.02~\text{ TeV }} \). JHEP 02, 102 (2019). https://doi.org/10.1007/JHEP02(2019)102. arXiv:1809.01404 [hep-ex]
LHCb Collaboration, R. Aaij et al., Measurement of \(\Xi _{{\rm c}}^{+}\) production in \({\text{ p--Pb }}\) collisions at \(\sqrt{s_{\text{ NN }}} =8.16\) TeV at LHCb. Phys. Rev. C 109, 044901 (2024). https://doi.org/10.1103/PhysRevC.109.044901. arXiv:2305.06711 [hep-ex]
LHCb Collaboration, R. Aaij et al., Measurement of the \(\Lambda _{\text{ c }}^{+}\) to D\(^{0}\) production ratio in peripheral \(\text{ Pb--Pb } \) collisions at \( \sqrt{s_{\rm NN}} \) = 5.02 TeV. JHEP 06, 132 (2023). https://doi.org/10.1007/JHEP06(2023)132. arXiv:2210.06939 [hep-ex]. [Erratum: JHEP 05, 021 (2024)]
CMS Collaboration, A.M. Sirunyan et al., Production of \(\Lambda _{{\rm c}}^+\) baryons in proton-proton and lead-lead collisions at \(\sqrt{s_{\rm NN}}=\) 5.02 TeV. Phys. Lett. B 803, 135328, (2020).https://doi.org/10.1016/j.physletb.2020.135328. arXiv:1906.03322 [hep-ex]
PHENIX Collaboration, A. Adare et al., Heavy-quark production in \({{\rm pp}}\) and energy loss and flow of heavy quarks in Au + Au collisions at \({\sqrt{s_{\text{ NN }}}}=200\) GeV. Phys. Rev. C 84, 044905 (2011). https://doi.org/10.1103/PhysRevC.84.044905
STAR Collaboration, L. Adamczyk et al., Measurements of \(\rm D^0\) and \({{{\rm D}}^{*+}}\) production in \({{\rm pp}}\) collisions at \(\sqrt{s} = 200 {\rm GeV}\). Phys. Rev. D 86, 072013 (2012). https://doi.org/10.1103/PhysRevD.86.072013
STAR Collaboration, J. Adams et al., Open charm yields in d–Au collisions at \(\sqrt{s_{NN}} =200 {{\rm GeV}}\). Phys. Rev. Lett. 94, 062301 (2005). https://doi.org/10.1103/PhysRevLett.94.062301
C. Bierlich et al., Open charm production cross section from combined LHC experiments in pp collisions at \(\sqrt{s} = 5.02\) TeV. Eur. Phys. J. Plus 139, 593 (2024). https://doi.org/10.1140/epjp/s13360-024-05355-0. arXiv:2311.11426 [hep-ph]
ALICE Collaboration, S. Acharya et al., xPseudorapidity distributions of charged particles as a function of mid- and forward rapidity multiplicities in pp collisions at \(\sqrt{s}\) = 5.02, 7 and 13 TeV. Eur. Phys. J. C81, 630 (2021). https://doi.org/10.1140/epjc/s10052-021-09349-5. arXiv:2009.09434 [nucl-ex]
ALICE Collaboration, B. Abelev et al., Pseudorapidity density of charged particles in p–Pb collisions at \(\sqrt{s_{\rm NN}}=5.02\) TeV. Phys. Rev. Lett. 110, 032301 (2013). https://doi.org/10.1103/PhysRevLett.110.032301. arXiv:1210.3615 [nucl-ex]
ALICE Collaboration, S. Acharya et al., Observation of a multiplicity dependence in the \(p_{\rm T}\)-differential charm baryon-to-meson ratios in proton–proton collisions at \(\sqrt{s}=13\) TeV. Phys. Lett. B829, 137065 (2022). https://doi.org/10.1016/j.physletb.2022.137065. arXiv:2111.11948 [nucl-ex]
LHCb Collaboration, R. Aaij et al., Measurement of \(b\) hadron fractions in 13 TeV pp collisions. Phys. Rev. D 100, 031102 (2019). https://doi.org/10.1103/PhysRevD.100.031102. arXiv:1902.06794 [hep-ex]
LHCb Collaboration, R. Aaij et al., Study of the production of \(\Lambda _{{\rm b}}^0\) and \(\overline{{\rm B}}^0\) hadrons in pp collisions and first measurement of the \(\Lambda _{{\rm b}}^0\rightarrow {{\rm J}}/\psi {\rm pK}^-\) branching fraction. Chin. Phys. C 40, 011001 (2016). https://doi.org/10.1088/1674-1137/40/1/011001. arXiv:1509.00292 [hep-ex]
LHCb Collaboration, R. Aaij et al., Enhanced production of \(\Lambda _{{\rm b}}^{0}\) baryons in high-multiplicity pp collisions at \(\sqrt{s} = 13\) TeV. Phys. Rev. Lett. 132, 081901 (2024). https://doi.org/10.1103/PhysRevLett.132.081901. arXiv:2310.12278 [hep-ex]
LHCb Collaboration, R. Aaij et al., Measurement of \({{\rm B}}^+\), \({{\rm B}}^0\) and \(\Lambda _{\rm b}^0\) production in \({\text{ p--Pb }}\) collisions at \(\sqrt{s_{\text{ NN }}} =8.16\;{\rm TeV}\). Phys. Rev. D 99, 052011 (2019). https://doi.org/10.1103/PhysRevD.99.052011. arXiv:1902.05599 [hep-ex]
ALICE Collaboration, S. Acharya et al., Study of flavor dependence of the baryon-to-meson ratio in proton-proton collisions at \(\sqrt{s}=13\) TeV. Phys. Rev. D 108, 112003 (2023). https://doi.org/10.1103/PhysRevD.108.112003. arXiv:2308.04873 [hep-ex]
M. Arneodo, Nuclear effects in structure functions. Phys. Rep. Ser. 240, 301–393 (1994). https://doi.org/10.1016/0370-1573(94)90048-5
M. Klasen, H. Paukkunen, Nuclear parton distribution functions after the first decade of LHC data. Annu. Rev. Nucl. Part. Sci. (2024). https://doi.org/10.1146/annurev-nucl-102122-022747. arXiv:2311.00450 [hep-ph]
M. Lev, B. Petersson, Nuclear effects at large transverse momentum in a QCD parton model. Z. Phys. C Ser. 21, 155 (1983). https://doi.org/10.1007/BF01648792
B.Z. Kopeliovich, J. Nemchik, A. Schafer, A.V. Tarasov, Cronin effect in hadron production off nuclei. Phys. Rev. Lett. Ser. 88, 232303 (2002). https://doi.org/10.1103/PhysRevLett.88.232303. arXiv:hep-ph/0201010
ALICE Collaboration, S. Acharya et al., Measurement of prompt \({{\rm D}}^{0}\), \({\rm D}^+\), \({\rm D}^{*+}\), and \({\rm D}_{s}^+ \) production in \({\text{ p--Pb }}\) collisions at \({\sqrt{s_{\rm NN}} = 5.02~\text{ TeV }} \). JHEP12, 092 (2019). https://doi.org/10.1007/JHEP12(2019)092. arXiv:1906.03425 [nucl-ex]
P. Bozek, Collective flow in p–Pb and d–Pb collisions at TeV energies. Phys. Rev. C Ser. 85, 014911 (2012). https://doi.org/10.1103/PhysRevC.85.014911. arXiv:1112.0915 [hep-ph]
R.D. Weller, P. Romatschke, One fluid to rule them all: viscous hydrodynamic description of event-by-event central pp, p–Pb and Pb–Pb collisions at \(\sqrt{s}=5.02\) TeV. Phys. Lett. B 774, 351–356 (2017). https://doi.org/10.1016/j.physletb.2017.09.077. arXiv:1701.07145 [nucl-th]
J. Song, H.-H. Li, F.-L. Shao, New feature of low \(p_{{\rm T}}\) charm quark hadronization in pp collisions at \(\sqrt{s}=7\) TeV. Eur. Phys. J. C Ser. 78, 344 (2018). https://doi.org/10.1140/epjc/s10052-018-5817-x. arXiv:1801.09402 [hep-ph]
V. Minissale, S. Plumari, V. Greco, Charm hadrons in pp collisions at LHC energy within a coalescence plus fragmentation approach. Phys. Lett. B Ser. 821, 136622 (2021). https://doi.org/10.1016/j.physletb.2021.136622. arXiv:2012.12001 [hep-ph]
A. Beraudo, A. De Pace, D. Pablos, F. Prino, M. Monteno, M. Nardi, Heavy-flavor transport and hadronization in pp collisions. Phys. Rev. D Ser. 109, L011501 (2024). https://doi.org/10.1103/PhysRevD.109.L011501. arXiv:2306.02152 [hep-ph]
ALICE Collaboration, B. Abelev et al., Multiplicity dependence of pion, kaon, proton and lambda production in p–Pb collisions at \(\sqrt{s_{\rm NN}}\) = 5.02 TeV. Phys. Lett. B 728, 25–38 (2014). https://doi.org/10.1016/j.physletb.2013.11.020. arXiv:1307.6796 [nucl-ex]
J. Altmann et al., Towards the understanding of heavy quarks hadronization: from leptonic to heavy-ion collisions. arXiv:2405.19137 [hep-ph]
ALICE Collaboration, K. Aamodt et al., The ALICE experiment at the CERN LHC. JINST3, S08002 (2008). https://doi.org/10.1088/1748-0221/3/08/S08002
ALICE Collaboration, B. Abelev et al., Measurement of visible cross sections in proton–lead collisions at \(\sqrt{s_{\rm NN}}\) = 5.02 TeV in van der Meer scans with the ALICE detector. JINST 9, P11003 (2014). https://doi.org/10.1088/1748-0221/9/11/P11003. arXiv:1405.1849 [nucl-ex]
ALICE Collaboration, S. Acharya et al., \(\Lambda ^{+}_{{\rm c}}\) production and baryon-to-meson ratios in pp and p–Pb collisions at \(\sqrt{s_{\rm NN}}=5.02\,\text{ TeV } \) at the LHC. Phys. Rev. Lett. 127, 202301 (2021). https://doi.org/10.1103/PhysRevLett.127.202301. arXiv:2011.06078 [nucl-ex]
ALICE Collaboration, S. Acharya et al., Inclusive, prompt and non-prompt J/\(\psi \) production at midrapidity in p-Pb collisions at \( \sqrt{s_{\rm NN}} \) = 5.02 TeV. JHEP 06, 011 (2022). https://doi.org/10.1007/JHEP06(2022)011. arXiv:2105.04957 [nucl-ex]
S. Frixione, P. Nason, G. Ridolfi, A Positive-weight next-to-leading-order Monte Carlo for heavy flavour hadroproduction. JHEP Ser. 09, 126 (2007). https://doi.org/10.1088/1126-6708/2007/09/126. arXiv:0707.3088 [hep-ph]
T. Sjostrand, S. Mrenna, P.Z. Skands, PYTHIA 6.4 physics and manual. JHEP 05, 026 (2006). https://doi.org/10.1088/1126-6708/2006/05/026. arXiv:hep-ph/0603175
S. Dulat et al., New parton distribution functions from a global analysis of quantum chromodynamics. Phys. Rev. D Ser. 93, 033006 (2016). https://doi.org/10.1103/PhysRevD.93.033006. arXiv:1506.07443 [hep-ph]
K.J. Eskola, P. Paakkinen, H. Paukkunen, C.A. Salgado, EPPS16: nuclear parton distributions with LHC data. Eur. Phys. J. C Ser. 77, 163 (2017). https://doi.org/10.1140/epjc/s10052-017-4725-9. arXiv:1612.05741 [hep-ph]
H. Li, F. Shao, J. Song, R. Wang, Production of single-charm hadrons by quark combination mechanism in \({\text{ p--Pb }}\) collisions at \(\sqrt{s_{\rm NN}}=5.02\) TeV. Phys. Rev. C 97, 064915 (2018). https://doi.org/10.1103/PhysRevC.97.064915. arXiv:1712.08921 [hep-ph]
H. Li, F. Shao, J. Song, Production of light-flavor and single-charmed hadrons in pp collisions at \({\sqrt{s} = 5.02~\text{ TeV }}\) in an equal-velocity quark combination model. Chin. Phys. C 45, 113105 (2021). https://doi.org/10.1088/1674-1137/ac1ef9. arXiv:2103.14900 [hep-ph]
C. Bierlich, G. Gustafson, L. Lönnblad, H. Shah, The Angantyr model for heavy-ion collisions in PYTHIA8. JHEP Ser. 10, 134 (2018). https://doi.org/10.1007/JHEP10(2018)134. arXiv:1806.10820 [hep-ph]
M. Cacciari, S. Frixione, N. Houdeau, M.L. Mangano, P. Nason, G. Ridolfi, Theoretical predictions for charm and bottom production at the LHC. JHEP Ser. 10, 137 (2012). https://doi.org/10.1007/JHEP10(2012)137. arXiv:1205.6344 [hep-ph]
M. Cacciari, M.L. Mangano, P. Nason, Gluon PDF constraints from the ratio of forward heavy-quark production at the LHC at \(\sqrt{s}=7\) and 13 TeV. Eur. Phys. J. C Ser. 75, 610 (2015). https://doi.org/10.1140/epjc/s10052-015-3814-x. arXiv:1507.06197 [hep-ph]
P.M. Nadolsky, H.-L. Lai, Q.-H. Cao, J. Huston, J. Pumplin, D. Stump, W.-K. Tung, C.P. Yuan, Implications of CTEQ global analysis for collider observables. Phys. Rev. D Ser. 78, 013004 (2008). https://doi.org/10.1103/PhysRevD.78.013004. arXiv:0802.0007 [hep-ph]
M. Cacciari, M. Greco, P. Nason, The \({p_{\rm T}}\) spectrum in heavy flavor hadroproduction. JHEP Ser. 05, 007 (1998). https://doi.org/10.1088/1126-6708/1998/05/007. arXiv:hep-ph/9803400
PHENIX Collaboration, A. Adare et al., Heavy quark production in pp and energy loss and flow of heavy quarks in Au–Au collisions at \(\sqrt{s_{\rm NN}}=200\) GeV. Phys. Rev. C 84, 044905 (2011). https://doi.org/10.1103/PhysRevC.84.044905. arXiv:1005.1627 [nucl-ex]
STAR Collaboration, L. Adamczyk et al., Measurements of \({{{\rm D}}^0}\) and \({{{\rm D}}^{*+}}\) production in pp collisions at \(\sqrt{s} = 200\) GeV. Phys. Rev. D 86, 072013 (2012). https://doi.org/10.1103/PhysRevD.86.072013. arXiv:1204.4244 [nucl-ex]
D. d’Enterria, A.M. Snigirev, Triple parton scatterings in high-energy proton–proton collisions. Phys. Rev. Lett. Ser. 118, 122001 (2017). https://doi.org/10.1103/PhysRevLett.118.122001. arXiv:1612.05582 [hep-ph]
D. d’Enterria, A.M. Snigirev, Triple-parton scatterings in proton–nucleus collisions at high energies. Eur. Phys. J. C Ser. 78, 359 (2018). https://doi.org/10.1140/epjc/s10052-018-5687-2. arXiv:1612.08112 [hep-ph]
M. Czakon, P. Fiedler, A. Mitov, Total top-quark pair-production cross section at hadron colliders through \(O(\alpha ^4_{{\rm s}})\). Phys. Rev. Lett. Ser. 110, 252004 (2013). https://doi.org/10.1103/PhysRevLett.110.252004. arXiv:1303.6254 [hep-ph]
K. Kovarik et al., nCTEQ15—global analysis of nuclear parton distributions with uncertainties in the CTEQ framework. Phys. Rev. D Ser. 93, 085037 (2016). https://doi.org/10.1103/PhysRevD.93.085037. arXiv:1509.00792 [hep-ph]
A. Kusina, J.-P. Lansberg, I. Schienbein, H.-S. Shao, Gluon shadowing in heavy-flavor production at the LHC. Phys. Rev. Lett. Ser. 121, 052004 (2018). https://doi.org/10.1103/PhysRevLett.121.052004. arXiv:1712.07024 [hep-ph]
ALICE Collaboration, S. Acharya et al., Measurement of prompt \({{{\rm D}}_s^+}\)-meson production and azimuthal anisotropy in Pb–Pb collisions at \(\sqrt{s_{\rm NN}}=5.02\) TeV. Phys. Lett. B 827, 136986 (2022). https://doi.org/10.1016/j.physletb.2022.136986. arXiv:2110.10006 [nucl-ex]
STAR Collaboration, J. Adam et al., Observation of \({\rm D}_{\rm s}^{\pm }/{{\rm D}}^0\) enhancement in Au–Au collisions at \(\sqrt{s_{\rm NN}}\) = 200 GeV. Phys. Rev. Lett. 127, 092301 (2021). https://doi.org/10.1103/PhysRevLett.127.092301. arXiv:2101.11793 [hep-ex]
ALICE Collaboration, S. Acharya et al., Constraining hadronization mechanisms with \(\rm \Lambda _{\rm c}^{+}\)/D\(^0\) production ratios in Pb–Pb collisions at \(\sqrt{s_{\rm NN}} = 5.02\) TeV. Phys. Lett. B 839, 137796 (2023). https://doi.org/10.1016/j.physletb.2023.137796. arXiv:2112.08156 [nucl-ex]
STAR Collaboration, J. Adam et al., First measurement of \(\Lambda _{{\rm c}}\) baryon production in Au–Au collisions at \(\sqrt{s_{\rm NN}}\) = 200 GeV. Phys. Rev. Lett. 124, 172301 (2020). https://doi.org/10.1103/PhysRevLett.124.172301. arXiv:1910.14628 [nucl-ex]
ALICE Collaboration, S. Acharya et al., \(\psi \)(2S) suppression in Pb–Pb collisions at the LHC. Phys. Rev. Lett. 132, 042301 (2024). https://doi.org/10.1103/PhysRevLett.132.042301. arXiv:2210.08893 [nucl-ex]
ALICE Collaboration, S. Acharya et al., Measurements of inclusive J/\(\psi \) production at midrapidity and forward rapidity in \(\text{ Pb--Pb } \) collisions at \({\sqrt{s_{\rm NN}} = 5.02~\text{ TeV }} \). Phys. Lett. B 849, 138451 (2024). https://doi.org/10.1016/j.physletb.2024.138451. arXiv:2303.13361 [nucl-ex]
Acknowledgements
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. I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation (ANSL), State Committee of Science and World Federation of Scientists (WFS), Armenia; Austrian Academy of Sciences, Austrian Science Fund (FWF): [M 2467-N36] and Nationalstiftung für Forschung, Technologie und Entwicklung, Austria; Ministry of Communications and High Technologies, National Nuclear Research Center, Azerbaijan; Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Financiadora de Estudos e Projetos (Finep), Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) and Universidade Federal do Rio Grande do Sul (UFRGS), Brazil; Bulgarian Ministry of Education and Science, within the National Roadmap for Research Infrastructures 2020-2027 (object CERN), Bulgaria; Ministry of Education of China (MOEC), Ministry of Science & Technology of China (MSTC) and National Natural Science Foundation of China (NSFC), China; Ministry of Science and Education and Croatian Science Foundation, Croatia; Centro de Aplicaciones Tecnológicas y Desarrollo Nuclear (CEADEN), Cubaenergía, Cuba; Ministry of Education, Youth and Sports of the Czech Republic, Czech Republic; The Danish Council for Independent Research | Natural Sciences, the VILLUM FONDEN and Danish National Research Foundation (DNRF), Denmark; Helsinki Institute of Physics (HIP), Finland; Commissariat à l’Energie Atomique (CEA) and Institut National de Physique Nucléaire et de Physique des Particules (IN2P3) and Centre National de la Recherche Scientifique (CNRS), France; Bundesministerium für Bildung und Forschung (BMBF) and GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany; General Secretariat for Research and Technology, Ministry of Education, Research and Religions, Greece; National Research, Development and Innovation Office, Hungary; Department of Atomic Energy Government of India (DAE), Department of Science and Technology, Government of India (DST), University Grants Commission, Government of India (UGC) and Council of Scientific and Industrial Research (CSIR), India; National Research and Innovation Agency-BRIN, Indonesia; Istituto Nazionale di Fisica Nucleare (INFN), Italy; Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT) and Japan Society for the Promotion of Science (JSPS) KAKENHI, Japan; Consejo Nacional de Ciencia (CONACYT) y Tecnología, through Fondo de Cooperación Internacional en Ciencia y Tecnología (FONCICYT) and Dirección General de Asuntos del Personal Academico (DGAPA), Mexico; Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO), Netherlands; The Research Council of Norway, Norway; Pontificia Universidad Católica del Perú, Peru; Ministry of Science and Higher Education, National Science Centre and WUT ID-UB, Poland; Korea Institute of Science and Technology Information and National Research Foundation of Korea (NRF), Republic of Korea; Ministry of Education and Scientific Research, Institute of Atomic Physics, Ministry of Research and Innovation and Institute of Atomic Physics and Universitatea Nationala de Stiinta si Tehnologie Politehnica Bucuresti, Romania; Ministry of Education, Science, Research and Sport of the Slovak Republic, Slovakia; National Research Foundation of South Africa, South Africa; Swedish Research Council (VR) and Knut & Alice Wallenberg Foundation (KAW), Sweden; European Organization for Nuclear Research, Switzerland; Suranaree University of Technology (SUT), National Science and Technology Development Agency (NSTDA) and National Science, Research and Innovation Fund (NSRF via PMU-B B05F650021), Thailand; Turkish Energy, Nuclear and Mineral Research Agency (TENMAK), Turkey; National Academy of Sciences of Ukraine, Ukraine; Science and Technology Facilities Council (STFC), United Kingdom; National Science Foundation of the United States of America (NSF) and United States Department of Energy, Office of Nuclear Physics (DOE NP), United States of America. In addition, individual groups or members have received support from: Czech Science Foundation (grant no. 23-07499S), Czech Republic; European Research Council (grant no. 950692), European Union; ICSC-Centro Nazionale di Ricerca in High Performance Computing, Big Data and Quantum Computing, European Union-NextGenerationEU; Academy of Finland (Center of Excellence in Quark Matter) (Grant nos. 346327, 346328), Finland.
Author information
Authors and Affiliations
Consortia
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3.
About this article
Cite this article
ALICE Collaboration. Charm fragmentation fractions and \(\mathbf {c\overline{c}}\) cross section in p–Pb collisions at \(\mathbf {\sqrt{s _{NN}}=5.02 TeV}\). Eur. Phys. J. C 84, 1286 (2024). https://doi.org/10.1140/epjc/s10052-024-13394-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-024-13394-1