$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 20--40\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 40--60\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 60--100\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to charged $\pi$ ($\pi^{+} + \pi^{-}$) at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to charged $\pi$ ($\pi^{+} + \pi^{-}$) as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to K$^{-}$ at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to K$^{-}$ as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to its ground state particle $\Lambda$ ($\Lambda + \bar{\Lambda}$) at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to its ground state particle $\Lambda$ ($\Lambda + \bar{\Lambda}$) as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Expected and observed events in the 16 discovery regions along with the according control regions.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic R-hadron search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector R-hadron search.

p0-values and model-independent upper limits on cross-section x acceptance x efficiency for the 16 discovery regions.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic search for metastable gluino R-hadrons.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector direct-stau search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector chargino search.

Upper cross-section limit in gluino R-hadron search.

Upper cross-section limit in sbottom R-hadron search.

Upper cross-section limit in stop R-hadron search.

Upper cross-section limit in stau search.

Upper cross-section limit in chargino search.

Lower mass limit as function of gluino lifetime.

Acceptance x efficiency, acceptance and efficiency for the full range of simulated masses in the MS-agnostic R-hadron search.

Upper cross-section limit in meta-stable gluino R-hadron search.

Flavor composition of 800 GeV stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

ETmiss trigger efficiency as function of true ETmiss (EtmissTurnOn).

Single-muon trigger efficiency as function of $|\eta|$ and $\beta$ (SingleMuTurnOn).

Candidate reconstruction efficiency for ID+Calo selection (IDCaloEff).

Candidate reconstruction efficiency for loose selection (LooseEff).

Efficiency for a loose candidate to be promoted to a tight candidate (TightPromotionEff).

Resolution and average of reconstructed dE/dx mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for FullDet candidates.

Predicted ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<1.37$.

Measured ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Predicted ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.37$.

Predicted ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.37$.

Measured double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$.

Predicted double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$.

Measured double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<1.37$.

Predicted double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<1.37$.

Measured double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Predicted double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.37$.

Predicted double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.37$.

Acceptance x efficiency versus resonance mass for both lephad and hadhad channels in the RS bulk model with k/MPl = 1

Acceptance x efficiency versus resonance mass for both lephad and hadhad channels in the RS bulk model with k/MPl = 2

Acceptance x efficiency versus resonance mass for both lephad and hadhad channels in the scalar model

Upper limits on the production cross-section times the HH to bbtautau branching ratio for non-resonant HH at 95% CLS and their interpretation as multiples of the SM prediction

Upper limits on the production cross-section times the HH to bbtautau branching ratio divided by the SM prediction for non-resonant HH at 95% CL

Post-fit expected number of signal and background events and observed number of data events after applying the selection criteria and requiring exactly 2 b-tagged jets and assuming a background-only hypothesis

Post-fit expected number of signal and background events and observed number of data events in the last two bins of the non-resonant BDT score distribution of the SM signal after applying the selection criteria and requiring exactly 2 b-tagged jets and assuming a background-only hypothesis

Dijet mass spectra after the background only fit with the background prediction in the high-mass region with two b-tags.

Dijet mass spectra after the background only fit with the background prediction in the low-mass region with two b-tags.

The online b-tagging efficiency with respect to the offline b-tagging efficiency as a function of pT. The b-tagging online and offline working points correspond to an efficiency of 60% and 70%, respectively.

Observed and expected 95% credibility-level upper limits on the cross-section for the b* model in the high-mass region with inclusive b-jet selection.

Observed and expected 95% credibility-level upper limits on the cross-section times branching ratio for the SSM and leptophobic Z' models in the low- and high-mass region with two b-tags selection.

Observed and expected 95% credibility-level upper limits on the cross-section for DM Z' models in the low-mass region with two b-tags selection. The Z' is expected to decay to all five quark flavors other than the top quark and the mediator to SM quark coupling (gSM) equal to 0.1 is assumed.

Observed and expected 95% credibility-level upper limits on the cross-section times branching ratio for DM Z'->bb models in the high-mass region with two b-tags selection. The Z' is expected to decay to bb only and the mediator to SM quark coupling (gSM) equal to 0.25 is assumed.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances with intrinsic width smaller than the detector resolution. The width of the Gaussian reconstructed shape is dominated by the dijet mass resolution. The table shows the limits obtained from the high-mass inclusive one b-tag selection.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances with intrinsic width smaller than the detector resolution. The width of the Gaussian reconstructed shape is dominated by the dijet mass resolution. The table shows the limits obtained from the combined low- and high-mass two b-tags selection.

The mass distributions for the inclusive one b-tagged selection and two b-tagged selection using an integrated luminosity of 36.1 $fb^{-1}$. The inclusive one b-tagged Pythia8 MC distribution is normalized to the inclusive one b-tagged data. The two b-tagged Pythia8 MC distribution is normalized to the two b-tagged data. The systematic uncertainty band is dominated by the b-tagging scale factor and the b-jet energy scale.

Signal acceptance times efficiency for inclusive 1 b-tag and 2 b-tag categories as a function of the simulated b* and the Z' masses.

Signal acceptance times efficiency for 2 b-tag categories as a function of the simulated Z' masses.

The flavor composition of the simulated dijet background as a function of dijet mass before tagging. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

The flavor composition of the simulated dijet background as a function of dijet mass with inclusive one b-tag. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

The flavor composition of the simulated dijet background as a function of dijet mass with two b-tags. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 3%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 7%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 10%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 15%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 3%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 7%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 10%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 15%.

The ratio $|C/A|$ of amplitudes of $\omega$ and $\rho^0$ mesons.

$d\sigma/dy$ for exclusively photoproduced $\rho^0$ mesons in $Xn\,Xn$ events.

$d\sigma/dy$ for exclusively photoproduced $\rho^0$ mesons in the previous STAR analysis, Phys. Rev. C 77 034910 (2008).

$d\sigma/dy$ for exclusively photoproduced $\rho^0$ mesons in $1n\,1n$ events.

The $-t$ distribution for exclusive $\rho^0$ mesons in events with $Xn\,Xn$ mutual event dissociation.

The $-t$ distribution for exclusive $\rho^0$ mesons in events with $1n\,1n$ mutual event dissociation.

$d\sigma/dt$ for coherent $\rho^0$ photoproduction in $Xn\,Xn$ events

$d\sigma/dt$ for coherent $\rho^0$ photoproduction in $1n\,1n$ events

The target distribution in the transverse plane, the result of a two-dimensional Fourier transform (Hankel transform) of the $1n\, 1n$ diffraction pattern shown in Fig. 8. The integration is limied to the region $|t| < 0.06 (\textrm{GeV}/c)^2$. The uncertainty is estimated by changing the maximum $-t$ to 0.05, 0.07, and 0.09 $(\textrm{GeV}/c)^2$

The target distribution in the transverse plane, the result of a two-dimensional Fourier transform (Hankel transform) of the $Xn\, Xn$ diffraction pattern shown in Fig. 8. The integration is limied to the region $|t| < 0.06 (\textrm{GeV}/c)^2$. The uncertainty is estimated by changing the maximum $-t$ to 0.05, 0.07, and 0.09 $(\textrm{GeV}/c)^2$

$e^+$ $e^-$ invariant mass pair distributions of signal pairs compared to the inclusive unlike sign and reconstructed background pairs in 200 GeV Au+Au minimum-bias collisions.

Signal-to-background ratios in minimum bias p+p and Au+Au collisions.

Signal-to-background ratios in central p+p and Au+Au collisions.

Invariant mass spectrum in the STAR acceptance ($p_T^e$ > 0.2 GeV/c, |$\eta^e$| < 1, and |$y_{ee}$|) from SQRT($s_{NN}$) = 200 GeV Au+Au minimum-bias collisions.

Invariant mass spectrum in the STAR acceptance ($p_T^e$ > 0.2 GeV/c, |$\eta^e$| < 1, and |$y_{ee}$|) from SQRT($s_{NN}$) = 200 GeV Au+Au minimum-bias collisions.

Invariant mass spectra from SQRT($s_{NN}$ = 200 GeV Au+Au minimum-bias collisions in pT range 0 - 5 GeV/c. The ratio of dielectron yield over cocktail for different pT bins, and the comparison with model calculations.

Invariant mass spectra from SQRT($s_{NN}$ = 200 GeV Au+Au minimum-bias collisions in pT range 5 - 10 GeV/c. The ratio of dielectron yield over cocktail for different pT bins, and the comparison with model calculations.

Invariant mass spectra from SQRT($s_{NN}$ = 200 GeV Au+Au minimum-bias collisions in pT range 1 - 1.5 GeV/c. The ratio of dielectron yield over cocktail for different pT bins, and the comparison with model calculations.

Invariant mass spectra from SQRT($s_{NN}$ = 200 GeV Au+Au minimum-bias collisions in pT range 1.5 - 2 GeV/c. The ratio of dielectron yield over cocktail for different pT bins, and the comparison with model calculations.

Invariant mass spectra from SQRT($s_{NN}$ = 200 GeV Au+Au minimum-bias collisions integrated. The ratio of dielectron yield over cocktail for different pT bins, and the comparison with model calculations.

The integrated dielectron yield as a function of pair pT in invariant mass range 0 - 0.15 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

The integrated dielectron yield as a function of pair pT in invariant mass range 0.15 - 0.3 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

The integrated dielectron yield as a function of pair pT in invariant mass range 0.3 - 0.76 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

The integrated dielectron yield as a function of pair pT in invariant mass range 0.76 - 1.05 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

The integrated dielectron yield as a function of pair pT in invariant mass range 1.05 - 1.8 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

The integrated dielectron yield as a function of pair pT in invariant mass range 1.8 - 2.8 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

The integrated dielectron yield as a function of pair pT in invariant mass range 2.8 to 3.5 GeV/c^2 compared with cocktail. The ratio of dielectron yield over cocktail for different mass ranges as a function of pair pT.

Invariant mass spectra from SQRT($s_{NN}$) = 200 GeV Au+Au collisions in different centralities. The ratio of dielectron yield over cocktail for a centrality of 0-10%.

Invariant mass spectra from SQRT($s_{NN}$) = 200 GeV Au+Au collisions in different centralities. The ratio of dielectron yield over cocktail for a centrality of 10-40%.

Invariant mass spectra from SQRT($s_{NN}$) = 200 GeV Au+Au collisions in different centralities. The ratio of dielectron yield over cocktail for a centrality of 40-80%.

Invariant mass spectra from SQRT($s_{NN}$) = 200 GeV Au+Au collisions in different centralities. The ratio of dielectron yield over cocktail for minimum bias.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 0 - 0.15 GeV/c^2, and the comparison with model calculations.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 0.15 - 0.3 GeV/c^2, and the comparison with model calculations.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 0.3 - 0.76 GeV/c^2, and the comparison with model calculations.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 0.76 - 1.05 GeV/c^2, and the comparison with model calculations.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 1.05 - 1.8 GeV/c^2, and the comparison with model calculations.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 1.8 - 2.8 GeV/c^2, and the comparison with model calculations.

The integrated dielectron yield and ratio of dielectron yield over cocktail in different centralities for mass window 2.8 - 3.5 GeV/c^2, and the comparison with model calculations.

Dielectron invariant mass spectra from minimum-bias.

Dielectron invariant mass spectra from the most central (0-10%) collisions that we are able to achieve most statistics at present.

The ratio of the Npart-scaled dielectron yield between minimum-bias and the most central collisions.

Mass spectrum of the excess (data - cocktail) in the low-mass region in Au+Au minimum-bias collisions compared to model calculations.

The yields scaled by Npart for the $\rho$-like region with the cocktail subtracted.

The Omega-like region without cocktail subtraction as a function of Npart.

The Phi-like region without cocktail subtraction as a function of Npart.

Correlated charm contributions to the dielectron mass spectra for different assumptions of the correlation strength.

Slope parameters Teff versus invariant mass for dielectrons from charm hadron decays.

Omega meson invariant mass distribution from SQRT($s_{NN}$) = 200 GeV Au+Au minimum-bias collisions after subtraction of the combinatorial background using the mixed-event method.

Phi meson invariant mass distribution from SQRT($s_{NN}$) = 200 GeV Au+Au minimum-bias collisions after subtraction of the combinatorial background using the mixed-event method.

The widths and mass positions of the omega signal from data compared to the values from the PDG and the full Geant simulation.

The widths and mass positions of the phi signal from data compared to the values from the PDG and the full Geant simulation.

The efficiency and acceptance correction factor as function of pT for mid-rapidity omega and phi mesons.

The pT distributions of the omega meson invariant yields from SQRT($s_{NN}$) = 200 GeV Au+Au minimum-bias collisions.

The pT distributions of the phi meson invariant yields from SQRT($s_{NN}$) = 200 GeV Au+Au minimum-bias collisions.

Unlike-sign/like-sign pair acceptance difference correction factor with the PHENIX phi acceptance compared with the full acceptance.

Efficiency corrected invariant mass spectra calculated using the STAR data filtered with the PHENIX azimuthal angle acceptance. The data points are compared to cocktail simulations.

The same data points compared to theoretical model calculations of medium vector meson and QGP contributions.

Relative beauty contribution to the heavy-flavour electron yield obtained from the method based on the azimuthal correlation of heavy-flavour decay electrons and charged tracks.

The pT-differential inclusive cross section of electrons from beauty hadron decays using the electron-hadron azimuthal correlation method.

The pT-differential inclusive cross section of electrons from beauty hadron decays. A 1.9% normalization uncertainty not shown.

Inclusive bbbar production cross section per rapidity unit measured.

Total bbbar production cross section.