Measurement of $B^0$, $B^0_s$, $B^+$ and $\Lambda^0_b$ production asymmetries in 7 and 8 TeV proton-proton collisions

TL;DR

This paper reports measurements of production asymmetries for various $b$ hadrons in proton-proton collisions at 7 and 8 TeV, providing insights into hadronization processes and CP violation effects.

Contribution

The study presents the first comprehensive measurement of $B^0$, $B^0_s$, $B^+$, and $\Lambda^0_b$ production asymmetries at these energies as functions of transverse momentum and rapidity.

Findings

Measured asymmetries as functions of kinematic variables.

Determined overall production asymmetries integrated over phase space.

Provided data to test QCD models of hadronization.

Abstract

The $B^0$, $B^0_s$, $B^+$ and $\Lambda^0_b$ hadron production asymmetries are measured using a data sample corresponding to an integrated luminosity of 3.0 fb$^{-1}$, collected by the LHCb experiment in proton-proton collisions at centre-of-mass energies of 7 and 8 TeV. The measurements are performed as a function of transverse momentum and rapidity of the $b$ hadrons within the LHCb detector acceptance. The overall production asymmetries, integrated over transverse momentum and rapidity, are also determined.

Tables

Table 1. Table 1: Values of the various physical inputs used in the fits, as reported in Ref.  [ 33 ] .

Parameter	Value
$\Delta m_d\left[{\mathrm{ps}}^{-1}\right]$	$0.5065\pm 0.0019$
$\Delta m_s\left[{\mathrm{ps}}^{-1}\right]$	$17.757\pm 0.021$
${\Gamma }_d\left[{\mathrm{ps}}^{-1}\right]$	$0.6579\pm 0.0017$
${\Gamma }_s\left[{\mathrm{ps}}^{-1}\right]$	$0.6645\pm 0.0018$
$\Delta {\Gamma }_s\left[{\mathrm{ps}}^{-1}\right]$	$0.083\pm 0.006$
${\left|q/p\right|}_{B^0}$	$1.0007\pm 0.0009$
${\left|q/p\right|}_{B_s^0}$	$1.0038\pm 0.0021$

Table 2. Table 2: Values of signal yields, A P subscript 𝐴 P A_{\rm P} and A D subscript 𝐴 D A_{\rm D} obtained from global fits for the two neutral B 𝐵 B -meson decays under study.

Parameter	$\sqrt{s}$ = 7$\mathrm{TeV}$		$\sqrt{s}$ = 8$\mathrm{TeV}$
	$B^0\to J/\psi K^{\ast 0}$	$B_s^0\to D_s^{-}{\pi }^{+}$	$B^0\to J/\psi K^{\ast 0}$	$B_s^0\to D_s^{-}{\pi }^{+}$
$N_{\mathrm{sig}}$	$\mathrm{95\hspace{0.17em}122}\pm 369$	$\mathrm{16\hspace{0.17em}932}\pm 174$	$\mathrm{221\hspace{0.17em}973}\pm 569$	$\mathrm{36\hspace{0.17em}726}\pm 250$
$A_P$	$-0.0113\pm 0.0063$	$-0.0001\pm 0.0166$	$-0.0109\pm 0.0042$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.0081}\pm 0.0111$
$A_D$	$-0.0098\pm 0.0046$	$-0.0143\pm 0.0086$	$-0.0056\pm 0.0030$	$-0.0103\pm 0.0058$

Table 3. Table 3: Values of signal yields and raw asymmetries obtained from global fits in the case of the B + → J / ψ ​ K + → superscript 𝐵 𝐽 𝜓 superscript 𝐾 {{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}} decay.

Parameter	$\sqrt{s}$ = 7$\mathrm{TeV}$	$\sqrt{s}$ = 8$\mathrm{TeV}$
$N_{\mathrm{sig}}$	$\mathrm{265\hspace{0.17em}574}\pm 576$	$\mathrm{619\hspace{0.17em}800}\pm 908$
$A_{\mathrm{raw}}$	$-0.017\pm 0.002$	$-0.014\pm 0.001$

Table 4. Table 4: Absolute values of systematic uncertainties for integrated production asymmetries. The total systematic uncertainties are obtained by summing the individual contributions in quadrature.

	Uncertainty [$\sqrt{s}$ = 7$\mathrm{TeV}$]
Source	$A_P\left(B^{+}\right)$	$A_P\left(B^0\right)$	$A_P\left(B_s^0\right)$	$A_P\left({\Lambda }_b^0\right)$
Signal mass shape	0.0016	0.0005	0.0036	0.0024
Decay-time bias	0.0000	0.0000	0.0008	0.0004
$\Delta m_d$, $\Delta m_s$	0.0000	0.0001	0.0014	0.0007
Decay-time resolution	0.0000	0.0000	0.0026	0.0014
Final-state radiation	0.0000	0.0001	0.0000	0.0001
Decay-time reconstruction efficiency	0.0000	0.0001	0.0000	0.0001
Combinatorial background mass shape	0.0003	0.0000	0.0004	0.0003
Partially reconstructed background mass shape	0.0000	0.0000	0.0029	0.0015
$\Delta {\Gamma }_s$	0.0000	0.0000	0.0000	0.0000
$A_D\left(K^{+}\right)$	0.0018	0.0000	0.0000	0.0013
${\left|q/p\right|}_{B^0}$, ${\left|q/p\right|}_{B_s^0}$	0.0000	0.0009	0.0021	0.0013
Uncertainties from fragmentation fractions	0.0000	0.0000	0.0000	0.0058
Difference between ${\omega }_i$ or ${\omega }_i^{\mathrm{data}}$	0.0003	0.0003	0.0003	0.0003
Neglecting term with $A_P\left({\Xi }_b\right)$ in Eq. 3	0.0000	0.0000	0.0000	0.0071
Validity of $N_b=N_{\overline{b}}$ in each bin	0.0000	0.0000	0.0000	0.0032
$A_{CP}\left(B^{+}\to J/\psi K^{+}\right)$	0.0028	0.0000	0.0000	0.0028
$A_D(\overline{K}{}_{}{}^0)$	0.0001	0.0000	0.0000	0.0002
Total systematic uncertainty	0.0037	0.0011	0.0059	0.0108
	Uncertainty [$\sqrt{s}$ = 8$\mathrm{TeV}$]
Source	$A_P\left(B^{+}\right)$	$A_P\left(B^0\right)$	$A_P\left(B_s^0\right)$	$A_P\left({\Lambda }_b^0\right)$
Signal mass shape	0.0006	0.0004	0.0035	0.0021
Decay-time bias	0.0000	0.0000	0.0008	0.0004
$\Delta m_d$, $\Delta m_s$	0.0000	0.0001	0.0015	0.0008
Decay-time resolution	0.0000	0.0000	0.0028	0.0016
Final-state radiation	0.0000	0.0001	0.0001	0.0001
Decay-time reconstruction efficiency	0.0000	0.0001	0.0001	0.0001
Combinatorial background mass shape	0.0002	0.0000	0.0004	0.0003
Partially reconstructed background mass shape	0.0000	0.0000	0.0027	0.0015
$\Delta {\Gamma }_s$	0.0000	0.0000	0.0001	0.0001
$A_D\left(K^{+}\right)$	0.0014	0.0000	0.0000	0.0011
${\left|q/p\right|}_{B^0}$, ${\left|q/p\right|}_{B_s^0}$	0.0000	0.0009	0.0021	0.0014
Uncertainties from fragmentation fractions	0.0000	0.0000	0.0000	0.0025
Difference between ${\omega }_i$ or ${\omega }_i^{\mathrm{data}}$	0.0002	0.0003	0.0003	0.0003
Neglecting term with $A_P\left({\Xi }_b\right)$ in Eq. 3	0.0000	0.0000	0.0000	0.0046
Validity of $N_b=N_{\overline{b}}$ in each bin	0.0000	0.0000	0.0000	0.0033
$A_{CP}\left(B^{+}\to J/\psi K^{+}\right)$	0.0028	0.0000	0.0000	0.0027
$A_D(\overline{K}{}_{}{}^0)$	0.0001	0.0000	0.0000	0.0002
Total systematic uncertainty	0.0032	0.0010	0.0059	0.0076

Table 5. Table 5: Values for m 𝑚 m and q 𝑞 q and their correlation coefficient ( ρ 𝜌 \rho ) obtained from fits to the values reported in Tables  10 – 17 with a first order polynomial function (FOPF), A P ​ ( b ​ hadron ) = a ​ x + b subscript 𝐴 P 𝑏 hadron 𝑎 𝑥 𝑏 A_{\mathrm{P}}(b\,\mathrm{hadron})=a\,x+b with x = p T , y 𝑥 subscript 𝑝 T 𝑦 x=\mbox{$p_{\rm T}$},y . The label SL indicates the fit to the values with a straight line.

$p_T$	$\sqrt{s}$ = 7$\mathrm{TeV}$
	$B^{+}$	$B^0$	$B_s^0$	${\Lambda }_b^0$
$a\left[c/\mathrm{GeV}\cdot {10}^{-4}\right]$	$-3\pm 6$	$7\pm 14$	$-33\pm 26$	$10\pm 60$
$b\left[{10}^{-3}\right]$	$-1\pm 5$	$-3\pm 12$	$44\pm 37$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}7}\pm 50$
$\rho (m,q)$	$-0.59$	$-0.78$	$-0.89$	$-0.85$
fit ${\chi }^2/\mathrm{ndf}$ (SL)	1.05	1.67	1.67	2.03
fit ${\chi }^2/\mathrm{ndf}$ (FOPF)	0.95	1.50	1.68	1.34
$p_T$	$\sqrt{s}$ = 8$\mathrm{TeV}$
	$B^{+}$	$B^0$	$B_s^0$	${\Lambda }_b^0$
$a\left[c/\mathrm{GeV}\cdot {10}^{-4}\right]$	$-5\pm 4$	$12\pm 9$	$2\pm 18$	$30\pm 30$
$b\left[{10}^{-3}\right]$	$-5\pm 4$	$-18\pm 8$	$8\pm 25$	$12\pm 32$
$\rho (m,q)$	$-0.49$	$-0.78$	$-0.89$	$-0.84$
fit ${\chi }^2/\mathrm{ndf}$ (SL)	0.99	0.54	1.80	0.07
fit ${\chi }^2/\mathrm{ndf}$ (FOPF)	1.12	0.67	1.19	0.28
$y$	$\sqrt{s}$ = 7$\mathrm{TeV}$
	$B^{+}$	$B^0$	$B_s^0$	${\Lambda }_b^0$
$a\left[{10}^{-4}\right]$	$22\pm 45$	$-71\pm 141$	$-542\pm 469$	$810\pm 420$
$b\left[{10}^{-3}\right]$	$-12\pm 14$	$-12\pm 46$	$165\pm 153$	$-240\pm 130$
$\rho (m,q)$	$-0.96$	$-0.99$	$-0.99$	$-0.98$
fit ${\chi }^2/\mathrm{ndf}$ (SL)	1.36	2.79	0.48	0.86
fit ${\chi }^2/\mathrm{ndf}$ (FOPF)	1.17	2.36	0.91	2.29
$y$	$\sqrt{s}$ = 8$\mathrm{TeV}$
	$B^{+}$	$B^0$	$B_s^0$	${\Lambda }_b^0$
$a\left[{10}^{-4}\right]$	$-86\pm 29$	$-44\pm 100$	$-217\pm 321$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}470}\pm 280$
$b\left[{10}^{-3}\right]$	$19\pm 9$	$-4\pm 32$	$85\pm 105$	$-111\pm 90$
$\rho (m,q)$	$-0.93$	$-0.99$	$-0.98$	$-0.98$
fit ${\chi }^2/\mathrm{ndf}$ (SL)	1.10	0.86	0.42	1.35
fit ${\chi }^2/\mathrm{ndf}$ (FOPF)	2.43	0.75	0.44	2.23

Table 6. Table 6: Values of A P ​ ( B + ) subscript 𝐴 P superscript 𝐵 A_{\rm P}({{{B}^{+}}}) and A P ​ ( B 0 ) subscript 𝐴 P superscript 𝐵 0 A_{\rm P}({{B}^{0}}) in each kinematic bin for data collected in proton-proton collisions at centre-of-mass energy of 7 TeV TeV \mathrm{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second systematic.

$p_T$ [$\mathrm{GeV}/c$ ]	$y$	$A_P{\left(B^{+}\right)}_{\sqrt{s}=7\mathrm{TeV}}$	$A_P{\left(B^0\right)}_{\sqrt{s}=7\mathrm{TeV}}$
$(0.00,2.00)$	$(2.10,2.70)$	$0.0085\pm 0.0156\pm 0.0036$	$0.0722\pm 0.0770\pm 0.0010$
$(0.00,2.00)$	$(2.70,2.85)$	$-0.0014\pm 0.0191\pm 0.0036$	$-0.1108\pm 0.0815\pm 0.0020$
$(0.00,2.00)$	$(2.85,3.00)$	$0.0016\pm 0.0177\pm 0.0036$	$-0.0300\pm 0.0733\pm 0.0024$
$(0.00,2.00)$	$(3.00,3.15)$	$-0.0052\pm 0.0171\pm 0.0036$	$-0.0849\pm 0.0624\pm 0.0038$
$(0.00,2.00)$	$(3.15,3.30)$	$-0.0006\pm 0.0171\pm 0.0037$	$-0.0662\pm 0.0638\pm 0.0035$
$(0.00,2.00)$	$(3.30,3.70)$	$0.0107\pm 0.0110\pm 0.0040$	$0.0116\pm 0.0397\pm 0.0011$
$(0.00,2.00)$	$(3.70,4.50)$	$-0.0104\pm 0.0141\pm 0.0046$	$0.0702\pm 0.0462\pm 0.0013$
$(2.00,4.50)$	$(2.10,2.70)$	$0.0007\pm 0.0088\pm 0.0036$	$0.0691\pm 0.0392\pm 0.0017$
$(2.00,4.50)$	$(2.70,2.85)$	$-0.0171\pm 0.0112\pm 0.0036$	$0.0136\pm 0.0409\pm 0.0013$
$(2.00,4.50)$	$(2.85,3.00)$	$-0.0120\pm 0.0105\pm 0.0036$	$-0.0284\pm 0.0375\pm 0.0010$
$(2.00,4.50)$	$(3.00,3.15)$	$-0.0269\pm 0.0101\pm 0.0037$	$-0.0273\pm 0.0360\pm 0.0009$
$(2.00,4.50)$	$(3.15,3.30)$	$0.0043\pm 0.0102\pm 0.0038$	$0.0137\pm 0.0351\pm 0.0015$
$(2.00,4.50)$	$(3.30,3.70)$	$-0.0167\pm 0.0071\pm 0.0041$	$-0.0273\pm 0.0230\pm 0.0028$
$(2.00,4.50)$	$(3.70,4.50)$	$0.0053\pm 0.0098\pm 0.0045$	$-0.0269\pm 0.0279\pm 0.0013$
$(4.50,7.00)$	$(2.10,2.70)$	$0.0023\pm 0.0087\pm 0.0035$	$0.0597\pm 0.0329\pm 0.0039$
$(4.50,7.00)$	$(2.70,2.85)$	$-0.0002\pm 0.0120\pm 0.0037$	$-0.0177\pm 0.0404\pm 0.0010$
$(4.50,7.00)$	$(2.85,3.00)$	$0.0034\pm 0.0116\pm 0.0038$	$-0.0103\pm 0.0362\pm 0.0019$
$(4.50,7.00)$	$(3.00,3.15)$	$0.0092\pm 0.0115\pm 0.0039$	$-0.0696\pm 0.0372\pm 0.0016$
$(4.50,7.00)$	$(3.15,3.30)$	$-0.0092\pm 0.0120\pm 0.0042$	$-0.0444\pm 0.0359\pm 0.0015$
$(4.50,7.00)$	$(3.30,3.70)$	$-0.0168\pm 0.0088\pm 0.0044$	$-0.0214\pm 0.0234\pm 0.0010$
$(4.50,7.00)$	$(3.70,4.50)$	$0.0010\pm 0.0129\pm 0.0044$	$0.0192\pm 0.0316\pm 0.0013$
$(7.00,8.25)$	$(2.10,2.70)$	$0.0031\pm 0.0140\pm 0.0036$	$-0.0239\pm 0.0441\pm 0.0025$
$(7.00,8.25)$	$(2.70,2.85)$	$-0.0591\pm 0.0208\pm 0.0039$	$-0.2197\pm 0.0602\pm 0.0017$
$(7.00,8.25)$	$(2.85,3.00)$	$-0.0089\pm 0.0203\pm 0.0040$	$-0.0619\pm 0.0595\pm 0.0031$
$(7.00,8.25)$	$(3.00,3.15)$	$0.0016\pm 0.0213\pm 0.0043$	$-0.0151\pm 0.0590\pm 0.0047$
$(7.00,8.25)$	$(3.15,3.30)$	$-0.0205\pm 0.0222\pm 0.0044$	$-0.0037\pm 0.0566\pm 0.0039$
$(7.00,8.25)$	$(3.30,3.70)$	$0.0303\pm 0.0172\pm 0.0046$	$-0.0305\pm 0.0406\pm 0.0018$
$(7.00,8.25)$	$(3.70,4.50)$	$0.0603\pm 0.0259\pm 0.0047$	$-0.0348\pm 0.0516\pm 0.0018$
$(8.25,9.50)$	$(2.10,2.70)$	$-0.0134\pm 0.0157\pm 0.0037$	$-0.0442\pm 0.0477\pm 0.0087$
$(8.25,9.50)$	$(2.70,2.85)$	$-0.0099\pm 0.0246\pm 0.0039$	$-0.0506\pm 0.0652\pm 0.0028$
$(8.25,9.50)$	$(2.85,3.00)$	$-0.0112\pm 0.0246\pm 0.0042$	$-0.0611\pm 0.0674\pm 0.0043$
$(8.25,9.50)$	$(3.00,3.15)$	$-0.0613\pm 0.0251\pm 0.0044$	$-0.0015\pm 0.0695\pm 0.0024$
$(8.25,9.50)$	$(3.15,3.30)$	$0.0552\pm 0.0279\pm 0.0045$	$0.0219\pm 0.0731\pm 0.0014$
$(8.25,9.50)$	$(3.30,3.70)$	$-0.0038\pm 0.0216\pm 0.0046$	$-0.0621\pm 0.0478\pm 0.0032$
$(8.25,9.50)$	$(3.70,4.50)$	$0.0047\pm 0.0342\pm 0.0047$	$-0.0856\pm 0.0637\pm 0.0025$
$(9.50,10.75)$	$(2.10,2.70)$	$-0.0249\pm 0.0182\pm 0.0037$	$0.0408\pm 0.0525\pm 0.0089$
$(9.50,10.75)$	$(2.70,2.85)$	$-0.0113\pm 0.0292\pm 0.0041$	$0.0228\pm 0.0759\pm 0.0036$
$(9.50,10.75)$	$(2.85,3.00)$	$-0.0241\pm 0.0290\pm 0.0045$	$0.0102\pm 0.0904\pm 0.0017$
$(9.50,10.75)$	$(3.00,3.15)$	$0.0267\pm 0.0318\pm 0.0045$	$-0.0586\pm 0.0847\pm 0.0023$
$p_T$ [$\mathrm{GeV}/c$ ]	$y$	$A_P{\left(B^{+}\right)}_{\sqrt{s}=7\mathrm{TeV}}$	$A_P{\left(B^0\right)}_{\sqrt{s}=7\mathrm{TeV}}$
$(9.50,10.75)$	$(3.15,3.30)$	$0.0118\pm 0.0352\pm 0.0048$	$-0.0577\pm 0.0775\pm 0.0012$
$(9.50,10.75)$	$(3.30,3.70)$	$-0.0164\pm 0.0281\pm 0.0048$	$0.0624\pm 0.0577\pm 0.0016$
$(9.50,10.75)$	$(3.70,4.50)$	$-0.0605\pm 0.0411\pm 0.0049$	$-0.0328\pm 0.0946\pm 0.0021$
$(10.75,12.00)$	$(2.10,2.70)$	$-0.0200\pm 0.0206\pm 0.0038$	$0.0154\pm 0.0636\pm 0.0023$
$(10.75,12.00)$	$(2.70,2.85)$	$-0.0068\pm 0.0344\pm 0.0044$	$-0.0104\pm 0.1017\pm 0.0044$
$(10.75,12.00)$	$(2.85,3.00)$	$-0.0017\pm 0.0362\pm 0.0045$	$0.0179\pm 0.0849\pm 0.0040$
$(10.75,12.00)$	$(3.00,3.15)$	$-0.0181\pm 0.0411\pm 0.0047$	$0.1481\pm 0.0890\pm 0.0024$
$(10.75,12.00)$	$(3.15,3.30)$	$-0.0239\pm 0.0441\pm 0.0047$	$0.0478\pm 0.0835\pm 0.0025$
$(10.75,12.00)$	$(3.30,3.70)$	$0.0058\pm 0.0362\pm 0.0048$	$0.0377\pm 0.0731\pm 0.0037$
$(10.75,12.00)$	$(3.70,4.50)$	$0.0485\pm 0.0547\pm 0.0051$	$0.1058\pm 0.1181\pm 0.0018$
$(12.00,15.00)$	$(2.10,2.70)$	$0.0059\pm 0.0174\pm 0.0039$	$-0.0071\pm 0.0446\pm 0.0039$
$(12.00,15.00)$	$(2.70,2.85)$	$0.0210\pm 0.0321\pm 0.0046$	$0.0264\pm 0.0924\pm 0.0042$
$(12.00,15.00)$	$(2.85,3.00)$	$0.0092\pm 0.0334\pm 0.0062$	$0.0230\pm 0.0775\pm 0.0046$
$(12.00,15.00)$	$(3.00,3.15)$	$-0.0267\pm 0.0386\pm 0.0050$	$-0.1190\pm 0.0791\pm 0.0040$
$(12.00,15.00)$	$(3.15,3.30)$	$-0.0516\pm 0.0420\pm 0.0046$	$0.1330\pm 0.0909\pm 0.0029$
$(12.00,15.00)$	$(3.30,3.70)$	$0.0071\pm 0.0349\pm 0.0052$	$0.0469\pm 0.0588\pm 0.0021$
$(12.00,15.00)$	$(3.70,4.50)$	$0.0748\pm 0.0542\pm 0.0049$	$-0.1026\pm 0.0854\pm 0.0031$
$(15.00,30.00)$	$(2.10,2.70)$	$0.0116\pm 0.0188\pm 0.0040$	$0.0703\pm 0.0456\pm 0.0014$
$(15.00,30.00)$	$(2.70,2.85)$	$-0.0763\pm 0.0401\pm 0.0046$	$-0.0009\pm 0.0748\pm 0.0034$
$(15.00,30.00)$	$(2.85,3.00)$	$-0.0541\pm 0.0458\pm 0.0047$	$-0.0550\pm 0.0755\pm 0.0049$
$(15.00,30.00)$	$(3.00,3.15)$	$-0.0449\pm 0.0512\pm 0.0046$	$-0.1637\pm 0.0925\pm 0.0026$
$(15.00,30.00)$	$(3.15,3.30)$	$0.0011\pm 0.0599\pm 0.0073$	$0.0456\pm 0.1119\pm 0.0018$
$(15.00,30.00)$	$(3.30,3.70)$	$0.0089\pm 0.0502\pm 0.0048$	$-0.0193\pm 0.0777\pm 0.0027$
$(15.00,30.00)$	$(3.70,4.50)$	$-0.0662\pm 0.0827\pm 0.0186$	$0.1690\pm 0.1332\pm 0.0030$

Table 7. Table 7: Values of A P ​ ( B + ) subscript 𝐴 P superscript 𝐵 A_{\rm P}({{{B}^{+}}}) and A P ​ ( B 0 ) subscript 𝐴 P superscript 𝐵 0 A_{\rm P}({{B}^{0}}) in each kinematic bin for data collected in proton-proton collisions at centre-of-mass energy of 8 TeV TeV \mathrm{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second systematic.

$p_T$ [$\mathrm{GeV}/c$ ]	$y$	$A_P{\left(B^{+}\right)}_{\sqrt{s}=8\mathrm{TeV}}$	$A_P{\left(B^0\right)}_{\sqrt{s}=8\mathrm{TeV}}$
$(0.00,2.00)$	$(2.10,2.70)$	$-0.0178\pm 0.0097\pm 0.0031$	$0.0068\pm 0.0537\pm 0.0009$
$(0.00,2.00)$	$(2.70,2.85)$	$-0.0027\pm 0.0126\pm 0.0031$	$-0.0735\pm 0.0719\pm 0.0017$
$(0.00,2.00)$	$(2.85,3.00)$	$0.0093\pm 0.0120\pm 0.0031$	$0.0503\pm 0.0628\pm 0.0011$
$(0.00,2.00)$	$(3.00,3.15)$	$0.0005\pm 0.0119\pm 0.0031$	$0.0086\pm 0.0549\pm 0.0034$
$(0.00,2.00)$	$(3.15,3.30)$	$-0.0230\pm 0.0119\pm 0.0033$	$0.0817\pm 0.0617\pm 0.0016$
$(0.00,2.00)$	$(3.30,3.70)$	$-0.0120\pm 0.0080\pm 0.0033$	$0.0668\pm 0.0367\pm 0.0009$
$(0.00,2.00)$	$(3.70,4.50)$	$-0.0077\pm 0.0103\pm 0.0037$	$-0.0419\pm 0.0453\pm 0.0010$
$(2.00,4.50)$	$(2.10,2.70)$	$0.0050\pm 0.0054\pm 0.0031$	$-0.0192\pm 0.0234\pm 0.0013$
$(2.00,4.50)$	$(2.70,2.85)$	$-0.0076\pm 0.0073\pm 0.0031$	$-0.0070\pm 0.0291\pm 0.0009$
$(2.00,4.50)$	$(2.85,3.00)$	$0.0009\pm 0.0070\pm 0.0031$	$-0.0088\pm 0.0278\pm 0.0009$
$(2.00,4.50)$	$(3.00,3.15)$	$-0.0046\pm 0.0069\pm 0.0032$	$-0.0213\pm 0.0271\pm 0.0009$
$(2.00,4.50)$	$(3.15,3.30)$	$-0.0018\pm 0.0070\pm 0.0032$	$-0.0635\pm 0.0260\pm 0.0012$
$(2.00,4.50)$	$(3.30,3.70)$	$-0.0081\pm 0.0049\pm 0.0034$	$-0.0169\pm 0.0174\pm 0.0009$
$(2.00,4.50)$	$(3.70,4.50)$	$-0.0133\pm 0.0067\pm 0.0036$	$-0.0131\pm 0.0203\pm 0.0009$
$(4.50,7.00)$	$(2.10,2.70)$	$-0.0045\pm 0.0054\pm 0.0031$	$-0.0074\pm 0.0192\pm 0.0028$
$(4.50,7.00)$	$(2.70,2.85)$	$-0.0002\pm 0.0077\pm 0.0031$	$-0.0440\pm 0.0264\pm 0.0027$
$(4.50,7.00)$	$(2.85,3.00)$	$-0.0019\pm 0.0075\pm 0.0032$	$0.0315\pm 0.0235\pm 0.0028$
$(4.50,7.00)$	$(3.00,3.15)$	$-0.0107\pm 0.0076\pm 0.0033$	$-0.0203\pm 0.0233\pm 0.0020$
$(4.50,7.00)$	$(3.15,3.30)$	$-0.0175\pm 0.0078\pm 0.0034$	$-0.0248\pm 0.0234\pm 0.0010$
$(4.50,7.00)$	$(3.30,3.70)$	$-0.0241\pm 0.0059\pm 0.0035$	$-0.0254\pm 0.0159\pm 0.0010$
$(4.50,7.00)$	$(3.70,4.50)$	$-0.0101\pm 0.0087\pm 0.0036$	$-0.0015\pm 0.0213\pm 0.0010$
$(7.00,8.25)$	$(2.10,2.70)$	$-0.0052\pm 0.0086\pm 0.0031$	$0.0080\pm 0.0276\pm 0.0028$
$(7.00,8.25)$	$(2.70,2.85)$	$-0.0177\pm 0.0131\pm 0.0033$	$-0.0383\pm 0.0390\pm 0.0014$
$(7.00,8.25)$	$(2.85,3.00)$	$-0.0083\pm 0.0132\pm 0.0033$	$-0.0543\pm 0.0382\pm 0.0025$
$(7.00,8.25)$	$(3.00,3.15)$	$0.0065\pm 0.0134\pm 0.0035$	$-0.0575\pm 0.0377\pm 0.0012$
$(7.00,8.25)$	$(3.15,3.30)$	$-0.0055\pm 0.0144\pm 0.0040$	$-0.0120\pm 0.0379\pm 0.0013$
$(7.00,8.25)$	$(3.30,3.70)$	$-0.0003\pm 0.0111\pm 0.0036$	$-0.0089\pm 0.0268\pm 0.0044$
$(7.00,8.25)$	$(3.70,4.50)$	$-0.0300\pm 0.0168\pm 0.0036$	$-0.0486\pm 0.0364\pm 0.0022$
$(8.25,9.50)$	$(2.10,2.70)$	$-0.0038\pm 0.0097\pm 0.0031$	$-0.0215\pm 0.0286\pm 0.0017$
$(8.25,9.50)$	$(2.70,2.85)$	$-0.0070\pm 0.0153\pm 0.0033$	$0.0710\pm 0.0415\pm 0.0013$
$(8.25,9.50)$	$(2.85,3.00)$	$-0.0228\pm 0.0157\pm 0.0034$	$0.0123\pm 0.0395\pm 0.0010$
$(8.25,9.50)$	$(3.00,3.15)$	$-0.0236\pm 0.0164\pm 0.0037$	$0.0747\pm 0.0411\pm 0.0023$
$(8.25,9.50)$	$(3.15,3.30)$	$-0.0252\pm 0.0182\pm 0.0042$	$-0.0533\pm 0.0459\pm 0.0025$
$(8.25,9.50)$	$(3.30,3.70)$	$-0.0036\pm 0.0141\pm 0.0037$	$0.0152\pm 0.0299\pm 0.0009$
$(8.25,9.50)$	$(3.70,4.50)$	$-0.0293\pm 0.0220\pm 0.0037$	$-0.0063\pm 0.0448\pm 0.0034$
$(9.50,10.75)$	$(2.10,2.70)$	$0.0060\pm 0.0109\pm 0.0032$	$0.0022\pm 0.0324\pm 0.0022$
$(9.50,10.75)$	$(2.70,2.85)$	$-0.0011\pm 0.0183\pm 0.0036$	$0.0429\pm 0.0491\pm 0.0050$
$(9.50,10.75)$	$(2.85,3.00)$	$0.0122\pm 0.0182\pm 0.0036$	$0.0513\pm 0.0509\pm 0.0021$
$(9.50,10.75)$	$(3.00,3.15)$	$0.0067\pm 0.0204\pm 0.0037$	$-0.0898\pm 0.0499\pm 0.0059$
$(9.50,10.75)$	$(3.15,3.30)$	$-0.0462\pm 0.0233\pm 0.0037$	$-0.0220\pm 0.0494\pm 0.0034$
$(9.50,10.75)$	$(3.30,3.70)$	$-0.0290\pm 0.0181\pm 0.0037$	$-0.0204\pm 0.0353\pm 0.0013$
$(9.50,10.75)$	$(3.70,4.50)$	$-0.0243\pm 0.0273\pm 0.0037$	$-0.0849\pm 0.0509\pm 0.0026$
$p_T$ [$\mathrm{GeV}/c$ ]	$y$	$A_P{\left(B^{+}\right)}_{\sqrt{s}=8\mathrm{TeV}}$	$A_P{\left(B^0\right)}_{\sqrt{s}=8\mathrm{TeV}}$
$(10.75,12.00)$	$(2.10,2.70)$	$0.0191\pm 0.0128\pm 0.0032$	$0.0034\pm 0.0355\pm 0.0056$
$(10.75,12.00)$	$(2.70,2.85)$	$-0.0562\pm 0.0220\pm 0.0034$	$-0.0193\pm 0.0593\pm 0.0026$
$(10.75,12.00)$	$(2.85,3.00)$	$0.0172\pm 0.0233\pm 0.0037$	$0.0198\pm 0.0628\pm 0.0066$
$(10.75,12.00)$	$(3.00,3.15)$	$-0.0080\pm 0.0262\pm 0.0044$	$-0.0056\pm 0.0565\pm 0.0012$
$(10.75,12.00)$	$(3.15,3.30)$	$0.0162\pm 0.0282\pm 0.0038$	$-0.0638\pm 0.0582\pm 0.0040$
$(10.75,12.00)$	$(3.30,3.70)$	$-0.0393\pm 0.0233\pm 0.0037$	$0.0205\pm 0.0454\pm 0.0083$
$(10.75,12.00)$	$(3.70,4.50)$	$0.0317\pm 0.0353\pm 0.0038$	$0.0139\pm 0.0709\pm 0.0009$
$(12.00,15.00)$	$(2.10,2.70)$	$0.0067\pm 0.0106\pm 0.0032$	$-0.0364\pm 0.0278\pm 0.0010$
$(12.00,15.00)$	$(2.70,2.85)$	$-0.0232\pm 0.0195\pm 0.0035$	$-0.0007\pm 0.0525\pm 0.0026$
$(12.00,15.00)$	$(2.85,3.00)$	$0.0171\pm 0.0211\pm 0.0047$	$0.0255\pm 0.0467\pm 0.0010$
$(12.00,15.00)$	$(3.00,3.15)$	$0.0065\pm 0.0241\pm 0.0046$	$0.0080\pm 0.0521\pm 0.0017$
$(12.00,15.00)$	$(3.15,3.30)$	$-0.0101\pm 0.0273\pm 0.0038$	$-0.0019\pm 0.0491\pm 0.0021$
$(12.00,15.00)$	$(3.30,3.70)$	$-0.0214\pm 0.0219\pm 0.0039$	$-0.0526\pm 0.0373\pm 0.0045$
$(12.00,15.00)$	$(3.70,4.50)$	$-0.0511\pm 0.0340\pm 0.0038$	$-0.0494\pm 0.0605\pm 0.0027$
$(15.00,30.00)$	$(2.10,2.70)$	$-0.0203\pm 0.0115\pm 0.0033$	$0.0217\pm 0.0267\pm 0.0012$
$(15.00,30.00)$	$(2.70,2.85)$	$-0.0340\pm 0.0252\pm 0.0036$	$-0.0204\pm 0.0491\pm 0.0038$
$(15.00,30.00)$	$(2.85,3.00)$	$-0.0231\pm 0.0277\pm 0.0055$	$0.0878\pm 0.0520\pm 0.0020$
$(15.00,30.00)$	$(3.00,3.15)$	$0.0347\pm 0.0317\pm 0.0037$	$0.0120\pm 0.0534\pm 0.0016$
$(15.00,30.00)$	$(3.15,3.30)$	$-0.0064\pm 0.0379\pm 0.0068$	$0.0153\pm 0.0626\pm 0.0025$
$(15.00,30.00)$	$(3.30,3.70)$	$-0.0221\pm 0.0311\pm 0.0042$	$-0.0647\pm 0.0434\pm 0.0013$
$(15.00,30.00)$	$(3.70,4.50)$	$-0.0987\pm 0.0496\pm 0.0063$	$0.0394\pm 0.0777\pm 0.0042$

Table 8. Table 8: Values of A P ​ ( B s 0 ) subscript 𝐴 P subscript superscript 𝐵 0 𝑠 A_{\rm P}({{B}^{0}_{s}}) and A P ​ ( Λ b 0 ) subscript 𝐴 P subscript superscript Λ 0 𝑏 A_{\rm P}({{\mathchar 28931\relax}^{0}_{b}}) in each kinematic bin for data collected in proton-proton collisions at centre-of-mass energy of 7 TeV TeV \mathrm{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second systematic.

$p_T$ [$\mathrm{GeV}/c$ ]	$y$	$A_P{\left(B_s^0\right)}_{\sqrt{s}=7\mathrm{TeV}}$	$A_P{\left({\Lambda }_b^0\right)}_{\sqrt{s}=7\mathrm{TeV}}$
$(2.00,7.00)$	$(2.10,3.00)$	$0.0166\pm 0.0632\pm 0.0125$	$-0.0892\pm 0.0508\pm 0.0214$
$(2.00,7.00)$	$(3.00,3.30)$	$0.0311\pm 0.0773\pm 0.0151$	$0.0507\pm 0.0539\pm 0.0208$
$(2.00,7.00)$	$(3.30,4.50)$	$-0.0833\pm 0.0558\pm 0.0132$	$0.0849\pm 0.0401\pm 0.0188$
$(7.00,9.50)$	$(2.10,3.00)$	$0.0364\pm 0.0479\pm 0.0068$	$0.1374\pm 0.0697\pm 0.0313$
$(7.00,9.50)$	$(3.00,3.30)$	$0.0206\pm 0.0682\pm 0.0127$	$0.0138\pm 0.0913\pm 0.0298$
$(7.00,9.50)$	$(3.30,4.50)$	$0.0058\pm 0.0584\pm 0.0089$	$0.0466\pm 0.0770\pm 0.0347$
$(9.50,12.00)$	$(2.10,3.00)$	$-0.0039\pm 0.0456\pm 0.0121$	$-0.0128\pm 0.0985\pm 0.0367$
$(9.50,12.00)$	$(3.00,3.30)$	$0.1095\pm 0.0723\pm 0.0179$	$-0.0848\pm 0.1379\pm 0.0452$
$(9.50,12.00)$	$(3.30,4.50)$	$0.1539\pm 0.0722\pm 0.0212$	$-0.1523\pm 0.1414\pm 0.0488$
$(12.00,30.00)$	$(2.10,3.00)$	$-0.0271\pm 0.0336\pm 0.0061$	$-0.0720\pm 0.1248\pm 0.0465$
$(12.00,30.00)$	$(3.00,3.30)$	$-0.0542\pm 0.0612\pm 0.0106$	$0.3291\pm 0.2299\pm 0.0918$
$(12.00,30.00)$	$(3.30,4.50)$	$-0.0586\pm 0.0648\pm 0.0150$	$-0.0571\pm 0.2162\pm 0.0800$

Table 9. Table 9: Values of A P ​ ( B s 0 ) subscript 𝐴 P subscript superscript 𝐵 0 𝑠 A_{\rm P}({{B}^{0}_{s}}) and A P ​ ( Λ b 0 ) subscript 𝐴 P subscript superscript Λ 0 𝑏 A_{\rm P}({{\mathchar 28931\relax}^{0}_{b}}) in each kinematic bin for data collected in proton-proton collisions at centre-of-mass energy of 8 TeV TeV \mathrm{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second systematic.

$p_T$ [$\mathrm{GeV}/c$ ]	$y$	$A_P{\left(B_s^0\right)}_{\sqrt{s}=8\mathrm{TeV}}$	$A_P{\left({\Lambda }_b^0\right)}_{\sqrt{s}=8\mathrm{TeV}}$
$(2.00,7.00)$	$(2.10,3.00)$	$0.0412\pm 0.0416\pm 0.0150$	$0.0032\pm 0.0318\pm 0.0139$
$(2.00,7.00)$	$(3.00,3.30)$	$-0.0241\pm 0.0574\pm 0.0079$	$0.0929\pm 0.0392\pm 0.0171$
$(2.00,7.00)$	$(3.30,4.50)$	$0.0166\pm 0.0391\pm 0.0092$	$0.0437\pm 0.0284\pm 0.0173$
$(7.00,9.50)$	$(2.10,3.00)$	$0.0482\pm 0.0320\pm 0.0067$	$0.0069\pm 0.0434\pm 0.0169$
$(7.00,9.50)$	$(3.00,3.30)$	$0.0983\pm 0.0470\pm 0.0155$	$0.0076\pm 0.0589\pm 0.0259$
$(7.00,9.50)$	$(3.30,4.50)$	$-0.0430\pm 0.0386\pm 0.0079$	$0.1053\pm 0.0524\pm 0.0252$
$(9.50,12.00)$	$(2.10,3.00)$	$0.0067\pm 0.0303\pm 0.0063$	$-0.0512\pm 0.0594\pm 0.0215$
$(9.50,12.00)$	$(3.00,3.30)$	$-0.1283\pm 0.0503\pm 0.0171$	$0.2355\pm 0.0877\pm 0.0399$
$(9.50,12.00)$	$(3.30,4.50)$	$-0.0500\pm 0.0460\pm 0.0104$	$0.1531\pm 0.0838\pm 0.0320$
$(12.00,30.00)$	$(2.10,3.00)$	$-0.0012\pm 0.0222\pm 0.0050$	$0.0453\pm 0.0762\pm 0.0300$
$(12.00,30.00)$	$(3.00,3.30)$	$0.0421\pm 0.0416\pm 0.0162$	$-0.0934\pm 0.1377\pm 0.0493$
$(12.00,30.00)$	$(3.30,4.50)$	$0.0537\pm 0.0447\pm 0.0124$	$0.3173\pm 0.1411\pm 0.0655$

Table 10. Table 10: Values of the production asymmetries in bins of p T subscript 𝑝 T p_{\rm T} , integrated over y 𝑦 y , for B + superscript 𝐵 {{B}^{+}} and B 0 superscript 𝐵 0 {B}^{0} mesons for data collected in proton-proton collisions at the centre-of-mass energy of 7 TeV TeV \mathrm{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second systematic. The uncertainties among the bins are correlated due to the external inputs: A C ​ P ​ ( B + → J / ψ ​ K + ) subscript 𝐴 𝐶 𝑃 → superscript 𝐵 𝐽 𝜓 superscript 𝐾 A_{{C\!P}}({{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}) and A D ( K ¯ ) 0 A_{\rm D}({{\kern 1.99997pt\overline{\kern-1.99997ptK}{}}{}^{0}}) for A P ​ ( B + ) subscript 𝐴 P superscript 𝐵 A_{\rm P}({{{B}^{+}}}) , and | q / p | 𝑞 𝑝 |q/p| for A P ​ ( B 0 ) subscript 𝐴 P superscript 𝐵 0 A_{\rm P}({{B}^{0}}) .

$p_T$ [GeV/$c$]	$A_P{\left(B^{+}\right)}_{\sqrt{s}=7\mathrm{TeV}}$	$A_P{\left(B^0\right)}_{\sqrt{s}=7\mathrm{TeV}}$
$(0.00,2.00)$	$0.0015\pm 0.0067\pm 0.0036$	$0.0215\pm 0.0297\pm 0.0025$
$(2.00,4.50)$	$-0.0050\pm 0.0040\pm 0.0037$	$0.0123\pm 0.0163\pm 0.0078$
$(4.50,7.00)$	$-0.0010\pm 0.0045\pm 0.0038$	$0.0124\pm 0.0150\pm 0.0042$
$(7.00,8.25)$	$0.0083\pm 0.0080\pm 0.0041$	$-0.0440\pm 0.0219\pm 0.0012$
$(8.25,9.50)$	$-0.0078\pm 0.0096\pm 0.0039$	$-0.0476\pm 0.0248\pm 0.0038$
$(9.50,10.75)$	$-0.0220\pm 0.0114\pm 0.0044$	$0.0155\pm 0.0297\pm 0.0056$
$(10.75,12.00)$	$-0.0045\pm 0.0138\pm 0.0043$	$0.0404\pm 0.0357\pm 0.0040$
$(12.00,15.00)$	$0.0107\pm 0.0124\pm 0.0053$	$-0.0050\pm 0.0269\pm 0.0035$
$(15.00,30.00)$	$-0.0146\pm 0.0150\pm 0.0065$	$0.0333\pm 0.0298\pm 0.0077$

Table 11. Table 11: Values of the production asymmetries in bins of y 𝑦 y , integrated over p T subscript 𝑝 T p_{\rm T} , for B + superscript 𝐵 {{B}^{+}} and B 0 superscript 𝐵 0 {B}^{0} mesons for data collected in proton-proton collisions at the centre-of-mass energy of 7 TeV TeV \mathrm{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second systematic. The uncertainties among the bins are correlated due to the external inputs: A C ​ P ​ ( B + → J / ψ ​ K + ) subscript 𝐴 𝐶 𝑃 → superscript 𝐵 𝐽 𝜓 superscript 𝐾 A_{{C\!P}}({{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}) and A D ( K ¯ ) 0 A_{\rm D}({{\kern 1.99997pt\overline{\kern-1.99997ptK}{}}{}^{0}}) for A P ​ ( B + ) subscript 𝐴 P superscript 𝐵 A_{\rm P}({{{B}^{+}}}) , and | q / p | 𝑞 𝑝 |q/p| for A P ​ ( B 0 ) subscript 𝐴 P superscript 𝐵 0 A_{\rm P}({{B}^{0}}) .

$y$	$A_P{\left(B^{+}\right)}_{\sqrt{s}=7\mathrm{TeV}}$	$A_P{\left(B^0\right)}_{\sqrt{s}=7\mathrm{TeV}}$
$(2.10,2.70)$	$0.0007\pm 0.0047\pm 0.0036$	$0.0488\pm 0.0205\pm 0.0017$
$(2.70,2.85)$	$-0.0131\pm 0.0064\pm 0.0036$	$-0.0366\pm 0.0232\pm 0.0027$
$(2.85,3.00)$	$-0.0063\pm 0.0061\pm 0.0037$	$-0.0251\pm 0.0213\pm 0.0010$
$(3.00,3.15)$	$-0.0125\pm 0.0061\pm 0.0039$	$-0.0478\pm 0.0203\pm 0.0017$
$(3.15,3.30)$	$-0.0009\pm 0.0063\pm 0.0039$	$-0.0130\pm 0.0203\pm 0.0018$
$(3.30,3.70)$	$-0.0060\pm 0.0044\pm 0.0043$	$-0.0143\pm 0.0133\pm 0.0017$
$(3.70,4.50)$	$0.0041\pm 0.0062\pm 0.0046$	$0.0044\pm 0.0173\pm 0.0045$

Table 12. Table 12: Values of the production asymmetries in bins of p T subscript 𝑝 T p_{\rm T} , integrated over y 𝑦 y , for B + superscript 𝐵 {{B}^{+}} and B 0 superscript 𝐵 0 {B}^{0} mesons for data collected in proton-proton collisions at the centre-of-mass energy of 8 TeV TeV \mathrm{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second systematic. The uncertainties among the bins are correlated due to the external inputs: A C ​ P ​ ( B + → J / ψ ​ K + ) subscript 𝐴 𝐶 𝑃 → superscript 𝐵 𝐽 𝜓 superscript 𝐾 A_{{C\!P}}({{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}) and A D ( K ¯ ) 0 A_{\rm D}({{\kern 1.99997pt\overline{\kern-1.99997ptK}{}}{}^{0}}) for A P ​ ( B + ) subscript 𝐴 P superscript 𝐵 A_{\rm P}({{{B}^{+}}}) , and | q / p | 𝑞 𝑝 |q/p| for A P ​ ( B 0 ) subscript 𝐴 P superscript 𝐵 0 A_{\rm P}({{B}^{0}}) .

$p_T$ [GeV/$c$]	$A_P{\left(B^{+}\right)}_{\sqrt{s}=8\mathrm{TeV}}$	$A_P{\left(B^0\right)}_{\sqrt{s}=8\mathrm{TeV}}$
$(0.00,2.00)$	$-0.0105\pm 0.0045\pm 0.0031$	$0.0065\pm 0.0230\pm 0.0017$
$(2.00,4.50)$	$-0.0033\pm 0.0026\pm 0.0031$	$-0.0188\pm 0.0103\pm 0.0009$
$(4.50,7.00)$	$-0.0093\pm 0.0029\pm 0.0032$	$-0.0111\pm 0.0092\pm 0.0011$
$(7.00,8.25)$	$-0.0094\pm 0.0051\pm 0.0033$	$-0.0192\pm 0.0141\pm 0.0015$
$(8.25,9.50)$	$-0.0126\pm 0.0061\pm 0.0033$	$0.0015\pm 0.0155\pm 0.0009$
$(9.50,10.75)$	$-0.0073\pm 0.0073\pm 0.0034$	$-0.0156\pm 0.0177\pm 0.0013$
$(10.75,12.00)$	$0.0036\pm 0.0090\pm 0.0034$	$0.0017\pm 0.0210\pm 0.0027$
$(12.00,15.00)$	$-0.0082\pm 0.0079\pm 0.0035$	$-0.0270\pm 0.0171\pm 0.0009$
$(15.00,30.00)$	$-0.0251\pm 0.0095\pm 0.0040$	$0.0137\pm 0.0177\pm 0.0009$

Table 13. Table 13: Values of the production asymmetries in bins of y 𝑦 y , integrated over p T subscript 𝑝 T p_{\rm T} , for B + superscript 𝐵 {{B}^{+}} and B 0 superscript 𝐵 0 {B}^{0} mesons for data collected in proton-proton collisions at the centre-of-mass energy of 8 TeV TeV \mathrm{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second systematic. The uncertainties among the bins are correlated due to the external inputs: A C ​ P ​ ( B + → J / ψ ​ K + ) subscript 𝐴 𝐶 𝑃 → superscript 𝐵 𝐽 𝜓 superscript 𝐾 A_{{C\!P}}({{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}) and A D ( K ¯ ) 0 A_{\rm D}({{\kern 1.99997pt\overline{\kern-1.99997ptK}{}}{}^{0}}) for A P ​ ( B + ) subscript 𝐴 P superscript 𝐵 A_{\rm P}({{{B}^{+}}}) , and | q / p | 𝑞 𝑝 |q/p| for A P ​ ( B 0 ) subscript 𝐴 P superscript 𝐵 0 A_{\rm P}({{B}^{0}}) .

$y$	$A_P{\left(B^{+}\right)}_{\sqrt{s}=8\mathrm{TeV}}$	$A_P{\left(B^0\right)}_{\sqrt{s}=8\mathrm{TeV}}$
$(2.10,2.70)$	$-0.0023\pm 0.0029\pm 0.0031$	$-0.0082\pm 0.0128\pm 0.0012$
$(2.70,2.85)$	$-0.0080\pm 0.0041\pm 0.0031$	$-0.0237\pm 0.0173\pm 0.0009$
$(2.85,3.00)$	$0.0003\pm 0.0040\pm 0.0032$	$0.0148\pm 0.0159\pm 0.0015$
$(3.00,3.15)$	$-0.0038\pm 0.0040\pm 0.0032$	$-0.0140\pm 0.0151\pm 0.0009$
$(3.15,3.30)$	$-0.0123\pm 0.0042\pm 0.0034$	$-0.0193\pm 0.0158\pm 0.0021$
$(3.30,3.70)$	$-0.0138\pm 0.0030\pm 0.0034$	$-0.0029\pm 0.0103\pm 0.0010$
$(3.70,4.50)$	$-0.0144\pm 0.0042\pm 0.0037$	$-0.0201\pm 0.0137\pm 0.0010$

Table 14. Table 14: Values of the production asymmetries in bins of p T subscript 𝑝 T p_{\rm T} , integrated over y 𝑦 y , for the B s 0 subscript superscript 𝐵 0 𝑠 {B}^{0}_{s} meson and the Λ b 0 subscript superscript Λ 0 𝑏 {\mathchar 28931\relax}^{0}_{b} baryon for data collected in proton-proton collisions at the centre-of-mass energy of 7 TeV TeV \mathrm{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second systematic. The uncertainties among the bins are correlated due to the external inputs: A C ​ P ​ ( B + → J / ψ ​ K + ) subscript 𝐴 𝐶 𝑃 → superscript 𝐵 𝐽 𝜓 superscript 𝐾 A_{{C\!P}}({{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}) , A D ( K ¯ ) 0 A_{\rm D}({{\kern 1.99997pt\overline{\kern-1.99997ptK}{}}{}^{0}}) , | q / p | B 0 subscript 𝑞 𝑝 superscript 𝐵 0 |q/p|_{{{B}^{0}}} and | q / p | B s 0 subscript 𝑞 𝑝 subscript superscript 𝐵 0 𝑠 |q/p|_{{{B}^{0}_{s}}} for A P ​ ( Λ b 0 ) subscript 𝐴 P subscript superscript Λ 0 𝑏 A_{\rm P}({{\mathchar 28931\relax}^{0}_{b}}) , and | q / p | B s 0 subscript 𝑞 𝑝 subscript superscript 𝐵 0 𝑠 |q/p|_{{{B}^{0}_{s}}} for A P ​ ( B s 0 ) subscript 𝐴 P subscript superscript 𝐵 0 𝑠 A_{\rm P}({{B}^{0}_{s}}) .

$p_T$ [GeV/$c$]	$A_P{\left(B_s^0\right)}_{\sqrt{s}=7\mathrm{TeV}}$	$A_P{\left({\Lambda }_b^0\right)}_{\sqrt{s}=7\mathrm{TeV}}$
$(2.0,7.0)$	$-0.0166\pm 0.0393\pm 0.0082$	$-0.0130\pm 0.0311\pm 0.0133$
$(7.0,9.5)$	$0.0247\pm 0.0334\pm 0.0050$	$0.0948\pm 0.0476\pm 0.0211$
$(9.5,12.0)$	$0.0566\pm 0.0349\pm 0.0096$	$-0.0596\pm 0.0722\pm 0.0262$
$(12.0,30.0)$	$-0.0382\pm 0.0273\pm 0.0054$	$-0.0146\pm 0.0985\pm 0.0369$

Table 15. Table 15: Values of the production asymmetries in bins of y 𝑦 y , integrated over p T subscript 𝑝 T p_{\rm T} , for the B s 0 subscript superscript 𝐵 0 𝑠 {B}^{0}_{s} meson and the Λ b 0 subscript superscript Λ 0 𝑏 {\mathchar 28931\relax}^{0}_{b} baryon for data collected in proton-proton collisions at the centre-of-mass energy of 7 TeV TeV \mathrm{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second systematic. The uncertainties among the bins are correlated due to the external inputs: A C ​ P ​ ( B + → J / ψ ​ K + ) subscript 𝐴 𝐶 𝑃 → superscript 𝐵 𝐽 𝜓 superscript 𝐾 A_{{C\!P}}({{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}) , A D ( K ¯ ) 0 A_{\rm D}({{\kern 1.99997pt\overline{\kern-1.99997ptK}{}}{}^{0}}) , | q / p | B 0 subscript 𝑞 𝑝 superscript 𝐵 0 |q/p|_{{{B}^{0}}} and | q / p | B s 0 subscript 𝑞 𝑝 subscript superscript 𝐵 0 𝑠 |q/p|_{{{B}^{0}_{s}}} for A P ​ ( Λ b 0 ) subscript 𝐴 P subscript superscript Λ 0 𝑏 A_{\rm P}({{\mathchar 28931\relax}^{0}_{b}}) , and | q / p | B s 0 subscript 𝑞 𝑝 subscript superscript 𝐵 0 𝑠 |q/p|_{{{B}^{0}_{s}}} for A P ​ ( B s 0 ) subscript 𝐴 P subscript superscript 𝐵 0 𝑠 A_{\rm P}({{B}^{0}_{s}}) .

$y$	$A_P{\left(B_s^0\right)}_{\sqrt{s}=7\mathrm{TeV}}$	$A_P{\left({\Lambda }_b^0\right)}_{\sqrt{s}=7\mathrm{TeV}}$
$(2.1,3.0)$	$0.0151\pm 0.0445\pm 0.0088$	$-0.0511\pm 0.0399\pm 0.0168$
$(3.0,3.3)$	$0.0296\pm 0.0566\pm 0.0111$	$0.0514\pm 0.0448\pm 0.0171$
$(3.3,4.5)$	$-0.0554\pm 0.0432\pm 0.0101$	$0.0638\pm 0.0348\pm 0.0160$

Table 16. Table 16: Values of the production asymmetries in bins of p T subscript 𝑝 T p_{\rm T} , integrated over y 𝑦 y , for the B s 0 subscript superscript 𝐵 0 𝑠 {B}^{0}_{s} meson and the Λ b 0 subscript superscript Λ 0 𝑏 {\mathchar 28931\relax}^{0}_{b} baryon for data collected in proton-proton collisions at the centre-of-mass energy of 8 TeV TeV \mathrm{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second systematic. The uncertainties among the bins are correlated, due to the external inputs: A C ​ P ​ ( B + → J / ψ ​ K + ) subscript 𝐴 𝐶 𝑃 → superscript 𝐵 𝐽 𝜓 superscript 𝐾 A_{{C\!P}}({{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}) , A D ( K ¯ ) 0 A_{\rm D}({{\kern 1.99997pt\overline{\kern-1.99997ptK}{}}{}^{0}}) , | q / p | B 0 subscript 𝑞 𝑝 superscript 𝐵 0 |q/p|_{{{B}^{0}}} and | q / p | B s 0 subscript 𝑞 𝑝 subscript superscript 𝐵 0 𝑠 |q/p|_{{{B}^{0}_{s}}} , for A P ​ ( Λ b 0 ) subscript 𝐴 P subscript superscript Λ 0 𝑏 A_{\rm P}({{\mathchar 28931\relax}^{0}_{b}}) and | q / p | B s 0 subscript 𝑞 𝑝 subscript superscript 𝐵 0 𝑠 |q/p|_{{{B}^{0}_{s}}} for A P ​ ( B s 0 ) subscript 𝐴 P subscript superscript 𝐵 0 𝑠 A_{\rm P}({{B}^{0}_{s}}) .

$p_T$ [GeV/$c$]	$A_P{\left(B_s^0\right)}_{\sqrt{s}=8\mathrm{TeV}}$	$A_P{\left({\Lambda }_b^0\right)}_{\sqrt{s}=8\mathrm{TeV}}$
$(2.0,7.0)$	$0.0235\pm 0.0264\pm 0.0083$	$0.0292\pm 0.0200\pm 0.0096$
$(7.0,9.5)$	$0.0257\pm 0.0223\pm 0.0049$	$0.0367\pm 0.0302\pm 0.0127$
$(9.5,12.0)$	$-0.0286\pm 0.0230\pm 0.0053$	$0.0442\pm 0.0437\pm 0.0164$
$(12.0,30.0)$	$0.0187\pm 0.0186\pm 0.0049$	$0.0902\pm 0.0612\pm 0.0253$

Table 17. Table 17: Values of the production asymmetries in bins of y 𝑦 y , integrated over p T subscript 𝑝 T p_{\rm T} , for the B s 0 subscript superscript 𝐵 0 𝑠 {B}^{0}_{s} meson and the Λ b 0 subscript superscript Λ 0 𝑏 {\mathchar 28931\relax}^{0}_{b} baryon for data collected in proton-proton collisions at the centre-of-mass energy of 8 TeV TeV \mathrm{\,Te\kern-1.00006ptV} . The first uncertainties are statistical and the second systematic. The uncertainties among the bins are correlated, due to the external inputs: A C ​ P ​ ( B + → J / ψ ​ K + ) subscript 𝐴 𝐶 𝑃 → superscript 𝐵 𝐽 𝜓 superscript 𝐾 A_{{C\!P}}({{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}) , A D ( K ¯ ) 0 A_{\rm D}({{\kern 1.99997pt\overline{\kern-1.99997ptK}{}}{}^{0}}) , | q / p | B 0 subscript 𝑞 𝑝 superscript 𝐵 0 |q/p|_{{{B}^{0}}} and | q / p | B s 0 subscript 𝑞 𝑝 subscript superscript 𝐵 0 𝑠 |q/p|_{{{B}^{0}_{s}}} , for A P ​ ( Λ b 0 ) subscript 𝐴 P subscript superscript Λ 0 𝑏 A_{\rm P}({{\mathchar 28931\relax}^{0}_{b}}) and | q / p | B s 0 subscript 𝑞 𝑝 subscript superscript 𝐵 0 𝑠 |q/p|_{{{B}^{0}_{s}}} for A P ​ ( B s 0 ) subscript 𝐴 P subscript superscript 𝐵 0 𝑠 A_{\rm P}({{B}^{0}_{s}}) .

$y$	$A_P{\left(B_s^0\right)}_{\sqrt{s}=8\mathrm{TeV}}$	$A_P{\left({\Lambda }_b^0\right)}_{\sqrt{s}=8\mathrm{TeV}}$
$(2.1,3.0)$	$0.0364\pm 0.0290\pm 0.0103$	$0.0028\pm 0.0247\pm 0.0107$
$(3.0,3.3)$	$-0.0078\pm 0.0413\pm 0.0063$	$0.0792\pm 0.0317\pm 0.0138$
$(3.3,4.5)$	$0.0055\pm 0.0298\pm 0.0070$	$0.0682\pm 0.0242\pm 0.0142$

Equations

$$A_{\mathrm{P}}\left(x\right)$$

$$S\left(t,\psi ,\xi \right)$$

$${\Omega }_{\pm }^{\xi }\equiv {\delta }_{+1\xi }\left(1-A_{\mathrm{P}}\right){\left|\frac{q}{p}\right|}^{1-\psi }\pm {\delta }_{-1\xi }\left(1+A_{\mathrm{P}}\right){\left|\frac{q}{p}\right|}^{-1-\psi },$$

$$A_{\mathrm{D}}\equiv \frac{{\epsilon }_{\bar{f}}-{\epsilon }_f}{{\epsilon }_{\bar{f}}+{\epsilon }_f}.$$

$$A_{CP}\equiv \frac{\mathcal{B}\left(\overline{B}{}_{(s)}^0\to \bar{f}\right)-\mathcal{B}\left(B_{(s)}^0\to f\right)}{\mathcal{B}\left(\overline{B}{}_{(s)}^0\to \bar{f}\right)+\mathcal{B}\left(B_{(s)}^0\to f\right)}$$

$$A_{\mathrm{raw}}\equiv \frac{N(B^{-}\to J/\psi K^{-})-N(B^{+}\to J/\psi K^{+})}{N(B^{-}\to J/\psi K^{-})+N(B^{+}\to J/\psi K^{+})},$$

$$A_{\mathrm{raw}}(B^{+}\to J/\psi K^{+})=A_{\mathrm{P}}(B^{+})+A_{\mathrm{D}}(K^{+})+A_{CP}(B^{+}\to J/\psi K^{+}),$$

$$A_{\mathrm{P}}({\text{mathchar}28931}_b^0)=-\left[\frac{f_u}{f_{{\text{mathchar}28931}_b^0}}{A}_{\mathrm{P}}(B^{+})+\frac{f_d}{f_{{\text{mathchar}28931}_b^0}}{A}_{\mathrm{P}}(B^0)+\frac{f_s}{f_{{\text{mathchar}28931}_b^0}}{A}_{\mathrm{P}}(B_s^0)+\frac{f_c}{f_{{\text{mathchar}28931}_b^0}}{A}_{\mathrm{P}}(B_c^{+})+\frac{f_{\mathrm{other}}}{f_{{\text{mathchar}28931}_b^0}}{A}_{\mathrm{P}}(\mathrm{other})\right]$$

$$A_{\mathrm{P}}({\Lambda }_b^0)\simeq -\left[\frac{f_u}{f_{{\text{mathchar}28931}_b^0}}{A}_{\mathrm{P}}(B^{+})+\frac{f_d}{f_{{\text{mathchar}28931}_b^0}}{A}_{\mathrm{P}}(B^0)+\frac{f_s}{f_{{\text{mathchar}28931}_b^0}}{A}_{\mathrm{P}}(B_s^0)\right].$$

$$A_{\mathrm{P}}=\sum _i\frac{N_i}{{\epsilon }_i}{A}_{\mathrm{P},i}/\sum _i\frac{N_i}{{\epsilon }_i}$$

$$N_i=\mathcal{L}{\sigma }_{b\bar{b}}2f_q\mathcal{B}{F}_i{\epsilon }_i$$

$$A_{\mathrm{P}}=\sum _i{\omega }_i{A}_{\mathrm{P},i}.$$

$$g(m)\propto {\int }_0^{+\infty }(m^{\prime }{)}^{s}G\left(m+m^{\prime };\mu \right)\mathrm{d}{m}^{\prime }$$

$$S\left(t,\psi \right)$$

$$ϵ\left(t\right)=\frac{1}{2}\left[1-\mathrm{erf}\left(p_1{t}^{p_2}\right)\right]\left(1+p_{3}t\right),$$

$$ϵ\left(t\right)=\frac{1}{2}\left[1-\frac{1}{2}\mathrm{erf}\left(\frac{n_1-t}{t}\right)-\frac{1}{2}\mathrm{erf}\left(\frac{n_2-t}{t}\right)\right]\left(1+n_{3}t\right),$$

$${\Lambda }_{\pm }\equiv \left(1-A_{\mathrm{P}}\right){\left|\frac{q}{p}\right|}^{1-\psi }\pm \left(1+A_{\mathrm{P}}\right){\left|\frac{q}{p}\right|}^{-1-\psi },$$

$$\delta z_i=\frac{N_{i,b}-N_{i,\overline{b}}}{N_{i,b}+N_{i,\overline{b}}}\cdot \frac{1}{f_{{\text{mathchar}28931}_b^0}}.$$

$$\frac{f_{\mathrm{other}}}{f_{{\text{mathchar}28931}_b^0}}{A}_{\mathrm{P}}(\mathrm{other})\simeq \frac{f_{{\Xi }_b}}{f_{{\text{mathchar}28931}_b^0}}{A}_{\mathrm{P}}({\Xi }_b).$$

$${\omega }_i^{\mathrm{data}}=\frac{N_i}{{\epsilon }_i^{\mathrm{total}}}/\sum _i\frac{N_i}{{\epsilon }_i^{\mathrm{total}}}$$

$$A_{\mathrm{P}}(B^{+}{)}_{\sqrt{s}=7\mathrm{Te}\mathrm{V}}$$

$$A_{\mathrm{P}}(B^{+}{)}_{\sqrt{s}=8\mathrm{Te}\mathrm{V}}$$

$$A_{\mathrm{P}}(B^0{)}_{\sqrt{s}=7\mathrm{Te}\mathrm{V}}$$

$$A_{\mathrm{P}}(B^0{)}_{\sqrt{s}=8\mathrm{Te}\mathrm{V}}$$

$$A_{\mathrm{P}}(B_s^0{)}_{\sqrt{s}=7\mathrm{Te}\mathrm{V}}$$

$$A_{\mathrm{P}}(B_s^0{)}_{\sqrt{s}=8\mathrm{Te}\mathrm{V}}$$

$$A_{\mathrm{P}}({\text{mathchar}28931}_b^0{)}_{\sqrt{s}=7\mathrm{Te}\mathrm{V}}$$

$$A_{\mathrm{P}}({\text{mathchar}28931}_b^0{)}_{\sqrt{s}=8\mathrm{Te}\mathrm{V}}$$

Full text

EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN)

​​​

CERN-EP-2017-036

LHCb-PAPER-2016-062

24 March 2017

Measurement of ${{B}^{0}}$, ${{B}^{0}_{s}}$, ${{{B}^{+}}}$ and ${{\mathchar 28931\relax}^{0}_{b}}$ production asymmetries in 7 and 8 TeV proton-proton collisions

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

The ${{B}^{0}}$, ${{B}^{0}_{s}}$, ${{{B}^{+}}}$ and ${{\mathchar 28931\relax}^{0}_{b}}$ hadron production asymmetries are measured using a data sample corresponding to an integrated luminosity of 3.0$\mbox{\,fb}^{-1}$, collected by the LHCb experiment in proton-proton collisions at centre-of-mass energies of 7 and 8$\mathrm{\,Te\kern-1.00006ptV}$. The measurements are performed as a function of transverse momentum and rapidity of the ${b}$ hadrons within the LHCb detector acceptance. The overall production asymmetries, integrated over transverse momentum and rapidity, are also determined.

Published in Phys. Lett. B 774 (2017) 139-158

© CERN on behalf of the LHCb collaboration, licence CC-BY-4.0.

1 Introduction

The production rates of ${b}$ and ${\overline{b}}$ hadrons are not expected to be identical in proton-proton collisions, as ${b}$ and ${\overline{b}}$ quarks, produced in a hard scattering at the partonic level, might have different probabilities for coalescing with ${u}$ or ${d}$ valence quarks from the beam remnant. As a consequence, the production rates of ${{B}^{+}}$ and ${B}^{0}$ mesons may exceed those of ${{B}^{-}}$ and ${\kern 1.79993pt\overline{\kern-1.79993ptB}{}}{}^{0}$, and $b$ baryons can be produced more abundantly than $\bar{b}$ baryons. In the case of ${B}^{0}_{s}$ and ${\kern 1.79993pt\overline{\kern-1.79993ptB}{}}{}^{0}_{s}$ the production rates depend on the values of the other production asymmetries as no valence strange quark is present within the colliding protons and $b$ and $\overline{b}$ quarks are predominantly produced in pairs.

The LHCb detector, thanks to its unique geometry as a forward spectrometer, is particularly suited to measure such asymmetries, as they are expected to be enhanced at forward rapidities and small transverse momenta. Other subtle effects of quantum chromodynamics, beyond the coalescence of ${b}$ quarks and light valence quarks, may also contribute [1, 2, 3].

The measurements of hadron production asymmetries are of primary importance, not only for the understanding of the production mechanisms, but also for enabling precise measurements of ${C\!P}$ violation in $c$ and $b$ hadrons at the LHC. Indeed, observed asymmetries must be corrected for production effects to obtain the $C\!P$ asymmetries in the decays. Simulations that model the non-perturbative fragmentation of ${b}$ quarks in proton-proton collisions at LHC energies predict asymmetries generally up to a few percent [4, 5]. Production asymmetries of ${B}^{0}$ and ${B}^{0}_{s}$ mesons have been measured by LHCb at a centre-of-mass energy of 7$\mathrm{\,Te\kern-1.00006ptV}$, excluding values larger than a few percent [6]. The LHCb collaboration has also searched for possible production asymmetries for ${{D}^{+}}$ and ${{D}^{+}_{s}}$ mesons, finding the integrated ${D}^{+}$ production asymmetry different from zero at approximately three standard deviations [7, 8]. In the $b$-baryon sector, the LHCb collaboration measured the sum of the ${\mathchar 28931\relax}^{0}_{b}$–${\kern 1.00006pt\overline{\kern-1.00006pt\mathchar 28931\relax}}{}^{0}_{b}$ production asymmetry and the $C\!P$ asymmetry in the ${{\mathchar 28931\relax}^{0}_{b}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}p{{K}^{-}}$ decay [9], finding evidence for a dependence on ${\mathchar 28931\relax}^{0}_{b}$ rapidity.

This paper reports measurements of the production asymmetries, $A_{\mathrm{P}}\left({{{B}^{+}}}\right)$, $A_{\mathrm{P}}\left({{B}^{0}}\right)$ and $A_{\mathrm{P}}\left({{B}^{0}_{s}}\right)$, measured using ${{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}$, ${{B}^{0}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{*0}}$ and ${{B}^{0}_{s}}\rightarrow{{D}^{-}_{s}}{{\pi}^{+}}$ decays. In addition, a measurement of $A_{\mathrm{P}}\left({{\mathchar 28931\relax}^{0}_{b}}\right)$, determined indirectly from the other asymmetries, is presented. Hereafter, ${K}^{*0}$ is used to refer to the $K^{*}(892)^{0}$ and the inclusion of charge-conjugate decay modes is implied throughout, except when referring to the production asymmetries, which are defined as

[TABLE]

where $\sigma$ denotes the inclusive production cross-section in a given region of phase space. The data sample, collected by LHCb in proton-proton collisions, corresponds to an integrated luminosity of 1.0$\mbox{\,fb}^{-1}$ at a centre-of-mass energy of 7$\mathrm{\,Te\kern-1.00006ptV}$, and 2.0$\mbox{\,fb}^{-1}$ at 8$\mathrm{\,Te\kern-1.00006ptV}$. The measurements are performed as a function of both the component of the momentum transverse to the beam ($p_{\rm T}$) and the rapidity ($y$) of the hadrons within the LHCb detector acceptance, and are then integrated over the ranges 0<\mbox{p_{\rm T}}<30$${\mathrm{\,Ge\kern-1.00006ptV\!/}c} and/or $2.1<y<4.5$ for ${{B}^{+}}$ and ${B}^{0}$ decays, and 2<\mbox{p_{\rm T}}<30$${\mathrm{\,Ge\kern-1.00006ptV\!/}c} and/or $2.1<y<4.5$ for ${B}^{0}_{s}$ and ${\mathchar 28931\relax}^{0}_{b}$ decays. The ranges in $p_{\rm T}$ are not identical due to different trigger requirements between decays with and without muons in the final states. This analysis improves the previous one performed on ${B}^{0}$ and ${B}^{0}_{s}$ production asymmetries [6], using a larger data sample and a finer binning scheme for investigating the dependence on $p_{\rm T}$ and $y$. In addition, new measurements of ${{B}^{+}}$ and ${\mathchar 28931\relax}^{0}_{b}$ production asymmetries have been included. Unlike in the previous analysis, the ${{B}^{0}}\rightarrow{{D}^{-}}{{\pi}^{+}}$ decay is not considered, as it has been found not to improve the precision on the ${B}^{0}$ measurement.

2 Detector, trigger and simulation

The LHCb detector [10] 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 proton-proton interaction region, a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about $4{\rm\,Tm}$, and three stations of silicon-strip detectors and straw drift tubes placed downstream of the magnet.

The tracking system provides a measurement of momentum, $p$, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200${\mathrm{\,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/\mbox{p_{\rm T}}){\,\upmu\rm m}, where $p_{\rm T}$ is measured in ${\mathrm{\,Ge\kern-1.00006ptV\!/}c}$. Different types of charged hadrons are distinguished using information from two ring-imaging Cherenkov detectors. Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillating-pad 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. The trigger [11] 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.

For ${{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}$ and ${{B}^{0}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{*0}}$ decays, the data are collected by using the hardware muon trigger, which requires a single muon with large transverse momentum (from p_{\mathrm{T}}>1.4$${\mathrm{\,Ge\kern-1.00006ptV\!/}c} to p_{\mathrm{T}}>1.8$${\mathrm{\,Ge\kern-1.00006ptV\!/}c}) or a pair of muons with a large product of their transverse momenta (from $\sqrt{p_{\mathrm{T,1}}\,p_{\mathrm{T,2}}}>1.3{\mathrm{\,Ge\kern-1.00006ptV\!/}c}$ to $\sqrt{p_{\mathrm{T,1}}\,p_{\mathrm{T,2}}}>1.6{\mathrm{\,Ge\kern-1.00006ptV\!/}c}$), depending on the data-taking conditions. For ${{B}^{0}_{s}}\rightarrow{{D}^{-}_{s}}{{\pi}^{+}}$ decays, data are collected using the hadronic hardware trigger, which requires at least one cluster in the hadronic calorimeter with a transverse energy greater than 3.5$$\mathrm{\,Ge\kern-1.00006ptV} or 3.7$$\mathrm{\,Ge\kern-1.00006ptV}, depending on the data-taking period. The output is then processed by the software trigger. In the case of ${{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}$ and ${{B}^{0}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{*0}}$ decays, ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}$ mesons consistent with coming from the decay of a ${b}$-hadron are selected by requiring that their decay products form a displaced vertex and have large IPs at the PV with respect to which the $B$ candidate has the smallest $\chi^{2}_{\rm IP}$. The quantity $\chi^{2}_{\rm IP}$ is defined as the difference in the vertex-fit $\chi^{2}$ of a given PV reconstructed with and without the particle under consideration. The ${{B}^{0}_{s}}\rightarrow{{D}^{-}_{s}}{{\pi}^{+}}$ decays are selected by requiring a two- or three-track secondary vertex with a significant displacement from all PVs. At least one charged particle must have a transverse momentum \mbox{p_{\rm T}}>1.7$${\mathrm{\,Ge\kern-1.00006ptV\!/}c} and be inconsistent with originating from a PV. A multivariate algorithm [12] is used for the identification of secondary vertices consistent with the decay of a $b$ hadron.

Simulated events are used to determine the signal selection efficiency as a function of $p_{\mathrm{T}}$ and $y$, and to study the modelling of the decay-time resolution, the reconstruction efficiency as function of the decay time and the shape of the invariant mass distribution of partially reconstructed background. In the simulation, proton-proton collisions are generated using Pythia [13, 4] with a specific LHCb configuration [14]. Decays of hadronic particles are described by EvtGen [15], in which final-state radiation is generated using Photos [16]. The interaction of the generated particles with the detector, and its response, are implemented using the Geant4 toolkit [17, *Agostinelli:2002hh] as described in Ref. [19].

3 Methodology

The asymmetries $A_{\mathrm{P}}\left({{B}^{0}}\right)$ and $A_{\mathrm{P}}\left({{B}^{0}_{s}}\right)$ are measured by means of a time-dependent analysis of ${{B}^{0}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{*0}}$ decays, with ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}\rightarrow{\mu^{-}}{\mu^{+}}$ and ${{K}^{*0}}\rightarrow{{K}^{-}}{{\pi}^{+}}$, and ${{B}^{0}_{s}}\rightarrow{{D}^{-}_{s}}{{\pi}^{+}}$ decays, with ${{D}^{-}_{s}}\rightarrow{{K}^{+}}{{K}^{-}}{{\pi}^{-}}$. The decay rate to a flavour-specific final state $f$ of a $B^{0}_{(s)}$ meson with average decay width $\Gamma_{{d}({s})}$ can be written as

[TABLE]

where $\Delta m_{{d}({s})}\equiv m_{{d}({s}),\,\rm H}-m_{{d}({s}),\,\rm L}$ and $\Delta\Gamma_{{d}({s})}\equiv\Gamma_{{d}({s}),\,\rm L}-\Gamma_{{d}({s}),\,\rm H}$ are the mass and width differences of the $B^{0}_{(s)}$–${\kern 1.79993pt\overline{\kern-1.79993ptB}{}}{}^{0}_{(s)}$ system mass eigenstates. The subscripts H and L denote the heavy and light eigenstates, respectively. The symbol $\psi$ is the tag of the final state, which assumes the values $\psi=1$ if the final state is $f$ and $\psi=-1$ if the final state is the $C\!P$ conjugate $\bar{f}$, while $\xi$ indicates the tag of the initial flavour of the $B^{0}_{(s)}$ meson, which takes the values $\xi=1$ for $B^{0}_{(s)}$ and $\xi=-1$ for ${\kern 1.79993pt\overline{\kern-1.79993ptB}{}}{}^{0}_{(s)}$. The terms $\Omega_{+}^{\xi}$ and $\Omega_{-}^{\xi}$ are defined as

[TABLE]

where $p$ and $q$ are complex parameters entering the definition of the two mass eigenstates of the effective Hamiltonian of the $B^{0}_{(s)}$ system, $|B_{\mathrm{H}}\rangle=p|B^{0}_{(s)}\rangle-q|{\kern 1.79993pt\overline{\kern-1.79993ptB}{}}{}^{0}_{(s)}\rangle$ and $|B_{\mathrm{L}}\rangle=p|B^{0}_{(s)}\rangle+q|{\kern 1.79993pt\overline{\kern-1.79993ptB}{}}{}^{0}_{(s)}\rangle$, and $\delta_{ij}$ is the Kronecker delta. The symbol $A_{\rm D}$ represents the detection asymmetry of the final state, defined in terms of the $f$ and $\bar{f}$ detection efficiencies, $\varepsilon$, as

[TABLE]

The direct $C\!P$ asymmetry $A_{{C\!P}}$ is defined as

[TABLE]

where the symbol $\mathcal{B}$ stands for the branching fraction of the decay considered.

The asymmetry $A_{\mathrm{P}}\left({{{B}^{+}}}\right)$ is measured by means of a time-integrated analysis of ${{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}$ decays, with ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}\rightarrow{\mu^{+}}{\mu^{-}}$, starting from the raw asymmetry defined as

[TABLE]

where $N$ denotes the observed yields. The raw asymmetry can be written, up to ${\cal O}(10^{-6})$ corrections, as

[TABLE]

where $A_{\rm D}({{K}^{+}})$ is the ${K}^{+}$ detection asymmetry, measured by means of charm control samples as in Ref. [20], and $A_{C\!P}({{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}})$ is the $C\!P$ asymmetry in the decay, measured by BaBar, Belle and D0 [21, 22, 23]. An improved measurement of the $C\!P$ asymmetry was also made recently by LHCb [24], using an independent data sample selected with different trigger requirements. The $A_{\mathrm{P}}$ values obtained from Eq. 3 and Eq. 2 are detector-independent quantities only if measured in kinematic regions where the reconstruction efficiencies are constant. To account for the dependence of the production asymmetries on the kinematics of the ${{B}^{+}}$, ${B}^{0}$ and ${B}^{0}_{s}$ mesons, each data sample is divided into bins of ($p_{\rm T}$, $y$), and the measurement is performed for each bin. Figure 1 shows the distribution of ($p_{\rm T}$, $y$) for ${{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}$, ${{B}^{0}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{*0}}$ and ${{B}^{0}_{s}}\rightarrow{{D}^{-}_{s}}{{\pi}^{+}}$ decays, where the background components are subtracted using the sPlot technique [25] and the definition of the various kinematic bins is overlaid. For the ${{B}^{+}}$ and ${B}^{0}$ decays a common set of bins is used, defined in Table LABEL:tab:resultsBpB02011 of the Appendix, and in the case of the ${B}^{0}_{s}$ decay, the binning scheme is reported in Table 8.

In proton-proton collisions at the LHC, ${b}$ and ${\overline{b}}$ quarks are predominantly pair-produced via strong interaction processes. This leads to a relation between the ${{\mathchar 28931\relax}^{0}_{b}}$ production asymmetry and the other ${b}$-hadron production asymmetries, namely

[TABLE]

where $f_{u}$, $f_{d}$, $f_{s}$, $f_{c}$, $f_{{{\mathchar 28931\relax}^{0}_{b}}}$ and $f_{\rm other}$ are the fragmentation fractions of a $b$ quark hadronizing into weakly-decaying ${{B}^{+}}$, ${B}^{0}$, ${B}^{0}_{s}$, ${B}_{c}^{+}$ mesons, ${\mathchar 28931\relax}^{0}_{b}$ baryons and all the other ${b}$-baryon species. The ratios of the fragmentation fractions, $f_{u}/f_{{{\mathchar 28931\relax}^{0}_{b}}}$, $f_{d}/f_{{{\mathchar 28931\relax}^{0}_{b}}}$ and $f_{s}/f_{{{\mathchar 28931\relax}^{0}_{b}}}$ are taken from LHCb measurements reported in Refs. [26, 27]. Their dependence on $p_{\rm T}$ and $y$ is taken into account. The terms $(f_{c}/f_{{{\mathchar 28931\relax}^{0}_{b}}})\cdot A_{P}({{B}_{c}^{+}})$ and $(f_{\rm other}/f_{{{\mathchar 28931\relax}^{0}_{b}}})\cdot A_{P}(\rm other)$ are of the order of $3\cdot 10^{-5}$ and $2\cdot 10^{-3}$, respectively. This is estimated assuming that the value of $A_{P}({{B}_{c}^{+}})$ and $A_{P}(\rm other)$ are of the same order as the $B$-meson production asymmetries ($\simeq 10^{-2}$) and taking the values of $f_{c}/f_{{{\mathchar 28931\relax}^{0}_{b}}}$ and $f_{\rm other}/f_{{{\mathchar 28931\relax}^{0}_{b}}}$ from simulation. Neglecting these terms, the ${\mathchar 28931\relax}^{0}_{b}$ production asymmetry can be measured using the approximate relation

[TABLE]

Possible small deviations from this approximation, due in particular to contributions from other $b$ baryons, are taken into account in the evaluation of systematic uncertainties.

3.1 Integrated production asymmetries

In addition to the measurements in bins, integrated production asymmetries, where efficiency corrections have been applied, are also provided. The integration of the $A_{\rm P}$ values is performed in the ranges 0<\mbox{p_{\rm T}}<30$${\mathrm{\,Ge\kern-1.00006ptV\!/}c} and $2.1<y<4.5$ for the ${{{B}^{+}}}$ and ${B}^{0}$ decays and in the ranges 2<p_{\mathrm{T}}<30$${\mathrm{\,Ge\kern-1.00006ptV\!/}c} and $2.1<y<4.5$ for the ${{B}^{0}_{s}}$ and ${\mathchar 28931\relax}^{0}_{b}$ decays. The integrated value of $A_{\rm P}$ is given by

[TABLE]

where the index $i$ runs over the bins, $N_{i}$ is the number of observed signal events in the $i$-th bin and $\varepsilon_{i}$ is the efficiency defined as the number of selected events divided by the number of produced events in the $i$-th bin. The signal yield in each bin can be expressed as

[TABLE]

where $\mathcal{L}$ is the integrated luminosity, $\sigma_{b\bar{b}}$ is the $b\bar{b}$ cross section, $f_{q}$ is the fragmentation fraction for quark flavour $q$, with $q\in\{{u},{d},{s}\}$, $F_{i}$ stands for the fraction of the $b$ hadrons produced in the $i$-th bin and $\mathcal{B}$ is the branching fraction of the ${b}$-hadron decay being considered. By substituting $N_{i}/\varepsilon_{i}$ from Eq. 5 into Eq. 4, the integrated value of $A_{\rm{P}}$ becomes

[TABLE]

where $\omega_{i}=F_{i}/\sum_{i}F_{i}$. The $\omega_{i}$ values are determined using simulated events, generated with proton-proton collisions at the centre-of-mass energies of 7 and 8$\mathrm{\,Te\kern-1.00006ptV}$.

4 Data set and event selections

The selections of ${{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}$ and ${{B}^{0}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{*0}}$ decays are based on the reconstruction of ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}\rightarrow{\mu^{-}}{\mu^{+}}$ decays combined with either a track identified as a kaon or with a ${K}^{*0}$ decaying to ${{K}^{+}}{{\pi}^{-}}$. The ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}$ candidates are formed from two oppositely charged tracks originating from a common vertex, identified as muons with \mbox{p_{\rm T}}>500$${\mathrm{\,Me\kern-1.00006ptV\!/}c}. The ${{K}^{*0}}$ candidates are formed from two oppositely charged tracks, one identified as a kaon and the other as a pion, originating from the same vertex. They are required to have \mbox{p_{\rm T}}>1$${\mathrm{\,Ge\kern-1.00006ptV\!/}c} and the ${K}^{+}$ ${\pi}^{-}$ invariant mass in the range 826\text{--}966$${\mathrm{\,Me\kern-1.00006ptV\!/}c^{2}}. The invariant mass of ${B}^{0}$ and ${{B}^{+}}$ candidates, calculated constraining the two muon candidates to have the known ${J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}$ mass, is required to be in the range 5150\text{--}5450$${\mathrm{\,Me\kern-1.00006ptV\!/}c^{2}}. The proper decay time of the $B$-meson candidate is calculated from a fit that constrains the candidate to originate from the PV with the smallest $\chi^{2}_{\rm IP}$ with respect to the $B$ candidate. Only $B$-meson candidates with a decay time greater than 0.2 ps are retained. This lower bound on the decay time rejects a large fraction of the combinatorial background.

In the case of ${{B}^{0}_{s}}\rightarrow{{D}^{-}_{s}}{{\pi}^{+}}$ decays, the ${D}^{-}_{s}$ candidates are reconstructed using the ${{K}^{+}}{{K}^{-}}{{\pi}^{-}}$ decay channel. Requirements are applied to the ${D}^{-}_{s}$ decay products before combining them to form a common vertex, namely the scalar $p_{\rm T}$ sum of the tracks must exceed 1.8${\mathrm{\,Ge\kern-1.00006ptV\!/}c}$ and the largest distance of closest approach between all possible pairs of tracks must be less than 0.5 mm. The ${D}^{-}_{s}$ candidates are then required to be significantly detached from the PV and to have the invariant mass within the range 1949\text{--}1989$${\mathrm{\,Me\kern-1.00006ptV\!/}c^{2}}. Each ${D}^{-}_{s}$ candidate is subsequently combined with a second charged pion, referred to as the accompanying pion in the following, to form the ${B}$-meson decay vertex. The sum of the $p_{\rm T}$ values of the ${D}^{-}_{s}$ and accompanying ${\pi}^{+}$ must be larger than 5${\mathrm{\,Ge\kern-1.00006ptV\!/}c}$ and the decay time of ${B}$-meson candidates must be greater than 0.2 ps. Furthermore, the cosine of the angle between the $B$-meson candidate momentum vector and the vector connecting the PV and $B$-meson candidate vertex is required to be larger than 0.999.

Stringent particle identification criteria are required to be satisfied for the kaons and pions forming the ${{K}^{*0}}$ and ${D}^{-}_{s}$ candidates, the kaon from the ${{B}^{+}}$ decay and the accompanying pion, in order to reduce to a negligible level the background from other ${B}$-meson decays with a misidentified kaon or pion, and from ${\mathchar 28931\relax}^{0}_{b}$ decays with a misidentified proton.

A final selection is applied using a multivariate analysis method based on a Boosted Decision Tree [28, *Roe], where the variables used in the selection are: the $p_{\rm T}$ and the IP of the $B$ decay products, the flight distance and the IP of the $B$ candidate, and, in the case of ${B}^{0}_{s}$, the flight distance of the ${D}^{-}_{s}$ meson. The multivariate selection is trained using simulated events as a proxy for the signal, and $B$-meson candidates from data selected in the upper mass sidebands to represent the background.

5 Fit model

For each signal and background component, the invariant mass distribution of all $B$ candidates, and, in the case of $B^{0}_{(s)}$, the decay time, is modelled by defining appropriate probability density functions (PDFs). Two categories of background are considered: the combinatorial background, due to the random association of tracks, and the partially reconstructed background, due to decays with a topology similar to that of the signal, but with one or more particles not reconstructed. The latter is only relevant for ${{B}^{0}_{s}}\rightarrow{{D}^{-}_{s}}{{\pi}^{+}}$ decays.

5.1 Invariant mass parameterization

The signal component for $B$ mesons is modelled by convolving a sum of two Gaussian functions with a function parameterizing the final-state QED radiation (FSR). The PDF of the invariant mass, $m$, is given by the convolution

[TABLE]

where $G$ is the sum of two Gaussian functions with different widths and common mean $\mu$ that represents the ${B}$-meson mass. The parameter $s$ governs the amount of FSR, and using simulation is found to be $s=-0.9966\pm 0.0005$ for the ${{B}^{+}}$ decay, $s=-0.9945\pm 0.0003$ for the ${B}^{0}$ decay and $s=-0.9832\pm 0.0004$ for the ${B}^{0}_{s}$ decay. The invariant mass shape of the combinatorial background is well described by an exponential PDF.

Regarding the ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}$ invariant mass spectrum, common parameters are used for both ${{B}^{+}}$ and ${{B}^{-}}$ mesons. In the case of the ${{D}^{-}_{s}}{{\pi}^{+}}$ spectrum, a background component due to partially reconstructed ${B}^{0}_{s}$ decays is also present in the low invariant mass region. The contributions with the highest branching fractions are from the ${{B}^{0}_{s}}\rightarrow D_{s}^{*-}{{\pi}^{+}}$ decay, with $D_{s}^{*-}\rightarrow{{D}^{-}_{s}}\gamma$ or $D_{s}^{*-}\rightarrow{{D}^{-}_{s}}{{\pi}^{0}}$, where the $\gamma$ or $\pi^{0}$ is not reconstructed, and from the ${{B}^{0}_{s}}\rightarrow D_{s}^{-}\rho^{+}$ decay, with $\rho^{+}\rightarrow{{\pi}^{+}}{{\pi}^{0}}$, where the ${{\pi}^{0}}$ is not reconstructed. The partially reconstructed components are parameterized by means of a kernel estimation technique [30] based on invariant mass distributions obtained from simulated events, where the same selection applied to data is used and differences in invariant mass resolution between data and simulation are taken into account. The yields are obtained from the fits.

In the case of the ${{B}^{0}_{s}}\rightarrow{{D}^{-}_{s}}{{\pi}^{+}}$ decay, an irreducible background component due to the ${{B}^{0}}\rightarrow{{D}^{+}_{s}}{{\pi}^{-}}$ decay is also present. This component is accounted for in the fits using the same parameterization adopted for the signal, where the mean values of the two signal PDFs are separated by the difference in the known masses between ${B}^{0}$ and ${B}^{0}_{s}$ mesons [31] and the production asymmetry is fixed to the ${B}^{0}$ measured value. The yield of this component is fixed according to the known branching fraction [31].

5.2 Decay time parameterization

Starting from Eq. 3 and summing over $\xi$, the decay rate to a flavour-specific final state of a neutral ${B}$ meson is parameterized by the convolution

[TABLE]

where $R(t)$ is a function describing the decay-time resolution, as discussed in Sec. 5.3, and $\epsilon(t)$ is the reconstruction efficiency as a function of the decay time determined from simulation and parameterized for the ${B}^{0}$ decay by

[TABLE]

and for the ${{B}^{0}_{s}}$ decay by

[TABLE]

where erf is the error function, and $p_{i}$ and $n_{i}$ are parameters determined from simulation. The terms $\Lambda_{+}$ and $\Lambda_{-}$ are defined as

[TABLE]

and the term $A_{{C\!P}}A_{\rm D}$ is neglected, as $A_{\rm D}$ is ${\cal O}(10^{-2})$ [32] and $A_{{C\!P}}$ is very small for the decays under study. For this reason, it is only possible for the fit to determine the sum of $A_{\rm D}$ and $A_{{C\!P}}$, but not their individual values.

The decay-time PDF of the combinatorial background is studied using events from a high invariant mass window where the signal is not present, namely in the range 5310\text{--}5340$${\mathrm{\,Me\kern-1.00006ptV\!/}c^{2}} for ${{B}^{0}}\!\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{*0}}$ and 5450\text{--}5900$${\mathrm{\,Me\kern-1.00006ptV\!/}c^{2}} for ${{B}^{0}_{s}}\rightarrow{{D}^{-}_{s}}{{\pi}^{+}}$ decays. The partially reconstructed component for the ${{B}^{0}_{s}}\rightarrow{{D}^{-}_{s}}{{\pi}^{+}}$ decay is determined from simulated events.

5.3 Decay time resolution

The decay-time resolutions of ${{B}^{0}}$ and ${{B}^{0}_{s}}$ mesons are estimated by studying the decay time of fake $B$ candidates, formed from a ${D}^{-}$ decaying to ${K}^{+}$ ${\pi}^{-}$ ${\pi}^{-}$ and a pion track, both coming from the same PV. These $B$ candidates are called fake, as the probability to form a real decay with this technique is negligible. In order to avoid the introduction of biases in the decay-time measurements, the accompanying pion is selected with requirements on momentum and $p_{\rm T}$, rather than on IP. The decay-time distribution of these fake $B$ candidates yields an estimate of the decay-time resolution of a real decay. This method is verified by means of simulated events, both for signal and fake $B$ decays. The resolution model, $R(t)$, consisting of a sum of three Gaussian functions with zero mean and three different widths, characterized by an average width of 49 fs, is used. The resolution is found to be overestimated by about 4 fs and to be dependent on the decay time. Taking these effects into account, an uncertainty of 8 fs on the average width is considered as a systematic uncertainty. It is estimated from simulation that the measurement of the decay time is biased by no more than 2 fs, and this effect is also accounted for as a systematic uncertainty.

6 Determination of the production asymmetries

The production asymmetries are determined by means of unbinned ($B^{0}_{(s)}$) and binned (${{B}^{+}}$) maximum likelihood fits, for each kinematic bin, to the invariant mass (${{B}^{+}}$) and invariant mass and decay time ($B^{0}_{(s)}$) distributions, using the models described in the previous section. The models are validated with a series of fits to the mass and lifetime distributions of events obtained from pseudoexperiments. No evidence of biases on central values nor on the uncertainty is found. Furthermore, a global fit to the total sample of selected candidates is performed for each of the three decay modes to validate the fitting model on data. In the case of the time-dependent analysis, the mass differences $\Delta m_{{d}}$ and $\Delta m_{{s}}$, the mixing parameters $|q/p|_{{{B}^{0}}}$ and $|q/p|_{{{B}^{0}_{s}}}$, the average decay widths $\Gamma_{{d}}$ and $\Gamma_{{s}}$, and the width difference $\Delta\Gamma_{{s}}$ are fixed to the central values of the measurements reported in Table 1. The width difference $\Delta\Gamma_{{d}}$ is fixed to zero.

As already mentioned, for small values of $A_{{C\!P}}$ and $A_{\rm D}$, the decay rate is to first order only sensitive to the sum of these two quantities. For this reason, $A_{{C\!P}}$ is fixed to zero and $A_{\rm D}$ is left as a free parameter in the fits and hence measured from data, oppositely to the ${{B}^{+}}$ case, where an external input is necessary for $A_{\rm D}$. It is verified that the choice of different $A_{{C\!P}}$ values, up to the few percent level, leads to negligible variations of $A_{\rm P}$. Figures 2 and 3 show the ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}$, ${{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}}{{\pi}^{-}}$ and ${{D}^{-}_{s}}{{\pi}^{+}}$ invariant mass distributions and, in the case of the neutral $B$ meson, the time distributions with the results of the global fits overlaid, for data recorded at centre-of-mass energies of 7 and 8$\mathrm{\,Te\kern-1.00006ptV}$.

Figure 4 shows the raw asymmetries for neutral $B$-meson decays, defined as the ratio between the difference and the sum of the overall decay-time distributions, as a function of the decay time for candidates in the signal invariant mass regions, defined as the ranges 5250\text{--}5310$${\mathrm{\,Me\kern-1.00006ptV\!/}c^{2}} for ${B}^{0}$ decay and 5290\text{--}5450$${\mathrm{\,Me\kern-1.00006ptV\!/}c^{2}} for ${B}^{0}_{s}$ decay. The results of the global fits are overlaid.

The signal yields, $A_{\rm P}$ values and detection asymmetries obtained from the global fits are reported in Table 2, for the neutral $B$-meson decays, while the signal yield and $A_{\rm raw}$ for the ${{B}^{+}}$ decay are reported in Table 3. The $A_{\rm P}$ values obtained from the time-dependent global fits, reported here for illustrative purposes, are detector-independent quantities only if efficiency corrections as a function of $p_{\rm T}$ and $y$ are applied. An accurate knowledge of the decay-time resolution is important for the ${{B}^{0}_{s}}\rightarrow{{D}^{-}_{s}}{{\pi}^{+}}$ decay, due to the fast oscillation of the ${B}^{0}_{s}$ meson. For this reason the decay-time resolution is determined using the method previously described, applied to candidates in each ($p_{\rm T}$, $y$) bin.

According to Eq. 2, the measurement of $A_{\rm P}({{{B}^{+}}})$ requires knowledge of the $C\!P$ asymmetry $A_{{C\!P}}({{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}})$ and $A_{\rm D}({{K}^{+}})$. The value recently measured by LHCb [24] with an independent data set is used for the former and corresponds to $A_{{C\!P}}({{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}})=\left(0.09\pm 0.27\,\mathrm{\,(stat)}\pm 0.07\,\mathrm{\,(syst)}\right)\times 10^{-2}$. The measurement of the kaon detection asymmetry is obtained from ${D}$-meson decays produced directly in proton-proton collisions, using the same technique reported in Ref. [20]. It consists of measuring raw asymmetries from the two decay modes, ${{D}^{+}}\rightarrow{{K}^{-}}{{\pi}^{+}}{{\pi}^{+}}$ and ${{D}^{+}}\rightarrow{{K}^{0}_{\rm\scriptscriptstyle S}}{{\pi}^{+}}$ with ${{K}^{0}_{\rm\scriptscriptstyle S}}\rightarrow{{\pi}^{+}}{{\pi}^{-}}$, to obtain the ${K}^{+}$ ${\pi}^{-}$ detection asymmetry, $A_{\rm D}({{K}^{+}}{{\pi}^{-}})$, in each ($p_{\rm T}$, $y$) bin of the ${{B}^{+}}$ mesons. Additionally, $A_{\rm D}({{K}^{+}})$ is obtained by subtracting from $A_{\rm D}({{K}^{+}}{{\pi}^{-}})$ the pion detection asymmetry, $A_{\rm D}({{\pi}^{-}})$, measured by means of a sample of partially and fully reconstructed $D^{*+}\rightarrow{{D}^{0}}({{K}^{-}}{{\pi}^{+}}{{\pi}^{-}}{{\pi}^{+}}){{\pi}^{+}}$ decays, as described in Ref. [8]. It is estimated that the pion detection asymmetries across the various ${{B}^{+}}$ meson bins of ($p_{\rm T}$, $y$) are in the range $0\text{--}0.2$ %. Finally, the detection asymmetry of the ${{K}^{0}}\rightarrow{{\pi}^{+}}{{\pi}^{-}}$ final state, $A_{\mathrm{D}}({{K}^{0}})$, measured by LHCb to be $A_{\mathrm{D}}({{K}^{0}})=(0.054\pm 0.014)\%$ [20], is also subtracted.

7 Systematic uncertainties

Several sources of systematic uncertainty are considered. They are evaluated for each kinematic bin and for each decay mode. For the invariant mass model, the effects of the uncertainty on the shapes of all components (signals, combinatorial and partially reconstructed backgrounds) are investigated. For the decay-time model, systematic effects related to the decay-time resolution and reconstruction efficiency are studied. The effects of the uncertainties on the external inputs used in the fits, reported in Table 1, are evaluated by repeating the fits with each parameter varied by $\pm 1$ standard deviation ($\sigma$). Alternative decay-time parameterizations of the background components are also considered. To estimate the contribution of each single source, the fit is repeated for each ($p_{\rm T}$, $y$) bin after having modified the baseline fit model. The shifts from the relevant baseline values are taken as the systematic uncertainties. A detailed description follows. To estimate a systematic uncertainty related to the parameterization of final-state radiation effects on the signal mass distributions, the parameter $s$ of Eq. 6 is varied by $\pm 1\sigma$ of the corresponding value obtained from fits to simulated decays. A systematic uncertainty related to the invariant mass resolution model is estimated by repeating the fit using a simplified model with a single Gaussian function. The systematic uncertainty related to the parameterization of the mass distribution for the combinatorial background is investigated by replacing the exponential function with a linear function. Concerning the partially reconstructed background, a systematic uncertainty is assessed by repeating the fits while excluding the low invariant mass region, applying the requirement m>5330$${\mathrm{\,Me\kern-1.00006ptV\!/}c^{2}} to ${{B}^{0}_{s}}\rightarrow{{D}^{-}_{s}}{{\pi}^{+}}$ candidates. To estimate the uncertainty related to the parameterization of signal decay-time reconstruction efficiency, different functions are considered. Effects of inaccuracies in the knowledge of the decay-time resolution are estimated by rescaling the widths of the baseline model to obtain an average resolution width differing by $\pm 8$ fs. The impact of the small bias in the reconstructed decay time is assessed by introducing a corresponding bias of $\pm 2$ fs in the decay-time resolution model. The determination of the systematic uncertainties related to the $|q/p|$ input values requires special treatment, as $A_{\rm P}$ is correlated with $|q/p|$. For this reason, any variation of $|q/p|$ produces the same shift of $A_{\rm P}$ in each of the kinematic bins. Such a correlation is taken into account when integrating over $p_{\rm T}$ and $y$. The values of the systematic uncertainties related to the knowledge of $|q/p|$ are 0.0009 in the case of $A_{\rm P}({{B}^{0}})$ and 0.0021 in the case of $A_{\rm P}({{B}^{0}_{s}})$. For the ${{B}^{+}}$ decay, the uncertainties on $A_{{C\!P}}({{{B}^{+}}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}})$ and $A_{\rm D}({{K}^{+}})$ are considered as systematic uncertainties. They introduce correlations among the bins that are considered when the integrated results are calculated.

The ${\mathchar 28931\relax}^{0}_{b}$ production asymmetry is calculated, in each kinematic bin, assuming that the number of produced hadrons of any species in the $i$-th bin containing a ${b}$ quark, $N_{i,\,{b}}$, is equal to the number of produced hadrons containing a ${\overline{b}}$ quark, $N_{i,\,{\overline{b}}}$, i.e. relying on Eq. 3. This assumption is strictly valid in the full phase space, but not necessarily in a specific bin. In the event that $N_{i,\,{b}}\neq N_{i,\,{\overline{b}}}$, $A_{\rm P}({{\mathchar 28931\relax}^{0}_{b}})$ is biased by the quantity

[TABLE]

Values for $\delta z_{i}$ are studied using simulated events. Systematic uncertainties on $A_{\rm P}({{\mathchar 28931\relax}^{0}_{b}})$, in each kinematic bin, are assigned as half of the maximum variation from zero of the quantities $\delta z_{i}\pm\sigma(\delta z_{i})$, where $\sigma(\delta z_{i})$ is the related uncertainty.

The term $(f_{c}/f_{{{\mathchar 28931\relax}^{0}_{b}}})\cdot A_{\rm P}({{B}_{c}^{+}})$, estimated to be $3\cdot 10^{-5}$, can be safely neglected, while a systematic uncertainty related to neglecting the term $(f_{\rm other}/f_{{{\mathchar 28931\relax}^{0}_{b}}})\cdot A_{\rm P}(\rm other)$ has to be assessed. Amongst all other ${b}$ baryons, the production rate of $\Xi_{b}$ baryons is estimated from the simulation (which well reproduces the ${{B}^{+}}$, ${B}^{0}$, ${B}^{0}_{s}$ and ${\mathchar 28931\relax}^{0}_{b}$ fragmentation fractions) to be dominant, corresponding to about 1% of all $b$-hadron species produced in the primary collisions. On this basis, the neglected term can be evaluated as

[TABLE]

The value of $A_{\rm P}(\Xi_{b})$ is found to be double that of $A_{\rm P}({{\mathchar 28931\relax}^{0}_{b}})$ in the simulation. A systematic uncertainty on $A_{\rm P}({{\mathchar 28931\relax}^{0}_{b}})$ is obtained by assuming $A_{\rm P}(\Xi_{b})=2\,A_{\rm P}({{\mathchar 28931\relax}^{0}_{b}})$.

The dominant systematic uncertainties for the ${{B}^{+}}$ and ${B}^{0}$ cases are related to the measured value of $A_{{C\!P}}({{B}^{+}}\rightarrow{{J\mskip-3.0mu/\mskip-2.0mu\psi\mskip 2.0mu}}{{K}^{+}})$ and to $|q/p|_{{{B}^{0}}}$, respectively. The systematic uncertainty associated with the signal mass shape is the main source for the ${{B}^{0}_{s}}$ case, while it is the one related to neglecting the term $f_{\rm other}/f_{{{\mathchar 28931\relax}^{0}_{b}}}\cdot A_{\rm P}(\rm other)$ in Eq. 3 in the case of the ${{\mathchar 28931\relax}^{0}_{b}}$ decay. All the systematic uncertainties are summed in quadrature for each kinematic bin. Their values are reported, together with the final measurements, in Tables LABEL:tab:resultsBpB02011–9 in the Appendix.

When the integrated results are calculated, all the systematic uncertainties estimated for each bin are propagated according to Eq. 4 and correlations among the bins are taken into account. An additional systematic uncertainty is considered by studying how the integrated values vary in the case that the values of $\omega_{i}$ are measured using a data driven approach. In this case $\omega_{i}^{\mathrm{data}}$ is measured as

[TABLE]

where $\varepsilon_{i}^{\mathrm{total}}$ is the total reconstruction efficiency, obtained as a combination of the selection efficiency, determined from simulation, and PID and trigger efficiencies, measured from data. Differences in the central values between $A_{\mathrm{P}}$ calculated using either $\omega_{i}$ or $\omega_{i}^{\mathrm{data}}$ are found to be small for all the decay modes. Table 4 summarizes systematic uncertainties associated with the integrated measurements.

8 Results and conclusions

Using a data sample corresponding to an integrated luminosity of 3.0$\mbox{\,fb}^{-1}$, the ${{B}^{+}}$, ${B}^{0}$ and ${B}^{0}_{s}$ hadron production asymmetries have been determined independently for each ($p_{\rm T}$,$y$) bin and then combined using Eq. 3 to derive the ${\mathchar 28931\relax}^{0}_{b}$ production asymmetry. Tables LABEL:tab:resultsBpB02011–9, in the Appendix, report the final results.

The ${{B}^{+}}$, ${B}^{0}$, ${B}^{0}_{s}$ and ${\mathchar 28931\relax}^{0}_{b}$ hadron production asymmetries are also determined integrating over $p_{\rm T}$ or $y$, in the range 0<\mbox{p_{\rm T}}<30$${\mathrm{\,Ge\kern-1.00006ptV\!/}c} and $2.1<y<4.5$ for ${{B}^{+}}$ and ${B}^{0}$ decays, and in the range 2<\mbox{p_{\rm T}}<30$${\mathrm{\,Ge\kern-1.00006ptV\!/}c} and $2.1<y<4.5$ for ${B}^{0}_{s}$ and ${\mathchar 28931\relax}^{0}_{b}$ decay. The corresponding numerical values are reported in Tables 10–17, in the Appendix, and in Figs. 5–8, where the results of the fits with a constant and a first-order polynomial function are also shown. Table 5 reports the values of the fit parameters. No evidence for any dependence is observed. Finally, integrating over both $p_{\rm T}$ and $y$, the $b$-hadron production asymmetries are found to be

[TABLE]

All the results are consistent with zero within $2.5$ standard deviations. The results of this analysis supersede the previous LHCb results of Ref. [6]. These measurements, once integrated using appropriate weights for any reconstructed ${{B}^{+}}$, ${B}^{0}$, ${B}^{0}$, ${\mathchar 28931\relax}^{0}_{b}$ decay in LHCb, can be used to determine effective production asymmetries, as inputs for $C\!P$ violation measurements with the LHCb data.

Acknowledgements

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 (The Netherlands); MNiSW and NCN (Poland); MEN/IFA (Romania); MinES and FASO (Russia); MinECo (Spain); SNSF and SER (Switzerland); NASU (Ukraine); STFC (United Kingdom); NSF (USA). We acknowledge the computing resources that are provided by CERN, IN2P3 (France), KIT and DESY (Germany), INFN (Italy), SURF (The 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), Conseil Général de Haute-Savoie, Labex ENIGMASS and OCEVU, Région Auvergne (France), RFBR and Yandex LLC (Russia), GVA, XuntaGal and GENCAT (Spain), Herchel Smith Fund, The Royal Society, Royal Commission for the Exhibition of 1851 and the Leverhulme Trust (United Kingdom).

Appendix