The first search for narrow resonances decaying to three well-separated hadronic jets is presented. The search uses proton-proton collision data corresponding to an integrated luminosity of 138 fb$^{-1}$ at $\sqrt{s}$ = 13 TeV, collected at the CERN LHC. No significant deviations from the background predictions are observed between 1.75-9.00 TeV. The results provide the first mass limits on a right-handed boson Z$_{\mathrm{R}}$ decaying to three gluons, an excited quark decaying via a vector boson to three quarks, as well as updated limits on a Kaluza-Klein gluon decaying via a radion to three gluons.

Using proton-proton collision data corresponding to an integrated luminosity of 140 fb$^{-1}$ collected by the CMS experiment at $\sqrt{s}$ = 13 TeV, the $\Lambda_\text{b}^0$$\to$ J/$\psi\Xi^-$K$^+$ decay is observed for the first time, with a statistical significance exceeding 5 standard deviations. The relative branching fraction, with respect to the $\Lambda_\text{b}^0$$\to$$\psi$(2S)$\Lambda$ decay, is measured to be $\mathcal{B}$($\Lambda_\text{b}^0$$\to$ J/$\psi\Xi^-$K$^+$)/$\mathcal{B}$( $\Lambda_\text{b}^0$$\to$$\psi$(2S)$\Lambda$) = [3.38 $\pm$ 1.02 $\pm$ 0.61 $\pm$ 0.03]%, where the first uncertainty is statistical, the second is systematic, and the third is related to the uncertainties in $\mathcal{B}$($\psi$(2S) $\to$ J/$\psi\pi^+\pi^-$) and $\mathcal{B}$($\Xi^-$ $\to$ $\Lambda\pi^-$).

The first evidence for the standard model production of a top quark in association with a W boson and a Z boson is reported. The measurement is performed in multilepton final states, where the Z boson is reconstructed via its decays to electron or muon pairs and the W boson decays either to leptons or hadrons. The analysed data were recorded by the CMS experiment at the CERN LHC in 2016-2018 in proton-proton collisions at $\sqrt{s}$ = 13 TeV, and correspond to an integrated luminosity of 138 fb$^{-1}$. The measured cross section is 354 $\pm$ 54 (stat) $\pm$ 95 (syst) fb, and corresponds to a statistical significance of 3.4 standard deviations.

A search for exotic decays of the Higgs boson (H) with a mass of 125 GeV to a pair of light pseudoscalars $\mathrm{a}_1$ is performed in final states where one pseudoscalar decays to two b quarks and the other to a pair of muons or $\tau$ leptons. A data sample of proton-proton collisions at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$ recorded with the CMS detector is analyzed. No statistically significant excess is observed over the standard model backgrounds. Upper limits are set at 95% confidence level (CL) on the Higgs boson branching fraction to $\mu\mu$bb and to $\tau\tau$bb, via a pair of $\mathrm{a}_1$s. The limits depend on the pseudoscalar mass $m_{\mathrm{a}_1}$ and are observed to be in the range (0.17-3.3) $\times$ 10$^{-4}$ and (1.7-7.7) $\times$ 10$^{2}$ in the $\mu\mu$bb and $\tau\tau$bb final states, respectively. In the framework of models with two Higgs doublets and a complex scalar singlet (2HDM+S), the results of the two final states are combined to determine model-independent upper limits on the branching fraction $\mathcal{B}$(H $\to$ $\mathrm{a}_1\mathrm{a}_1$ $\to$ $\ell\ell$bb) at 95% CL, with $\ell$ being a muon or a $\tau$ lepton. For different types of 2HDM+S, upper bounds on the branching fraction $\mathcal{B}$(H $\to$ $\mathrm{a}_1\mathrm{a}_1$) are extracted from the combination of the two channels. In most of the Type II 2HDM+S parameter space, $\mathcal{B}($H $\to$ $\mathrm{a}_1\mathrm{a}_1$) values above 0.23 are excluded at 95% CL for $m_{\mathrm{a}_1}$ values between 15 and 60 GeV.

A search for Higgs boson pair (HH) production with one Higgs boson decaying to two bottom quarks and the other to two W bosons are presented. The search is done using proton-proton collisions data at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$ recorded by the CMS detector at the LHC from 2016 to 2018. The final states considered include at least one leptonically decaying W boson. No evidence for the presence of a signal is observed and corresponding upper limits on the HH production cross section are derived. The limit on the inclusive cross section of the nonresonant HH production, assuming that the distributions of kinematic observables are as expected in the standard model (SM), is observed (expected) to be 14 (18) times the value predicted by the SM, at 95% confidence level. The limits on the cross section are also presented as functions of various Higgs boson coupling modifiers, and anomalous Higgs boson coupling scenarios. In addition, limits are set on the resonant HH production via spin-0 and spin-2 resonances within the mass range 250-900 GeV.

A measurement is presented of the primary Lund jet plane (LJP) density in inclusive jet production in proton-proton collisions. The analysis uses 138 fb$^{-1}$ of data collected by the CMS experiment at $\sqrt{s}$ = 13 TeV. The LJP, a representation of the phase space of emissions inside jets, is constructed using iterative jet declustering. The transverse momentum $k_\mathrm{T}$ and the splitting angle $\Delta R$ of an emission relative to its emitter are measured at each step of the jet declustering process. The average density of emissions as function of $\ln(k_\mathrm{T}$/GeV) and $\ln(R/\Delta R)$ is measured for jets with distance parameters $R$ = 0.4 or 0.8, transverse momentum $p_\mathrm{T} \gt$ 700 GeV, and rapidity $\vert y\vert \lt $ 1.7. The jet substructure is measured using the charged-particle tracks of the jet. The measured distributions, unfolded to the level of stable particles, are compared with theoretical predictions from simulations and with perturbative quantum chromodynamics calculations. Due to the ability of the LJP to factorize physical effects, these measurements can be used to improve different aspects of the physics modeling in event generators.

A search is presented for a third-generation leptoquark (LQ) coupled exclusively to a $\tau$ lepton and a b quark. The search is based on proton-proton collision data at a center-of-mass energy of 13 TeV recorded with the CMS detector, corresponding to an integrated luminosity of 138 fb$^{-1}$. Events with $\tau$ leptons and a varying number of jets originating from b quarks are considered, targeting the single and pair production of LQs, as well as nonresonant $t$-channel LQ exchange. An excess is observed in the data with respect to the background expectation in the combined analysis of all search regions. For a benchmark LQ mass of 2 TeV and an LQ-b-$\tau$ coupling strength of 2.5, the excess reaches a local significance of up to 2.8 standard deviations. Upper limits at the 95% confidence level are placed on the LQ production cross section in the LQ mass range 0.5-2.3 TeV, and up to 3 TeV for $t$-channel LQ exchange. Leptoquarks are excluded below masses of 1.22-1.88 TeV for different LQ models and varying coupling strengths up to 2.5. The study of nonresonant $\tau\tau$ production through $t$-channel LQ exchange allows lower limits on the LQ mass of up to 2.3 TeV to be obtained.

A search for Higgs boson pair (HH) production in association with a vector boson V (W or Z boson) is presented. The search is based on proton-proton collision data at a center-of-mass energy of 13 TeV, collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. All hadronic and leptonic decays of V bosons are used. The leptons considered are electrons, muons, and neutrinos. The HH production is searched for in the $\mathrm{b\bar{b}b\bar{b}}$ decay channel. An observed (expected) upper limit at 95% confidence level of VHH production cross section is set at 294 (124) times the standard model prediction. Constraints are also set on the modifiers of the Higgs boson trilinear self-coupling, $\kappa_{\lambda}$, assuming $\kappa_{2\mathrm{V}}$ = 1 and vice versa on the coupling of two Higgs bosons with two vector bosons, $\kappa_{2\mathrm{V}}$. The observed (expected) 95% confidence intervals of these coupling modifiers are -37.7 $\lt$ $\kappa_{\lambda}$ $\lt$ 37.2 (-30.1 $\lt$ $\kappa_{\lambda}$ $\lt$ 8.9) and -12.2 $\lt$ $\kappa_{2\mathrm{V}}$ $\lt$ 13.5 (-7.2 $\lt$ $\kappa_{2\mathrm{V}}$ $\lt$ 8.9), respectively.

Diboson production in association with jets is studied in the fully leptonic final states, pp $\to$ (Z$\gamma^*$)(Z/$\gamma^*$)+jets $\to$ 2$\ell$2$\ell'$+jets, ($\ell,\ell'$ = e or $\mu$) in proton-proton collisions at a center-of-mass energy of 13 TeV. The data sample corresponds to an integrated luminosity of 138 fb$^{-1}$ collected with the CMS detector at the LHC. Differential distributions and normalized differential cross sections are measured as a function of jet multiplicity, transverse momentum $p_\mathrm{T}$, pseudorapidity $\eta$, invariant mass and $\Delta\eta$ of the highest-$p_\mathrm{T}$ and second-highest-$p_\mathrm{T}$ jets, and as a function of invariant mass of the four-lepton system for events with various jet multiplicities. These differential cross sections are compared with theoretical predictions that mostly agree with the experimental data. However, in a few regions we observe discrepancies between the predicted and measured values. Further improvement of the predictions is required to describe the ZZ+jets production in the whole phase space.

The observation of the production of four top quarks in proton-proton collisions is reported, based on a data sample collected by the CMS experiment at a center-of-mass energy of 13 TeV in 2016-2018 at the CERN LHC and corresponding to an integrated luminosity of 138 fb$^{-1}$. Events with two same-sign, three, or four charged leptons (electrons and muons) and additional jets are analyzed. Compared to previous results in these channels, updated identification techniques for charged leptons and jets originating from the hadronization of b quarks, as well as a revised multivariate analysis strategy to distinguish the signal process from the main backgrounds, lead to an improved expected signal significance of 4.9 standard deviations above the background-only hypothesis. Four top quark production is observed with a significance of 5.6 standard deviations, and its cross section is measured to be 17.7 $^{+3.7}_{-3.5}$ (stat) $^{+2.3}_{-1.9}$ (syst) fb, in agreement with the available standard model predictions.

A search for W' bosons decaying to a top and a bottom quark in final states including an electron or a muon is performed with the CMS detector at the LHC. The analyzed data correspond to an integrated luminosity of 138 fb$^{-1}$ of proton-proton collisions at a center-of-mass energy of 13 Tev. Good agreement with the standard model expectation is observed and no evidence for the existence of the W' boson is found over the mass range examined. The largest observed deviation from the standard model expectation is found for a W' boson mass ($m_\mathrm{W'}$) hypothesis of 3.8 TeV with a relative decay width of 1%, with a local (global) significance of 2.6 (2.0) standard deviations. Upper limits on the production cross sections of W' bosons decaying to a top and a bottom quark are set. Left- and right-handed W' bosons with $m_\mathrm{W'}$ below 3.9 and 4.3 TeV, respectively, are excluded at the 95% confidence level, under the assumption that the new particle has a narrow decay width. Limits are also set for relative decay widths up to 30%. These are the most stringent limits to date on this W' boson decay channel.

A search is performed for charged-lepton flavor violating processes in top quark (t) production and decay. The data were collected by the CMS experiment from proton-proton collisions at a center-of-mass energy of 13 TeV and correspond to an integrated luminosity of 138 fb$^{-1}$. The selected events are required to contain one opposite-sign electron-muon pair, a third charged lepton (electron or muon), and at least one jet of which no more than one is associated with a bottom quark. Boosted decision trees are used to distinguish signal from background, exploiting differences in the kinematics of the final states particles. The data are consistent with the standard model expectation. Upper limits at 95% confidence level are placed in the context of effective field theory on the Wilson coefficients, which range between 0.024-0.424 TeV$^{-2}$ depending on the flavor of the associated light quark and the Lorentz structure of the interaction. These limits are converted to upper limits on branching fractions involving up (charm) quarks, t$\to$e$\mu$u (t$\to$e$\mu$c), of 0.032 (0.498)$\times$10$^{-6}$, 0.022 (0.369)$\times$10$^{-6}$, and 0.012 (0.216)$\times$10$^{-6}$ for tensor-like, vector-like, and scalar-like interactions, respectively.

Measurements of inclusive and normalized differential cross sections of the associated production of top quark-antiquark and bottom quark-antiquark pairs, ttbb, are presented. The results are based on data from proton-proton collisions collected by the CMS detector at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The cross sections are measured in the lepton+jets decay channel of the top quark pair, using events containing exactly one isolated electron or muon and at least five jets. Measurements are made in four fiducial phase space regions, targeting different aspects of the ttbb process. Distributions are unfolded to the particle level through maximum likelihood fits, and compared with predictions from several event generators. The inclusive cross section measurements of this process in the fiducial phase space regions are the most precise to date. In most cases, the measured inclusive cross sections exceed the predictions with the chosen generator settings. The only exception is when using a particular choice of dynamic renormalization scale, $\mu_\mathrm{R}=\frac{1}{2} \prod_{i=\mathrm{t, \bar{t}, b, \bar{b}}} m_{\mathrm{T},i}^{1/4}$, where $m_{\mathrm{T},i}^2=m_i^2+p^2_{\mathrm{T},i}$ are the transverse masses of top and bottom quarks. The differential cross sections show varying degrees of compatibility with the theoretical predictions, and none of the tested generators with the chosen settings simultaneously describe all the measured distributions.

A search is presented for fractionally charged particles with charge below 1$e$, using their small energy loss in the tracking detector as a key variable to observe a signal. The analyzed data set corresponds to an integrated luminosity of 138 fb$^{-1}$ of proton-proton collisions collected at $\sqrt{s}$ = 13 TeV in 2016-2018 at the CERN LHC. This is the first search at the LHC for new particles with charges between $e/$3 and $e$. Masses up to 640 GeV and charges as low as $e/$3 are excluded at 95% confidence level. These are the most stringent limits to date for the considered Drell-Yan-like production mode.

A search for an exotic decay of the Higgs boson to a Z boson and a light pseudoscalar particle (a), decaying to a pair of leptons and a pair of photons, respectively, is presented. The search is based on proton-proton collision data at a center-of-mass energy of $\sqrt{s}$ = 13 TeV, collected with the CMS detector and corresponding to an integrated luminosity of 138 fb$^{-1}$. The analysis probes pseudoscalar masses $m_\mathrm{a}$ between 1 and 30 GeV, leading to two pairs of well-isolated leptons and photons. Upper limits at 95% confidence level are set on the Higgs boson production cross section times its branching fraction to two leptons and two photons. The observed (expected) limits are in the range of 1.1-17.8 (1.7-17.9) fb within the probed $m_\mathrm{a}$ interval. An excess of data above the expected standard model background with a local (global) significance of 2.6 (1.3) standard deviations is observed for a mass hypothesis of $m_\mathrm{a}$ = 3 GeV. Limits on models involving axion-like particles, formulated as an effective field theory, are also reported.

A search for the lepton-flavor violating decay of the Higgs boson and potential additional Higgs bosons with a mass in the range 110-160 GeV to an e$^{\pm}\mu^{\mp}$ pair is presented. The search is performed with a proton-proton collision dataset at a center-of-mass energy of 13 TeV collected by the CMS experiment at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. No excess is observed for the Higgs boson. The observed (expected) upper limit on the e$^{\pm}\mu^{\mp}$ branching fraction for it is determined to be 4.4 (4.7) $\times$ 10$^{-5}$ at 95% confidence level, the most stringent limit set thus far from direct searches. The largest excess of events over the expected background in the full mass range of the search is observed at an e$^{\pm}\mu^{\mp}$ invariant mass of approximately 146 GeV with a local (global) significance of 3.8 (2.8) standard deviations.

A measurement of the Higgs boson (H) production via vector boson fusion (VBF) and its decay into a bottom quark-antiquark pair ($\mathrm{b\bar{b}}$) is presented using proton-proton collision data recorded by the CMS experiment at $\sqrt{s}$ = 13 TeV and corresponding to an integrated luminosity of 90.8 fb$^{-1}$. Treating the gluon-gluon fusion process as a background and constraining its rate to the value expected in the standard model (SM) within uncertainties, the signal strength of the VBF process, defined as the ratio of the observed signal rate to that predicted by the SM, is measured to be $\mu^\text{qqH}_\mathrm{Hb\bar{b}}$ = 1.01 $^{+0.55}_{-0.46}$. The VBF signal is observed with a significance of 2.4 standard deviations relative to the background prediction, while the expected significance is 2.7 standard deviations. Considering inclusive Higgs boson production and decay into bottom quarks, the signal strength is measured to be $\mu^\text{incl.}_\mathrm{Hb\bar{b}}$ = 0.99 $^{+0.48}_{-0.41}$, corresponding to an observed (expected) significance of 2.6 (2.9) standard deviations.

A measurement of the dijet production cross section is reported based on proton-proton collision data collected in 2016 at $\sqrt{s}$ = 13 TeV by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of up to 36.3 fb$^{-1}$. Jets are reconstructed with the anti-$k_\mathrm{T}$ algorithm for distance parameters of $R$ = 0.4 and 0.8. Cross sections are measured double-differentially (2D) as a function of the largest absolute rapidity $\lvert y_\text{max}\rvert$ of the two jets with the highest transverse momenta $p_\mathrm{T}$ and their invariant mass $m_{1,2}$, and triple-differentially (3D) as a function of the rapidity separation $y^*$, the total boost $y_\mathrm{b}$, and either $m_{1,2}$ or the average $p_\mathrm{T}$ of the two jets. The cross sections are unfolded to correct for detector effects and are compared with fixed-order calculations derived at next-to-next-to-leading order in perturbative quantum chromodynamics. The impact of the measurements on the parton distribution functions and the strong coupling constant at the mass of the Z boson is investigated, yielding a value of $\alpha_\mathrm{S}$ = 0.1179 $\pm$ 0.0019.

Measurements of the inclusive and differential fiducial cross sections for the Higgs boson production in the H → ZZ → 4ℓ (ℓ = e, μ) decay channel are presented. The results are obtained from the analysis of proton-proton collision data recorded by the CMS experiment at the CERN LHC at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{−1}$. The measured inclusive fiducial cross section is 2.73 ± 0.26 fb, in agreement with the standard model expectation of 2.86 ± 0.1 fb. Differential cross sections are measured as a function of several kinematic observables sensitive to the Higgs boson production and decay to four leptons. A set of double-differential measurements is also performed, yielding a comprehensive characterization of the four leptons final state. Constraints on the Higgs boson trilinear coupling and on the bottom and charm quark coupling modifiers are derived from its transverse momentum distribution. All results are consistent with theoretical predictions from the standard model.

The observation of WW$\gamma$ production in proton-proton collisions at a center-of-mass energy of 13 TeV with an integrated luminosity of 138 fb$^{-1}$ is presented. The observed (expected) significance is 5.6 (5.1) standard deviations. Events are selected by requiring exactly two leptons (one electron and one muon) of opposite charge, moderate missing transverse momentum, and a photon. The measured fiducial cross section for WW$\gamma$ is 5.9 $\pm$ 0.8 (stat) $\pm$ 0.8 (syst) $\pm$ 0.7 (modeling) fb, in agreement with the next-to-leading order quantum chromodynamics prediction. The analysis is extended with a search for the associated production of the Higgs boson and a photon, which is generated by a coupling of the Higgs boson to light quarks. The result is used to constrain the Higgs boson couplings to light quarks.

A search for direct production of low-mass dimuon resonances is performed using $\sqrt{s}$ = 13 TeV proton-proton collision data collected by the CMS experiment during the 2017-2018 operation of the CERN LHC with an integrated luminosity of 96.6 fb$^{-1}$. The search exploits a dedicated high-rate trigger stream that records events with two muons with transverse momenta as low as 3 GeV but does not include the full event information. The search is performed by looking for narrow peaks in the dimuon mass spectrum in the ranges of 1.1-2.6 GeV and 4.2-7.9 GeV. No significant excess of events above the expectation from the standard model background is observed. Model-independent limits on production rates of dimuon resonances within the experimental fiducial acceptance are set. Competitive or world's best limits are set at 90% confidence level for a minimal dark photon model and for a scenario with two Higgs doublets and an extra complex scalar singlet (2HDM+S). Values of the squared kinetic mixing coefficient $\varepsilon^2$ in the dark photon model above 10$^{-6}$ are excluded over most of the mass range of the search. In the 2HDM+S, values of the mixing angle $\sin(\theta_\text{H})$ above 0.08 are excluded over most of the mass range of the search with a fixed ratio of the Higgs doublets vacuum expectation $\tan\beta$ = 0.5.

The polarization of $\tau$ leptons is measured using leptonic and hadronic $\tau$ lepton decays in Z $\to$$\tau^+\tau^-$ events in proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded by CMS at the CERN LHC with an integrated luminosity of 36.3 fb$^{-1}$. The measured $\tau^-$ polarization at the Z boson mass pole is $\mathcal{P}_{\tau}$(Z) = $-$0.144 $\pm$ 0.006 (stat) $\pm$ 0.014 (syst) = $-$0.144 $\pm$ 0.015, in good agreement with the measurement of the $\tau$ lepton asymmetry parameter of $A_{\tau}$ = 0.1439 $\pm$ 0.0043 = $-\mathcal{P}_{\tau}$(Z) at LEP. The $\tau$ polarization depends on the ratio of the vector to axial-vector couplings of the $\tau$ leptons in the neutral current expression, and thus on the effective weak mixing angle $\sin^{2}\theta_\mathrm{W}^{\text{eff}}$, independently of the Z boson production mechanism. The obtained value $\sin^{2}\theta_\mathrm{W}^{\text{eff}}$ = 0.2319 $\pm$ 0.0008 (stat) $\pm$ 0.0018 (syst) = 0.2319 $\pm$ 0.0019 is in good agreement with measurements at e$^+$e$^-$ colliders.

The first search for scalar leptoquarks produced in $\tau$-lepton-quark collisions is presented. It is based on a set of proton-proton collision data recorded with the CMS detector at the LHC at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$. The reconstructed final state consists of a jet, significant missing transverse momentum, and a $\tau$ lepton reconstructed through its hadronic or leptonic decays. Limits are set on the product of the leptoquark production cross section and branching fraction and interpreted as exclusions in the plane of the leptoquark mass and the leptoquark-$\tau$-quark coupling strength.

A search for high-mass dimuon resonance production in association with one or more b quark jets is presented. The study uses proton-proton collision data collected with the CMS detector at the LHC corresponding to an integrated luminosity of 138 fb$^{-1}$ at a center-of-mass energy of 13 TeV. Model-independent limits are derived on the number of signal events with exactly one or more than one b quark jet. Results are also interpreted in a lepton-flavor-universal model with Z$'$ boson couplings to a bb quark pair ($g_\mathrm{b}$), an sb quark pair ($g_\mathrm{b}\delta_\mathrm{bs}$), and any same-flavor charged lepton ($g_\ell$) or neutrino pair ($g_\nu$), with $\left|g_{\nu}\right| = \left|g_\ell\right|$. For a Z$'$ boson with a mass $m_{\mathrm{Z}'}$ = 350 GeV (2 TeV) and $\left|\delta_\mathrm{bs}\right|$$\lt$ 0.25, the majority of the parameter space with 0.0057 $\lt$$\left|g_\ell\right|$$\lt$ 0.35 (0.25 $\lt$$\left|g_\ell\right|$$\lt$ 0.43) and 0.0079 $\lt$$\left|g_\mathrm{b}\right|$$\lt$ 0.46 (0.34 $\lt$$\left|g_\mathrm{b}\right|$$\lt$ 0.57) is excluded at 95% confidence level. Finally, constraints are set on a Z$'$ model with parameters consistent with low-energy b $\to$ s$\ell\ell$ measurements. In this scenario, most of the allowed parameter space is excluded for a Z$'$ boson with 350 $\lt m_{\mathrm{Z}'}$ $\lt$ 500 GeV, while the constraints are less stringent for higher $m_{\mathrm{Z}'}$ hypotheses. This is the first dedicated search at the LHC for a high-mass dimuon resonance produced in association with multiple b quark jets, and the constraints obtained on models with this signature are the most stringent to date.

A search for the rare $\eta$$\to$$\mu^+\mu^-\mu^+\mu^-$ double-Dalitz decay is performed using a sample of proton-proton collisions, collected by the CMS experiment at the CERN LHC with high-rate muon triggers in 2017-2018 and corresponding to an integrated luminosity of 101 fb$^{-1}$. A signal having a statistical significance well in excess of 5 standard deviations is observed. Using the $\eta$$\to$$\mu^+ \mu^-$ decay as normalization, the branching fraction $\mathcal{B}(\eta$$\to$$\mu^+\mu^-\mu^+\mu^-)$ = [5.0 $\pm$ 0.8 (stat) $\pm$ 0.7 (syst) $\pm$ 0.7 ($\mathcal{B}_{2\mu}$)] $\times$ 10$^{-9}$ is measured, where the last term is the uncertainty in the normalization channel branching fraction. This work achieves an improved precision of over five orders of magnitude compared to previous results, leading to the first measurement of this branching fraction, which is found to agree with theoretical predictions.