Inclusive cross-sections for top-quark pair production in association with charm quarks are measured with proton-proton collision data at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 140 fb$^{-1}$, collected with the ATLAS experiment at the LHC between 2015 and 2018. The measurements are performed by requiring one or two charged leptons (electrons and muons), two $b$-tagged jets, and at least one additional jet in the final state. A custom flavor-tagging algorithm is employed for the simultaneous identification of $b$-jets and $c$-jets. In a fiducial phase space that replicates the acceptance of the ATLAS detector, the cross-sections for $t\bar{t}+ {\geq} 2c$ and $t\bar{t}+1c$ production are measured to be $1.28^{+0.27}_{-0.24}\;\text{pb}$ and $6.4^{+1.0}_{-0.9}\;\text{pb}$, respectively. The measurements are primarily limited by uncertainties in the modeling of inclusive $t\bar{t}$ and $t\bar{t}+b\bar{b}$ production, in the calibration of the flavor-tagging algorithm, and by data statistics. Cross-section predictions from various $t\bar{t}$ simulations are largely consistent with the measured cross-section values, though all underpredict the observed values by 0.5 to 2.0 standard deviations. In a phase-space volume without requirements on the $t\bar{t}$ decay products and the jet multiplicity, the cross-section ratios of $t\bar{t}+ {\geq} 2c$ and $t\bar{t}+1c$ to total $t\bar{t}+\text{jets}$ production are determined to be $(1.23 \pm 0.25) \%$ and $(8.8 \pm 1.3) \%$.
Measured cross-section values in the fiducial phase space and inclusive volume for the various $t\bar{t}+jets$ categories.
Post-fit agreement between data and MC prediction for $SR_{\mathrm{loose}}^{1\ell5j}$ signal region, which uses the invariant mass of the two geometrically closest c-tagged jets, $m_{\mathit{cc}}^{\mathrm{min}\Delta R}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.
Post-fit agreement between data and MC prediction for the $SR_{\mathrm{tight}}^{1\ell5j}$ signal region, which uses the invariant mass of the two geometrically closest jets tagged with c@11%, $m_{\mathit{cc}}^{\mathrm{min}\Delta R}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.
A search is presented for a heavy resonance decaying into a Z boson and a Higgs (H) boson. The analysis is based on data from proton-proton collisions at a centre-of-mass energy of 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$, recorded with the CMS experiment in the years 2016-2018. Resonance masses between 1.4 and 5 TeV are considered, resulting in large transverse momenta of the Z and H bosons. Final states that result from Z boson decays to pairs of electrons, muons, or neutrinos are considered. The H boson is reconstructed as a single large-radius jet, recoiling against the Z boson. Machine-learning flavour-tagging techniques are employed to identify decays of a Lorentz-boosted H boson into pairs of charm or bottom quarks, or into four quarks via the intermediate H $\to$ WW* and ZZ* decays. The analysis targets H boson decays that were not generally included in previous searches using the H $\to$$\mathrm{b\bar{b}}$ channel. Compared with previous analyses, the sensitivity for high resonance masses is improved significantly in the channel where at most one b quark is tagged.
The product of signal acceptance and efficiency for signal events as a function of $m_{Z'}$ for the charged-lepton and neutrino channels in the SR. The efficiency is calculated with respect to Z boson decays to charged leptons and neutrinos for the charged-lepton and neutrino channels, respectively. For comparison, the results from the $\leq$ 1 b category of the previous CMS search in the ZH channel are shown as dashed lines.
Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the muon channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.
Distributions in $m_{Z'}^{T}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the invisible channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.
We report a measurement of the $e^+e^- \to \pi^+\pi^-\pi^0$ cross section in the energy range from 0.62 to 3.50 GeV using an initial-state radiation technique. We use an $e^+e^-$ data sample corresponding to 191 $\text{fb}^{-1}$ of integrated luminosity, collected at a center-of-mass energy at or near the $\Upsilon{(4S)}$ resonance with the Belle II detector at the SuperKEKB collider. Signal yields are extracted by fitting the two-photon mass distribution in $e^+e^- \to \pi^+\pi^-\pi^0\gamma$ events, which involve a $\pi^0 \to \gamma\gamma$ decay and an energetic photon radiated from the initial state. Signal efficiency corrections with an accuracy of 1.6% are obtained from several control data samples. The uncertainty on the cross section at the $\omega$ and $\phi$ resonances is dominated by the systematic uncertainty of 2.2%. The resulting cross sections in the 0.62-1.80 GeV energy range yield $ a_\mu^{3\pi} = [48.91 \pm 0.23~(\mathrm{stat}) \pm 1.07~(\mathrm{syst})] \times 10^{-10} $ for the leading-order hadronic vacuum polarization contribution to the muon anomalous magnetic moment. This result differs by $2.5$ standard deviations from the most precise current determination.
Energy bin range ($\sqrt{s'}$), number of events after unfolding ($N_{\mathrm{unf}}$), corrected efficiency ($\varepsilon$), and cross section ($\sigma_{3\pi}$) for $e^{+}e^{-} \to \pi^{+} \pi^{-} \pi^{0}$ in energy range 0.62--1.05~GeV. The two uncertainties in the cross section are the statistical and systematic contributions. The statistical uncertainties for the unfolding and cross section are square roots of the diagonal components of the unfolding covariance matrix. The image shows Figure 23 in the PRD paper, and the points with error bars indicate the cross section in the table.
Energy bin range ($\sqrt{s'}$), number of events after unfolding ($N_{\mathrm{unf}}$), corrected efficiency ($\varepsilon$), and cross section ($\sigma_{3\pi}$) for $e^{+}e^{-} \to \pi^{+} \pi^{-} \pi^{0}$ in energy range 1.05--3.50~GeV. The two uncertainties in the cross section are the statistical and systematic contributions. The statistical uncertainties for the unfolding and cross section are square roots of the diagonal components of the unfolding covariance matrix. The image shows Figure 23 in the PRD paper, and the points with error bars indicate the cross section in the table.
The statistic covariance matrix for the $e^+e^- \to \pi^+ \pi^- \pi^0$ cross section measurement at the Belle II. The 212 x 212 matrix of the energy ranges from 0.62 to 3.50 GeV. This covariance matrix, obtained by propagating the covariance matrix in the unfolding procedure, shows the total statistical uncertainties for the cross section results.
This paper reports the observation of top-quark pair production in proton-lead collisions in the ATLAS experiment at the Large Hadron Collider. The measurement is performed using 165 nb$^{-1}$ of $p$+Pb data collected at $\sqrt{s_\mathrm{NN}}=8.16$ TeV in 2016. Events are categorised in two analysis channels, consisting of either events with exactly one lepton (electron or muon) and at least four jets, or events with two opposite-charge leptons and at least two jets. In both channels at least one $b$-tagged jet is also required. Top-quark pair production is observed with a significance over five standard deviations in each channel. The top-quark pair production cross-section is measured to be $\sigma_{t\bar{t}}= 58.1\pm 2.0\;\mathrm{(stat.)\;^{+4.8}_{-4.4} \;\mathrm{(syst.)}}\;\mathrm{nb}$, with a total uncertainty of 9%. In addition, the nuclear modification factor is measured to be $R_{p\mathrm{A}} = 1.090\pm0.039\;(\mathrm{stat.})\;^{+0.094}_{-0.087}\;(\mathrm{syst.})$. The measurements are found to be in good agreement with theory predictions involving nuclear parton distribution functions.
The figure shows the pre-fit distribution of events as a function of $H_{\mathrm{T}}^{\ell,j} = \sum_{\ell,j} p_{T}^{\ell,j}$, scalar sum of $p_T$ for all jets and leptons in the $\ell+$jets channel, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, with an integrated luminosity of 165 nb$^{-1}$. The data correspond to the $1\ell 1b$ $e$+jets channel in a pre-fit configuration. The stacked histograms represent different processes contributing to the event yield, including top quark pair production ($t\bar{t}$), single top, $W$ boson production with $b$, $c$, and light quarks, $Z$ boson production with $b$, $c$, and light quarks, diboson, and fake lepton backgrounds.
The figure shows the post-fit distribution of events as a function of $H_{\mathrm{T}}^{\ell,j} = \sum_{\ell,j} p_{T}^{\ell,j}$, scalar sum of $p_T$ for all jets and leptons in the $\ell+$jets channel, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, with an integrated luminosity of 165 nb$^{-1}$. The data correspond to the $1\ell 1b$ $e$+jets channel in a pre-fit configuration. The stacked histograms represent different processes contributing to the event yield, including top quark pair production ($t\bar{t}$), single top, $W$ boson production with $b$, $c$, and light quarks, $Z$ boson production with $b$, $c$, and light quarks, diboson, and fake lepton backgrounds.
The figure shows the pre-fit distribution of events as a function of $H_{\mathrm{T}}^{\ell,j} = \sum_{\ell,j} p_{T}^{\ell,j}$, scalar sum of $p_T$ for all jets and leptons in the $\ell+$jets channel, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, with an integrated luminosity of 165 nb$^{-1}$. The data correspond to the $1\ell 1b$ $\mu$+jets channel in a pre-fit configuration. The stacked histograms represent different processes contributing to the event yield, including top quark pair production ($t\bar{t}$), single top, $W$ boson production with $b$, $c$, and light quarks, $Z$ boson production with $b$, $c$, and light quarks, diboson, and fake lepton backgrounds.
A search is presented for the pair production of new heavy resonances, each decaying into a top quark (t) or antiquark and a gluon (g). The analysis uses data recorded with the CMS detector from proton-proton collisions at a center-of-mass energy of 13 TeV at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. Events with one muon or electron, multiple jets, and missing transverse momentum are selected. After using a deep neural network to enrich the data sample with signal-like events, distributions in the scalar sum of the transverse momenta of all reconstructed objects are analyzed in the search for a signal. No significant deviations from the standard model prediction are found. Upper limits at 95% confidence level are set on the product of cross section and branching fraction squared for the pair production of excited top quarks in the $\mathrm{t^*}$ $\to$ tg decay channel. The upper limits range from 0.12 pb to 0.8 fb for a $\mathrm{t^*}$ with spin-1/2 and from 0.015 pb to 1.0 fb for a $\mathrm{t^*}$ with spin-3/2. These correspond to mass exclusion limits up to 1050 and 1700 GeV for spin-1/2 and spin-3/2 $\mathrm{t^*}$ particles, respectively. These are the most stringent limits to date on the existence of $\mathrm{t^*}$ $\to$ tg resonances.
Expected and observed 95% CL upper limits on the product of the $t^{*} \overline{t}^{*}$ production cross section and the branching fraction squared $BR^2(t^{*} \rightarrow tg)$ for a spin-1/2 $t^{*}$ as a function of $m_{t^{*}}$. The inner (green) and outer (yellow) bands give the central probability intervals containing 68 and 95% of the expected upper limits under the background-only hypothesis. The cross section predicted by theory, following an EFT approach, is shown in blue, assuming $BR(t^{*} \rightarrow tg)=1$.
Expected and observed 95% CL upper limits on the product of the $t^{*} \overline{t}^{*}$ production cross section and the branching fraction squared $BR^2(t^{*} \rightarrow tg)$ for a spin-3/2 $t^{*}$ as a function of $m_{t^{*}}$. The inner (green) and outer (yellow) bands give the central probability intervals containing 68 and 95% of the expected upper limits under the background-only hypothesis. The cross section predicted by theory, following an EFT approach, is shown in blue, assuming $BR(t^{*} \rightarrow tg)=1$. The results of the previous CMS analysis, using data corresponding to an integrated luminosity of 35.9 $fb^{-1}$, are shown in red.
Distributions in $S_T$ in the SR for the muon channel, after a background-only fit to the data. The signal distributions are scaled to the cross section predicted by the theory. The hatched bands show the post-fit uncertainty band, combining all sources of uncertainty. The ratio of data to the background predictions is shown in the panels below the distributions.
A search is presented for the pair production of new heavy resonances, each decaying into a top quark or antiquark and a gluon. The analysis uses data recorded with the CMS detector from proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 $\text{fb}^{-1}$. Events with one muon or electron, multiple jets, and missing transverse momentum are selected. After using a deep neural network to enrich the data sample with signal-like events, distributions in the scalar sum of the transverse momenta of all reconstructed objects are analyzed in search for a signal. No significant deviations from the standard model predictions are found. Upper limits at $95\%$ confidence level are set on the product of cross section times branching fraction squared for the pair production of two excited top quarks in the $\mathrm{t}^{*} \to \mathrm{t}\mathrm{g}$ decay channel. The upper limits range from 0.12 pb to 0.8 fb for a $\text{t}^{*}$ with spin-1/2 and from 0.015 pb to 1.0 fb for a $\text{t}^{*}$ with spin-3/2. This corresponds to mass exclusion limits up to 1050 and 1700 GeV for spin-1/2 and spin-3/2 $\mathrm{t}^{*}$ particles, respectively.
Expected and observed 95% CL upper limits on the product of the $t^{*} \overline{t}^{*}$ production cross section and the branching fraction squared $BR^2(t^{*} \rightarrow tg)$ for a spin-1/2 $t^{*}$ as a function of $m_{t^{*}}$. The inner (green) and outer (yellow) bands give the central probability intervals containing 68 and 95% of the expected upper limits under the background-only hypothesis. The cross section predicted by theory, following an EFT approach, is shown in blue, assuming $BR(t^{*} \rightarrow tg)=1$.
Expected and observed 95% CL upper limits on the product of the $t^{*} \overline{t}^{*}$ production cross section and the branching fraction squared $BR^2(t^{*} \rightarrow tg)$ for a spin-1/2 $t^{*}$ as a function of $m_{t^{*}}$. The inner (green) and outer (yellow) bands give the central probability intervals containing 68 and 95% of the expected upper limits under the background-only hypothesis. The cross section predicted by theory, following an EFT approach, is shown in blue, assuming $BR(t^{*} \rightarrow tg)=1$.
Expected and observed 95% CL upper limits on the product of the $t^{*} \overline{t}^{*}$ production cross section and the branching fraction squared $BR^2(t^{*} \rightarrow tg)$ for a spin-3/2 $t^{*}$ as a function of $m_{t^{*}}$. The inner (green) and outer (yellow) bands give the central probability intervals containing 68 and 95% of the expected upper limits under the background-only hypothesis. The cross section predicted by theory, following an EFT approach, is shown in blue, assuming $BR(t^{*} \rightarrow tg)=1$. The results of the previous CMS analysis, using data corresponding to an integrated luminosity of 35.9 $fb^{-1}$, are shown in red.
We report the results of a test of the compositeness of the heaviest charged lepton, $\tau$, using data collected in proton-proton collisions at a center of mass energy of 13 TeV at the Large Hadron Collider at CERN in 2016-2018, corresponding to an integrated luminosity of $138\,\mathrm{fb}^{-1}$. We seek tau lepton pair production in which one tau lepton is produced in an excited state and decays to a ground state tau lepton and photon. The analysis selects events with two isolated tau lepton decay candidates plus a high energy photon. The mass of the excited tau lepton is reconstructed using the missing transverse momentum in the event, assuming the momentum of the neutrinos from each tau lepton decay lie along the visible tau lepton decay product direction. We observe no excess of events above the standard model background prediction and use this null result to set lower bounds on the excited tau lepton mass. For a compositeness scale $\Lambda$ equal to the excited tau lepton mass (10 TeV), we exclude excited tau leptons with masses below 4700 GeV (2800 GeV).
Signal process cutflow for all taustar hypothesis masses, channels, and years. Cut descriptions: - numLeps : 1e0mu for ETau, 0e1mu For muTau, 0e0mu for TauTau - havePair : 1e + 1tauh for ETau, 1mu + 1tauh for MuTau, 2tauh for TauTau - trigger : passed appropriate high level trigger (HLT) for channel and year - m_vis : Visible mass of pair is >= 100 GeV - Zveto : No Z->ee or Z->mumu pairs - bveto : No btagged jets in event (only for ETau and MuTau) - gamma20GeV : have a photon with pt>20GeV - gamma100GeV : have a photon with pt>100GeV
Background-only-fit total covariance matrix for all channels and regions. Values are extracted from the Higgs Combine tool's fitDiagnostics option. Please see associated image ("fitRegions.pdf") for guide on region and channel setup.
Background-only-fit total covariance matrix for all channels and regions. Values are extracted from the Higgs Combine tool's fitDiagnostics option. Please see associated image ("fitRegions.pdf") for guide on region and channel setup.
The Large Hadron Collider at CERN, delivering proton-proton collisions at much higher energies and far higher luminosities than previous machines, has enabled a comprehensive programme of measurements of the standard model (SM) processes by the CMS experiment. These unprecedented capabilities facilitate precise measurements of the properties of a wide array of processes, the most fundamental being cross sections. The discovery of the Higgs boson and the measurement of its mass became the keystone of the SM. Knowledge of the mass of the Higgs boson allows precision comparisons of the predictions of the SM with the corresponding measurements. These measurements span the range from one of the most copious SM processes, the total inelastic cross section for proton-proton interactions, to the rarest ones, such as Higgs boson pair production. They cover the production of Higgs bosons, top quarks, single and multibosons, and hadronic jets. Associated parameters, such as coupling constants, are also measured. These cross section measurements can be pictured as a descending stairway, on which the lowest steps represent the rarest processes allowed by the SM, some never seen before.
Cross sections of selected high-energy processes measured by the CMS experiment. Measurements performed at different LHC pp collision energies are marked by unique symbols and the coloured bands indicate the combined statistical and systematic uncertainty of the measurement. Grey bands indicate the uncertainty of the corresponding SM theory predictions. Shaded hashed bars indicate the excluded cross section region for a production process with the measured 95% CL upper limit on the process indicated by the solid line of the same colour.
Summary of production cross section measurements involving top quarks. Measurements performed at different LHC pp collision energies are marked by unique symbols and the coloured bands indicate the combined statistical and systematic uncertainty of the measurement. Grey bands indicate the uncertainty of the corresponding SM theory predictions. Shaded hashed bars indicate the excluded cross section region for a production process with the measured 95% C.L. upper limit on the process indicated by the solid line of the same colour.
Summary of measurements of jet cross sections and electroweak processes in association with jets. Measurements performed at different LHC pp collision energies are marked by unique symbols and the coloured bands indicate the combined statistical and systematic uncertainty of the measurement. Grey bands indicate the uncertainty of the corresponding SM theory predictions. Shaded hashed bars indicate the excluded cross section region for a production process with the measured 95% C.L. upper limit on the process indicated by the solid line of the same colour. Versions of these plots in pdf format with links to the publications can be found at https://cms-results.web.cern.ch/cms-results/public-results/publications/SMP-23-004/.
This paper presents a measurement of the production cross-section of a $Z$ boson in association with $b$- or $c$-jets, in proton-proton collisions at $\sqrt{s} = 13$ TeV with the ATLAS experiment at the Large Hadron Collider using data corresponding to an integrated luminosity of 140 fb$^{-1}$. Inclusive and differential cross-sections are measured for events containing a $Z$ boson decaying into electrons or muons and produced in association with at least one $b$-jet, at least one $c$-jet, or at least two $b$-jets with transverse momentum $p_\textrm{T} > 20$ GeV and rapidity $|y| < 2.5$. Predictions from several Monte Carlo generators based on next-to-leading-order matrix elements interfaced with a parton-shower simulation, with different choices of flavour schemes for initial-state partons, are compared with the measured cross-sections. The results are also compared with novel predictions, based on infrared and collinear safe jet flavour dressing algorithms. Selected $Z + \ge 1 c$-jet observables, optimized for sensitivity to intrinsic-charm, are compared with benchmark models with different intrinsic-charm fractions.
Figure 6(left) of the article. Measured fiducial cross sections for events with $Z (\rightarrow ll) \ge 1 $ b-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.
Figure 6(right) of the article. Measured fiducial cross sections for events with $Z (\rightarrow ll) \ge 2 $ b-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.
Figure 7 of the article. Measured fiducial cross sections for events with $Z (\rightarrow ll) \ge 1 $ c-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.
This article reports on the inclusive production cross section of several quarkonium states, $\mathrm{J}/\psi$, $\psi {\rm (2S)}$, $\Upsilon\rm(1S)$, $\Upsilon\rm(2S)$, and $\Upsilon\rm(3S)$, measured with the ALICE detector at the LHC, in pp collisions at $\sqrt{s} = 5.02$ TeV. The analysis is performed in the dimuon decay channel at forward rapidity ($2.5 < y < 4$). The integrated cross sections and transverse-momentum ($p_{\rm T}$) and rapidity ($y$) differential cross sections for $\mathrm{J}/\psi$, $\psi {\rm (2S)}$, $\Upsilon\rm(1S)$, and the $\psi {\rm (2S)}$-to-$\mathrm{J}/\psi$ cross section ratios are presented. The integrated cross sections, assuming unpolarized quarkonia, are: $\sigma_{\mathrm{J}/\psi}$($p_{\rm T}<20$ GeV/c) = 5.88 $\pm$ 0.03 $\pm$ 0.34 $\mu$b, $\sigma_{\psi {\rm (2S)}}$($p_{\rm T}<12$ GeV/c) = 0.87 $\pm$ 0.06 $\pm$ 0.10 $\mu$b, $\sigma_{\Upsilon\rm(1S)}$($p_{\rm T}<15$ GeV/c) = 45.5 $\pm$ 3.9 $\pm$ 3.5 nb, $\sigma_{\Upsilon\rm(2S)}$($p_{\rm T}<15$ GeV/c) = 22.4 $\pm$ 3.2 $\pm$ 2.7 nb, and $\sigma_{\Upsilon\rm(3S)}$($p_{\rm T}<15$ GeV/c) = 4.9 $\pm$ 2.2 $\pm$ 1.0 nb, where the first (second) uncertainty is the statistical (systematic) one. For the first time, the cross sections of the three $\Upsilon$ states, as well as the $\psi {\rm (2S)}$ one as a function of $p_{\rm T}$ and $y$, are measured at $\sqrt{s} = 5.02$ TeV at forward rapidity. These measurements also significantly extend the $\mathrm{J}/\psi$$p_{\rm T}$ reach and supersede previously published results. A comparison with ALICE measurements in pp collisions at $\sqrt{s} = 2.76$, 7, 8, and 13 TeV is presented and the energy dependence of quarkonium production cross sections is discussed. Finally, the results are compared with the predictions from several production models.
Differential production cross sections of J/$\psi$ as a function of $p_{\rm T}$.
Differential production cross sections of J/$\psi$ as a function of rapidity.
Rapidity Differential production cross sections of J/$\psi$ in the range 0.3 $\leq$ p$_{\rm T}$ < 2 GeV/c