Measurement of psi(2S) production cross-sections in proton-proton collisions at 7 and 13 TeV

TL;DR

This paper measures the production cross-sections of psi(2S) mesons in proton-proton collisions at 7 and 13 TeV, providing new data that aligns well with theoretical predictions and enhances understanding of quarkonium production.

Contribution

First measurements of psi(2S) production cross-sections at 7 and 13 TeV in specific kinematic regions using LHCb data, improving the experimental dataset for quarkonium production.

Findings

Measured cross-sections at 13 TeV: prompt 1.43 μb, from-b 0.43 μb.

Measured cross-sections at 7 TeV: prompt 0.47 μb, from-b 0.13 μb.

Results agree with theoretical models.

Abstract

The cross-sections of $\psi(2S)$ meson production in proton-proton collisions at $\sqrt{s}=13~\mathrm{TeV}$ are measured with a data sample collected by the LHCb detector corresponding to an integrated luminosity of $275~p\mathrm{b}^{-1}$. The production cross-sections for prompt $\psi(2S)$ mesons and those for $\psi(2S)$ mesons from $b$-hadron decays ($\psi{(2S)}\mathrm{-from-}b$) are determined as functions of the transverse momentum, $p_{\mathrm{T}}$, and the rapidity, $y$, of the $\psi(2S)$ meson in the kinematic range $2<p_{\mathrm{T}}<20~\mathrm{GeV}/c$ and $2.0<y<4.5$. The production cross-sections integrated over this kinematic region are \begin{equation*} \begin{split} \sigma(\mbox{prompt }\psi(2S),13~\mathrm{TeV}) &= {1.430 \pm 0.005(\mathrm{stat}) \pm 0.099 (\mathrm{syst})\mu\mathrm{b}},\\ \sigma(\psi(2S)\mathrm{-from-}b,13~\mathrm{TeV})&={0.426 \pm 0.002(\mathrm{stat}) \pm0.030 (\mathrm{syst})\mu\mathrm{b}}. \end{split} \end{equation*} A new measurement of $\psi(2S)$ production cross-sections in $pp$ collisions at $\sqrt{s}=7~\mathrm{TeV}$ is also performed using data collected in 2011, corresponding to an integrated luminosity of $614~{p\mathrm{b}^{-1}}$.The integrated production cross-sections in the kinematic range $3.5<p_{\mathrm{T}}<14~\mathrm{GeV}/c$ and $2.0<y<4.5$ are \begin{equation*} \begin{split} \sigma(\mbox{prompt }\psi(2S),7~\mathrm{TeV}) &={0.471 \pm0.001 (\mathrm{stat}) \pm 0.025 (\mathrm{syst})\mu\mathrm{b}},\\ \sigma(\psi(2S)\mathrm{-from-}b,7~\mathrm{TeV}) &={0.126\pm0.001 (\mathrm{stat}) \pm0.008 (\mathrm{syst})\mu\mathrm{b}}. \end{split} \end{equation*} All results show reasonable agreement with theoretical calculations.

Tables

Table 1. Table 1: Systematic uncertainties on the ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} cross-section measurements. The uncertainty from the t z subscript 𝑡 𝑧 t_{z} fit only affects the ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 {\psi{(2S)}} -from- b 𝑏 b result. Uncertainties labelled with “ ∗ ∗ \ast ” are correlated between kinematic bins.

Source	$13\text{ TeV}$ (%)	$7\text{ TeV}$ (%)
Signal mass shape∗	$0.0$–$4.1$	$0.0$–$8.5$
Radiative tail∗	$1.0$	$1.0$
Tracking∗	$(0.1$–$2.4)\oplus (2\times 0.8)$	$(0.7$–$3.0)\oplus (2\times 0.4)$
Muon ID∗	$(0.1$–$1.1)\oplus (0.1$–$4.6)$	$(0.7$–$8.9)\oplus (0.4$–$5.4)$
Trigger∗	$0.1$–$9.3$	$0.0$–$4.4$
Kinematic spectrum	$0.0$–$2.0$	$0.0$–$4.9$
Luminosity∗	$3.9$	$1.7$
$ℬ{\left(\psi \left(2S\right)\to e^{+}{e}^{-}\right)}^{\ast }$	$2.2$	$2.2$
Simulated sample size (prompt $\psi \left(2S\right)$)	$0.7$–$11.5$	$1.3$–$13.1$
Simulated sample size ($\psi \left(2S\right)$-from-$b$)	$0.8$–$5.7$	$1.2$–$9.5$
$t_z$ fit∗ ($\psi \left(2S\right)$-from-$b$ only)	$0.1$–$8.4$	$0.1$–$9.2$

Table 2. Table 2: Double-differential production cross-sections (in  nb /( GeV / c GeV 𝑐 {\mathrm{Ge\kern-1.00006ptV\!/}c} )) of prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} mesons at 13 ​  TeV 13  TeV 13\text{\,Te\kern-1.00006ptV} in bins of ( p T subscript 𝑝 T p_{\mathrm{T}} , y 𝑦 y ). The first uncertainties are statistical, the second are the uncorrelated systematic uncertainties between bins, and the last are the correlated systematic uncertainties between bins. Adjacent bins with large statistical uncertainties have been merged.

$p_T$ ($\mathrm{GeV}/c$)
2–3	$232.86\pm 5.80\pm 5.39\pm 18.24$	$228.46\pm 3.53\pm 1.65\pm 17.70$	$198.36\pm 2.86\pm 1.34\pm 15.33$	$167.07\pm 2.19\pm 1.30\pm 12.93$	$133.13\pm 1.99\pm 1.62\pm 10.42$
3–4	$192.93\pm 3.96\pm 4.02\pm 9.74$	$166.62\pm 2.36\pm 1.27\pm 8.22$	$147.39\pm 1.90\pm 1.07\pm 7.25$	$120.42\pm 1.49\pm 0.96\pm 6.00$	$85.03\pm 1.42\pm 1.09\pm 4.37$
4–5	$124.73\pm 2.38\pm 2.70\pm 9.15$	$108.56\pm 1.41\pm 0.88\pm 7.92$	$94.79\pm 1.09\pm 0.74\pm 6.90$	$76.73\pm 0.93\pm 0.68\pm 5.59$	$57.41\pm 1.00\pm 0.82\pm 4.28$
5–6	$77.27\pm 1.46\pm 1.78\pm 4.89$	$66.05\pm 0.82\pm 0.59\pm 4.16$	$58.12\pm 0.65\pm 0.52\pm 3.65$	$47.25\pm 0.44\pm 0.48\pm 2.98$	$34.35\pm 0.66\pm 0.58\pm 2.29$
6–7	$45.13\pm 0.91\pm 1.11\pm 4.73$	$41.43\pm 0.53\pm 0.43\pm 4.34$	$35.20\pm 0.43\pm 0.37\pm 3.69$	$28.70\pm 0.39\pm 0.36\pm 3.01$	$19.22\pm 0.47\pm 0.38\pm 2.08$
7–8	$28.87\pm 0.62\pm 0.77\pm 1.95$	$25.38\pm 0.35\pm 0.32\pm 1.70$	$21.00\pm 0.30\pm 0.27\pm 1.41$	$15.97\pm 0.27\pm 0.25\pm 1.07$	$12.23\pm 0.32\pm 0.31\pm 0.91$
8–9	$17.52\pm 0.42\pm 0.49\pm 1.19$	$15.39\pm 0.25\pm 0.24\pm 1.04$	$12.75\pm 0.22\pm 0.20\pm 0.86$	$9.86\pm 0.20\pm 0.19\pm 0.67$	$6.08\pm 0.21\pm 0.18\pm 0.47$
9–10	$10.94\pm 0.29\pm 0.34\pm 0.56$	$9.63\pm 0.18\pm 0.18\pm 0.48$	$7.46\pm 0.16\pm 0.14\pm 0.38$	$5.90\pm 0.15\pm 0.15\pm 0.30$	$3.81\pm 0.16\pm 0.14\pm 0.26$
10–11	$7.66\pm 0.23\pm 0.28\pm 0.50$	$5.98\pm 0.14\pm 0.13\pm 0.39$	$4.84\pm 0.13\pm 0.11\pm 0.31$	$3.83\pm 0.12\pm 0.11\pm 0.25$	$2.47\pm 0.11\pm 0.11\pm 0.19$
11–12	$4.25\pm 0.16\pm 0.17\pm 0.23$	$3.86\pm 0.11\pm 0.10\pm 0.21$	$3.16\pm 0.10\pm 0.09\pm 0.17$	$2.37\pm 0.09\pm 0.08\pm 0.13$	$1.67\pm 0.10\pm 0.09\pm 0.13$
12–13	$3.09\pm 0.13\pm 0.16\pm 0.27$	$2.44\pm 0.13\pm 0.07\pm 0.21$	$2.02\pm 0.08\pm 0.07\pm 0.18$	$1.39\pm 0.07\pm 0.06\pm 0.12$	$0.74\pm 0.06\pm 0.05\pm 0.07$
13–14	$1.63\pm 0.09\pm 0.08\pm 0.12$	$1.57\pm 0.07\pm 0.06\pm 0.12$	$1.26\pm 0.06\pm 0.05\pm 0.09$	$0.91\pm 0.05\pm 0.04\pm 0.07$	$0.58\pm 0.06\pm 0.04\pm 0.05$
14–15	$1.42\pm 0.08\pm 0.08\pm 0.07$	$1.26\pm 0.06\pm 0.05\pm 0.06$	$0.81\pm 0.05\pm 0.04\pm 0.04$	$0.59\pm 0.04\pm 0.04\pm 0.03$	$0.26\pm 0.03\pm 0.02\pm 0.02$
15–16	$1.18\pm 0.08\pm 0.09\pm 0.06$	$0.86\pm 0.05\pm 0.04\pm 0.04$	$0.58\pm 0.04\pm 0.03\pm 0.03$	$0.41\pm 0.03\pm 0.02\pm 0.02$
16–17	$0.81\pm 0.06\pm 0.07\pm 0.04$	$0.59\pm 0.04\pm 0.03\pm 0.03$	$0.39\pm 0.03\pm 0.03\pm 0.02$		$0.15\pm 0.01\pm 0.02\pm 0.01$
17–18	$0.54\pm 0.05\pm 0.06\pm 0.03$	$0.42\pm 0.03\pm 0.03\pm 0.02$	$0.21\pm 0.01\pm 0.01\pm 0.01$	$0.15\pm 0.01\pm 0.01\pm 0.01$
18–19	$0.29\pm 0.02\pm 0.02\pm 0.01$	$0.28\pm 0.02\pm 0.02\pm 0.02$
19–20

Table 3. Table 3: Double-differential production cross-sections (in  nb / ( GeV / c )  nb GeV 𝑐 \text{\,nb}/({\mathrm{Ge\kern-1.00006ptV\!/}c}) ) of ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} -from- b 𝑏 b mesons at 13 ​  TeV 13  TeV 13\text{\,Te\kern-1.00006ptV} in bins of ( p T subscript 𝑝 T p_{\mathrm{T}} , y 𝑦 y ). The first uncertainties are statistical, the second are the uncorrelated systematic uncertainties between bins, and the last are the correlated systematic uncertainties between bins. Adjacent bins with large statistical uncertainties have been merged.

$p_T$ ($\mathrm{GeV}/c$)
2–3	$69.99\pm 1.94\pm 1.65\pm 5.51$	$61.15\pm 1.14\pm 0.48\pm 4.78$	$51.47\pm 0.97\pm 0.40\pm 3.99$	$37.86\pm 0.87\pm 0.35\pm 2.95$	$24.04\pm 0.94\pm 0.36\pm 1.89$
3–4	$56.37\pm 1.45\pm 1.17\pm 3.05$	$47.89\pm 0.83\pm 0.36\pm 2.45$	$39.58\pm 0.70\pm 0.31\pm 2.45$	$29.92\pm 0.64\pm 0.27\pm 1.57$	$19.49\pm 0.70\pm 0.29\pm 1.01$
4–5	$37.99\pm 1.00\pm 0.78\pm 2.80$	$34.38\pm 0.58\pm 0.25\pm 2.51$	$28.42\pm 0.48\pm 0.21\pm 2.08$	$21.25\pm 0.45\pm 0.19\pm 1.55$	$12.89\pm 0.50\pm 0.20\pm 0.97$
5–6	$27.93\pm 0.72\pm 0.60\pm 1.77$	$23.18\pm 0.40\pm 0.18\pm 1.46$	$19.18\pm 0.34\pm 0.15\pm 1.24$	$14.62\pm 0.24\pm 0.14\pm 0.94$	$7.94\pm 0.35\pm 0.14\pm 0.56$
6–7	$16.99\pm 0.49\pm 0.37\pm 1.78$	$15.74\pm 0.29\pm 0.13\pm 1.65$	$12.33\pm 0.24\pm 0.11\pm 1.29$	$8.91\pm 0.22\pm 0.10\pm 0.94$	$4.83\pm 0.25\pm 0.09\pm 0.52$
7–8	$12.29\pm 0.37\pm 0.28\pm 0.83$	$10.14\pm 0.21\pm 0.09\pm 0.68$	$7.77\pm 0.18\pm 0.08\pm 0.53$	$5.75\pm 0.17\pm 0.07\pm 0.39$	$3.41\pm 0.18\pm 0.08\pm 0.26$
8–9	$9.05\pm 0.29\pm 0.20\pm 0.62$	$6.84\pm 0.16\pm 0.07\pm 0.46$	$5.24\pm 0.14\pm 0.06\pm 0.36$	$3.56\pm 0.13\pm 0.05\pm 0.24$	$2.09\pm 0.13\pm 0.05\pm 0.16$
9–10	$6.12\pm 0.22\pm 0.13\pm 0.32$	$4.85\pm 0.13\pm 0.06\pm 0.24$	$3.60\pm 0.11\pm 0.05\pm 0.19$	$2.53\pm 0.10\pm 0.04\pm 0.13$	$1.50\pm 0.10\pm 0.05\pm 0.10$
10–11	$4.17\pm 0.16\pm 0.11\pm 0.27$	$3.38\pm 0.11\pm 0.05\pm 0.22$	$2.52\pm 0.09\pm 0.04\pm 0.16$	$1.63\pm 0.08\pm 0.03\pm 0.11$	$0.87\pm 0.06\pm 0.03\pm 0.07$
11–12	$2.90\pm 0.13\pm 0.07\pm 0.16$	$2.38\pm 0.09\pm 0.04\pm 0.13$	$1.66\pm 0.08\pm 0.03\pm 0.11$	$1.07\pm 0.06\pm 0.03\pm 0.06$	$0.65\pm 0.06\pm 0.03\pm 0.06$
12–13	$2.09\pm 0.11\pm 0.06\pm 0.19$	$1.68\pm 0.07\pm 0.03\pm 0.15$	$1.22\pm 0.06\pm 0.03\pm 0.11$	$0.78\pm 0.05\pm 0.02\pm 0.07$	$0.41\pm 0.05\pm 0.02\pm 0.04$
13–14	$1.39\pm 0.08\pm 0.04\pm 0.11$	$1.24\pm 0.06\pm 0.03\pm 0.09$	$0.76\pm 0.05\pm 0.02\pm 0.06$	$0.61\pm 0.04\pm 0.02\pm 0.05$	$0.29\pm 0.04\pm 0.02\pm 0.03$
14–15	$1.18\pm 0.07\pm 0.04\pm 0.06$	$0.82\pm 0.05\pm 0.02\pm 0.04$	$0.71\pm 0.05\pm 0.02\pm 0.04$	$0.42\pm 0.04\pm 0.02\pm 0.02$	$0.16\pm 0.02\pm 0.01\pm 0.02$
15–16	$0.84\pm 0.06\pm 0.03\pm 0.04$	$0.70\pm 0.05\pm 0.02\pm 0.04$	$0.42\pm 0.03\pm 0.01\pm 0.03$	$0.24\pm 0.02\pm 0.01\pm 0.02$
16–17	$0.62\pm 0.05\pm 0.02\pm 0.03$	$0.54\pm 0.04\pm 0.02\pm 0.03$	$0.34\pm 0.03\pm 0.01\pm 0.02$		$0.07\pm 0.01\pm 0.00\pm 0.01$
17–18	$0.63\pm 0.05\pm 0.03\pm 0.04$	$0.41\pm 0.03\pm 0.02\pm 0.02$	$0.21\pm 0.01\pm 0.01\pm 0.01$	$0.09\pm 0.01\pm 0.00\pm 0.01$
18–19	$0.32\pm 0.03\pm 0.01\pm 0.02$	$0.28\pm 0.02\pm 0.01\pm 0.02$
19–20

Table 4. Table 4: Double-differential production cross-section in  nb /( GeV / c GeV 𝑐 {\mathrm{Ge\kern-1.00006ptV\!/}c} ) of prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} mesons at 7 ​  TeV 7  TeV 7\text{\,Te\kern-1.00006ptV} in bins of ( p T subscript 𝑝 T p_{\mathrm{T}} , y 𝑦 y ). The first uncertainty is statistical, the second is the uncorrelated systematic uncertainties shared between bins and the last is the correlated systematic uncertainties.

$p_T$ ($\text{ GeV/}c$)
3.5–4	$94.34\pm 2.49\pm 5.59\pm 8.65$	$109.85\pm 1.38\pm 2.50\pm 5.04$	$97.07\pm 0.99\pm 1.90\pm 4.33$	$77.87\pm 0.78\pm 1.51\pm 3.11$	$57.90\pm 0.80\pm 1.57\pm 2.09$
4–5	$73.08\pm 1.21\pm 4.35\pm 4.40$	$75.42\pm 0.66\pm 1.32\pm 3.53$	$65.66\pm 0.47\pm 0.85\pm 2.64$	$53.33\pm 0.40\pm 0.88\pm 2.23$	$36.34\pm 0.40\pm 0.80\pm 1.30$
5–6	$44.86\pm 0.71\pm 2.47\pm 2.24$	$43.26\pm 0.38\pm 1.04\pm 4.14$	$37.21\pm 0.30\pm 0.58\pm 1.54$	$31.96\pm 0.27\pm 0.86\pm 1.31$	$18.98\pm 0.26\pm 0.93\pm 0.65$
6–7	$25.99\pm 0.45\pm 1.64\pm 2.58$	$24.94\pm 0.25\pm 0.63\pm 1.20$	$22.42\pm 0.21\pm 0.58\pm 0.86$	$17.86\pm 0.19\pm 0.61\pm 0.69$	$10.33\pm 0.17\pm 0.39\pm 0.50$
7–8	$16.58\pm 0.30\pm 1.26\pm 1.04$	$15.21\pm 0.17\pm 0.41\pm 0.66$	$12.82\pm 0.14\pm 0.40\pm 0.50$	$9.98\pm 0.13\pm 0.33\pm 0.44$	$5.61\pm 0.12\pm 0.30\pm 0.26$
8–9	$7.03\pm 0.14\pm 1.16\pm 0.35$	$9.02\pm 0.12\pm 0.39\pm 0.41$	$7.16\pm 0.10\pm 0.31\pm 0.33$	$5.05\pm 0.09\pm 0.28\pm 0.23$	$2.92\pm 0.08\pm 0.18\pm 0.15$
9–10	$6.10\pm 0.14\pm 0.63\pm 0.28$	$5.26\pm 0.09\pm 0.26\pm 0.23$	$4.63\pm 0.08\pm 0.20\pm 0.18$	$3.16\pm 0.07\pm 0.20\pm 0.13$	$1.89\pm 0.07\pm 0.20\pm 0.12$
10–11	$3.89\pm 0.11\pm 0.36\pm 0.23$	$3.03\pm 0.06\pm 0.17\pm 0.14$	$2.50\pm 0.05\pm 0.17\pm 0.11$	$1.93\pm 0.06\pm 0.18\pm 0.07$	$1.04\pm 0.05\pm 0.15\pm 0.06$
11–12	$2.56\pm 0.08\pm 0.29\pm 0.13$	$2.30\pm 0.06\pm 0.35\pm 0.10$	$1.50\pm 0.04\pm 0.16\pm 0.13$	$1.46\pm 0.05\pm 0.25\pm 0.10$	$0.38\pm 0.02\pm 0.08\pm 0.02$
12–13	$2.42\pm 0.10\pm 0.90\pm 0.11$	$1.12\pm 0.04\pm 0.09\pm 0.05$	$0.84\pm 0.03\pm 0.12\pm 0.04$	$0.58\pm 0.03\pm 0.08\pm 0.05$	$0.35\pm 0.03\pm 0.10\pm 0.04$
13–14	$0.85\pm 0.04\pm 0.16\pm 0.04$	$0.88\pm 0.03\pm 0.12\pm 0.11$	$0.60\pm 0.03\pm 0.09\pm 0.04$	$0.33\pm 0.02\pm 0.06\pm 0.02$	$0.07\pm 0.01\pm 0.02\pm 0.01$

Table 5. Table 5: Double-differential production cross-section in  nb /( GeV / c GeV 𝑐 {\mathrm{Ge\kern-1.00006ptV\!/}c} ) of ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} -from- b 𝑏 b mesons at 7 ​  TeV 7  TeV 7\text{\,Te\kern-1.00006ptV} in bins of ( p T subscript 𝑝 T p_{\mathrm{T}} , y 𝑦 y ). The first uncertainties are statistical, the second are the uncorrelated systematic uncertainties shared between bins and the last are the correlated systematic uncertainties.

$p_T$ ($\text{ GeV/}c$)
3.5–4	$27.65\pm 1.10\pm 1.64\pm 2.86$	$24.88\pm 0.55\pm 0.44\pm 1.66$	$21.06\pm 0.42\pm 0.35\pm 0.98$	$16.01\pm 0.36\pm 0.43\pm 0.65$	$9.62\pm 0.36\pm 0.32\pm 0.53$
4–5	$20.85\pm 0.54\pm 1.12\pm 1.36$	$19.60\pm 0.29\pm 0.26\pm 1.14$	$15.58\pm 0.22\pm 0.19\pm 0.74$	$11.05\pm 0.19\pm 0.18\pm 0.46$	$6.58\pm 0.18\pm 0.18\pm 0.23$
5–6	$13.31\pm 0.36\pm 0.96\pm 0.93$	$13.04\pm 0.19\pm 0.19\pm 1.31$	$9.80\pm 0.15\pm 0.13\pm 0.43$	$6.90\pm 0.13\pm 0.12\pm 0.28$	$4.00\pm 0.13\pm 0.14\pm 0.14$
6–7	$8.77\pm 0.24\pm 0.45\pm 0.93$	$8.15\pm 0.14\pm 0.12\pm 0.40$	$6.20\pm 0.11\pm 0.09\pm 0.26$	$4.53\pm 0.10\pm 0.08\pm 0.18$	$2.01\pm 0.09\pm 0.08\pm 0.14$
7–8	$6.24\pm 0.18\pm 0.33\pm 0.41$	$4.96\pm 0.10\pm 0.08\pm 0.25$	$4.02\pm 0.08\pm 0.07\pm 0.16$	$2.64\pm 0.07\pm 0.05\pm 0.15$	$1.38\pm 0.07\pm 0.06\pm 0.07$
8–9	$3.68\pm 0.12\pm 0.19\pm 0.23$	$3.40\pm 0.08\pm 0.07\pm 0.16$	$2.69\pm 0.07\pm 0.05\pm 0.14$	$1.61\pm 0.05\pm 0.04\pm 0.08$	$0.76\pm 0.05\pm 0.04\pm 0.04$
9–10	$3.14\pm 0.11\pm 0.21\pm 0.19$	$2.15\pm 0.06\pm 0.04\pm 0.09$	$1.33\pm 0.05\pm 0.04\pm 0.12$	$1.00\pm 0.04\pm 0.03\pm 0.04$	$0.37\pm 0.03\pm 0.02\pm 0.04$
10–11	$1.96\pm 0.08\pm 0.12\pm 0.12$	$1.46\pm 0.05\pm 0.04\pm 0.07$	$1.11\pm 0.04\pm 0.04\pm 0.07$	$0.66\pm 0.03\pm 0.02\pm 0.03$	$0.32\pm 0.03\pm 0.02\pm 0.02$
11–12	$1.28\pm 0.06\pm 0.07\pm 0.07$	$1.00\pm 0.04\pm 0.03\pm 0.05$	$0.75\pm 0.03\pm 0.03\pm 0.07$	$0.42\pm 0.02\pm 0.02\pm 0.03$	$0.18\pm 0.02\pm 0.02\pm 0.01$
12–13	$0.80\pm 0.04\pm 0.05\pm 0.05$	$0.73\pm 0.03\pm 0.02\pm 0.04$	$0.51\pm 0.03\pm 0.02\pm 0.02$	$0.28\pm 0.02\pm 0.01\pm 0.03$	$0.09\pm 0.02\pm 0.01\pm 0.01$
13–14	$0.66\pm 0.04\pm 0.07\pm 0.04$	$0.51\pm 0.03\pm 0.02\pm 0.06$	$0.33\pm 0.02\pm 0.01\pm 0.02$	$0.17\pm 0.02\pm 0.01\pm 0.01$	$0.06\pm 0.01\pm 0.01\pm 0.01$

Table 6. Table 6: Differential production cross-sections d ​ σ / d ​ p T d 𝜎 d subscript 𝑝 T \mathrm{d}\sigma/\mathrm{d}p_{\mathrm{T}} (in  nb /( GeV / c GeV 𝑐 {\mathrm{Ge\kern-1.00006ptV\!/}c} )) of prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} and ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} -from- b 𝑏 {b} mesons at 13 ​  TeV 13  TeV 13\text{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

$p_T$ ($\mathrm{GeV}/c$)	Prompt $\psi \left(2S\right)$	$\psi \left(2S\right)$-from-$b$
2–3	$479.94\pm 3.97\pm 3.08\pm 37.31$	$122.26\pm 1.38\pm 0.91\pm 9.56$
3–4	$356.19\pm 2.69\pm 2.29\pm 17.79$	$96.62\pm 1.02\pm 0.66\pm 5.27$
4–5	$231.12\pm 1.64\pm 1.56\pm 16.92$	$67.47\pm 0.71\pm 0.45\pm 4.95$
5–6	$141.52\pm 0.98\pm 1.04\pm 8.99$	$46.43\pm 0.49\pm 0.34\pm 2.98$
6–7	$84.84\pm 0.65\pm 0.68\pm 8.92$	$29.40\pm 0.35\pm 0.21\pm 3.10$
7–8	$51.72\pm 0.44\pm 0.48\pm 3.52$	$19.68\pm 0.26\pm 0.16\pm 1.34$
8–9	$30.80\pm 0.30\pm 0.32\pm 2.11$	$13.38\pm 0.20\pm 0.12\pm 0.92$
9–10	$18.87\pm 0.22\pm 0.23\pm 0.99$	$9.29\pm 0.16\pm 0.08\pm 0.49$
10–11	$12.39\pm 0.17\pm 0.18\pm 0.82$	$6.29\pm 0.12\pm 0.07\pm 0.42$
11–12	$7.65\pm 0.13\pm 0.13\pm 0.44$	$4.33\pm 0.10\pm 0.05\pm 0.26$
12–13	$4.84\pm 0.11\pm 0.10\pm 0.43$	$3.09\pm 0.08\pm 0.04\pm 0.28$
13–14	$2.98\pm 0.07\pm 0.06\pm 0.23$	$2.14\pm 0.06\pm 0.03\pm 0.17$
14–20	$1.13\pm 0.02\pm 0.02\pm 0.06$	$0.90\pm 0.02\pm 0.01\pm 0.05$

Table 7. Table 7: Differential production cross-sections d ​ σ / d ​ y d 𝜎 d 𝑦 \mathrm{d}\sigma/\mathrm{d}y (in  nb ) of prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} and ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} -from- b 𝑏 {b} mesons at 13 ​  TeV 13  TeV 13\text{\,Te\kern-1.00006ptV} per rapidity unit. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

$y$	Prompt $\psi \left(2S\right)$	$\psi \left(2S\right)$-from-$b$
2.0–2.5	$751.4\pm 7.7\pm 7.6\pm 51.8$	$251.2\pm 2.8\pm 2.3\pm 17.6$
2.5–3.0	$679.1\pm 4.6\pm 2.4\pm 46.7$	$215.9\pm 1.6\pm 0.7\pm 15.0$
3.0–3.5	$588.7\pm 3.7\pm 2.0\pm 40.4$	$175.8\pm 1.4\pm 0.6\pm 12.7$
3.5–4.0	$482.3\pm 2.9\pm 1.9\pm 33.2$	$129.6\pm 1.2\pm 0.5\pm 9.1$
4.0–4.5	$357.8\pm 2.8\pm 2.3\pm 25.6$	$79.0\pm 1.4\pm 0.5\pm 5.7$

Table 8. Table 8: Differential cross-sections d ​ σ / d ​ p T d 𝜎 d subscript 𝑝 T \mathrm{d}\sigma/\mathrm{d}p_{\mathrm{T}} (in  nb /( GeV / c GeV 𝑐 {\mathrm{Ge\kern-1.00006ptV\!/}c} )) of prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} and ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} -from- b 𝑏 b mesons at 7 ​  TeV 7  TeV 7\text{\,Te\kern-1.00006ptV} , integrated over y 𝑦 y between 2.0 and 4.5. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

$p_T$ ($\text{ GeV/}c$)	Prompt $\psi \left(2S\right)$	$\psi \left(2S\right)$-from-$b$
3.5-4	$218.52\pm 1.61\pm 3.38\pm 11.61$	$49.61\pm 0.70\pm 0.91\pm 3.34$
4–5	$151.91\pm 0.78\pm 2.39\pm 7.05$	$36.83\pm 0.35\pm 0.60\pm 1.97$
5–6	$88.13\pm 0.47\pm 1.51\pm 4.94$	$23.52\pm 0.24\pm 0.50\pm 1.55$
6–7	$50.77\pm 0.30\pm 0.99\pm 2.92$	$14.83\pm 0.16\pm 0.24\pm 0.95$
7–8	$30.10\pm 0.21\pm 0.73\pm 1.45$	$9.62\pm 0.12\pm 0.18\pm 0.52$
8–9	$15.59\pm 0.12\pm 0.65\pm 0.74$	$6.08\pm 0.09\pm 0.11\pm 0.32$
9–10	$10.52\pm 0.10\pm 0.38\pm 0.47$	$3.99\pm 0.07\pm 0.11\pm 0.24$
10–11	$6.20\pm 0.08\pm 0.25\pm 0.30$	$2.75\pm 0.05\pm 0.07\pm 0.16$
11–12	$4.10\pm 0.06\pm 0.27\pm 0.24$	$1.81\pm 0.04\pm 0.04\pm 0.12$
12–13	$2.65\pm 0.06\pm 0.46\pm 0.15$	$1.20\pm 0.03\pm 0.03\pm 0.08$
13–14	$1.37\pm 0.03\pm 0.12\pm 0.11$	$0.86\pm 0.03\pm 0.04\pm 0.07$

Table 9. Table 9: Differential cross-sections d ​ σ / d ​ y d 𝜎 d 𝑦 \mathrm{d}\sigma/\mathrm{d}y (in  nb ) of prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} and ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} -from- b 𝑏 b mesons at 7 ​  TeV 7  TeV 7\text{\,Te\kern-1.00006ptV} , integrated over p T subscript 𝑝 T p_{\mathrm{T}} between 3.5 and 14  GeV​/ ​ c  GeV​/ 𝑐 \text{\,Ge\kern-1.00006ptV\!/}c . The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

$y$	Prompt $\psi \left(2S\right)$	$\psi \left(2S\right)$-from-$b$
2.0–2.5	$230.5\pm 2.0\pm 6.3\pm 15.7$	$74.5\pm 0.9\pm 1.8\pm 5.8$
2.5–3.0	$235.4\pm 1.1\pm 2.3\pm 13.1$	$67.4\pm 0.5\pm 0.4\pm 4.4$
3.0–3.5	$203.9\pm 0.8\pm 1.6\pm 8.5$	$52.8\pm 0.4\pm 0.3\pm 2.5$
3.5–4.0	$164.6\pm 0.7\pm 1.7\pm 6.8$	$37.3\pm 0.3\pm 0.3\pm 1.6$
4.0–4.5	$106.9\pm 0.7\pm 1.6\pm 4.2$	$20.6\pm 0.3\pm 0.3\pm 1.0$

Table 10. Table 10: Fractions of ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} -from- b 𝑏 {b} (in %) at 13 ​  TeV 13  TeV 13\text{\,Te\kern-1.00006ptV} in bins of ( p T , y ) subscript 𝑝 T 𝑦 (p_{\mathrm{T}},y) of ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} mesons. The uncertainties are statistical only. The systematic uncertainties are negligible. Adjacent bins with large statistical uncertainty have been merged.

$p_T$ ($\mathrm{GeV}/c$)
2–3	$23.0\pm 0.6$	$20.8\pm 0.4$	$20.4\pm 0.4$	$18.1\pm 0.4$	$14.8\pm 0.5$
3–4	$22.4\pm 0.6$	$22.2\pm 0.4$	$21.0\pm 0.4$	$19.5\pm 0.4$	$18.4\pm 0.6$
4–5	$23.0\pm 0.6$	$23.9\pm 0.4$	$22.8\pm 0.4$	$21.4\pm 0.4$	$18.8\pm 0.7$
5–6	$25.9\pm 0.6$	$25.5\pm 0.4$	$24.4\pm 0.4$	$22.8\pm 0.5$	$18.5\pm 0.7$
6–7	$27.1\pm 0.7$	$27.1\pm 0.5$	$25.6\pm 0.5$	$23.3\pm 0.5$	$20.2\pm 1.0$
7–8	$29.8\pm 0.8$	$28.7\pm 0.5$	$26.4\pm 0.6$	$25.8\pm 0.7$	$22.1\pm 1.0$
8–9	$32.8\pm 0.9$	$30.6\pm 0.6$	$28.8\pm 0.7$	$26.6\pm 0.8$	$25.3\pm 1.4$
9–10	$34.2\pm 1.0$	$32.8\pm 0.8$	$32.4\pm 0.9$	$29.7\pm 1.0$	$28.7\pm 1.7$
10–11	$35.5\pm 1.2$	$35.4\pm 0.9$	$33.7\pm 1.1$	$29.7\pm 1.2$	$26.9\pm 1.8$
11–12	$38.5\pm 1.5$	$37.6\pm 1.1$	$34.8\pm 1.4$	$30.6\pm 1.6$	$28.1\pm 2.4$
12–13	$40.7\pm 1.7$	$38.7\pm 1.4$	$36.4\pm 1.5$	$34.3\pm 2.0$	$31.7\pm 3.4$
13–14	$42.1\pm 2.1$	$41.8\pm 1.6$	$35.9\pm 1.9$	$37.6\pm 1.0$	$31.2\pm 3.9$
14–15	$44.2\pm 2.1$	$39.4\pm 1.8$	$44.7\pm 2.2$	$36.4\pm 2.8$	$39.2\pm 1.7$
15–16	$43.4\pm 2.5$	$42.5\pm 2.1$	$41.5\pm 1.0$	$36.8\pm 2.6$
16–17	$43.6\pm 2.8$	$45.9\pm 2.5$	$44.8\pm 1.3$		$29.0\pm 3.9$
17–18	$51.7\pm 3.1$	$47.1\pm 3.0$	$47.2\pm 2.5$	$34.7\pm 3.6$
18–19	$49.7\pm 3.2$	$49.1\pm 2.5$
19–20

Table 11. Table 11: Fractions of ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} -from- b 𝑏 {b} (in %) at 7 ​  TeV 7  TeV 7\text{\,Te\kern-1.00006ptV} in bins of ( p T , y ) subscript 𝑝 T 𝑦 (p_{\mathrm{T}},y) of ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} mesons. The uncertainties are statistical only. The systematic uncertainties are negligible.

$p_T$ ($\mathrm{GeV}/c$)
3.5–4	$21.11\pm 0.79$	$17.91\pm 0.37$	$17.26\pm 0.33$	$16.52\pm 0.34$	$13.93\pm 0.50$
4–5	$21.73\pm 0.53$	$20.21\pm 0.28$	$18.42\pm 0.25$	$16.95\pm 0.27$	$15.47\pm 0.40$
5–6	$22.15\pm 0.55$	$22.40\pm 0.31$	$20.08\pm 0.29$	$17.89\pm 0.32$	$16.49\pm 0.50$
6–7	$24.12\pm 0.61$	$24.49\pm 0.37$	$21.45\pm 0.35$	$20.21\pm 0.41$	$16.18\pm 0.65$
7–8	$27.34\pm 0.68$	$25.04\pm 0.44$	$23.81\pm 0.44$	$21.64\pm 0.52$	$18.48\pm 0.83$
8–9	$28.44\pm 0.80$	$27.26\pm 0.54$	$26.81\pm 0.57$	$23.68\pm 0.70$	$19.12\pm 1.10$
9–10	$30.46\pm 0.08$	$29.40\pm 0.68$	$21.50\pm 0.69$	$24.95\pm 0.88$	$17.28\pm 1.38$
10–11	$32.97\pm 1.08$	$31.67\pm 0.83$	$28.14\pm 0.87$	$27.30\pm 1.17$	$24.10\pm 1.77$
11–12	$34.86\pm 1.33$	$33.31\pm 1.02$	$30.98\pm 1.16$	$28.55\pm 1.44$	$27.33\pm 2.10$
12–13	$35.56\pm 1.62$	$37.85\pm 1.31$	$33.30\pm 1.44$	$28.50\pm 1.78$	$20.81\pm 3.01$
13–14	$38.79\pm 1.93$	$36.68\pm 1.51$	$34.21\pm 1.80$	$30.97\pm 2.50$	$36.20\pm 6.19$

Table 12. Table 12: Ratios of production cross-sections at 13 ​  TeV 13  TeV 13\text{\,Te\kern-1.00006ptV} between ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} mesons and J / ψ 𝐽 𝜓 {J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu} mesons in bins of p T subscript 𝑝 T p_{\mathrm{T}} for prompt production and for those from b 𝑏 {b} -hadron decays integrated in the rapidity range 2.0 < y < 4.5 2.0 𝑦 4.5 2.0<y<4.5 . The statistical and systematic uncertainties are added in quadrature.

$p_T$ ($\mathrm{GeV}/c$)	Prompt	From $b$-hadron decays
2–3	$0.135\pm 0.010$	$0.241\pm 0.018$
3–4	$0.158\pm 0.007$	$0.263\pm 0.013$
4–5	$0.182\pm 0.012$	$0.278\pm 0.019$
5–6	$0.201\pm 0.011$	$0.317\pm 0.018$
6–7	$0.225\pm 0.023$	$0.311\pm 0.031$
7–8	$0.259\pm 0.016$	$0.341\pm 0.021$
8–9	$0.271\pm 0.017$	$0.374\pm 0.024$
9–10	$0.284\pm 0.012$	$0.365\pm 0.018$
10–11	$0.310\pm 0.019$	$0.370\pm 0.025$
11–12	$0.305\pm 0.017$	$0.357\pm 0.023$
12–13	$0.315\pm 0.029$	$0.380\pm 0.037$
13–14	$0.294\pm 0.026$	$0.361\pm 0.033$

Table 13. Table 13: Ratios of production cross-sections between ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} mesons and J / ψ 𝐽 𝜓 {J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu} mesons at 13 ​  TeV 13  TeV 13\text{\,Te\kern-1.00006ptV} in bins of y 𝑦 y for prompt production and for those from b 𝑏 {b} -hadron decays integrated in the transverse momentum range 2 < p T < 14 ​  GeV​/ ​ c 2 subscript 𝑝 T 14  GeV​/ 𝑐 2<p_{\mathrm{T}}<14\text{\,Ge\kern-1.00006ptV\!/}c . The statistical and systematic uncertainties are added in quadrature.

$y$	Prompt	From $b$-hadron decays
2.0–2.5	$0.175\pm 0.011$	$0.288\pm 0.018$
2.5–3.0	$0.168\pm 0.010$	$0.289\pm 0.017$
3.0–3.5	$0.164\pm 0.009$	$0.279\pm 0.017$
3.5–4.0	$0.159\pm 0.010$	$0.261\pm 0.017$
4.0–4.5	$0.150\pm 0.012$	$0.237\pm 0.019$

Table 14. Table 14: Ratios of production cross-sections between 13 ​  TeV 13  TeV 13\text{\,Te\kern-1.00006ptV} and 7 ​  TeV 7  TeV 7\text{\,Te\kern-1.00006ptV} in bins of p T subscript 𝑝 T p_{\mathrm{T}} for prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} and ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} -from- b 𝑏 {b} mesons integrated in the rapidity range 2.0 < y < 4.5 2.0 𝑦 4.5 2.0<y<4.5 . The statistical and systematic uncertainties are added in quadrature.

$p_T$ ($\mathrm{GeV}/c$)	Prompt $\psi \left(2S\right)$	$\psi \left(2S\right)$-from-$b$
3.5–4	$1.63\pm 0.09$	$1.95\pm 0.13$
4–5	$1.52\pm 0.12$	$1.83\pm 0.15$
5–6	$1.61\pm 0.11$	$1.97\pm 0.15$
6–7	$1.67\pm 0.19$	$1.98\pm 0.23$
7–8	$1.72\pm 0.12$	$2.05\pm 0.14$
8–9	$1.97\pm 0.13$	$2.20\pm 0.15$
9–10	$1.79\pm 0.08$	$2.33\pm 0.14$
10–11	$2.00\pm 0.13$	$2.28\pm 0.17$
11–12	$1.87\pm 0.12$	$2.39\pm 0.15$
12–13	$1.83\pm 0.20$	$2.57\pm 0.24$
13–14	$2.18\pm 0.25$	$2.49\pm 0.27$

Table 15. Table 15: Ratios of cross-sections between measurements at 13 ​  TeV 13  TeV 13\text{\,Te\kern-1.00006ptV} and 7 ​  TeV 7  TeV 7\text{\,Te\kern-1.00006ptV} in different bins of y 𝑦 y for prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} and ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} -from- b 𝑏 {b} mesons integrated in the transverse momentum range 3.5 < p T < 14 ​  GeV​/ ​ c 3.5 subscript 𝑝 T 14  GeV​/ 𝑐 3.5<p_{\mathrm{T}}<14\text{\,Ge\kern-1.00006ptV\!/}c . The statistical and systematic uncertainties are added in quadrature.

$y$	Prompt $\psi \left(2S\right)$	$\psi \left(2S\right)$-from-$b$
2.0–2.5	$1.81\pm 0.14$	$2.00\pm 0.17$
2.5–3.0	$1.54\pm 0.11$	$1.89\pm 0.14$
3.0–3.5	$1.54\pm 0.10$	$1.94\pm 0.13$
3.5–4.0	$1.54\pm 0.10$	$2.03\pm 0.13$
4.0–4.5	$1.69\pm 0.10$	$2.17\pm 0.14$

Table 16. Table 16: Multiplicative scaling factors needed to obtain the prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} differential cross-sections from unpolarised cross-section measurements at 13 ​  TeV 13  TeV 13\text{\,Te\kern-1.00006ptV} as reported in Table 2 under the assumption of fully transverse polarisation ( α = + 1 𝛼 1 \alpha=+1 ).

$p_T$ ($\mathrm{GeV}/c$)
2–3	$1.084$	$1.081$	$0.997$	$0.984$	$0.993$
3–4	$1.202$	$1.047$	$0.964$	$0.950$	$0.948$
4–5	$1.164$	$1.034$	$0.961$	$0.937$	$0.921$
5–6	$1.138$	$1.016$	$0.959$	$0.940$	$0.922$
6–7	$1.100$	$1.012$	$0.970$	$0.939$	$0.923$
7–8	$1.110$	$0.996$	$0.965$	$0.937$	$0.915$
8–9	$1.120$	$0.993$	$0.962$	$0.943$	$0.937$
9–10	$1.037$	$0.985$	$0.962$	$0.940$	$0.941$
10–11	$1.070$	$0.941$	$0.941$	$0.960$	$0.910$
11–12	$1.022$	$0.952$	$0.958$	$0.926$	$0.925$
12–13	$1.096$	$0.940$	$0.932$	$0.944$	$0.897$
13–14	$1.007$	$0.964$	$0.949$	$0.932$	$0.924$
14–15	$1.060$	$0.944$	$0.939$	$0.924$	$0.909$
15–16	$0.932$	$0.940$	$0.910$	$0.931$
16–17	$1.051$	$0.932$	$0.930$		$0.850$
17–18	$1.080$	$0.946$	$0.924$	$0.911$
18–19	$0.969$	$0.934$
19–20

Table 17. Table 17: Multiplicative scaling factors needed to obtain the prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} differential cross-sections from unpolarised cross-section measurements at 13 ​  TeV 13  TeV 13\text{\,Te\kern-1.00006ptV} as reported in Table 2 under the assumption of fully longitudinal polarisation ( α = − 1 𝛼 1 \alpha=-1 ).

$p_T$ ($\mathrm{GeV}/c$)
2–3	$0.846$	$0.870$	$0.987$	$1.016$	$1.008$
3–4	$0.754$	$0.909$	$1.039$	$1.079$	$1.092$
4–5	$0.777$	$0.934$	$1.053$	$1.092$	$1.123$
5–6	$0.806$	$0.950$	$1.051$	$1.091$	$1.132$
6–7	$0.829$	$0.958$	$1.040$	$1.069$	$1.129$
7–8	$0.843$	$0.977$	$1.041$	$1.079$	$1.102$
8–9	$0.865$	$0.992$	$1.054$	$1.073$	$1.120$
9–10	$0.908$	$1.004$	$1.061$	$1.066$	$1.129$
10–11	$0.882$	$1.024$	$1.068$	$1.098$	$1.096$
11–12	$0.944$	$1.035$	$1.077$	$1.069$	$1.114$
12–13	$0.941$	$1.050$	$1.079$	$1.076$	$1.126$
13–14	$0.992$	$1.069$	$1.067$	$1.100$	$1.115$
14–15	$0.966$	$1.076$	$1.085$	$1.114$	$1.143$
15–16	$0.964$	$1.082$	$1.142$	$1.102$
16–17	$1.011$	$1.099$	$1.102$		$1.013$
17–18	$0.992$	$1.074$	$1.114$	$1.140$
18–19	$1.037$	$1.094$
19–20

Table 18. Table 18: Multiplicative scaling factors needed to obtain the prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} differential cross-sections from unpolarised cross-section measurements at 13 ​  TeV 13  TeV 13\text{\,Te\kern-1.00006ptV} as reported in Table 2 under the assumption of 20 % longitudinal polarisation ( α = − 0.2 𝛼 0.2 \alpha=-0.2 ), corresponding to a conservative limit of the ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} polarisation measured at 7 ​  TeV 7  TeV 7\text{\,Te\kern-1.00006ptV}   [ 55 ] .

$p_T$ ($\mathrm{GeV}/c$)
2–3	$0.975$	$0.976$	$0.997$	$1.001$	$0.998$
3–4	$0.955$	$0.983$	$1.005$	$1.013$	$1.012$
4–5	$0.955$	$0.990$	$1.010$	$1.012$	$1.012$
5–6	$0.963$	$0.988$	$1.007$	$1.014$	$1.018$
6–7	$0.960$	$0.990$	$1.008$	$1.004$	$1.020$
7–8	$0.976$	$0.991$	$1.006$	$1.007$	$1.004$
8–9	$0.993$	$0.998$	$1.010$	$1.008$	$1.025$
9–10	$0.973$	$0.998$	$1.013$	$1.003$	$1.033$
10–11	$0.981$	$0.984$	$1.004$	$1.030$	$0.999$
11–12	$0.986$	$0.996$	$1.019$	$0.996$	$1.016$
12–13	$1.023$	$0.995$	$1.004$	$1.010$	$1.005$
13–14	$1.005$	$1.021$	$1.008$	$1.014$	$1.016$
14–15	$1.024$	$1.010$	$1.011$	$1.016$	$1.019$
15–16	$0.948$	$1.011$	$1.020$	$1.014$
16–17	$1.049$	$1.013$	$1.013$		$0.929$
17–18	$1.086$	$1.009$	$1.015$	$1.019$
18–19	$1.004$	$1.013$
19–20

Table 19. Table 19: Multiplicative scaling factors needed to obtain the prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} differential cross-sections from unpolarised cross-section measurements at 7 ​  TeV 7  TeV 7\text{\,Te\kern-1.00006ptV} as reported in Table 4 under the assumption of fully transverse polarisation ( α = + 1 𝛼 1 \alpha=+1 ).

$p_T$ ($\mathrm{GeV}/c$)
3.5–4	$1.087$	$1.288$	$1.164$	$1.111$	$1.079$
4–5	$1.435$	$1.260$	$1.157$	$1.103$	$1.067$
5–6	$1.397$	$1.231$	$1.156$	$1.102$	$1.066$
6–7	$1.362$	$1.220$	$1.144$	$1.102$	$1.071$
7–8	$1.305$	$1.186$	$1.142$	$1.097$	$1.072$
8–9	$1.287$	$1.184$	$1.114$	$1.088$	$1.074$
9–10	$1.250$	$1.150$	$1.124$	$1.086$	$1.071$
10–11	$1.350$	$1.218$	$1.188$	$1.162$	$1.106$
11–12	$1.347$	$1.209$	$1.091$	$1.143$	$1.085$
12–13	$1.246$	$1.222$	$1.127$	$1.139$	$1.103$
13–14	$1.253$	$1.162$	$1.181$	$1.115$	$1.057$

Table 20. Table 20: Multiplicative scaling factors needed to obtain the prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} differential cross-sections from unpolarised cross-section measurements at 7 ​  TeV 7  TeV 7\text{\,Te\kern-1.00006ptV} as reported in Table 4 under the assumption of fully longitudinal polarisation ( α = − 1 𝛼 1 \alpha=-1 ).

$p_T$ ($\mathrm{GeV}/c$)
3.5–4	$0.858$	$0.673$	$0.784$	$0.847$	$0.893$
4–5	$0.583$	$0.697$	$0.792$	$0.854$	$0.908$
5–6	$0.596$	$0.716$	$0.788$	$0.851$	$0.906$
6–7	$0.616$	$0.726$	$0.802$	$0.851$	$0.893$
7–8	$0.627$	$0.734$	$0.780$	$0.837$	$0.874$
8–9	$0.647$	$0.734$	$0.814$	$0.850$	$0.877$
9–10	$0.676$	$0.769$	$0.808$	$0.846$	$0.873$
10–11	$0.743$	$0.821$	$0.833$	$0.868$	$0.896$
11–12	$0.745$	$0.845$	$0.940$	$0.879$	$0.882$
12–13	$0.769$	$0.820$	$0.874$	$0.861$	$0.908$
13–14	$0.802$	$0.872$	$0.834$	$0.922$	$0.947$

Table 21. Table 21: Multiplicative scaling factors needed to obtain the prompt ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} differential cross-sections from unpolarised cross-section measurements at 7 ​  TeV 7  TeV 7\text{\,Te\kern-1.00006ptV} as reported in Table 4 under the assumption of 20 % longitudinal polarisation ( α = − 0.2 𝛼 0.2 \alpha=-0.2 ), corresponding to a conservative limit of the ψ ​ ( 2 ​ S ) 𝜓 2 𝑆 \psi{(2S)} polarisation measured at 7 ​  TeV 7  TeV 7\text{\,Te\kern-1.00006ptV}   [ 55 ]

$p_T$ ($\mathrm{GeV}/c$)
3.5–4	$0.980$	$0.944$	$0.967$	$0.979$	$0.986$
4–5	$0.922$	$0.949$	$0.969$	$0.980$	$0.989$
5–6	$0.923$	$0.951$	$0.967$	$0.979$	$0.988$
6–7	$0.928$	$0.953$	$0.969$	$0.978$	$0.986$
7–8	$0.917$	$0.945$	$0.957$	$0.969$	$0.978$
8–9	$0.928$	$0.944$	$0.963$	$0.972$	$0.980$
9–10	$0.930$	$0.953$	$0.963$	$0.969$	$0.977$
10–11	$1.013$	$1.014$	$1.007$	$1.017$	$1.008$
11–12	$1.016$	$1.023$	$1.020$	$1.013$	$0.988$
12–13	$0.991$	$1.013$	$1.003$	$1.001$	$1.010$
13–14	$1.011$	$1.018$	$1.003$	$1.022$	$1.005$

Equations

$$\begin{array}{rl}\sigma (\text{mbox}prompt\psi (2S),13\text{TeV})& =1.430\pm 0.005\text{(stat)}\pm 0.099\text{(syst)}\text{upmu}\text{b},\\ \sigma (\psi (2S)\text{mbox}-from-b,13\text{TeV})& =0.426\pm 0.002\text{(stat)}\pm 0.030\text{(syst)}\text{upmu}\text{b}.\end{array}$$

$$\begin{array}{rl}\sigma (\text{mbox}prompt\psi (2S),7\text{TeV})& =0.471\pm 0.001\text{(stat)}\pm 0.025\text{(syst)}\text{upmu}\text{b},\\ \sigma (\psi (2S)\text{mbox}-from-b,7\text{TeV})& =0.126\pm 0.001\text{(stat)}\pm 0.008\text{(syst)}\text{upmu}\text{b}.\end{array}$$

$$t_z=\frac{\left(z_{\psi (2S)}-z_{\mathrm{PV}}\right)\times M_{\psi (2S)}}{p_z},$$

$$\frac{{\mathrm{d}}^2\sigma }{\mathrm{d}y\mathrm{d}{p}_{\mathrm{T}}}=\frac{N(p_{\mathrm{T}},y)}{{\epsilon }_{\mathrm{tot}}(p_{\mathrm{T}},y)\times {\mathcal{L}}_{\text{int}}\times \mathcal{B}\times \Delta y\times \Delta p_{\mathrm{T}}},$$

$$\begin{array}{rl}\sigma (\text{mbox}prompt\psi (2S),13\text{TeV})& =1.430\pm 0.005\text{(stat)}\pm 0.099\text{(syst)}\text{upmu}\text{b},\\ \sigma (\psi (2S)\text{mbox}-from-b,13\text{TeV})& =0.426\pm 0.002\text{(stat)}\pm 0.030\text{(syst)}\text{upmu}\text{b}.\end{array}$$

$$\begin{array}{rl}\sigma (\text{mbox}prompt\psi (2S),7\text{TeV})& =0.471\pm 0.001\text{(stat)}\pm 0.025\text{(syst)}\text{upmu}\text{b},\\ \sigma (\psi (2S)\text{mbox}-from-b,7\text{TeV})& =0.126\pm 0.001\text{(stat)}\pm 0.008\text{(syst)}\text{upmu}\text{b}.\end{array}$$

$$\xi \equiv \frac{{\alpha }_{4\pi }(\psi (2S))}{{\alpha }_{4\pi }(J/\psi )}=1.402,$$

$$\frac{\mathcal{B}(b\to \psi (2S)X)}{\mathcal{B}(b\to J/\psi X)}=\xi \frac{\sigma (\psi (2S)\text{mbox}-from-b,13\mathrm{TeV})}{\sigma (J/\psi \text{mbox}-from-b,13\mathrm{TeV})}.$$

$$\frac{\mathcal{B}(b\to \psi (2S)X)}{\mathcal{B}(b\to J/\psi X)}=0.265\pm 0.002(\mathrm{stat})\pm 0.015(\mathrm{syst})\pm 0.006(\mathcal{B}),$$

$$\mathcal{B}(b\to \psi (2S)X)=(3.08\pm 0.02(\mathrm{stat})\pm 0.18(\mathrm{syst})\pm 0.27(\mathcal{B}))\times 10^{-3}.$$

Full text

EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN)

​​​

CERN-EP-2019-150

LHCb-PAPER-2018-049

28 February 2020

**Measurement of ${\psi{(2S)}}$ production cross-sections in proton-proton collisions at $\sqrt{s}=7$ and **$13\text{\,Te\kern-2.07413ptV}$

LHCb collaboration†††Authors are listed at the end of this paper.

The cross-sections of $\psi{(2S)}$ meson production in proton-proton collisions at $\sqrt{s}=13\text{\,Te\kern-1.00006ptV}$ are measured with a data sample collected by the LHCb detector corresponding to an integrated luminosity of $275\text{\,pb}^{-1}$. The production cross-sections for prompt $\psi{(2S)}$ mesons and those for $\psi{(2S)}$ mesons from $b$-hadron decays (${\psi{(2S)}}\mbox{-from-}{b}$) are determined as functions of the transverse momentum, $p_{\mathrm{T}}$, and the rapidity, $y$, of the $\psi{(2S)}$ meson in the kinematic range $2<p_{\mathrm{T}}<20\text{\,Ge\kern-1.00006ptV\!/}c$ and $2.0<y<4.5$. The production cross-sections integrated over this kinematic region are

[TABLE]

A new measurement of $\psi{(2S)}$ production cross-sections in $pp$ collisions at $\sqrt{s}=7\text{\,Te\kern-1.00006ptV}$ is also performed using data collected in 2011, corresponding to an integrated luminosity of $614\text{\,pb}^{-1}$. The integrated production cross-sections in the kinematic range $3.5<p_{\mathrm{T}}<14\text{\,Ge\kern-1.00006ptV\!/}c$ and $2.0<y<4.5$ are

[TABLE]

All results show reasonable agreement with theoretical calculations.

Published in Eur. Phys. J. C 80 (2020) 185

© 2024 CERN for the benefit of the LHCb collaboration. CC-BY-4.0 licence.

1 Introduction

The study of hadronic production of heavy quarkonia can provide important information about quantum chromodynamics (QCD). The production of heavy quark pairs, ${Q{\overline{Q}}}$, can be calculated with perturbative QCD, while the hadronisation of ${Q{\overline{Q}}}$ pairs into heavy quarkonia is nonperturbative and must be determined using input from experimental results. Heavy-quarkonium production therefore probes both perturbative and nonperturbative aspects of QCD by providing stringent tests of theoretical models. Knowledge of hadronic production of heavy quarkonium has been significantly improved in the past forty years [1, 2], but the mechanism behind it is still not fully understood. Colour-singlet model calculations [3, 4, 5, 6, 7, 8, 9] require that the intermediate ${Q{\overline{Q}}}$ state is colourless and has the same $J^{PC}$ quantum numbers as those of the outgoing quarkonium state. In the nonrelativistic QCD (NRQCD) approach [10, 11, *Cho:1995ce], intermediate ${Q{\overline{Q}}}$ states with all possible colour-spin-parity quantum numbers have nonzero probability to be transformed into the desired quarkonium. The transition probability of a ${Q{\overline{Q}}}$ pair into the quarkonium state is described by a long-distance matrix element (LDME), which is assumed to be universal and can be determined from experimental data.

In high-energy proton-proton ($pp$) collisions, charmonium states can be produced directly from hard collisions of partons inside the protons, through the feed-down from excited states, or via weak decays of ${b}$ hadrons. The first two contributions, which cannot be distinguished experimentally, are referred to as prompt production; while the third component can be separated from prompt production by exploiting the lifetime of $b$-hadrons. For prompt ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}$ production the feed-down contribution is large, mostly from radiative decays of ${\chi_{{c}J}}$ ($J=0,1,2$) mesons. This complicates the comparison between theoretical calculations and experimental results. On the contrary, the feed-down contribution to ${\psi{(2S)}}$ mesons is negligible [13], thus theoretical calculations can be directly compared with measurements.

The studies of heavy quarkonium production are crucial to separate the contributions of single parton scattering (SPS) [14] and double parton scattering (DPS) [15] to multiple-quarkonium production. Multiple-quarkonium production through the SPS process shares the same LDMEs as the single quarkonium production, thus providing a new method to test the theoretical calculations. The DPS process can reveal the transverse profile of partons inside the proton. Further theoretical and experimental works provide deeper insights on how to interpret the production mechanism of multiple quarkonia. In particular, additional data help in improving the precision of LDME determination.

The differential cross-sections of inclusive $\psi{(2S)}$ meson production in ${p\overline{p}}$ collisions at centre-of-mass energies of $\sqrt{s}=1.8$ and $1.96\text{\,Te\kern-1.00006ptV}$ were measured by the CDF experiment at the Fermilab Tevatron Collider [16, 17], and in $pp$ collisions at $\sqrt{s}=7\text{\,Te\kern-1.00006ptV}$ [18, 19, 20, 21, 22, 23], $8\text{\,Te\kern-1.00006ptV}$ [23], and $13\text{\,Te\kern-1.00006ptV}$ [24] with LHC data. This paper presents measurements of $\psi{(2S)}$ production cross-sections in ${pp}$ collisions using a data sample collected by LHCb in 2015 (2011) corresponding to an integrated luminosity of $275\pm 11\text{\,pb}^{-1}$ at $\sqrt{s}=13\text{\,Te\kern-1.00006ptV}$ ($614\pm 11\text{\,pb}^{-1}$ at $\sqrt{s}=7\text{\,Te\kern-1.00006ptV}$). The ${\psi{(2S)}}$ mesons from prompt production are abbreviated as “prompt ${\psi{(2S)}}$”, while those from ${b}$-hadron decays are abbreviated as “${\psi{(2S)}}$-from-${b}$”. The $\psi{(2S)}$ mesons are reconstructed through their decay mode ${{\psi{(2S)}}\rightarrow{\mu^{+}\mu^{-}}}$. The double-differential production cross-sections of prompt $\psi{(2S)}$ and $\psi{(2S)}$-from-$b$ as functions of transverse momentum $p_{\mathrm{T}}$ and rapidity $y$ and their integrated production cross-sections are measured, assuming zero polarisation of the $\psi{(2S)}$ meson. The kinematic region of the measurement at $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) is $2<p_{\mathrm{T}}<20\text{\,Ge\kern-1.00006ptV\!/}c$ ($3.5<p_{\mathrm{T}}<14\text{\,Ge\kern-1.00006ptV\!/}c$) and $2.0<y<4.5$. Compared to the previous LHCb measurement at $7\text{\,Te\kern-1.00006ptV}$ using 2010 data [19], the new analysis at $7\text{\,Te\kern-1.00006ptV}$ has several advantages: the 2011 data sample is much larger than that in the previous measurement corresponding to an integrated luminosity of $36\text{\,pb}^{-1}$, the previous measurement did not provide the $\psi{(2S)}$ production cross-section as a function of the rapidity $y$ because of the limited sample size, and the same final state and offline selection criteria as those in the $13\text{\,Te\kern-1.00006ptV}$ measurement are used in the new $7\text{\,Te\kern-1.00006ptV}$ measurement. This guarantees that the maximum number of systematic uncertainties cancel in the cross-section ratio between $13\text{\,Te\kern-1.00006ptV}$ and $7\text{\,Te\kern-1.00006ptV}$, which is measured in the present analysis. This represents a more stringent test of the theoretical models, since many of the experimental and theoretical uncertainties cancel. Finally, the $\psi{(2S)}$ meson differential production cross-sections are compared with those of the ${J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}$ meson at $\sqrt{s}=13\text{\,Te\kern-1.00006ptV}$ [25].

2 Detector and simulation

The LHCb detector [26, 27] is a single-arm forward spectrometer covering the pseudorapidity range $2<\eta<5$, designed for the study of particles containing $b$ or $c$ quarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector surrounding the $pp$ interaction region [28], a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about $4{\mathrm{\,Tm}}$, and three stations of silicon-strip detectors and straw drift tubes [29, 30] placed downstream of the magnet. The tracking system provides a measurement of the momentum, $p$, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200$\text{\,Ge\kern-1.00006ptV\!/}c$. The minimum distance of a track to a primary vertex (PV), the impact parameter (IP), is measured with a resolution of $(15+29/p_{\mathrm{T}})\,\upmu\text{m}$, where $p_{\mathrm{T}}$ is in $\text{\,Ge\kern-1.00006ptV\!/}c$. Different types of charged hadrons are distinguished using information from two ring-imaging Cherenkov detectors [31]. Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillating-pad (SPD) and preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers [32]. The online event selection is performed by a trigger [33], which consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies a full event reconstruction.

Simulated samples are used to evaluate the $\psi{(2S)}$ detection efficiency. In the simulation, ${pp}$ collisions are generated using Pythia 8 [34, 35] with a specific LHCb configuration [36]. Decays of unstable particles are described by EvtGen [37], in which final-state radiation is generated using Photos [38]. Both the leading-order colour-singlet and colour-octet contributions are included in the generated prompt charmonium states [36, 39]. These states are generated with zero polarisation. The interaction of the generated particles with the detector, and its response, are implemented using the Geant4 toolkit [40, *Agostinelli:2002hh] as described in Ref. [42].

3 Selection of $\psi{(2S)}$ candidates

The decay channel ${\psi{(2S)}}\rightarrow{\mu^{+}\mu^{-}}$ is used in the measurements of the $\psi{(2S)}$ production cross-sections at both $13\text{\,Te\kern-1.00006ptV}$ and $7\text{\,Te\kern-1.00006ptV}$. The same strategy is used for both analyses, except for different trigger requirements. The hardware trigger selects events that contain two tracks consistent with muon hypotheses, and the product of the transverse momenta of the two muons is required to be greater than $(1.3\text{\,Ge\kern-1.00006ptV\!/}c)^{2}$. At the software trigger stage the two muons are required to be oppositely charged, to have good track quality, to form a good-quality vertex, and to each have a momentum larger than $6\text{\,Ge\kern-1.00006ptV\!/}c$. The invariant mass of the $\psi{(2S)}$ candidates is required to be within the range $3566<m_{{\mu^{+}\mu^{-}}}<3806\text{\,Me\kern-1.00006ptV\!/}c^{2}$. The transverse momentum of each muon is required to be larger than $0.3\text{\,Ge\kern-1.00006ptV\!/}c$ ($0.5\text{\,Ge\kern-1.00006ptV\!/}c$) and that of the $\psi{(2S)}$ candidate is required to be larger than $2\text{\,Ge\kern-1.00006ptV\!/}c$ ($3.5\text{\,Ge\kern-1.00006ptV\!/}c$) for the $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) data trigger. Due to the different triggers, the $\psi{(2S)}$ candidates are selected in different $p_{\mathrm{T}}$ ranges. For the $13\text{\,Te\kern-1.00006ptV}$ data taking, an alignment and calibration of the detector is performed in near real-time [43] and updated constants are made available for the trigger.

To suppress the background associated to random combination of tracks (combinatorial) more stringent criteria are applied offline on the ${\psi{(2S)}}$ vertex fit quality, the muon kinematics and particle identification requirements. Each muon must have $p_{\mathrm{T}}>1.2\text{\,Ge\kern-1.00006ptV\!/}c$ and $2.0<\eta<4.9$. At least one PV should be reconstructed in the event from at least four tracks in the vertex detector.

For events with more than one PV, the ${\psi{(2S)}}$ candidate is associated to the PV for which the difference in the $\chi^{2}$ of the PV fit with and without the ${\psi{(2S)}}$ candidate is the smallest. This is equivalent to select the vertex with respect to which the signal candidate has the smallest impact parameter, compared to resolution. Simulation studies have shown that, using the above procedure, the fraction of candidates associated to the wrong PV is $0.3\%$, which is negligible. To select $\psi{(2S)}$ candidates, additional requirements on the pseudo decay time, $t_{z}$, $|t_{z}|<10\text{\,ps}$, and its uncertainty, $\sigma_{t_{z}}$, $\sigma_{t_{z}}<0.3\text{\,ps}$, are applied. The pseudo decay time $t_{z}$ is defined as

[TABLE]

where $z_{\psi{(2S)}}$ ($z_{\mathrm{PV}}$) is the $z$ coordinate of the reconstructed ${\psi{(2S)}}$ decay vertex (the PV), $p_{z}$ is the $z$-component of the measured ${\psi{(2S)}}$ momentum, and $M_{{\psi{(2S)}}}$ is the world average ${\psi{(2S)}}$ mass [13]. The $z$-axis is the direction of the proton beam pointing downstream into the LHCb acceptance [26]. The pseudo decay time defined above provides a good approximation of the $b$-hadron decay time [44] and is used to separate prompt $\psi{(2S)}$ and $\psi{(2S)}$-from-$b$ candidates.

4 Cross-section determination

The double-differential production cross-section for prompt $\psi{(2S)}$ or $\psi{(2S)}$-from-${b}$ in a given $(p_{\mathrm{T}},y)$ bin is defined as

[TABLE]

where $N(p_{\mathrm{T}},y)$ is the signal yield, ${\varepsilon_{\mathrm{tot}}}(p_{\mathrm{T}},y)$ is the total detection efficiency of the ${{\psi{(2S)}}\rightarrow{\mu^{+}\mu^{-}}}$ decay evaluated independently for prompt $\psi{(2S)}$ or $\psi{(2S)}$-from-${b}$ in the given $(p_{\mathrm{T}},y)$ bin, $\mathcal{L}_{\text{int}}$ is the integrated luminosity, ${\mathcal{B}}$ is the branching fraction of the decay ${{\psi{(2S)}}\rightarrow{\mu^{+}\mu^{-}}}$, and $\Delta p_{\mathrm{T}}=1\text{\,Ge\kern-1.00006ptV\!/}c$ and $\Delta y=0.5$ are the bin widths. The integrated luminosity is determined using the beam-gas imaging and, for the $7\text{\,Te\kern-1.00006ptV}$ data, also the van der Meer scan methods [45]. Assuming lepton universality in electromagnetic decays, ${\mathcal{B}}({{\psi{(2S)}}\rightarrow{e^{+}e^{-}}})=(7.89\pm 0.17)\times 10^{-3}$ [46] is used in Eq. 2, taking advantage of the much smaller uncertainty compared to the ${{\psi{(2S)}}\rightarrow{\mu^{+}\mu^{-}}}$ decay. The difference of the two branching fractions introduced by the mass difference between electrons and muons is negligible.

The yields of prompt $\psi{(2S)}$ and $\psi{(2S)}$-from-$b$ candidates in each $(p_{\mathrm{T}},y)$ bin are determined from a two-dimensional extended unbinned maximum-likelihood fit to the distributions of the invariant mass, $m_{{\mu^{+}\mu^{-}}}$, and $t_{z}$ of the $\psi{(2S)}$ candidates. The correlation between $m_{{\mu^{+}\mu^{-}}}$ and $t_{z}$ is found to be negligible. The invariant-mass distribution of the signal candidates in each bin is described by the sum of two Crystal Ball (CB) functions [47] with a common mean value and different widths. The parameters of the power-law tails, the relative fractions and the difference between the widths of the two CB functions are fixed to values obtained from simulation, leaving the mean value and the width of one of the CB functions as free parameters. The invariant-mass distribution of the combinatorial background is described by an exponential function with the slope parameter free to vary in the fit. The $t_{z}$ distribution of prompt $\psi{(2S)}$ mesons is described by a Dirac $\delta$ function at $t_{z}=0$, and that of $\psi{(2S)}$-from-$b$ by an exponential function, both convolved with the sum of two Gaussian functions. A $\psi{(2S)}$ candidate can also be associated to a wrong PV, resulting in a long tail component in the $t_{z}$ distribution. This shape is modelled from data by calculating $t_{z}$ between the $\psi{(2S)}$ candidate from a given event and the closest PV in the next event of the sample. The background $t_{z}$ distribution is parametrised with an empirical function based on the observed shape of the $t_{z}$ distribution in the $\psi{(2S)}$ mass sidebands ($3566<m_{{\mu^{+}\mu^{-}}}<3620\text{\,Me\kern-1.00006ptV\!/}c^{2}$ and $3750<m_{{\mu^{+}\mu^{-}}}<3806\text{\,Me\kern-1.00006ptV\!/}c^{2}$). It is parametrised as the combination of a Dirac $\delta$ function and the sum of five exponential functions, three for positive $t_{z}$ and two for negative $t_{z}$. This sum is convolved with the sum of two Gaussian functions. All parameters of the background $t_{z}$ distribution are fixed to values determined from the $\psi{(2S)}$ mass sidebands independently in each $(p_{\mathrm{T}},y)$ bin. Figure 1 shows as an example the $m_{{\mu^{+}\mu^{-}}}$ and $t_{z}$ distributions in the kinematic bin corresponding to $5<p_{\mathrm{T}}<6\text{\,Ge\kern-1.00006ptV\!/}c$ and $2.5<y<3.0$ for the $13\text{\,Te\kern-1.00006ptV}$ data sample. The one-dimensional projections of the fit result are also presented. The total signal yields of prompt $\psi{(2S)}$ and $\psi{(2S)}$-from-$b$ in the kinematic range for the $13\text{\,Te\kern-1.00006ptV}$ sample are $(440.7\pm 1.2)\times 10^{3}$ and $(140.0\pm 0.5)\times 10^{3}$, and for the $7\text{\,Te\kern-1.00006ptV}$ sample are $(433.9\pm 0.9)\times 10^{3}$ and $(115.1\pm 0.4)\times 10^{3}$, respectively.

The total efficiency, ${\varepsilon_{\mathrm{tot}}}$, in each kinematic bin is determined as the product of the geometrical acceptance of the detector and the efficiencies of particle reconstruction, event selection, muon identification and trigger requirements. The detector acceptance, selection and trigger efficiencies are calculated using simulated samples in each $(p_{\mathrm{T}},y)$ bin, independently for prompt $\psi{(2S)}$ and $\psi{(2S)}$-from-$b$. The trigger efficiencies are also validated using data, as explained in Section 5. The track reconstruction and the muon-identification efficiencies are evaluated using simulated samples and calibrated with data. The efficiencies of prompt $\psi{(2S)}$ and those of $\psi{(2S)}$-from-$b$ are very similar.

5 Systematic uncertainties

A variety of sources of systematic uncertainty are studied as described below and are summarised in Table 1. For the uncertainties that vary in kinematic bins, the largest uncertainties always appear in the bins with small sample sizes.

The uncertainty related to the modelling of the signal mass shape is studied by replacing the baseline model with a kernel-density estimated distribution [48] obtained from the simulated sample in each kinematic bin. In order to account for the resolution difference between data and simulation, a Gaussian function is used to smear the shape of the distribution in simulation. The relative difference of the signal yield in each kinematic bin, $0.0$–$4.1\%$ ($0.0$–$8.5\%$) for the $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) sample, is taken as the systematic uncertainty due to signal mass shape.

Due to the presence of final-state radiation in the ${\psi{(2S)}}\rightarrow{\mu^{+}\mu^{-}}$ decay, a fraction of $\psi{(2S)}$ candidates fall outside the mass window used to determine the signal yields. The efficiency of the selection of the mass window is estimated using simulated samples, and the imperfect modelling of the radiation is studied by comparisons of the radiative tails between simulation and data, from which an uncertainty of $1.0\%$ is assigned to the cross-sections in all kinematic bins.

The track detection efficiencies are determined from a simulated sample in each $(p_{\mathrm{T}},y)$ bin of the $\psi{(2S)}$ meson, and are corrected by using ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}\rightarrow{\mu^{+}\mu^{-}}$ decays reconstructed in a control data sample and in simulation. These efficiencies are calculated as functions of $p$ and $\eta$ with a tag-and-probe approach [49]. The uncertainties due to the finite size of the control samples are propagated to the results using a large number of pseudoexperiments. In each pseudoexperiment, a new efficiency-correction ratio in each $(p_{\mathrm{T}},y)$ bin is generated according to a Gaussian distribution where the original ratio and its uncertainty are used as the Gaussian mean and standard deviation, respectively. The contribution to the systematic uncertainty in each kinematic bin of $\psi{(2S)}$ mesons varies from $0.1\%$ ($0.7\%$) to $2.4\%$ ($3.0\%$) for the $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) data sample. The distribution of the number of SPD hits in simulation is weighted to match that in data to correct the effect of the detector occupancy. As a crosscheck the number of tracks is used as an alternative weighting variable. The tracking efficiencies are found to be different when different variables are used. Therefore, an additional systematic uncertainty of $0.8\%$ ($0.4\%$) per muon track is assigned for the $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) sample.

The muon identification efficiency is determined from simulation and calibrated with a data sample of ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}\rightarrow{\mu^{+}\mu^{-}}$ decays. The statistical uncertainty due to the finite size of the calibration sample is propagated to the final results using pseudoexperiments. The resulting uncertainties vary from $0.1\%$ ($0.7\%$) to $1.1\%$ ($8.9\%$) in different $(p_{\mathrm{T}},y)$ bins for the $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) sample. The uncertainty related to the kinematic binning scheme of the calibration samples is studied by changing the size and the boundaries of the $p$ and $\eta$ bins, and number of SPD hits. This leads to systematic uncertainties of $0.1$–$4.6\%$ ($0.4$–$5.4\%$) for the $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) sample.

The trigger efficiency is determined from simulated samples. To estimate the systematic uncertainty, a tag-and-probe method is used to estimate the trigger efficiencies in each ($p_{\mathrm{T}},y$) bin of a $\psi{(2S)}$ data sample that is independent of the detection of $\psi{(2S)}$ signals [33]. The same procedure is applied to the simulated $\psi{(2S)}$ samples, and the relative difference of efficiencies between data and simulation in each kinematic bin, $0.1$–$9.3\%$ ($0.0$–$4.4\%$) for the $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) sample, is taken as a systematic uncertainty.

The $p_{\mathrm{T}}$ and $y$ distributions of $\psi{(2S)}$ mesons in simulation and in data could be different within each kinematic bin due to the finite bin size, causing differences in efficiencies. The possible discrepancy is studied by weighting the kinematic distribution in simulation to match that in data. All efficiencies are recalculated, and the relative differences of the total efficiencies between the new and the nominal results, which are found to be in the range $0.0$–$2.0\%$ ($0.0$–$4.9\%$) for the $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) sample, are taken as systematic uncertainties.

The integrated luminosity is determined using the beam-gas imaging method for the $13\text{\,Te\kern-1.00006ptV}$ data sample, and by a combination of the beam-gas imaging and van der Meer scan methods [45] for the $7\text{\,Te\kern-1.00006ptV}$ data sample. The uncertainty associated with the luminosity determination is $3.9\%$ ($1.7\%$) for the $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) sample. The uncertainty of the branching fraction of the ${{\psi{(2S)}}\rightarrow{e^{+}e^{-}}}$ decay, $2.2\%$, is taken as a systematic uncertainty [46]. The limited size of the simulated sample in each bin leads to uncertainties of $0.7$–$11.5\%$ ($1.3$–$13.1\%$) for prompt $\psi{(2S)}$ and $0.8$–$5.7\%$ ($1.2$–$9.5\%$) for $\psi{(2S)}$-from-${b}$ for the $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) sample, and are smaller than or comparable with the data statistical uncertainty in each bin.

There are sources of systematic uncertainties that are related to the $t_{z}$ variable, the effects of which are notable for $\psi{(2S)}$-from-$b$ and are negligible for prompt $\psi{(2S)}$. The modelling of the $t_{z}$ resolution is modified by adding a third Gaussian to the nominal resolution model. The variation in the $\psi{(2S)}$-from-$b$ fraction $F_{b}$ is found to be negligible. An alternative method is adopted to estimate the systematic uncertainty due to the modelling of the background $t_{z}$ distribution. In this method, the background distribution is obtained with the sPlot technique [50] using the invariant mass as the discriminating variable. The $t_{z}$ distribution is then parametrised for the two-dimensional fits to obtain the fraction $F_{b}$. The relative difference of $F_{b}$ in each kinematic bin between the two methods is taken as a systematic uncertainty. The total systematic uncertainty on the ${\psi{(2S)}}\text{-from-}b$ cross-section related to the $t_{z}$ fit model is $0.1$–$8.4\%$ ($0.1$–$9.2\%$) for the $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) sample.

6 Results

6.1 Production cross-sections

The double-differential production cross-sections for prompt $\psi{(2S)}$ and $\psi{(2S)}$-from-$b$ are measured as functions of $p_{\mathrm{T}}$ and $y$ assuming no polarisation of $\psi{(2S)}$ mesons. The results are shown in Figs. 2 and 3, respectively. The corresponding values are listed in Tables 2, 3, 4, and 5 in Appendix A.

By integrating the double-differential results over $y$ in the range $2.0<y<4.5$, the differential production cross-sections of prompt $\psi{(2S)}$ and $\psi{(2S)}$-from-$b$ as functions of $p_{\mathrm{T}}$ are shown in Fig. 4. The results of prompt $\psi{(2S)}$ production are compared with the theoretical calculations based on NRQCD [51], and those of $\psi{(2S)}$-from-$b$ are compared with the fixed-order-plus-next-leading-logarithm (FONLL) calculations [52]. The differential cross-section as function of $y$ at $13\text{\,Te\kern-1.00006ptV}$ ($7\text{\,Te\kern-1.00006ptV}$) is obtained by integrating the double-differential results over $p_{\mathrm{T}}$ in the range $2<p_{\mathrm{T}}<20\text{\,Ge\kern-1.00006ptV\!/}c$ ($3.5<p_{\mathrm{T}}<14\text{\,Ge\kern-1.00006ptV\!/}c$). The results are presented in Fig. 5. The theoretical calculations based on FONLL are shown for $\psi{(2S)}$-from-$b$. The NRQCD calculations are omitted since they are not reliable in the low $p_{\mathrm{T}}$ region [53]. The values of the differential cross-sections are shown in Tables 6, 7, 8, and 9 in Appendix A. In the NRQCD calculations, only the dominant uncertainties associated with the LDMEs are considered [51]. The FONLL calculations include the uncertainty due to $b$-quark mass and the scales of renormalisation and factorisation. The NRQCD calculations show reasonable agreement with experimental data for $p_{\mathrm{T}}>7\text{\,Ge\kern-1.00006ptV\!/}c$. The FONLL calculations agree well with the measurements. The production cross-sections of prompt $\psi{(2S)}$ and $\psi{(2S)}$-from-$b$ integrated in the kinematic range $2.0<y<4.5$ and $2<p_{\mathrm{T}}<20\text{\,Ge\kern-1.00006ptV\!/}c$ at $13\text{\,Te\kern-1.00006ptV}$, are measured to be:

[TABLE]

The production cross-sections of prompt $\psi{(2S)}$ and $\psi{(2S)}$-from-$b$ integrated in the kinematic range $2.0<y<4.5$ and $3.5<p_{\mathrm{T}}<14\text{\,Ge\kern-1.00006ptV\!/}c$ at $7\text{\,Te\kern-1.00006ptV}$, are measured to be:

[TABLE]

As mentioned above, these results are obtained under the assumption of zero polarisation of $\psi{(2S)}$ mesons. Possible polarisation of $\psi{(2S)}$ meson would affect the detection efficiency. This effect is studied for extreme cases of fully transverse and fully longitudinal polarisation corresponding to the parameter $\alpha$ be equal to $+1$ or $-1$, respectively, within the helicity frame [54] approach. Also the polarisation case of $\alpha=-0.2$, corresponding to a conservative limit of the $\psi{(2S)}$ polarisation measured at $7\text{\,Te\kern-1.00006ptV}$ [55], is considered. Resulting scaling factors for prompt $\psi{(2S)}$ production cross-sections are listed in Appendix B in Tables 16, 17 and 18 for $13\text{\,Te\kern-1.00006ptV}$ results, and in Tables 19, 20 and 21 for $7\text{\,Te\kern-1.00006ptV}$ results, respectively.

6.2 Fraction of $\psi{(2S)}$-from-${b}$ mesons

The fraction of $\psi{(2S)}$-from-${b}$ is $F_{b}\equiv N_{b}/(N_{b}+N_{p})$, where $N_{p}$ is the efficiency-corrected signal yield of prompt $\psi{(2S)}$ and $N_{b}$ is that of $\psi{(2S)}$-from-${b}$. The fractions $F_{b}$ as functions of $p_{\mathrm{T}}$ and $y$ are shown in Fig. 6. The corresponding values are presented in Table 10 in Appendix A. Only statistical uncertainties are shown owing to the cancellation of most systematic contributions, except for that due to the $t_{z}$ fit, which is negligible. For each $y$ bin, the fraction increases with increasing $p_{\mathrm{T}}$ of the $\psi{(2S)}$ mesons. For each $p_{\mathrm{T}}$ bin, the fraction decreases with increasing $y$ of the $\psi{(2S)}$ mesons.

6.3 Comparison with ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}$ results at $13\text{\,Te\kern-1.00006ptV}$

The production cross-sections of $\psi{(2S)}$ mesons at $13\text{\,Te\kern-1.00006ptV}$ are compared with those of ${J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}$ mesons measured by LHCb at $\sqrt{s}=13\text{\,Te\kern-1.00006ptV}$ in the range $0<p_{\mathrm{T}}<14\text{\,Ge\kern-1.00006ptV\!/}c$ and $2.0<y<4.5$ [25], where the ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}$ meson is also assumed to be produced with zero polarisation. The ratio, $R_{{\psi{(2S)}}/{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}}$, of the differential production cross-sections in the common range between prompt $\psi{(2S)}$ and prompt ${J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}$ mesons is shown in Fig. 7 as a function of $p_{\mathrm{T}}$ ($y$) integrated over $2.0<y<4.5$ ($2<p_{\mathrm{T}}<14\text{\,Ge\kern-1.00006ptV\!/}c$). The NRQCD calculation of $R_{{\psi{(2S)}}/{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}}$ for prompt productions [51] is also shown. The ratio of production cross-sections between $\psi{(2S)}$-from-$b$ and ${J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}$-from-$b$ is shown in Fig. 8 as a function of $p_{\mathrm{T}}$ ($y$) integrated over $2.0<y<4.5$ ($2<p_{\mathrm{T}}<14\text{\,Ge\kern-1.00006ptV\!/}c$). The FONLL calculations [56] are compared to the measured values. To calculate these ratios from the measured cross-sections of $\psi{(2S)}$ and ${J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}$ mesons, the systematic uncertainties due to the luminosity, the tracking correction, and the fit model are considered to be fully correlated. All other uncertainties are assumed to be uncorrelated. The numerical results of the measured ratios are listed in Tables 12 and 13 in Appendix A. The FONLL prediction agrees well with the experimental data for the production cross-section ratio between $\psi{(2S)}$-from-$b$ and ${J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}$ mesons from $b$-hadron decays, while the NRQCD predictions show reasonable agreement with the measurements for prompt $\psi{(2S)}$ and prompt ${J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}$.

6.4 Comparison between $13\text{\,Te\kern-1.00006ptV}$ and $7\text{\,Te\kern-1.00006ptV}$

The production cross-sections of ${\psi{(2S)}}$ mesons in $pp$ collisions at $13\text{\,Te\kern-1.00006ptV}$ and $7\text{\,Te\kern-1.00006ptV}$ are compared by means of their ratio, $R_{13/7}$. Figures 9 and 10 show the ratios as functions of $p_{\mathrm{T}}$ integrated over $2.0<y<4.5$ and as functions of $y$ integrated over $3.5<p_{\mathrm{T}}<14\text{\,Ge\kern-1.00006ptV\!/}c$ for prompt $\psi{(2S)}$ and $\psi{(2S)}$-from-$b$. The NRQCD (FONLL) calculations of $R_{13/7}$ for prompt $\psi{(2S)}$ ($\psi{(2S)}$-from-$b$) are also shown in the left (right) panel for comparison. Both FONLL and NRQCD predictions on $R_{13/7}$ agree well with the corresponding experimental data. The measured ratios are also presented in Tables 14 and 15 in Appendix A.

For both the theoretical calculations and the experimental measurements, some of the uncertainties in the ratio cancel, which allows for a more precise comparison to theory. In the calculation of these ratios from the measured $\psi{(2S)}$ production cross-sections at $13\text{\,Te\kern-1.00006ptV}$ and $7\text{\,Te\kern-1.00006ptV}$ the systematic uncertainty related to the branching fraction is cancelled. The uncertainties due to the luminosity, the fit model and the tracking correction are partially correlated. Other uncertainties are assumed to be uncorrelated.

6.5 Measurement of the inclusive ${b}\rightarrow{\psi{(2S)}}X$ branching fraction

The reported results of the cross-section $\sigma({\psi{(2S)}}\mbox{-from-}{b},13\text{\,Te\kern-1.00006ptV})$, in combination with the previous results about ${J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}$ production [25], can be used to determine the inclusive branching fraction $\mathcal{B}(b\rightarrow\psi(2S)X)$. To achieve this, both results must be extrapolated to the full phase space, as they are measured only for a limited range of phase space. The extrapolation factors $\alpha_{4\pi}({\psi{(2S)}})$ and $\alpha_{4\pi}({{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}})$ are determined with LHCb-tuned versions of Pythia 8 [34] for the $\psi{(2S)}$ and of Pythia 6 [35] for the ${J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}$. The factors $\alpha_{4\pi}({\psi{(2S)}})$ and $\alpha_{4\pi}({{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}})$ are found to be 7.29 and 5.20, respectively. In the ratio of the two factors,

[TABLE]

most of the theoretical uncertainties are expected to cancel. Alternatively, the correction factor $\xi$ can be obtained using FONLL calculations which uses different parton distribution functions. The values of $\xi$ obtained from the two methods differ by 2.89 %.

With the definition of the ratio $\xi$, the $\mathcal{B}(b\rightarrow{\psi{(2S)}}X)$ branching fraction can be obtained from the ratio

[TABLE]

By inserting the value $\sigma({{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}\mbox{-from-}{b},\,\mathrm{13~{}TeV})=2.25\pm 0.01(\mathrm{stat})\pm 0.14(\mathrm{syst})\,\upmu\text{b}$ [25] and the value of $\xi$, the ratio of the branching fractions is

[TABLE]

where possible correlations between uncertainties originating from ${\psi{(2S)}}\text{-from-}b$ and ${J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}$-from-$b$, respectively, are taken into account. The last uncertainty is from the uncertainty of the branching fractions $\mathcal{B}({\psi{(2S)}}\rightarrow e^{+}e^{-})$ and $\mathcal{B}({{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}\rightarrow\mu^{+}\mu^{-})$. Using the known value $\mathcal{B}(b\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}X)=(1.16\pm 0.10)\times 10^{-2}$ [13], one obtains

[TABLE]

This result is in agreement with the world-average value [13]. The $\mathcal{B}(b\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}X)$ uncertainty dominates the total uncertainty from the branching fractions.

7 Conclusions

The production cross-sections of $\psi{(2S)}$ mesons in proton-proton collisions at a centre-of-mass energy of $13\text{\,Te\kern-1.00006ptV}$ are reported with a data sample corresponding to an integrated luminosity of $275\pm 11\text{\,pb}^{-1}$, collected by the LHCb detector in 2015. The double-differential cross-sections, as functions of $p_{\mathrm{T}}$ and $y$ of the $\psi{(2S)}$ meson in the range of $2<p_{\mathrm{T}}<20\text{\,Ge\kern-1.00006ptV\!/}c$ and $2.0<y<4.5$, are determined for prompt $\psi{(2S)}$ mesons and $\psi{(2S)}$ mesons from ${b}$-hadron decays. A new measurement of the branching fraction $\mathcal{B}(b\rightarrow\psi(2S)X)$ is presented, which is in agreement with the world average [13]. The measured prompt $\psi{(2S)}$ production cross-section as a function of transverse momentum is in good agreement in the high $p_{\mathrm{T}}$ region with theoretical calculations in the NRQCD framework. Theoretical predictions based on the FONLL calculations describe well the measured cross-sections for $\psi{(2S)}$ mesons from ${b}$-hadron decays.

A new measurement of $\psi{(2S)}$ production cross-sections at $7\text{\,Te\kern-1.00006ptV}$ is performed using the 2011 data sample corresponding to an integrated luminosity of $614\pm 11\text{\,pb}^{-1}$. The new result provides a significantly reduced uncertainty compared to the previous independent LHCb result [19].

The cross section ratios between $13\text{\,Te\kern-1.00006ptV}$ and $7\text{\,Te\kern-1.00006ptV}$ show reasonable agreement with theoretical calculations.

Acknowledgements

We thank Kuang-Ta Chao and Yan-Qing Ma for frequent and interesting discussions on the production of ${\psi{(2S)}}$ mesons. We express our gratitude to our colleagues in the CERN accelerator departments for the excellent performance of the LHC. We thank the technical and administrative staff at the LHCb institutes. We acknowledge support from CERN and from the national agencies: CAPES, CNPq, FAPERJ and FINEP (Brazil); MOST and NSFC (China); CNRS/IN2P3 (France); BMBF, DFG and MPG (Germany); INFN (Italy); NWO (Netherlands); MNiSW and NCN (Poland); MEN/IFA (Romania); MSHE (Russia); MinECo (Spain); SNSF and SER (Switzerland); NASU (Ukraine); STFC (United Kingdom); DOE NP and NSF (USA). We acknowledge the computing resources that are provided by CERN, IN2P3 (France), KIT and DESY (Germany), INFN (Italy), SURF (Netherlands), PIC (Spain), GridPP (United Kingdom), RRCKI and Yandex LLC (Russia), CSCS (Switzerland), IFIN-HH (Romania), CBPF (Brazil), PL-GRID (Poland) and OSC (USA). We are indebted to the communities behind the multiple open-source software packages on which we depend. Individual groups or members have received support from AvH Foundation (Germany); EPLANET, Marie Skłodowska-Curie Actions and ERC (European Union); ANR, Labex P2IO and OCEVU, and Région Auvergne-Rhône-Alpes (France); Key Research Program of Frontier Sciences of CAS, CAS PIFI, and the Thousand Talents Program (China); RFBR, RSF and Yandex LLC (Russia); GVA, XuntaGal and GENCAT (Spain); the Royal Society and the Leverhulme Trust (United Kingdom).

Appendices

Appendix A Result tables

[FIGURE:]

[FIGURE:]

[FIGURE:]

[FIGURE:]

Appendix B Scaling factors for alternative polarisation scenarios