Fig. 6b: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Away and Toward regions after the subtraction of Number density $N_{\rm ch}$ and $\Sigma p_{\rm T}$ distributions in the transverse region for p-Pb collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. 7a: $<p_{\rm T}>$ as a function of $p_{\rm T}^{\rm trig}$ in Toward (left) and Away (right) regions after the subtraction of transverse region for pp collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. 7b: $<p_{\rm T}>$ as a function of $p_{\rm T}^{\rm trig}$ in Toward (left) and Away (right) regions after the subtraction of transverse region for p-Pb collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A1: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 0.15 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A2: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 1.0 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A3: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 0.15 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A4: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 1.0 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Mean $p_T$ values of the $\Upsilon(1$S$)$, $\Upsilon(2$S$)$, and $\Upsilon(3$S$)$ states with $p_T\,>\,0\,GeV$ and $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $0<p_T<5\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $5<p_T<7\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $7<p_T<9\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $9<p_T<11\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $11<p_T<15\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $15<p_T<20\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $20<p_T<50\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $0<p_T<5\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $5<p_T<7\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $7<p_T<9\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $9<p_T<11\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $11<p_T<15\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $15<p_T<20\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $20<p_T<50\,$GeV range

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of "forward" track multiplicity $N_{track}^{\Delta\phi}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$. Forward tracks are those with momentum direction in $\Delta\phi\,<\,\pi/3$ w.r.t. the $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of "transverse" track multiplicity $N_{track}^{\Delta\phi}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$. Transverse tracks are those with momentum direction in $\pi/3\,<\,\Delta\phi\,<\,2\pi/3$ w.r.t. the $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of "backward" track multiplicity $N_{track}^{\Delta\phi}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$. Backward tracks are those with momentum direction in $\Delta\phi\,>\,2\pi/3$ w.r.t. the $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,=\,0$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,=\,1$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,=\,2$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,>\,2$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.00\,<\,S_T\,<\,0.55$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.55\,<\,S_T\,<\,0.70$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.70\,<\,S_T\,<\,0.85$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.85\,<\,S_T\,<\,1.00$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Efficiency-corrected multiplicity bins used in the $\Upsilon(n$S$)$ ratio analysis and the corresponding mean number of charged particle tracks.

Observed event yields and estimated backgrounds in the low- and high-mass control regions. The uncertainty in the background yield is the sum in quadrature of the statistical and systematic uncertainties.

Summary of the relative systematic uncertainties in heavy Majorana neutrino signal yields and the background from prompt same-sign leptons, both estimated from simulation. The relative systematic uncertainties assigned to the data-driven backgrounds for the misidentified lepton background and mismeasured charge background are also shown. The uncertainties are given for the low-mass selections.

Summary of the relative systematic uncertainties in heavy Majorana neutrino signal yields and the background from prompt same-sign leptons, both estimated from simulation. The relative systematic uncertainties assigned to the data-driven backgrounds for the misidentified lepton background and mismeasured charge background are also shown. The uncertainties are given for the high-mass selections.

Summary of contributions to the systematic uncertainty related to the prompt same-sign leptons background, misidentified lepton background, and mismeasured charge background on the total background uncertainty for the case of $m_{\rm N}$ = 100 and 500 GeV.

Observed event yields and estimated backgrounds after the application of all selection, except for the final optimization. The background predictions from prompt same-sign leptons (Prompt bkgd.), misidentified leptons (Misid. bkgd.), mismeasured charge (Charge mismeas. bkgd.) and the total background (Total bkgd.) are shown together withthe number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

Dielectron channel results after final optimization. The background predictions from prompt same-sign leptons, misidentified leptons, and mismeasured charge are shown along with the total background estimate and the number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

Electron-muon channel results after final optimization. The background predictions from prompt same-sign leptons and misidentified leptons are shown along with the total background estimate and the number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

Nuclear modicifation factor vs. N_coll $J\psi$ with|y|<1 and pT < 3 GeV/c in d Au collisions.

The invariant cross section versus transverse momentum for J/$\psi$ with $|y|\lt 1$ in $p+p$ collisions.

The invariant yield versus transverse momentum for $J/\Psi$ with |y| < 1 in 0 - 100 percent central d+Au collisions.

Fully corrected average charged particle scalar Sum(pT) per unit of pseudorapidity and per radian as a function of the leading track-jet transverse momentum for proton-proton collisions at a centre-of-mass energy of 2.76 TeV in the TransMAX region.

Fully corrected average charged particle multiplicity per unit of pseudorapidity and per radian as a function of the leading track-jet transverse momentum for proton-proton collisions at a centre-of-mass energy of 2.76 TeV in the TransMIN region.

Fully corrected average charged particle scalar Sum(pT) per unit of pseudorapidity and per radian as a function of the leading track-jet transverse momentum for proton-proton collisions at a centre-of-mass energy of 2.76 TeV in the TransMIN region.

Fully corrected average charged particle multiplicity per unit of pseudorapidity and per radian as a function of the leading track-jet transverse momentum for proton-proton collisions at a centre-of-mass energy of 2.76 TeV in the TransDIF region.

Fully corrected average charged particle scalar Sum(pT) per unit of pseudorapidity and per radian as a function of the leading track-jet transverse momentum for proton-proton collisions at a centre-of-mass energy of 2.76 TeV in the TransDIF region.

The dihadron correlated yield normalized per radian per unit of pseudorapidity as function of $\Delta\eta$ in d+Au collisions on the away side (|$\Delta\phi$ - $\pi$| < $\pi$/3). Shown is the high FTPC-Au activity data. Trigger and associated particles have 1 < $p_T$ < 3 GeV/c and |$\eta$| < 1.

The dihadron correlated yield normalized per radian per unit of pseudorapidity as function of $\Delta\eta$ in d+Au collisions on the away side (|$\Delta\phi$ - $\pi$| < $\pi$/3). Shown is the high-activity data after subtracting the unscaled low-activity data. Trigger and associated particles have 1 < $p_T$ < 3 GeV/c and |$\eta$| < 1.

The dihadron correlated yield normalized per radian per unit of pseudorapidity as function of $\Delta\eta$ in d+Au collisions on the away side (|$\Delta\phi$ - $\pi$| < $\pi$/3). Shown is the high-activity data after subtracting the unscaled low-activity data. Trigger and associated particles have 1 < $p_T$ < 3 GeV/c and |$\eta$| < 1.

The dihadron correlated yield normalized per radian per unit of pseudorapidity as function of $\Delta\eta$ in d+Au collisions on the near (|$\Delta\phi$| < $\pi$/3) side. Shown is the high-activity data after subtracting the unscaled low-activity data. Trigger and associated particles have 1 < $p_T$ < 3 GeV/c and |$\eta$| < 1.

The dihadron correlated yield normalized per radian per unit of pseudorapidity as function of $\Delta\eta$ in d+Au collisions on the near (|$\Delta\phi$| < $\pi$/3, solid circles) and away side (|$\Delta\phi$ - $\pi$| < $\pi$/3, open circles). Shown are the (a) low and (b) high FTPC-Au activity data, and the high-activity data after subtracting the (c) unscaled. Trigger and associated particles have 1 < $p_T$ < 3 GeV/c and |$\eta$| < 1.

(a) The near-side jetlike correlated yield obtained from Gaussian fit as in Fig. 1 as function of the uncorrected dN/d$\eta$ at midrapidity measured in the TPC. Two event selections are used: FTPC-Au multiplicity and ZDC-Au energy. The curve is the result from a HIJING calculation.

(a) The near-side jetlike correlated yield obtained from Gaussian fit as in Fig. 1 as function of the uncorrected dN/d$\eta$ at midrapidity measured in the TPC. Two event selections are used: FTPC-Au multiplicity and ZDC-Au energy. The curve is the result from a HIJING calculation.

$p_T^{(a)}$ dependence. The ratio of the correlated yields in high over low FTPC-Au multiplicity events as a function of $p_T^{(a)}$ $p_T^{(t)}$ where $p_T^{(a)}$ $p_T^{(t)}$ is fixed.

$p_T^{(a)}$ dependence. The ratio of the correlated yields in high over low FTPC-Au multiplicity events as a function of $p_T^{(a)}$ $p_T^{(t)}$ where $p_T^{(a)}$ $p_T^{(t)}$ is fixed.

(a) The dihadron correlated yield normalized per radian per unit of pseudorapidity as a function of $\Delta\phi$ in d+Au collisions at high (0-20%) FTPC-Au multiplicities. Trigger and associated particles are 1 < $p_T$ < 3 GeV/c within 0.5 < |$\Delta\eta$| < 0.7. ZYAM positions are indicated with arrows.

(a) The dihadron correlated yield normalized per radian per unit of pseudorapidity as a function of $\Delta\phi$ in d+Au collisions at low (40-100%) FTPC-Au multiplicities. Trigger and associated particles are 1 < $p_T$ < 3 GeV/c within 0.5 < |$\Delta\eta$| < 0.7. ZYAM positions are indicated with arrows.multiplicity versus $\Delta\phi$.

The raw associated yield at high FTPC-Au multiplicity minus the unscaled ZYAM-subtracted correlated yields at low FTPC-Au multiplicity versus $\Delta\phi$.

The raw associated yield at high FTPC-Au multiplicity minus the scaled ZYAM-subtracted correlated yields at low FTPC-Au multiplicity versus $\Delta\phi$.

Inclusive jet $A_{LL}$ vs. parton jet $p_T$ for 0.5<|eta|<1.0.

The uncertainties on the parton jet pT values are strongly correlated across the bins. The point-to-point statistical and systematic uncertainties on the measured A_LL values given above are only weakly correlated among the bins. The correlation matrix for the quadrature sum of the point-to-point statistical and systematic uncertainties is given below. In that matrix, the order of the rows and columns matches the bin numbering above. In addition to the point-to-point statistical and systematic uncertainties, there are two systematic uncertainties that are common to all the measured A_LL values (and not included in the systematic uncertainty numbers quoted above): (a) +/-0.0005 vertical shift uncertainty from relative luminosity and (b) +/-6.5% vertical scale uncertainty from polarization.

$A_{LL}$ model predictions for |eta|<0.5.

$A_{LL}$ model predictions for 0.5<|eta|<1.0.

Mean $p_T$, leading in-jet charged particle.

Mean $p_T$, charged particle jets, $p^{ch.jet}_T > 5$ GeV, $|\eta^{ch.jet}| < 1.9$.

Charged jet rate, $p^\text{ch.jet}_T > 5$ GeV, $|\eta^{ch.jet}| < 1.9$.

Charged jet rate, $p^\text{ch.jet}_T > 30$ GeV, $|\eta^{ch.jet}| < 1.9$.

Jet $p_T$ spectrum, $|\eta^{ch.jet}| < 1.9$, $10 < N_\text{ch} \le 30$.

Jet $p_T$ spectrum, $|\eta^{ch.jet}| < 1.9$, $30 < N_\text{ch} \le 50$.

Jet $p_T$ spectrum, $|\eta^{ch.jet}| < 1.9$, $50 < N_\text{ch} \le 80$.

Jet $p_T$ spectrum, $|\eta^{ch.jet}| < 1.9$, $80 < N_\text{ch} \le 110$.

Jet $p_T$ spectrum, $|\eta^{ch.jet}| < 1.9$, $110 < N_\text{ch} \le 140$.

Intrajet ring $p_{T}$ density, $10 < N_\text{ch} \le 30$.

Intrajet ring $p_{T}$ density, $30 < N_\text{ch} \le 50$.

Intrajet ring $p_{T}$ density, $50 < N_\text{ch} \le 80$.

Intrajet ring $p_{T}$ density, $80 < N_\text{ch} \le 110$.

Intrajet ring $p_{T}$ density, $110 < N_\text{ch} \le 140$.

(Color online.) $J/\psi$ RAA versus pT for several centrality bins for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The statistical (systematic) uncertainties are shown with vertical bars (open boxes). The filled boxes about unity on the right show the size of the normalization uncertainty. PHENIX low-$p_{T}$ $J/\psi$ results [6] and STAR high-$p_{T}$ results in Cu+Cu collisions [15] are shown for comparison. The curves are from the predictions by Model I (Liu et al.) [11] and Model II (Zhao et al.) [12].

(Color online) $R_{AA}$ versus $N_{part}$ for high-$p_{T}$ $J/\psi$ in comparison with a) low-$p_{T}$ $J/\psi$ data from PHENIX [6], and high-$p_{T}$ $\pi^{\pm}$ from STAR [41, 42], and b) high-$p_{T}$ $J/\psi$ from CMS [8]. The statistical (systematic) uncertainties are shown in vertical bars (boxes). The shaded green band about unity shows the systematic uncertainties from $N_{bin}$ and the box about unity on the right shows the $R_{AA}$ normalization uncertainty from the statistical and global systematic uncertainties of the $p+p$ reference data. The percentage in the top horizontal axis of (a) and (b) refer to the collision centrality for high-$p_{T}$ $J/\psi$ data at the RHIC and LHC, respectively.

Toward $\sum p_\perp$ density vs $p_\perp^{\mu\mu}$.

Transverse $\sum p_\perp$ density vs $p_\perp^{\mu\mu}$.

Away $\sum p_\perp$ density vs $p_\perp^{\mu\mu}$.

Toward $\langle p_\perp \rangle$ vs $p_\perp^{\mu\mu}$.

Transverse $\langle p_\perp \rangle$ vs $p_\perp^{\mu\mu}$.

Away $\langle p_\perp \rangle$ vs $p_\perp^{\mu\mu}$.

Towards $+$ transverse $N_\text{chg}$ density vs $m_{\mu\mu}$, $p_\perp^{\mu\mu} < 5$ GeV.

Towards $+$ transverse $\sum p_\perp$ density vs $m_{\mu\mu}$, $p_\perp^{\mu\mu} < 5$ GeV.

Towards $+$ transverse $\langle p_\perp \rangle$ vs $m_{\mu\mu}$, $p_\perp^{\mu\mu} < 5$ GeV.

Toward $N_\text{chg}$.

Transverse $N_\text{chg}$.

Away $N_\text{chg}$.

Toward $p_\perp$.

Transverse $p_\perp$.

Away $p_\perp$.

Transverse $N_\text{chg}$, $p_\perp(\mu\mu) < 5$ GeV.

Transverse $p_\perp$, $p_\perp(\mu\mu) < 5$ GeV.