A search for new heavy particles manifested as resonances in two-jet final states is presented. The data were produced in 7 TeV proton-proton collisions by the Large Hadron Collider (LHC) and correspond to an integrated luminosity of 315 nb^-1 collected by the ATLAS detector. No resonances were observed. Upper limits were set on the product of cross section and signal acceptance for excited-quark (q*) production as a function of q* mass. These exclude at the 95% CL the q* mass interval 0.30 < mq* < 1.26 TeV, extending the reach of previous experiments.

The first measurements from proton-proton collisions recorded with the ATLAS detector at the LHC are presented. Data were collected in December 2009 using a minimum-bias trigger during collisions at a centre-of-mass energy of 900 GeV. The charged-particle multiplicity, its dependence on transverse momentum and pseudorapidity, and the relationship between mean transverse momentum and charged-particle multiplicity are measured for events with at least one charged particle in the kinematic range |eta|<2.5 and pT>500 MeV. The measurements are compared to Monte Carlo models of proton-proton collisions and to results from other experiments at the same centre-of-mass energy. The charged-particle multiplicity per event and unit of pseudorapidity at eta = 0 is measured to be 1.333 +/- 0.003 (stat.) +/- 0.040 (syst.), which is 5-15% higher than the Monte Carlo models predict.

Measurements are presented from proton-proton collisions at centre-of-mass energies of sqrt(s) = 0.9, 2.36 and 7 TeV recorded with the ATLAS detector at the LHC. Events were collected using a single-arm minimum-bias trigger. The charged-particle multiplicity, its dependence on transverse momentum and pseudorapidity and the relationship between the mean transverse momentum and charged-particle multiplicity are measured. Measurements in different regions of phase-space are shown, providing diffraction-reduced measurements as well as more inclusive ones. The observed distributions are corrected to well-defined phase-space regions, using model-independent corrections. The results are compared to each other and to various Monte Carlo models, including a new AMBT1 PYTHIA 6 tune. In all the kinematic regions considered, the particle multiplicities are higher than predicted by the Monte Carlo models. The central charged-particle multiplicity per event and unit of pseudorapidity, for tracks with pT >100 MeV, is measured to be 3.483 +- 0.009 (stat) +- 0.106 (syst) at sqrt(s) = 0.9 TeV and 5.630 +- 0.003 (stat) +- 0.169 (syst) at sqrt(s) = 7 TeV.

A measurement of the cross section for the inclusive production of isolated prompt photons in pp collisions at a centre-of-mass energy sqrt(s) = 7TeV is presented. The measurement covers the pseudorapidity ranges |eta|<1.37 and 1.52<|eta|<1.81 in the transverse energy range 15 < E_T <100 GeV. The results are based on an integrated luminosity of 880 nb-1, collected with the ATLAS detector at the Large Hadron Collider. Photon candidates are identified by combining information from the calorimeters and from the inner tracker. Residual background in the selected sample is estimated from data based on the observed distribution of the transverse isolation energy in a narrow cone around the photon candidate. The results are compared to predictions from next-to-leading order perturbative QCD calculations.

Jet cross sections have been measured for the first time in proton-proton collisions at a centre-of-mass energy of 7 TeV using the ATLAS detector. The measurement uses an integrated luminosity of 17 nb-1 recorded at the Large Hadron Collider. The anti-kt algorithm is used to identify jets, with two jet resolution parameters, R = 0.4 and 0.6. The dominant uncertainty comes from the jet energy scale, which is determined to within 7% for central jets above 60 GeV transverse momentum. Inclusive single-jet differential cross sections are presented as functions of jet transverse momentum and rapidity. Dijet cross sections are presented as functions of dijet mass and the angular variable $\chi$. The results are compared to expectations based on next-to-leading-order QCD, which agree with the data, providing a validation of the theory in a new kinematic regime.

Using a data sample of proton-proton collisions at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 140 fb$^{-1}$ collected by the CMS experiment in 2016-2018, the B$^0_\mathrm{s}$ $\to$ X(3872)$\phi$ decay is observed. Decays into J/$\psi\,\pi^+\pi^-$ and K$^+$K$^-$ are used to reconstruct, respectively, the X(3872) and $\phi$. The ratio of the product of branching fractions $\mathcal{B}($B$^0_\mathrm{s}$ $\to$ X(3872)$\phi)\,\mathcal{B}($X(3872) $\to$ J$/\psi\,\pi^+\pi^-)$ to the product $\mathcal{B}($B$^0_\mathrm{s}$ $\to$ $\psi$(2S)$\phi)\,\mathcal{B}(\psi$(2S) $\to$ J/$\psi\,\pi^+\pi^-)$ is measured to be (2.21 $\pm$ 0.29 (stat) $\pm$ 0.17 (syst))%. The ratio $\mathcal{B}($B$^0_\mathrm{s}$ $\to$ X(3872)$\phi) / \mathcal{B}($B$^{0}$ $\to$ X(3872)K$^0)$ is found to be consistent with one, while the ratio $\mathcal{B}($B$^0_\mathrm{s}$ $\to$ X(3872)$\phi) / \mathcal{B}($B$^+$ $\to$ X(3872)K$^+)$ is two times smaller. This suggests a difference in the production dynamics of the X(3872) in B$^0$ and B$^0_\mathrm{s}$ meson decays compared to B$^+$. The reported observation may shed new light on the nature of the X(3872) particle.

The measurements of the Higgs boson (H) production cross sections performed by the CMS Collaboration in the four-lepton (4$\ell$, $\ell$ = e, $\mu$) final state at a center-of-mass energy $\sqrt{s}$ = 13.6 TeV are presented. These measurements are based on data collected with the CMS detector at the CERN LHC in 2022, corresponding to an integrated luminosity of 34.7 fb$^{-1}$. Cross sections are measured in a fiducial region closely matching the experimental acceptance, both inclusively and differentially, as a function of the transverse momentum and the absolute value of the rapidity of the four-lepton system. The H $\to$ ZZ $\to$ 4$\ell$ inclusive fiducial cross section is measured to be 2.89 $^{+0.53}_{-0.49}$ (stat) $^{+0.29}_{-0.21}$ (syst) fb, in agreement with the standard model expectation of 3.09 $^{+0.27}_{-0.24}$ fb.

Dijet angular distributions from the first LHC pp collisions at center-of-mass energy sqrt(s) = 7 TeV have been measured with the ATLAS detector. The dataset used for this analysis represents an integrated luminosity of 3.1 pb-1. Dijet $\chi$ distributions and centrality ratios have been measured up to dijet masses of 2.8 TeV, and found to be in good agreement with Standard Model predictions. Analysis of the $\chi$ distributions excludes quark contact interactions with a compositeness scale $\Lambda$ below 3.4 TeV, at 95% confidence level, significantly exceeding previous limits.

The structure of nucleons is multidimensional and depends on the transverse momenta, spatial geometry, and polarization of the constituent partons. Such a structure can be studied using high-energy photons produced in ultraperipheral heavy-ion collisions. The first measurement of the azimuthal angular correlations of exclusively produced events with two jets in photon-lead interactions at large momentum transfer is presented, a process that is considered to be sensitive to the underlying nuclear gluon polarization. This study uses a data sample of ultraperipheral lead-lead collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, corresponding to an integrated luminosity of 0.38 nb$^{-1}$, collected with the CMS experiment at the LHC. The measured second harmonic of the correlation between the sum and difference of the two jet momenta is found to be positive, and rising, as the dijet momentum increases. A well-tuned model that has been successful at describing a wide range of proton scattering data from the HERA experiments fails to describe the observed correlations, suggesting the presence of gluon polarization effects.

The Higgs boson (H) trilinear self-coupling, $\lambda_3$, is constrained via its measured properties and limits on the HH pair production using the proton-proton collision data collected by the CMS experiment at $\sqrt{s}$ = 13 TeV. The combination of event categories enriched in single-H and HH events is used to measure $\kappa_\lambda$, defined as the value of $\lambda_3$ normalized to its standard model prediction, while simultaneously constraining the Higgs boson couplings to fermions and vector bosons. Values of $\kappa_\lambda$ outside the interval $-$1.2 $\lt$$\kappa_\lambda$$\lt$ 7.5 are excluded at 2$\sigma$ confidence level, which is compatible with the expected range of $-$2.0 $\lt$$\kappa_\lambda$$\lt$ 7.7 under the assumption that all other Higgs boson couplings are equal to their standard model predicted values. Relaxing the assumption on the Higgs couplings to fermions and vector bosons the observed (expected) $\kappa_\lambda$ interval is constrained to be within $-$1.4 $\lt$$\kappa_\lambda$$\lt$ 7.8 ($-$2.3 $\lt$$\kappa_\lambda$$\lt$ 7.8) at 2$\sigma$ confidence level.

In July 2012, the ATLAS and CMS Collaborations at the CERN Large Hadron Collider announced the observation of a Higgs boson at a mass of around 125 GeV. Ten years later, and with the data corresponding to the production of 30 times larger number of Higgs bosons, we have learnt much more about the properties of the Higgs boson. The CMS experiment has observed the Higgs boson in numerous fermionic and bosonic decay channels, established its spin-parity quantum numbers, determined its mass and measured its production cross sections in various modes. Here the CMS Collaboration reports the most up-to-date combination of results on the properties of the Higgs boson, including the most stringent limit on the cross section for the production of a pair of Higgs bosons, on the basis of data from proton-proton collisions at a centre-of-mass energy of 13 TeV. Within the uncertainties, all these observations are compatible with the predictions of the standard model of elementary particle physics. Much evidence points to the fact that the standard model is a low-energy approximation of a more comprehensive theory. Several of the standard model issues originate in the sector of Higgs boson physics. An order of magnitude larger number of Higgs bosons, expected to be examined over the next fifteen years, will help deepen our understanding of this crucial sector.

The $\Xi^-_\mathrm{b} \pi^+ \pi^-$ invariant mass spectrum is investigated with an event sample of proton-proton collisions at $\sqrt{s} = $ 13 TeV, collected by the CMS experiment at the LHC in 2016-2018 and corresponding to an integrated luminosity of 140 fb$^{-1}$. The ground state $\Xi^-_\mathrm{b}$ is reconstructed via its decays to J$/\psi \Xi^-$ and J$/\psi \Xi^-\Lambda$K$^-$. A narrow resonance, labeled $\Xi_\mathrm{b}$(6100)$^-$, is observed at a $\Xi^-_\mathrm{b} \pi^+ \pi^-$ invariant mass of 6100.3 $\pm$ 0.2 (stat) $\pm$ 0.1 (syst) $\pm$ 0.6 ($\Xi^-_\mathrm{b}$) MeV, where the last uncertainty reflects the precision of the $\Xi^-_\mathrm{b}$ baryon mass. The upper limit on the $\Xi_\mathrm{b}$(6100)$^-$ natural width is determined to be 1.9 MeV at 95% confidence level. Following analogies with the established excited $\Xi_\mathrm{c}$ baryon states, the new $\Xi_\mathrm{b}$(6100)$^-$ resonance and its decay sequence are consistent with the orbitally excited $\Xi^-_\mathrm{b}$ baryon, with spin and parity quantum numbers $J^P=$ 3/2$^-$.

Protons consist of three valence quarks, two up-quarks and one down-quark, held together by gluons and a sea of quark-antiquark pairs. Collectively, quarks and gluons are referred to as partons. In a proton-proton collision, typically only one parton of each proton undergoes a hard scattering - referred to as single-parton scattering - leaving the remainder of each proton only slightly disturbed. Here, we report the study of double- and triple-parton scatterings through the simultaneous production of three J/$\psi$ mesons, which consist of a charm quark-antiquark pair, in proton-proton collisions recorded with the CMS experiment at the Large Hadron Collider. We observed this process - reconstructed through the decays of J/$\psi$ mesons into pairs of oppositely charged muons - with a statistical significance above five standard deviations. We measured the inclusive fiducial cross section to be 272 $^{+141}_{-104}$ (stat) $\pm$ 17 (syst) fb, and compared it to theoretical expectations for triple-J/$\psi$ meson production in single-, double- and triple-parton scattering scenarios. Assuming factorization of multiple hard-scattering probabilities in terms of single-parton scattering cross sections, double- and triple-parton scattering are the dominant contributions for the measured process.

The first observation of the decay $\Xi^-_\mathrm{b}$$\to$$\psi$(2S)$\Xi^-$ and measurement of the branching ratio of $\Xi^-_\mathrm{b}$$\to$$\psi$(2S)$\Xi^-$ to $\Xi^-_\mathrm{b}$$\to$ J/$\psi$$\Xi^-$ are presented. The J/$\psi$ and $\psi$(2S) mesons are reconstructed using their dimuon decay modes. The results are based on proton-proton colliding beam data from the LHC collected by the CMS experiment at $\sqrt{s}$ = 13 TeV in 2016-2018, corresponding to an integrated luminosity of 140 fb$^{-1}$. The branching fraction ratio is measured to be $\mathcal{B}$($\Xi^-_\mathrm{b}$$\to$$\psi$(2S)$\Xi^-$) / $\mathcal{B}$($\Xi^-_\mathrm{b}$$\to$ J/$\psi$$\Xi^-$) = 0.84 $^{+0.21}_{-0.19}$ (stat) $\pm$ 0.10 (syst) $\pm$ 0.02 ($\mathcal{B}$), where the last uncertainty comes from the uncertainties in the branching fractions of the charmonium states. New measurements of the $\Xi_\mathrm{b}$(5945)$^{0}$ baryon mass and natural width are also presented, using the $\Xi_\mathrm{b}^-\pi^+$ final state, where the $\Xi^-_\mathrm{b}$ baryon is reconstructed through the decays J/$\psi \Xi^-$, $\psi$(2S)$\Xi^-$, J/$\psi \Lambda$K$^-$, and J/$\psi \Sigma^0$K$^-$. Finally, the fraction of the $\Xi^-_\mathrm{b}$ baryons produced from $\Xi_\mathrm{b}$(5945)$^{0}$ decays is determined.

This paper reports searches for heavy resonances decaying into $ZZ$ or $ZW$ using data from proton--proton collisions at a centre-of-mass energy of $\sqrt{s}=13$ TeV. The data, corresponding to an integrated luminosity of 36.1 fb$^{-1}$, were recorded with the ATLAS detector in 2015 and 2016 at the Large Hadron Collider. The searches are performed in final states in which one $Z$ boson decays into either a pair of light charged leptons (electrons and muons) or a pair of neutrinos, and the associated $W$ boson or the other $Z$ boson decays hadronically. No evidence of the production of heavy resonances is observed. Upper bounds on the production cross sections of heavy resonances times their decay branching ratios to $ZZ$ or $ZW$ are derived in the mass range 300--5000 GeV within the context of Standard Model extensions with additional Higgs bosons, a heavy vector triplet or warped extra dimensions. Production through gluon--gluon fusion, Drell--Yan or vector-boson fusion are considered, depending on the assumed model.

Jet shapes have been measured in inclusive jet production in proton-proton collisions at sqrt(s) = 7 TeV using 3 pb^{-1} of data recorded by the ATLAS experiment at the LHC. Jets are reconstructed using the anti-kt algorithm with transverse momentum 30 GeV < pT < 600 GeV and rapidity in the region |y| < 2.8. The data are corrected for detector effects and compared to several leading-order QCD matrix elements plus parton shower Monte Carlo predictions, including different sets of parameters tuned to model fragmentation processes and underlying event contributions in the final state. The measured jets become narrower with increasing jet transverse momentum and the jet shapes present a moderate jet rapidity dependence. Within QCD, the data test a variety of perturbative and non-perturbative effects. In particular, the data show sensitivity to the details of the parton shower, fragmentation, and underlying event models in the Monte Carlo generators. For an appropriate choice of the parameters used in these models, the data are well described.

The first observation of the concurrent production of two J/$\psi$ mesons in proton-nucleus collisions is presented. The analysis is based on a proton-lead (pPb) data sample recorded at a nucleon-nucleon center-of-mass energy of 8.16 TeV by the CMS experiment at the CERN LHC and corresponding to an integrated luminosity of 174.6 nb$^{-1}$. The two J/$\psi$ mesons are reconstructed in their $\mu^+\mu^-$ decay channels with transverse momenta $p_\mathrm{T}$$\gt$ 6.5 GeV and rapidity $\lvert y \rvert$$\lt$ 2.4. Events where one of the J/$\psi$ mesons is reconstructed in the dielectron channel are also considered in the search. The pPb $\to$ J/$\psi$J/$\psi$ + X process is observed with a significance of 5.3 standard deviations. The measured inclusive fiducial cross section, using the four-muon channel alone, is $\sigma$(pPb $\to$ J/$\psi$J/$\psi$ + X) = 22.0 $\pm$ 8.9 (stat) $\pm$ 1.5 (syst) nb. A fit of the data to the expected rapidity separation for pairs of J/$\psi$ mesons produced in single (SPS) and double (DPS) parton scatterings yields $\sigma^{\mathrm{pPb} \to \mathrm{J}/\psi\mathrm{J}/\psi +\mathrm{X}}_\text{SPS}$ = 16.5 $\pm$ 10.8 (stat) $\pm$ 0.1 (syst) nb and $\sigma^{\mathrm{pPb} \to \mathrm{J}/\psi\mathrm{J}/\psi + \mathrm{X}}_\text{DPS}$ = 5.4 $\pm$ 6.2 (stat) $\pm$ 0.4 (syst) nb, respectively. This latter result can be transformed into a lower bound on the effective DPS cross section, closely related to the squared average interparton transverse separation in the collision, of $\sigma_\text{eff}$$\gt$ 1.0 mb at 95% confidence level.

A pioneering machine-learning-based flavor-tagging algorithm combining same-side and opposite-side tagging is used to obtain the equivalent of 27$\,$000 tagged B$^0_\mathrm{s}$$\to$ J/$\psi\, \phi$(1020) decays from pp collisions at $\sqrt{s} =$ 13 TeV, collected by the CMS experiment and corresponding to an integrated luminosity of 96.5 fb$^{-1}$. A time- and flavor-dependent angular analysis of the $\mu^+\mu^-$K$^+$K$^-$ final state is used to measure parameters of the $\mathrm{B}^0_\mathrm{s}$-$\overline{\mathrm{B}}^0_\mathrm{s}$ system. The weak phase is measured to be $\phi_\mathrm{s}$ = $-$73 $\pm$ 23 (stat) $\pm$ 7 (syst) mrad, which, combined with a $\sqrt{s}$ = 8 TeV CMS result, gives $\phi_\mathrm{s}$ = $-$74 $\pm$ 23 mrad. This value differs from zero by 3.2 standard deviations, providing evidence for $CP$ violation in B$^0_\mathrm{s}$$\to$ J/$\psi\,\phi$(1020) decays. All measured physics parameters are found to agree with standard model predictions where available.

A search for the standard model Higgs boson decaying to a charm quark-antiquark pair, H $\to$$\mathrm{c\bar{c}}$, produced in association with a leptonically decaying V (W or Z) boson is presented. The search is performed with proton-proton collisions at $\sqrt{s}$ = 13 TeV collected by the CMS experiment, corresponding to an integrated luminosity of 138 fb$^{-1}$. Novel charm jet identification and analysis methods using machine learning techniques are employed. The analysis is validated by searching for Z $\to$$\mathrm{c\bar{c}}$ in VZ events, leading to its first observation at a hadron collider with a significance of 5.7 standard deviations. The observed (expected) upper limit on $\sigma$(VH)$ \mathcal{B}$(H $\to$$\mathrm{c\bar{c}}$) is 0.94 (0.50 $^{+0.22}_{-0.15}$) pb at 95% confidence level (CL), corresponding to 14 (7.6 $^{+3.4}_{-2.3}$) times the standard model prediction. For the Higgs-charm Yukawa coupling modifier, $\kappa_\mathrm{c}$, the observed (expected) 95% CL interval is 1.1 $\lt$ $\vert\kappa_\mathrm{C}\vert$ $\lt$ 5.5 ($\vert\kappa_\mathrm{c}\vert$ $\lt$ 3.4), the most stringent constraint to date.

The polarizations of prompt and non-prompt J$/\psi$ and $\psi$(2S) mesons are measured in proton-proton collisions at $\sqrt{s}$ = 13 TeV, using data samples collected by the CMS experiment in 2017 and 2018, corresponding to a total integrated luminosity of 103.3 fb$^{-1}$. Based on the analysis of the dimuon decay angular distributions in the helicity frame, the polar anisotropy, $\lambda_\theta$, is measured as a function of the transverse momentum, $p_\mathrm{T}$, of the charmonium states, in the 25-120 and 20-100 GeV ranges for the J$/\psi$ and $\psi$(2S), respectively. The non-prompt polarizations agree with predictions based on the hypothesis that, for $p_\mathrm{T}$$\gtrsim$ 25 GeV, the non-prompt J$/\psi$ and $\psi$(2S) are predominantly produced in two-body B meson decays. The prompt results clearly exclude strong transverse polarizations, even for $p_\mathrm{T}$ exceeding 30 times the J$/\psi$ mass, where $\lambda_\theta$ tends to an asymptotic value around 0.3. Taken together with previous measurements, by CMS and LHCb at $\sqrt{s}$ = 7 TeV, the prompt polarizations show a significant variation with $p_\mathrm{T}$, at low $p_\mathrm{T}$.

The J/$\psi$$\to$$\mu^+\mu^-\mu^+\mu^-$ decay has been observed with a statistical significance in excess of five standard deviations. The analysis is based on an event sample of proton-proton collisions at a center-of-mass energy of 13 TeV, collected by the CMS experiment in 2018 and corresponding to an integrated luminosity of 33.6 fb$^{-1}$. Normalizing to the J/$\psi$$\to$$\mu^+\mu^-$ decay mode leads to a branching fraction [10.1 $^{+3.3}_{-2.7}$ (stat) $\pm$ 0.4 (syst)] $\times$ 10$^{-7}$, a value that is consistent with the standard model prediction.

This paper presents a measurement of the triple-differential cross section for the Drell--Yan process $Z/\gamma^*\rightarrow \ell^+\ell^-$ where $\ell$ is an electron or a muon. The measurement is performed for invariant masses of the lepton pairs, $m_{\ell\ell}$, between $46$ and $200$ GeV using a sample of $20.2$ fb$^{-1}$ of $pp$ collisions data at a centre-of-mass energy of $\sqrt{s}=8$ TeV collected by the ATLAS detector at the LHC in 2012. The data are presented in bins of invariant mass, absolute dilepton rapidity, $|y_{\ell\ell}|$, and the angular variable $\cos\theta^{*}$ between the outgoing lepton and the incoming quark in the Collins--Soper frame. The measurements are performed in the range $|y_{\ell\ell}|<2.4$ in the muon channel, and extended to $|y_{\ell\ell}|<3.6$ in the electron channel. The cross sections are used to determine the $Z$ boson forward-backward asymmetry as a function of $|y_{\ell\ell}|$ and $m_{\ell\ell}$. The measurements achieve high-precision, below the percent level in the pole region, excluding the uncertainty in the integrated luminosity, and are in agreement with predictions. These precision data are sensitive to the parton distribution functions and the effective weak mixing angle.

We present an observation of photon-photon production of $\tau$ lepton pairs in ultraperipheral lead-lead collisions. The measurement is based on a data sample with an integrated luminosity of 404 $\mu$b$^{-1}$ collected by the CMS experiment at a nucleon-nucleon center-of-mass energy of 5.02 TeV. The $\gamma\gamma$$\to$$\tau^+\tau^-$ process is observed for $\tau\tau$ events with a muon and three charged hadrons in the final state. The measured fiducial cross section is $\sigma(\gamma\gamma$$\to$$\tau^+\tau^-)$ = 4.8 $\pm$ 0.6 (stat) $\pm$ 0.5 (syst) $\mu$b, in agreement with leading-order QED predictions. Using $\sigma(\gamma\gamma$$\to$$\tau^+\tau^-)$, we estimate a model-dependent value of the anomalous magnetic moment of the $\tau$ lepton of $a_\tau$ = 0.001 $^{+0.055}_{-0.089}$.

An analysis is presented based on models of the intrinsic transverse momentum (intrinsic $k_\mathrm{T}$) of partons in nucleons by studying the dilepton transverse momentum in Drell-Yan events. Using parameter tuning in event generators and existing data from fixed-target experiments and from hadron colliders, our investigation spans three orders of magnitude in center-of-mass energy and two orders of magnitude in dilepton invariant mass. The results show an energy-scaling behavior of the intrinsic $k_\mathrm{T}$ parameters, independent of the dilepton invariant mass at a given center-of-mass energy.

A full set of optimized observables is measured in an angular analysis of the decay B$^0$$\to$ K$^*$(892)$^0\mu^+\mu^-$ using a sample of proton-proton collisions at $\sqrt{s}$ = 13 TeV, collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 140 fb$^{-1}$. The analysis is performed in six bins of the squared invariant mass of the dimuon system, $q^2$, over the range 1.1 $\lt$$q^2$$\lt$ 16 GeV$^2$. The results are among the most precise experimental measurements of the angular observables for this decay and are compared to a variety of predictions based on the standard model. Some of these predictions exhibit tension with the measurements.