Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 4. (Color online.) Top panel: The number of photons divided by $\left\langle N_{\text{part}}\right\rangle / 2$ as a function of average number of participating nucleons for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s} _{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$ for $-3.7<\eta<-2.3 .$ Errors shown are systematic only and include those for $\left\langle N_{\text {part }}\right\rangle .$ Results from HIJING are shown as lines (solid for $\mathrm{Au}+\mathrm{Au}$ and dashed for $\mathrm{Cu}+\mathrm{Cu}$ ). Bottom panel: Same as above, for both photons and charged particles from PHOBOS [10]. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 4. (Color online.) Top panel: The number of photons divided by $\left\langle N_{\text{part}}\right\rangle / 2$ as a function of average number of participating nucleons for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s} _{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$ for $-3.7<\eta<-2.3 .$ Errors shown are systematic only and include those for $\left\langle N_{\text {part }}\right\rangle .$ Results from HIJING are shown as lines (solid for $\mathrm{Au}+\mathrm{Au}$ and dashed for $\mathrm{Cu}+\mathrm{Cu}$ ). Bottom panel: Same as above, for both photons and charged particles from PHOBOS [10]. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 5. (Color online.) Photon pseudorapidity distributions normalized by the average number of participating nucleon pairs for different collision centralities are plotted as a function of pseudorapidity shifted by the beam rapidity (-5.36 for $200 \mathrm{GeV}$ and -4.19 for $62.4 \mathrm{GeV}$ ) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ collisions at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$. Errors are systematic only, statistical errors are negligible in comparison. For clarity of presentation, results for only four centralities are shown. The $\mathrm{Cu}+$ Cu data are shifted by 0.1 unit in $\eta-y_{\text {beam }} .$ The solid line is a second order polynomial fit to the data (see text for details).

Fig. 5. (Color online.) Photon pseudorapidity distributions normalized by the average number of participating nucleon pairs for different collision centralities are plotted as a function of pseudorapidity shifted by the beam rapidity (-5.36 for $200 \mathrm{GeV}$ and -4.19 for $62.4 \mathrm{GeV}$ ) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ collisions at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$. Errors are systematic only, statistical errors are negligible in comparison. For clarity of presentation, results for only four centralities are shown. The $\mathrm{Cu}+$ Cu data are shifted by 0.1 unit in $\eta-y_{\text {beam }} .$ The solid line is a second order polynomial fit to the data (see text for details).

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Photon conversion point radius distribution from real data and MC simulation in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Conversion probability correction factor for $\pi^0$ as a function of $p_{T}$

The overall detection efficiency of $\pi^0$ from embedding studt in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The overall detection efficiency of $\pi^0$ from embedding studt in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The overall detection efficiency of $\pi^0$ from embedding studt in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Invariant yield of STAR $\pi^0$ as a function of $p_{T}$ at midrapidity for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Invariant yield of STAR $\pi^0$ as a function of $p_{T}$ at midrapidity for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Invariant yield of STAR $\pi^0$ as a function of $p_{T}$ at midrapidity for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The ratios of STAR $\pi^0$ spectra over $\pi^{\pm}$ from STAR and $\pi^0$ from PHENIX for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The nuclear modification factor $R_{CP}$ as a functon of $p_{T}$ of STAR $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The nuclear modification factor $R_{AA}$ as a functon of $p_{T}$ of STAR $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Sample correlations in $\Delta\phi$ ($|\Delta\eta|<$ 1.78) for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-95% $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV, 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV, and 0-12% central Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. The data are averaged between positive and negative $\Delta\phi$ and reflected in the plot.

Background subtracted sample correlations for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ on the near-side for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-95% $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV, 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV and 0-12% central Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. The dependence of the jet-like correlation is shown as a function of $\Delta\eta$ ($|\Delta\phi|<$ 0.78). The data are averaged between positive and negative $\Delta\eta$.

Background subtracted sample correlations for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ on the near-side for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-95% $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV, 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV and 0-12% central Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. The dependence of the jet-like correlation is shown as a function of $\Delta\phi$ ($|\Delta\eta|<$ 1.78). The data are averaged between positive and negative $\Delta\phi$.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $p_T^{\mathrm{trigger}}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-95% $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV, 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV and 0-12% central Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. The 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $p_T^{\mathrm{associated}}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ for 0-60% Cu+Cu and 0-80% Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV and 0-95% $d$+Au, 0-60% Cu+Cu, 0-12% Au+Au and 40-80% Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of the widths in $\Delta\phi$ on $p_T^{\mathrm{trigger}}$ for 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-95% $d$+Au, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV and $\sqrt{s_{NN}}$ = 200 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, and 0-12% and 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\phi$ on $p_T^{\mathrm{associated}}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ for 0-95% $d$+Au, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV and $\sqrt{s_{NN}}$ = 200 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, and 0-12% and 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-95% $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $p_T^{\mathrm{trigger}}$ for 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-95% $d$+Au, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV and $\sqrt{s_{NN}}$ = 200 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, and 0-12% and 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $p_T^{\mathrm{associated}}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ for 0-95% $d$+Au, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV and $\sqrt{s_{NN}}$ = 200 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, and 0-12% and 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-95% $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of ridge yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of ridge yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of ridge yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of ridge yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Ratio of the ridge and jet-like yields on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Ratio of the ridge and jet-like yields on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Ratio of the ridge and jet-like yields on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Ratio of the ridge and jet-like yields on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Transverse momentum distribution for the second leading lepton for events in the SR2 signal region for DATA and SM predictions.

Transverse momentum distribution for the third leading lepton for events in the SR1 signal region for DATA and SM predictions.

Transverse momentum distribution for the third leading lepton for events in the SR2 signal region for DATA and SM predictions.

Missing transverse energy for events in the SR1 signal region for DATA and SM predictions.

Missing transverse energy for events in the SR2 signal region for DATA and SM predictions.

Invariant mass of the same-flavour-opposite-sign (SFOS) lepton pair for events in the SR1 signal region for DATA and SM predictions.

The Cross Section in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Cross Section in signal region SR1 for the SUSY simplified model grid.

The Number of generated Events in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Number of generated Events in signal region SR1 for the SUSY simplified model grid.

The Efficiency in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Efficiency in signal region SR1 for the SUSY simplified model grid.

The Acceptance in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Acceptance in signal region SR1 for the SUSY simplified model grid.

The Acceptance*Efficiency in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Acceptance*Efficiency in signal region SR1 for the SUSY simplified model grid.

The Systematic Uncertainty of the data (excluding the Monte Carlo) in signal region SR1 for the SUSY pMSSM model with M1=100 GeV grid.

The Systematic Uncertainty of the data (excluding the Monte Carlo) in signal region SR1 for the SUSY simplified model grid.

CL values for the pMSSM with M1=100 GeV model grid for the SR1 signal region.

CL values for the simplified model model grid for the SR1 signal region.

Numbers of observed and background events for SR3L low for each bin of the distribution in Meff. The table corresponds to Fig. 4(d). The statistical and systematic uncertainties are combined for the predicted numbers.

Numbers of observed and background events for SR3L high for each bin of the distribution in Meff. The table corresponds to Fig. 4(e). The statistical and systematic uncertainties are combined for the predicted numbers.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The efficiencies are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The efficiencies are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

The efficiencies are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into b s and gluinos decay into t stop (see Fig. 5d in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The efficiencies are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The efficiencies are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The acceptances (in percent, %) are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The acceptances (in percent, %) are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

The acceptances (in percent, %) are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The acceptances (in percent, %) are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The acceptances (in percent, %) are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The limits on observed cross section are calculated for all simplified models. The simplified models are for direct pair production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct pair-production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct pair production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The limits on observed cross sections are calculated for all simplified models. The simplified models are for direct pair production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The signal event yields are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The signal event yields are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop)-20 GeV.

The signal event yields are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The signal event yields are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

Experimental uncertainties on the signal event yields are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

Experimental uncertainties on the signal event yields are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

Experimental uncertainties on the signal event yields are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

Experimental uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the values are given for the five signal regions and their combination.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

Statistical uncertainties on the signal event yields are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, and mu>0.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

Statistical uncertainties on the signal event yields are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

Statistical uncertainties on the signal event yields are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the values are given for the five signal regions and their combination. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

Statistical uncertainties on the signal event yields are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W ^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluino), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The confidence levels are calculated for all simplified models. For each model, the observed and expected values are given. The simplified model is for direct production of gluinos that decay into t tbar t tbar chi1^0 chi1^0 (see Fig. 5a in the paper). This particular model assumes that top quark is much heavier than gluino.

The confidence levels are calculated for all simplified models. For each model, the observed and expected values are given. The simplified model is for direct production of squarks that decay into two steps into q q W Z W Z chi1^0 chi1^0 (see Fig. 6c in the paper).

The confidence levels are calculated for all simplified models. For each model, the values are given for the five signal regions and their combination. The simplified model is for direct pair-production of gluinos that decay via a two-step process into q q q q W Z W Z chi1^0 chi1^0 (see Fig. 6b in the paper).

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of gluinos that decay via sleptons into q q q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6d in the paper).

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct pair-production of gluinos. A gluino decays into t stop. Consequently, a top squark squark decays into b chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 5b in the paper). This particular model assumes that m(stop) < m(gluion), m(chi1^0)=6 GeV, and m(chi1^(+-))=118 GeV.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of gluinos. A gluino decays into t c chi1^0 (see Fig. 5c in the paper). This particular model assumes that m(chi1^0) = m(stop) - 20 GeV.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7b in the paper). This particular model assumes that m(chi1^0)=2(chi1^0).

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of bottom squarks. A bottom squark decays into t chi1^(+-) and chi1^(+-) --> W^(+-) chi1^0 (see Fig. 7a in the paper). This particular model assumes that m(chi1^0)=60 GeV.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of squarks. Squarks decay into q q l l (l l) chi1^0 chi1^0 + neutrinos (see Fig. 6e in the paper).

The confidence levels are calculated for all GMSB models (see Fig. 8c in the paper). For each model, the expected and observed values are given. The model assumes mmess=250 TeV, m5=3, mu>0, and Cgrav=1.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of gluinos and top squarks. Top squarks undergo R-parity violating decays into bs and gluinos decay into t stop (see Fig. 5d in the paper).

The confidence levels are calculated for all mSUGRA/CMSSM models with bRPV (see Fig. 8b in the paper). For each model, the expected and observed values are given. The model assumes tan(beta)=30, A0=2m0, mu>0, and bRPV.

The confidence levels are calculated for all simplified extra dimension model (see Fig. 8d in the paper). For each model, the expected and observed values are given.

The confidence levels are calculated for all simplified models. For each model, the expected and observed values are given. The simplified model is for direct production of gluinos that decay into q q q q W W chi1^0 chi1^0 (see Fig. 6a in the paper).

The confidence levels are calculated for all mSUGRA models (see Fig. 8a in the paper). For each model, the expected and observed values are given. The model assumes tan(beta)=30, A0=2m0, and mu>0.

Observed and expected numbers of events with exactly three leptons with the scalar sum of jet transverse momentum values HT < 200 GeV. "On-Z" refers to events with at least one E+ E- or MU+ MU- (OSSF) pair with dilepton mass between 75 and 105 GeV, while "Off-Z" refers to events with one or two OSSF pairs, none of which fall in this mass range. The OSSFN designation refers to the number of E+ E- and MU+ MU- pairs in the event. Search channels binned in ET have been combined into coarse ET bins for the purposes of presentation.

The product of acceptance and efficiency for the calculation of the 95% C.L. limits for the natural Higgsino NLSP scenario with BR(NEUTRALINO1 --> HIGGS GRAVITINO) = BR(NEUTRALINO1 --> Z GRAVITINO) = 0.5 for strong plus electroweak production.

The product of acceptance and efficiency for the calculation of the 95% C.L. limits for the natural Higgsino NLSP scenario with BR(NEUTRALINO1 --> HIGGS GRAVITINO) = 1.0 for strong plus electroweak production.

The product of acceptance and efficiency for the calculation of the 95% C.L. limits for the natural Higgsino NLSP scenario with BR(NEUTRALINO1 --> Z0 GRAVITINO) = 1.0 for strong plus electroweak production.

The product of acceptance and efficiency for the calculation of the 95% C.L. limit for the natural Higgsino NLSP scenario with strong plus electroweak production for BR(NEUTRALINO1 --> HIGGS GRAVITINO) = BR(NEUTRALINO1 --> Z0 GRAVITINO) = 0.5. The chargino-Higgsino mass is fixed at 150 GeV.

The product of acceptance and efficiency for the calculation of the 95% C.L. limit for the natural Higgsino NLSP scenario with strong plus electroweak production for BR(NEUTRALINO1 --> HIGGS GRAVITINO) = BR(NEUTRALINO1 --> Z0 GRAVITINO) = 0.5. The chargino-Higgsino mass is fixed at 350 GeV.

The product of acceptance and efficiency for the calculation of the 95% C.L. limits for the natural Higgsino NLSP scenario with strong plus electroweak production for BR(NEUTRALINO1 --> HIGGS GRAVITINO) = 1.0. The chargino-Higgsino mass is fixed at 150 GeV.

The product of acceptance and efficiency for the calculation of the 95% C.L. limits for the natural Higgsino NLSP scenario with strong plus electroweak production for BR(NEUTRALINO1 --> HIGGS GRAVITINO) = 1.0. The chargino-Higgsino mass is fixed at 350 GeV.

The product of acceptance and efficiency for the calculation of the 95% C.L. limits for the natural Higgsino NLSP scenario with strong plus electroweak production for BR(NEUTRALINO1 --> Z0 GRAVITINO) = 1.0. The chargino-Higgsino mass is fixed at 150 GeV.

The product of acceptance and efficiency for the calculation of the 95% C.L. limits for the natural Higgsino NLSP scenario with strong plus electroweak production for BR(NEUTRALINO1 --> Z0 GRAVITINO) = 1.0. The chargino-Higgsino mass is fixed at 350 GeV.

The product of acceptance and efficiency for the calculation of the 95% C.L. limit on BR(NEUTRALINO1 --> HIGGS GRAVITINO) for the natural Higgsino NLSP scenario with fixed chargino-Higgsino mass of 150 GeV assuming BR(NEUTRALINO1 --> HIGGS GRAVITINO) + BR(NEUTRALINO1 --> Z GRAVITINO) = 1.0 and strong plus weak production.

The product of acceptance and efficiency for the calculation of the 95% C.L. limit on BR(NEUTRALINO1 --> HIGGS GRAVITINO) for the natural Higgsino NLSP scenario with fixed chargino-Higgsino mass of 250 GeV assuming BR(NEUTRALINO1 --> HIGGS GRAVITINO) + BR(NEUTRALINO1 --> Z0 GRAVITINO) = 1.0 and strong plus weak production.

The product of acceptance and efficiency for the calculation of the 95% C.L. limit on BR(NEUTRALINO1 --> HIGGS GRAVITINO) for the natural Higgsino NLSP scenario with fixed chargino-Higgsino mass of 350 GeV assuming BR(NEUTRALINO1 --> HIGGS GRAVITINO) + BR(NEUTRALINO1 --> Z0 GRAVITINO) = 1.0 and strong plus weak production.

The product of acceptance and efficiency for the calculation of the 95% C.L. limits for the slepton co-NLSP model in the gluino-chargino (wino-like chargino) mass plane.

The product of acceptance and efficiency for the calculation of the 95% C.L. limits for the stau-NLSP and stau-NNLSP scenario in the degenerate smuon- and selectron- stau mass plane.

The product of acceptance and efficiency for the calculation of the 95% C.L. limits for the T1tttt model in the gluino-neutralino mass plane.

The product of acceptance and efficiency for the calculation of the 95% C.L. limits for the T6ttWW model in the sbottom-chargino mass plane.

Observed 95% C.L. exclusion contour for the T1tttt scenario.

Observed - 1sig(theoretical) 95% C.L. exclusion contour for the T1tttt scenario.

Observed + 1sig(theoretical) 95% C.L. exclusion contour for the T1tttt scenario.

Expected 95% C.L. exclusion contour for the T1tttt scenario.

Expected - 1sig 95% C.L. exclusion contour for the T1tttt scenario.

Expected + 1sig 95% C.L. exclusion contour for the T1tttt scenario.

Observed 95% C.L. exclusion contour for the T6ttWW scenario.

Observed - 1sig(theoretical) 95% C.L. exclusion contour for the T6ttWW scenario.

Observed + 1sig(theoretical) 95% C.L. exclusion contour for the T6ttWW scenario.

Expected 95% C.L. exclusion contour for the T6ttWW scenario.

Expected - 1sig 95% C.L. exclusion contour for the T6ttWW scenario.

Expected + 1sig 95% C.L. exclusion contour for the T6ttWW scenario.

Observed and expected 95% C.L. limits for the natural Higgsino NLSP scenario with BR(NEUTRALINO1 --> HIGGS GRAVITINO) = BR(NEUTRALINO1 --> Z GRAVITINO) = 0.5 for strong plus electroweak production.

Observed and expected 95% C.L. limits for the natural Higgsino NLSP scenario with BR(NEUTRALINO1 --> HIGGS GRAVITINO) = 1.0 for strong plus electroweak production. Successive instances of SIG(OBS) = SIG(EXP) =0 have been rebinned over M(STOP) to reduce the size of the table.

Observed and expected 95% C.L. limits for the natural Higgsino NLSP scenario with BR(NEUTRALINO1 --> Z0 GLUINO) = 1.0 for strong plus electroweak production.

Observed 95% C.L. limits on BR(NEUTRALINO1 --> HIGGS GRAVITINO) for the natural Higgsino NLSP scenario with fixed chargino-Higgsino mass of 150 GeV assuming BR(NEUTRALINO1 --> HIGGS GRAVITINO) + BR(NEUTRALINO1 --> Z GRAVITINO) = 1.0 and strong plus weak production.

Observed 95% C.L. limits on BR(NEUTRALINO1 --> HIGGS GRAVITINO) for the natural Higgsino NLSP scenario with fixed chargino-Higgsino mass of 250 GeV assuming BR(NEUTRALINO1 --> HIGGS GRAVITINO) + BR(NEUTRALINO1 --> Z GRAVITINO) = 1.0 and strong plus weak production.

Observed 95% C.L. limits on BR(NEUTRALINO1 --> HIGGS GRAVITINO) for the natural Higgsino NLSP scenario with fixed chargino-Higgsino mass of 350 GeV assuming BR(NEUTRALINO1 --> HIGGS GRAVITINO) + BR(NEUTRALINO1 --> Z GRAVITINO) = 1.0 and strong plus weak production.

Observed and expected 95% C.L. limits for the slepton co-NLSP model in the gluino-chargino (wino-like chargino) mass plane are shown.

Observed and expected 95% C.L. limits for the stau-NLSP and stau-NNLSP scenario in the degenerate smuon- and selectron- stau mass plane are shown. Successive instances of SIG(OBS) = SIG(EXP) =0 have been rebinned over M(SELECTRON) to reduce the size of the table.

Observed and expected 95% C.L. limits for the T1tttt model in the gluino-neutralino mass plane are shown.

Observed and expected 95% C.L. limits for the T6ttWW model in the sbottom-chargino mass plane are given.

The effective mass distribution in VR2Z.

The ETmiss distribution in SR0noZa.

The effective mass distribution in SR0noZa.

The ETmiss distribution in SR1noZa.

The effective mass distribution in SR1noZa.

The ETmiss distribution in SR2noZa.

The effective mass distribution in SR2noZa.

The ETmiss distribution in SR0noZb.

The effective mass distribution in SR0noZb.

The ETmiss distribution in SR1noZb.

The effective mass distribution in SR1noZb.

The ETmiss distribution in SR2noZb.

The effective mass distribution in SR2noZb.

The ETmiss distribution in SR0Z.

The effective mass distribution in SR0Z.

The ETmiss distribution in SR1Z.

The effective mass distribution in SR1Z.

The ETmiss distribution in SR2Z.

The effective mass distribution in SR2Z.

Observed 95% CL exclusion contour for the RPV chargino NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV chargino NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV chargino NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV chargino NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV chargino NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV chargino NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV chargino NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV chargino NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the RPV gluino NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV gluino NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV gluino NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV gluino NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV gluino NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV gluino NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV gluino NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV gluino NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV Lslepton NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV Rslepton NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_121 != 0.

Expected 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_121 != 0.

Observed 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_122 != 0.

Expected 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_122 != 0.

Observed 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_133 != 0.

Expected 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_133 != 0.

Observed 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_233 != 0.

Expected 95% CL exclusion contour for the RPV sneutrino NLSP model with lambda_233 != 0.

Observed 95% CL exclusion contour for the R-slepton RPC model.

Expected 95% CL exclusion contour for the R-slepton RPC model.

Observed and expected 95% CL cross-section upper limits for the Stau RPC model, together with the theoretically predicted cross-section.

Observed and expected 95% CL cross-section upper limits for the Z RPC model, together with the theoretically predicted cross-section.

Observed 95% CL exclusion contour for the GGM tan beta = 1.5 model.

Expected 95% CL exclusion contour for the GGM tan beta = 1.5 model.

Observed 95% CL exclusion contour for the GGM tan beta = 30 model.

Expected 95% CL exclusion contour for the GGM tan beta = 30 model.

Observed 95% CL cross-section upper limit for the RPV chargino NLSP models with lambda_121 != 0 and lambda_122 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV chargino NLSP models with lambda_133 != 0 and lambda_233 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV gluino NLSP models with lambda_121 != 0 and lambda_122 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV gluino NLSP models with lambda_133 != 0 and lambda_233 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV Lslepton NLSP models with lambda_121 != 0 and lambda_122 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV Lslepton NLSP models with lambda_133 != 0 and lambda_233 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV Rslepton NLSP models with lambda_121 != 0 and lambda_122 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV Rslepton NLSP models with lambda_133 != 0 and lambda_233 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV sneutrino NLSP models with lambda_121 != 0 and lambda_122 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the RPV sneutrino NLSP models with lambda_133 != 0 and lambda_233 != 0, and the selection of Z-veto signal regions used to set limits in these models. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bba' means that the regions SR0noZb, SR1noZb and SR2noZa were used, in addition to the three Z-rich regions (SR0-2Z).

Observed 95% CL cross-section upper limit for the R-slepton RPC model, and the selection of Z-veto signal regions used to set limits in this model. The combination of regions used is ordered by the minimum number of hadronic taus required. For example, ``bbb' means that the regions SR0noZb, SR1noZb and SR2noZb were used, in addition to the three Z-rich regions (SR0-2Z). For the RPC stau and Z models, the ``aaa' combination of regions was used throughout.

Performance of the SR0noZa selection in the R-slepton RPC model: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Performance of the SR0noZb selection in the RPV chargino NLSP model with lambda_121 != 0: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Performance of the SR1noZa selection in the RPV sneutrino NLSP model with lambda_233 != 0: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Performance of the SR1noZb selection in the RPV gluino NLSP model with lambda_133 != 0: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Performance of the SR2noZa selection in the RPV sneutrino NLSP model with lambda_233 != 0: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Performance of the SR2noZb selection in the RPV gluino NLSP model with lambda_133 != 0: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Performance of the SR0Z selection in the GGM tan beta = 30 model: number of generated signal events; total signal cross-section; acceptance; efficiency; total experimental systematic uncertainty, not including Monte Carlo statistics; observed CL using this region alone; expected CL using this region alone.

Cut flows for a representative selection of SUSY signal points in the Z-veto signal regions. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the larger of the two masses. The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise.

Cut flows for a representative selection of SUSY signal points in the Z-rich signal regions. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the larger of the two masses (or the value of mu in the case of GGM models). The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise.

Cut flows by lepton channel for a representative selection of SUSY signal points in the SR0noZa signal region. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the larger of the two masses. The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise. The RPC R-slepton model is used, with (m2,m1) = (450,300) GeV.

Cut flows by lepton channel for a representative selection of SUSY signal points in the SR1noZb signal region. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the larger of the two masses. The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise. The RPV gluino NLSP model is used, with lambda_133 != 0 and (m2,m1) = (800,400) GeV.

Cut flows by lepton channel for a representative selection of SUSY signal points in the SR0Z signal region. In each case, m2 and m1 refer to the axes of the plots in Sec. XI, where m2 is the value of mu. The number of events expected for a luminosity of 20.3 fb-1 is quoted at each step of the selection. The preselection requires four baseline leptons, at least two of which are light leptons; the signal lepton selection is made at the ``Lepton Multiplicity' stage. ``Event Cleaning' refers to the selection criteria applied to remove non-collision backgrounds and detector noise. The GGM tan beta = 30 model is used, with (m2,m1) = (200,1000) GeV.

Expected and observed $am_{T2}$ distribution for bCd_high1 SR, before applying the $am_{T2}>200$ GeV requirement. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Expected and observed leading b-jet $p_T$ distribution for bCd_high2 SR, before applying the leading b-jet $p_T>170$ GeV requirement. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Expected and observed $E_T^{miss}$ distribution for tNbC_mix SR, before applying the $E_T^{miss}>270$ GeV requirement. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Expected and observed lepton $p_T$ distribution for bCa_low SR. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Expected and observed lepton $p_T$ distribution for bCa_med SR. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Expected and observed $am_T2$ distribution for bCb_med1 SR. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Expected and observed $am_T2$ distribution for bCb_high SR. The uncertainty includes statistical and all experimental systematic uncertainties. The last bin includes overflows.

Best expected signal region for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1$ three-body scenario ($\tilde t_1\to bW\chi^0_1$). This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1$ four-body scenario ($\tilde t_1\to bff'\chi^0_1$). This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=150$ GeV. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=106$ GeV. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+5$ GeV. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\tilde t_1}-10$ GeV. This mapping is used for the final combined exclusion limits.

Best expected signal region for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\tilde t_1}=300$ GeV. This mapping is used for the final combined exclusion limits.

Upper limits on the model cross-section for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Observed exclusion contour for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Expected exclusion contour for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Upper limit on signal events for the $\tilde t_1$ three-body scenario ($\tilde t_1\to bW\chi^0_1$).

Observed exclusion contour for the $\tilde t_1$ three-body scenario ($\tilde t_1\to bW\chi^0_1$).

Expected exclusion contour for the $\tilde t_1$ three-body scenario ($\tilde t_1\to bW\chi^0_1$).

Upper limit on signal events for the $\tilde t_1$ four-body scenario ($\tilde t_1\to bff'\chi^0_1$).

Observed exclusion contour for the $\tilde t_1$ four-body scenario ($\tilde t_1\to bff'\chi^0_1$).

Expected exclusion contour for the $\tilde t_1$ four-body scenario ($\tilde t_1\to bff'\chi^0_1$).

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=150$ GeV.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=150$ GeV.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=150$ GeV.

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=106$ GeV.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=106$ GeV.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=106$ GeV.

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+5$ GeV.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+5$ GeV.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+5$ GeV.

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\tilde t_1}-10$ GeV.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\tilde t_1}-10$ GeV.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\tilde t_1}-10$ GeV.

Upper limit on signal events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\tilde t_1}=300$ GeV.

Observed exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\tilde t_1}=300$ GeV.

Expected exclusion contour for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\tilde t_1}=300$ GeV.

Acceptance of tN_diag SR ($E_T^{miss}>150$ GeV, $m_T>140$ GeV) for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of tN_med SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of tN_boost SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCb_med2 SR ($am_{T2}>250$ GeV, $m_T>60$ GeV) for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCc_diag SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCd_bulk SR ($am_{T2}>175$ GeV, $m_T>120$ GeV) for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCd_high1 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCd_high2 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCa_med for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCa_low for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCb_med1 for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of bCb_high for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of 3-body SR ($80<am_{T2}<90$ GeV, $m_T>120$ GeV) for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$). The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance of tNbC_mix SR for the asymmetric scenario ($\tilde t_1$, $\tilde t_1\to t\chi^0_1$, b $\chi^\pm_1$) with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Efficiency of tN_diag SR ($E_T^{miss}>150$ GeV, $m_T>140$ GeV) for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of tN_med SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of tN_boost SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCb_med2 SR ($am_{T2}>250$ GeV, $m_T>60$ GeV) for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCc_diag SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCd_bulk SR ($am_{T2}>175$ GeV, $m_T>120$ GeV) for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCd_high1 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCd_high2 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCa_med for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCa_low for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCb_med1 for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of bCb_high for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of 3-body SR ($80<am_{T2}<90$ GeV, $m_T>120$ GeV) for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$). The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency of tNbC_mix SR for the asymmetric scenario ($\tilde t_1$, $\tilde t_1\to t\chi^0_1$, b $\chi^\pm_1$) with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Number of generated events for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Number of generated events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Number of generated events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV; $E_T^{miss}$(gen)$>60$ GeV.

Number of generated events for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV; $E_T^{miss}$(gen)$>250$ GeV.

Number of generated events for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$).

Number of generated events for the asymmetric scenario ($\tilde t_1$, $\tilde t_1\to t\chi^0_1$, b $\chi^\pm_1$) with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Cross-section for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Cross-section for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Cross-section for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Cross-section for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$).

Cross-section for the asymmetric scenario ($\tilde t_1$, $\tilde t_1\to t\chi^0_1$, b $\chi^\pm_1$) with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Combined experimental systematic uncertainty of expected tN_diag SR yields for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$, using the 2 highest $E_T^{miss}$ and 2 highest $m_T$ bins.

Combined experimental systematic uncertainty of expected tN_med SR yields for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Combined experimental systematic uncertainty of expected tN_boost SR yields for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Combined experimental systematic uncertainty of expected bCb_med2 SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$, using the 2 highest $am_{T2}$ and 2 highest $m_T$ bins.

Combined experimental systematic uncertainty of expected bCc_diag SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Combined experimental systematic uncertainty of expected bCd_bulk SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$, using the 2 highest $am_{T2}$ and 2 highest $m_T$ bins.

Combined experimental systematic uncertainty of expected bCd_high1 SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Combined experimental systematic uncertainty of expected bCd_high2 SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Combined experimental systematic uncertainty of expected bCa_med SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Combined experimental systematic uncertainty of expected bCa_low SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Combined experimental systematic uncertainty of expected bCb_med1 SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Combined experimental systematic uncertainty of expected bCb_high SR yields for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Combined experimental systematic uncertainty of expected 3-body SR yields for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$), using the 2 lowest $am_{T2}$ and 2 highest $m_T$ bins.

Combined experimental systematic uncertainty of expected tNbC_mix SR yields for the asymmetric scenario ($\tilde t_1$, $\tilde t_1\to t\chi^0_1$, b $\chi^\pm_1$) with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed CLs in tN_diag SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Observed CLs in tN_med SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Observed CLs in tN_boost SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Observed CLs in bCb_med2 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed CLs in bCc_diag SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed CLs in bCd_bulk SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed CLs in bCd_high1 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed CLs in bCd_high2 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Observed CLs in bCa_med SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Observed CLs in bCa_low SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Observed CLs in bCb_med1 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Observed CLs in bCb_high SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Observed CLs in 3-body SR for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$).

Observed CLs in tNbC_mix SR for the mixed scenario (50% $\tilde t_1\to t\chi^0_1$, 50% $\tilde t_1\to b\chi^0_1$).

Expected CLs in tN_diag SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Expected CLs in tN_med SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Expected CLs in tN_boost SR for the $\tilde t_1\to t\chi^0_1$ scenario with $m_{\tilde t_1}>m_t+m_{\chi^0_1}$.

Expected CLs in bCb_med2 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Expected CLs in bCc_diag SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Expected CLs in bCd_bulk SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Expected CLs in bCd_high1 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Expected CLs in bCd_high2 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=2\times m_{\chi^0_1}$.

Expected CLs in bCa_med SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Expected CLs in bCa_low SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Expected CLs in bCb_med1 SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Expected CLs in bCb_high SR for the $\tilde t_1\to b\chi^\pm_1$ scenario with $m_{\chi^\pm_1}=m_{\chi^0_1}+20$ GeV.

Expected CLs in 3-body SR for the 3-body scenario ($\tilde t_1\to b W\chi^0_1$).

Expected CLs in tNbC_mix SR for the mixed scenario (50% $\tilde t_1\to t\chi^0_1$, 50% $\tilde t_1\to b\chi^\pm_1$).

Measured differential cross section with associated uncertainties as a function of absolute rapidity of diphoton system. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of multiplicity of jets with transverse momentum pT(jet) > 30 GeV. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of multiplicity of jets with transverse momentum pT(jet) > 30 GeV. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of multiplicity of jets with transverse momentum pT(jet) > 30 GeV. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of multiplicity of jets with transverse momentum pT(jet) > 30 GeV. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of the leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of the leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of absolute rapidity of the leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of absolute rapidity of the leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of scalar transverse momentum sum of all jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of scalar transverse momentum sum of all jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of second leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of second leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of rapidity separation of leading two jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of rapidity separation of leading two jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of azimuthal angle between diphoton and dijet systems. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of azimuthal angle between diphoton and dijet systems. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of azimuthal angle between the two leading jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of azimuthal angle between the two leading jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of diphoton thrust pT. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of diphoton thrust pT. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of rapidity separation of the two photons. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of rapidity separation of the two photons. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of tau of highest-tau jet (see paper for description). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of tau of highest-tau jet (see paper for description). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of scalar sum of tau for all jets (see paper for description). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of scalar sum of tau for all jets (see paper for description). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of absolute rapidity of second leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of absolute rapidity of second leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of third leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of third leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of dijet invariant mass. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of dijet invariant mass. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of the combined diphoton and dijet system. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of the combined diphoton and dijet system. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame in bins of diphoton transverse momentum. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame in bins of diphoton transverse momentum. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of diphoton transverse momentum in jet multiplicity bins. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of diphoton transverse momentum in jet multiplicity bins. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of leading jet transverse momentum for exclusive one-jet events. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of leading jet transverse momentum for exclusive one-jet events. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Diphoton kinematic acceptances in percent for gluon fusion for each fiducial region/variable bin studied in this paper, defined as the probability to fulfil the diphoton kinematic criteria: pT(2gam) < 0.35 (0.25) for the leading (subleading) photon and |eta(gam)|<2.37. The factors are evaluated using the POWHEG event generator with MPI modelling and hadronisation turned off. Consistent results for the diphoton variables are obtained by HRes 2.2. Uncertainties are taken from PDF variations. QCD scale varaitions have a negligible impact on these factors.

Diphoton kinematic acceptances in percent for gluon fusion for each fiducial region/variable bin studied in this paper, defined as the probability to fulfil the diphoton kinematic criteria: pT(2gam) < 0.35 (0.25) for the leading (subleading) photon and |eta(gam)|<2.37. The factors are evaluated using the POWHEG event generator with MPI modelling and hadronisation turned off. Consistent results for the diphoton variables are obtained by HRes 2.2. Uncertainties are taken from PDF variations. QCD scale varaitions have a negligible impact on these factors.

Isolation efficiencies in percent for gluon fusion H -> 2gam for each fiducial region/variable bin measured in this analysis. The isolation efficiency is defined as the probability for both photons to fulfil the isolation criteria (ETiso < 14 GeV as described in the text) for events that pass the diphoton kinematic criteria. Uncertainties are assigned in the same way as for the non-perturbative correction factors: by varying the fragmentation and underlying event modelling. These factors can be multiplied by the kinematic acceptance factors (see table~ ef{tab:fid_acceptance}) to extrapolate an inclusive gluon fusion Higgs prediction to the fiducial volume used in this analysis.

Isolation efficiencies in percent for gluon fusion H -> 2gam for each fiducial region/variable bin measured in this analysis. The isolation efficiency is defined as the probability for both photons to fulfil the isolation criteria (ETiso < 14 GeV as described in the text) for events that pass the diphoton kinematic criteria. Uncertainties are assigned in the same way as for the non-perturbative correction factors: by varying the fragmentation and underlying event modelling. These factors can be multiplied by the kinematic acceptance factors (see table~ ef{tab:fid_acceptance}) to extrapolate an inclusive gluon fusion Higgs prediction to the fiducial volume used in this analysis.

Non-perturbative correction factors in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Uncertainties are evaluated by deriving these factors using different generators and tunes as described in the text.

Non-perturbative correction factors in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Uncertainties are evaluated by deriving these factors using different generators and tunes as described in the text.

Fiducial cross sections in fb from WH, ZH and ttH combined, in each variable bin and fiducial region.

Fiducial cross sections in fb from WH, ZH and ttH combined, in each variable bin and fiducial region.

Statistical bin-to-bin correlations between the 24 bins used in the analysis described in Phys.Lett. B753 (2016) 69-85 (arXiv:1508.02507). See Figure 3 of this paper for the bins used.

Selection efficiencies in % for the signal in the top squark search, estimated from the simulation.

The final $S_{\text{T}}$ distribution for the leptoquark search with the $\text{e}\tau_{\text{h}}$ and $\mu\tau_{\text{h}}$ channels combined. The last bin contains the overflow events.

The final $S_{\text{T}}$ distribution for the top squark search with the $\text{e}\tau_{\text{h}}$ and $\mu\tau_{\text{h}}$ channels combined. The last bin contains the overflow events.

The expected and observed combined upper limits on the third-generation LQ pair production cross section $\sigma$ times the square of the branching fraction, $\mathcal{B}^2$, at 95% CL, as a function of the LQ mass. These limits also apply to top squarks decaying directly via the coupling $\lambda^{\prime}_{333}$.

The expected and observed combined upper limits on the top squark pair production cross section $\sigma$ times the square of the branching fraction, $\mathcal{B}^2$, at 95% CL, as a function of the top squark mass. These limits apply to top squarks with a chargino-mediated decay through the coupling $\lambda^{\prime}_{3jk}$.

The expected and observed 95% CL upper limits on the branching fraction for the leptoquark decay to a tau lepton and a bottom quark, as a function of the leptoquark mass. A search for top squark pair production [EPJC 73(2013)2677] has the same kinematic signature as the leptoquark decay to a tau neutrino and a top quark. This search is reinterpreted to provide the expected and observed 95% CL upper limits for low values of $\mathcal{B}$, assuming the leptoquark only decays to third-generation fermions.