FONLL Heavy Quark Production

JavaScript is disabled in your browser. Please enable it in order to use this form.

Calculate total or single inclusive differential cross sections for heavy quark production (charm or bottom) at ppbar or pp colliders, with cuts on p_T and y or eta. Non perturbative fragmentation into heavy hadrons and their subsequent decays into other final states (leptons, quarkonia, ...) can also be included.

NB. Use of the "back" button in the browser to come back to the form after having produced a set of results can occasionally lead to unexpected behaviour, and is deprecated. Reloading the form from scratch is safer and preferable, until this bug is fixed.

NEWS:

• 28 Aug 2024 - Addition of 16.8 GeV CM energy
• 06 Aug 2024 - Addition of 5.36 TeV CM energy and of NNPDF30nlo_as0118 set for 5.5 TeV
• 05 Dec 2022 - Changed default PDF set to NNPDF30nlo_as0118, where available. Added 5.03 TeV with this PDF set.
• 03 Aug 2022 - Addition of 13.6 TeV CM energy
• 11 Feb 2022 - Addition of new CM energies, 8.16 TeV and 68.6 GeV
• 11 Oct 2018 - Addition of double-differential output options ds/dpt/dy and ds/dpt/deta

Be sure to check the older news, as well as the list of bugs discovered so far. Feedback to matteo.cacciari@cern.ch is welcome.

Collider: pp, 16.8 GeV    [ptmax ≤ 3 GeV] RHIC (pp, 62.4 GeV)    [ptmax ≤ 12 GeV] pp, 68.6 GeV    [ptmax ≤ 12 GeV] RHIC (pp, 200 GeV)    [ptmax ≤ 30 GeV] RHIC (pp, 500 GeV)    [ptmax ≤ 50 GeV] Tevatron (ppbar, 900 GeV)    [ptmax ≤ 50 GeV] Tevatron Run 2 (ppbar, 1960 GeV)    [ptmax ≤ 100 GeV] LHC (pp, 2.75 TeV)    [ptmax ≤ 300 GeV] LHC (pp, 4+1.58 TeV)    [ptmax ≤ 300 GeV] LHC (pp, 5.03 TeV)    [ptmax ≤ 300 GeV] LHC (pp, 5.36 TeV)    [ptmax ≤ 300 GeV] LHC (pp, 5.5 TeV)    [ptmax ≤ 300 GeV] LHC (pp, 7 TeV)    [ptmax ≤ 300 GeV] LHC (pp, 8 TeV)    [ptmax ≤ 300 GeV] LHC (pp, 8.16 TeV)    [ptmax ≤ 300 GeV] LHC (pp, 13 TeV)    [ptmax ≤ 300 GeV] LHC (pp, 13.6 TeV)    [ptmax ≤ 300 GeV] LHC (pp, 14 TeV)    [ptmax ≤ 300 GeV]

PDFs: CTEQ6.6 NNPDF30nlo_as0118

Perturbative order: FONLL NLO

Heavy quark: charm bottom

Hadronic final state: bare quark B hadron

Weak decay: - B hadron to electron B to D to electron B to J/psi B to psi(2S) B to D B to D*

Cross section type: Total cross section (within cuts) dsigma/dpt dsigma/dy dsigma/deta 1/pi dsigma/dpt^2/dy|_y=0 dsigma/dpt/dy dsigma/dpt/deta

central prediction only uncertainty range from scales only uncertainty range from scales and masses

Output all scales

Include PDFs uncertainties

ptmin (GeV)

ptmax (GeV)

y(eta)min

y(eta)max

Use y Use eta

npoints

Change these defaults only if you know what you are doing

Non-pert. Fragm. Funct.: default Kartvelishvili Peterson BCFY for D* BCFY for D0 BCFY for D+

Non-pert. parameter =

Use default N=5 moment

Fragmentation Fraction c->D =

Fragmentation Fraction b->B =

These fragmentation fraction values are generic. They can (and should) be set manually for a specific hadron.

BR(D->l) =

BR(B->l) =

BR(B->D->l) =

BR(B->D) =

BR(B->D*) =

BR(B->J/psi) =

BR(B->psi(2S)) =

Set you inputs, click submit, and be patient. The run time may be as long as a few minutes depending on your choices (e.g. a total cross section with full uncertainty band, or many points of a differential cross section).
Display results Email results to

Submit   Reset

Notes
Central values: m_b = 4.75 GeV for bottom, m_c = 1.5 GeV for charm, μ_R = μ_F = μ₀ = sqrt(m² + p_T²)
Scales uncertainties: μ₀/2 < μ_R, μ_F < 2μ₀ with 1/2 < μ_R/μ_F < 2.
Mass uncertainties: m_b = 4.5, 5.0 GeV for bottom, m_c = 1.3, 1.7 GeV for charm, summed in quadrature to scales uncertainties.
PDFs uncertainties: calculated according to the individual PDF set recipe, and summed in quadrature to scales and mass uncertainties.

No fragmentation fractions (unless specified above) are included for the heavy quark -> heavy hadron fragmentation.
This means that all heavy quarks are hadronised as if they fragmented into the chosen heavy hadron.
To construct the proper mixing, the correct fragmentation fraction (FF) must be provided and the results summed separately.
The D0 and D+ already include feedown from D*.
The correct branching ratios (BR) for weak decays into leptons and other hadrons are instead provided by default.

Last modified:    © 2007-2024 Matteo Cacciari