Signals consistent with the B$^+_\mathrm{c}$(2S) and B$^{*+}_\mathrm{c}$(2S) states are observed in proton-proton collisions at $\sqrt{s} =$ 13 TeV, in an event sample corresponding to an integrated luminosity of 140 fb$^{-1}$, collected by the CMS experiment during the 2016, 2017, and 2018 LHC running periods. These excited $\bar{\mathrm{b}}$c states are observed in the B$^+_\mathrm{c}\pi^+\pi^-$ invariant mass spectrum, with the ground state B$^+_\mathrm{c}$ reconstructed through its decay to J/$\psi\,\pi^+$. The two states are well resolved from each other and are observed with a significance exceeding five standard deviations. The mass of the B$^+_\mathrm{c}$(2S) meson is measured to be 6871.0 $\pm$ 1.2 (stat) $\pm$ 0.8 (syst) $\pm$ 0.8 (B$^+_\mathrm{c}$) MeV, where the last term corresponds to the uncertainty in the world-average B$^+_\mathrm{c}$ mass.

A measurement is presented of electroweak (EW) production of a W boson in association with two jets in proton-proton collisions at $\sqrt{s} =$ 13 TeV. The data sample was recorded by the CMS Collaboration at the LHC and corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The measurement is performed for the $\ell\nu$jj final state (with $\ell\nu$ indicating a lepton-neutrino pair, and j representing the quarks produced in the hard interaction) in a kinematic region defined by invariant mass $m_\mathrm{jj}$ $>$ 120 GeV and transverse momenta $p_\mathrm{T j}$ $>$ 25 GeV. The cross section of the process is measured in the electron and muon channels yielding $\sigma_\mathrm{EW}$(Wjj) = 6.23 $\pm$ 0.12 (stat) $\pm$ 0.61 (syst) pb per channel, in agreement with leading-order standard model predictions. The additional hadronic activity of events in a signal-enriched region is studied, and the measurements are compared with predictions. The final state is also used to perform a search for anomalous trilinear gauge couplings. Limits on anomalous trilinear gauge couplings associated with dimension-six operators are given in the framework of an effective field theory. The corresponding 95% confidence level intervals are $-$2.3 $<$ $c_{\mathrm{WWW}}/\Lambda^2$ $<$ 2.5 TeV$^{-2}$, $-$8.8 $<$ $c_{\mathrm{W}}/\Lambda^2$ $<$ 16 TeV$^{-2}$, and $-$45 $<$ $c_{\mathrm{B}}/\Lambda^2$ $<$ 46 TeV$^{-2}$. These results are combined with the CMS EW Zjj analysis, yielding the constraint on the $c_{\mathrm{WWW}}$ coupling: $-$1.8 $<$ $c_{\mathrm{WWW}}/\Lambda^2$ $<$ 2.0 TeV$^{-2}$.

A search for anomalous electroweak production of WW, WZ, and ZZ boson pairs in association with two jets in proton-proton collisions at $\sqrt{s} =$ 13 TeV at the LHC is reported. The data sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$ collected with the CMS detector. Events are selected by requiring two jets with large rapidity separation and invariant mass, one or two leptons (electrons or muons), and a W or Z boson decaying hadronically. No excess of events with respect to the standard model background predictions is observed and constraints on the structure of quartic vector boson interactions in the framework of dimension-8 effective field theory operators are reported. Stringent limits on parameters of the effective field theory operators are obtained. The observed 95% confidence level limits for the S0, M0, and T0 operators are $-$2.7 $<$ f$_{\mathrm{S0}}/ \Lambda^{4}$ $<$ 2.7, $-$1.0 $<$ f$_{\mathrm{M0}}/ \Lambda^{4}$ $<$ 1.0, and $-$0.17 $<$ f$_{\mathrm{T0}}/ \Lambda^{4}$ $<$ 0.16, in units of TeV$^{-4}$. Constraints are also reported on the product of the cross section and branching fraction for vector boson fusion production of charged Higgs bosons as a function of mass from 600 to 2000 GeV. The results are interpreted in the context of the Georgi-Machacek model.

Measurements of the top quark polarization and top quark pair ($\mathrm{t\bar{t}}$) spin correlations are presented using events containing two oppositely charged leptons (e$^+$e$^-$, e$^\pm\mu^\mp$, or $\mu^+\mu^-$) produced in proton-proton collisions at a center-of-mass energy of 13 TeV. The data were recorded by the CMS experiment at the LHC in 2016 and correspond to an integrated luminosity of 35.9 fb$^{-1}$. A set of parton-level normalized differential cross sections, sensitive to each of the independent coefficients of the spin-dependent parts of the $\mathrm{t\bar{t}}$ production density matrix, is measured for the first time at 13 TeV. The measured distributions and extracted coefficients are compared with standard model predictions from simulations at next-to-leading-order (NLO) accuracy in quantum chromodynamics (QCD), and from NLO QCD calculations including electroweak corrections. All measurements are found to be consistent with the expectations of the standard model. The normalized differential cross sections are used in fits to constrain the anomalous chromomagnetic and chromoelectric dipole moments of the top quark to $-$0.24 $<C_\text{tG}/\Lambda^{2}$ $<$ 0.07 TeV$^{-2}$ and $-$0.33 $< C^{I}_\text{tG}/\Lambda^{2}$ $<$ 0.20 TeV$^{-2}$, respectively, at 95% confidence level.

For the first time, a search for the rare decay of the W boson to three charged pions has been performed. Proton-proton collision data recorded by the CMS experiment at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 77.3 fb$^{-1}$, have been analyzed. No significant excess is observed above the background expectation. An upper limit of 1.01 $\times$ 10$^{-6}$ is set at 95% confidence level on the branching fraction of the W boson to three charged pions. This provides a strong motivation for theoretical calculations of this branching fraction.

A search for nonresonant excesses in the invariant mass spectra of electron and muon pairs is presented. The analysis is based on data from proton-proton collisions at a center-of-mass energy of 13 TeV recorded by the CMS experiment in 2016, corresponding to a total integrated luminosity of 36 fb$^{-1}$. No significant deviation from the standard model is observed. Limits are set at 95% confidence level on energy scales for two general classes of nonresonant models. For a class of fermion contact interaction models, lower limits ranging from 20 to 32 TeV are set on the characteristic compositeness scale $\Lambda$. For the Arkani-Hamed, Dimopoulos, and Dvali model of large extra dimensions, the first results in the dilepton final state at 13 TeV are reported, and values of the ultraviolet cutoff parameter $\Lambda_\mathrm{T}$ below 6.9 TeV are excluded. A combination with recent CMS diphoton results improves this exclusion to $\Lambda_\mathrm{T}$ below 7.7 TeV, providing the most sensitive limits to date in nonhadronic final states.

A search for the pair production of heavy vector-like partners T and B of the top and bottom quarks has been performed by the CMS experiment at the CERN LHC using proton-proton collisions at $\sqrt{s} =$ 13 TeV. The data sample was collected in 2016 and corresponds to an integrated luminosity of 35.9 fb$^{-1}$. Final states studied for $\mathrm{T\overline{T}}$ production include those where one of the T quarks decays via T $\to$ tZ and the other via T $\to$ bW, tZ, or tH, where H is a Higgs boson. For the $\mathrm{B\overline{B}}$ case, final states include those where one of the B quarks decays via B $\to$ bZ and the other B $\to$ tW, bZ, or bH. Events with two oppositely charged electrons or muons, consistent with coming from the decay of a Z boson, and jets are investigated. The number of observed events is consistent with standard model background estimations. Lower limits at 95% confidence level are placed on the masses of the T and B quarks for a range of branching fractions. Assuming 100% branching fractions for T $\to$ tZ, and B $\to$ bZ, T and B quark mass values below 1280 and 1130 GeV, respectively, are excluded.

A search for supersymmetry is presented based on events with at least one photon, jets, and large missing transverse momentum produced in proton-proton collisions at a center-of-mass energy of 13 TeV. The data correspond to an integrated luminosity of 35.9 fb$^{-1}$ and were recorded at the LHC with the CMS detector in 2016. The analysis characterizes signal-like events by categorizing the data into various signal regions based on the number of jets, the number of b-tagged jets, and the missing transverse momentum. No significant excess of events is observed with respect to the expectations from standard model processes. Limits are placed on the gluino and top squark pair production cross sections using several simplified models of supersymmetric particle production with gauge-mediated supersymmetry breaking. Depending on the model and the mass of the next-to-lightest supersymmetric particle, the production of gluinos with masses as large as 2120 GeV and the production of top squarks with masses as large as 1230 GeV are excluded at 95% confidence level.

A top quark mass measurement is performed using 35.9 fb$^{-1}$ of LHC proton-proton collision data collected with the CMS detector at $\sqrt{s} =$ 13 TeV. The measurement uses the $\mathrm{t\overline{t}}$ all-jets final state. A kinematic fit is performed to reconstruct the decay of the $\mathrm{t\overline{t}}$ system and suppress the multijet background. Using the ideogram method, the top quark mass ($m_\mathrm{t}$) is determined, simultaneously constraining an additional jet energy scale factor (JSF). The resulting value of $m_\mathrm{t}$ = 172.34 $\pm$ 0.20 (stat+JSF) $\pm$ 0.70 (syst) GeV is in good agreement with previous measurements. In addition, a combined measurement that uses the $\mathrm{t\overline{t}}$ lepton+jets and all-jets final states is presented, using the same mass extraction method, and provides an $m_\mathrm{t}$ measurement of 172.26 $\pm$ 0.07 (stat+JSF) $\pm$ 0.61 (syst) GeV. This is the first combined $m_\mathrm{t}$ extraction from the lepton+jets and all-jets channels through a single likelihood function.

A general search is presented for a low-mass $\tau^-\tau^+$ resonance produced in association with a bottom quark. The search is based on proton-proton collision data at a center-of-mass energy of 13 TeV collected by the CMS experiment at the LHC, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The data are consistent with the standard model expectation. Upper limits at 95% confidence level on the cross section times branching fraction are determined for two signal models: a light pseudoscalar Higgs boson decaying to a pair of $\tau$ leptons produced in association with bottom quarks, and a low-mass boson X decaying to a $\tau$-lepton pair that is produced in the decay of a bottom-like quark B such that B $\to$ bX. Masses between 25 and 70 GeV are probed for the light pseudoscalar boson with upper limits ranging from 250 to 44 pb. Upper limits from 20 to 0.3 pb are set on B masses between 170 and 450 GeV for X boson masses between 20 and 70 GeV.

A search for supersymmetric particles produced in the vector boson fusion topology in proton-proton collisions is presented. The search targets final states with one or zero leptons, large missing transverse momentum, and two jets with a large separation in rapidity. The data sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$ of proton-proton collisions at $\sqrt{s} =$ 13 TeV collected in 2016 with the CMS detector at the LHC. The observed dijet invariant mass and lepton-neutrino transverse mass spectra are found to be consistent with the standard model predictions. Upper limits are set on the cross sections for chargino ($\widetilde\chi_1^\pm$) and neutralino ($\widetilde\chi_2^0$) production with two associated jets. For a compressed mass spectrum scenario in which the $\widetilde\chi_1^\pm$ and $\widetilde\chi_2^0$ decays proceed via a light slepton and the mass difference between the lightest neutralino $\widetilde\chi_1^0$ and the mass-degenerate particles $\widetilde\chi_1^\pm$ and $\widetilde\chi_2^0$ is 1 (30) GeV, the most stringent lower limit to date of 112 (215) GeV is set on the mass of these latter two particles.

A search for long-lived particles decaying into jets is presented. Data were collected with the CMS detector at the LHC from proton-proton collisions at a center-of-mass energy of 13 TeV in 2016, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The search examines the distinctive topology of displaced tracks and secondary vertices. The selected events are found to be consistent with standard model predictions. For a simplified model in which long-lived neutral particles are pair produced and decay to two jets, pair production cross sections larger than 0.2 fb are excluded at 95% confidence level for a long-lived particle mass larger than 1000 GeV and proper decay lengths between 3 and 130 mm. Several supersymmetry models with gauge-mediated supersymmetry breaking or $R$-parity violation, where pair-produced long-lived gluinos or top squarks decay to several final-state topologies containing displaced jets, are also tested. For these models, in the mass ranges above 200 GeV, gluino masses up to 2300-2400 GeV and top squark masses up to 1350-1600 GeV are excluded for proper decay lengths approximately between 10 and 100 mm. These are the most restrictive limits to date on these models.

A search is presented for charged Higgs bosons in the H$^{\pm}$ $\to$ $\tau^{\pm}\nu_\tau$ decay mode in the hadronic final state and in final states with an electron or muon. The search is based on proton-proton collision data recorded by the CMS experiment in 2016 at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The results agree with the background expectation from the standard model. Upper limits at 95% confidence level are set on the production cross section times branching fraction to $\tau^{\pm}\nu_\tau$ for an H$^{\pm}$ in the mass range of 80 GeV to 3 TeV, including the region near the top quark mass. The observed limit ranges from 6 pb at 80 GeV to 5 fb at 3 TeV. The limits are interpreted in the context of the minimal supersymmetric standard model $m_\mathrm{h}^\mathrm{mod-}$ scenario.

A search for exotic decays of the Higgs boson to a pair of light pseudoscalar particles a$_1$ is performed under the hypothesis that one of the pseudoscalars decays to a pair of opposite sign muons and the other decays to b$\overline{\mathrm{b}}$. Such signatures are predicted in a number of extensions of the standard model (SM), including next-to-minimal supersymmetry and two-Higgs-doublet models with an additional scalar singlet. The results are based on a data set of proton-proton collisions corresponding to an integrated luminosity of 35.9 fb$^{-1}$, accumulated with the CMS experiment at the CERN LHC in 2016 at a centre-of-mass energy of 13 TeV. No statistically significant excess is observed with respect to the SM backgrounds in the search region for pseudoscalar masses from 20 GeV to half of the Higgs boson mass. Upper limits at 95% confidence level are set on the product of the production cross section and branching fraction, $\sigma_{\mathrm{h}}\mathcal{B}$(h $\to$ a$_1$ a$_1$ $\to$ $\mu^+\mu^-\mathrm{b}\bar{\mathrm{b}}$), ranging from 5 to 33 fb, depending on the pseudoscalar mass. Corresponding limits on the branching fraction, assuming the SM prediction for $\sigma_{\mathrm{h}}$, are (1$-$7)$\times$ 10$^{-4}$.

The observation of the standard model (SM) Higgs boson decay to a pair of bottom quarks is presented. The main contribution to this result is from processes in which Higgs bosons are produced in association with a W or Z boson (VH), and are searched for in final states including 0, 1, or 2 charged leptons and two identified bottom quark jets. The results from the measurement of these processes in a data sample recorded by the CMS experiment in 2017, comprising 41.3 fb$^{-1}$ of proton-proton collisions at $\sqrt{s} =$ 13 TeV, are described. When combined with previous VH measurements using data collected at $\sqrt{s}=$ 7, 8, and 13 TeV, an excess of events is observed at $m_\mathrm{H} =$ 125.09 GeV with a significance of 4.8 standard deviations, where the expectation for the SM Higgs boson is 4.9. The corresponding measured signal strength is 1.01 $\pm$ 0.22. The combination of this result with searches by the CMS experiment for H $\to\mathrm{b\overline{b}}$ in other production processes yields an observed (expected) significance of 5.6 (5.5) standard deviations and a signal strength of 1.04 $\pm$ 0.20.

A search is presented for pairs of light pseudoscalar bosons, in the mass range from 4 to 15 GeV, produced from decays of the 125 GeV Higgs boson. The decay modes considered are final states that arise when one of the pseudoscalars decays to a pair of tau leptons, and the other one either into a pair of tau leptons or muons. The search is based on proton-proton collisions collected by the CMS experiment in 2016 at a center-of-mass energy of 13 TeV that correspond to an integrated luminosity of 35.9 fb${-1}$. The 2$\mu$2$\tau$ and 4$\tau$ channels are used in combination to constrain the product of the Higgs boson production cross section and the branching fraction into 4$\tau$ final state, $\sigma\mathcal{B}$, exploiting the linear dependence of the fermionic coupling strength of pseudoscalar bosons on the fermion mass. No significant excess is observed beyond the expectation from the standard model. The observed and expected upper limits at 95% confidence level on $\sigma\mathcal{B}$, relative to the standard model Higgs boson production cross section, are set respectively between 0.022 and 0.23 and between 0.027 and 0.19 in the mass range probed by the analysis.

A measurement is presented of the associated production of a single top quark and a W boson in proton-proton collisions at $\sqrt{s}=$ 13 TeV by the CMS Collaboration at the CERN LHC. The data collected corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The measurement is performed using events with one electron and one muon in the final state along with at least one jet originated from a bottom quark. A multivariate discriminant, exploiting the kinematic properties of the events, is used to separate the signal from the dominant $\mathrm{t\overline{t}}$ background. The measured cross section of 63.1 $\pm$ 1.8 (stat) $\pm$ 6.4 (syst) $\pm$ 2.1 (lumi) pb is in agreement with the standard model expectation.

The production of a Z boson, decaying to two charged leptons, in association with jets in proton-proton collisions at a centre-of-mass energy of 13 TeV is measured. Data recorded with the CMS detector at the LHC are used that correspond to an integrated luminosity of 2.19 fb$^{-1}$. The cross section is measured as a function of the jet multiplicity and its dependence on the transverse momentum of the Z boson, the jet kinematic variables (transverse momentum and rapidity), the scalar sum of the jet momenta, which quantifies the hadronic activity, and the balance in transverse momentum between the reconstructed jet recoil and the Z boson. The measurements are compared with predictions from four different calculations. The first two merge matrix elements with different parton multiplicities in the final state and parton showering, one of which includes one-loop corrections. The third is a fixed-order calculation with next-to-next-to-leading order accuracy for the process with a Z boson and one parton in the final state. The fourth combines the fully differential next-to-next-to-leading order calculation with next-to-next-to-leading logarithm resummation and parton showering.

A search for the Higgs boson decaying to two oppositely charged muons is presented using data recorded by the CMS experiment at the CERN LHC in 2016 at a center-of-mass energy $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Data are found to be compatible with the predicted background. For a Higgs boson with a mass of 125.09 GeV, the 95% confidence level observed (background-only expected) upper limit on the production cross section times branching fraction to a pair of muons is found to be 3.0 (2.5) times the standard model expectation. In combination with data recorded at center-of-mass energies $\sqrt{s} =$ 7 and 8 TeV, the background-only expected upper limit improves to 2.2 times the standard model value with a standard model expected significance of 1.0 standard deviations. The corresponding observed upper limit is 2.9 with an observed significance of 0.9 standard deviations. This corresponds to an observed upper limit on the standard model Higgs boson branching fraction to muons of 6.4 $\times$ 10$^{-4}$ and to an observed signal strength of 1.0 $\pm$ 1.0 (stat) $\pm$ 0.1 (syst).

A search for Higgs bosons that decay into a bottom quark-antiquark pair and are accompanied by at least one additional bottom quark is performed with the CMS detector. The data analyzed were recorded in proton-proton collisions at a centre-of-mass energy of $\sqrt{s} =$ 13 TeV at the LHC, corresponding to an integrated luminosity of 35.7 fb$^{-1}$. The final state considered in this analysis is particularly sensitive to signatures of a Higgs sector beyond the standard model, as predicted in the generic class of two Higgs doublet models (2HDMs). No signal above the standard model background expectation is observed. Stringent upper limits on the cross section times branching fraction are set for Higgs bosons with masses up to 1300 GeV. The results are interpreted within several MSSM and 2HDM scenarios.

Bose-Einstein correlations of charged hadrons are measured over a broad multiplicity range, from a few particles up to about 250 reconstructed charged hadrons in proton-proton collisions at $\sqrt{s} =$ 13 TeV. The results are based on data collected using the CMS detector at the LHC during runs with a special low-pileup configuration. Three analysis techniques with different degrees of dependence on simulations are used to remove the non-Bose-Einstein background from the correlation functions. All three methods give consistent results. The measured lengths of homogeneity are studied as functions of particle multiplicity as well as average pair transverse momentum and mass. The results are compared with data from both CMS and ATLAS at $\sqrt{s} =$ 7 TeV, as well as with theoretical predictions.

A search for an exotic decay of the Higgs boson to a pair of light pseudoscalar bosons is performed for the first time in the final state with two b quarks and two $\tau$ leptons. The search is motivated in the context of models of physics beyond the standard model (SM), such as two Higgs doublet models extended with a complex scalar singlet (2HDM+S), which include the next-to-minimal supersymmetric SM (NMSSM). The results are based on a data set of proton-proton collisions corresponding to an integrated luminosity of 35.9 fb$^{-1}$, accumulated by the CMS experiment at the LHC in 2016 at a center-of-mass energy of 13 TeV. Masses of the pseudoscalar boson between 15 and 60 GeV are probed, and no excess of events above the SM expectation is observed. Upper limits between 3 and 12% are set on the branching fraction $\mathcal{B}$(h $\to$ aa $\to$ 2$\tau$2b) assuming the SM production of the Higgs boson. Upper limits are also set on the branching fraction of the Higgs boson to two light pseudoscalar bosons in different 2HDM+S scenarios. Assuming the SM production cross section for the Higgs boson, the upper limit on this quantity is as low as 20% for a mass of the pseudoscalar of 40 GeV in the NMSSM.

Measurements of inclusive isolated-photon and photon+jet production in proton-proton collisions at $\sqrt{s} =$ 13 TeV are presented. The analysis uses data collected by the CMS experiment in 2015, corresponding to an integrated luminosity of 2.26 fb$^{-1}$. The cross section for inclusive isolated photon production is measured as a function of the photon transverse energy in a fiducial region. The cross section for photon+jet production is measured as a function of the photon transverse energy in the same fiducial region with identical photon requirements and with the highest transverse momentum jet. All measurements are in agreement with predictions from next-to-leading-order perturbative QCD.

A study is presented of anomalous HVV interactions of the Higgs boson, including its $CP$ properties. The study uses Higgs boson candidates produced mainly in vector boson fusion and gluon fusion that subsequently decay to a pair of $\tau$ leptons. The data were recorded by the CMS experiment at the LHC in 2016 at a center-of-mass energy of 13 TeV and correspond to an integrated luminosity of 35.9 fb$^{-1}$. A matrix element technique is employed for the analysis of anomalous interactions. The results are combined with those from the H $\to 4\ell$ decay channel presented earlier, yielding the most stringent constraints on anomalous Higgs boson couplings to electroweak vector bosons expressed as effective cross section fractions and phases: the $CP$-violating parameter $f_{a3}\cos(\phi_{a3})$ $=$ $(0.00 \pm 0.27) \times 10^{-3}$ and the $CP$-conserving parameters $f_{a2}\cos(\phi_{a2})$ $=$ $(0.08 ^{+1.04}_{-0.21}) \times 10^{-3}$, $f_{\Lambda1}\cos(\phi_{\Lambda1})$ $=$ $(0.00 ^{+0.53}_{-0.09}) \times 10^{-3}$, and $f_{\Lambda1}^{\mathrm{Z}\gamma}\cos(\phi_{\Lambda1}^{\mathrm{Z}\gamma})$ $=$ $(0.0 ^{+1.1}_{-1.3}) \times 10^{-3}$. The current dataset does not allow for precise constraints on $CP$ properties in the gluon fusion process. The results are consistent with standard model expectations.

A search for a Higgs boson decaying into a pair of electrons or muons and a photon is described. Higgs boson decays to a Z boson and a photon (H $\to$ Z$\gamma\to\ell\ell\gamma$, $\ell =$ e or $\mu$), or to two photons, one of which has an internal conversion into a muon pair (H $\to\gamma^{*}\gamma\to\mu\mu\gamma$) were considered. The analysis is performed using a data set recorded by the CMS experiment at the LHC from proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. No significant excess above the background prediction has been found. Limits are set on the cross section for a standard model Higgs boson decaying to opposite-sign electron or muon pairs and a photon. The observed limits on cross section times the corresponding branching fractions vary between 1.4 and 4.0 (6.1 and 11.4) times the standard model cross section for H $\to\gamma^{*}\gamma\to\mu\mu\gamma$ (H $\to$ Z$\gamma\to\ell\ell\gamma$) in the 120-130 GeV mass range of the $\ell\ell\gamma$ system. The H $\to\gamma^*\gamma\to\mu\mu\gamma$ and H $\to$ Z$\gamma\to\ell\ell\gamma$ analyses are combined for $m_\mathrm{H} =$ 125 GeV, obtaining an observed (expected) 95% confidence level upper limit of 3.9 (2.0) times the standard model cross section.