Predicted and observed yields as a function of $m_{4\ell}$ in the VBS-Suppressed region. Overflow events are included in the last bin of the distribution.

Relative systematic uncertainties in the unfolded differential cross section for 4ljj production in the VBS-Enhanced region as a function of m_{jj}

Relative systematic uncertainties in the unfolded differential cross section for 4ljj production in the VBS-Enhanced region as a function of m_{4\ell}

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $m_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $m_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $m_{4\ell}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $m_{4\ell}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $ p_{T,4l}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $ p_{T,4l}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $\Delta\phi_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $\Delta\phi_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $|\Delta y_{jj}|$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $|\Delta y_{jj}|$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $cos(\theta^{*}_{12})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $cos(\theta^{*}_{12})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $cos(\theta^{*}_{34})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $cos(\theta^{*}_{34})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $p_{T,jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $p_{T,jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $p_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $p_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $S_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $S_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Statistical correlation between different bins of the unfolded cross-section for all the measured observables in the VBS-Enhanced region.

Statistical correlation between different bins of the unfolded cross-section for all the measured observables in the VBS-Suppressed region.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the signal region, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $W$+jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Expected distributions of the MEM discriminant $P(S|X)$ in the SR, for the s-channel single-top signal, and for the $t\bar{t}$ and $W$+jets backgrounds, for MEM discriminant values larger than $2.0\times10^{-4}$. Each distribution is normalised to unity. The binning is the same as the optimised binning used in the signal extraction fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $W$+jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 3-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 4-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty bands indicate the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR before the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore not included. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$, after subtraction of all backgrounds. The fitted distribution for the simulation of the signal is shown together with the post-fit uncertainty in the backgrounds. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Pre-fit and post-fit event yields in the SR, for MEM discriminant values larger than $2.0\times10^{-4}$. The central value of the event yield for each process is calculated by summing the values of the discriminant bin contents, using the nominal expected yield for the pre-fit value, and the best-fit estimate for the post-fit value. The error includes statistical and systematic uncertainties summed in quadrature. All sources of systematic uncertainties are included, taking into account correlations and anti-correlations in the post-fit case. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore only included in the post-fit uncertainties.

Observed impact of the different sources of uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. The 'Systematic uncertainties' category combines all sources of systematic uncertainties. The statistical uncertainty is obtained by repeating the fit after having fixed all nuisance parameters, including the $t\bar{t}$ and $W$+jets normalisation factors. 'Total' gives the total uncertainty on the measurement.

Observed impact of the different sources of $t\bar{t}$ modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and 'ME/PS matching' to the matching of the ME to the parton shower.

Observed impact of the different sources of s-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of t-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of $tW$ modelling uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and '$t\bar{t}$ overlap' to the algorithm removing the overlap between $tW$ and $t\bar{t}$ production at NLO, as described in Section 7 in the paper.

Observed impact of the different sources of PDF uncertainties on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the ${E}_{T}^{miss}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the $m_{T}^{W}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Nuisance parameters ranked according to their post-fit impacts on the best-fit value of the ratio $\mu$ of the measured cross-section to the predicted cross-section. In the figure, only the 20 nuisance parameters with the largest post-fit impacts are shown. The empty (solid) blue rectangles illustrate the pre-fit (post-fit) impact on $\mu$, corresponding to the upper axis. The pre-fit (post-fit) impact of each nuisance parameter, $\Delta\mu$, is calculated as the difference in the fitted value of $\mu$ between the nominal fit and the fit when fixing the corresponding nuisance parameter to $\hat{\theta}\pm\Delta\theta$ ($\hat{\theta}\pm\Delta\hat{\theta}$), where $\hat{\theta}$ is the best-fit value of the nuisance parameter and $\Delta\theta$ ($\Delta\hat{\theta}$) is its pre-fit (post-fit) uncertainty. Several systematic uncertainties are split into different nuisance parameters, which are indicated by NP. JES (JER) indicates jet energy scale (resolution), and $\gamma$ indicates a nuisance parameter associated to the MC statistics in one of the 18 bins numbered from 0 to 17. The black points show the best-fit values of the nuisance parameters, with the error bars representing the post-fit uncertainties. Each nuisance parameter is shown wrt. its nominal value, $\theta_0$, and in units of its pre-fit uncertainty, except the free-floating normalisation factors of the $t\bar{t}$ and $W$+jets backgrounds, and the parameters associated to the MC statistics in each bin, for which the post-fit values and uncertainties are shown.

Correlation matrix of the nuisance parameters and of the ratio $\mu$ of the measured cross-section to the predicted cross-section. The correlations are given after the fit to data. In the figure, only the parameters which have a correlation of at least 0.2 with any other parameter are shown.

Distribution of the MEM discriminant $P(S|X)$ in the SR for MEM discriminant values larger than $2.0\times10^{-4}$, for the collision data used for the measurement, and for 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

Measured values of the signal cross-section and of the $t\bar{t}$ and $W$+jets normalisation factors, obtained by statistical-only fits to the collision data used for the measurement, and to 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The central values and their statistical uncertainties are obtained by repeating the fit after having fixed all nuisance parameters, except the $t\bar{t}$ and $W$+jets normalisation factors, which are let free-floating (unlike for the statistical uncertainty on the cross-section quoted in the paper). The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$.

Signal region detector-level distribution for the observable $m_{e\mu}$.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{e\mu}$.

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

Parton dijet $M_{inv}$ vs $A_{LL}$ values with associated uncertainties, for topology C.

Parton dijet $M_{inv}$ vs $A_{LL}$ values with associated uncertainties, for topology D.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements.The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology A). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology B). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology A). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology A) with forward-central dijet measurements (topology B). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology A) with forward-central dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology A) with forward-central dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology B). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology B) with forward-central dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology B) with forward-central dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology C) with forward-central dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

Parton inclusive-jet $p_T$ and $A_{LL}$ values with associated uncertainties for jet-$\eta$ region $0.5<|\eta|<1$.

Parton inclusive-jet $p_T$ and $A_{LL}$ values with associated uncertainties for jet-$\eta$ region $|\eta|<0.5$.

Parton dijet invariant mass $M_{inv}$ and $A_{LL}$ values with associated uncertainties for the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology.

Parton dijet invariant mass $M_{inv}$ and $A_{LL}$ values with associated uncertainties for the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology.

Parton inclusive-jet $p_T$, $x_T$, and $A_{LL}$ values with associated uncertainties for jet-$\eta$ region $|\eta|<1$.

The correlation matrix for the point-to-point uncertainties for inclusive jet measurements with jets in the $0.5<|\eta|<1$ region. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties for inclusive jet measurements with jets in the $|\eta|<0.5$ region. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties for inclusive jet measurements with jets in the $|\eta|<1$ region. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<0.5$ region and jets in the $0.5<|\eta|<1$ region. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $0.5<|\eta|<1$ region with dijet measurements with the the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $0.5<|\eta|<1$ region with dijet measurements with the the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<0.5$ region with dijet measurements with the the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<0.5$ region with dijet measurements with the the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<1$ region with dijet measurements with the the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<1$ region with dijet measurements with the the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties for dijet measurements with the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties for dijet measurements with the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling dijet measurements with the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology and the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

Out projection of raw 3D LCMS KS KS correlation function for 1.0 < kT < 1.5 GeV/c bin

Side projection of raw 3D LCMS KS KS correlation function for 1.0 < kT < 1.5 GeV/c bin

Long projection of raw 3D LCMS KS KS correlation function for 1.0 < kT < 1.5 GeV/c bin

Rout vs. mT for PI+- PI+- for centrality 0-10%

Rside vs. mT for PI+- PI+- for centrality 0-10%

Rlong vs. mT for PI+- PI+- for centrality 0-10%

Rout vs. mT for PI+- PI+- for centrality 10-30%

Rside vs. mT for PI+- PI+- for centrality 10-30%

Rlong vs. mT for PI+- PI+- for centrality 10-30%

Rout vs. mT for PI+- PI+- for centrality 30-50%

Rside vs. mT for PI+- PI+- for centrality 30-50%

Rlong vs. mT for PI+- PI+- for centrality 30-50%

Rout vs. mT for PI+- PI+- for centrality 0-5%

Rside vs. mT for PI+- PI+- for centrality 0-5%

Rlong vs. mT for PI+- PI+- for centrality 0-5%

Ratio Rout/Rside vs. mT for PI+- PI+- for centrality 0-5%

Ratio Rout/Rside vs. mT for PI+- PI+- for centrality 0-10%

Ratio Rout/Rside vs. mT for PI+- PI+- for centrality 10-30%

Ratio Rout/Rside vs. mT for PI+- PI+- for centrality 30-50%

Rlong^2 vs. mT for PI+- PI+- for centrality 0-5%

C(k*)/(linear baseline fit) for K0s K+ and all kT bin

C(k*)/(linear baseline fit) for K0s K+ and kT < 0.675 GeV/c bin

C(k*)/(linear baseline fit) for K0s K+ and kT > 0.675 GeV/c bin

C(k*)/(linear baseline fit) for K0s K- and all kT bin

C(k*)/(linear baseline fit) for K0s K- and kT < 0.675 GeV/c bin

C(k*)/(linear baseline fit) for K0s K- and kT > 0.675 GeV/c bin

R vs. kT for K0s K- with Achasov1 parameters

lambda vs. kT for K0s K- with Achasov1 parameters

R vs. kT for K0s K+- with Achasov1 parameters and linear fit

lambda vs. kT for K0s K+- with Achasov1 parameters and linear fit

The azimuthal dependence of $R_{out}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20-30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $R_{side}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $R_{side}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $R_{side}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $R_{side}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $R_{long}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $R_{long}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $R_{long}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $R_{long}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $\lambda$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $\lambda$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $\lambda$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $\lambda$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $R_{os}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $R_{os}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $R_{os}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The azimuthal dependence of $R_{os}^{2}$ as a function of $\Delta\varphi=\varphi_{\mathrm{pair}}-\Psi_{\mathrm EP,2}$ for the centrality 20--30% and different $k_{\mathrm{T}}$ ranges.

The average radii $R_{out,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{out,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{out,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{out,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{side,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{side,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{side,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{side,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{long,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{long,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{long,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{long,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{os,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{os,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{os,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

The average radii $R_{os,0}^{2}$ as a function of centrality for different $k_{\mathrm{T}}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm out,2}/R^{2}_{\mathrm side,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm out,2}/R^{2}_{\mathrm side,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm out,2}/R^{2}_{\mathrm side,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm out,2}/R^{2}_{\mathrm side,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm side,2}/R^{2}_{\mathrm side,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm side,2}/R^{2}_{\mathrm side,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm side,2}/R^{2}_{\mathrm side,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm side,2}/R^{2}_{\mathrm side,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm long,2}/R^{2}_{\mathrm long,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm long,2}/R^{2}_{\mathrm long,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm long,2}/R^{2}_{\mathrm long,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm long,2}/R^{2}_{\mathrm long,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm os,2}/R^{2}_{\mathrm side,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm os,2}/R^{2}_{\mathrm side,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm os,2}/R^{2}_{\mathrm side,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

Amplitudes of the relative radius oscillations~$R^{2}_{\mathrm os,2}/R^{2}_{\mathrm side,0}$ versus centrality for different $k_{\mathrm T}$ ranges.

An estimate of freeze-out eccentricity $2 R^{2}_{\mathrm side,2}/R^{2}_{\mathrm side,0}$ for different $k_{\mathrm T}$ ranges vs. initial state eccentricity from Monte Carlo Glauber model for six centrality ranges, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, and 40-50%.

An estimate of freeze-out eccentricity $2 R^{2}_{\mathrm side,2}/R^{2}_{\mathrm side,0}$ for different $k_{\mathrm T}$ ranges vs. initial state eccentricity from Monte Carlo Glauber model for six centrality ranges, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, and 40-50%.

An estimate of freeze-out eccentricity $2 R^{2}_{\mathrm side,2}/R^{2}_{\mathrm side,0}$ for different $k_{\mathrm T}$ ranges vs. initial state eccentricity from Monte Carlo Glauber model for six centrality ranges, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, and 40-50%.

An estimate of freeze-out eccentricity $2 R^{2}_{\mathrm side,2}/R^{2}_{\mathrm side,0}$ for different $k_{\mathrm T}$ ranges vs. initial state eccentricity from Monte Carlo Glauber model for six centrality ranges, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, and 40-50%.

$\Delta\eta$ integrated projections of correlation functions for (a) $\pi^+\pi^-$, (b) $\rm K^+K^-$, (c) $\mathrm{p\overline{p}}$ pairs.

$\Delta\eta$ integrated projections of correlation functions for combined pairs of (a) $\pi^+\pi^+ + \pi^-\pi^-$, (b) $\rm K^+K^+ + K^-K^-$, (c) $\mathrm{pp} + \rm \overline{p}\overline{p}$ pairs.

$\Delta\eta$ integrated projections of correlation functions for (a) $\Lambda\Lambda+\overline{\Lambda}\overline{\Lambda}$, (b)$\Lambda\overline{\Lambda}$ pairs.

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $5 < p_{\rm T} < 8~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ ($-0.96 < y_{\rm cms} < 0.04$) and charged particles with $p_{\rm T} > 0.3~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp (p-Pb) collisions at $\sqrt{s} = 7~{\rm TeV}$ ($\sqrt{s_{\rm NN}} = 5.02~{\rm TeV}$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $8 < p_{\rm T} < 16~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ ($-0.96 < y_{\rm cms} < 0.04$) and charged particles with $p_{\rm T} > 0.3~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp (p-Pb) collisions at $\sqrt{s} = 7~{\rm TeV}$ ($\sqrt{s_{\rm NN}} = 5.02~{\rm TeV}$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $3 < p_{\rm T} < 5~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ and charged particles with $0.3 < p_{\rm T} < 1~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp collisions at $\sqrt{s} = 7~{\rm TeV}$.

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $5 < p_{\rm T} < 8~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ ($-0.96 < y_{\rm cms} < 0.04$) and charged particles with $0.3 < p_{\rm T} < 1~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp (p-Pb) collisions at $\sqrt{s} = 7~{\rm TeV}$ ($\sqrt{s_{\rm NN}} = 5.02~{\rm TeV}$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $8 < p_{\rm T} < 16~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ ($-0.96 < y_{\rm cms} < 0.04$) and charged particles with $0.3 < p_{\rm T} < 1~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp (p-Pb) collisions at $\sqrt{s} = 7~{\rm TeV}$ ($\sqrt{s_{\rm NN}} = 5.02~{\rm TeV}$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $3 < p_{\rm T} < 5~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ and charged particles with $p_{\rm T} > 1~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp collisions at $\sqrt{s} = 7~{\rm TeV}$.

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $5 < p_{\rm T} < 8~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ ($-0.96 < y_{\rm cms} < 0.04$) and charged particles with $p_{\rm T} > 1~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp (p-Pb) collisions at $\sqrt{s} = 7~{\rm TeV}$ ($\sqrt{s_{\rm NN}} = 5.02~{\rm TeV}$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $8 < p_{\rm T} < 16~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ ($-0.96 < y_{\rm cms} < 0.04$) and charged particles with $p_{\rm T} > 1~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp (p-Pb) collisions at $\sqrt{s} = 7~{\rm TeV}$ ($\sqrt{s_{\rm NN}} = 5.02~{\rm TeV}$).

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = 7$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = 7$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = 7$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = 7$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $p_{\rm T} > 1$ GeV/$c$, in pp collisions at $\sqrt{s} = 7$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $p_{\rm T} > 1$ GeV/$c$, in pp collisions at $\sqrt{s} = 7$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $3 < p_{\rm T} < 5~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ and charged particles with $p_{\rm T} > 0.3~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp collisions at $\sqrt{s} = 7~{\rm TeV}$. The points are shown after the subtraction of the baseline of the distribution, which is defined as the physical minimum of the distribution and estimated from the weighted average of the points in the transverse region ($\pi/4 < |\Delta\varphi| < \pi/2$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $5 < p_{\rm T} < 8~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ and charged particles with $p_{\rm T} > 0.3~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp collisions at $\sqrt{s} = 7~{\rm TeV}$. The points are shown after the subtraction of the baseline of the distribution, which is defined as the physical minimum of the distribution and estimated from the weighted average of the points in the transverse region ($\pi/4 < |\Delta\varphi| < \pi/2$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $8 < p_{\rm T} < 16~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ and charged particles with $p_{\rm T} > 0.3~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp collisions at $\sqrt{s} = 7~{\rm TeV}$. The points are shown after the subtraction of the baseline of the distribution, which is defined as the physical minimum of the distribution and estimated from the weighted average of the points in the transverse region ($\pi/4 < |\Delta\varphi| < \pi/2$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $3 < p_{\rm T} < 5~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ and charged particles with $0.3 < p_{\rm T} < 1~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp collisions at $\sqrt{s} = 7~{\rm TeV}$. The points are shown after the subtraction of the baseline of the distribution, which is defined as the physical minimum of the distribution and estimated from the weighted average of the points in the transverse region ($\pi/4 < |\Delta\varphi| < \pi/2$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $5 < p_{\rm T} < 8~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ and charged particles with $0.3 < p_{\rm T} < 1~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp collisions at $\sqrt{s} = 7~{\rm TeV}$. The points are shown after the subtraction of the baseline of the distribution, which is defined as the physical minimum of the distribution and estimated from the weighted average of the points in the transverse region ($\pi/4 < |\Delta\varphi| < \pi/2$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $8 < p_{\rm T} < 16~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ and charged particles with $0.3 < p_{\rm T} < 1~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp collisions at $\sqrt{s} = 7~{\rm TeV}$. The points are shown after the subtraction of the baseline of the distribution, which is defined as the physical minimum of the distribution and estimated from the weighted average of the points in the transverse region ($\pi/4 < |\Delta\varphi| < \pi/2$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $3 < p_{\rm T} < 5~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ and charged particles with $p_{\rm T} > 1~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp collisions at $\sqrt{s} = 7~{\rm TeV}$. The points are shown after the subtraction of the baseline of the distribution, which is defined as the physical minimum of the distribution and estimated from the weighted average of the points in the transverse region ($\pi/4 < |\Delta\varphi| < \pi/2$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $5 < p_{\rm T} < 8~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ and charged particles with $p_{\rm T} > 1~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp collisions at $\sqrt{s} = 7~{\rm TeV}$. The points are shown after the subtraction of the baseline of the distribution, which is defined as the physical minimum of the distribution and estimated from the weighted average of the points in the transverse region ($\pi/4 < |\Delta\varphi| < \pi/2$).

Azimuthal correlation of D mesons (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) with $8 < p_{\rm T} < 16~{\rm GeV}/c$ and $|y_{\rm cms}| < 0.5$ and charged particles with $p_{\rm T} > 1~{\rm GeV}/c$ for $|\Delta\eta| = |\eta_{\rm ch}-\eta_{\rm D}| < 1$ measured in pp collisions at $\sqrt{s} = 7~{\rm TeV}$. The points are shown after the subtraction of the baseline of the distribution, which is defined as the physical minimum of the distribution and estimated from the weighted average of the points in the transverse region ($\pi/4 < |\Delta\varphi| < \pi/2$).

Near-side-peak associated yield values measured in pp collisions at $\sqrt{s} = 7$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $p_{\rm T} > 0.3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side-peak associated yield values measured in pp collisions at $\sqrt{s} = 7$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side-peak associated yield values measured in pp collisions at $\sqrt{s} = 7$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $p_{\rm T} > 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side-peak width values measured in pp collisions at $\sqrt{s} = 7$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $p_{\rm T} > 0.3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side-peak width values measured in pp collisions at $\sqrt{s} = 7$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side-peak width values measured in pp collisions at $\sqrt{s} = 7$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $p_{\rm T} > 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Baseline values measured in pp collisions at $\sqrt{s} = 7$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $p_{\rm T} > 0.3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Baseline values measured in pp collisions at $\sqrt{s} = 7$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Baseline values measured in pp collisions at $\sqrt{s} = 7$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $p_{\rm T} > 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side-peak associated yield values measured in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $p_{\rm T} > 0.3$ GeV/$c$. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side-peak associated yield values measured in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side-peak associated yield values measured in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $p_{\rm T} > 1$ GeV/$c$. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side-peak width values measured in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $p_{\rm T} > 0.3$ GeV/$c$. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side-peak width values measured in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side-peak width values measured in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, for correlations of D mesons (D$^{0}$, D$^{+}$, D$^{*+}$ average) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated track $p_{\rm T} > 1$ GeV/$c$. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta| =|\eta_{\rm ch}-\eta_{\rm D}|< 1$.