Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH3 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 and ISR starting scale SpaceShower:pT0Ref = 1 GeV at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH3 and ISR starting scale SudakovCommon:pTmin=0.7 GeV at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-proton center-of-mass energy 38.8 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-proton center-of-mass energy 8000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-nucleon center-of-mass energy 8160 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-nucleon center-of-mass energy 13000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 38.8 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 8000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-nucleon center-of-mass energy 8160 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 13000 GeV versus the dilepton mass range of the process.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Self-Normalized $J/\psi$ yields in the same arms before and after dimuon subtraction and the comparison with PYTHIA 8 simulations

Self-normalized $\psi(2S)$ to $J/\psi$ yields ratio as a function of Self-normalized $N_{ch}$ in the same arms

Self-normalized $\psi(2S)$ to $J/\psi$ yields ratio as a function of Self-normalized $N_{ch}$ in opposite arms

Distributions of $m_T$ in the $W^{+}$ signal selection for mu final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_T$ in the $W^{-}$ signal selection for e final states for the pp collisions at $\sqrt{s}=$ 5TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_T$ in the $W^{-}$ signal selection for mu final states for the pp collisions at $\sqrt{s}=$ 5TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_T$ in the $W^{-}$ signal selection for e final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_T$ in the $W^{-}$ signal selection for mu final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_{ll}$ in the Z signal selection for ee final states for the pp collisions at $\sqrt{s}=$ 5TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to l\nu$, and diboson processes.

Distributions of $m_{ll}$ in the Z signal selection for mumu final states for the pp collisions at $\sqrt{s}=$ 5TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to l\nu$, and diboson processes.

Distributions of $m_{ll}$ in the Z signal selection for ee final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to l\nu$, and diboson processes.

Distributions of $m_{ll}$ in the Z signal selection for mumu final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to l\nu$, and diboson processes.

Figure 4 Fiducial cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes at 13TeV

Figure 5 Fiducial cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes at 5.02TeV

Figure 6 Fiducial cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes between 13TeV and 5.02TeV

Figure 7 Total cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes at 13TeV

Figure 8 Total cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes at 5.02TeV

Figure 9 Total cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes between 13TeV and 5.02TeV

Summary of theoretical predictions of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at ppbar collisions.

Summary of experimental measurements of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at UA1 and UA2.

Summary of experimental measurements of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at D0.

Summary of experimental measurements of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at CDF.

Summary of theoretical predictions of the total inclusive cross sections for W boson production times leptonic branching fractions at pp collisions.

Summary of experimental measurements of the total inclusive cross sections for Z boson production times leptonic branching fractions at pp collisions.

Summary of experimental measurements of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at CMS.

Values of $N^{EXP}_{coll}$ versus $p_{T}$ as defined in Eq. 2 for three d+Au event classes, (a) $0\%–100\%$, (b) $0\%–5\%$ and (c) $60\%–88\%$. Also shown are fits to the data (solid lines) and the corresponding value $N^{GL}_{coll}$ (dashed lines).

Panels (d) to (f) show the nuclear modification factors $R^{\pi^{0}}_{dAu,EXP}$, calculated with Eq. 3, for the same event selections as in (a) to (c) together with fits to the data, (d) $0\%–100\%$, (e) $0\%–5\%$ and (f) $60\%–88\%$.

The Ratio $N^{EXP}_{coll}/N^{GL}_{coll}$ as a function of $N^{EXP}_{coll}$. Horizontal and vertical bars are the statistical uncertainties. The values for $0\%–100\%$ d+Au collisions are represented by a solid blue line, with the statistical uncertainty given as band. The scale uncertainties that are common to all data points are shown for the $0\%–100\%$ value.

Average $R^{\pi^{0}}_{dAu,EXP}$ as a function of $N^{EXP}_{coll}$. Horizontal and vertical bars are the statistical uncertainties. The values for $0\%–100\%$ d+Au collisions are represented by a solid blue line, with the statistical uncertainty given as band. The scale uncertainties(16.5%) that are common to all data points are shown for the $0\%–100\%$ value.

Average $R^{\pi^{0}}_{dAu,EXP}$ as a function of $N^{EXP}_{coll}$. Horizontal and vertical bars are the statistical uncertainties. The values for $0\%–100\%$ d+Au collisions are represented by a solid blue line, with the statistical uncertainty given as band. The scale uncertainties($16.5\%$) that are common to all data points are shown for the $0\%–100\%$ value. The systematic uncertainties on $N^{EXP}_{coll}$, $\approx16\%$, are dominated by uncertainties on the p+p data set and thus are a common scale uncertainty for all d+Au event classes.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in pp collisions at $\sqrt{s} = 5020~\mathrm{GeV}$.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in Pb-p collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in pp collisions at various LHC energies.

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of collision energy in pp collisions.

Centrality-dependent pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$ (Centrality estimator - CL1).

Centrality-dependent pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$ (Centrality estimator - ZNA).

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of $\langle N_{\rm part} \rangle$ for CL1 centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (high-$p_{\rm T}$) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (Pb-side) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm part} \rangle$ for CL1 centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (high-$p_{\rm T}$) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (Pb-side) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm \gamma} \rangle$ for CL1 centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm \gamma} \rangle$ for ZNA (high-$p_{\rm T}$) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm \gamma} \rangle$ for ZNA (Pb-side) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

The $\phi$-meson integrated nuclear modification factors $R_{AB}$ measured as a function of $N_{part}$

The comparison of $\phi$ meson (a) $v_2$ and (b) $v_2/(\epsilon N_{part}^{1/3})$ measured as a function of $p_T$ in Cu+Au collisions at $\sqrt{s}$ = 200 GeV and U+U collisions at $\sqrt{s}$ = 193 GeVat midrapidity (|η| < 0.35)

The comparison of elliptic flow (a) $v_2$ values measured for $\phi$ mesons as a function of $p_T$ in $0\%–20\%$ Cu+Au collisions

The comparison of elliptic flow (b) $v_2$ values measured for $\phi$ mesons as a function of $KE_T$ in $0\%–20\%$ Cu+Au collisions

The comparison of elliptic flow (c) $v_2/n_q$ values measured for $\phi$ mesons as a function of $KE_T/n_q$ in $0\%–20\%$ Cu+Au collisions

The comparison of elliptic flow (a) $v_2$ values measured for $\phi$ mesons as a function of $p_T$ in $20\%–60\%$ Cu+Au collisions

The comparison of elliptic flow (b) $v_2$ values measured for $\phi$ mesons as a function of $KE_T$ in $20\%–60\%$ Cu+Au collisions

The comparison of elliptic flow (c) $v_2/n_q$ values measured for $\phi$ mesons as a function of $KE_T/n_q$ in $20\%–60\%$ Cu+Au collisions

The comparison of elliptic flow (a) $v_2$ values measured for $\phi$ mesons as a function of $p_T$ in $0\%–50\%$ U+U collisions

The comparison of elliptic flow (b) $v_2$ values measured for $\phi$ mesons as a function of $KE_T$ in $0\%–50\%$ U+U collisions

The comparison of elliptic flow (c) $v_2/n_q$ values measured for $\phi$ mesons as a function of $KE_T/n_q$ in $0\%–50\%$ U+U collisions

Invariant yield of direct photons, every 20% centrality

Invariant yield of nonprompt photons, every 20% centrality

$T_{eff}$ vs charged particle multiplicity

Integrated direct photon yield of 1-5GeV vs charged particle multiplicity

Integrated nonprompt direct photon yield of different $p_T$ integration window vs charged particle multiplicity

Scaling factor $\alpha$ vs charged particle multiplicity

Direct photon spectra for minimum bias (0-86%) Au+Au collision at $\sqrt{s_{NN}} = 62.4$ GeV (a) and $39$ GeV (b). For $62.4$ GeV also centrality bins of 0-20% (c) and 20-40% (d) are shown. Data points are shown with statistical and systematic uncertainties, unless the central value is negative (arrows) or is consistent with zero within the statistical uncertainties (arrows with data point). In these cases upper limit with CL = 95$%$ are given.

Inverse slopes, $T_{eff}$ , obtained from fitting the combined data from central collisions shown in Fig. 11 is compared to the fit results of the individual data sets at $\sqrt{s_{NN}} = 62.4$ GeV, $200$ GeV, and $2760$ GeV. Also included is the value for $\sqrt{s_{NN}} = 39$ GeV obtained from fitting the min. bias data set in the lower $p_T$ range.

Integrated invariant direct-photon yields vs. charged particle multiplicity for $p_{T}$ integrated from (a) 0.4 GeV/c, (b) 1.0 GeV/c, (c) 1.5 GeV/c, and (d) 2.0 to 5.0 GeV/c for all available A+A datasets.

Integrated invariant direct-photon yields vs. charged particle multiplicity for $p_{T}$ integrated from (a) 0.4 GeV/c, (b) 1.0 GeV/c, (c) 1.5 GeV/c, and (d) 2.0 to 5.0 GeV/c for all available A+A datasets.

Integrated invariant direct-photon yields vs. charged particle multiplicity for $p_{T}$ integrated from (a) 0.4 GeV/c, (b) 1.0 GeV/c, (c) 1.5 GeV/c, and (d) 2.0 to 5.0 GeV/c for all available A+A datasets.

Integrated invariant direct-photon yields vs. charged particle multiplicity for $p_{T}$ integrated from (a) 0.4 GeV/c, (b) 1.0 GeV/c, (c) 1.5 GeV/c, and (d) 2.0 to 5.0 GeV/c for all available A+A datasets.

Integrated direct-photon yields from A+A collisions for $p_{T,min}$ of $5$ GeV/c (a) and $8$ GeV/c. The representation is the same as in Fig. 13. Also shown are the results from pQCD calculations scaled by $N_{coll}$

Integrated direct-photon yields from A+A collisions for $p_{T,min}$ of $5$ GeV/c (a) and $8$ GeV/c. The representation is the same as in Fig. 13. Also shown are the results from pQCD calculations scaled by $N_{coll}$

The $\alpha$ values extracted using fits to integrated direct photon yields. The dashed line gives the average $\alpha$ value for the 4 lower $p_{T,min}$ points. Also shown is a model calculation for $\alpha$ discussed in the text.