Correlation matrix of all the parameters including nuisance parameters used in the profile likelihood unfolding fit on $A_{C}$.

Target asymmetry for $p\pi^0$ as a function of the polar angle for bins of the incident photon energy in the range of $E_\gamma$ = 650-2600 MeV.

Target asymmetry for $p\pi^0$ as a function of the $p\pi0$ invariant mass for bins of the incident photon energy in the range of $E_\gamma$ = 650-2600 MeV.

Target asymmetry for $p\pi^0$ as a function of the $\phi^*$ angle for bins of the incident photon energy in the range of $E_\gamma$ = 650-2600 MeV.

Double polarization observables for $\pi^0\pi^0$ as a function of the polar angle for bins of the incident photon energy in the range of $E_\gamma$ = 650-950 MeV.

Double polarization observables for $\pi^0\pi^0$ as a function of the $\pi^0\pi0$ invariant mass for bins of the incident photon energy in the range of $E_\gamma$ = 650-950 MeV.

Double polarization observables for $\pi^0\pi^0$ as a function of the polar angle for bins of the incident photon energy in the range of $E_\gamma$ = 650-950 MeV.

Double polarization observables for $p\pi^0$ as a function of the $p\pi0$ invariant mass for bins of the incident photon energy in the range of $E_\gamma$ = 650-950 MeV.

Double polarization observables for $\pi^0\pi^0$ as a function of the $\phi^*$ angle for bins of the incident photon energy in the range of $E_\gamma$ = 650-950 MeV.

Double polarization observables for $p\pi^0$ as a function of the $\phi^*$ angle for bins of the incident photon energy in the range of $E_\gamma$ = 650-950 MeV.

Target asymmetries for $\pi^0\pi^0$ as a function of $\cos\vartheta$, invariant mass and $\phi^*$ for bins of the incident photon energy in the range of $E_\gamma$ = 650-800 MeV.

Target asymmetries for $\pi^0\pi^0$ as a function of $\cos\vartheta$, invariant mass and $\phi^*$ for bins of the incident photon energy in the range of $E_\gamma$ = 800-950 MeV.

Target asymmetries for $\pi^0\pi^0$ as a function of $\cos\vartheta$, invariant mass and $\phi^*$ for bins of the incident photon energy in the range of $E_\gamma$ = 950-1100 MeV.

Target asymmetries for $\pi^0\pi^0$ as a function of $\cos\vartheta$, invariant mass and $\phi^*$ for bins of the incident photon energy in the range of $E_\gamma$ = 1100-1250 MeV.

Target asymmetries for $p\pi^0$ as a function of $\cos\vartheta$, invariant mass and $\phi^*$ for bins of the incident photon energy in the range of $E_\gamma$ = 650-800 MeV.

Target asymmetries for $p\pi^0$ as a function of $\cos\vartheta$, invariant mass and $\phi^*$ for bins of the incident photon energy in the range of $E_\gamma$ = 800-950 MeV.

Target asymmetries for $p\pi^0$ as a function of $\cos\vartheta$, invariant mass and $\phi^*$ for bins of the incident photon energy in the range of $E_\gamma$ = 950-1100 MeV.

Target asymmetries for $p\pi^0$ as a function of $\cos\vartheta$, invariant mass and $\phi^*$ for bins of the incident photon energy in the range of $E_\gamma$ = 1100-1250 MeV.

Double polarization observables for $\pi^0\pi^0$ as a function of $\cos\vartheta$, invariant mass and $\phi^*$ for bins of the incident photon energy in the range of $E_\gamma$ = 650-800 MeV.

Double polarization observables for $\pi^0\pi^0$ as a function of $\cos\vartheta$, invariant mass and $\phi^*$ for bins of the incident photon energy in the range of $E_\gamma$ = 800-950 MeV.

Double polarization observables for $p\pi^0$ as a function of $\cos\vartheta$, invariant mass and $\phi^*$ for bins of the incident photon energy in the range of $E_\gamma$ = 650-800 MeV.

Double polarization observables for $p\pi^0$ as a function of $\cos\vartheta$, invariant mass and $\phi^*$ for bins of the incident photon energy in the range of $E_\gamma$ = 800-950 MeV.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1500-1600 MeV.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 25° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 1° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV, allowing for all proton and photon laboratory angles.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1170 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1180 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1200 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1210 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1220 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1400-1420 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1460 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1460-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1500-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1540-1600 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{\pi^0\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\pi^0}$, with a cut on the structure in the $M^2_{p\eta}$ invariant mass distribution.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Differential cross section $d\sigma/dM_{\pi^0\eta}$ for $E_\gamma$ = 1540-1600 MeV.

$M_{p\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$M_{\pi^0\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$\theta_{p\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

$\theta_{\pi^0\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons with a Soeding-inspired analytic function --- statistical correlations only. Only one half of the (symmetric) matrix is stored.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons, differential in the dipion mass and in bins of $W$. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons, differential in the dipion mass and in bins of $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, in bins of the photon-proton energy $W$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross sections off protons as a function of energy. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross sections off protons as a function of energy --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic $\rho^0(770)$ photoproduction cross section off protons as a function of energy (various experiments)

Fit of elastic $\rho^0(770)$ photoproduction cross sections off protons as a function of energy (various experiments) --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction cross section off protons, double-differential in the dipion mass and the momentum transfer $\vert t\vert$. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction cross section off protons, double-differential in the dipion mass and the momentum transfer $\vert t\vert$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single-differential in the momentum transfer $\vert t\vert$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single differential in the momentum transfer $\vert t\vert$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction double-differenial cross section off protons, as a function of the dipion mass and the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$ The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction double-differenial cross section off protons, as a function of the dipion mass and the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single-differential in the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross section off protons, single differential in the momentum transfer $\vert t\vert$ and in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Regge fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ and photon-proton energy $W$. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Regge fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ and photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of $b$-slopes in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$.

Fit of $b$-slopes in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of Pomeron trajectories in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$.

Fit of Pomeron trajectories in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Measured $\gamma+c$ fiducial differential cross section as a function of $E_\text{T}^\gamma$ for $|\eta^\gamma|<1.37$.

Measured $\gamma+c$ fiducial differential cross section as a function of $E_\text{T}^\gamma$ for $1.56<|\eta^\gamma|<2.37$.

Measured ratio of the $\gamma+b$ fiducial differential cross section as a function of $E_\text{T}^\gamma$ for $|\eta^\gamma|<1.37$ to that for $1.56<|\eta^\gamma|<2.37$.

Measured ratio of the $\gamma+c$ fiducial differential cross section as a function of $E_\text{T}^\gamma$ for $|\eta^\gamma|<1.37$ to that for $1.56<|\eta^\gamma|<2.37$.

Statistical correlation between the $\gamma+b$ and the $\gamma+c$ cross sections in a given $E_\text{T}^\gamma$ bin and $|\eta^\gamma|$ region. The two cross sections are correlated as the heavy flavour fractions are extracted simultaneously from a template fit, performed in each $E_\text{T}^\gamma$ bin and separately for the two $|\eta^\gamma|$ regions.

Signed shifts of the individual systematic uncertainties on the $\gamma+b$ cross section for $|\eta^\gamma|<1.37$. The numbers after the name of the uncertainty source refer to the individual component in that uncertainty. Each bin of the MC statistical uncertainty is independent of any other bin. The first four components of the photon energy scale uncertainty are specific to this $|\eta^\gamma|$ region and are independent of the components in the other region. The region is indicated as part of their name to indicate the independence between the $|\eta^\gamma|$ regions. The uncertainties on the prompt photon modelling, non-perturbative QCD models and particle-level migration effects are only varied once and not up and down by their nature, but are symmetrised for the final results. Only uncertainties which have at least a 1% variation in at least one bin of the $\gamma+b$ and $\gamma+c$ cross section measurements, including the ratios, are listed. The others are summed in quadrature and listed as a single entry.

Signed shifts of the individual systematic uncertainties on the $\gamma+b$ cross section for $1.56<|\eta^\gamma|<2.37$. The numbers after the name of the uncertainty source refer to the individual component in that uncertainty. Each bin of the MC statistical uncertainty is independent of any other bin. The first four components of the photon energy scale uncertainty are specific to this $|\eta^\gamma|$ region and are independent of the components in the other region. The region is indicated as part of their name to indicate the independence between the $|\eta^\gamma|$ regions. The uncertainties on the prompt photon modelling, non-perturbative QCD models and particle-level migration effects are only varied once and not up and down by their nature, but are symmetrised for the final results. Only uncertainties which have at least a 1% variation in at least one bin of the $\gamma+b$ and $\gamma+c$ cross section measurements, including the ratios, are listed. The others are summed in quadrature and listed as a single entry.

Signed shifts of the individual systematic uncertainties on the $\gamma+c$ cross section for $|\eta^\gamma|<1.37$. The numbers after the name of the uncertainty source refer to the individual component in that uncertainty. Each bin of the MC statistical uncertainty is independent of any other bin. The first four components of the photon energy scale uncertainty are specific to this $|\eta^\gamma|$ region and are independent of the components in the other region. The region is indicated as part of their name to indicate the independence between the $|\eta^\gamma|$ regions. The uncertainties on the prompt photon modelling, non-perturbative QCD models and particle-level migration effects are only varied once and not up and down by their nature, but are symmetrised for the final results. Only uncertainties which have at least a 1% variation in at least one bin of the $\gamma+b$ and $\gamma+c$ cross section measurements, including the ratios, are listed. The others are summed in quadrature and listed as a single entry.

Signed shifts of the individual systematic uncertainties on the $\gamma+c$ cross section for $1.56<|\eta^\gamma|<2.37$. The numbers after the name of the uncertainty source refer to the individual component in that uncertainty. Each bin of the MC statistical uncertainty is independent of any other bin. The first four components of the photon energy scale uncertainty are specific to this $|\eta^\gamma|$ region and are independent of the components in the other region. The region is indicated as part of their name to indicate the independence between the $|\eta^\gamma|$ regions. The uncertainties on the prompt photon modelling, non-perturbative QCD models and particle-level migration effects are only varied once and not up and down by their nature, but are symmetrised for the final results. Only uncertainties which have at least a 1% variation in at least one bin of the $\gamma+b$ and $\gamma+c$ cross section measurements, including the ratios, are listed. The others are summed in quadrature and listed as a single entry.

$W$ dependence of the differential cross section ${\rm d}\sigma/{\rm d}|\cos\theta^*|$ in five $|\cos\theta^*|$ bins for $Q^2$=6.89 GeV$^2$.

$W$ dependence of the differential cross section ${\rm d}\sigma/{\rm d}|\cos\theta^*|$ in five $|\cos\theta^*|$ bins for $Q^2$=8.92 GeV$^2$.

$W$ dependence of the differential cross section ${\rm d}\sigma/{\rm d}|\cos\theta^*|$ in five $|\cos\theta^*|$ bins for $Q^2$=10.93 GeV$^2$.

$W$ dependence of the differential cross section ${\rm d}\sigma/{\rm d}|\cos\theta^*|$ in five $|\cos\theta^*|$ bins for $Q^2$=13.37 GeV$^2$.

$W$ dependence of the differential cross section ${\rm d}\sigma/{\rm d}|\cos\theta^*|$ in five $|\cos\theta^*|$ bins for $Q^2$=17.23 GeV$^2$.

$W$ dependence of the differential cross section ${\rm d}\sigma/{\rm d}|\cos\theta^*|$ in five $|\cos\theta^*|$ bins for $Q^2$=24.25 GeV$^2$.

$W$ dependence of the integrated cross section in nine $Q^2$ bins.