The fraction of the direct photon component as a function of $p_T$.

The fraction of the direct photon component as a function of $p_T$.

Invariant yield (Au+Au) of direct photons as a function of $p_T$.

Invariant yield (Au+Au) of direct photons as a function of $p_T$.

Invariant yield (Au+Au) of direct photons as a function of $p_T$.

Invariant cross section ($p$+$p$) of direct photons as a function of $p_T$.

Observed and expected upper limits at 95% CL on the production cross-section times branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed and expected upper limits at 95% CL on the production cross-section times branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed and expected upper limits at 95% CL on the production cross-section times branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

IRing2 for HT2>=500 GeV, NJets>=5

IRing2 for HT2>=1000 GeV, NJets>=2

IRing2 for HT2>=1000 GeV, NJets>=3

IRing2 for HT2>=1000 GeV, NJets>=4

IRing2 for HT2>=1000 GeV, NJets>=5

IRing2 for HT2>=1500 GeV, NJets>=2

IRing2 for HT2>=1500 GeV, NJets>=3

IRing2 for HT2>=1500 GeV, NJets>=4

IRing2 for HT2>=1500 GeV, NJets>=5

IRing128 for HT2>=500 GeV, NJets>=2

IRing128 for HT2>=500 GeV, NJets>=3

IRing128 for HT2>=500 GeV, NJets>=4

IRing128 for HT2>=500 GeV, NJets>=5

IRing128 for HT2>=1000 GeV, NJets>=2

IRing128 for HT2>=1000 GeV, NJets>=3

IRing128 for HT2>=1000 GeV, NJets>=4

IRing128 for HT2>=1000 GeV, NJets>=5

IRing128 for HT2>=1500 GeV, NJets>=2

IRing128 for HT2>=1500 GeV, NJets>=3

IRing128 for HT2>=1500 GeV, NJets>=4

IRing128 for HT2>=1500 GeV, NJets>=5

ICyl16 for HT2>=500 GeV, NJets>=2

ICyl16 for HT2>=500 GeV, NJets>=3

ICyl16 for HT2>=500 GeV, NJets>=4

ICyl16 for HT2>=500 GeV, NJets>=5

ICyl16 for HT2>=1000 GeV, NJets>=2

ICyl16 for HT2>=1000 GeV, NJets>=3

ICyl16 for HT2>=1000 GeV, NJets>=4

ICyl16 for HT2>=1000 GeV, NJets>=5

ICyl16 for HT2>=1500 GeV, NJets>=2

ICyl16 for HT2>=1500 GeV, NJets>=3

ICyl16 for HT2>=1500 GeV, NJets>=4

ICyl16 for HT2>=1500 GeV, NJets>=5

IRing2 covariance for HT2>=500 GeV, NJets>=2 (Table 1)

IRing2 covariance for HT2>=500 GeV, NJets>=3 (Table 2)

IRing2 covariance for HT2>=500 GeV, NJets>=4 (Table 3)

IRing2 covariance for HT2>=500 GeV, NJets>=5 (Table 4)

IRing2 covariance for HT2>=1000 GeV, NJets>=2 (Table 5)

IRing2 covariance for HT2>=1000 GeV, NJets>=3 (Table 6)

IRing2 covariance for HT2>=1000 GeV, NJets>=4 (Table 7)

IRing2 covariance for HT2>=1000 GeV, NJets>=5 (Table 8)

IRing2 covariance for HT2>=1500 GeV, NJets>=2 (Table 9)

IRing2 covariance for HT2>=1500 GeV, NJets>=3 (Table 10)

IRing2 covariance for HT2>=1500 GeV, NJets>=4 (Table 11)

IRing2 covariance for HT2>=1500 GeV, NJets>=5 (Table 12)

IRing128 covariance for HT2>=500 GeV, NJets>=2 (Table 13)

IRing128 covariance for HT2>=500 GeV, NJets>=3 (Table 14)

IRing128 covariance for HT2>=500 GeV, NJets>=4 (Table 15)

IRing128 covariance for HT2>=500 GeV, NJets>=5 (Table 16)

IRing128 covariance for HT2>=1000 GeV, NJets>=2 (Table 17)

IRing128 covariance for HT2>=1000 GeV, NJets>=3 (Table 18)

IRing128 covariance for HT2>=1000 GeV, NJets>=4 (Table 19)

IRing128 covariance for HT2>=1000 GeV, NJets>=5 (Table 20)

IRing128 covariance for HT2>=1500 GeV, NJets>=2 (Table 21)

IRing128 covariance for HT2>=1500 GeV, NJets>=3 (Table 22)

IRing128 covariance for HT2>=1500 GeV, NJets>=4 (Table 23)

IRing128 covariance for HT2>=1500 GeV, NJets>=5 (Table 24)

ICyl16 covariance for HT2>=500 GeV, NJets>=2 (Table 25)

ICyl16 covariance for HT2>=500 GeV, NJets>=3 (Table 26)

ICyl16 covariance for HT2>=500 GeV, NJets>=4 (Table 27)

ICyl16 covariance for HT2>=500 GeV, NJets>=5 (Table 28)

ICyl16 covariance for HT2>=1000 GeV, NJets>=2 (Table 29)

ICyl16 covariance for HT2>=1000 GeV, NJets>=3 (Table 30)

ICyl16 covariance for HT2>=1000 GeV, NJets>=4 (Table 31)

ICyl16 covariance for HT2>=1000 GeV, NJets>=5 (Table 32)

ICyl16 covariance for HT2>=1500 GeV, NJets>=2 (Table 33)

ICyl16 covariance for HT2>=1500 GeV, NJets>=3 (Table 34)

ICyl16 covariance for HT2>=1500 GeV, NJets>=4 (Table 35)

ICyl16 covariance for HT2>=1500 GeV, NJets>=5 (Table 36)

IRing2 covariance, complete

1-IRing128 covariance, complete

1-ICyl16 covariance, complete

Pre-fit comparison between data and background in the baseline SR for two of the variables used as input for the BSM pBDT: HT.

Data and post-fit background comparison for the distributions of the discriminant variables fitted in the control regions obtained with the background-only fit: CR HF e.

Data and post-fit background comparison for the distributions of the discriminant variables fitted in the control regions obtained with the background-only fit: CR HF mu.

Data and post-fit background comparison for the distributions of the discriminant variables fitted in the control regions obtained with the background-only fit: CR Conv.

Data and post-fit background comparison for the distributions of the discriminant variables fitted in the control regions obtained with the background-only fit: CR ttW.

Data and post-fit background comparison for the distributions of the discriminant variables fitted in the control regions obtained with the background-only fit: CR lowBDT.

Data and post-fit background comparison obtained with the background-only fit in the validation region defined for ttW+jets events: number of jets.

Data and post-fit background comparison obtained with the background-only fit in the validation region defined for ttW+jets events: SM BDT.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 400 GeV.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 1000 GeV.

Observed (black solid line) and expected (black dashed line) 95% CL upper limits on the ttH/A cross-section times branching fraction of H/A->tt as a function of mH/A.

Observed (red line) and expected (black dashed line) exclusion regions at 95% CL in the tan(beta) versus mass plane assuming that both, a heavy scalar H and pseudo-scalar A, contribute to the tttt final state and have the same mass mH=mA.

Observed (red line) and expected (black dashed line) exclusion regions at 95% CL in the tan(beta) versus mass plane assuming that only the scalar H contributes to the tttt final state.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 500 GeV.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 600 GeV.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 700 GeV.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 800 GeV.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 900 GeV.

Observed (red line) and expected (black dashed line) exclusion regions at 95% CL in the tan(beta) versus mass plane assuming that only the pseudo-scalar A contributes to the tttt final state.

The product of acceptance and efficiency (A×ε), defined as the number of signal events satisfying the full set of selection criteria of the BSM SR, divided by the total number of generated signal events, for the tt̄H(→ tt̄ ) process in a type-II 2HDM for tanβ = 1 for the different heavy Higgs mass hypotheses.

2 sigma band of expected combination limits at 95% CL in the $(m_{a},m_{A})$ plane under the assumption of $sin\theta$ = 0.35.

Observed combination limits at 95% CL in the $(m_{a},m_{A})$ plane under the assumption of $sin\theta$ = 0.7.

Expected combination limits at 95% CL in the $(m_{a},m_{A})$ plane under the assumption of $sin\theta$ = 0.7.

1 sigma band of expected combination limits at 95% CL in the $(m_{a},m_{A})$ plane under the assumption of $sin\theta$ = 0.7.

2 sigma band of expected combination limits at 95% CL in the $(m_{a},m_{A})$ plane under the assumption of $sin\theta$ = 0.7.

Observed limits of H(inv) at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35.

Expected limits of H(inv) at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35.

Observed limits of monoh(bb) at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35.

Expected limits of monoh(bb) at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35.

Observed limits of $h\to aa \to \mu\mu\tau\tau$ at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35. Data points included are corresponding to the lower bound of the exclusion contour. The region above the bound was excluded.

Observed limits of $h\to aa \to bbbb$ at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35. Data points included are corresponding to the lower bound of the exclusion contour. The region above the bound was excluded.

Observed limits of $h\to aa \to bb\mu\mu$ at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35. Data points included are corresponding to the lower bound of the exclusion contour. The region above the bound was excluded.

Observed limits of $h\to aa \to 4 \mu$ at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35. Data points included are corresponding to the lower bound of the exclusion contour. The region above the bound was excluded.

Observed limits of $h\to aa \to 4 \mu$ at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35. Data points included are corresponding to the lower bound of the exclusion contour. The region above the bound was excluded.

Observed limits of $h\to aa \to 4 \mu$ at 95% CL in the $(m_{a},m_{\chi})$ plane under the assumption of $m_{A}$ = 1.2 TeV, $tan\beta$ = 1.0 and $sin\theta$ = 0.35. Data points included are corresponding to the lower bound of the exclusion contour. The region above the bound was excluded.

Observed combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

Expected combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

1 sigma band of expected combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

2 sigma band of expected combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

Observed combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

Expected combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

1 sigma band of expected combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

2 sigma band of expected combination limits at 95% CL in the $(m_{A},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

Observed combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

Expected combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

1 sigma band of expected combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

2 sigma band of expected combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.35.

Observed combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

Expected combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

1 sigma band of expected combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

2 sigma band of expected combination limits at 95% CL in the $(m_{a},tan\beta)$ plane under the assumption of $sin\theta$ = 0.7.

Observed and expected combination limits at 95% CL as a function of $sin\theta$ for the low-mass hypothesis $m_{A}$ = 0.6 TeV, $m_{a}$ = 200 GeV.

Observed and expected combination limits at 95% CL as a function of $sin\theta$ for the high-mass hypothesis $m_{A}$ = 1 TeV, $m_{a}$ = 350 GeV.

Observed and expected combination limits at 95% CL as a function of $m_{\chi}$ following the parameter choices of $m_{A}$ = 1 TeV, $m_{a}$ = 400 GeV.

Direct $\gamma$ invariant yields as a function of transverse momentum for 9 centrality selections and minimum bias Au+AU collisions at $\sqrt{s_{NN}}$ = 200 GeV. Data with no errors represents 90% confidence level upper limit. The bin range is not an uncertainty in the x-axis because the actual uncertainty by having the finite bin width is corrected for by the bin-shift correction. These bins were constructed using the corrected finite values as centers.

Ratio of Au+Au yield to $p$+$p$ yield normalized by the number of binary nucleon collisions as a function of centrality given by $N_{part}$ for direct $\gamma$ and $\pi^0$ yields integrated above 6 GeV/$c$. Data with no errors represents 90% confidence level upper limit.