Combined analysis of charm-quark fragmentation-fraction measurements

A summary of measurements of the fragmentation of charm quarks into a specific hadron is given. Measurements performed in photoproduction and deep inelastic scattering in e±p\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^{\pm }p$$\end{document}, pp and e+e-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-$$\end{document} collisions are compared, using up-to-date branching ratios. Within uncertainties, all measurements agree, supporting the hypothesis that fragmentation is independent of the specific production process. Averages of the fragmentation fractions over all measurements are presented. The average has significantly reduced uncertainties compared to individual measurements.

The fragmentation process is soft and hence can not be calculated with the techniques of perturbative QCD (pQCD). Therefore, these measurements are a necessary ingredient for any QCD prediction of charm-hadron production. In this context, it is important to validate the hypothesis that fragmentation fractions are universal, i.e. independent of the hard production mechanism. Thus, once precisely measured in one experiment, they can be applied in any reaction. Another important check is that the sum of fragmentation fractions of all known weakly decaying charm hadrons is equal to unity, thus checking if all weakly decaying states are known.
To achieve these goals, a comparison of fragmentationfraction measurements obtained in different production regimes is performed using a combination of individual a e-mail: andrii.verbytskyi@mpp.mpg.de measurements. Due to independent data sets and different detector types and constructions, the experimental statistical and systematic uncertainties in most cases can be treated as uncorrelated between measurements. However, a careful treatment of correlated uncertainties due to common usage of branching ratios and theory inputs is essential, as in many measurements these are one of the leading uncertainty sources. In the past several combinations of fragmentationfraction data were performed with fewer inputs: the summary of the charm fragmentation fractions in e + e − at the Z resonance [25], the combination of e + e − measurements [26,27] as well as the combination of e + e − and e ± p measurements [28]. Compared to those, the present analysis extends to a larger set of measurements, in particular the final measurement in PHP by the ZEUS experiment at HERA [15], the pp measurements from LHCb [19,20], ALICE [21][22][23], ATLAS [24] and the + c measurements from the BABAR experiment [8]. It uses the up-to-date branching-ratio values [25,29,30], treats correlations of branching-ratio uncertainties and recent theory predictions with reduced uncertainties [31,32] as input.

Combination procedure
The combination of the measurements used in the present analysis is based on numerical χ 2 minimisation with respect to observables of interest. The numerical minimisation was performed with the MINUIT package [35] and the procedure for calculation of χ 2 itself is outlined below. For a set of m measurements and corresponding expectation values calculated from n parameters, a column-vector of the residuals R(1 × m) is calculated as a difference of a measurement and the corresponding expectation. The covariance matrix, V (m × m), is calculated as where U i stands for an uncorrelated uncertainty of ith residual, C i,k stands for the correlated uncertainty of source k of the ith measurement and the sum runs over all sources of correlated uncertainties. The χ 2 is then calculated as The correlated uncertainties are treated multiplicatively in the construction of the covariance matrix, i.e. the relative uncertainties are used to scale the corresponding expectation values instead of the measurement. This avoids the bias for normalisation uncertainties, such as branching ratio uncertainties, which are the main correlated uncertainties considered in the presented analysis. The statistical and uncorrelated systematic uncertainties are treated additively. Data sets and their systematic uncertainties are assumed to be independent between experiments. In addition, most of the measurements do not contain the information about a potential correlation between cross-section values for different charm hadrons. Therefore, in the following all experimental uncertainties are treated as uncorrelated, unless otherwise stated. Uncertainties on the combined values of the fragmentation fractions are determined using the Hessian method with the criterion χ 2 = 1. 1 The evaluated total uncertainties on the free parameters comprise experimental, branching ratio and uncertainties of the λ parameter.
The combination of all the measurements is obtained imposing the normalisation constraint on the sum of all ground state hadrons by adding an additional "measurement" of S calculated from Eq. (6) with an uncertainty on λ and the corresponding prediction S = 1. In order to keep the main result with the normalisation constraint as model independent as possible, no theory inputs on the charm cross-section are used in such a combination, and the fragmentation fractions are calculated according to Eq. (5). For the same reason, any measurements that require theoretical inputs for conversion into cross-sections or fragmentation fractions and do not have these inputs in the original publications are also excluded from the main combination. However, such data are included in a more constrained combination. In the following, treatment of such measurements will be discussed case-by-case in the relevant sections.
The quantities commonly used as Monte Carlo generator parameters, were calculated from the fit results with the full error propagation and taking into account the correlation between parameters.

Charm-quark fragmentation into hadrons in e + e − collisions
Measurements of charm-hadron-production cross-sections in e + e − collisions hadrons were based on the differential Footnote 1 continued to the total uncertainties on the diagonal. The correlated uncertainties considered are those related to λ in Eq. (5) and branching ratios. The residuals are obtained subtracting from the measurements cross-section expectations calculated from n free parameters in the fit, which could be fragmentation fractions, total charm cross-sections, kinematic factors, etc. The details of this calculation are outlined in the each section. With all these components at hand the χ 2 can be evaluated and iteratively numerically minimised with respect to the free parameters. momentum spectrum dσ (e + e − → H )/dx p . The extrapolation to the total cross-section was made in the original papers using a theoretical fragmentation function (e.g. Bowler [36] or Peterson [37]). 2 As mentioned before, the precise predictions of the total charm-production cross-section in e + e − allow calculation of the fragmentation fractions without constraints on the sum of fractions. This way the used hypothesis about the sum of fragmentation fractions (Eq. (6)) can be verified.

Charm-quark fragmentation fractions from measurements at B-factories
The B-factories provided many results on charm-hadron production around the ϒ resonances, which can be used for the calculation of the charm-quark fragmentation fractions in hadrons (see Table 2). The results of the CLEO [1,2] and ARGUS [3][4][5] experiments are represented as a product of the charm-hadron cross-sections times decay branching ratios, σ (e + e − → H ) · B(H → daughters). The BELLE experiment [7] provided measurements of σ (e + e − → H ). The BABAR experiment [8] provided a measurement of an average number of + c → pK − π + decays per hadronic event where R c = σ (e + e − →cc) σ (e + e − →hadrons) is the average number of charmquark pairs per hadronic event. A prediction of R c is needed to use this measurement as an input, therefore, as discussed in Sect. 2.3, it is used only in the case of σ (e + e − → cc) fixed to a theoretical prediction.
For the calculation of the charm-quark fragmentation fractions a fit procedure is used as described in Sect. 2. The total charm-quark-production cross-section is calculated as described in Appendix A. The fit parameters are the fragmentation fractions. The obtained results are given in the middle column of Table 3. The sum of the charm-quark fragmentation fractions into weakly decaying states calculated according to Eq. (6), S ϒ = 0.9701 ± 0.0284, is consistent with unity. 2 A proper extrapolation procedure requires only the hadrons produced directly in fragmentation to be used in the fits. The hadrons produced in decays of excited charm hadrons should be treated separately. In many cases the limited precision of the measurements makes this requirement hard to follow and the decay part of the meson production is treated together with the fragmentation part. In the cases where the contribution of hadrons from decays is comparable with the contribution of the direct production in fragmentation, e.g. for D 0 and D + , the joint treatment could bias the results.  [6,8]. The numbers are given as in the original publications with uncertainties related to branching ratios omitted. The first uncertainty is statistical and the second is systematic 7 .55 ± 0.20 ± 0.34 [6] 231 ± 2 ± 92 [7] D * s → φπ + γ 5.8 ± 0.7 ± 0.5 [6] The combination is also done according to Eq. (2) and imposing the constraint S ϒ − 1 = 0, to be consistent with the definition used for e ± p and pp data. The fit parameters are the fragmentation fractions and the total charm crosssection. The centre-of-mass energy dependence of the charmquark cross-section is accounted for, according to formulae in Appendix A taking the total charm-quark cross-section at a centre-of-mass energy √ s = 10.5 GeV as a reference. The results are given in Table 3 (right column). In this approach, the precise BABAR measurement of + c production [8] is not included in the combination since it requires usage of the R c theoretical prediction. The latter has an influence on other fragmentation-fraction results.

Charm-quark fragmentation fractions from measurements at LEP
The LEP collider provided many results on the charmhadron production. The most valuable for the studies of fragmentation are results obtained from Z decays. Most of those results are represented in the form of fraction of charm events multiplied by branching ratios Table 4). In addition, ALEPH [11], DELPHI [13] and OPAL [10] provided measurements of f (c → D * + ) from the fits of fragmentation functions (see Table 4).
For the calculation of charm-quark fragmentation fractions, a fit procedure is used, as described in Sect. 2. The theoretically calculated value that is used, [32], is in agreement with the experimental world average 0.1721 ± 0.003 [25]. The fit parameters are the fragmentation fractions. The results are given in the middle column of Table 5. The sum of the charm-quark fragmentation fractions into weakly decaying states calculated according to Eq. (6), S Z = 0.9292 ± 0.0261, differs from unity by 2.7 standard deviations.
The combination is also done using Eq. (2), and imposing the constraint S Z −1 = 0, to be consistent with the definition Table 3 Average of charm-quark fragmentation fractions in hadrons in e + e − collisions around √ s = 10.5 GeV. The quantities S, R u/d , P d V and γ s were recalculated from the fit results taking into account correlation of fit parameters. The value of minimised χ 2 and the number degrees of freedom of the fit n dof are given as well used for e ± p and pp data. The fit parameters are the fragmentation fractions and the The results, given in Table 5 (right column), are in good agreement with Ref. [27].

Charm-quark fragmentation into hadrons in e ± p collisions
The charm-hadron-production cross-sections at HERA were measured in a restricted fiducial phase space. The extraction of the charm-quark fragmentation fractions requires a special treatment, as described in detail in Appendix B. The approach followed in this analysis is similar to the one originally used by the ZEUS collaboration [14].

Charm-quark fragmentation fractions from measurements in DIS
Charm-quark fragmentation fractions in DIS in e ± p collisions are calculated from ZEUS and H1 measurements given in Table 6. For the calculation of charm-quark fragmentation fractions a fit procedure is used as it is described in Sect. 2. The free parameters in the fit are the charm fragmentation fractions and pairs of variables Here, σ (c) i is the total charm crosssection in e ± p , while κ i is the kinematic factor for decays from higher states (see Appendix B). The parameter κ is fixed to one for the lowp T measurements in Ref. [17] since the whole p T kinematic space was covered. The present beauty contributions in Ref.

Charm-quark fragmentation fractions from measurements in PHP
Charm-quark fragmentation fractions in PHP in e ± p collisions were calculated from measurements of ZEUS collaboration and given in Table 8. For the update of the latest ZEUS measurement [15] to the decay branching ratios from Table 1, the measured fragmentation fractions are first trans- Table 7 Average of charm-quark fragmentation fractions in e ± p collisions in DIS. The quantities S, R u/d , P d V and γ s are recalculated from the fit results taking into account correlation of fit parameters. The value of minimised χ 2 and the number degrees of freedom of the fit n dof are given as well formed into total charm-hadron cross-sections according to the formulae in Appendix B and only then used in the calculations. In this procedure, the kinematic factor for decays from higher states, κ, is set to 1, since the total phase space is considered from the fragmentation fraction definition, and the σ (c) value cancels out in the procedure. For the calculation of charm-quark fragmentation fractions, a fit procedure is used as it is described in Sect. 2. The free parameters in the fit are the charm fragmentation fractions and pairs of variables σ (c) i | i=1...2 and κ i | i=1...2 for each set of measurement. Here, σ (c) i is the total charm crosssection in e ± p , while κ i is the kinematic factor for decays from higher states (see Appendix B). The sum of charm fragmentation fractions S ep PHP is constrained to unity. The results of the averaging procedure are given in Table 9. The Table 6 Measurements of charm-hadron-production cross-sections in DIS in e ± p collisions. The numbers are given as in the original publications. The first uncertainty is statistical, the second is systematical and the third one corresponds to the branching ratio  obtained fragmentation fractions are in agreement with those obtained in the original publications [14,15].

Charm-quark fragmentation into hadrons in pp collisions
The ALICE experiment measured fiducial cross-sections of D + s [22] and differential p T cross-sections of D 0 , D + and D * + mesons [21,23] at √ s = 2.76 TeV and √ s = 7 TeV. With an integration of the differential cross-sections of D 0 , D + and D * + from Ref. [21] and D 0 from Ref. [23], a coherent set of measurements in the kinematic range 2 < p T < 12 GeV, |y| < 0.5 has been constructed (see Table 10) for the √ s = 2.76 TeV and √ s = 7 TeV. The LHCb experiment provided measurements of charm-hadron cross-sections at √ s = 7 TeV [19] and at √ s = 13 TeV [20]. The ATLAS experiment recently measured the production cross-sections of D * + , D + and D + s mesons at √ s = 7 TeV [24] in the kinematic range 3.5 < p T < 20 GeV, |η| < 2.1.
The measurements together with the correlation matrix for LHCb √ s = 7 TeV are given in Table 10. For the calculation of charm-quark fragmentation fractions a fit procedure is used as it is described in Sect. 2. The free parameters in the fit are: the charm fragmentation fractions, the fiducial cross-sections for LHCb, ALICE and ATLAS measurements in kinematic regions given in Table 10 and corresponding κ parameters.
The constraint S pp − 1 = 0 is imposed. As the Refs. [21,23] do not provide detailed decomposition of the systematic uncertainties, for every bin all systematic uncertainties were conservatively assumed to be fully correlated. For all of these measurements we assume the statistical and systematic uncertainties uncorrelated and luminosity uncertainties -fully correlated within a set of measurements at a given value of √ s. A set of orthogonal fully correlated uncertainties was obtained from the covariance matrix of the √ s = 7 TeV Table 9 Average of charm-quark fragmentation fractions in hadrons in e ± p collisions in photoproduction. The quantities S, R u/d , P d V and γ s are recalculated from the fit results taking into account correlation of fit parameters. The value of minimised χ 2 and the number degrees of freedom of the fit n dof are given as well LHCb measurements with an eigenvector decomposition. The obtained uncertainties are later treated in the same way as other correlated sources in the combination. The Ref. [20] does not contain the correlation matrix for the measurements at √ s = 13 TeV, therefore simplified correlations between measurements were calculated as follows. All of the measurements include 3.9 % fully correlated uncertainty related to luminosity included in the Ref. [20] to the systematic uncertainty. The systematic uncertainties also include the uncertainties on the branching ratios, which were treated correlated with other branching-ratio uncertainties. The remaining systematic uncertainty were treated as fully uncorrelated for different measurements with the same p T cuts and fully correlated for the same measurements with different p T cuts. The statistical uncertainties of σ (D 0 ) p T <8 GeV , σ (D + ) p T <8 GeV were split in two parts, which correspond to p T < 1 GeV and p T > 1 GeV regions. The later were considered fully correlated to the statistical uncertainties Table 10 Measurements of charm-hadron-production cross-sections in pp collisions. For the measurements from Ref. [19] the total uncertainty is given. For the rest of measurements the first uncertainty is statistical, the second is systematic. p T is given in GeV

Decay
LHCb [19] L H C b [ 20] A L I C E [ Table 11 Average of charm-quark fragmentation fractions in pp collisions. The quantities S, R u/d , P d V and γ s are recalculated from the fit results taking into account correlation of fit parameters. The value of minimised χ 2 and the number degrees of freedom of the fit n dof are given as well 0.0639 ± 0.0122 χ 2 6.9 n dof 7 S 1.0000 ± 0.0005 For the LHCb and ATLAS measurements of the D + s → K + K − π + in the limited mass windows M(K + K − ) the following approach was used. The branching ratios were obtained from the integrals over the M(K + K − ) line shape that was parametrised as in Ref. [34] with the total D + s → K + K − π + signal normalised to B(D + s → K + K − π + ) = 5.45 ± 0.17 % [25]. The results are given in Table 1.
The results of the fit are reported in Table 11. In addition to the values of the fragmentation fractions, the fit delivers the inclusive charm-production cross sections in the corresponding fiducial regions, which have particular interest. Therefore, the values of these cross-sections obtained in the global combination with better precision are discussed below.

Selection of measurements for the extraction of fragmentation fractions
The selection of the measurements for the extraction of fragmentation fractions was done according a set of criteria explained below. First, the selection is limited to the measurements obtained in the collisions of high energy particle beams as it assures an absence of possible matter effects and the charm quark production mechanism in these environments is well understood. The measurements of charm-hadron production in proton-meson, proton-nucleon and nucleon-nucleon collisions [38][39][40][41][42][43] were omitted as those provide results in very specific production environment and energy ranges which cannot be easily compared to the results in other experiments.
The third criterion of the selection is the availability of sufficient measurements in the given physical environment needed for the extraction procedure. Several results on charm production in e ± p collisions (e.g. Ref. [56]) and pp collisions (e.g. Refs. [42,43,57]) do not contain enough simultaneous measurements of hadron production and, therefore, cannot be treated independently and/or constrain the results of the combination.

The global combination
To check the consistency of the data from different production regimes and also to extract the charm-quark fragmentation fractions with high precision, all input measurements introduced in the previous sections are used together to produce a global combination. As discussed in Sect. 3.1, the + c measurement by the BABAR experiment [8] is not included while obtaining the combined result. The free parameters of the fit are the charm-quark fragmentation fractions and pairs of variables σ (c) i | i=1...5 and κ i | i=1...5 for three DIS and two PHP sets of measurements, cc had , σ (e + e − → c) at √ s = 10.5 GeV, and the fiducial charm-quark crosssection and κ parameters in pp collisions, at √ s = 7 TeV and √ s = 13 TeV, corresponding to the phase space of the measurements. The constraint on the sum of the cross-sections of the weakly decaying charm states, S, is imposed in the combination, i.e. the prediction for the total charm cross-sections in e + e − collisions is not used, in order to minimise model dependence of the averaging procedure. The result of averaging e + e − , e ± p and pp data, with the constraint S = 1 is presented in the middle column of Table 12 and is shown in Fig. 1. The correlations between the fitted parameters are given in Table 13. The input data are in very good agreement with χ 2 /n dof = 65.6/64. The result of the combination has significantly reduced uncertainties compared to individual measurements.
As an alternative, the combination is also performed using both the constraint on S as well as theoretical predictions of charm production in e + e − collisions and Z decays, i.e. also allows to include the precise BABAR measurement of + c production [8] using the R c calculation as described in Appendix A, which significantly affects the averaged value of f ( + c ). The result of the averaging procedure with this approach is given in the right column of Table 12 for completeness. The result is more model dependent than the default combination, but has a higher precision. At the same time, the result visibly differs from the result of the default procedure. This may partially be traced to the value S Z = 0.9292 ± 0.0261 for the accurate LEP measurements,  Fig. 2 The values of R u/d , P d V and γ s in different experiments with the S constraint. The global combination with the S constraint is shown with the shaded band. Combinations of included data in different production regimes are shown with various full symbols. Data that were not included in the combination [57,58] are shown with open symbols. Note, that the latter are quoted from the original papers, i.e. without correction to the up-to-date branching ratios and with no branching ratio uncertainty, if not given in the source which differs markedly from 1 (see Sect. 3.2). This difference is also reflected in the larger χ 2 /n dof value compared to the default global combination. The difference in the f (c → + c ) precision is to a large extent due to inclusion of the precise BABAR data [8].
The extracted R u/d , P d V and γ s factors are provided in Table 12 and shown in Fig. 2. The combined data are also compared to recent measurements [57,58] that were not included in the combination. In particular, R u/d = 1.097 ± 0.035 is in fair agreement with the isospin invariance hypothesis R u/d = 1 within 2.7 standard deviations. The values of the σ ( pp → c) cross-sections, obtained in the global fit (see Table 14) are consistent with those obtained in the original analysis, but have significantly reduced uncertainties. The consistent treatment of the LHCb and ALICE measurements in the combination procedure allows unbiased calculation of the ratio of the inclusive fiducial charm-quark production cross-sections:  [59]. The R 13/7 value is visibly higher than the theoretical prediction R 13/7 (theory) = 1.39 +0.12 −0.29 [60].

Excited states
In addition to the average fragmentation fractions for the ground, L = 0, states, some fragmentation fractions for the excited, L = 1 charm hadrons are calculated. The measurements used for the calculations are shown in Table 15. The unpublished measurement of f (c → D + s1 ) from Ref. [61] was not used. The fragmentation fractions were not updated to the most recent branching ratios, as the difference between the used branching ratios and the newest is negligible in comparison to statistical and systematical uncertainties of the measurements, and is well below the given numerical precision of the individual measurements in Table 15.
The averages are calculated with an assumption of fully uncorrelated statistical and systematical uncertainties. The results of the averaging procedure are given in Table 16. The strangeness-suppression factor for L = 1, J = 1 + charm mesons is calculated neglecting D(2430) 0 contribution and assuming D + 1 is 1 + state: .

Summary
A summary of measurements of the fragmentation of charm quarks into a specific charm hadron is given. The analysis includes data collected in photoproduction and deep inelastic scattering in e ± p collisions and well as e + e − and pp data. Measurements in different production regimes agree within uncertainties, supporting the hypothesis that fragmentation proceeds independent of the specific production process. Averages of the fragmentation fractions are presented. The global average has significantly reduced uncertainties compared to individual measurements. In addition, the hypothesis that the sum of fragmentation fractions of all known weakly decaying charm hadrons is equal to unity is checked to hold within 3 standard deviations using the e + e − data.
σ (e + e − → l + l − ) = 4α 2 (s)π/3s (8) is the total cross-section of massless charged lepton pair production. In this work, the calculations of the r q (s) were done according to Ref. [31] at the reference energy of √ s = 10.5 GeV and assuming the c quark is the heavy one. The constants used for the calculations in Eqs. (7) and (8)  can be compared with the existing measurements and predictions. It was found that R had is in agreement with the direct measurement from CLEO below √ s = 10.56 GeV R had,CLEO = 3.591 ± 0.003 ± 0.067 ± 0.049 [68] and R c is in agreement with the CLEO Monte-Carlo based estimation R c,CLEO = 0.37 ± 0.05 [1]. For all the theoretically calculated values, the uncertainties of calculations are negligible.

Appendix B: Extraction of the charm-quark fragmentation fractions from the measurements in the restricted phase space
Often the production cross-sections of charm hadrons are measured in some restricted (e.g. in transverse momentum and pseudorapidity) kinematical region v and cannot be used directly (i.e. without extrapolation to the full kinematical space) in the Eq. (2). To avoid the extrapolation and obtain unbiased charm-quark fragmentation fractions, the following approach is used.
The full cross-section of a charm hadron, σ (H ) ∈v can be split into cross-section of direct production σ (H ) dir,∈v and the contribution from the decays of heavier charm states H * , where k H * →H < 1 is a fraction of H * → H decays with H in v. The lack of experimental data in this analysis allows to consider only the heavier charm states that are giving the largest contribution to the total cross-sections. For this reason all D 0 and D + produced not in D * + and D * 0 decays are considered as produced directly. Because of similar kinematics of D * → D decays it is also assumed that k D * →D = k. The measurements on heavier charm baryons and charm strange mesons are absent, but all of those decay dominantly to + c and D + s so we treat all + c and D + s as produced directly. Within a common phase space, the impact of the different masses of the hadrons on the fragmentation process is neglected.
With the assumptions above, we have the equations: In the full kinematical space: In general, to solve the system the measurements of D * 0 production are needed. However, these can be avoided with an assumption of isospin invariance: The last two systems are the working equations for the calculation of the charm fragmentation fractions from the crosssection measurements in the restricted phase space.