Cross section of D-mesons production via the c-quark fragmentation process at next-to-leading order accuracy in pp collisions at the 7TeV LHC

In this paper, the differential cross section of charmed mesons production (Ds±\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_{s}^{\pm }$$\end{document}, D0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^{0}$$\end{document}, D∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^{*}$$\end{document} and D∗±\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^{*\pm }$$\end{document} and D±\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^{\pm }$$\end{document}) are computed as a function of transverse momentum pT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_{T}$$\end{document} and the total cross sections. The hadronization mechanism is considered via a fragmentation process of a c heavy quark to a light-heavy meson according to Suzuki’s model. The two-cutoff phase-space slicing method has been applied to calculate the phase-space integral of two- and three-body subprocess cross sections and fragmentation function (FF) at next-to-leading order (NLO) accuracy. By using Monte Carlo (MC) techniques, we try to properly estimate the ratio of neutral to charged D-meson production Ru/d\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{u/d}$$\end{document} and the strangeness-suppression factor γs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma _{s}$$\end{document}. Our numerical results show that the NLO corrections of the FF have significant effect on D-mesons production and enhance the cross section ratio at NLO respect to leading order (LO) about 50%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$50\%$$\end{document}. A comparison indicates that our cross section results are in agreement (within uncertainties) with available data from ATLAS and ALICE and theoretical predictions of the fixed-order next-to-leading-logarithm (FONLL) approach and NLO parton-shower MC simulation (MC@NLO) in high-energy proton–proton (pp) collisions at s=7TeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=7TeV$$\end{document}.


Introduction
The charmed mesons are the ground states of the c,q (which q can be the light quarks of u, d and s). They are abundantly produced in high-energy colliders. The D ± and D 0 pseudoscalar (P) mesons are binding systems of cd and cū with spin 0, respectively, and the D * ± and D * are their corresponding vector (V) states with spin 1. Carrying a heavy flavor with a light quark, they make the unique heavy mesons in the standard model. Thus, the study of them has attracted a e-mail: dadfar@phys.usb.ac.ir (corresponding author) particular attention of experimenters and theoreticians, particularly, ever since observation of the first charm quark in 1976 and subsequent detection of the D mesons [1].
Almost all these data have been compared with different theoretical predictions, including the fixed-order next-toleading-logarithm (FONLL) approach in Refs. [19,20], the general-mass variable-flavor-number scheme (GM-VFNS) [21,22] calculations, and the NLO quantum chromodynamic (NLO QCD) calculations matched with an MC@NLO simulation [23]. The cross section of D mesons based on the k t factorization approach has been published in the LO and NLO accuracies at LHC energy and compared to experimental results [24].
In the highp T region, the main mechanism for producing a hadron is parton fragmentation. The theoretical predictions of the heavy hadron production cross section, in this region, are based on the factorization approach [19,20]. In this way, the cross section of hadron production is computed as a convolution of three terms: the partonic hard scattering cross section (σ ), the parton distribution functions (PDFs) of the incoming protons, and the FF, D L O (z), that give the probability for the splitting of a parton with high transverse momentum into the desired hadron plus other partons. In fact, the factorization theorem allows us to separate the short-distance, high-energy parton production, and the long-distance fragmentation process. Thus, it provides necessary tools to derive a cross section of hadron production in a systematic manner. The short-distance part is computed through a perturbative calculation [25], and the PDFs and the FFs can be parameterized based on experimental data. The PDFs of the proton are prevalently parameterized by different predictions such as MSTW [26] and MMHT [27].
The typical formula of a differential cross section for producing a hadron h in the collision of hadrons of A and B can be factorized at LO in α s as: where the ψ (4) functions have been calculated by Ellis and Sexton in Ref. [28], and A and B are the incoming hadrons. In Eq. (1), p = {p A , p B , p 1 , p 2 } and a = {a A , a B , a 1 , a 2 } denote the four vectors and flavor indices of partons, respectively, s = ( p A + p B ) 2 and w(a) are the number of spins and colors for particle a, x A and x B are momentum fractions, z is the energy fraction of the produced hadron, and Different theoretical models have been proposed for the mechanism of heavy-light meson production via a fragmentation process at LO and NLO accuracies. Each model contains different parameters which can be effective in the theoretical results of the FF. The LO FF of a heavy quark to light-heavy meson was first calculated by Peterson et al. [29]. The function obtained in Peterson's model is proportional to E −2 , where E is the difference between the final state energy and initial heavy quark energy. Collins and Spillers [30] evaluated the FF at a typical value of k 2 T . In their model, the probability of the transition is related to the hadronic structure function. In Bratten's model, the PQCD fragmentation function was calculated from the general gauge-invariant definition by integrating over the invariant mass s of the fragmenting quark. The FF of Suzuki's model [31] is calculated in Feynman gauge by neglecting Fermi motion, integrating over the transverse momentum of final parton and averaging over the spin of the initial heavy quark states and sum over the spin of the final states. In Ref. [32], the FF is calculated at LO accuracy for producing a free on-shell QQ state D Q→h according initial conditions of Suzuki's model by considering different spin states of the initial quark.
The experimental observations indicate that the heavy quark fragmentation into the D mesons is not well modeled in theoretical calculations. In order to achieve a better theoretical estimation of the charmed mesons production cross section, here we consider the FF at the NLO accuracy in the factorization approach as well asσ and the PDFs at NLO accuracy.
The FFs of a heavy quark into a light-heavy meson [33] and heavy quarkonia [34] have been computed at NLO accuracy by adopting the assumptions of Suzuki's model at LO accuracy [31]. The NLO corrections have been applied on the right side of the fragmentation process c( p 2 ) → cq(k) + q(k 3 ) in which q can be light quarks (u, d, s). In this process, is the four momentum of the fragmenting parton, k = z( √ s 1 , 0, √ s 1 , 0) is the four momentum of the final meson and is the four momentum of the outcoming quark. In the real corrections, a gluon with the four momentumk = (x √ s 1 , 0, x √ s 1 , 0) is added on the right side of this process so that it carries momentum fraction x of the initial parton's momentum, and this fraction is zero in the rest of the NLO corrections. The conservation of four energy-momentum yields to the relation p 2 =k + k 3 + k, so in the rest frame of the Qq binding state (the same as Ref. [31]), the minimum squared energy of the fragmenting quark is obtained as , where m q and m Q are the mass of the light quark and the heavy quark, respectively, and M = m q + m Q .
The FF is obtained for this process at the NLO accuracy as: where |M| 2 is the squared amplitude, and d 3 and d 2 are the three-body and two-body phase spaces of the final states in the real and virtual corrections. The phase-space slicing technique with two cutoffs [35], collinear and soft cutoff δ c and δ s , respectively, has been applied for calculating the phase-space integral and the cancellation of soft and collinear singularities. The dimensional regularization is adopted to regularize the infrared (IR) and ultraviolet (UV) divergences, and it is done in d = 4 − 2ε dimensional space-time. The final-state hard collinear and soft singularities appearing in the FF are canceled with the addition of the LO diagrams to the renormalized one-loop real and virtual diagrams. More details of the NLO calculations by using dimensional regularization technique and two slicing method have been completely explained in Refs. [33][34][35][36].
After cancellation of all the divergences, the total FF for this process at the NLO accuracy is obtained as follows: where D L O can be referred to the FF at LO for vector, V , and pseudoscalar, P, states [33,34]. Suzuki's FFs for pseudoscalar and vector state are respectively D L O P and In the above equations, N is a normalization constant, a = m q is the mass invariant for producing q + Qq according to the kinematic of Ref. [33], and μ 0 = m Q + 2m q is the initial scale of factorization. F(z), H (z) and t 12 = −2k 1 · p 2 are the functions of z as: All terms of order and the soft and collinear cutoffs ,δ s and δ c , are neglectable. The final function D N L O (z, μ 0 ) depends on z and a factorization scale μ 0 . Since the FF generally contains terms like ln(μ 2 /m 2 Q ), in order to correctly take into account the possible large logarithmic terms, it should be evolved from μ 0 to a higher and arbitrary scale μ f by using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations [37,38].
Here, we would like to derive the cross section at NLO accuracy for pseudoscalar and vector D mesons production in the 1S wave state at a center of mass energy √ s = 7T eV via the fragmentation process of a heavy c quark. By doing this, we indicate the importance of the NLO corrections of FF in the production cross section of these light-heavy mesons.
In the remainder of the present paper, in Sect. 2, we briefly present the calculation of the production cross section using the factorization approach and two slicing method. In Sect. 3, the numerical evaluation of the cross section is performed for vector and pseudoscalar states of D mesons produced via the fragmentation process of the c quark at NLO accuracy. In this section, our numerical results are presented and compared with experimental data and theoretical predictions in the lowand high-transverse-momentum p T region.

Inclusive production cross section at NLO
The calculations of the production cross section of hadrons are typically performed with a combination of analytical and numerical integration techniques. In the course of the NLO calculations, UV singularities appear in loop integrals where the momenta go to infinity. They are removed through the renormalization tools. Soft or IR divergences often arise in the calculations owing to the existence of a massless field, such as a gluon in QCD theory. They appear in both loop and phase-space integrals and are found in the low-energy region where the integration momenta go to zero. The calculation of these singularities should be made in such a way that no information of observable physics quantities is lost in terms. There are several techniques for handling the cancellations of divergences in the higher-order corrections, such as MC techniques. In this technique, the phase-space slicing method can be served to perform NLL calculations. In this method, two cutoff parameters are used to separate the phase-space regions of the soft and collinear singularities from the nonsingular regions. This method would be useful for the transition from partonic to hadronic states and also study of a wide variety of high-energy processes. On this basis, we use this method as tools for calculation of the production cross section of D mesons at pp collision. For the case of LO, the differential cross section is calculated through Eq. (1). This equation involves lowest order two-to-two subprocesses, dσ 2→2 , which have two-body final states. In addition, the FF obtained at the LO is served to calculate the dσ L O . For the calculation of the NLO case through factorization theorems, we should consider the following procedure: as a first step, dσ should be calculated at the NLO that is independent of the species of the produced hadron so it can be applied in the pQCD factorization formula to product a heavy quark (Q) or heavy antiquark (Q). These calculations involve higher-order twoto-three subprocesses, dσ 2→3 , which lead to both two-and three-body final states. The calculation of the squared matrix elements of subprocesses dσ 2→2 and dσ 2→3 is performed based on the VEGAS algorithms described in Refs. [39,40]. The dominant parton processes leading to the production of a heavy quark-antiquark pair in final states are: The subprocess cross section for n body processes at NLO accuracy is written as: where dσ N L O 2→n contains virtual and real corrections: and dσ Real involves the soft and collinear corrections of cross section as: Here, dσ Real and dσ V ir can be precisely found by adapting two-and three-body differential cross sections of Refs. [35,36]. The second step is to apply the NLO FF calculated in Ref. [33] to the total differential cross section. Thus, at NLO accuracy, the general structure of calculation of the total differential cross section is written as: where dσ L O is the LO cross section, dσ N L O A+B→i+h and dσ N L O A+B→i+ j+h denote the NLO cross sections of two-and three-body cases, respectively. Our strategy in calculation of the NLO cross section is regularization of the poles appearing in the amplitude by employing the two cutoff phase-space slicing method.
In the two-body process, i and j are the outcoming partons where the j parton fragments into h hadron. The differential cross section for the two-body process is factorized to the FF as follows: where the contribution of NLO corrections is written as: where G H O and G CT are the higher-order part and factorization counterterm part, respectively, which are a function of ψ 4 and ψ 6 , and their structures are found in Ref. [36]. To remove singularities appearing in dσ N L O A+B→i+h , we introduce the mass factorization scale for fragmentation function as: where μ f is the final-state factorization scale, and P + j j 's are the regulated splitting functions which can be found in Ref. [37]. By replacing Eq. (16) with D(z) in Eq. (14), two body processes are obtained as follows: The two-to-three contribution to the inclusive cross section for producing a hadron h via fragmentation process can be written as: where The expressions for the three-body squared matrix elements may be found in Ref. [28]. The above equation can be split into three terms as in which the three-body phase space,( p A + p B → p i + p j + p k ), is divided into three parts, soft, collinear and noncollinear, with two parameters δ c and δ s . The phase-space integration of the non-collinear region, which is free of singularities, is performed numerically, while in the soft and collinear regions, the phase-space integrations are done analytically in D = 4 − 2 dimensions. The soft terms are as follows: where φ is the usual flux factor which depends on the incident particle masses m A , m B and the partonic center of momentum energy squared, s. Here, the k parton is considered as the third particle. It should be noted that the condition p 3 being soft is defined by δ s > where d 2 is the two-body phase-space differential of i and j partons, and p 3 is the gloun's four momentum in the p A + p B rest frame as: and dp 3−2 3 = dp 0 3 P 0 where Thus, we have: So the soft contribution is written as: where dσ 0 n,m is as follows: where ψ 0 nm is the square the color-connected Born amplitude, and φ is the usual flux factor which depends on the partonic energy squared s and the incident particle masses m A and m B . Integration over their soft region is defined by p 0 3 < δ s √ 2 p A · p B /2. Therefore, the contribution of the soft cross section is as follows: where A s i 's (with i = 0, 1, 2) are in accordance with Eqs. (101-103) in Ref. [35].
The collinear singularities associated with the nonfragmenting parton in the final state correspond to gluon or quark emission. For the collinear region, the integral limit is k 0 3 > δ s √ 2 p A · p B /2. The collinear contribution of the differential cross section at NLO is written as follows: where where s jk 2 p j · p k and P i j are unregulated splitting functions related to usual Altarelli-Parisi splitting kernels, which are as follows Here, the j parton splits into a j and k collinear pair. Now, by inserting Eqs. (32,33) into Eq. (31) and integrating over ds jk , we have: The limits of integration on y for the splitting functions P qq and P gg are δ s ≤ y ≤ 1 − δ s , and for the splitting functions P gq and P qg , are 0 ≤ y ≤ 1 − δ s . It should be noted that a method for evaluating the angular integral of the soft and collinear differential cross section has been mentioned in Ref. [35]. The structure of the finite three-body differential cross section dσ f inite is presented as: where dρ involves the momenta differential in four spacetime dimensions and some kinematic factors. After subtracting the soft and collinear terms from the first term, the singularities are canceled, and the finite contribution is defined as: where μ f is the factorization scale, andf andD functions are modified PDFs and FFs which are given as follows: To obtainD, f should be replaced by D in the above equation. The integration on phase space contains the logarithms in the soft cutoff δ s and single poles in d = 4 − 2ε. The finalstate differential cross section is dependent on the soft and collinear cutoffs through the boundaries imposed on their phase space (δ c δ s ). In the approximation where the cutoffs are small, terms of order δ s can be neglected. The mass singularities are factorized into the FF, and the remaining collinear singularities are factorized and absorbed into the PDFs f a i (x i ) [41,42], so the final result is finite in a d = 4 dimension. The sum of the soft, collinear and ultraviolet renormalized terms is finite, and the final result for physics observables are independent of the cutoffs. In the course of calculations, the average over initial quark spin and the sum over final spin states were performed. In order to make the theoretical predictions for differential and total cross section, the Mathematica FeynCalc package [43,44] was used to carry out the color and Dirac traces, and the Apart package [45] and the FIRE package [46] were used to do partial fractions, and the master integrals were calculated using the LoopTools package [47]. There are Fortran codes corresponding to [28,48] which calculate the one-loop scalar integrals and the matrix elements of amplitude. Here, we attempt to make more reliable predictions for the cross section production of D mesons by using the NLO corrections on the FF and compare them with the presently available results.

The numerical results and conclusion
In this paper, we have computed the production cross section σ , and the differential cross section dσ/dp T , as a function of p T , of charmed mesons (D 0 , D ± , D pm s and vector mesons D * and D * ± ) via the fragmentation process at NLO accuracy in a pQCD approach. The NLO FF calculations of c → cq + q, which q can be u, d and s, have been completely done in our previews work in Refs. [33,34]. Actually, our present results are complementary of previous work. In the course of calculations, we adopted the two-loop formula of the running strong coupling constant α s as follows: where β 0 = 11 − 2/3n f and β 1 = 102 − 38/3n f are the one-and two-loop coefficients of the β function in QCD, respectively. We obtain n f =3 = 0.387 GeV and n f =4 = 0.337 GeV according to α s (m Z ) = 0.118 (with m Z = 91.187 GeV) [49]. In calculations of the production cross section of the hadron at NLO accuracy, the factorization scale μ f should be larger than μ 0 , so it is defined from the order of p T as μ f = m 2 Q + p 2 T,Q . We set the renormalization scale, μ r , to be μ r = 2m q , which is the minimum value of the invariant mass of the virtual gluon.
In numerical calculations of the cross section, the Vegas program was used for integrating over the phase space [40], and the PDFs were adopted from the MSTW predictions at LO and NLO accuracies [26]. Here, the numerical results of cross sections are presented and compared with the ALICE and ATLAS data at energy scale √ s = 7T eV . The ATLAS experimental data have been reported in the range of transverse momentums in a lowp T region (3.5 < p T < 20 GeV) and highp T region (20 < p T < 100 GeV) [50], and the ALICE experimental data have been presented in a region of 2 < p T < 20 GeV. It must be stated that the range of pseudorapidity for the ATLAS data is |η| < 2.1, and the range of rapidity for the ALICE data is |y| < 0.5.
By using the mentioned parameters, the cross sections of D mesons production were obtained numerically at LO and NLO accuracies. In Fig. 1, the differential cross section distributions of the D ± s mesons in c quark fragmentation are presented at LO and NLO accuracies.
The differential cross sections for D 0 and D * mesons are shown in Fig. 2. The D mesons production threshold cannot be lower than approximately 2 GeV. Thus, in low p T from 0 to 2 GeV, the dσ/d p T threshold increases and reaches a maximum at p T = 2 GeV and then decreases towards p T = 100 GeV [50]. The sharp fall-off of the differential cross section distribution is expected in the framework of our study. This distribution for vector meson D * is slightly Fig. 1 The differential cross sections of mesons D ± s as a function of p T compared to ALICE data Fig. 2 The differential cross sections of D 0 and D * mesons as a function of p T in |y| < 0.5 larger than the one for pseudoscalar meson D 0 at both LO and NLO accuracies, which may be due to their mass and spin alignment.
In Fig. 3, our results for D * ± and D ± are compared with ATLAS data at 3.5 < p T < 100 GeV and |η| < 2.1. As it is seen, by increasing the p T distribution, the cross section at the NLO accuracy gets closer to the data than the calculations in LO. Actually, the compatibility of NLO calculations with the data is seen especially in 2 < p T < 8 GeV for D * ± and D ± mesons. The transverse momentum distribution of the D 0 meson is shown in Fig. 4 at 2 < p T < 20 GeV and rapidity |y| = 0.5 and compared with ALICE data and the KMR model [51] for Peterson parameterization with c = 0.02. According to this figure, a comparison shows cross sections at the NLO accuracy have good agreement with ALICE data [4], and it also corresponds to the upper limit of KMR uncer- Fig. 3 The differential cross sections of D * ± (up) and D ± (down) mesons as a function of p T compared to ATLAS data (symbols) tainties. In the bottom panel of Figs. 3 and 4, the ratios of experimental data to our results at both LO and NLO are displayed as a function of p T . The ratio of our results at NLO to ALICE data for D 0 and D ± s mesons gets near to 1 approximately at 2 < p T < 8 GeV. These ratios indicate that our results are in good agreement with the experimental data in the lowp T region. For mesons D ± and D * ± , the downward trend of this ratio is also observed in the range of 8-40 GeV. As can be seen in Figs. 1, 2, 3 and 4, by applying the NLO corrections into the calculations, the differential cross section increases for the all states.
In Table 1, we present our calculations of D 0 and D * mesons cross sections at both LO and NLO accuracies at |η| < 2.1 and at transverse momentums between 3.5 < p T < 20 GeV and 20 < p T < 100 GeV. In Table 2, the cross section values of D ± s , D ± and D * ± are presented in the lowp T region (3.5 < p T < 20 GeV) and highp T region (20 < p T < 100 GeV). In this table, our calculations at LO and NLO accuracies are compared with experimental data from ATLAS and also with theoretical results from FNOLL Fig. 4 The differential cross section of D 0 compared with ALICE data and the KMR model at 2 < p T < 20 GeV Table 1 The production cross sections of D 0 and D * mesons with |η| < 2.1 at 3.5 < p T < 20 GeV. The uncertainties in our results are just due to the quark's mass  [19,20,52] and MC@NLO [23,53] predictions. Therefore, it can be concluded that our results are comparable to these data and predictions within their uncertainties. The uncertainties in our results are just due to the mass of quarks. At a specific collision energy, the mass of the light and heavy quarks has a significant effect on not only the fragmentation probabilities, but also the value of μ and the PDFs and subsequently the production cross section. Here, the effective mass of the heavy quark (c) and the light quarks (u, d, s) and their uncertainties are considered as follows: It should be pointed out that the effective mass of a heavy quark is set to nonzero in the partonic cross sections, the light quarks are considered massless in the proton, and the nonzero values of the effective mass of light quarks only enter through the initial conditions of the FFs. Our numerical results in Tables 1 and 2 show that the NLO QCD corrections increase the lowest-order cross section by about 50% in a factorization approach. Based on the obtained values in Tables 1 and 2, the ratio of neutral to charged D mesons production at both the LO and NLO accuracies can be calculated as: This result agrees with the measured R u/d value in e − p interactions and the e − e + collider which are approximately unity. It is consistent with isospin invariance, which implies that u and d quarks are produced equally in charm fragmentation [8].
At a certain energy, the production probability inversely depends on quark mass, so a decrease in quark mass will increase the production probability, and vice versa. Since the s-quark is heavier than u-and d-quarks, the D ± s mesons cross section becomes smaller than the other D mesons. According to our results in Table 2, at high p T , the strangenesssuppression factor of 0.24 is obtained, which agrees with measured data [50].
where Br(D * ± → D 0 π + ) = 0.677 ± 0.005 [54] is the branching fraction of the D * ± → D 0 π + decay. In summary, we have derived the NLO cross section for a fragmentation process of a c quark to charmed mesons (D ± s , D 0 , D * , D * ± and D ± ) using a factorization approach. In this paper, the phase-space integrals of two-and three-body subprocess cross sections have been calculated at NLO accuracy Table 2 The lowp T , 3.5 < p T < 20 GeV, and highp T , 20 < p T < 100 GeV, cross sections of D * ± , D ± and D ± s production with |η| < 2. by applying the two-cutoff phase-space slicing method. After applying the NLO subprocess cross section and the NLO FF, we have presented the results of a hadronic production cross section in pp collision at high and low transverse momenta.
Our results indicate a good overall agreement with the data, and their discrepancy may be explained as an indication of uncalculated QCD effects. It should be emphasized again that we have considered only uncertainties due to the quark mass.
It should be noted that we have not followed the widespread habit of estimating theory uncertainties by varying factorization and renormalization scales by factors of 2 up and down. Of course, there is also the option to extend the calculations and present numerical results for these additional sources of uncertainties.
Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: We have used the LHC data for the D mesons cross section in this article.] 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://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .