Measurement of ${\rm D^0}$, ${\rm D^+}$, ${\rm D^{*+}}$ and ${{\rm D^+}_s}$ production in pp collisions at $\mathbf{\sqrt{{\textit s}}~=~5.02~TeV}$ with ALICE

The measurements of the production of prompt ${\rm D^0}$, ${\rm D^+}$, ${\rm D^{*+}}$, and ${{\rm D^+}_s}$ mesons in proton--proton (pp) collisions at $\sqrt{s}=5.02$ TeV with the ALICE detector at the Large Hadron Collider (LHC) are reported. D mesons were reconstructed at mid-rapidity ($|y|<0.5$) via their hadronic decay channels ${\rm D}^0 \to {\rm K}^-\pi^+$, ${\rm D}^+\to {\rm K}^-\pi^+\pi^+$, ${\rm D}^{*+} \to {\rm D}^0 \pi^+ \to {\rm K}^- \pi^+ \pi^+$, ${\rm D^{+}_{s}\to \phi\pi^+\to K^{+} K^{-} \pi^{+}}$, and their charge conjugates. The production cross sections were measured in the transverse momentum interval $0<p_{\rm T}<36~\mathrm{GeV}/c$ for ${\rm D^0}$, $1<p_{\rm T}<36~\mathrm{GeV}/c$ for ${\rm D^+}$ and ${\rm D^{*+}}$, and in $2<p_{\rm T}<24~\mathrm{GeV}/c$ for ${{\rm D^+}_s}$ mesons. Thanks to the higher integrated luminosity, an analysis in finer $p_{\rm T}$ bins with respect to the previous measurements at $\sqrt{s}=7$ TeV was performed, allowing for a more detailed description of the cross-section $p_{\rm T}$ shape. The measured $p_{\rm T}$-differential production cross sections are compared to the results at $\sqrt{s}=7$ TeV and to four different perturbative QCD calculations. Its rapidity dependence is also tested combining the ALICE and LHCb measurements in pp collisions at $\sqrt{s}=5.02$ TeV. This measurement will allow for a more accurate determination of the nuclear modification factor in p-Pb and Pb-Pb collisions performed at the same nucleon-nucleon centre-of-mass energy.


Introduction
The study of the production of hadrons containing heavy quarks, i.e. charm and beauty, in proton-proton (pp) collisions at LHC energies is a sensitive test of Quantum Chromodynamics (QCD) calculations with the factorisation approach. In this scheme, the transverse momentum (p T ) differential production cross sections of hadrons containing charm or beauty quarks are calculated as a convolution of three terms: (i) the parton distribution functions (PDFs) of the incoming protons, (ii) the partonic scattering cross section, calculated as a perturbative series in powers of the strong coupling constant α s , and (iii) the fragmentation function, which parametrises the non-perturbative evolution of a heavy quark into a given species of heavy-flavour hadron. Factorisation is implemented in terms of the squared momentum transfer Q 2 (collinear factorisation) [1] or of the partonic transverse momentum k T [2]. At LHC energies, calculations based on collinear factorisation are available in the general-mass variable-flavournumber scheme, GM-VFNS [3][4][5][6], and in the fixed order plus next-to-leading logarithms approach, FONLL [7,8], both of them having next-to-leading order (NLO) accuracy with all-order resummation of next-to-leading logarithms. Within the k T -factorisation framework, heavy-flavour production crosssection calculations exist only at leading order (LO) approximation in α s [2,9,10]. These calculations describe within uncertainties the production cross sections of D and B mesons measured in pp and pp collisions in different kinematic regions at centre-of-mass energies from 0.2 to 13 TeV (see e.g. Ref. [11] and references therein). In the case of charm production, pQCD calculations have significantly larger uncertainties, dominated by the choice of the scales of the perturbative calculation (e.g. the factorisation and renormalisation scales), than the experimental ones [12][13][14][15][16][17][18][19][20][21][22][23]. However, as pointed out in Ref. [24], in the ratios of cross sections at different LHC energies and in different rapidity intervals the uncertainty due to choice of the factorisation and renormalisation scales becomes subdominant with respect to the uncertainty on the PDFs, thus making the measurement sensitive to the gluon PDF at small Bjorken-x values. A precise measurement of the D-meson production cross sections down to p T = 0 can therefore provide important constraints to pQCD calculations and to low-x gluon PDFs. Furthermore, D-meson measurements in pp collisions represent an essential reference for the study of effects induced by cold and hot strongly-interacting matter in the case of proton-nucleus and nucleus-nucleus collisions (see e.g. the recent reviews [11,25,26]).
In this article, the measurements of the p T -differential production cross sections of prompt D 0 , D + , D * + , and D + s mesons (as average of particles and anti-particles) in pp collisions at the centre-of-mass energy √ s = 5.02 TeV are reported together with their ratios. The measurements are performed at mid-rapidity (|y| < 0.5) in the transverse momentum intervals 0 < p T < 36 GeV/c for D 0 mesons, 1 < p T < 36 GeV/c for D + and D * + mesons, and 2 < p T < 24 GeV/c for D + s mesons. The p T -integrated D-meson production cross sections per unit of rapidity is also reported for each D-meson specie. The ratios of the D 0 , D + , and D * + -meson production cross sections measured at √ s = 7 TeV [27] and √ s = 5.02 TeV are presented as well, and compared to FONLL calculations. Finally, the ratios of D 0 -meson production cross sections at mid-and forward rapidity are also reported, using the measurements done at forward rapidity by the LHCb collaboration in pp collisions at √ s = 5.02 TeV [22].

Experimental apparatus and data sample
The ALICE experimental apparatus is composed of a set of detectors for particle reconstruction and identification at mid-rapidity, embedded in a large solenoidal magnet that provides a B = 0.5 T field parallel to the beams. It also includes a forward muon spectrometer and various forward and backward detectors for triggering and event characterisation. A complete description and an overview of their typical performance in pp, p-Pb, and Pb-Pb collisions is presented in Refs. [28,29].
The tracking and particle identification capabilities of the ALICE central barrel detectors were exploited to reconstruct the D-meson decay products at mid-rapidity. The Inner Tracking System (ITS), consisting 3 Data analysis 3

.1 Analysis with D-meson decay vertex reconstruction
The D mesons and their charge conjugates were reconstructed via the decay channels D 0 → K − π + (with branching ratio, BR = 3.89 ± 0.04%), D + → K − π + π + (BR = 8.98 ± 0.28%), D * + → D 0 π + → K − π + π + (BR = 2.63 ± 0.03%), and D + s → φ π + → K + K − π + (BR = 2.27 ± 0.08%) [30]. The analysis was based on the reconstruction of decay vertices displaced from the interaction vertex, exploiting the separation of a few hundred µm induced by the weak decays of D 0 , D + , and D + s mesons (cτ ≃ 123, 312, and 150 µm, respectively [30]). The D 0 , D + , and D + s candidates were built combining pairs or triplets of tracks with the proper charge, |η| < 0.8, p T > 0.3 GeV/c, at least 70 associated TPC space points, χ 2 /ndf < 2 in the TPC (where ndf is the number of degrees of freedom involved in the track fit procedure), and at least one hit in either of the two layers of the SPD. The D * + candidates were defined by the combination of D 0 candidates with tracks reconstructed with at least two points in the ITS, including at least one in the SPD, and p T > 80 MeV/c. As a consequence of these track selection criteria, the acceptance for D mesons decreases rapidly for |y| > 0.5 at low p T and for |y| > 0.8 for p T > 5 GeV/c. Therefore, only D-meson candidates within a fiducial acceptance region, |y| < y fid (p T ), were selected. The y fid (p T ) factor was defined as a second-order polynomial function, increasing from 0.5 to 0.8 in the transverse momentum range 0 < p T < 5 GeV/c, and a constant term, y fid = 0.8, for p T > 5 GeV/c.
In order to reduce the combinatorial background and to increase the signal-over-background ratio (S/B), geometrical selections on the D 0 , D + , and D + s -meson decay topology were applied. In the D * + → D 0 π + case, the decay vertex cannot be resolved from the primary vertex and geometrical selections were applied on the secondary vertex topology of the produced D 0 mesons. The selection requirements, tuned to provide a large statistical significance for the signal and to keep the selection efficiency as high as possible, were mainly based on the displacement of the tracks from the primary vertex (d 0 ), the distance between the D-meson decay vertex and the primary vertex (decay length, L), and the pointing of the reconstructed D-meson momentum to the primary vertex. Additional selection criteria, already introduced in Refs. [27,31], were applied to D + and D + s candidates. These selections reject both combinatorial background and D mesons from beauty-hadron decays (up to 50% at high p T ), denoted as "feed-down" in the following. For the D + s -candidate selection, one of the two pairs of opposite-sign tracks was required to be within ±10 MeV/c 2 with respect to the PDG world average of the φ meson [30]. D-meson production in pp collisions at √ s = 5.02 TeV ALICE Collaboration Further reduction of the combinatorial background was obtained by applying particle identification (PID) to the decay tracks, except for the soft-pion track coming from D * + → D 0 π + decays. Pions and kaons were identified requiring three standard deviations (3 σ ) compatibility between the measured and the expected signals for both the TPC dE/dx and the time-of-flight. Tracks without TOF hits were identified using only the TPC information with a 3 σ selection, except for the decay products of D + s candidates with p T < 6 GeV/c, for which a 2 σ selection was needed.
The D-meson raw yields, including both particles and antiparticles, were obtained from binned maximum likelihood fits to the invariant-mass (M) distributions of D 0 , D + , and D + s candidates and to the mass difference ∆M = M(Kππ) − M(Kπ) distributions of D * + candidates, in the transverse-momentum intervals 0.5 < p T < 36 GeV/c for D 0 mesons, 1 < p T < 36 GeV/c for D + and D * + mesons, and 2 < p T < 24 GeV/c for D + s mesons. The signal extraction was performed in finer p T bins with respect the previous measurements at √ s = 7 TeV [27], allowing for a more detailed description of the crosssection p T shape. The fit function was composed by a Gaussian for the description of the signal and by an exponential term for the background of D 0 , D + , and D + s candidates, and by a threshold function for D * + candidates [27]. The contribution of signal candidates present in the invariant-mass distribution of the D 0 meson with the wrong decay-particle mass assignment (reflections) was modelled from the fit with the sum of two Gaussian functions of the simulated reflected invariant-mass distributions and included in the fit function (about 2% − 3% of the raw signal, depending on p T ). For the M(KKπ) distribution, an additional Gaussian was used to describe the signal of the decay D + → K + K − π + , with a branching ratio of (9.51 ± 0.34) · 10 −3 [30], present on the left side of the D + s -meson signal. Figure 1 shows the invariant mass (mass-difference) distributions together with the result of the fits, in 1.5 < p T < 2 GeV/c, 16 < p T < 24 GeV/c, 7 < p T < 7.5 GeV/c, and 3 < p T < 4 GeV/c intervals for D 0 , D + , D * + , and D + s candidates, respectively. The statistical significance of the observed signals, S/ (S + B), varies from 4 to 28, depending on the meson species and on the p T interval. The S/B values obtained applying the selections described above are 0.01-1.85 for D 0 , 0.5-2.2 for D + , 0.3-4.2 for D * + , and 0.3-2.2 for D + s mesons, depending on p T .
The p T -differential cross section of prompt D mesons in each p T interval was computed as: The raw yield values (sum of particles and antiparticles, N D+D,raw ) were divided by a factor of two and multiplied by the prompt fraction f prompt to obtain the charged-averaged yields of prompt D mesons. Furthermore, they were divided by the acceptance-times-efficiency of prompt D mesons (Acc × ε) prompt , the BR of the decay channel, the width of the p T interval (∆p T ), the correction factor for the rapidity coverage c ∆y , and the integrated luminosity L int = N ev /σ MB , where N ev is the number of analysed events and σ MB = (50.9 ± 0.9) mb is the cross section for the MB trigger condition [32].
The (Acc × ε) correction was obtained simulating pp collisions with the PYTHIA 6.4.25 event generator [33] (Perugia-11 tune [34]), and propagating the generated particles through the apparatus using GEANT3 [35]. Each simulated PYTHIA pp event contained a cc or bb pair, and D mesons were forced to decay into the hadronic channels of interest for the analysis. The luminous region distribution and the conditions of all the ALICE detectors in terms of active channels, gain, noise level and alignment, and their evolution with time during the data taking, were taken into account in the simulations. Figure 2 shows the (Acc× ε) as a function of p T for prompt and feed-down D 0 , D + , D * + , and D + s mesons within the fiducial acceptance region. The average larger displacement from the primary vertex of beauty hadrons due to their long lifetime (cτ ≈ 500 µm [30]) results in a more efficient selection of feed-down D mesons compared to prompt D mesons in most of the p T intervals.
The correction factor for the rapidity acceptance c ∆y was computed with the PYTHIA 6.4.25 event gener-D-meson production in pp collisions at √ s = 5.02 TeV ALICE Collaboration   s candidates and charge conjugates in 1.5 < p T < 2 GeV/c, 16 < p T < 24 GeV/c, 7 < p T < 7.5 GeV/c, and 3 < p T < 4 GeV/c intervals, respectively. The blue solid lines show the total fit functions as described in the text and the red dashed lines are the combinatorial-background terms. In case of D 0 , the grey dashed line represents the combinatorial background with the contribution of the reflections. The values of the mean (µ) and the width (σ ) of the signal peak are reported together with the signal counts (S) and the signal over background ratio (S/B) in the mass interval (µ − 3σ , µ + 3σ ). The reported uncertainties are only the statistical uncertainties from the fit. ator with Perugia-11 tune. It was defined as the ratio between the generated D-meson yield in ∆y = 2 y fid , and that in |y| < 0.5. It was checked that calculations of the c ∆y correction factor based on FONLL pQCD calculations [8] or on the assumption of uniform D-meson rapidity distribution in |y| < y fid would give the same result, because both in PYTHIA and in FONLL the D-meson yield is uniform within 1% in the range |y| < 0.8.
The f prompt fraction was calculated similarly to previous measurements (see e.g. Refs. [27, 31]) using the beauty-hadron production cross sections from FONLL calculations [7,36], the beauty hadron → D + X decay kinematics from the EvtGen package [37], and the efficiencies for feed-down D mesons reported in Fig. 2. The values of f prompt range between 0.8 and 0.96 depending on D-meson species and p T . D-meson production in pp collisions at √ s = 5.02 TeV ALICE Collaboration

Analysis without D-meson decay vertex reconstruction
A different analysis method, not based on geometrical selections of the displaced decay-vertex topology, was developed for the two-body decay D 0 → K − π + (and its charge conjugate) in order to extend the measurement of the cross section down to p T = 0 [19]. Indeed, the poor track impact parameter resolution at very low p T and the small Lorentz boost limit the effectiveness of the selections based on the displaced decay-vertex topology. Furthermore, geometrical selections based on the displacement of the D 0 -meson decay vertex tend to enhance the contribution of feed-down D mesons, increasing the related systematic uncertainty. This alternative analysis technique is mainly based on particle identification and on the estimation and subtraction of the combinatorial background.
The D 0 candidates were formed combining pairs of kaons and pions tracks with opposite charge sign (UnLike Sign, ULS), |η| < 0.8, and p T > 0.3 GeV/c. Track selection and pion and kaon identification were performed with the same strategy used in the analysis with decay-vertex reconstruction described in Section 3.1. The resulting D 0 and D 0 candidates were selected by applying the same fiducial acceptance D-meson production in pp collisions at √ s = 5.02 TeV ALICE Collaboration  selection |y| < y fid (p T ) adopted for the analysis with decay-vertex reconstruction. The invariant-mass distribution of Kπ pairs was obtained in fourteen transverse momentum intervals, in the range 0 < p T < 12 GeV/c. The background distribution was estimated with the track-rotation technique. For each D 0 (and D 0 ) candidate, up to 19 combinatorial-background-like candidates were created by rotating the kaon track by different angles in the range between π 10 and 19π 10 radians in azimuth. The left hand panel of Fig. 3 shows the invariant-mass distribution of ULS Kπ pairs together with that of the background estimated with the track-rotation technique in the interval 0 < p T < 0.5 GeV/c.
After subtracting the background distribution from the ULS Kπ invariant-mass distribution, the D 0meson raw signal (sum of particle and antiparticle contributions) was extracted from the resulting distribution via a fit to the background-subtracted invariant-mass distribution, as reported in Fig. 3 (right panel) for the interval 0 < p T < 0.5 GeV/c. In the fit function, the signal was modelled with a Gaussian term, while the residual background with second-order polynomial function. The statistical significance of the signal extracted in 0 The (Acc × ε) correction factors of prompt and feed-down D 0 mesons were determined from the same Monte Carlo simulations, as those used for the analyses with decay-vertex reconstruction. The (Acc × ε) obtained with the two different analyses are compared in Fig. 4. For the analysis that does not exploit the selections on the D 0 -meson decay vertex, the efficiency is higher by a factor of about 30 (3) at low (high) p T and is almost independent of p T . The mild increase with the increasing p T is mainly determined by the geometrical acceptance of the detector. Unlike in the analysis with decay-vertex reconstruction, the efficiency is the same for prompt D 0 and for feed-down D 0 , as expected when no selection is made on the displacement of the D 0 -meson decay vertex from the interaction point.
The prompt fraction to the D 0 -meson raw yield, f prompt , was estimated with the same FONLL-based approach used for the analysis with decay-vertex. The resulting f prompt values decrease with increasing p T , from a value of about 0.95 for p T < 4 GeV/c to about 0.90 in the interval 8 < p T < 12 GeV/c and are larger compared to the analysis with decay-vertex reconstruction, due to the fact that the feed-down component is not enhanced by the topological selection criteria.

Measurement of the fraction of prompt D mesons
The fractions of prompt D 0 and D + s mesons in the raw yields, f prompt , were measured exploiting the different shapes for the distributions of the transverse-plane impact parameter to the primary vertex (d 0 ) D-meson production in pp collisions at √ s = 5.02 TeV ALICE Collaboration In this function, S and B are the signal raw yield and background in the selected invariant-mass range, constrained from the values obtained from the invariant-mass fit; are the functions describing the impact-parameter distributions of prompt and feed-down D mesons and background, respectively. The function F prompt is a detector resolution term modelled with a Gaussian and a symmetric exponential term. The function F feed-down is the convolution of a sum of two symmetric exponential functions (F feed-down true ), which describe the intrinsic impact-parameter distribution of secondary D mesons from beauty-hadron decays, and the detector resolution term (F prompt ). All the parameters of the F prompt and F feed-down true functions were fixed in the data fit to the values obtained by fitting the distributions from Monte Carlo simulations, except for the Gaussian width of the detector-resolution term, which was kept free in order to compensate a possible not perfect description of the impactparameter resolution. The distribution describing the combinatorial background was parametrised by fitting the impact-parameter distribution background candidates in the side bands of the signal peak in the invariant-mass distributions, with a function composed by a Gaussian and a symmetric exponential terms (F backgr ). Figure 5 (left) shows examples of fits to the impact-parameter distributions of D 0 and D + s mesons in the transverse-momentum intervals 3 < p T < 4 GeV/c and 5 < p T < 6 GeV/c, respectively. For this study, wider p T intervals were adopted compared to the analysis, due to the poor quality of the fit when reducing the sample. The D 0 candidates used in the impact-parameter fit were selected with the same criteria described in Section 3.1. For the D + s mesons, the impact-parameter selection, used to extract the raw yield from the invariant-mass distribution, was not applied for this study. In this case, the prompt fraction, f prompt , was obtained by integrating the functions obtained from the fit in the restricted impact-parameter range used in the analysis. Entries and charge conj. background); (ii) the uncertainty on the signal and background yields extracted from the invariant-mass fits; and (iii) the consistency of the procedure, evaluated with a Monte Carlo closure test. These uncertainties were estimated with the procedures described in Ref. [19]. The total systematic uncertainty on f prompt with the data-driven approach ranges, depending on p T , between 1% and 9% for the D 0 meson, and between about 4% and 17% for the D + s meson. The prompt fractions in the raw yields of D 0 and D + s mesons measured with the data-driven method are compared to the one calculated with the FONLL-based approach in the right panels of Fig. 5 and found to be compatible within uncertainties. For the interval 24 < p T < 36 GeV/c (16 < p T < 24 GeV/c), given the poor precision of the impact-parameter fit, it was not possible to determine the data-driven prompt fraction for the D 0 (D + s ) meson.
D-meson production in pp collisions at √ s = 5.02 TeV ALICE Collaboration

Systematic uncertainties
Systematic uncertainties on the D-meson cross sections were estimated considering the following sources: (i) extraction of the raw yield from the invariant-mass distributions; (ii) track reconstruction efficiency; (iii) D-meson selection efficiency; (iv) PID efficiency; (v) the shape of the p T spectrum generated for D mesons in the simulation; (vi) subtraction of the feed-down from beauty-hadron decays. In addition, the uncertainties on the branching ratios and on the integrated luminosity were considered. A summary of the systematic uncertainties is reported in Table 1 for different p T intervals.
The systematic uncertainties on the raw yield extraction were evaluated by repeating the fits several hundred times varying the fit interval and the functional form of the background fit function. The same strategy was performed using a bin-counting method, in which the signal yield was obtained by integrating the invariant-mass distribution after subtracting the background, estimated from a fit to the side-bands only. The systematic uncertainty was defined as the RMS of the distribution of the signal yields obtained from all these variations and ranges between 1% and 9% depending on the D-meson species and p T interval. This includes for the D 0 mesons a contribution of about 1% obtained by varying the ratio of the integral of the reflections to the integral of the signal and the shape of the templates used in the invariant-mass fits. For the background estimation of the D 0 -meson analysis without decay-vertex reconstruction with the track-rotation technique, different configurations of the rotation angle were used. In addition, three alternative approaches were tested to estimate the background distribution: like-sign (LS) pairs, event mixing, and side-band fit [19]. The raw yield values obtained subtracting these alternative background distributions were found to be consistent with those from the default configuration of the track-rotation method within the uncertainty estimated by varying the fit conditions and therefore no additional systematic uncertainty was assigned.
The systematic uncertainty on the track reconstruction efficiency has two different contributions. The first one is estimated by varying the track-quality selection criteria and the second one is estimated by comparing the probability to match the tracks from the TPC to the ITS hits in data and simulation (matching efficiency). To obtain the matching efficiency, the abundances of primary and secondary particles in data were estimated via template fits to the track impact-parameter distributions, where the relative abundances in the simulation were weighted to match those in data [27,38]. The estimated uncertainty, a quadratic sum of the two contributions, depends on the D-meson p T and it ranges from 3% to 5% for the two-body decay of D 0 mesons and from 3.5% to 7% for the three-body decays of D + , D * + , and D + s mesons. The systematic uncertainty on the D-meson selection efficiency originates from imperfections in the simulation of the D-meson decay kinematics and topology and of the resolutions and alignments of detectors in the simulation. For the analyses with decay-vertex reconstruction, the systematic uncertainty was estimated by repeating the analysis with different sets of selection criteria, resulting in a significant modification of the efficiencies, raw yield, and background values. The systematic uncertainties are largest at low p T (up to 5%), where the efficiencies are low and vary steeply with p T , because of the tighter geometrical selections. For the D + s meson, for which more stringent selection criteria were used, slightly larger uncertainties were estimated, ranging from 5% at high p T to 8% at low p T . In the case of the D 0 -meson analysis without decay-vertex reconstruction, the stability of the corrected yield was tested against variations of the single-track p T selection and no systematic effect was observed.
To estimate the uncertainty on the PID selection efficiency, the analysis was repeated without PID selection for the three non-strange D-meson species and D + s mesons with p T > 6 GeV/c. The resulting cross sections were found to be compatible with those obtained with the PID selection and therefore no systematic uncertainty was assigned. For D + s mesons with p T < 6 GeV/c and the D 0 -meson analysis without decay-vertex reconstruction, an analysis without applying PID selections could not be performed due to the insufficient statistical significance of the signal. The systematic uncertainty for low-p T D + s mesons  was therefore estimated by comparing the pion and kaon PID selection efficiencies in the data and in the simulation and combining the observed differences using the D + s -meson decay kinematics [31]. A 3% systematic uncertainty was assigned for 4 < p T < 6 GeV/c, and 2.5% for p T < 4 GeV/c. For the D 0 -meson analysis without decay-vertex reconstruction, compatible cross sections were obtained when using more stringent PID criteria. Based on this result and on the fact that the PID selections are the same as used in the analysis with decay-vertex reconstruction, no uncertainty due to PID was assigned.
The systematic uncertainty due to the generated D-meson p T shape was estimated by using FONLL as an alternative generator with respect to PYTHIA to simulate the D-meson p T distribution [15], and was found to be 0%-5% for p T < 3 GeV/c and negligible at higher p T . The p T shape of both considered distributions were found to be compatible with the measured one within uncertainties. Finally, the systematic uncertainty on the subtraction of feed-down from beauty-hadron decays (i.e. the calculation of the f prompt fraction) was estimated by varying the FONLL parameters (b-quark mass, factorisation, and renormalisation scales) as prescribed in Ref. [8]. It ranges between +1.0 −1.2 % and +4.4 −6.3 % depending on the D-meson species and p T interval.
The total systematic uncertainties were obtained as the square-root of the quadratic sum of all the contributions described previously. The systematic uncertainty on PID, tracking, and selection efficiencies are mainly correlated among the different p T intervals, while the raw-yield extraction uncertainty is mostly uncorrelated. The total systematic uncertainty varies from 6.5%-10.0%, 6.5%-10.5%, 5.4%-11.3%, and 8.7%-12.1% for the D 0 , D + , D * + , and D + s mesons, respectively. The p T -differential cross sections have an additional systematic uncertainty due to the uncertainties on the integrated luminosity [32] and on the branching ratios of the considered D-meson decays [30].

Transverse momentum-differential cross sections
The p T -differential production cross section for prompt D 0 mesons in |y| < 0.5 in pp collisions at √ s = 5.02 TeV was obtained from the analyses with and without decay-vertex reconstruction. The two results are compared in Fig. 6 with the inset showing their ratio in the common p T range. In all the figures in this section, the vertical error bars represent the statistical uncertainties and the systematic uncertainties are depicted as boxes around the data points. In each p T interval the symbols are positioned horizontally at the center of the bin and the horizontal bars represents the width of the p T interval. The two results for prompt D 0 -meson cross section are found to be consistent within statistical uncertainties, which are independent between the two measurements because of their very different signal-to-background ratios and efficiencies. The most precise measurement of the prompt D 0 -meson production cross section is obtained using the results of the analysis without decay-vertex reconstruction in the interval 0 < p T < D-meson production in pp collisions at √ s = 5.02 TeV ALICE Collaboration 1 GeV/c and those of the analysis with decay-vertex reconstruction for p T > 1 GeV/c.
The p T -differential cross sections for prompt D 0 , D + , D * + , and D + s -meson production in |y| < 0.5 are depicted in Fig. 7. The prompt D 0 -meson p T -differential cross section is compatible with the one measured by the CMS collaboration at the same centre-of-mass energy in |y| < 1 and 2 < p T < 100 GeV/c [20].
In Figs. 8, 9, 10 and 11 the measured prompt D 0 , D + , D * + , D + s -meson p T -differential cross sections are compared with results of perturbative QCD calculations performed with different schemes: FONLL [7,8] (not available for the D + s meson), two calculations using the GM-VNFS framework with different prescriptions to regulate the divergences at small transverse momentum, dubbed as GM-VFNS(modµ R,F ) [39,40] and GM-VFNS(SACOT-m T ) [6], and a calculation based on k T -factorisation [41]. In particular, the GM-VFNS(mod-µ R,F ) calculations were performed with a different choice of the factorisation and renormalisation scales µ F and µ R with respect to the GM-VFNS predictions of Ref. [5] that were compared in Ref.
[27] to the cross sections measured at √ s = 7 TeV. With this modification of QCD scale, the calculations could be extended to lower p T . In GM-VFNS(SACOT-m T ), the divergences of the heavy-quark PDFs and light-parton fragmentation functions at low p T are regulated by the heavy-quark mass, thus allowing the calculation of the D-meson cross section down to p T = 0. Note also that the authors of the k T -factorisation calculations changed the treatment of the running strong coupling constant α S and the gluon distributions [41], with respect to the predictions shown in Ref. [27]. In GM-VFNS(mod-µ R,F ) the value of charm mass is set to 1.3 GeV/c 2 , while in FONLL, GM-VFNS(SACOT-m T ) and k T -factorisation predictions the mass is set to 1.5 GeV/c 2 . The four frameworks utilise different sets of PDFs (CTEQ6.6 [42], CTEQ14 [43], NNPDF3.1 [44] and MMHT2014 [45] for The analysis without decay-vertex reconstruction provides also a direct measurement of the inclusive D 0 -meson cross section because no selections are applied on the decay topology, which alter the fraction of prompt and feed-down D mesons. The inclusive D 0 -meson cross section is shown in Fig. 12 and compared with results from FONLL calculations [7,8] with the B → D + X decay kinematics from the EvtGen package [37]. The contributions of prompt D 0 -meson poduction from FONLL and D 0 mesons from B-meson decays from FONLL+EvtGen are also shown separately. The measured cross sections are described by the calculation within the theoretical uncertainties, with the central value of the prediction lying below the data in all the p T intervals, similarly to what observed for prompt D mesons. The mean p T of prompt D 0 mesons, p T , was evaluated for p T >0 with a fit of the prompt D 0 -meson cross section, that is measured down to p T = 0, using a power-law function, as was done in Ref.  Figure 8: p T -differential production cross sections for prompt D 0 meson compared to pQCD calculations: FONLL [7,8], GM-VFNS(mod-µ R,F ) [39,40], GM-VFNS(SACOT-m T ) [6], and k T -factorisation [41]. The ratios of the data to the theoretical predictions are shown in the lower part of each panel.
The systematic uncertainty on the p T was estimated as described in Refs. [19,27]. The contributions due to the correlated and uncorrelated systematic uncertainties on the measured p T -differential cross section were taken into account separately and the contribution due to the choice of the fit function has D-meson production in pp collisions at √ s = 5.02 TeV ALICE Collaboration  Figure 9: p T -differential production cross sections for prompt D + meson compared to pQCD calculations: FONLL [7,8], GM-VFNS(mod-µ R,F ) [39,40], GM-VFNS(SACOT-m T ) [6], and k T -factorisation [41]. The ratios of the data to the theoretical predictions are shown in the lower part of each panel.
been estimated by comparing results obtained using different functions and using a method based on direct calculations of p T from the data points.

D-meson cross-section ratios
The ratios of the p T -differential cross sections of prompt D 0 , D + , D * + , and D + s mesons in pp collisions at √ s = 5.02 TeV are reported in Fig. 13. In the evaluation of the systematic uncertainties on these ratios, the sources of correlated and uncorrelated systematic effects were treated separately. In particular, the contributions of the yield extraction and cut efficiency were considered as uncorrelated, while those of the feed-down from beauty-hadron decays and the tracking efficiency were treated as fully correlated among D-meson production in pp collisions at √ s = 5.02 TeV ALICE Collaboration  Figure 10: p T -differential production cross sections for prompt D * + meson compared to pQCD calculations: FONLL [7,8], GM-VFNS(mod-µ R,F ) [39,40], GM-VFNS(SACOT-m T ) [6], and k T -factorisation [41]. The ratios of the data to the theoretical predictions are shown in the lower part of each panel.
the different D-meson species. The measured D-meson cross-section ratios do not show a significant p T dependence within the experimental uncertainties, thus suggesting no discernible difference between the fragmentation functions of charm quarks to pseudoscalar (D 0 , D + , and D + s ) and vector (D * + ) mesons and to strange and non-strange mesons. The results are compatible within uncertainties with the ratios measured in pp collisions at √ s = 7 TeV [27] 1 .
To study the evolution of prompt D-meson production with the centre-of-mass energy of the collision, 1 The cross section for D 0 and D + mesons in pp collisions at √ s = 7 TeV were updated with respect to Ref.
the ratios of the production cross sections in pp collisions at √ s = 7 TeV [27] and √ s = 5.02 TeV were computed for D 0 , D + , D * + and D + s mesons. The systematic uncertainties on the measured ratios were obtained treating the contribution originating from the subtraction of the feed-down from beauty-hadron decays as correlated, while all the other systematic uncertainties on the cross sections were propagated as uncorrelated between the measurements at the two different energies, except for the uncertainty on the BR, which cancels out in the ratio. The results for D 0 , D + , D * + and D + s are compared in Fig. 14 TeV, from the analysis without decay-vertex reconstruction, compared to FONLL pQCD calculations [7,8] with the B → D + X decay kinematics from the EvtGen package [37] (grey boxes). The contributions of prompt D 0 from FONLL (red) and D 0 from B-meson decays from FONLL+EvtGen (blue) are also shown separately. The vertical error bars and the empty boxes represent the statistical and systematic uncertainties, respectively. the magnitude of the theoretical uncertainties comparable with those of the data.
The rapidity dependence of D 0 -meson production in pp collisions at √ s = 5 TeV can be studied from the ratios between our measurements at midrapidity and the LHCb results in different y intervals at forward rapidity [22]. The precise measurement of the D 0 -meson cross section down to p T = 0 presented in this paper, when analysed together with other results at different centre-of-mass energies and rapidities, can provide sensitivity to the gluon PDF at small values of Bjorken-x (10 −4 -10 −5 ) [24]. In Fig. 15 the ratios of the D 0 -meson production cross sections per unit of rapidity measured with ALICE at mid-rapidity (|y| < 0.5) and by the LHCb collaboration in three rapidity intervals at forward rapidity 2 < y < 2.5 (left panel), 3 < y < 3.5 (middle panel), 4 < y < 4.5 (right panel) [22] are shown as a function of p T . The error bars and boxes represent the uncertainty obtained from the propagation of the statistical and systematic uncertainties, respectively, from the p T -differential cross sections. The systematic uncertainties, including the one on the luminosity determination, were treated as uncorrelated between the ALICE and LHCb results, except for the uncertainty on the BR, which cancels out in the ratio. The central values and the uncertainties of the FONLL calculations are evaluated as described in Ref. [27]. The measured ratios are compatible with FONLL calculations, shown as red boxes in Fig. 15. Nevertheless the comparison seems to hint at a different slope in data with respect to FONLL, since at low (high) p T the data tend to stay above (below) the FONLL central values, in all rapidity intervals.

Transverse momentum-integrated cross sections and ratios
The visible production cross sections of prompt D mesons were evaluated by integrating the p T -differential cross sections over the narrower p T intervals of the D + , D * + , and D + s -meson measurements, in the measured p T range. The results are reported in Table 2. The systematic uncertainty was evaluated by propa-D-meson production in pp collisions at √ s = 5.02 TeV ALICE Collaboration
The ratios of the p T -integrated yields of the different D-meson species were computed from the cross sections integrated over the common p T range. The systematic uncertainties on the ratios were computed treating the BR, yield extraction and cut efficiency uncertainties as uncorrelated among the different species and the other sources as correlated. The results are reported in Table 3.  given D-meson species was computed using the FONLL central parameters to evaluate the ratio between the total production cross section in |y| < 0.5 and that in the experimentally covered phase space. It was verified that the extrapolation factors computed with FONLL were compatible with those resulting from GM-VFNS calculations. The systematic uncertainty on the extrapolation factor was estimated as proposed in Ref. [8], considering sources due to (i) the CTEQ6.6 PDFs uncertainties [42], (ii) the variation of the charm-quark mass and (iii) the renormalisation and factorisation scales in the FONLL calculation. For D 0 mesons, for which the measurement extends down to p T = 0, the extrapolation factor accounts only for the very small contribution of D mesons with p T > 36 GeV/c and therefore its value is very close to unity with negligible uncertainty. The FONLL predictions are not available for D + s mesons, hence in this case the central value of the extrapolation factor was computed as described in Ref.
[27], combining the prediction based on the p T -differential cross section of charm quarks from FONLL, the fractions f (c → D + s ) and f (c → D * + s ) from ALEPH [46], and the fragmentation functions from Ref. [47], which have one parameter, r, that was set to 0.1 as done in FONLL [48]. An additional contribution to the systematic uncertainty was assigned based on the envelope of the results obtained using the FONLL p Tdifferential cross sections of non-strange D mesons to compute the D + s -meson extrapolation factor. The computed extrapolation factors and the prompt D-meson production cross sections per unit of rapidity dσ /dy in |y| < 0.5, are presented in Table 4.
In Ref.
[27], the cc production cross section per unit of rapidity at mid-rapidity (|y| < 0.5) and the total charm production cross sections in pp collisions at √ s = 7 TeV were reported. They were computed from the prompt D 0 -meson production cross section, which was divided by the fraction of charm quarks hadronising into D 0 mesons, f (c → D 0 ) = 0.542 ± 0.024, derived in Ref. [49] by averaging the measurements in e + e − collisions at LEP. However, recent measurements of the Λ + c baryon production cross section in pp collisions at √ s = 7 TeV and in p-Pb collisions at √ s = 5.02 TeV [50] show a significant enhancement of the Λ + c /D 0 ratio for p T >1 GeV/c as compared to the values measured in e + e − and ep collisions at lower centre-of-mass energies. This suggests that the fragmentation fractions of charm quarks into charmed baryons in pp collisions at LHC energies might differ significantly from the LEP results reported in Ref. [49] and that measurements of charmed-baryon production cross sections in pp collisions at √ s = 5.02 TeV are needed for an accurate calculation of the charm production cross section.

Summary
We have reported the measurement of the inclusive p T -differential production cross sections of prompt D 0 , D + , D * + , and D + s mesons at mid-rapidity (|y| < 0.5) in pp collisions at a centre-of-mass energy of √ s = 5.02 TeV, obtained with the data collected at the end of 2017 with the ALICE detector. The measurement was performed in the transverse-momentum range 0 < p T < 36 GeV/c for D 0 , 1 < p T < 36 GeV/c for D + and D * + , and 2 < p T < 24 GeV/c for D + s mesons. It is measured in finer p T bins with respect to the previous measurements at √ s = 7 TeV [27], providing a more detailed description of the cross-section p T shape. The results were compared and found compatible, within uncertainties, with different pQCD calculations performed with different schemes: FONLL [7,8], two calculations D-meson production in pp collisions at √ s = 5.02 TeV ALICE Collaboration using the GM-VNFS framework with different prescriptions [6,39,40], and a calculation based on k Tfactorisation [41]. The ratios of D 0 -meson production cross sections measured with ALICE and LHCb in different rapidity intervals were compatible with FONLL calculations, showing a hint of different slope in data with respect to theoretical predictions. The ratios of the cross sections of D 0 , D + , and D * + mesons at √ s = 7 TeV [27] and √ s = 5.02 TeV are described consistently by FONLL pQCD calculations. The ratios of the p T -differential cross sections of D 0 , D + , D * + , and D + s mesons were found to be compatible within uncertainties with the D-meson cross-section ratios measured in pp collisions at √ s = 7 TeV [27]. The new measurement will allow for a more accurate determination of the nuclear modification factor R pA in p-Pb collisions and R AA in Pb-Pb collisions at √ s NN = 5.02 TeV, due to the larger statistics available and since it is performed at the same centre-of-mass energy of the other collision systems. D-meson production in pp collisions at √ s = 5.02 TeV ALICE Collaboration     [32] ALICE Collaboration Collaboration, "ALICE 2017 luminosity determination for pp collisions at √ s = 5 TeV",. http://cds.cern.ch/record/2648933.