Correlations between event-by-event fluctuations of anisotropic flow harmonics are measured in Pb-Pb collisions at a center-of-mass energy per nucleon pair of 5.02 TeV, as recorded by the ALICE detector at the LHC. This study presents correlations up to the hexagonal flow harmonic, $v_6$, which was measured for the first time. The magnitudes of these higher-order correlations are found to vary as a function of collision centrality and harmonic order. These measurements are compared to viscous hydrodynamic model calculations with EKRT initial conditions and to the iEBE-VISHNU model with TRENTo initial conditions. The observed discrepancies between the data and the model calculations vary depending on the harmonic combinations. Due to the sensitivity of model parameters estimated with Bayesian analyses to these higher-order observables, the results presented in this work provide new and independent constraints on the initial conditions and transport properties in theoretical models used to describe the system created in heavy-ion collisions.
NSC(5,2) vs centrality in Pb-Pb collisions at 5.02 TeV
NSC(5,3) vs centrality in Pb-Pb collisions at 5.02 TeV
NSC(6,2) vs centrality in Pb-Pb collisions at 5.02 TeV
This paper presents a new $τ$-lepton reconstruction and identification procedure at the ATLAS detector at the Large Hadron Collider, which leads to significantly improved performance in the case of physics processes where a highly boosted pair of $τ$-leptons is produced and one $τ$-lepton decays into a muon and two neutrinos ($τ_μ$), and the other decays into hadrons and one neutrino ($τ_{had}$). By removing the muon information from the signals used for reconstruction and identification of the $τ_{had}$ candidate in the boosted pair, the efficiency is raised to the level expected for an isolated $τ_{had}$. The new procedure is validated by selecting a sample of highly boosted $Z\rightarrowτ_μτ_{had}$ candidates from the data sample of $140$${fb}^{-1}$ of proton-proton collisions at $13$ TeV recorded with the ATLAS detector. Good agreement is found between data and simulation predictions in both the $Z\rightarrowτ_μτ_{had}$ signal region and in a background validation region. The results presented in this paper demonstrate the effectiveness of the $τ_{had}$ reconstruction with muon removal in enhancing the signal sensitivity of the boosted $τ_μτ_{had}$ channel at the ATLAS detector.
The distribution of the TauID jet RNN score for $\tau_\mathrm{had}^{\mu\mkern-10mu\backslash}$ in the SR. `$Z(\rightarrow\tau\tau)$+jets' represents the contributions from the signal process. `Top' represents the predicted contributions from the $t\bar{t}$, single-top-quark, and $tW$ processes. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic sources.
The distribution of the TauID jet RNN score for $\tau_\mathrm{had}^{\mu\mkern-10mu\backslash}$ in the VR. `$Z(\rightarrow\tau\tau)$+jets' represents the contributions from the signal process. `Top' represents the predicted contributions from the $t\bar{t}$, single-top-quark, and $tW$ processes. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic sources.
The distribution of the $p_\mathrm{T}{}_{\mu\mathrm{-had}}^\mathrm{col}$ in the SR. `$Z(\rightarrow\tau\tau)+\text{jets}$' represents the contributions from the signal process. `Diboson' indicates the contributions from $WW$, $WZ$, and $ZZ$ processes. `Top' represents the predicted contributions from the $t\bar{t}$, single-top-quark, and $tW$ processes. `Other' includes the contributions from the $Z(\rightarrow\ell\ell)$+jets, $W$+jets, and Higgs boson processes. The uncertainties shown include both statistical and systematic sources.
A search for exotic decays of the 125 GeV Higgs boson into a pair of new spin-0 particles, $H \to aa$, where one decays into a photon pair and the other into a $\tau$-lepton pair, is presented. Hadronic decays of the $\tau$-leptons are considered and reconstructed using a dedicated tagger for collimated $\tau$-lepton pairs. The search uses 140 fb$^{-1}$ of proton-proton collision data at a centre-of-mass energy of $\sqrt{s}=13$ TeV recorded between 2015 and 2018 by the ATLAS experiment at the Large Hadron Collider. The search is performed in the mass range of the $a$ boson between 10 GeV and 60 GeV. No significant excess of events is observed above the Standard Model background expectation. Model-independent upper limits at 95$\% $ confidence level are set on the branching ratio of the Higgs boson to the $\gamma\gamma\tau\tau$ final state, $\mathcal{B}(H\to aa\to \gamma\gamma\tau\tau)$, ranging from 0.2$\% $ to 2$\% $, depending on the $a$-boson mass hypothesis.
Distribution of the diphoton invariant mass for all events satisfying the analysis selections in the full Run 2 dataset.
Scan of the observed $p$-value as a function of $m_{a}$ for the background-only hypothesis.
The observed and expected ($\pm1\sigma$) upper limits at 95% CL on the branching ratio for $H\rightarrow aa\rightarrow \gamma\gamma\tau\tau$ as a function of the resonance mass hypothesis $m_{a}$.
A measurement of off-shell Higgs boson production in the $H^*\to ZZ\to 4\ell$ decay channel is presented. The measurement uses 140 fb$^{-1}$ of proton-proton collisions at $\sqrt{s}=13$ TeV collected by the ATLAS detector at the Large Hadron Collider and supersedes the previous result in this decay channel using the same dataset. The data analysis is performed using a neural simulation-based inference method, which builds per-event likelihood ratios using neural networks. The observed (expected) off-shell Higgs boson production signal strength in the $ZZ\to 4\ell$ decay channel at 68% CL is $0.87^{+0.75}_{-0.54}$ ($1.00^{+1.04}_{-0.95}$). The evidence for off-shell Higgs boson production using the $ZZ\to 4\ell$ decay channel has an observed (expected) significance of $2.5\sigma$ ($1.3\sigma$). The expected result represents a significant improvement relative to that of the previous analysis of the same dataset, which obtained an expected significance of $0.5\sigma$. When combined with the most recent ATLAS measurement in the $ZZ\to 2\ell 2\nu$ decay channel, the evidence for off-shell Higgs boson production has an observed (expected) significance of $3.7\sigma$ ($2.4\sigma$). The off-shell measurements are combined with the measurement of on-shell Higgs boson production to obtain constraints on the Higgs boson total width. The observed (expected) value of the Higgs boson width at 68% CL is $4.3^{+2.7}_{-1.9}$ ($4.1^{+3.5}_{-3.4}$) MeV.
Values of the test statistic $t_{\mu_{\mathrm{off-shell}}}$ assuming a single parameter of interest $\mu_{\mathrm{off-shell}}$ obtained with an Asimov dataset and with data in the $H^*\rightarrow ZZ\rightarrow 4\ell$ decay channel. The values from the histogram-based analysis (Phys. Lett. B 846 (2023) 138223) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.
Values of the test statistic $t_{\mu_{\mathrm{off-shell}}}$ assuming a single parameter of interest $\mu_{\mathrm{off-shell}}$ obtained with an Asimov dataset and with data in the $H^*\rightarrow ZZ\rightarrow 4\ell$ decay channel. The values with all nuisance parameters fixed at their best-fit values (stat-only) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.
Values of the test statistic $t_{\mu_{\mathrm{off-shell}}}$ assuming a single parameter of interest $\mu_{\mathrm{off-shell}}$ obtained with an Asimov dataset and with data when combining the $H^*\rightarrow ZZ\rightarrow 4\ell$ and $H^*\rightarrow ZZ\rightarrow 2\ell 2\nu$ decay channels. The values with all nuisance parameters fixed at their best-fit values (stat-only) are added for comparison. The 68% and 95% confidence intervals obtained from the Neyman construction are also added.
A search for decays of the Higgs boson into a $Z$ boson and a light resonance, with a mass of 0.5-3.5 GeV, is performed using the full 140 fb$^{-1}$ dataset of 13 TeV proton-proton collisions recorded by the ATLAS detector during Run 2 of the LHC. Leptonic decays of the $Z$ boson and hadronic decays of the light resonance are considered. The resonance can be interpreted as a $J/ψ$ or $η_c$ meson, an axion-like particle, or a light pseudoscalar in two-Higgs-doublet models. Due to its low mass, it would be produced with high boost and reconstructed as a single small-radius jet of hadrons. A neural network is used to correct the Monte Carlo simulation of the background in a data-driven way. Two additional neural networks are used to distinguish signal from background. A binned profile-likelihood fit is performed on the final-state invariant mass distribution. No significant excess of events relative to the expected background is observed, and upper limits at 95% confidence level are set on the Higgs boson's branching fraction to a $Z$ boson and a light resonance. The exclusion limit is ~10% for the lower masses, and increases for higher masses. Upper limits on the effective coupling $C^\text{eff}_{ZH}/Λ$ of an axion-like particle to a Higgs boson and $Z$ boson are also set at 95% confidence level, and range from 0.9 to 2 TeV$^{-1}$.
The angularity, for data, background (pre- and post-reweighting) and three $H\rightarrow Za$ signal hypotheses (for $a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection but not the classification NN requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.
The angularity, for data, background (pre- and post-reweighting) and three $H\rightarrow Za$ signal hypotheses (for $a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection but not the classification NN requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.
The modified energy correlation function, for data, background (pre- and post-reweighting) and three $H\rightarrow Za$ signal hypotheses (for $a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection but not the classification NN requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.
In this letter, we present the first measurement of direct photons at the transverse momentum of $ 1 < p_{\rm T} < 6$ GeV/$c$ at midrapidity $|η| < 0.8$ in inelastic and high-multiplicity proton--proton collisions at a centre-of-mass energy of $\sqrt{s} =$ 13 TeV. The fraction of virtual direct photons in the inclusive virtual photon spectrum is obtained from a fit to the dielectron invariant mass spectrum. In the limit of zero invariant mass, this fraction is equal to the relative contribution of real direct photons in the inclusive real photon spectrum. Contributions from decays of light-flavour neutral mesons are estimated using independent measurements in proton-proton collisions at the same energy and the same event class. For the first time at the LHC energies, a direct-photon signal is observed at low $p_{\rm T}$ in both inelastic and high-multiplicity event classes, with a significance of 3.2$σ$ and 1.9$σ$ in terms of standard deviations, correspondingly. The yield of direct photons in inelastic pp collisions is compared to perturbative QCD calculations. The integrated photon yield is studied as a function of charged-particle multiplicity and is compared to the results from other experiments and theoretical calculations. The results show a significant increase of direct-photon yield with charged-particle multiplicity.
The dielectron cross section in inelastic pp collisions at $\sqrt{s}$ = 13 TeV as a function of invariant mass for 2 < $p_{\rm T,ee}$ < 3 GeV/$c$.
The dielectron cross section in high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV as a function of invariant mass for 2 < $p_{\rm T,ee}$ < 3 GeV/$c$.
The direct-photon fraction r in inelastic pp collisions at $\sqrt{s}$ = 13 TeV as a function of transverse momentum for 1 < $p_{\rm T}$ < 6 GeV/$c$. r is the ratio of direct GAMMA to inclusive GAMMA.
Top-quark pair production is observed in lead-lead (Pb+Pb) collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV at the Large Hadron Collider with the ATLAS detector. The data sample was recorded in 2015 and 2018, amounting to an integrated luminosity of 1.9 nb$^{-1}$. Events with exactly one electron and one muon and at least two jets are selected. Top-quark pair production is measured with an observed (expected) significance of 5.0 (4.1) standard deviations. The measured top-quark pair production cross-section is $\sigma_{t\bar{t}} = 3.6\;^{+1.0}_{-0.9}\;\mathrm{(stat.)}\;^{+0.8}_{-0.5}\;\mathrm{(syst.)} ~\mathrm{\mu b}$, with a total relative uncertainty of 31%, and is consistent with theoretical predictions using a range of different nuclear parton distribution functions. The observation of this process consolidates the evidence of the existence of all quark flavors in the pre-equilibrium stage of the quark-gluon plasma at very high energy densities, similar to the conditions present in the early universe.
The figure shows the post-fit distribution of events as a function of the dilepton invariant mass ($m_{e\mu}$), in lead-lead (Pb+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV, with an integrated luminosity of 1.9 nb$^{-1}$. The data correspond to the SR1 (Signal Region 1 (SR\(_1\)):} Events with exactly one muon and one oppositely charged electron, a dilepton invariant mass \( m_{e\mu} \geq 30 \, \mathrm{GeV} \), at least two jets with \( p_T \geq 35 \, \mathrm{GeV} \), and a dilepton transverse momentum \( p_T^{e\mu} > 40 \, \mathrm{GeV} \). This region is expected to be signal-dominated) channel in a pre-fit configuration. The stacked histograms represent different processes contributing to the event yield.
The figure shows the post-fit distribution of events as a function of the dilepton invariant mass ($m_{e\mu}$), in lead-lead (Pb+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV, with an integrated luminosity of 1.9 nb$^{-1}$. The data correspond to the SR2 (Signal Region 2 (SR\(_2\)):} Events meeting the same criteria as SR\(_1\), but with a dilepton transverse momentum \( p_T^{e\mu} \leq 40 \, \mathrm{GeV} \). This region includes events with a lower \( p_T^{e\mu} \) and has a larger background contribution) channel in a pre-fit configuration. The stacked histograms represent different processes contributing to the event yield.
The impact of systematic uncertainties on the fitted signal-strength parameter $\hat{\mu}$ for the combined fit of all channels. Only the 10 most significant systematic uncertainties are shown and listed in decreasing order of their impact on $\mu$ on the $y$-axis. The empty (filled) blue/cyan boxes correspond to the pre-fit (post-fit) impact on $\mu$, referring to the upper $x$-axis. The impact of each systematic uncertainty, $\Delta \mu$, is calculated by comparing the nominal best-fit value of $\mu$ with the result of the fit when fixing the corresponding nuisance parameter $\theta$ to its best-fit value $\hat{\theta}$ shifted by its pre-fit (post-fit) uncertainties $\hat{\theta} \pm \Delta \theta(\hat{\theta} \pm \Delta \hat{\theta})$. The black points, which refer to the lower $x$-axis, show the pulls of the fitted nuisance parameters, i.e., the deviations of the fitted parameters $\hat{\theta}$ from their nominal values $\theta_0$, normalized to their nominal uncertainties $\Delta \theta$. The black lines show the post-fit uncertainties of the nuisance parameters, relative to their nominal uncertainties, which are indicated by the dashed lines.
A measurement of the $B^{0}$ meson lifetime and related properties using $B^0 \to J/ψK^{*0}$ decays in data from 13 TeV proton-proton collisions with an integrated luminosity of 140 fb$^{-1}$ recorded by the ATLAS detector at the LHC is presented. The measured effective lifetime is $$ τ= 1.5053 \pm 0.0012 ~\mathrm{(stat.)} \pm 0.0035 ~\mathrm{(syst.)~ps}. $$ The average decay width extracted from the effective lifetime, using parameters from external sources, is $$ Γ_d = 0.6639 \pm 0.0005 ~\mathrm{(stat.)} \pm 0.0016 ~\mathrm{(syst.)}\pm 0.0038 ~\textrm{(ext.)} \textrm{ ps}^{-1}, $$ where the uncertainties are statistical, systematic and from external sources. The earlier ATLAS measurement of $Γ_s$ in the $B^0_s \to J/ψϕ$ decay was used to derive a value for the ratio of the average decay widths $Γ_d$ and $Γ_s$ for $B^{0}$ and $B_s^{0}$ mesons respectively, of $$ \frac{Γ_d }{Γ_s } = 0.9905 \pm 0.0022 ~\textrm{(stat.)} \pm 0.0036 ~\textrm{(syst.)} \pm 0.0057 ~\textrm{(ext.)}. $$ The measured lifetime, average decay width and decay width ratio are in agreement with theoretical predictions and with measurements by other experiments. This measurement provides the most precise result of the effective lifetime of the $B^{0}$ meson to date.
The measured effective lifetime for the $B^0 \rightarrow J/\psi\,K^{*0}$ decay.
The measured average decay width $\Gamma_{d}\,$ extracted from the average lifetime.
The measured ratio $\Gamma_{d} / \Gamma_{s}\,$ of the average decay widths.
The distribution of angular correlations between prompt charm hadrons and primary charged particles in pp collisions is sensitive to the charm-quark hadronisation process. In this letter, charm-baryon correlations are measured for the first time by studying the azimuthal-angle difference between charged particles and prompt $Λ_{\rm c}^{+}$ baryons produced in pp collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV, with the ALICE detector. $Λ_{\rm c}^{+}$ baryons are reconstructed at midrapidity ($|y| < 0.5$) in the transverse-momentum interval $3 < p_{\rm T} < 16$ GeV/$c$, and correlated with charged particles with $p_{\rm T} > 0.3$ GeV/$c$ and pseudorapidity $|η| < 0.8$. For $3< p_{\rm T}^{Λ_{\rm c}^{+},{\rm D}} <5$ GeV/$c$, the comparison with published measurements of D-meson and charged-particle correlations in the same collision system hints at a larger number of low-momentum particles associated with $Λ_{\rm c}^{+}$-baryon triggers than with D-meson triggers, both in the collinear and opposite directions with respect to the trigger particle. These differences can be quantified by the comparison of the properties of the near- and away-side correlation peaks, and are not reproduced by predictions of various Monte Carlo event generators, generally underpredicting the associated particle yields at $p_{\rm T}^{\rm assoc}<1$ GeV/$c$. This tension between $Λ_{\rm c}^{+}$-baryon and D-meson associated peak yields could suggest a modified fragmentation of the charm quark, or a different hadronisation process, when a charm baryon is produced in the final state.
Azimuthal-correlation distributions of $\Lambda_\mathrm{c}^{+}$ baryons with associated particles in pp collisions at $\sqrt{s} = 13$ TeV, compared to the average of the azimuthal-correlation distributions of D$^{0}$, D$^{+}$, and D$^{*+}$ mesons with associated particles in the same collision system. The distributions are shown after the baseline subtraction. Charm hadrons are reported for $3 < p_{\rm T} < 5$ GeV/$c$, associated particles for $0.3 < p_{\rm T} < 1$ GeV/$c$. The rapidity range for the D mesons and $\Lambda_\mathrm{c}^{+}$ baryons is $|y^{\rm D,\Lambda_\mathrm{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.
Azimuthal-correlation distributions of $\Lambda_\mathrm{c}^{+}$ baryons with associated particles, i.e. primary charged particles, in pp collisions at $\sqrt{s} = 13$ TeV, compared to the average of the azimuthal-correlation distributions of D$^{0}$, D$^{+}$, and D$^{*+}$ mesons with associated particles in the same collision system. The distributions are shown after the baseline subtraction. Charm hadrons are reported for $5 < p_{\rm T} < 8$ GeV/$c$, associated particles for $0.3 < p_{\rm T} < 1$ GeV/$c$. The rapidity range for the D mesons and $\Lambda_\mathrm{c}^{+}$ baryons is $|y^{\rm D,\Lambda_\mathrm{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.
Azimuthal-correlation distributions of $\Lambda_\mathrm{c}^{+}$ baryons with associated particles, i.e. primary charged particles, in pp collisions at $\sqrt{s} = 13$ TeV, compared to the average of the azimuthal-correlation distributions of D$^{0}$, D$^{+}$, and D$^{*+}$ mesons with associated particles in the same collision system. The distributions are shown after the baseline subtraction. Charm hadrons are reported for $8 < p_{\rm T} < 16$ GeV/$c$, associated particles for $p_{\rm T} > 1$ GeV/$c$. The rapidity range for the D mesons and $\Lambda_\mathrm{c}^{+}$ baryons is $|y^{\rm D,\Lambda_\mathrm{c}^{+}}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D,\Lambda_{c}^{+}}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.
The momentum-differential invariant cross sections of ${\pi^{0}}$ and $\eta$ mesons are reported for pp collisions at $\sqrt{s}$ = 13 TeV at midrapidity ($|y|<0.8$). The measurement is performed in a broad transverse-momentum range of $0.2<p_{\rm T}<200$ GeV/$c$ and $0.4<p_{\rm T}<60$ GeV/$c$ for the ${\pi^{0}}$ and $\eta$, respectively, extending the $p_{\rm T}$ coverage of previous measurements. Transverse-mass-scaling violation of up to 60% at low transverse momentum has been observed, agreeing with measurements at lower collision energies. Transverse Bjorken $x$ ($x_{\rm T}$) scaling of the ${\pi^{0}}$ cross sections at LHC energies is fulfilled with a power-law exponent of $n$ = 5.01$\pm$0.05, consistent with values obtained for charged pions at similar collision energies. The data is compared to predictions from next-to-leading order perturbative QCD calculations, where the ${\pi^{0}}$ spectrum is best described using the CT18 parton distribution function and the NNFF1.0 or BDSS fragmentation function. Expectations from PYTHIA8 and EPOS LHC overestimate the spectrum for the ${\pi^{0}}$ and are not able to describe the shape and magnitude of the $\eta$ spectrum. The charged-particle multiplicity dependent ${\pi^{0}}$ and $\eta$ $p_{\rm T}$ spectra show the expected hardening with increasing multiplicity. This is demonstrated across a broad transverse-momentum range and up to events with a charged-particle multiplicity exceeding five times the mean value in minimum bias collisions. The absolute magnitude of the $\eta/\pi^{0}$ ratio shows a dependence on the charged-particle multiplicity for $p_{\rm T}<4$ GeV/$c$, qualitatively described by PYTHIA8 and EPOS LHC due to a rising contribution from feed-down of heavier particles to the ${\pi^{0}}$ spectrum.
Invariant differential cross section of the $\pi^{0}$ versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV.
Invariant differential cross section of the $\eta$ meson versus transverse momentum for pp collisions at $\sqrt{s}$ = 13 TeV.
The $\eta/\pi^{0}$ ratio as a function of $p_{\rm T}$ for pp collisions at $\sqrt{s}$ = 13 TeV.