A measurement of dielectron production in proton-proton (pp) collisions at $\sqrt{s} = 13$ TeV, recorded with the ALICE detector at the CERN LHC, is presented in this Letter. The data set was recorded with a reduced magnetic solenoid field. This enables the investigation of a kinematic domain at low dielectron invariant mass $m_{\rm ee}$ and pair transverse momentum $p_{\rm T,ee}$ that was previously inaccessible at the LHC. The cross section for dielectron production is studied as a function of $m_{\rm ee}$, $p_{\rm T,ee}$, and event multiplicity ${\rm d} N_{\rm ch}/{\rm d} \eta$. The expected dielectron rate from hadron decays, called hadronic cocktail, utilizes a parametrization of the measured $\eta/\pi^0$ ratio in pp and proton-nucleus (p-A) collisions, assuming that this ratio shows no strong dependence on collision energy at low transverse momentum. Comparison of the measured dielectron yield to the hadronic cocktail at $0.15<m_{\rm ee}<0.6$ GeV/$c^2$ and for $p_{\rm T,ee} < 0.4$ GeV/$c$ indicates an enhancement of soft dielectrons, reminiscent of the 'anomalous' soft-photon and -dilepton excess in hadron-hadron collisions reported by several experiments under different experimental conditions. The enhancement factor over the hadronic cocktail amounts to $1.61\pm 0.13\,(\rm{stat.})\pm 0.17\,(\rm{syst., data}) \pm 0.34\,(\rm{syst., cocktail})$ in the ALICE acceptance. Acceptance-corrected excess spectra in $m_{\rm ee}$ and $p_{\rm T,ee}$ are extracted and compared with calculations of dielectron production from hadronic bremsstrahlung and thermal radiation within a hadronic many-body approach.

Measurements of two- and multi-particle angular correlations in pp collisions at sqrt(s) = 5, 7, and 13 TeV are presented as a function of charged-particle multiplicity. The data, corresponding to integrated luminosities of 1.0 inverse picobarn (5 TeV), 6.2 inverse picobarns (7 TeV), and 0.7 inverse picobarns (13 TeV), were collected using the CMS detector at the LHC. The second-order (v[2]) and third-order (v[3]) azimuthal anisotropy harmonics of unidentified charged particles, as well as v[2] of K0 short and Lambda/anti-Lambda particles, are extracted from long-range two-particle correlations as functions of particle multiplicity and transverse momentum. For high-multiplicity pp events, a mass ordering is observed for the v[2] values of charged hadrons (mostly pions), K0 short, and Lambda/anti-Lambda, with lighter particle species exhibiting a stronger azimuthal anisotropy signal below pt of about 2 GeV/c. For 13 TeV data, the v[2] signals are also extracted from four- and six-particle correlations for the first time in pp collisions, with comparable magnitude to those from two-particle correlations. These observations are similar to those seen in pPb and PbPb collisions, and support the interpretation of a collective origin for the observed long-range correlations in high-multiplicity pp collisions.

The production cross-sections of $J/\psi$ mesons in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=5$ TeV are measured using a data sample corresponding to an integrated luminosity of $9.13\pm0.18~\text{pb}^{-1}$, collected by the LHCb experiment. The cross-sections are measured differentially as a function of transverse momentum, $p_{\text{T}}$, and rapidity, $y$, and separately for $J/\psi$ mesons produced promptly and from beauty hadron decays (nonprompt). With the assumption of unpolarised $J/\psi$ mesons, the production cross-sections integrated over the kinematic range $0<p_{\text{T}}<20~\text{GeV}/c$ and $2.0<y<4.5$ are $8.154\pm0.010\pm0.283~\mu\text{b}$ for prompt $J/\psi$ mesons and $0.820\pm0.003\pm0.034~\mu\text{b}$ for nonprompt $J/\psi$ mesons, where the first uncertainties are statistical and the second systematic. These cross-sections are compared with those at $\sqrt{s}=8$ TeV and $13$ TeV, and are used to update the measurement of the nuclear modification factor in proton-lead collisions for $J/\psi$ mesons at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\text{NN}}}=5$ TeV. The results are compared with theoretical predictions.

The production cross section of inclusive J/$\psi$ pairs in pp collisions at a centre-of-mass energy $\sqrt{s} = 13$ TeV is measured with ALICE. The measurement is performed for J/$\psi$ in the rapidity interval $2.5 < y < 4.0$ and for transverse momentum $p_{\rm T} > 0$. The production cross section of inclusive J/$\psi$ pairs is reported to be $10.3 \pm 2.3 {\rm (stat.)} \pm 1.3 {\rm (syst.)}$ nb in this kinematic interval. The contribution from non-prompt J/$\psi$ (i.e. originated from beauty-hadron decays) to the inclusive sample is evaluated. The results are discussed and compared with data.

Understanding the production mechanism of light (anti)nuclei is one of the key challenges of nuclear physics and has important consequences for astrophysics, since it provides an input for indirect dark-matter searches in space. In this paper, the latest results about the production of light (anti)nuclei in pp collisions at $\sqrt{s} = 13$ TeV are presented, focusing on the comparison with the predictions of coalescence and thermal models. For the first time, the coalescence parameters $B_2$ for deuterons and $B_3$ for helions are compared with parameter-free theoretical predictions that are directly constrained by the femtoscopic measurement of the source radius in the same event class. A fair description of the data with a Gaussian wave function is observed for both deuteron and helion, supporting the coalescence mechanism for the production of light (anti)nuclei in pp collisions. This method paves the way for future investigations of the internal structure of more complex nuclear clusters, including the hypertriton.

A precision measurement of the $Z$ boson production cross-section at $\sqrt{s} = 13$ TeV in the forward region is presented, using $pp$ collision data collected by the LHCb detector, corresponding to an integrated luminosity of 5.1 fb$^{-1}$. The production cross-section is measured using $Z\rightarrow\mu^+\mu^-$ events within the fiducial region defined as pseudorapidity $2.0<\eta<4.5$ and transverse momentum $p_{T}>20$ GeV/$c$ for both muons and dimuon invariant mass $60<M_{\mu\mu}<120$ GeV/$c^2$. The integrated cross-section is determined to be $\sigma (Z \rightarrow \mu^+ \mu^-)$ = 196.4 $\pm$ 0.2 $\pm$ 1.6 $\pm$ 3.9~pb, where the first uncertainty is statistical, the second is systematic, and the third is due to the luminosity determination. The measured results are in agreement with theoretical predictions within uncertainties.

Results are presented from a search for CP violation in top quark pair production, using proton-proton collisions at a center-of-mass energy of 13 TeV. The data used for this analysis consist of final states with two charged leptons collected by the CMS experiment, and correspond to an integrated luminosity of 35.9 fb$^{-1}$. The search uses two observables, $\mathcal{O}_1$ and $\mathcal{O}_3$, which are Lorentz scalars. The observable $\mathcal{O}_1$ is constructed from the four-momenta of the charged leptons and the reconstructed top quarks, while $\mathcal{O}_3$ consists of the four-momenta of the charged leptons and the b quarks originating from the top quarks. Asymmetries in these observables are sensitive to CP violation, and their measurement is used to determine the chromoelectric dipole moment of the top quark. The results are consistent with the expectation from the standard model.

Cross-sections for the production of a $Z$ boson in association with two photons are measured in proton$-$proton collisions at a centre-of-mass energy of 13 TeV. The data used correspond to an integrated luminosity of 139 fb$^{-1}$ recorded by the ATLAS experiment during Run 2 of the LHC. The measurements use the electron and muon decay channels of the $Z$ boson, and a fiducial phase-space region where the photons are not radiated from the leptons. The integrated $Z(\rightarrow\ell\ell)\gamma\gamma$ cross-section is measured with a precision of 12% and differential cross-sections are measured as a function of six kinematic variables of the $Z\gamma\gamma$ system. The data are compared with predictions from MC event generators which are accurate to up to next-to-leading order in QCD. The cross-section measurements are used to set limits on the coupling strengths of dimension-8 operators in the framework of an effective field theory.

We report on the inclusive J/$\psi$ production cross section measured at the CERN Large Hadron Collider in proton-proton collisions at a centre-of-mass energy $\sqrt{s}~=~13$ TeV. The J/$\psi$ mesons are reconstructed in the $\rm e^{+} e^{-}$ decay channel and the measurements are performed at midrapidity ($|y|<0.9$) in the transverse-momentum interval $0< p_{\rm T} <40$ GeV/$c$, using a minimum bias data sample corresponding to an integrated luminosity $L_{\text{int}} = 32.2~\text{nb}^{-1}$ and an Electromagnetic Calorimeter triggered data sample with $L_{\text{int}} = 8.3~\mathrm{pb}^{-1}$. The $p_{\rm T}$-integrated J/$\psi$ production cross section at midrapidity, computed using the minimum bias data sample, is $\text{d}\sigma/\text{d}y|_{y=0} = 8.97\pm0.24~(\text{stat})\pm0.48~(\text{syst})\pm0.15~(\text{lumi})~\mu\text{b}$. An approximate logarithmic dependence with the collision energy is suggested by these results and available world data, in agreement with model predictions. The integrated and $p_{\rm T}$-differential measurements are compared with measurements in pp collisions at lower energies and with several recent phenomenological calculations based on the non-relativistic QCD and Color Evaporation models.

The pseudorapidity density of charged particles with minimum transverse momentum ($p_{\rm T}$) thresholds of 0.15, 0.5, 1, and 2 GeV$/c$ is measured in pp collisions at the centre of mass energies of $\sqrt{s} =$ 5.02 and 13 TeV with the ALICE detector. The study is carried out for inelastic collisions with at least one primary charged particle having a pseudorapidity ($\eta$) within $\pm0.8$ and $p_{\rm T}$ larger than the corresponding threshold. In addition, measurements without $p_{\rm T}$-thresholds are performed for inelastic and non-single-diffractive events as well as for inelastic events with at least one charged particle having $|\eta|<1$ in pp collisions at $\sqrt{s} =$ 5.02 TeV for the first time at the LHC. These measurements are compared to the PYTHIA 6, PYTHIA 8, and EPOS-LHC models. In general, the models describe the $\eta$ dependence of particle production well. However, discrepancies are observed for the highest transverse momentum threshold ($p_{\rm T}>2 {\rm\ GeV}/c$), highlighting the importance of such measurements for tuning event generators. The new measurements agree within uncertainties with results from the ATLAS and CMS experiments obtained at $\sqrt{s} = 13$ TeV.

Differential cross-section measurements of $Z\gamma$ production in association with hadronic jets are presented, using the full 139 fb$^{-1}$ dataset of $\sqrt{s}=13$ TeV proton-proton collisions collected by the ATLAS detector during Run 2 of the LHC. Distributions are measured using events in which the $Z$ boson decays leptonically and the photon is usually radiated from an initial-state quark. Measurements are made in both one and two observables, including those sensitive to the hard scattering in the event and others which probe additional soft and collinear radiation. Different Standard Model predictions, from both parton-shower Monte Carlo simulation and fixed-order QCD calculations, are compared with the measurements. In general, good agreement is observed between data and predictions from MATRIX and MiNNLO$_\text{PS}$, as well as next-to-leading-order predictions from MadGraph5_aMC@NLO and Sherpa.

The production of $\phi$ mesons has been studied in pp collisions at LHC energies with the ALICE detector via the dimuon decay channel in the rapidity region $2.5 < y < 4$. Measurements of the differential cross section ${\rm d}^2\sigma/{\rm d}y {\rm d}p_{\rm T}$ are presented as a function of the transverse momentum ($p_{\rm T}$) at the center-of-mass energies $\sqrt{s}=5.02$, 8 and 13 TeV and compared with the ALICE results at midrapidity. The differential cross sections at $\sqrt{s}=5.02$ and 13 TeV are also studied in several rapidity intervals as a function of $p_{\rm T}$, and as a function of rapidity in three $p_{\rm T}$ intervals. A hardening of the $p_{\rm T}$-differential cross section with the collision energy is observed, while, for a given energy, $p_{\rm T}$ spectra soften with increasing rapidity and, conversely, rapidity distributions get slightly narrower at increasing $p_{\rm T}$. The new results, complementing the published measurements at $\sqrt{s}=2.76$ and 7 TeV, allow one to establish the energy dependence of $\phi$ meson production and to compare the measured cross sections with phenomenological models. None of the considered models manages to describe the evolution of the cross section with $p_{\rm T}$ and rapidity at all the energies.

The multiplicity dependence of jet production in pp collisions at the centre-of-mass energy of $\sqrt{s} = 13\ \mathrm{TeV}$ is studied for the first time. Jets are reconstructed from charged particles using the anti-$k_\mathrm{T}$ algorithm with resolution parameters $R$ varying from $0.2$ to $0.7$. The jets are measured in the pseudorapidity range $|\eta_{\rm jet}|< 0.9-R$ and in the transverse momentum range $5<p_\mathrm{T,jet}^{\rm ch}<140\ \mathrm{GeV}/c$. The multiplicity intervals are categorised by the ALICE forward detector V0. The $p_{\mathrm{T}}$ differential cross section of charged-particle jets are compared to leading order (LO) and next-to-leading order (NLO) perturbative quantum chromodynamics (pQCD) calculations. It is found that the data are better described by the NLO calculation, although the NLO prediction overestimates the jet cross section below $20\ \mathrm{GeV}/c$. The cross section ratios for different $R$ are also measured and compared to model calculations. These measurements provide insights into the angular dependence of jet fragmentation. The jet yield increases with increasing self-normalised charged-particle multiplicity. This increase shows only a weak dependence on jet transverse momentum and resolution parameter at the highest multiplicity. While such behaviour is qualitatively described by the present version of PYTHIA, quantitative description may require implementing new mechanisms for multi-particle production in hadronic collisions.

A search is presented for long-lived particles produced in pairs in proton-proton collisions at the LHC operating at a center-of-mass energy of 13 TeV. The data were collected with the CMS detector during the period from 2015 through 2018, and correspond to a total integrated luminosity of 140 fb$^{-1}$. This search targets pairs of long-lived particles with mean proper decay lengths between 0.1 and 100 mm, each of which decays into at least two quarks that hadronize to jets, resulting in a final state with two displaced vertices. No significant excess of events with two displaced vertices is observed. In the context of $R$-parity violating supersymmetry models, the pair production of long-lived neutralinos, gluinos, and top squarks is excluded at 95% confidence level for cross sections larger than 0.08 fb, masses between 800 and 3000 GeV, and mean proper decay lengths between 1 and 25 mm.

The yields and average transverse momenta of pions, kaons, and antiprotons produced at the Fermilab p¯p collider at s=300, 540, 1000, and 1800 GeV are presented and compared with data from the energies reached at the CERN collider. We also present data on the dependence of average transverse momentum 〈pt〉 and particle ratios as a function of charged particle density dNcdη; data for particle densities as high as six times the average value, corresponding to a Bjorken energy density 6 GeV/fm3, are reported. These data are relevant to the search for quark-gluon phase of QCD.

Two-particle azimuthal correlations are measured with the ALICE apparatus in pp collisions at $\sqrt{s} = 13$ TeV to explore strangeness- and multiplicity-related effects in the fragmentation of jets and the transition regime between bulk and hard production, probed with the condition that a strange meson (K$^{0}_{\rm S}$) or baryon ($\Lambda$) with transverse momentum $p_{\rm T} > 3$ GeV/c is produced. Azimuthal correlations between kaons or $\Lambda$ hyperons with other hadrons are presented at midrapidity for a broad range of the trigger ($3 < p_{\rm T}^{\rm trigg} < 20$ GeV/$c$) and associated particle $p_{\rm T}$ (1 GeV/$c$$< p_{\rm T}^{\rm assoc} < p_{\rm T}^{\rm trigg}$), for minimum-bias events and as a function of the event multiplicity. The near- and away-side peak yields are compared for the case of either K$^{0}_{\rm S}$ or $\Lambda$($\overline{\Lambda}$) being the trigger particle with that of inclusive hadrons (a sample dominated by pions). In addition, the measurements are compared with predictions from PYTHIA 8 and EPOS LHC event generators.

An inclusive search for long-lived exotic particles decaying to a pair of muons is presented. The search uses data collected by the CMS experiment at the CERN LHC in proton-proton collisions at $\sqrt{s}$ = 13 TeV in 2016 and 2018 and corresponding to an integrated luminosity of 97.6 fb$^{-1}$. The experimental signature is a pair of oppositely charged muons originating from a common secondary vertex spatially separated from the pp interaction point by distances ranging from several hundred $\mu$m to several meters. The results are interpreted in the frameworks of the hidden Abelian Higgs model, in which the Higgs boson decays to a pair of long-lived dark photons Z$_\mathrm{D}$, and of a simplified model, in which long-lived particles are produced in decays of an exotic heavy neutral scalar boson. For the hidden Abelian Higgs model with $m_\mathrm{Z_D}$ greater than 20 GeV and less than half the mass of the Higgs boson, they provide the best limits to date on the branching fraction of the Higgs boson to dark photons for $c\tau$(Z$_\mathrm{D}$) (varying with $m_\mathrm{Z_D}$) between 0.03 and ${\approx}$ 0.5 mm, and above ${\approx}$ 0.5 m. Our results also yield the best constraints on long-lived particles with masses larger than 10 GeV produced in decays of an exotic scalar boson heavier than the Higgs boson and decaying to a pair of muons.

A measurement of single top-quark production in the s-channel is performed in proton$-$proton collisions at a centre-of-mass energy of 13 TeV with the ATLAS detector at the CERN Large Hadron Collider. The dataset corresponds to an integrated luminosity of 139 fb$^{-1}$. The analysis is performed on events with an electron or muon, missing transverse momentum and exactly two $b$-tagged jets in the final state. A discriminant based on matrix element calculations is used to separate single-top-quark s-channel events from the main background contributions, which are top-quark pair production and $W$-boson production in association with jets. The observed (expected) signal significance over the background-only hypothesis is 3.3 (3.9) standard deviations, and the measured cross-section is $\sigma=8.2^{+3.5}_{-2.9}$ pb, consistent with the Standard Model prediction of $\sigma^{\mathrm{SM}}=10.32^{+0.40}_{-0.36}$ pb.

Measurements of differential cross sections are presented for inclusive isolated-photon production in $pp$ collisions at a centre-of-mass energy of 13 TeV provided by the LHC and using 139 fb$^{-1}$ of data recorded by the ATLAS experiment. The cross sections are measured as functions of the photon transverse energy in different regions of photon pseudorapidity. The photons are required to be isolated by means of a fixed-cone method with two different cone radii. The dependence of the inclusive-photon production on the photon isolation is investigated by measuring the fiducial cross sections as functions of the isolation-cone radius and the ratios of the differential cross sections with different radii in different regions of photon pseudorapidity. The results presented in this paper constitute an improvement with respect to those published by ATLAS earlier: the measurements are provided for different isolation radii and with a more granular segmentation in photon pseudorapidity that can be exploited in improving the determination of the proton parton distribution functions. These improvements provide a more in-depth test of the theoretical predictions. Next-to-leading-order QCD predictions from JETPHOX and SHERPA and next-to-next-to-leading-order QCD predictions from NNLOJET are compared to the measurements, using several parameterisations of the proton parton distribution functions. The measured cross sections are well described by the fixed-order QCD predictions within the experimental and theoretical uncertainties in most of the investigated phase-space region.

A search for pair-produced vector-like quarks using events with exactly one lepton ($e$ or $\mu$), at least four jets including at least one $b$-tagged jet, and large missing transverse momentum is presented. Data from proton-proton collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV, recorded by the ATLAS detector at the LHC from 2015 to 2018 and corresponding to an integrated luminosity of 139 fb$^{-1}$, are analysed. Vector-like partners $T$ and $B$ of the top and bottom quarks are considered, as is a vector-like $X$ with charge +5/3, assuming their decay into a $W$, $Z$, or Higgs boson and a third-generation quark. No significant deviations from the Standard Model expectation are observed. Upper limits on the production cross-section of $T$ and $B$ quark pairs as a function of their mass are derived for various decay branching ratio scenarios. The strongest lower limits on the masses are 1.59 TeV assuming mass-degenerate VLQs and branching ratios corresponding to the weak-isospin doublet model, and 1.47 TeV (1.46 TeV) for exclusive $T \rightarrow Zt$ ($B/X \rightarrow Wt$) decays. In addition, lower limits on the $T$ and $B$ quark masses are derived for all possible branching ratios.

In a special run of the LHC with $\beta^\star = 2.5~$km, proton-proton elastic-scattering events were recorded at $\sqrt{s} = 13~$TeV with an integrated luminosity of $340~\mu \textrm{b}^{-1}$ using the ALFA subdetector of ATLAS in 2016. The elastic cross section was measured differentially in the Mandelstam $t$ variable in the range from $-t = 2.5 \cdot 10^{-4}~$GeV$^{2}$ to $-t = 0.46~$GeV$^{2}$ using 6.9 million elastic-scattering candidates. This paper presents measurements of the total cross section $\sigma_{\textrm{tot}}$, parameters of the nuclear slope, and the $\rho$-parameter defined as the ratio of the real part to the imaginary part of the elastic-scattering amplitude in the limit $t \rightarrow 0$. These parameters are determined from a fit to the differential elastic cross section using the optical theorem and different parameterizations of the $t$-dependence. The results for $\sigma_{\textrm{tot}}$ and $\rho$ are \begin{equation*} \sigma_{\textrm{tot}}(pp\rightarrow X) = \mbox{104.7} \pm 1.1 \; \mbox{mb} , \; \; \; \rho = \mbox{0.098} \pm 0.011 . \end{equation*} The uncertainty in $\sigma_{\textrm{tot}}$ is dominated by the luminosity measurement, and in $\rho$ by imperfect knowledge of the detector alignment and by modelling of the nuclear amplitude.

A search for physics beyond the standard model (SM) in the final state with a hadronically decaying tau lepton and a neutrino is presented. This analysis is based on data recorded by the CMS experiment from proton-proton collisions at a center-of-mass energy of 13 TeV at the LHC, corresponding to a total integrated luminosity of 138 fb$^{=1}$. The transverse mass spectrum is analyzed for the presence of new physics. No significant deviation from the SM prediction is observed. Limits are set on the production cross section of a W' boson decaying into a tau lepton and a neutrino. Lower limits are set on the mass of the sequential SM-like heavy charged vector boson and the mass of a quantum black hole. Upper limits are placed on the couplings of a new boson to the SM fermions. Constraints are put on a nonuniversal gauge interaction model and an effective field theory model. For the first time, upper limits on the cross section of $t$-channel leptoquark (LQ) exchange are presented. These limits are translated into exclusion limits on the LQ mass and on its coupling in the $t$-channel. The sensitivity of this analysis extends into the parameter space of LQ models that attempt to explain the anomalies observed in B meson decays. The limits presented for the various interpretations are the most stringent to date. Additionally, a model-independent limit is provided.

Production of inclusive charmonia in pp collisions at center-of-mass energy of $\sqrt{s}$ = 13 TeV and p-Pb collisions at center-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}}$ = 8.16 TeV is studied as a function of charged-particle pseudorapidity density with ALICE. Ground and excited charmonium states (J/$\psi$, $\psi$(2S)) are measured from their dimuon decays in the interval of rapidity in the center-of-mass frame $2.5 < y_{\rm cms} < 4.0$ for pp collisions, and $2.03 < y_{\rm cms} < 3.53$ and $-4.46 < y_{\rm cms} < -2.96$ for p-Pb collisions. The charged-particle pseudorapidity density is measured around midrapidity ($|\eta|<1.0$). In pp collisions, the measured charged-particle multiplicity extends to about six times the average value, while in p-Pb collisions at forward (backward) rapidity a multiplicity corresponding to about three (four) times the average is reached. The $\psi$(2S) yield increases with the charged-particle pseudorapidity density. The ratio of $\psi$(2S) over J/$\psi$ yield does not show a significant multiplicity dependence in either colliding system, suggesting a similar behavior of J/$\psi$ and $\psi$(2S) yields with respect to charged-particle pseudorapidity density. Results for the $\psi$(2S) yield and its ratio with respect to J/$\psi$ agree with available model calculations.

The multiplicity dependence of the pseudorapidity density of charged particles in proton-proton (pp) collisions at centre-of-mass energies $\sqrt{s}$ = 5.02, 7 and 13 TeV measured by ALICE is reported. The analysis relies on track segments measured in the midrapidity range ($|\eta| < 1.5$). Results are presented for inelastic events having at least one charged particle produced in the pseudorapidity interval $|\eta|<1$ ($\mathrm{INEL}_{>0}$). The multiplicity dependence of the pseudorapidy density of charged particles is measured with mid and forward rapidity multiplicity estimators, the latter being less affected by autocorrelations. A detailed comparison with predictions from the PYTHIA 8 and EPOS LHC event generators is also presented. Both generators provide a good description of the data.

A search for pairs of dijet resonances with the same mass is conducted in final states with at least four jets. Results are presented separately for the case where the four jet production proceeds via an intermediate resonant state and for nonresonant production. The search uses a data sample corresponding to an integrated luminosity of 138 fb$^{-1}$ collected by the CMS detector in proton-proton collisions at $\sqrt{s}$ = 13 TeV. Model-independent limits, at 95% confidence level, are reported on the production cross section of four-jet and dijet resonances. These first LHC limits on resonant pair production of dijet resonances via high mass intermediate states are applied to a signal model of diquarks that decay into pairs of vector-like quarks, excluding diquark masses below 7.6 TeV for a particular model scenario. There are two events in the tails of the distributions, each with a four-jet mass of 8 TeV and an average dijet mass of 2 TeV, resulting in local and global significances of 3.9 and 1.6 standard deviations, respectively, if interpreted as a signal. The nonresonant search excludes pair production of top squarks with masses between 0.50 TeV to 0.77 TeV, with the exception of a small interval between 0.52 and 0.58 TeV, for supersymmetric $R$-parity-violating decays to quark pairs, significantly extending previous limits. Here, the most significant excess above the predicted background occurs at an average dijet mass of 0.95 TeV, for which the local and global significances are 3.6 and 2.5 standard deviations, respectively.