Inclusive \(B\) -meson production at small \(p_T\) in the general-mass variable-flavor-number scheme

The European Physical Journal C, Mar 2015

We calculate the cross section for the inclusive production of \(B\) mesons in \(pp\) and \(p\bar{p}\) collisions at next-to-leading order in the general-mass variable-flavor-number scheme and show that a suitable choice of factorization scales leads to a smooth transition to the fixed-flavor-number scheme. Our numerical results are in good agreement with data from the Tevatron and LHC experiments at small and at large transverse momenta.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-015-3360-6.pdf

Inclusive \(B\) -meson production at small \(p_T\) in the general-mass variable-flavor-number scheme

Eur. Phys. J. C Inclusive B-meson production at small pT in the general-mass variable-flavor-number scheme B. A. Kniehl 2 G. Kramer 2 I. Schienbein 1 H. Spiesberger 0 3 0 PRISMA Cluster of Excellence, Institut fur Physik, Johannes Gutenberg-Universitat , 55099 Mainz , Germany 1 LPSC, Universite Grenoble-Alpes , CNRS/IN2P3, 53 avenue des Martyrs, 38026 Grenoble , France 2 II. Institut fur Theoretische Physik, Universitat Hamburg , Luruper Chaussee 149, 22761 Hamburg , Germany 3 Department of Physics, Centre for Theoretical and Mathematical Physics, University of Cape Town , Rondebosch 7700 , South Africa We calculate the cross section for the inclusive production of B mesons in pp and p p collisions at nextto-leading order in the general-mass variable-flavor-number scheme and show that a suitable choice of factorization scales leads to a smooth transition to the fixed-flavor-number scheme. Our numerical results are in good agreement with data from the Tevatron and LHC experiments at small and at large transverse momenta. In the last 25 years, there has been much interest in the study of B-meson production in p p and pp collisions at hadron colliders, both experimentally and theoretically. First measurements were done in the late 1980s by the UA1 Collabo- - roaft-imonasast ethneerCgEyRoNf Sp p S collider [1] operating at a center S = 0.63 TeV. Later measurements were performed by the CDF and D0 Collaborations at the Fermilab Tevatron in p p collisions at S = 1.8 TeV [2 5] and S = 1.96 TeV [68]. Recently, the CMS [911], ATLAS [12,13], and LHCb [14,15] Collaborations at the CERN Large Hadron Collider (LHC) published their results for inclusive B+, B0, and Bs meson production in pp collisions at S = 7 TeV. In all these measurements, the differential cross sections d/d pT , where pT is the transverse momentum of the B meson, integrated over specific rapidity (y) regions, or d/dy integrated over pT pT,min were obtained. While the CMS and ATLAS data were obtained at central rapidities, LHCb performed measurements at forward rapidities, 2 < y < 4.5. Actually, only a few measurements were able to explore the small- pT range down to pT,min 0, namely those by CDF [6] and the two by LHCb [14,15]. In all other data, the cut pT,min > 5 GeV was imposed. A unified theoretical description of these data, which covers both the very small and the large pT range, requires special efforts, which we shall undertake in this paper. On the theoretical side, it is generally accepted that for the treatment of B-meson production at small pT values, as well as for the calculation of the integrated cross section including the small- pT range, i.e., with pT of order O(mb) and below, where mb is the b-quark mass, one should use the so-called massive scheme or fixed-flavor-number scheme (FFNS) [16 20], in which the number of active quark flavors in the initial state is limited to n f = 4, and the b quark appears with explicit mass dependence only in the final state. In this case, the b quark is treated as a heavy particle in the final state and not as a parton in the initial (anti)proton. In the large- pT region, characterized by pT mb, the socalled massless scheme or zero-mass variable-flavor-number scheme (ZM-VFNS) [2132] is considered to be appropriate. This is the conventional parton model approach, where the b quark is considered massless like any other parton. The b quark is also treated as an incoming parton coming from the (anti)proton leading to additional contributions from hardscattering subprocesses besides those with u, d, s, and c quarks or the gluon (g) in the initial state. Although this approach can be used as soon as the factorization scales I and F associated with the initial- and final-state singularities are above the starting scale of the parton distribution functions (PDFs) and fragmentation functions (FFs), its predictions are reliable only in the large- pT region, pT mb, where terms of the order of m2b/ pT2 can be safely neglected. A next-to-leading-order (NLO) calculation in this scheme automatically resums leading and next-to-leading logarithmic terms. At the same time, all non-logarithmic terms through O(s ) relative to the Born approximation are retained for mb = 0. With the conventional choice of renormalization and factorization scales, R = I = F = mT with mT = pT2 + m2b, the results are dominated by contributions from the b-quark PDF down to pT 0. On the other hand, the general-mass variable-flavornumber scheme (GM-VFNS) [3348] provides a theoretical framework which combines the FFNS and the ZM-VFNS. The ZM-VFNS is extended into the intermediate- pT range by retaining the mass-dependent terms of the FFNS. With the conventional choice of scales indicated above, the results in this scheme are also dominated by the contributions of the b-quark PDF down to pT 0. As a consequence, there is no smooth transition from the GM-VFNS at large pT values to the FFNS in the small- pT range, and the GM-VFNS fails to describe the small- pT CDF and LHCb data. The dominance of contributions with b quarks in the initial state at small pT values is linked to the fact that this part is treated in the massless scheme, as a calculation of the b-quark-initiated subprocesses in a scheme with massive partons (like the ACOT scheme [49]) is not available for hadroproduction.1 The cross section with massless partons is, however, divergent for pT 0. For a realistic description, we thus have to find a way to eliminate or modify this contribution in the small- pT region. In this paper, we shall develop an approach to modify the GM-VFNS in such a way that it matches the FFNS with the exact mb dependence by a suitable choice of I and F . We shall study how these modifications can lead to a better agreement with presently available experimental data at small pT values. The content of this paper is as follows. In Sect. 2, we introduce our strategy for the transition to the FFNS and compare our predictions with the CDF [6] and LHCb [15] data. We shall also present results to be compared with measurements by the ATLAS Collaboration [13]. Our conclusions are presented in Sect. 3. 2 Small- pT results and comparisons with data In this section, we shall discuss a viable unified framework for theoretical predictions of inclusive B-meson production at small and large pT values. We shall compare with the cross section distributions d/d pT measured by CDF [6] and LHCb [15]. Throughout this paper, we take the b-quark pole mass to be mb = 4.5 GeV, evaluate s(n f )(R ) at NLO with n f = 4 and (M4)S = 328 MeV if R < mb and with n f = 5 and (M5)S = 226 MeV if R > mb, and use the CTEQ6.6M proton PDFs [51] unless otherwise stated. We start with results to be compared with the CDF data [6]. In Fig. 1, we show NLO predictions in the FFNS with n f = 4. The full line shows the result for the default choice of 1 For deep inelastic scattering, heavy-quark-initiated processes at NLO with massive quarks have been considered in Ref. [50]. FFNS S = 1.96 TeV -1.0 y 1.0 Tevatron Data default scale = (pT2+mb2) 1 0 2.5 5 7.5 10 12.5 15 17.5 20 22.5 25 pT (GeV) Fig. 1 NLO predictions for d/d pT of p p B+ + X with S = 1.96 TeV and 1 y 1 in the FFNS are compared with CDF data [6]. The dashed lines represent the theoretical uncertainty estimated by varying R and I up and down by a factor of 2 about the default choice R = I = mT scales, i = i mT with i = 1 for i = R, I , while the dashed lines represent an estimate of the theoretical error obtained in the usual way, by varying i up and down by a factor of 2. We take the transition of b and b quarks to the observed B mesons into account by using the branching fraction B(b B) = 39.8 % [52] as an overall normalization factor. The prediction in the FFNS agrees with the CDF data quite well, within experimental errors, up to pT 15 GeV. Beyond this value of pT , the FFNS starts to overestimate the data, as has been shown already in our previous publication [42]. In the FFNS, there is no need for FFs. However, a F independent FF might be introduced on phenomenological grounds and because of theoretical considerations to guarantee a proper matching between the schemes with n f = 4 and n f = 5. In the left panel of Fig. 2, we show results obtained using the F -independent Peterson FF [53] with parameter = 104. We find only marginal differences with respect to the case where a constant branching fraction is used. Note that there are no g, q, q B transitions in the FFNS. In addition to uncertainties from scale variations, there are also uncertainties due to errors in the input. We postpone the discussion of errors in the parametrizations of the PDFs to when we present predictions for the LHCb experiments, in Fig. 8 below, but instead show the influence of mb variations on the default prediction for the Tevatron measurement in the right panel of Fig. 2. At small pT values, the uncertainty is comparable in size with the scale uncertainty, but it is negligible for pT 2mb. In Ref. [42], we presented detailed comparisons of FFNS and GM-VFNS results with CDF data for B+ production FFNS; =0.0001 i= i mT mb=4.5 GeV 2.5 5 7.5 10 12.5 15 17.5 20 22.5 25 1 pT (GeV) FFNS; No FF; i=mT mb=4.3 GeV mb=5.0 GeV Fig. 2 NLO predictions for d/d pT of p p B+ + X with S = 1.96 TeV and 1 y 1 in the FFNS are compared with CDF data [6]. Left panel evaluation using the Peterson FF with = 104 and varying R and I by a factor of 2 up and down about the default choice R = I = mT . Right panel uncertainties from variations of mb d /dpT (nb/GeV) p p (B++B-) X GM-VFNS FFNS COMP S = 7.0 TeV 2.0 y 4.5 LHCb Data default scale = (pT2+mb2) Fig. 3 NLO predictions for d/d pT of pp B + X with S = 7 TeV and 2 y 4.5 in the FFNS (dashed line) and the GM-VFNS (solid line) are compared with LHCb data [15] (see, e.g., Figs. 7 and 8 in Ref. [42]). For calculations in the GM-VFNS, we use the F -dependent FFs described in Ref. [42]. In this case, there are also small contributions due to transitions from light quarks and gluons to B mesons. With the default choice R = I = F = 1, the GM-VFNS predictions diverge for pT 0, in obvious disagreement with the data. We notice, however, that the FFNS and GM-VFNS predictions approach each other at around pT 20 GeV, i.e. 4 to 5 times mb. In Fig. 3, we show a similar comparison of the FFNS and GM-VFNS predictions for the case of the LHCb data [15]. Although these data correspond to much higher S values and to different y ranges compared with the previous results in Ref. [42], we observe similar qualitative behaviors of the FFNS and GM-VFNS predictions and a transition point at about the same value of pT , namely pT 20 GeV. One should expect that a correct treatment of kinematic constraints due to the finite heavy-hadron mass is important, in particular at small pT values. However, there is no unique prescription to take into account the finite masses of the heavy quark and the heavy hadron at the same time. A prescription for including the heavy-hadron mass in connection with massless quarks based on the light-cone scaling variable was introduced in Ref. [54]. Here, we propose a generalization to the massive-quark case. This amounts to scaling the partonic cross section for the production of the massive quark a, with mass ma , energy Ea , and three-momentum pa , which initiates the formation of the heavy hadron H , with mass MH , energy E H = M H2 + pT2 cosh y, and threemomentum pH = M H2 sinh2 y + pT2 cosh2 y, as M H2 z2ma2 RaH = 1 (E H + pH )2 z2ma2 , where z = (E H + pH )/(Eh + ph ) is the light-cone scaling variable. In contrast to naive expectations, one finds RH < 1 corresponding to a slight enhancement of the cross section. In Fig. 4, we show the result of a calculation where this correction factor is taken into account. The effect is small mainly because it is only the difference between the b-quark and B-meson masses that enters. Therefore, also the additional suppression from tighter phase space limits is numerically not relevant. The poor small- pT behavior of the GM-VFNS calculation cannot be remedied by such a naive treatment of phase space restrictions. The results shown in Figs. 1, 2, 3 and 4 suggest that the predictions of the FFNS and the GM-VFNS are very simdefault scale = (pT2+mb2) Fig. 4 NLO predictions for d/d pT of p p B+ + X with S = 1.96 TeV and 1 y 1 in the GM-VFNS are compared with CDF data [6]. Besides the default prediction for scale parameters i = 1 (solid line), a theoretical-error band encompassed by the predictions for i = 1/2 (upper dashed line) and i = 2 (lower dashed line) is shown. The effect of including quark and hadron mass corrections via the phase space in the default prediction is also shown, for MH = 5.28 GeV (dotted line) ilar in the range 15 GeV pT 20 GeV. In this range, both theories are justified , with the FFNS naturally extending to smaller pT values and the GM-VFNS to larger pT values. Therefore, one could designate a transition point t at pT 20 GeV, say, where one switches from one scheme to the other [55, 56]. For practical purposes, this would require the knowledge of the PDFs in the scheme with n f = 4 active flavors up to I = t [57]. A naive prescription to combine the two schemes could be to use matching functions (t2 pT2 ) and [1 (t2 pT2 )] to multiply the FFNS and GM-VFNS cross sections, respectively. Such a prescription will, however, lead to a discontinuity in the combined cross section at pT = t , a property that is certainly unphysical. One could try to invent different ways to combine the two schemes by introducing some smooth matching function, e.g., G( pT ) = pT2 /( pT2 + c2m2) b with c = 5 [58], and imposing, schematically, the prescription = G( pT )GMVFNS + [1 G( pT )]FFNS. However, such an approach would introduce a new parameter, and it is unclear how theoretical uncertainties related to the choice of this parameter should be estimated. It will be helpful to take a closer look at the structure of the partonic cross sections in the GM-VFNS. They involve subtraction terms needed to separate the large logarithms ln( pT2 /m2) at fixed relative O (s ). For example, the cross b section of the partonic subprocess gg gbb contains terms with large logarithms which can be written, schematically, d /dpT sub (nb/GeV) GM-VFNS S = 1.96 TeV default scale = (pT2+mb2) 0 2.5 5 7.5 10 12.5 15 17.5 20 22.5 25 pT (GeV) Fig. 5 Subtraction term for the gg channel evaluated with zero (solid line) and finBite+(d+asXhewdliitnhe)mb value in the NLO prediction for d/d pT of p p S = 1.96 TeV and 1 y 1 in the GM-VFNS as convolutions f (g b) (gb gb) and (gg gg) d(g b) with a partonic PDF f (g b) and a partonic FF d(g b) (see Ref. [39] for definitions and a precise formulation of the corresponding cross sections). These splitting functions are absorbed in the PDFs and FFs, respectively. The large logarithms subtracted in this way are added back by calculating the cross sections of the gb gb and gg gg subprocesses with a b-quark PDF and a g B FF, respectively. There are similar subtraction terms and corresponding contributions involving PDFs and FFs of light quarks and the gluon in other channels as well. In Fig. 5, we show numerical results for the subtraction terms needed in the gg channel (see Eqs. (45)(53) in Ref. [39]), using a linear scale for better visibility. Their contribution is small, but not negligible, even at large pT values. Specifically, the evaluations using the correct prescription with mb = 0 (dashed line) and the approximate one with mb = 0 (full line) are compared with each other. We observe from this that such mass effects are small and cannot be responsible for the unphysical increase of the total result for d/d pT towards small pT values. It is, of course, unavoidable that the subtracted terms differ from those added back in the PDFs and FFs. The subtractions are obtained at fixed order, O (s ) in our case, while the PDFs and FFs contain the large logarithms resummed to all orders. Formally, the differences are of higher order in s . However, these higher-order terms are folded with cross sections calculated in the ZM-VFNS, and these cross sections are singular for pT 0. Therefore, it is not surprising that the contributions with b quarks in the initial state dominate at small pT values. Fig. 6 NLO predictions for d/d pT of p p B+ + X with S = 1.96 TeV and 1 y 1 in the GM-VFNS are compared with CDF data [6]. Left panel evaluations with R = 1 and I = F = 0.5 (solid line), 0.4 (lower dashed line), and 0.6 (upper dashed line). Right panel evaluations with I = F = 0.5 and R = 1 (solid line), 0.5 (upper dashed line), and 2 (lower dashed line). If I,F < 1, then I,F = mb is GM-VFNS S = 1.96 TeV -1.0 y 1.0 Tevatron Data default scale = (pT2+mb2) I and F frozen below =mb default scale = (pT2+mb2) I and F frozen below =mb Obviously, NLO cross sections of b-quark-initiated subprocesses that are convoluted with PDFs and FFs evolved at NLO contain terms singular for pT 0 at one order beyond the subtracted terms. In fact, it would be a major task to derive the missing next-to-next-to-leading-order (NNLO) subtraction terms. With such an extended version of the GM-VFNS including NNLO subtractions, but still using partonic cross sections derived in the ZM-VFNS, the problem would be shifted to one order higher, but it remains to be seen whether the required cancellations of singular terms can be obtained with the required numerical precision. In addition, strictly speaking, the NNLO subtraction terms only make sense in combination with the fixed-order calculation at NNLO. Otherwise, at large pT values, the NNLO subtraction terms and their NLO fixed-order counterparts do not cancel, and the ZM-VFNS is not recovered. In turn, one could argue that b-quark-initiated processes evaluated with b-quark PDFs and FFs should be included in the GM-VFNS at LO only. The predictions thus obtained [59] exhibit better agreement with data in the medium- pT range, between 2 and 7 GeV or so, but switching off these NLO terms is again not sufficient to completely eliminate a singular behavior for pT 0. Moreover, and most importantly, numerical evaluations show that NLO corrections in the zero-mass part of the GM-VFNS are essential at large pT values. Therefore, we do not follow this option either. Instead, we try to exploit the freedom offered by the presence of R , I , and F , parameters that are present anyway. Their values are not determined by theory, but some choice has to be made, based on some reasonable but ad-hoc physical argument. In fact, a judicious choice of scales can lead to a suppression of the potentially dangerous contributions from initial-state b quarks. This exploits the fact that all commonly available PDF fits assume that the b-quark PDF is zero below some starting scale, usually chosen to be I = mb. The same is true for the FFs: the FF for the b B transition vanishes for F < mb. Therefore, with I,F = I,F m T , a value I,F < 1 will render the b-quark PDF and FF zero for pT < mb 1/I2,F 1. In Fig. 6, we show results obtained with such scale choices and compare them with CDF data [6]. Indeed, values for I and F of about 1/2 lead to the required suppression of bquark-initiated contributions. Specifically, in the left panel of Fig. 6, we choose R = 1 and I = F = 0.4, 0.5, 0.6. In all these cases, there is a turn-over near pT = 2.5 GeV, and the agreement with the CDF data is reasonably good. We impose the additional constraint that I and F are not allowed to take values I , F < mb, i.e. the DGLAP evolutions of the PDFs and FFs are frozen below this scale.2 This explains the slight bumps that occur in the pT distribution at pT = 6.0, 7.8, and 10.3 GeV for I = F = 0.6, 0.5, and 0.4, respectively. Obviously, the freedom in the choice of the default values of I and F leads to a moderate extra uncertainty close to the transition region, but the results in the small- pT range, pT 5 GeV, and for large pT values stay unaffected. Taking now R = 1, I = F = 0.5 as the new default scale choice, we estimate the theoretical errors in the usual way by varying the scale parameters by a factor 2 up and down about the default scale choice. In fact, it turns out that the variation of R is the dominant source of the theoretical uncertainties, and we simplify the subsequent calculations by only considering variations of R in the range from 0.5 2 In fact, we have to freeze I slightly below mb, I = Cmb with C < 1, since the b-quark PDF parametrization of the CTEQ6.6M set, which we use here, vanishes only strictly below threshold. For our numerical evaluations, we choose C = 0.99. Fig. 7 NLO predictions for d/d pT of pp B + X with S = 7 TeV and 2 y 4.5 in the FFNS (left panel) and the GM-VFNS (right panel) are compared with LHCb data [15]. The default scale choice is R = 1 and I = F = 0.5 (solid lines), and the theoretical-error bands are obtained by varying R by a factor of 2 up (lower dashed lines) and down (upper dashed lines). If I,F < 1, then I,F = mb is put for d /dpT (nb/GeV) Fig. 8 NLO predictions for d/d pT of pp B + X with S = 7 TeV and 2 y 4.5 in the GM-VFNS with default scale choice are compared with LHCb data [15]. The PDF error is estimated by using the CT10 (solid line) [60], MSTW (dashed line) [61], HERAPDF1.5(NLO) (lower dotted line) [62], and NNPDF 2.3 (upper dotted line) [63,64] sets to 2. Note that we do not introduce an extra prescription to freeze R below mb because, first, the choice of R is not related to switching off b-quark-initiated subprocesses and, second, full variations of R are needed to obtain realistic estimates of the theoretical uncertainty. The resulting error band is shown in the right panel of Fig. 6. We emphasize that the freedom in choosing specific values of I and F as default does not introduce a large additional uncertainty, as may be understood by comparing the left and right panels of Fig. 6. In Fig. 7, we compare the LHCb data [15] with our theoretical predictions implemented with the scale conventions introduced above. In the left panel, we see again nice agreedefault scale = (pT2+mb2) ment between the FFNS predictions and the experimental data. The right panel tells us that the data are not quite as well described by the GM-VFNS predictions, but the agreement is quite acceptable, as the data points are covered by the error band for pT 2.5 GeV. In this pT range, the agreement of the data with predictions obtained in the FONLL scheme [32] is quite similar [15]. In Fig. 8, we show a comparison of results obtained using different PDF parametrizations. Specifically, we use the CT10 [60], MSTW [61], HERA-PDF1.5(NLO) [62], and NNPDF 2.3 [63, 64] PDF sets. The CTEQ6.6M [51] PDFs adopted elsewhere in this paper yield a result that is very similar to the one obtained using the CT10 PDFs. We observe that there are differences in the small- pT range, pT 4 GeV, which are somewhat larger than the experimental errors. We should, therefore, expect that B-meson production data at the LHC will help us to further constrain the PDFs. In particular, there is sensitivity to the gluon PDF in this kinematic range. For completeness, we also consider the production of B0 and Bs0 mesons. Appropriate experimental data were published by the LHCb Collaboration in Ref. [15]. In Fig. 9, we present comparisons with NLO predictions in the GMVFNS using the scale setting and theoretical-error estimation prescriptions described above to find good agreement. Finally, in Fig. 10, we compare B+-meson production data taken by the ATLAS Collaboration [13] with our NLO GMVFNS predictions. These data extend into the very-large- pT range, where we expect the GM-VFNS to be quite appropriate. Indeed, we find good agreement, except for the lowest pT bin, 913 GeV, and for central rapidities, where the data are slightly overestimated. Because of the large pT values probed, the adjustment of scales to match to the FFNS as described above is not an issue here. Fig. 9 NLO predictions for d/d pT of pp B0/B0 + X (left panel) and 0 wp pithBs0/Bs + X (right panel) S = 7 TeV and 2 y 4.5 in the GM-VFNS are compared with LHCb data [15]. The default predictions and theoretical-error bands are evaluated as in the right panel of Fig. 6 Fig. 10 NLO predictions for d/d pT of pp B+ + X with S = 7 TeV and 0 y 0.5 (upper left panel), 0.5 y 1 (upper right panel), 1 y 1.5 (lower left panel), and 1.5 y 2.25 (lower right panel) in the GM-VFNS are compared with ATLAS data [13]. The default predictions and theoretical-error bands are evaluated as in the right panel of Fig. 6 d /dpT (nb/GeV) p p (B0+B0) X default scale = (pT2+mb2) I and F frozen below =mb default scale = (pT2+mb2) I and F frozen below =mb 3 Conclusions Any theoretical prediction for hadronic collisions within perturbative QCD requires the factorization of initial- and final-state singularities. This unavoidably introduces factorization scale parameters, which cannot be predicted from theory. Their choice must be based on physical arguments. We exploited the freedom in this choice to find a prescription that extends the reliability of NLO predictions in the GM-VFNS down to small pT values. With scale parameters I < 1, it is possible to eliminate contributions from the heavy quark in the initial state. These contributions, dominated by the subprocess gb gb, are treated in the parton model with zero quark masses and would lead to an unphysical increase of the cross section for pT 0. We showed, however, that, with a judicious choice of R , I , and F , one can switch off these contributions in the small- pT range. Our prescription leads to a modified GM-VFNS yielding results in good agreement with recent data from the Tevatron and LHC experiments. Acknowledgments We thank I. Bierenbaum for help in checking some of our numerical calculations and J. Rojo for clarifying an issue related to the heavy-quark threshold used in the NNPDF PDFs [63,64]. This work was supported in part by the German Federal Ministry for Education and Research BMBF through Grant No. 05H12GUE. Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited. Funded by SCOAP3 / License Version CC BY 4.0. 1. C. Albajar et al., UA1 Collaboration, Phys. Lett. B 213 , 405 ( 1988 ) 2. F. Abe et al., CDF Collaboration, Phys. Rev. Lett . 75 , 1451 ( 1995 ). hep-ex/ 9503013 3. D.E. Acosta et al., CDF Collaboration, Phys. Rev. D 66 , 052005 ( 2002 ) 4. S. Abachi et al., D0 Collaboration, Phys. Rev. Lett . 74 , 3548 ( 1995 ) 5. B. Abbott et al., D0 Collaboration, Phys. Rev. Lett . 84 , 5478 ( 2000 ). hep-ex/ 9907029 6. D. Acosta et al., CDF Collaboration, Phys. Rev. D 71 , 032001 ( 2005 ). hep-ex/ 0412071 7. A. Abulencia et al., CDF Collaboration, Phys. Rev. D 75 , 012010 ( 2007 ). hep-ex/ 0612015 8. T. Aaltonen et al., CDF Collaboration, Phys. Rev. D 79 , 092003 ( 2009 ). arXiv:0903.2403 [hep-ex] 9. V. Khachatryan et al., CMS Collaboration, Phys. Rev. Lett . 106 , 112001 ( 2011 ). arXiv:1101.0131 [hep-ex] 10. S. Chatrchyan et al., CMS Collaboration, Phys. Rev. Lett . 106 , 252001 ( 2011 ). arXiv:1104.2892 [hep-ex] 11. S. Chatrchyan et al., CMS Collaboration, Phys. Rev. D 84 , 052008 ( 2011 ). arXiv:1106.4048 [hep-ex] 12. G. Aad et al., ATLAS Collaboration , Nucl. Phys . B 864 , 341 ( 2012 ). arXiv:1206.3122 [hep-ex] 13. G. Aad et al., ATLAS Collaboration , JHEP 1310 , 042 ( 2013 ). arXiv:1307.0126 [hep-ex] 14. R. Aaij et al., LHCb Collaboration , JHEP 1204 , 093 ( 2012 ). arXiv:1202.4812 [hep-ex] 15. R. Aaij et al., LHCb Collaboration , JHEP 1308 , 117 ( 2013 ). arXiv:1306.3663 [hep-ex] 16. P. Nason , S. Dawson , R.K. Ellis , Nucl. Phys . B 303 , 607 ( 1988 ) 17. P. Nason , S. Dawson , R.K. Ellis , Nucl. Phys . B 327 , 49 ( 1989 ). [Erratum-ibid . B 335 , 260 ( 1990 )] 18. W. Beenakker , H. Kuijf , W.L. van Neerven , J. Smith, Phys. Rev. D 40 , 54 ( 1989 ) 19. W. Beenakker , W.L. van Neerven , R. Meng , G.A. Schuler , J. Smith, Nucl. Phys . B 351 , 507 ( 1991 ) 20. I. Bojak , M. Stratmann , Phys. Rev. D 67 , 034010 ( 2003 ). hep-ph/ 0112276 21. M. Cacciari , M. Greco , Nucl. Phys . B 421 , 530 ( 1994 ). hep-ph/ 9311260 22. B.A. Kniehl , M. Krmer , G. Kramer , M. Spira , Phys. Lett . B 356 , 539 ( 1995 ). hep-ph/ 9505410 23. M. Cacciari , M. Greco , B.A. Kniehl , M. Krmer , G. Kramer , M. Spira , Nucl. Phys . B 466 , 173 ( 1996 ). hep-ph/ 9512246 24. J. Binnewies , B.A. Kniehl , G. Kramer , Z. Phys . C 76 , 677 ( 1997 ). hep-ph/ 9702408 25. B.A. Kniehl , G. Kramer , M. Spira , Z. Phys . C 76 , 689 ( 1997 ). hep-ph/ 9610267 26. J. Binnewies , B.A. Kniehl , G. Kramer , Phys. Rev. D 58 , 014014 ( 1998 ). hep-ph/ 9712482 27. J. Binnewies , B.A. Kniehl , G. Kramer , Phys. Rev. D 58 , 034016 ( 1998 ). hep-ph/ 9802231 28. B.A. Kniehl , G. Kramer , Phys. Rev. D 60 , 014006 ( 1999 ). hep-ph/ 9901348 29. B.A. Kniehl, in Proceedings of the 14th Topical Conference on Hadron Collider Physics: Hadron Collider Physics 2002 , Karlsruhe, Germany, 2002 , ed. by M. Erdmann, Th. Mller (Springer, Berlin, 2003 ), p. 161 . hep-ph/ 0211008 30. B.A. Kniehl , G. Kramer , Phys. Rev. D 71 , 094013 ( 2005 ). hep-ph/ 0504058 31. B.A. Kniehl , G. Kramer , Phys. Rev. D 74 , 037502 ( 2006 ). hep-ph/ 0607306 32. M. Cacciari , S. Frixione , N. Houdeau , M.L. Mangano , P. Nason , G. Ridolfi , JHEP 1210 , 137 ( 2012 ). arXiv: 1205 .6344 [hep-ph] 33. M. Buza , Y. Matiounine , J. Smith , R. Migneron , W.L. van Neerven, Nucl. Phys . B 472 , 611 ( 1996 ). hep-ph/ 9601302 34. M. Buza , Y. Matiounine , J. Smith , W.L. van Neerven, Eur. Phys. J. C 1 , 301 ( 1998 ). hep-ph/ 9612398 35. G. Kramer , H. Spiesberger , Eur. Phys. J. C 22 , 289 ( 2001 ). hep-ph/ 0109167 36. G. Kramer , H. Spiesberger , Eur. Phys. J. C 28 , 495 ( 2003 ). hep-ph/ 0302081 37. G. Kramer , H. Spiesberger , Eur. Phys. J. C 38 , 309 ( 2004 ). hep-ph/ 0311062 38. B.A. Kniehl , G. Kramer , I. Schienbein , H. Spiesberger , Phys. Rev. D 71 , 014018 ( 2005 ). hep-ph/ 0410289 39. B.A. Kniehl , G. Kramer , I. Schienbein , H. Spiesberger , Eur. Phys. J. C 41 , 199 ( 2005 ). hep-ph/ 0502194 40. B.A. Kniehl , G. Kramer , I. Schienbein , H. Spiesberger , A.I.P. Conf . Proc. 792 , 867 ( 2005 ). hep-ph/ 0507068 41. B.A. Kniehl , G. Kramer , I. Schienbein , H. Spiesberger , Phys. Rev. Lett . 96 , 012001 ( 2006 ). hep-ph/ 0508129 42. B.A. Kniehl , G. Kramer , I. Schienbein , H. Spiesberger , Phys. Rev. D 77 , 014011 ( 2008 ). arXiv: 0705 .4392 [hep-ph] 43. T. Kneesch , B.A. Kniehl , G. Kramer , I. Schienbein , Nucl. Phys . B 799 , 34 ( 2008 ). arXiv: 0712 .0481 [hep-ph] 44. B.A. Kniehl , G. Kramer , I. Schienbein , H. Spiesberger , Phys. Rev. D 79 , 094009 ( 2009 ). arXiv: 0901 .4130 [hep-ph] 45. B.A. Kniehl , G. Kramer , I. Schienbein , H. Spiesberger , Eur. Phys. J. C 62 , 365 ( 2009 ). arXiv: 0902 .3166 [hep-ph] 46. G. Kramer , H. Spiesberger , Phys. Lett . B 679 , 223 ( 2009 ). arXiv: 0906 .2533 [hep-ph] 47. B.A. Kniehl , G. Kramer , I. Schienbein , H. Spiesberger , Phys. Rev. D 84 , 094026 ( 2011 ). arXiv: 1109 .2472 [hep-ph] 48. B.A. Kniehl , G. Kramer , I. Schienbein , H. Spiesberger , Eur. Phys. J. C 72 , 2082 ( 2012 ). arXiv: 1202 .0439 [hep-ph] 49. M.A.G. Aivazis , J.C. Collins , F.I. Olness , W.-K. Tung , Phys. Rev. D 50 , 3102 ( 1994 ). hep-ph/ 9312319 50. S. Kretzer , I. Schienbein , Phys. Rev. D 58 , 094035 ( 1998 ). hep-ph/ 9805233 51. P.M. Nadolsky , H.-L. Lai , Q.-H. Cao , J. Huston , J. Pumplin , D. Stump , W.-K. Tung , C.-P. Yuan , CTEQ Collaboration, Phys. Rev. D 78 , 013004 ( 2008 ). arXiv: 0802 .0007 [hep-ph] 52. K. Nakamura et al., Particle Data Group , J. Phys . G 37 , 075021 ( 2010 ) 53. C. Peterson , D. Schlatter , I. Schmitt , P.M. Zerwas , Phys. Rev. D 27 , 105 ( 1983 ) 54. S. Albino , B.A. Kniehl , G. Kramer , Nucl. Phys . B 803 , 42 ( 2008 ). arXiv: 0803 .2768 [hep-ph] 55. J.C. Collins, Phys. Rev. D 58 , 094002 ( 1998 ). hep-ph/ 9806259 56. W.- K. Tung , S. Kretzer , C. Schmidt , J. Phys . G 28 , 983 ( 2002 ). hep-ph/ 0110247 57. A. Kusina , F.I. Olness , I. Schienbein , T. Jeo , K. Kovark , T. Stavreva , and J.Y. Yu , Phys. Rev. D 88 , 074032 ( 2013 ) arXiv:1306.6553 [hep-ph] 58. M. Cacciari , M. Greco , P. Nason , JHEP 9805 , 007 ( 1998 ). hep-ph/ 9803400 59. F.I. Olness , R.J. Scalise , W.-K. Tung , Phys. Rev. D 59 , 014506 ( 1999 ). hep-ph/ 9712494 60. H.-L. Lai , M. Guzzi , J. Huston , Z. Li , P.M. Nadolsky , J. Pumplin , C.-P. Yuan , CTEQ Collaboration, Phys. Rev. D 82 , 074024 ( 2010 ). arXiv: 1007 .2241 [hep-ph] 61. A.D. Martin , W.J. Stirling , R.S. Thorne , G. Watt , Eur. Phys. J. C 64 , 653 ( 2009 ). arXiv: 0905 .3531 [hep-ph] 62. https://wiki-zeuthen.desy.de/HERAFitter 63. R.D. Ball , V. Bertone , F. Cerutti , L. Del Debbio , S. Forte , A. Guffanti , J.I. Latorre , J. Rojo , M. Ubiali , NNPDF Collaboration, Nucl. Phys . B 849 , 296 ( 2011 ). arXiv: 1101 .1300 [hep-ph] 64. R.D. Ball , V. Bertone , F. Cerutti , L. Del Debbio , S. Forte , A. Guffanti , J.I. Latorre , J. Rojo , M. Ubiali , NNPDF Collaboration, Nucl. Phys . B 855 , 153 ( 2012 ). arXiv: 1107 .2652 [hep-ph]


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-015-3360-6.pdf

B. A. Kniehl, G. Kramer, I. Schienbein, H. Spiesberger. Inclusive \(B\) -meson production at small \(p_T\) in the general-mass variable-flavor-number scheme, The European Physical Journal C, 2015, 140, DOI: 10.1140/epjc/s10052-015-3360-6