Production cross sections of $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) states decaying into $\mu^+\mu^-$ in proton-lead (pPb) collisions are reported using data collected by the CMS experiment at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV. A comparison is made with corresponding cross sections obtained with pp data measured at the same collision energy and scaled by the Pb nucleus mass number. The nuclear modification factor for $\Upsilon$(1S) is found to be $R_\mathrm{pPb}(\Upsilon(1S))$ = 0.806 $\pm$ 0.024 (stat) $\pm$ 0.059 (syst). Similar results for the excited states indicate a sequential suppression pattern, such that $R_\mathrm{pPb}(\Upsilon(1S))$$\gt$$R_\mathrm{pPb}(\Upsilon(2S))$$\gt$$R_\mathrm{pPb}(\Upsilon(3S))$. The suppression is much less pronounced in pPb than in PbPb collisions, and independent of transverse momentum $p_\mathrm{T}^\Upsilon$ and center-of-mass rapidity $y_\mathrm{CM}^\Upsilon$ of the individual $\Upsilon$ state in the studied range $p_\mathrm{T}^\Upsilon$$\lt$ 30 GeV$/c$ and $\vert y_\mathrm{CM}^\Upsilon\vert$$\lt$ 1.93. Models that incorporate sequential suppression of bottomonia in pPb collisions are in better agreement with the data than those which only assume initial-state modifications.

Double-parton scattering is investigated using events with a Z boson and jets. The Z boson is reconstructed using only the dimuon channel. The measurements are performed with proton-proton collision data recorded by the CMS experiment at the LHC at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$ collected in the year 2016. Differential cross sections of Z + $\geq$ 1 jet and Z + $\geq$ 2 jets are measured with transverse momentum of the jets above 20 GeV and pseudorapidity $|\eta|$$\lt$ 2.4. Several distributions with sensitivity to double-parton scattering effects are measured as functions of the angle and the transverse momentum imbalance between the Z boson and the jets. The measured distributions are compared with predictions from several event generators with different hadronization models and different parameter settings for multiparton interactions. The measured distributions show a dependence on the hadronization and multiparton interaction simulation parameters, and are important input for future improvements of the simulations.

Differential cross sections for the Drell-Yan process, including Z boson production, using the dimuon decay channel are measured in proton-lead (pPb) collisions at a nucleon-nucleon centre-of-mass energy of 8.16 TeV. A data sample recorded with the CMS detector at the LHC is used, corresponding to an integrated luminosity of 173 nb$^{-1}$. The differential cross section as a function of the dimuon mass is measured in the range 15-600 GeV, for the first time in proton-nucleus collisions. It is also reported as a function of dimuon rapidity over the mass ranges 15-60 GeV and 60-120 GeV, and ratios for the p-going over the Pb-going beam directions are built. In both mass ranges, the differential cross sections as functions of the dimuon transverse momentum $p_\mathrm{T}$ and of a geometric variable $\phi^*$ are measured, where $\phi^*$ highly correlates with $p_\mathrm{T}$ but is determined with higher precision. In the Z mass region, the rapidity dependence of the data indicate a modification of the distribution of partons within a lead nucleus as compared to the proton case. The data are more precise than predictions based upon current models of parton distributions.

Angular distributions of the decay B$^+$$\to$ K$^*$(892)$^+\mu^+\mu^-$ are studied using events collected with the CMS detector in $\sqrt{s} =$ 8 TeV proton-proton collisions at the LHC, corresponding to an integrated luminosity of 20.0 fb$^{-1}$. The forward-backward asymmetry of the muons and the longitudinal polarization of the K$^*$(892)$^+$ meson are determined as a function of the square of the dimuon invariant mass. These are the first results from this exclusive decay mode and are in agreement with a standard model prediction.

The second-order Fourier coefficients ($v_2$) characterizing the azimuthal distribution of $\Upsilon$(1S) and $\Upsilon$(2S) mesons arising from PbPb collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV are studied. The $\Upsilon$ mesons are reconstructed in their dimuon decay channel, as measured by the CMS detector. The data set corresponds to an integrated luminosity of 1.7 nb$^{-1}$. The scalar product method is used to extract the $v_2$ coefficients of the azimuthal distribution. Results are reported for the rapidity range $|y|$$\lt$ 2.4, with the transverse momentum 0 $\lt$$p_\mathrm{T}$$\lt$ 50 GeV/$c$, and in three centrality ranges of 10-30%, 30-50% and 50-90%. In contrast to the J/$\psi$ mesons, the measured $v_2$ values for the $\Upsilon$ mesons are found to be consistent with zero.

The cross section of top quark pair production is measured in the $\mathrm{t\bar{t}}\to (\ell\nu_{\ell})(\tau_\mathrm{h}\nu_{\tau})\mathrm{b\bar{b}}$ final state, where $\tau_\mathrm{h}$ refers to the hadronic decays of the $\tau$ lepton, and $\ell$ is either an electron or a muon. The data sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$ collected in proton-proton collisions at $\sqrt{s}=$ 13 TeV with the CMS detector. The measured cross section is $\sigma_{\mathrm{t\bar{t}}} =$ 781 $\pm$ 7 (stat) $\pm$ 62 (syst) $\pm$ 20 (lum) pb, and the ratio of the partial width $\Gamma($t$\to\tau\nu_{\tau}$b) to the total decay width of the top quark is measured to be 0.1050 $\pm$ 0.0009 (stat) $\pm$ 0.0071 (syst). This is the first measurement of the $\mathrm{t\bar{t}}$ production cross section in proton-proton collisions at $\sqrt{s}=$ 13 TeV that explicitly includes $\tau$ leptons. The ratio of the cross sections in the $\ell\tau_\mathrm{h}$ and $\ell\ell$ final states yields a value $R_{\ell\tau_\mathrm{h}/\ell\ell}=$ 0.973 $\pm$ 0.009 (stat) $\pm$ 0.066 (syst), consistent with lepton universality.

A measurement of the production cross section of top quark pairs in association with two b jets ($\mathrm{t\bar{t}}\mathrm{b\bar{b}}$) is presented using data collected in proton-proton collisions at $\sqrt{s} =$ 13 TeV by the CMS detector at the LHC corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The cross section is measured in the all-jet decay channel of the top quark pair by selecting events containing at least eight jets, of which at least two are identified as originating from the hadronization of b quarks. A combination of multivariate analysis techniques is used to reduce the large background from multijet events not containing a top quark pair, and to help discriminate between jets originating from top quark decays and other additional jets. The cross section is determined for the total phase space to be 5.5 $\pm$ 0.3 (stat)${}^{+1.6}_{-1.3}$ (syst) pb and also measured for two fiducial $\mathrm{t\bar{t}}\mathrm{b\bar{b}}$ definitions. The measured cross sections are found to be larger than theoretical predictions by a factor of 1.5-2.4, corresponding to 1-2 standard deviations.

Two related searches for phenomena beyond the standard model (BSM) are performed using events with hadronic jets and significant transverse momentum imbalance. The results are based on a sample of proton-proton collisions at a center-of-mass energy of 13 TeV, collected by the CMS experiment at the LHC in 2016-2018 and corresponding to an integrated luminosity of 137 fb$^{-1}$. The first search is inclusive, based on signal regions defined by the hadronic energy in the event, the jet multiplicity, the number of jets identified as originating from bottom quarks, and the value of the kinematic variable $M_\mathrm{T2}$ for events with at least two jets. For events with exactly one jet, the transverse momentum of the jet is used instead. The second search looks in addition for disappearing tracks produced by BSM long-lived charged particles that decay within the volume of the tracking detector. No excess event yield is observed above the predicted standard model background. This is used to constrain a range of BSM models that predict the following: the pair production of gluinos and squarks in the context of supersymmetry models conserving $R$-parity, with or without intermediate long-lived charginos produced in the decay chain; the resonant production of a colored scalar state decaying to a massive Dirac fermion and a quark; or the pair production of scalar and vector leptoquarks each decaying to a neutrino and a top, bottom, or light-flavor quark. In most of the cases, the results obtained are the most stringent constraints to date.

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.

The production of W$^\pm$ bosons is studied in proton-lead (pPb) collisions at a nucleon-nucleon centre-of-mass energy of $\sqrt{s_\mathrm{NN}} =$ 8.16 TeV. Measurements are performed in the W$^\pm$ $\to$ $\mu^\pm\nu_\mu$ channel using a data sample corresponding to an integrated luminosity of 173.4 $\pm$ 8.7 nb$^{-1}$, collected by the CMS Collaboration at the LHC. The number of positively and negatively charged W bosons is determined separately in the muon pseudorapidity region in the laboratory frame $|\eta^\mu_\mathrm{lab}|$ $<$ 2.4 and transverse momentum $p_\mathrm{T}^\mu$ $>$ 25 GeV/$c$. The W$^\pm$ boson differential cross sections, muon charge asymmetry, and the ratios of W$^\pm$ boson yields for the proton-going over the Pb-going beam directions are reported as a function of the muon pseudorapidity in the nucleon-nucleon centre-of-mass frame. The measurements are compared to the predictions from theoretical calculations based on parton distribution functions (PDFs) at next-to-leading-order. The results favour PDF calculations that include nuclear modifications and provide constraints on the nuclear PDF global fits.

Measurements of the differential cross section for the Drell-Yan process, based on proton-proton collision data at a centre-of-mass energy of 13 TeV, collected by the CMS experiment, are presented. The data correspond to an integrated luminosity of 2.8 (2.3) fb$^{-1}$ in the dimuon (dielectron) channel. The total and fiducial cross section measurements are presented as a function of dilepton invariant mass in the range 15 to 3000 GeV, and compared with the perturbative predictions of the standard model. The measured differential cross sections are in good agreement with the theoretical calculations.

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

Three of the most significant measured deviations from standard model predictions, the enhanced decay rate for B $\to$ D$^{(*)}\tau\nu$, hints of lepton universality violation in B $\to$ K$^{(*)}\ell\ell$ decays, and the anomalous magnetic moment of the muon, can be explained by the existence of leptoquarks (LQs) with large couplings to third-generation quarks and masses at the TeV scale. The existence of these states can be probed at the LHC in high energy proton-proton collisions. A novel search is presented for pair production of LQs coupled to a top quark and a muon using data at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$, recorded by the CMS experiment. No deviation from the standard model prediction has been observed and scalar LQs decaying exclusively into t$\mu$ are excluded up to masses of 1420 GeV. The results of this search are combined with those from previous searches for LQ decays into t$\tau$ and b$\nu$, which excluded scalar LQs below masses of 900 and 1080 GeV. Vector LQs are excluded up to masses of 1190 GeV for all possible combinations of branching fractions to t$\mu$, t$\tau$ and b$\nu$. With this analysis, all relevant couplings of LQs with an electric charge of -1/3 to third-generation quarks are probed for the first time.

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.

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

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.

The cross sections for $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) production in lead-lead (PbPb) and proton-proton (pp) collisions at $\sqrt{s_{_\mathrm{NN}}} =$ 5.02 TeV have been measured using the CMS detector at the LHC. The nuclear modification factors, R$_\mathrm{AA}$, derived from the PbPb-to-pp ratio of yields for each state, are studied as functions of meson rapidity and transverse momentum, as well as PbPb collision centrality. The yields of all three states are found to be significantly suppressed, and compatible with a sequential ordering of the suppression, R$_\mathrm{AA}$($\Upsilon$(1S)) $>$ R$_\mathrm{AA}$($\Upsilon$(2S)) $>$ R$_\mathrm{AA}$($\Upsilon$(3S)) . The suppression of $\Upsilon$(1S) is larger than that seen at $\sqrt{s_{_\mathrm{NN}}} =$ 2.76 TeV, although the two are compatible within uncertainties. The upper limit on the R$_\mathrm{AA}$ of $\Upsilon$(3S) integrated over $p_\mathrm{T}$ and rapidity is 0.094 at 95% confidence level, which is the strongest suppression observed for any hadron species in heavy ion collisions to date.

A search for exotic Higgs boson decays to light pseudoscalars in the final state of two muons and two $\tau$ leptons is performed using proton-proton collision data recorded by the CMS experiment at the LHC at a center-of-mass energy of 13 TeV in 2016, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Masses of the pseudoscalar boson between 15.0 and 62.5 GeV are probed, and no significant excess of data is observed above the prediction of the standard model. Upper limits are set on the branching fraction of the Higgs boson to two light pseudoscalar bosons in different types of two-Higgs-doublet models extended with a complex scalar singlet.

The nuclear modification factors of J/$\psi$ and $\psi$(2S) mesons are measured in PbPb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV. The analysis is based on PbPb and pp data samples collected by CMS at the LHC in 2015, corresponding to integrated luminosities of 464 $\mu$b$^{-1}$ and 28 pb$^{-1}$, respectively. The measurements are performed in the dimuon rapidity range of $|y| <$ 2.4 as a function of centrality, rapidity, and transverse momentum (p$_\mathrm{T}$) from p$_\mathrm{T}=$ 3 GeV/$c$ in the most forward region and up to 50 GeV/$c$. Both prompt and nonprompt (coming from b hadron decays) mesons are observed to be increasingly suppressed with centrality, with a magnitude similar to the one observed at $\sqrt{s_{\mathrm{NN}}}=$ 2.76 TeV for the two J/$\psi$ meson components. No dependence on rapidity is observed for either prompt or nonprompt J/$\psi$ mesons. An indication of a lower prompt J/$\psi$ meson suppression at p$_\mathrm{T} >$ 25 GeV/$c$ is seen with respect to that observed at intermediate p$_\mathrm{T}$. The prompt $\psi$(2S) meson yield is found to be more suppressed than that of the prompt J/$\psi$ mesons in the entire p$_\mathrm{T}$ range.

Differential production cross sections of J/$\psi$ and $\psi$(2S) charmonium and $\Upsilon$(nS) (n = 1, 2, 3) bottomonium states are measured in proton-proton collisions at $\sqrt{s} =$ 13 TeV, with data collected by the CMS detector at the LHC, corresponding to an integrated luminosity of 2.3 fb$^{-1}$ for the J/$\psi$ and 2.7 fb$^{-1}$ for the other mesons. The five quarkonium states are reconstructed in the dimuon decay channel, for dimuon rapidity $|y| <$ 1.2. The double-differential cross sections for each state are measured as a function of $y$ and transverse momentum, and compared to theoretical expectations. In addition, ratios are presented of cross sections for prompt $\psi$(2S) to J/$\psi$, $\Upsilon$(2S) to $\Upsilon$(1S), and $\Upsilon$(3S) to $\Upsilon$(1S) production.

Four-lepton production in proton-proton collisions, $\mathrm{pp}\to (\mathrm{Z}/ \gamma^*)(\mathrm{Z}/\gamma^*) \to 4\ell$, where $\ell = \mathrm{e}$ or $\mu$, is studied at a center-of-mass energy of 13 TeV with the CMS detector at the LHC. The data sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The ZZ production cross section, $\sigma(\mathrm{pp} \to \mathrm{Z}\mathrm{Z}) = 17.2 \pm 0.5\text{ (stat) }\pm 0.7\text{ (syst) }\pm 0.4(\mathrm{theo}) \pm 0.4\text{ (lumi)}$ pb, measured using events with two opposite-sign, same-flavor lepton pairs produced in the mass region $60 < m_{\ell^+\ell^-} < $120 GeV, is consistent with standard model predictions. Differential cross sections are measured and are well described by the theoretical predictions. The Z boson branching fraction to four leptons is measured to be $\mathcal{B}(\mathrm{Z}\to 4\ell) = 4.8 \pm 0.2\text{ (stat) }\pm 0.2\text{ (syst) } \pm 0.1\text{ (theo) }\pm 0.1\text{ (lumi) }\times 10^{-6}$ for events with a four-lepton invariant mass in the range 80 $ < m_{4\ell} < $ 100 GeV and a dilepton mass $m_{\ell\ell} > $4 GeV for all opposite-sign, same-flavor lepton pairs. The results agree with standard model predictions. The invariant mass distribution of the four-lepton system is used to set limits on anomalous ZZZ and ZZ$\gamma$ couplings at 95% confidence level: $-0.0012 < f_4^\mathrm{Z} < 0.0010$, $-0.0010 < f_5^\mathrm{Z} < 0.0013$, $-0.0012 < f_4^{\gamma} < 0.0013$, $-0.0012 < f_5^{\gamma} < 0.0013$.

A measurement of the differential cross sections for a W boson produced in association with jets in the muon decay channel is presented. The measurement is based on 13 TeV proton-proton collision data corresponding to an integrated luminosity of 2.2 inverse femtobarns, recorded by the CMS detector at the LHC. The cross sections are reported as functions of jet multiplicity, jet transverse momentum pT, jet rapidity, the scalar pT sum of the jets, and angular correlations between the muon and the jet for different jet multiplicities. The measured cross sections are in agreement with predictions that include multileg leading-order (LO) and next-to-LO matrix element calculations interfaced with parton showers, as well as a next-to-next-to-LO calculation for the W boson and one jet production.

The differential cross section for inclusive particle production as a function of energy in proton-proton collisions at a center-of-mass energy of 13 TeV is measured in the very forward region of the CMS detector. The measurement is based on data collected with the CMS apparatus at the LHC, and corresponds to an integrated luminosity of 0.35 inverse microbarns. The energy is measured in the CASTOR calorimeter, which covers the pseudorapidity region -6.6 < eta < -5.2. The results are given as a function of the total energy deposited in CASTOR, as well as of its electromagnetic and hadronic components. The spectra are sensitive to the modeling of multiparton interactions in pp collisions, and provide new constraints for hadronic interaction models used in collider and in high energy cosmic ray physics.

The cross section of top quark-antiquark pair production in proton-proton collisions at sqrt(s) = 13 TeV is measured by the CMS experiment at the LHC, using data corresponding to an integrated luminosity of 2.2 inverse femtobarns. The measurement is performed by analyzing events in which the final state includes one electron, one muon, and two or more jets, at least one of which is identified as originating from hadronization of a b quark. The measured cross section is 815 +/- 9 (stat) +/- 38 (syst) +/- 19 (lumi) pb, in agreement with the expectation from the standard model.