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.

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.

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.

The collective behavior of K$^0_\mathrm{S}$ and $\Lambda/\bar{\Lambda}$ strange hadrons is studied by measuring the elliptic azimuthal anisotropy ($v_2$) using the scalar-product and multiparticle correlation methods. Proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV and lead-lead (PbPb) collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV collected by the CMS experiment at the LHC are investigated. Nonflow effects in the pPb collisions are studied by using a subevent cumulant analysis and by excluding events where a jet with transverse momentum greater than 20\GeV is present. The strange hadron $v_2$ values extracted in \pPb collisions via the four- and six-particle correlation method are found to be nearly identical, suggesting the collective behavior. Comparisons of the pPb and PbPb results for both strange hadrons and charged particles illustrate how event-by-event flow fluctuations depend on the system size.

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.

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.

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 first observation of top quark production in proton-nucleus collisions is reported using proton-lead data collected by the CMS experiment at the CERN LHC at a nucleon-nucleon center-of-mass energy of sqrt(s[NN]) = 8.16 TeV. The measurement is performed using events with exactly one isolated electron or muon and at least four jets. The data sample corresponds to an integrated luminosity of 174 inverse nanobarns. The significance of the tt-bar signal against the background-only hypothesis is above five standard deviations. The measured cross section is sigma[tt-bar] = 45 +/- 8 nb, consistent with predictions from perturbative quantum chromodynamics.

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.

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.

The central exclusive production of charged-hadron pairs in pp collisions at a centre-of-mass energy of 13\TeV is examined, based on data collected in a special high-$\beta^*$ run of the LHC. The nonresonant continuum processes are studied with the invariant mass of the centrally produced two-pion system in the resonance-free region, $m_{\pi^+\pi^-}$$\lt$ 0.7 GeV or $m_{\pi^+\pi^-}$$\gt$ 1.8 GeV. Differential cross sections as functions of the azimuthal angle between the surviving protons, squared exchanged four-momenta, and $m_{\pi^+\pi^-}$ are measured in a wide region of scattered proton transverse momenta, between 0.2 and 0.8 GeV, and for pion rapidities $\lvert y\rvert$$\lt$ 2. A rich structure of interactions related to double-pomeron exchange is observed. A parabolic minimum in the distribution of the two-proton azimuthal angle is observed for the first time. It can be interpreted as an effect of additional pomeron exchanges between the protons from the interference between the bare and the rescattered amplitudes. After model tuning, various physical quantities are determined that are related to the pomeron cross section, proton-pomeron and meson-pomeron form factors, pomeron trajectory and intercept, and coefficients of diffractive eigenstates of the proton.

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.

The inclusive jet cross section is measured as a function of jet transverse momentum $p_\mathrm{T}$ and rapidity $y$. The measurement is performed using proton-proton collision data at $\sqrt{s}$ = 5.02 TeV, recorded by the CMS experiment at the LHC, corresponding to an integrated luminosity of 27.4 pb$^{-1}$. The jets are reconstructed with the anti-$k_\mathrm{T}$ algorithm using a distance parameter of $R$ = 0.4, within the rapidity interval $\lvert y\rvert$$\lt$ 2, and across the kinematic range 0.06 $\lt$$p_\mathrm{T}$$\lt$ 1 TeV. The jet cross section is unfolded from detector to particle level using the determined jet response and resolution. The results are compared to predictions of perturbative quantum chromodynamics, calculated at both next-to-leading order and next-to-next-to-leading order. The predictions are corrected for nonperturbative effects, and presented for a variety of parton distribution functions and choices of the renormalization/factorization scales and the strong coupling $\alpha_\mathrm{S}$.

A search is presented for charged, long-lived supersymmetric particles in final states with one or more disappearing tracks. The search is based on data from proton-proton collisions at a center-of-mass energy of 13 TeV collected with the CMS detector at the CERN LHC between 2016 and 2018, corresponding to an integrated luminosity of 137 fb$^{-1}$. The search is performed over final states characterized by varying numbers of jets, b-tagged jets, electrons, and muons. The length of signal-candidate tracks in the plane perpendicular to the beam axis is used to characterize the lifetimes of wino- and higgsino-like charginos produced in the context of the minimal supersymmetric standard model. The d$E$/d$x$ energy loss of signal-candidate tracks is used to increase the sensitivity to charginos with a large mass and thus a small Lorentz boost. The observed results are found to be statistically consistent with the background-only hypothesis. Limits on the pair production cross section of gluinos and squarks are presented in the framework of simplified models of supersymmetric particle production and decay, and for electroweakino production based on models of wino and higgsino dark matter. The limits presented are the most stringent to date for scenarios with light third-generation squarks and a wino- or higgsino-like dark matter candidate capable of explaining the known dark matter relic density.

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.

A search for partonic collective effects inside jets produced in proton-proton collisions is performed via correlation measurements of charged constituents using the CMS detector at the CERN LHC. The analysis uses data collected at a center-of-mass energy of $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. Jets are reconstructed with the anti-$k_\mathrm{T}$ algorithm with a distance parameter of 0.8 and are required to have transverse momentum greater than 550 GeV and pseudorapidity $\lvert\eta\rvert$$\lt$ 1.6. Two-particle correlations among the charged constituents within the jets are studied as functions of the particles' azimuthal angle and pseudorapidity separations ($\Delta\phi^*$ and $\Delta\eta^*$) in a jet coordinate basis, where constituents' $\eta^*$, $\phi^*$ are defined relative to the direction of the jet. The correlation functions are studied in classes of in-jet charged-particle multiplicity up to $N_\text{ch}^\mathrm{j}$$\approx$ 100. Fourier harmonics are extracted from long-range azimuthal correlation functions to characterize azimuthal anisotropy for $\lvert\Delta\eta^*\rvert$$\gt$ 2. For low-multiplicity jets, the long-range elliptic anisotropic harmonic, $v^*_2$, is observed to decrease with $N_\text{ch}^\mathrm{j}$. This trend is well described by Monte Carlo event generators. However, a rising trend for $v^*_2$ emerges at $N_\text{ch}^\mathrm{j}$$\gtrsim$ 80, hinting at a possible onset of collective behavior, which is not reproduced by the models tested. This observation yields new insights into the dynamics of parton fragmentation processes in the vacuum.

Ultrarelativistic nuclear collisions create a strongly interacting state of hot and dense quark-gluon matter that exhibits a remarkable collective flow behavior with minimal viscous dissipation. To gain deeper insights into its intrinsic nature and fundamental degrees of freedom, we extracted the speed of sound in this medium created using lead-lead (PbPb) collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data were recorded by the CMS experiment at the CERN LHC and correspond to an integrated luminosity of 0.607 nb$^{-1}$. The measurement is performed by studying the multiplicity dependence of the average transverse momentum of charged particles emitted in head-on PbPb collisions. Our findings reveal that the speed of sound in this matter is nearly half the speed of light, with a squared value of 0.241 $\pm$ 0.002 (stat) $\pm$ 0.016 (syst) in natural units. The effective medium temperature, estimated using the mean transverse momentum, is 219 $\pm$ 8 (syst) MeV. The measured squared speed of sound at this temperature aligns precisely with predictions from lattice quantum chromodynamic (QCD) calculations. This result provides a stringent constraint on the equation of state of the created medium and direct evidence for a deconfined QCD phase being attained in relativistic nuclear collisions.

The study of (anti-)deuteron production in pp collisions has proven to be a powerful tool to investigate the formation mechanism of loosely bound states in high energy hadronic collisions. In this paper the production of (anti-)deuterons is studied as a function of the charged particle multiplicity in inelastic pp collisions at $\sqrt{s}=13$ TeV using the ALICE experiment. Thanks to the large number of accumulated minimum bias events, it has been possible to measure (anti-)deuteron production in pp collisions up to the same charged particle multiplicity ($\rm{d} N_{ch}/\rm{d}\eta\sim26$) as measured in p-Pb collisions at similar centre-of-mass energies. Within the uncertainties, the deuteron yield in pp collisions resembles the one in p-Pb interactions, suggesting a common formation mechanism behind the production of light nuclei in hadronic interactions. In this context the measurements are compared with the expectations of coalescence and Statistical Hadronisation Models (SHM).

A search is presented for baryon number violating interactions in top quark production and decay. The analysis uses data from proton-proton collisions at a center-of-mass energy of 13 TeV, collected with the CMS detector at the LHC with an integrated luminosity of 138 fb$^{-1}$. Candidate events are selected by requiring two oppositely-charged leptons (electrons or muons) and exactly one jet identified as originating from a bottom quark. Multivariate discriminants are used to separate the signal from the background. No significant deviation from the standard model prediction is observed. Upper limits are placed on the strength of baryon number violating couplings. For the first time the production of single top quarks via baryon number violating interactions is studied. This allows the search to set the most stringent constraints to date on the branching fraction of the top quark decay to a lepton, an up-type quark (u or c), and a down-type quark (d, s, or b). The results improve the previous bounds by three to six orders of magnitude based on the fermion flavor combination of the baryon number violating interactions.

A search for long-lived particles (LLPs) decaying in the CMS muon detectors is presented. A data sample of proton-proton collisions at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$ recorded at the LHC in 2016-2018, is used. The decays of LLPs are reconstructed as high multiplicity clusters of hits in the muon detectors. In the context of twin Higgs models, the search is sensitive to LLP masses from 0.4 to 55 GeV and a broad range of LLP decay modes, including decays to hadrons, $\tau$ leptons, electrons, or photons. No excess of events above the standard model background is observed. The most stringent limits to date from LHC data are set on the branching fraction of the Higgs boson decay to a pair of LLPs with masses below 10 GeV. This search also provides the best limits for various intervals of LLP proper decay length and mass. Finally, this search sets the first limits at the LHC on a dark quantum chromodynamic sector whose particles couple to the Higgs boson through gluon, Higgs boson, photon, vector, and dark-photon portals, and is sensitive to branching fractions of the Higgs boson to dark quarks as low as 2$\times$10$^{-3}$.

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.

Despite the f$_0$(980) hadron having been discovered half a century ago, the question about its quark content has not been settled: it might be an ordinary quark-antiquark ($\mathrm{q\bar{q}}$) meson, a tetraquark ($\mathrm{q\bar{q}q\bar{q}}$) exotic state, a kaon-antikaon ($\mathrm{K\bar{K}}$) molecule, or a quark-antiquark-gluon ($\mathrm{q\bar{q}g}$) hybrid. This paper reports strong evidence that the f$_0$(980) state is an ordinary $\mathrm{q\bar{q}}$ meson, inferred from the scaling of elliptic anisotropies ($v_2$) with the number of constituent quarks ($n_\mathrm{q}$), as empirically established using conventional hadrons in relativistic heavy ion collisions. The f$_0$(980) state is reconstructed via its dominant decay channel f$_0$(980) $\to$$\pi^+\pi^-$, in proton-lead collisions recorded by the CMS experiment at the LHC, and its $v_2$ is measured as a function of transverse momentum ($p_\mathrm{T}$). It is found that the $n_q$ = 2 ($\mathrm{q\bar{q}}$ state) hypothesis is favored over $n_q$ = 4 ($\mathrm{q\bar{q}q\bar{q}}$ or $\mathrm{K\bar{K}}$ states) by 7.7, 6.3, or 3.1 standard deviations in the $p_\mathrm{T}$$\lt$ 10, 8, or 6 GeV/$c$ ranges, respectively, and over $n_\mathrm{q}$ = 3 ($\mathrm{q\bar{q}g}$ hybrid state) by 3.5 standard deviations in the $p_\mathrm{T}$$\lt$ 8 GeV/$c$ range. This result represents the first determination of the quark content of the f$_0$(980) state, made possible by using a novel approach, and paves the way for similar studies of other exotic hadron candidates.

A combination of the results of several searches for the electroweak production of the supersymmetric partners of standard model bosons, and of charged leptons, is presented. All searches use proton-proton collision data at $\sqrt{s}$ = 13 TeV recorded with the CMS detector at the LHC in 2016-2018. The analyzed data correspond to an integrated luminosity of up to 137 fb$^{-1}$. The results are interpreted in terms of simplified models of supersymmetry. Two new interpretations are added with this combination: a model spectrum with the bino as the lightest supersymmetric particle together with mass-degenerate higgsinos decaying to the bino and a standard model boson, and the compressed-spectrum region of a previously studied model of slepton pair production. Improved analysis techniques are employed to optimize sensitivity for the compressed spectra in the wino and slepton pair production models. The results are consistent with expectations from the standard model. The combination provides a more comprehensive coverage of the model parameter space than the individual searches, extending the exclusion by up to 125 GeV, and also targets some of the intermediate gaps in the mass coverage.

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 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.

A search for the lepton flavor violating $\tau$$\to$ 3$\mu$ decay is performed using proton-proton collision events at a center-of-mass energy of 13 TeV collected by the CMS experiment at the LHC in 2017-2018, corresponding to an integrated luminosity of 97.7 fb$^{-1}$. Tau leptons produced in both heavy-flavor hadron and W boson decays are exploited in the analysis. No evidence for the decay is observed. The results of this search are combined with an earlier null result based on data collected in 2016 to obtain a total integrated luminosity of 131 fb$^{-1}$. The observed (expected) upper limits on the branching fraction $\mathcal{B}$($\tau$$\to$ 3$\mu$) at confidence levels of 90 and 95% are 2.9 $\times$ 10$^{-8}$ (2.4 $\times$ 10$^{-8}$) and 3.6 $\times$ 10$^{-8}$ (3.0 $\times$ 10$^{-8}$), respectively.

Ultrarelativistic heavy ion collisions recreate in the laboratory the thermodynamical conditions prevailing in the early universe up to 10$^{-6}$ seconds, thereby allowing the study of the quark-gluon plasma (QGP), a state of quantum chromodynamics (QCD) matter with deconfined partons. The top quark, the heaviest elementary particle known, is accessible in nucleus-nucleus collisions at the CERN LHC, and constitutes a novel probe of the QGP. Here, we report the first-ever evidence for the production of top quarks in nucleus-nucleus collisions, using lead-lead collision data at a nucleon-nucleon centre-of-mass energy of 5.02 TeV recorded by the CMS experiment. Two methods are used to measure the cross section for top quark pair production ($\sigma_\mathrm{t\bar{t}}$) via the decay into charged leptons (electrons or muons) and bottom quarks. One method relies on the leptonic information alone, and the second one exploits, in addition, the presence of bottom quarks. The measured cross sections, $\sigma_\mathrm{t\bar{t}} = $ 2.54 $^{+0.84}_{-0.74}$ and 2.03 $^{+0.71}_{-0.64}$ $\mu$b, respectively, are compatible with expectations from scaled proton-proton data and QCD predictions.

A search for the production of long-lived particles in proton-proton collisions at a center-of-mass energy of 13 TeV at the CERN LHC is presented. The search is based on data collected by the CMS experiment in 2016-2018, corresponding to a total integrated luminosity of 137 fb$^{-1}$. This search is designed to be sensitive to long-lived particles with mean proper decay lengths between 0.1 and 1000 $\mu$m, whose decay products produce a final state with at least one displaced vertex and missing transverse momentum. A machine learning algorithm, which improves the background rejection power by more than an order of magnitude, is applied to improve the sensitivity. The observation is consistent with the standard model background prediction, and the results are used to constrain split supersymmetry (SUSY) and gauge-mediated SUSY breaking models with different gluino mean proper decay lengths and masses. This search is the first CMS search that shows sensitivity to hadronically decaying long-lived particles from signals with mass differences between the gluino and neutralino below 100 GeV. It sets the most stringent limits to date for split-SUSY models and gauge-mediated SUSY breaking models with gluino proper decay length less than 6 $\mu$m.

A search is reported for near-threshold structures in the J/$\psi$J/$\psi$ invariant mass spectrum produced in proton-proton collisions at $\sqrt{s}$ = 13 TeV from data collected by the CMS experiment, corresponding to an integrated luminosity of 135 fb$^{-1}$. Three structures are found, and a model with quantum interference among these structures provides a good description of the data. A new structure is observed with a significance above 5 standard deviations at a mass of 6638 $^{+43}_{-38}$ (stat) $^{+16}_{-31}$ (syst) MeV. Another structure with even higher significance is found at a mass of 6847 $^{+44}_{-28}$ (stat) $^{+48}_{-20}$ (syst) MeV, which is consistent with the X(6900) resonance reported by the LHCb experiment and confirmed by the ATLAS experiment. Evidence for another new structure, with a local significance of 4.7 standard deviations, is found at a mass of 7134 $^{+48}_{-25}$ (stat) $^{+41}_{-15}$ (syst) MeV. Results are also reported for a model without interference, which does not fit the data as well and shows mass shifts up to 150 MeV relative to the model with interference.

A search for dark matter particles is performed using events with a pair of W bosons and large missing transverse momentum. Candidate events are selected by requiring one or two leptons ($\ell =$ electrons or muons). The analysis is based on proton-proton collision data collected at a center-of-mass energy of 13 TeV by the CMS experiment at the LHC and corresponding to an integrated luminosity of 138 fb$^{-1}$. No significant excess over the expected standard model background is observed in the $\ell\nu$qq and 2$\ell$2$\nu$ final states of the W$^+$W$^-$ boson pair. Limits are set on dark matter production in the context of a simplified dark Higgs model, with a dark Higgs boson mass above the W$^+$W$^-$ mass threshold. The dark matter phase space is probed in the mass range 100-300 GeV, extending the scope of previous searches. Current exclusion limits are improved in the range of dark Higgs masses from 160 to 250 GeV, for a dark matter mass of 200 GeV.

A study of the associated production of a Z boson and a charm quark jet (Z + c), and a comparison to production with a b quark jet (Z + b), in pp collisions at a centre-of-mass energy of 8 TeV are presented. The analysis uses a data sample corresponding to an integrated luminosity of 19.7 fb$^{-1}$, collected with the CMS detector at the CERN LHC. The Z boson candidates are identified through their decays into pairs of electrons or muons. Jets originating from heavy flavour quarks are identified using semileptonic decays of c or b flavoured hadrons and hadronic decays of charm hadrons. The measurements are performed in the kinematic region with two leptons with $p_{\rm T}^{\ell} > $ 20 GeV, ${|\eta^{\ell}|} < $ 2.1, 71 $ < m_{\ell\ell} < $ 111 GeV, and heavy flavour jets with $p_{\rm T}^{{\rm jet}} > $ 25 GeV and ${|\eta^{{\rm jet}}|} < $ 2.5. The Z + c production cross section is measured to be $\sigma({\mathrm{p}}{\mathrm{p}} \rightarrow \mathrm{Z} + \mathrm{c} + \mathrm{X}) {\cal B}(\mathrm{Z} \rightarrow \ell^+\ell^-) = $ 8.8 $ \pm $ 0.5 (stat) $ \pm $ 0.6 (syst) pb. The ratio of the Z + c and Z + b production cross sections is measured to be $\sigma({\mathrm{p}}{\mathrm{p}} \rightarrow \mathrm{Z} + \mathrm{c} + \mathrm{X}) / \sigma({\mathrm{p}}{\mathrm{p}} \rightarrow \mathrm{Z} + \mathrm{b} + \mathrm{X}) = $ 2.0 $ \pm $ 0.2 (stat) $ \pm $ 0.2 (syst). The Z + c production cross section and the cross section ratio are also measured as a function of the transverse momentum of the Z boson and of the heavy flavour jet. The measurements are compared with theoretical predictions.

For the first time a principle-component analysis is used to separate out different orthogonal modes of the two-particle correlation matrix from heavy ion collisions. The analysis uses data from sqrt(s[NN]) = 2.76 TeV PbPb and sqrt(s[NN]) = 5.02 TeV pPb collisions collected by the CMS experiment at the LHC. Two-particle azimuthal correlations have been extensively used to study hydrodynamic flow in heavy ion collisions. Recently it has been shown that the expected factorization of two-particle results into a product of the constituent single-particle anisotropies is broken. The new information provided by these modes may shed light on the breakdown of flow factorization in heavy ion collisions. The first two modes ("leading" and "subleading") of two-particle correlations are presented for elliptical and triangular anisotropies in PbPb and pPb collisions as a function of pt over a wide range of event activity. The leading mode is found to be essentially equivalent to the anisotropy harmonic previously extracted from two-particle correlation methods. The subleading mode represents a new experimental observable and is shown to account for a large fraction of the factorization breaking recently observed at high transverse momentum. The principle-component analysis technique has also been applied to multiplicity fluctuations. These also show a subleading mode. The connection of these new results to previous studies of factorization is discussed.

The production of four top quarks ($\mathrm{t\bar{t}t\bar{t}}$) is studied with LHC proton-proton collision data samples collected by the CMS experiment at a center-of-mass energy of 13 TeV, and corresponding to integrated luminosities of up to 138 fb$^{-1}$. Events that have no leptons (all-hadronic), one lepton, or two opposite-sign leptons (where lepton refers only to prompt electrons or prompt muons) are considered. This is the first $\mathrm{t\bar{t}t\bar{t}}$ measurement that includes the all-hadronic final state. The observed significance of the $\mathrm{t\bar{t}t\bar{t}}$ signal in these final states of 3.9 standard deviations (1.5 expected) provides evidence for $\mathrm{t\bar{t}t\bar{t}}$ production, with a measured cross section of 36 $^{+12}_{-11}$ fb. Combined with earlier CMS results in other final states, the signal significance is 4.0 standard deviations (3.2 expected). The combination returns an observed cross section of 17 $\pm$ 4 (stat) $\pm$ 3 (syst) fb, which is consistent with the standard model prediction.

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.

An inclusive search for long-lived exotic particles (LLPs) decaying to final states with a pair of muons is presented. The search uses data corresponding to an integrated luminosity of 36.6 fb$^{-1}$ collected by the CMS experiment from the proton-proton collisions at $\sqrt{s}$ = 13.6 TeV in 2022, the first year of Run 3 of the CERN LHC. The experimental signature is a pair of oppositely charged muons originating from a common vertex spatially separated from the proton-proton interaction point by distances ranging from several hundred $\mu$m to several meters. The sensitivity of the search benefits from new triggers for displaced dimuons developed for Run 3. The results are interpreted in the framework of the hidden Abelian Higgs model, in which the Higgs boson decays to a pair of long-lived dark photons, and of an $R$-parity violating supersymmetry model, in which long-lived neutralinos decay to a pair of muons and a neutrino. The limits set on these models are the most stringent to date in wide regions of lifetimes for LLPs with masses larger than 10 GeV.

A measurement of the decorrelation of azimuthal angles between the two jets with the largest transverse momenta is presented for seven regions of leading jet transverse momentum up to 2.2 TeV. The analysis is based on the proton-proton collision data collected with the CMS experiment at a centre-of-mass energy of 8 TeV corresponding to an integrated luminosity of 19.7 inverse femtobarns. The dijet azimuthal decorrelation is caused by the radiation of additional jets and probes the dynamics of multijet production. The results are compared to fixed-order predictions of perturbative quantum chromodynamics (QCD), and to simulations using Monte Carlo event generators that include parton showers, hadronization, and multiparton interactions. Event generators with only two outgoing high transverse momentum partons fail to describe the measurement, even when supplemented with next-to-leading-order QCD corrections and parton showers. Much better agreement is achieved when at least three outgoing partons are complemented through either next-to-leading-order predictions or parton showers. This observation emphasizes the need to improve predictions for multijet production.

A search for dark matter in events with a displaced nonresonant muon pair and missing transverse momentum is presented. The analysis is performed using an integrated luminosity of 138 fb$^{-1}$ of proton-proton (pp) collision data at a center-of-mass energy of 13 TeV produced by the LHC in 2016-2018. No significant excess over the predicted backgrounds is observed. Upper limits are set on the product of the inelastic dark matter production cross section $\sigma$(pp $\to$ A' $\to$$\chi_1$$\chi_2$) and the decay branching fraction $\mathcal{B}$($\chi_2$$\to$$\chi_1 \mu^+ \mu^-$), where A' is a dark photon and $\chi_1$ and $\chi_2$ are states in the dark sector with near mass degeneracy. This is the first dedicated collider search for inelastic dark matter.

A study of the anomalous couplings of the Higgs boson to vector bosons, including $CP$-violation effects, has been conducted using its production and decay in the WW channel. This analysis is performed on proton-proton collision data collected with the CMS detector at the CERN LHC during 2016-2018 at a center-of-mass energy of 13 TeV, and corresponds to an integrated luminosity of 138 fb$^{-1}$. The different-flavor dilepton (e$\mu$) final state is analyzed, with dedicated categories targeting gluon fusion, electroweak vector boson fusion, and associated production with a W or Z boson. Kinematic information from associated jets is combined using matrix element techniques to increase the sensitivity to anomalous effects at the production vertex. A simultaneous measurement of four Higgs boson couplings to electroweak vector bosons is performed in the framework of a standard model effective field theory. All measurements are consistent with the expectations for the standard model Higgs boson and constraints are set on the fractional contribution of the anomalous couplings to the Higgs boson production cross section.

A search for pair production of scalar and vector leptoquarks (LQs) each decaying to a muon and a bottom quark is performed using proton-proton collision data collected at $\sqrt{s}$ = 13 TeV with the CMS detector at the CERN LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. No excess above standard model expectation is observed. Scalar (vector) LQs with masses less than 1810 (2120) GeV are excluded at 95% confidence level, assuming a 100% branching fraction of the LQ decaying to a muon and a bottom quark. These limits represent the most stringent to date.

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 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.

Jet charge is an estimator of the electric charge of a quark, antiquark, or gluon initiating a jet. It is based on the momentum-weighted sum of the electric charges of the jet constituents. Measurements of three charge observables of the leading jet in transverse momentum pT are performed with dijet events. The analysis is carried out with data collected by the CMS experiment at the CERN LHC in proton-proton collisions at sqrt(s) = 8 TeV corresponding to an integrated luminosity of 19.7 inverse femtobarns. The results are presented as a function of the pT of the leading jet and compared to predictions from leading- and next-to-leading-order event generators combined with parton showers. Measured jet charge distributions, unfolded for detector effects, are reported, which expand on previous measurements of the jet charge average and standard deviation in pp collisions.

The production of Z bosons in pPb collisions at sqrt(s[NN]) = 5.02 TeV is studied by the CMS experiment via the electron and muon decay channels. The inclusive cross section is compared to pp collision predictions, and found to scale with the number of elementary nucleon-nucleon collisions. The differential cross sections as a function of the Z boson rapidity and transverse momentum are measured. Though they are found to be consistent within uncertainty with theoretical predictions both with and without nuclear effects, the forward-backward asymmetry suggests the presence of nuclear effects at large rapidities. These results provide new data for constraining nuclear parton distribution functions.

A search for long-lived heavy neutrinos (N) in the decays of \PB mesons produced in proton-proton collisions at $\sqrt{s}$ = 13 TeV is presented. The data sample corresponds to an integrated luminosity of 41.6 fb$^{-1}$ collected in 2018 by the CMS experiment at the CERN LHC, using a dedicated data stream that enhances the number of recorded events containing B mesons. The search probes heavy neutrinos with masses in the range 1 $\lt$$m_\mathrm{N}$$\lt$ 3 GeV and decay lengths in the range 10$^{-2}$$\lt$$c\tau$$\lt$ 10$^{4}$ mm, where $\tau_\mathrm{N}$ is the N proper mean lifetime. Signal events are defined by the signature B $\to$$\ell_\mathrm{B}$NX; N $\to$$\ell^{\pm} \pi^{\mp}$, where the leptons $\ell_\mathrm{B}$ and $\ell$ can be either a muon or an electron, provided that at least one of them is a muon. The hadronic recoil system, X, is treated inclusively and is not reconstructed. No significant excess of events over the standard model background is observed in any of the $\ell^{\pm}\pi^{\mp}$ invariant mass distributions. Limits at 95% confidence level on the sum of the squares of the mixing amplitudes between heavy and light neutrinos, $\vert V_\mathrm{N}\vert^2$, and on $c\tau$ are obtained in different mixing scenarios for both Majorana and Dirac-like N particles. The most stringent upper limit $\vert V_\mathrm{N}\vert^2$ $\lt$ 2.0$\times$10$^{-5}$ is obtained at $m_\mathrm{N}$ = 1.95 GeV for the Majorana case where N mixes exclusively with muon neutrinos. The limits on $\vert V_\mathrm{N}\vert^2$ for masses 1 $\lt$ $m_\mathrm{N}$ $\lt$ 1.7 GeV are the most stringent from a collider experiment to date.

A search is described for the production of a pair of bottom-type vector-like quarks (B VLQs) with mass greater than 1000 GeV. Each B VLQ decays into a b quark and a Higgs boson, a b quark and a Z boson, or a t quark and a W boson. This analysis considers both fully hadronic final states and those containing a charged lepton pair from a Z boson decay. The products of the H $to$ bb boson decay and of the hadronic Z or W boson decays can be resolved as two distinct jets or merged into a single jet, so the final states are classified by the number of reconstructed jets. The analysis uses data corresponding to an integrated luminosity of 138 fb$^{-1}$ collected in proton-proton collisions at $\sqrt{s}$ = 13 TeV with the CMS detector at the LHC from 2016 to 2018. No excess over the expected background is observed. Lower limits are set on the B VLQ mass at 95% confidence level. These depend on the B VLQ branching fractions and are 1570 and 1540 GeV for 100% B $\to$ bH and 100% B $\to$ bZ, respectively. In most cases, the mass limits obtained exceed previous limits by at least 100 GeV.

A search is presented for new Higgs bosons in proton-proton (pp) collision events in which a same-sign top quark pair is produced in association with a jet, via the pp $\to$ tH/A $\to$ t$\mathrm{\bar{t}}$c and pp $\to$ tH/A $\to$ t$\mathrm{\bar{t}}$u processes. Here, H and A represent the extra scalar and pseudoscalar boson, respectively, of the second Higgs doublet in the generalized two-Higgs-doublet model (g2HDM). The search is based on pp collision data collected at a center-of-mass energy of 13 TeV with the CMS detector at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. Final states with a same-sign lepton pair in association with jets and missing transverse momentum are considered. New Higgs bosons in the 200-1000 GeV mass range and new Yukawa couplings between 0.1 and 1.0 are targeted in the search, for scenarios in which either H or A appear alone, or in which they coexist and interfere. No significant excess above the standard model prediction is observed. Exclusion limits are derived in the context of the g2HDM.

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.

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.

A search for ZZ and ZH production in the $\mathrm{b\bar{b}b\bar{b}}$ final state is presented, where H is the standard model (SM) Higgs boson. The search uses an event sample of proton-proton collisions corresponding to an integrated luminosity of 133 fb$^{-1}$ collected at a center-of-mass energy of 13 TeV with the CMS detector at the CERN LHC. The analysis introduces several novel techniques for deriving and validating a multi-dimensional background model based on control samples in data. A multiclass multivariate classifier customized for the $\mathrm{b\bar{b}b\bar{b}}$ final state is developed to derive the background model and extract the signal. The data are found to be consistent, within uncertainties, with the SM predictions. The observed (expected) upper limits at 95% confidence level are found to be 3.8 (3.8) and 5.0 (2.9) times the SM prediction for the ZZ and ZH production cross sections, respectively.

Using 1.8 fb-1 of pp collisions at a center-of-mass energy of 7 TeV recorded by the ATLAS detector at the Large Hadron Collider, we present measurements of the production cross sections of Upsilon(1S,2S,3S) mesons. Upsilon mesons are reconstructed using the di-muon decay mode. Total production cross sections for p_T<70 GeV and in the rapidity interval |Upsilon|<2.25 are measured to be 8.01+-0.02+-0.36+-0.31 nb, 2.05+-0.01+-0.12+-0.08 nb, 0.92+-0.01+-0.07+-0.04 nb respectively, with uncertainties separated into statistical, systematic, and luminosity measurement effects. In addition, differential cross section times di-muon branching fractions for Upsilon(1S), Upsilon(2S), and Upsilon(3S) as a function of Upsilon transverse momentum p_T and rapidity are presented. These cross sections are obtained assuming unpolarized production. If the production polarization is fully transverse or longitudinal with no azimuthal dependence in the helicity frame the cross section may vary by approximately +-20%. If a non-trivial azimuthal dependence is considered, integrated cross sections may be significantly enhanced by a factor of two or more. We compare our results to several theoretical models of Upsilon meson production, finding that none provide an accurate description of our data over the full range of Upsilon transverse momenta accessible with this dataset.