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.

A measurement of the inclusive cross section of top quark pair production in association with a Z boson using proton-proton collisions at a center-of-mass energy of 13 TeV at the LHC is performed. The data sample corresponds to an integrated luminosity of 77.5 fb$^{-1}$, collected by the CMS experiment during 2016 and 2017. The measurement is performed using final states containing three or four charged leptons (electrons or muons), and the Z boson is detected through its decay to an oppositely charged lepton pair. The production cross section is measured to be $\sigma(\mathrm{t\bar{t}Z})$ $=$ 0.95 $\pm$ 0.05 (stat) $\pm$ 0.06 (syst) pb. For the first time, differential cross sections are measured as functions of the transverse momentum of the Z boson and the angular distribution of the negatively charged lepton from the Z boson decay. The most stringent direct limits to date on the anomalous couplings of the top quark to the Z boson are presented, including constraints on the Wilson coefficients in the framework of the standard model effective field theory.

A search for exclusive two-photon production via photon exchange in proton-proton collisions, pp $\to$ p$\gamma\gamma$p with intact protons, is presented. The data correspond to an integrated luminosity of 9.4 fb$^{-1}$ collected in 2016 using the CMS and TOTEM detectors at a center-of-mass energy of 13 TeV at the LHC. Events are selected with a diphoton invariant mass above 350 GeV and with both protons intact in the final state, to reduce backgrounds from strong interactions. The events of interest are those where the invariant mass and rapidity calculated from the momentum losses of the forward-moving protons matches the mass and rapidity of the central, two-photon system. No events are found that satisfy this condition. Interpreting this result in an effective dimension-8 extension of the standard model, the first limits are set on the two anomalous four-photon coupling parameters. If the other parameter is constrained to its standard model value, the limits at 95% CL are $\lvert\zeta_1\rvert$ $\lt$ 2.9 $\times$ 10$^{-13}$ GeV$^{-4}$ and $\lvert\zeta_2\rvert$ $\lt$ 6.0 $\times$ 10$^{-13}$ GeV$^{-4}$.

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.

Many measurements at the LHC require efficient identification of heavy-flavour jets, i.e. jets originating from bottom (b) or charm (c) quarks. An overview of the algorithms used to identify c jets is described and a novel method to calibrate them is presented. This new method adjusts the entire distributions of the outputs obtained when the algorithms are applied to jets of different flavours. It is based on an iterative approach exploiting three distinct control regions that are enriched with either b jets, c jets, or light-flavour and gluon jets. Results are presented in the form of correction factors evaluated using proton-proton collision data with an integrated luminosity of 41.5 fb$^{-1}$ at $\sqrt{s}$ = 13 TeV, collected by the CMS experiment in 2017. The closure of the method is tested by applying the measured correction factors on simulated data sets and checking the agreement between the adjusted simulation and collision data. Furthermore, a validation is performed by testing the method on pseudodata, which emulate different miscalibration conditions. The calibrated results enable the use of the full distributions of heavy-flavour identification algorithm outputs, e.g. as inputs to machine-learning models. Thus, they are expected to increase the sensitivity of future physics analyses.

Multiparticle azimuthal correlations of prompt D$^0$ mesons are measured in PbPb collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV. For the first time, a four-particle cumulant method is used to extract the second Fourier coefficient of the azimuthal distribution ($v_2$) of D$^0$ mesons as a function of event centrality and the D$^0$ transverse momentum. The ratios of the four-particle $v_2$ values to previously measured two-particle cumulant results provide direct experimental access to event-by-event fluctuations of charm quark azimuthal anisotropies. These ratios are also found to be comparable to those of inclusive charged particles in the event. However, hints of deviations are seen in the most central and peripheral collisions. To investigate the origin of flow fluctuations in the charm sector, these measurements are compared with models implementing fluctuations of charm quark energy loss via collisional or radiative processes in the quark-gluon plasma. These models cannot quantitatively describe the data over the full transverse momentum and centrality ranges, although the calculations with collisional energy loss provide a better description of the data.

We search for new massive scalar particles X and Y through the resonant process X $\to$ YH $\to$$\mathrm{b\bar{b}b\bar{b}}$, where H is the standard model Higgs boson. Data from CERN LHC proton-proton collisions are used, collected at a centre-of-mass energy of 13 TeV in 2016-2018 and corresponding to an integrated luminosity of 138 fb$^{-1}$. The search is performed in mass ranges of 0.9-4 TeV for X and 60-600 GeV for Y, where both Y and H are reconstructed as Lorentz-boosted single large-area jets. The results are interpreted in the context of the next-to-minimal supersymmetric standard model and also in an extension of the standard model with two additional singlet scalar fields. The 95% confidence level upper limits for the production cross section vary between 0.1 and 150 fb depending on the X and Y masses, and represent a significant improvement over results from previous searches.

A new algorithm is presented to discriminate reconstructed hadronic decays of tau leptons ($\tau_\mathrm{h}$) that originate from genuine tau leptons in the CMS detector against $\tau_\mathrm{h}$ candidates that originate from quark or gluon jets, electrons, or muons. The algorithm inputs information from all reconstructed particles in the vicinity of a $\tau_\mathrm{h}$ candidate and employs a deep neural network with convolutional layers to efficiently process the inputs. This algorithm leads to a significantly improved performance compared with the previously used one. For example, the efficiency for a genuine $\tau_\mathrm{h}$ to pass the discriminator against jets increases by 10-30% for a given efficiency for quark and gluon jets. Furthermore, a more efficient $\tau_\mathrm{h}$ reconstruction is introduced that incorporates additional hadronic decay modes. The superior performance of the new algorithm to discriminate against jets, electrons, and muons and the improved $\tau_\mathrm{h}$ reconstruction method are validated with LHC proton-proton collision data at $\sqrt{s} =$ 13 TeV.

The first study of the shapes of jets arising from bottom (b) quarks in heavy ion collisions is presented. Jet shapes are studied using charged hadron constituents as a function of their radial distance from the jet axis. Lead-lead (PbPb) collision data at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV were recorded by the CMS detector at the LHC, with an integrated luminosity of 1.69 nb$^{-1}$. Compared to proton-proton collisions, a redistribution of the energy in b jets to larger distances from the jet axis is observed in PbPb collisions. This medium-induced redistribution is found to be substantially larger for b jets than for inclusive jets.

We present the first direct search for exotic Higgs boson decays H $\to$$\mathcal{A}\mathcal{A}$, $\mathcal{A}$$\to$$\gamma\gamma$ in events with two photonlike objects. The hypothetical particle $\mathcal{A}$ is a low-mass spin-0 particle decaying promptly to a merged diphoton reconstructed as a single photonlike object. We analyze the data collected by the CMS experiment at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 136 fb$^{-1}$. No excess above the estimated background is found. We set upper limits on the branching fraction $\mathcal{B}$(H $\to$$\mathcal{A}\mathcal{A}$$\to$ 4$\gamma$) of (0.9-3.3) $\times$ 10$^{-3}$ at 95% confidence level for masses of $\mathcal{A}$ in the range 0.1-1.2 GeV.