Example of $\Lambda$ signal for 0-40% most central events, over mixed-event background for the bin $-0.05 < y_{cm} < 0.05$ and reduced transverse masses between $100-150 MeV/c^{2}$.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $K^{0}_{S}$ for the 0-40% most central events.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $K^{0}_{S}$ for the 0-40% most central events. NOTE: The spectra are not scaled by $1/N_{Events}$! To compare the data, divide by $N_{Events} = 2.1997626 x 10^{9}$

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $\Lambda$ for the 0-40% most central events.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $\Lambda$ for the 0-40% most central events. NOTE: The spectra are not scaled by $1/N_{Events}$! To compare the data, divide by $N_{Events} = 2.1997626 x 10^{9}$

Rapidity distribution of $K^{0}_{S}$

Rapidity distribution of $K^{0}_{S}$

Rapidity distribution of $\Lambda$

Rapidity distribution of $\Lambda$

Compilation of mid-rapidity yields for central Au+Au collisions as a function of $\sqrt{s_{NN}}$ of $K^{0}_{S}$ and $\Lambda$.

Compilation of mid-rapidity yields for central Au+Au collisions as a function of $\sqrt{s_{NN}}$ of $K^{0}_{S}$ and $\Lambda$.

Multiplicities per mean number of participants $Mult/\langle A_{Part}\rangle$ as a function of $\langle A_{Part}\rangle$

Multiplicities per mean number of participants $Mult/\langle A_{Part}\rangle$ as a function of $\langle A_{Part}\rangle$

Multiplicities per mean number of participants $Mult/\langle A_{Part}\rangle$ as a function of $\langle A_{Part}\rangle$ for $K^{0}_{S}$

Multiplicities per mean number of participants $Mult/\langle A_{Part}\rangle$ as a function of $\langle A_{Part}\rangle$ for $K^{0}_{S}$

Multiplicities per mean number of participants $Mult/\langle A_{Part}\rangle$ as a function of $\langle A_{Part}\rangle$ for $\Lambda$

Multiplicities per mean number of participants $Mult/\langle A_{Part}\rangle$ as a function of $\langle A_{Part}\rangle$ for $\Lambda$

Comparison of the shape of the rapidity distribution of $K^{0}_{S}$ to various transport model versions.

Comparison of the shape of the rapidity distribution of $K^{0}_{S}$ to various transport model versions.

Comparison of the shape of the rapidity distribution of $\Lambda$ to various transport model versions.

Comparison of the shape of the rapidity distribution of $\Lambda$ to various transport model versions.

Comparison of the shape of the $p_{t}$-spectra for $\pm0.15$ rapidity units around mid-rapidity of $K^{0}_{S}$ to various transport model versions.

Comparison of the shape of the $p_{t}$-spectra for $\pm0.15$ rapidity units around mid-rapidity of $K^{0}_{S}$ to various transport model versions.

Comparison of the shape of the $p_{t}$-spectra for $\pm0.15$ rapidity units around mid-rapidity of $\Lambda$ to various transport model versions.

Comparison of the shape of the $p_{t}$-spectra for $\pm0.15$ rapidity units around mid-rapidity of $\Lambda$ to various transport model versions.

Differential yield of $K^{0}_{S}$ as function of rapdidity and reduced transverse mass for the four centrality classes.

Differential yield of $K^{0}_{S}$ as function of rapdidity and reduced transverse mass for the four centrality classes. NOTE: The spectra are not scaled by $1/N_{Events}$! To compare the data, divide by $N_{Events} = 5.64247975 x 10^{8} (0 - 10\%)$, $5.85743394 x 10^{8} (10 - 20\%)$, $5.76369451 x 10^{8} (20 - 30\%)$ or $4.73401752 x 10^{8} (30 - 40\%)$

Differential yield of $\Lambda$ as function of rapdidity and reduced transverse mass for the four centrality classes.

Differential yield of $\Lambda$ as function of rapdidity and reduced transverse mass for the four centrality classes. NOTE: The spectra are not scaled by $1/N_{Events}$! To compare the data, divide by $N_{Events} = 5.64247975 x 10^{8} (0 - 10\%)$, $5.85743394 x 10^{8} (10 - 20\%)$, $5.76369451 x 10^{8} (20 - 30\%)$ or $4.73401752 x 10^{8} (30 - 40\%)$

$K^{0}_{S}$ and $\Lambda$ multiplicities in full phase space and inverse slopes at mid-rapidity $T_{Eff}$ for a given centrality.

$K^{0}_{S}$ and $\Lambda$ multiplicities in full phase space and inverse slopes at mid-rapidity $T_{Eff}$ for a given centrality.

Value of the parameter $\alpha$ extracted for $K^{+}$, $K^{-}$, $K^{0}_{S}$, $\Lambda$ and $\Phi$, as displayed in Figure 5 - Global $Mult/\langle A_{Part}\rangle$ Fit

Value of the parameter $\alpha$ extracted for $K^{+}$, $K^{-}$, $K^{0}_{S}$, $\Lambda$ and $\Phi$, as displayed in Figure 5 - Global $Mult/\langle A_{Part}\rangle$ Fit

$K^{+}$ signal and the corresponding background fit for the region covering mid-rapidity and $m_{t}−m_{0}$ between 25 and 50 $MeV/c^{2}$.

$K^{+}$ signal and the corresponding background fit for the region covering mid-rapidity and $m_{t}−m_{0}$ between 25 and 50 $MeV/c^{2}$.

$K^{+}$ signal and the corresponding background fit for the region covering mid-rapidity and $m_{t}−m_{0}$ between 25 and 50 $MeV/c^{2}$.

$K^{-}$ signal and the corresponding background fit for the region covering mid-rapidity and $m_{t}−m_{0}$ between 50 and 75 $MeV/c^{2}$.

$K^{-}$ signal and the corresponding background fit for the region covering mid-rapidity and $m_{t}−m_{0}$ between 50 and 75 $MeV/c^{2}$.

$K^{-}$ signal and the corresponding background fit for the region covering mid-rapidity and $m_{t}−m_{0}$ between 50 and 75 $MeV/c^{2}$.

$K^{+}K^{-}$ invariant mass distribution for the mid-rapidity region and $m_{t}−m_{0}$ between 0 and 100 $MeV/c^{2}$ after subtraction of the background.

$K^{+}K^{-}$ invariant mass distribution for the mid-rapidity region and $m_{t}−m_{0}$ between 0 and 100 $MeV/c^{2}$ after subtraction of the background.

$K^{+}K^{-}$ invariant mass distribution for the mid-rapidity region and $m_{t}−m_{0}$ between 0 and 100 $MeV/c^{2}$ after subtraction of the background.

Rapidity distributions of $K^{-}$ and $\phi$ as well as extracted inverse slope parameters obtained from the Boltzmann fits to the the $m_{t}$ spectra.

Rapidity distributions of $K^{-}$ and $\phi$ as well as extracted inverse slope parameters obtained from the Boltzmann fits to the the $m_{t}$ spectra.

Rapidity distributions of $K^{-}$ and $\phi$ as well as extracted inverse slope parameters obtained from the Boltzmann fits to the the $m_{t}$ spectra.

$K^{-}$ transverse-mass spectrum around mid-rapidity

$K^{-}$ transverse-mass spectrum around mid-rapidity

$K^{-}$ transverse-mass spectrum around mid-rapidity

Multiplicity ratio $\phi/K^{-}$ as function of $\sqrt{s_{NN}}$ for central heavy-ion collisions.

Multiplicity ratio $\phi/K^{-}$ as function of $\sqrt{s_{NN}}$ for central heavy-ion collisions.

Multiplicity ratio $\phi/K^{-}$ as function of $\sqrt{s_{NN}}$ for central heavy-ion collisions.

Multiplicities and effective inverse slopes $T_{eff}$ at mid-rapidity for given centrality classes.

Multiplicities and effective inverse slopes $T_{eff}$ at mid-rapidity for given centrality classes.

Multiplicities and effective inverse slopes $T_{eff}$ at mid-rapidity for given centrality classes.

Multiplicity ratios for given centrality classes.

Multiplicity ratios for given centrality classes.

Multiplicity ratios for given centrality classes.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $K^{+}$ for the 0-40% most central events.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $K^{+}$ for the 0-40% most central events.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $K^{+}$ for the 0-40% most central events.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $K^{-}$ for the 0-40% most central events.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $K^{-}$ for the 0-40% most central events.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $K^{-}$ for the 0-40% most central events.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $\phi$ for the 0-40% most central events.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $\phi$ for the 0-40% most central events.

Reduced transverse mass ($m_{t}-m_{0}$) spectra of $\phi$ for the 0-40% most central events.

Rapidity distribution of $K^{+}$ for the 0-40% most central events.

Differential production cross-sections for prompt $D^{+} + D^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

Differential production cross-sections for prompt $D_{s}^{+} + D_{s}^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

Differential production cross-sections for prompt $D_{s}^{+} + D_{s}^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

Differential production cross-sections for prompt $D^{*+} + D^{*-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

Differential production cross-sections for prompt $D^{*+} + D^{*-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D^{0} + \bar{D}^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D^{0} + \bar{D}^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D^{+} + D^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D^{+} + D^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D_{s}^{+} + D_{s}^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D_{s}^{+} + D_{s}^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D^{*+} + D^{*-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D^{*+} + D^{*-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{*+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{*+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{*+}$ and $D^{+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{*+}$ and $D^{+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{*+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{*+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

Optimised angular observables evaluated by the unbinned maximum likelihood fit. The first uncertainties are statistical and the second systematic.

CP-averaged angular observables evaluated using the method of moments. The first uncertainties are statistical and the second systematic.

CP-asymmetries evaluated using the method of moments. The first uncertainties are statistical and the second systematic.

Optimised observables evaluated using the method of moments. The first uncertainties are statistical and the second systematic.

Zero-crossing points determined with an amplitude fit.

Likelihood correlation matrix $0.1 < q^2 < 0.98~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $1.1 < q^2 < 2.5~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $2.5 < q^2 < 4.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $4.0 <q^2< 6.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $6.0 < q^2 < 8.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $11.0 <q^2< 12.5 ~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $15.0 < q^2 < 17.0 ~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $17.0 <q^2< 19.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $1.1 <q^2< 6.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $15.0 <q^2< 19.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $0.1 < q^2 < 0.98~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $1.1 < q^2 < 2.5~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $2.5 < q^2 < 4.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $4.0 <q^2< 6.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $6.0 < q^2 < 8.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $11.0 <q^2< 12.5 ~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $15.0 <q^2< 17.0 ~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $17.0 <q^2< 19.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $1.1 <q^2< 6.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $15.0 <q^2< 19.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $0.1 < q^2 < 0.98~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $1.1 < q^2 < 2.5~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $2.5 < q^2 < 4.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $4.0 <q^2< 6.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $6.0 < q^2 < 8.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $11.0 <q^2< 12.5 ~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $15.0 <q^2< 17.0 ~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $17.0 <q^2< 19.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $1.1 <q^2< 6.0~{\rm GeV}^2/c^4$.

Likelihood correlation matrix $15.0 <q^2< 19.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $0.10 < q^2 < 0.98~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $1.1 < q^2 < 2.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $2.0 < q^2 < 3.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $3.0 < q^2 < 4.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $4.0 < q^2 < 5.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $5.0 < q^2 < 6.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $6.0 < q^2 < 7.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $7.0 < q^2 < 8.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $11.00 <q^2 < 11.75~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $11.75 <q^2 < 12.50~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $15.0 <q^2 < 16.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $16.0 <q^2 < 17.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $17.0 <q^2 < 18.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $18.0 <q^2 < 19.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $15.0 <q^2 < 19.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $0.10 < q^2 < 0.98~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $1.1 < q^2 < 2.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $2.0 < q^2 < 3.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $3.0 < q^2 < 4.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $4.0 < q^2 < 5.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $5.0 < q^2 < 6.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $6.0 < q^2 < 7.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $7.0 < q^2 < 8.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $11.00 <q^2 < 11.75~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $11.75 <q^2 < 12.50~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $15.0 <q^2 < 16.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $16.0 <q^2 < 17.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $17.0 <q^2 < 18.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $18.0 <q^2 < 19.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $15.0 <q^2 < 19.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $0.1 <q^2 < 0.98~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $1.1 <q^2 < 2.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $2.0 <q^2 < 3.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $3.0 <q^2 < 4.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $4.0 <q^2 < 5.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $5.0 <q^2 < 6.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $6.0 <q^2 < 7.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $7.0 <q^2 < 8.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $11.0 <q^2 < 11.75~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $11.75 <q^2 < 12.5~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $15.0 <q^2 < 16.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $16.0 <q^2 < 17.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $17.0 < q^2 < 18.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $18.0 < q^2 < 19.0~{\rm GeV}^2/c^4$.

Bootstrap correlation matrix $15.0 < q^2 < 19.0~{\rm GeV}^2/c^4$.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} p_T(D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(\Upsilon(1S))}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(\Upsilon(1S))}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(D^0)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{\pi}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \left| \Delta \phi \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{\pi}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \left| \Delta \phi \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \left| \Delta y \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \left| \Delta y \right| }$ for $2<y(\Upsilon(1S))<4.5$, for $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} p_T(\Upsilon(1S)D^0) }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} p_T(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(\Upsilon(1S)D^0) }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \mathcal{A}_T }$, where $\mathcal{A}_T = \frac{p_T(\Upsilon(1S))-p_T(D^0)}{p_T(\Upsilon(1S))-p_T(D^0)}$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \mathcal{A}_T }$, where $\mathcal{A}_T = \frac{p_T(\Upsilon(1S))-p_T(D^0)}{p_T(\Upsilon(1S))-p_T(D^0)}$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} m(\Upsilon(1S)D^0) }$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} m(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Double-differential yields, $\frac{d^{2}n}{x_{_F}p_{_T}}$ and $f_n(x_{_F},p_{T})$, for $x_{_F}>0$.

$p_{_T}$ integrated yield $\frac{dn}{dy}$, the inverse slope parameter $T$ and the mean transverse mass $\langle m_{_T}\rangle - m_{\Lambda}$.

$p_{_T}$ integrated yield $\frac{dn}{x_{_F}}$ and the invariant cross section $F(x_{_F})$.