A search for new massive particles decaying into a pair of Higgs bosons in proton-proton collisions at a center-of-mass energy of 13 TeV is presented. Data were collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The search is performed for resonances with a mass between 0.8 and 3.5 TeV using events in which one Higgs boson decays into a bottom quark pair and the other decays into two W bosons that subsequently decay into a lepton, a neutrino, and a quark pair. The Higgs boson decays are reconstructed with techniques that identify final state quarks as substructure within boosted jets. The data are consistent with standard model expectations. Exclusion limits are placed on the product of the cross section and branching fraction for generic spin-0 and spin-2 massive resonances. The results are interpreted in the context of radion and bulk graviton production in models with a warped extra spatial dimension. These are the best results to date from searches for an HH resonance decaying to this final state, and they are comparable to the results from searches in other channels for resonances with masses below 1.5 TeV.

A statistical combination of searches for heavy resonances decaying to pairs of bosons or leptons is presented. The data correspond to an integrated luminosity of 35.9 fb$^{-1}$ collected during 2016 by the CMS experiment at the LHC in proton-proton collisions at a center-of-mass energy of 13 TeV. The data are found to be consistent with expectations from the standard model background. Exclusion limits are set in the context of models of spin-1 heavy vector triplets and of spin-2 bulk gravitons. For mass-degenerate W' and Z' resonances that predominantly couple to the standard model gauge bosons, the mass exclusion at 95% confidence level of heavy vector bosons is extended to 4.5 TeV as compared to 3.8 TeV determined from the best individual channel. This excluded mass increases to 5.0 TeV if the resonances couple predominantly to fermions.

A search is presented for a charged Higgs boson heavier than the top quark, produced in association with a top quark, or with a top and a bottom quark, and decaying into a top-bottom quark-antiquark pair. The search is performed using proton-proton collision data collected by the CMS experiment at the LHC at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Events are selected by the presence of a single isolated charged lepton (electron or muon) or an opposite-sign dilepton (electron or muon) pair, categorized according to the jet multiplicity and the number of jets identified as originating from b quarks. Multivariate analysis techniques are used to enhance the discrimination between signal and background in each category. The data are compatible with the standard model, and 95% confidence level upper limits of 9.6-0.01 pb are set on the charged Higgs boson production cross section times branching fraction to a top-bottom quark-antiquark pair, for charged Higgs boson mass hypotheses ranging from 200 GeV to 3 TeV. The upper limits are interpreted in different minimal supersymmetric extensions of the standard model.

A search for decays of the Higgs and Z boson to pairs of J/$\psi$ or $\Upsilon$(nS) (n=1, 2, 3) mesons, with their subsequent decay to $\mu^+\mu^-$ pairs, is presented. The analysis uses data from proton-proton collisions at $\sqrt{s}=$ 13 TeV, collected with the CMS detector at the LHC in 2017 and corresponding to an integrated luminosity of 37.5 fb$^{-1}$. While an observation of such a decay with this sample would indicate the presence of physics beyond the standard model, no significant excess is observed. Upper limits at 95% confidence level are placed on the branching fractions of these decays. In the J/$\psi$ pair channel, the limits are 1.8$\times$10$^{-3}$ and 2.2$\times$10$^{-6}$ for the Higgs and Z boson, respectively, while in the combined $\Upsilon$(nS) pair channel, the limits are 1.4$\times$ 10$^{-3}$ and 1.5$\times$10$^{-6}$, respectively, when the mesons from the Higgs and Z boson decay are assumed to be unpolarized. When fully longitudinal and transverse polarizations are considered the limits reduce by about 22-29% and increase by about 10-13%, respectively.

A measurement is presented of differential cross sections for $t$-channel single top quark and antiquark production in proton-proton collisions at a centre-of-mass energy of 13 TeV by the CMS experiment at the LHC. From a data set corresponding to an integrated luminosity of 35.9 fb$^{-1}$, events containing one muon or electron and two or three jets are analysed. The cross section is measured as a function of the top quark transverse momentum ($p_\mathrm{T}$), rapidity, and polarisation angle, the charged lepton $p_\mathrm{T}$ and rapidity, and the $p_\mathrm{T}$ of the W boson from the top quark decay. In addition, the charge ratio is measured differentially as a function of the top quark, charged lepton, and W boson kinematic observables. The results are found to be in agreement with standard model predictions using various next-to-leading-order event generators and sets of parton distribution functions. Additionally, the spin asymmetry, sensitive to the top quark polarisation, is determined from the differential distribution of the polarisation angle at parton level to be 0.440 $\pm$ 0.070, in agreement with the standard model prediction.

A search is presented for three additional operators that would lead to anomalous WW$\gamma$ or WWZ couplings with respect to those in the standard model. They are constrained by studying events with two vector bosons; a W boson decaying to e$\nu$ or $\mu\nu$, and a W or Z boson decaying hadronically, reconstructed as a single, massive, large-radius jet. The search uses a data set of proton-proton collisions at a centre-of-mass energy of 13 TeV, recorded by the CMS experiment at the CERN LHC in 2016, and corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Using the reconstructed diboson invariant mass, 95% confidence intervals are obtained for the anomalous coupling parameters of $-1.58< c_\mathrm{WWW}/\Lambda^2$ $<$ 1.59 TeV$^{-2}$, $-$2.00 $<$ $c_\mathrm{W}/\Lambda^2$ $<$ 2.65 TeV$^{-2}$, and $-$8.78 $<$ $c_\mathrm{B}/\Lambda^2$ $<$ 8.54 TeV$^{-2}$, in agreement with standard model expectations of zero for each parameter. These are the strictest bounds on these parameters to date.

A measurement is presented of the cross section for electroweak production of a Z boson and a photon in association with two jets (Z$\gamma$jj) in proton-proton collisions. The Z boson candidates are selected through their decay into a pair of electrons or muons. The process of interest, electroweak Z$\gamma$jj production, is isolated by selecting events with a large dijet mass and a large pseudorapidity gap between the two jets. The measurement is based on data collected at the CMS experiment at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The observed significance of the signal is 3.9 standard deviations, where a significance of 5.2 standard deviations is expected in the standard model. These results are combined with published results by CMS at $\sqrt{s} =$ 8 TeV, which leads to observed and expected respective significances of 4.7 and 5.5 standard deviations. From the 13 TeV data, a value is obtained for the signal strength of electroweak Z$\gamma$jj production and bounds are given on quartic vector boson interactions in the framework of dimension-eight effective field theory operators.

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.

Central exclusive and semiexclusive production of $\pi^+\pi^-$ pairs is measured with the CMS detector in proton-proton collisions at the LHC at center-of-mass energies of 5.02 and 13 TeV. The theoretical description of these nonperturbative processes, which have not yet been measured in detail at the LHC, poses a significant challenge to models. The two pions are measured and identified in the CMS silicon tracker based on specific energy loss, whereas the absence of other particles is ensured by calorimeter information. The total and differential cross sections of exclusive and semiexclusive central $\pi^+\pi^-$ production are measured as functions of invariant mass, transverse momentum, and rapidity of the $\pi^+\pi^-$ system in the fiducial region defined as transverse momentum $p_\mathrm{T}(\pi)$ $>$ 0.2 GeV and pseudorapidity $|\eta(\pi)|$ $<$ 2.4. The production cross sections for the four resonant channels f$_0(500)$, $\rho^0(770)$, f$_0(980)$, and f$_2(1270)$ are extracted using a simple model. These results represent the first measurement of this process at the LHC collision energies of 5.02 and 13 TeV.

Measurements are presented of the differential cross sections for Z bosons produced in proton-proton collisions at $\sqrt{s} =$ 13 TeV and decaying to muons and electrons. The data analyzed were collected in 2016 with the CMS detector at the LHC and correspond to an integrated luminosity of 35.9 fb$^{-1}$. The measured fiducial inclusive product of cross section and branching fraction agrees with next-to-next-to-leading order quantum chromodynamics calculations. Differential cross sections of the transverse momentum $p_\mathrm{T}$, the optimized angular variable $\phi^*_\eta$, and the rapidity of lepton pairs are measured. The data are corrected for detector effects and compared to theoretical predictions using fixed order, resummed, and parton shower calculations. The uncertainties of the measured normalized cross sections are smaller than 0.5% for $\phi^*_\eta <$ 0.5 and for $p_\mathrm{T}^\mathrm{Z} <$ 50 GeV.

The differential cross sections of $\Lambda_\mathrm{c}^+$ baryon production are measured via the exclusive decay channel $\Lambda_\mathrm{c}^+ \to $pK$^-\pi^+$, as a function of transverse momentum ($p_\mathrm{T}$) in proton-proton (pp) and lead-lead (PbPb) collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV with the CMS detector at the LHC. The measurement is performed within the $\Lambda_\mathrm{c}^+$ rapidity interval $|y|<$1.0 in the $p_\mathrm{T}$ range of 5-20 GeV/$c$ in pp and 10-20 GeV/$c$ in PbPb collisions. The observed yields of $\Lambda_\mathrm{c}^+$ for $p_\mathrm{T}$ of 10-20 GeV/$c$ suggest a possible suppression in central PbPb collisions compared to pp collisions. The $\Lambda_\mathrm{c}^+/$D$^0$ production ratio in pp collisions is compared to theoretical models. In PbPb collisions, this ratio is consistent with the result from pp collisions in their common $p_\mathrm{T}$ range.

Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in PbPb collisions at nucleon-nucleon center-of-mass energies $\sqrt{s_\mathrm{NN}} =$ 2.76 and 5.02 TeV with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.

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.

The fiducial cross section for $\Upsilon$(1S) pair production in proton-proton collisions at a center-of-mass energy of 13 TeV in the region where both $\Upsilon$(1S) mesons have an absolute rapidity below 2.0 is measured to be 79 $\pm$ 11 (stat) $\pm$ 6 (syst) $\pm$ 3 ($\mathcal{B}$) pb assuming the mesons are produced unpolarized. The last uncertainty corresponds to the uncertainty in the $\Upsilon$(1S) meson dimuon branching fraction. The measurement is performed in the final state with four muons using proton-proton collision data collected in 2016 by the CMS experiment at the LHC, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. This process serves as a standard model reference in a search for narrow resonances decaying to $\Upsilon$(1S)$\mu^+\mu^-$ in the same final state. Such a resonance could indicate the existence of a tetraquark that is a bound state of two b quarks and two $\bar{\mathrm{b}}$ antiquarks. The tetraquark search is performed for masses in the vicinity of four times the bottom quark mass, between 17.5 and 19.5 GeV, while a generic search for other resonances is performed for masses between 16.5 and 27 GeV. No significant excess of events compatible with a narrow resonance is observed in the data. Limits on the production cross section times branching fraction to four muons via an intermediate $\Upsilon$(1S) resonance are set as a function of the resonance mass.

A study of the production of prompt J/$\psi$ mesons contained in jets in proton-proton collisions at $\sqrt{s} =$ 8 TeV is presented. The analysis is based on data corresponding to an integrated luminosity of 19.1 fb$^{-1}$ collected with the CMS detector at the LHC. For events with at least one observed jet, the angular separation between the J/$\psi$ meson and the jet is used to test whether the J/$\psi$ meson is part of the jet. The analysis shows that most prompt J/$\psi$ mesons with energy above 15 GeV and rapidity $|y|<$ 1 are contained in jets with pseudorapidity $|\eta_{\text{jet}}|$ $<$ 1. The differential distributions of the probability to have a J/$\psi$ meson contained in a jet as a function of jet energy for a fixed J/$\psi$ energy fraction are compared to a theoretical model using the fragmenting jet function approach. The data agree best with fragmenting jet function calculations that use a long-distance matrix element parameter set in which prompt J/$\psi$ mesons are predicted to be unpolarized. This technique demonstrates a new way to test predictions for prompt J/$\psi$ production using nonrelativistic quantum chromodynamics.

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.

This paper reports on a search for an extended scalar sector of the standard model, where a new CP-even (odd) boson decays to a Z boson and a lighter CP-odd (even) boson, and the latter further decays to a b quark pair. The Z boson is reconstructed via its decays to electron or muon pairs. The analysed data were recorded in proton-proton collisions at a center-of-mass energy $\sqrt{s} = $ 13 TeV, collected by the CMS experiment at the LHC during 2016, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Data and predictions from the standard model are in agreement within the uncertainties. Upper limits at 95% confidence level are set on the production cross section times branching fraction, with masses of the new bosons up to 1000 GeV. The results are interpreted in the context of the two-Higgs-doublet model.

The momentum-weighted sum of the electric charges of particles inside a jet, known as jet charge, is sensitive to the electric charge of the particle initiating the parton shower. This paper presents jet charge distributions in $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV lead-lead (PbPb) and proton-proton (pp) collisions recorded with the CMS detector at the LHC. These data correspond to integrated luminosities of 404 $\mu$b$^{-1}$ and 27.4 pb$^{-1}$ for PbPb and pp collisions, respectively. Leveraging the sensitivity of the jet charge to fundamental differences in the electric charges of quarks and gluons, the jet charge distributions from simulated events are used as templates to extract the quark- and gluon-like jet fractions from data. The modification of these jet fractions is examined by comparing pp and PbPb data as a function of the overlap of the colliding Pb nuclei (centrality). This measurement tests the color charge dependence of jet energy loss due to interactions with the quark-gluon plasma. No significant modification between different centrality classes and with respect to pp results is observed in the extracted fractions of quark- and gluon-like jet fractions.

The WZ production cross section is measured in proton-proton collisions at a centre-of-mass energy $\sqrt{s} =$ 13 TeV using data collected with the CMS detector, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The inclusive cross section is measured to be $\sigma_{\text{tot}}$(pp $\to$ WZ$)$ = 48.09 $^{+1.00}_{-0.96}$ (stat) $^{+0.44}_{-0.37}$ (theo) $^{+2.39}_{-2.17}$ (syst) $\pm$ 1.39 (lumi) pb, resulting in a total uncertainty of $-$2.78/$+$2.98 pb. Fiducial cross section and ratios of charge-dependent cross section measurements are provided. Differential cross section measurements are also presented with respect to three variables: the Z boson transverse momentum $p_\mathrm{T}$, the leading jet $p_\mathrm{T}$, and the $m$(WZ) variable, defined as the invariant mass of the system composed of the three leptons and the missing transverse momentum. Differential measurements with respect to the W boson $p_\mathrm{T}$, separated by charge, are also shown. Results are consistent with standard model predictions, favouring next-to-next-to-leading-order predictions over those at next-to-leading order. Constraints on anomalous triple gauge couplings are derived via a binned maximum likelihood fit to the $m$(WZ) variable.

A study of excited $\Lambda_\mathrm{b}^0$ baryons is reported, based on a data sample collected in 2016-2018 with the CMS detector at the LHC in proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of up to 140 fb$^{-1}$. The existence of four excited $\Lambda_\mathrm{b}^0$ states: $\Lambda_\mathrm{b}$(5912)$^0$, $\Lambda_\mathrm{b}$(5920)$^0$, $\Lambda_\mathrm{b}$(6146)$^0$, and $\Lambda_\mathrm{b}$(6152)$^0$ in the $\Lambda_\mathrm{b}^0\pi^+\pi^-$ mass spectrum is confirmed, and their masses are measured. The $\Lambda_\mathrm{b}^0\pi^+\pi^-$ mass distribution exhibits a broad excess of events in the region of 6040-6100 MeV, whose origin cannot be discerned with the present data.

A search is presented for the production of a Higgs boson in association with a single top quark, based on data collected in 2016 by the CMS experiment at the LHC at a center-of-mass energy of 13 TeV, which corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The production cross section for this process is highly sensitive to the absolute values of the top quark Yukawa coupling, $y_t$, the Higgs boson coupling to vector bosons, $g_\mathrm{HVV}$, and, uniquely, to their relative sign. Analyses using multilepton signatures, targeting H $\to$ WW, H $\to$ $\tau\tau$, and H $\to$ ZZ decay modes, and signatures with a single lepton and a $\mathrm{b\overline{b}}$ pair, targeting the H $\to$ $\mathrm{b\overline{b}}$ decay, are combined with a reinterpretation of a measurement in the H $\to$ $\gamma\gamma$ channel to constrain $y_\mathrm{t}$. For a standard model-like value of $g_\mathrm{HVV}$, the data favor positive values of $y_\mathrm{t}$ and exclude values of $y_\mathrm{t}$ below about $-$0.9 $y_\mathrm{t}^\mathrm{SM}$.

Measurements are presented of the single-diffractive dijet cross section and the diffractive cross section as a function of the proton fractional momentum loss $\xi$ and the four-momentum transfer squared $t$. Both processes pp$\to$pX and pp$\to$Xp, ie with the proton scattering to either side of the interaction point, are measured, where X includes at least two jets; the results of the two processes are averaged. The analyses are based on data collected simultaneously with the CMS and TOTEM detectors at the LHC in proton-proton collisions at $\sqrt{s}$= 8 TeV during a dedicated run with $\beta^{\ast} =$ 90 m at low instantaneous luminosity and correspond to an integrated luminosity of 37.5 nb$^{-1}$. The single-diffractive dijet cross section $\sigma^\mathrm{pX}_{\mathrm{jj}}$, in the kinematic region $\xi \lt$ 0.1, 0.03 $\lt |$t$| \lt 1$ GeV$^2$, with at least two jets with transverse momentum $p_\mathrm{T} >$ 40 GeV, and pseudorapidity $|\eta| \lt$ 4.4, is 21.7$\pm$0.9 (stat)$^{+3.0}_{-3.3}$ (syst) $\pm$ 0.9 (lum) nb. The ratio of the single-diffractive to inclusive dijet yields, normalised per unit of $\xi$, is presented as a function of $x$, the longitudinal momentum fraction of the proton carried by the struck parton. The ratio in the kinematic region defined above, for $x$ values in the range $-$2.9 $\leq \log_{10} x \leq$$-$1.6, is $R = (\sigma^\mathrm{pX}_{\mathrm{jj}}/\Delta\xi)/\sigma_{\mathrm{jj}} =$ 0.025$\pm$0.001 (stat) $\pm$ 0.003 (syst), where $\sigma^\mathrm{pX}_{\mathrm{jj}}$ and $\sigma_{\mathrm{jj}}$ are the single-diffractive and inclusive dijet cross sections, respectively. The results are compared with predictions from models of diffractive and nondiffractive interactions. Monte Carlo predictions based on the HERA diffractive parton distribution functions agree well with the data when corrected for the effect of soft rescattering between the spectator partons.

We present the first study of charged-hadron production associated with jets originating from b quarks in proton-proton collisions at a center-of-mass energy of 5.02 TeV. The data sample used in this study was collected with the CMS detector at the CERN LHC and corresponds to an integrated luminosity of 27.4 pb$^{-1}$. To characterize the jet substructure, the differential jet shapes, defined as the normalized transverse momentum distribution of charged hadrons as a function of angular distance from the jet axis, are measured for b jets. In addition to the jet shapes, the per-jet yields of charged particles associated with b jets are also quantified, again as a function of the angular distance with respect to the jet axis. Extracted jet shape and particle yield distributions for b jets are compared with results for inclusive jets, as well as with the predictions from the PYTHIA and HERWIG++ event generators.

A search for dark matter produced in association with top quarks in proton-proton collisions at a center-of-mass energy of 13 TeV is presented. The data set used corresponds to an integrated luminosity of 35.9 fb$^{-1}$ recorded with the CMS detector at the LHC. Whereas previous searches for neutral scalar or pseudoscalar mediators considered dark matter production in association with a top quark pair only, this analysis also includes production modes with a single top quark. The results are derived from the combination of multiple selection categories that are defined to target either the single top quark or the top quark pair signature. No significant deviations with respect to the standard model predictions are observed. The results are interpreted in the context of a simplified model in which a scalar or pseudoscalar mediator particle couples to a top quark and subsequently decays into dark matter particles. Scalar and pseudoscalar mediator particles with masses below 290 and 300 GeV, respectively, are excluded at 95 % confidence level, assuming a dark matter particle mass of 1 GeV and mediator couplings to fermions and dark matter particles equal to unity.

The polarizations of promptly produced $\chi_\mathrm{c1}$ and $\chi_\mathrm{c2}$ mesons are studied using data collected by the CMS experiment at the LHC, in proton-proton collisions at $\sqrt{s} = $ 8 TeV. The $\chi_\mathrm{c}$ states are reconstructed via their radiative decays $\chi_\mathrm{c}$ $\to$ $\mathrm{J}/\psi\, \gamma$, with the photons being measured through conversions to e$^+$e$^-$, which allows the two states to be well resolved. The polarizations are measured in the helicity frame, through the analysis of the $\chi_\mathrm{c2}$ to $\chi_\mathrm{c1}$ yield ratio as a function of the polar or azimuthal angle of the positive muon emitted in the $\mathrm{J}/\psi$ $\to$ $\mu^+\mu^-$ decay, in three ranges of $\mathrm{J}/\psi$ transverse momentum. While no differences are seen between the two states in terms of azimuthal decay angle distributions, they are observed to have significantly different polar anisotropies. The measurement favors a scenario where at least one of the two states is strongly polarized along the helicity quantization axis, in agreement with nonrelativistic quantum chromodynamics predictions. This is the first measurement of significantly polarized quarkonia produced at high transverse momentum.