Diboson production in association with jets is studied in the fully leptonic final states, pp $\to$ (Z$\gamma^*$)(Z/$\gamma^*$)+jets $\to$ 2$\ell$2$\ell'$+jets, ($\ell,\ell'$ = e or $\mu$) in proton-proton collisions at a center-of-mass energy of 13 TeV. The data sample corresponds to an integrated luminosity of 138 fb$^{-1}$ collected with the CMS detector at the LHC. Differential distributions and normalized differential cross sections are measured as a function of jet multiplicity, transverse momentum $p_\mathrm{T}$, pseudorapidity $\eta$, invariant mass and $\Delta\eta$ of the highest-$p_\mathrm{T}$ and second-highest-$p_\mathrm{T}$ jets, and as a function of invariant mass of the four-lepton system for events with various jet multiplicities. These differential cross sections are compared with theoretical predictions that mostly agree with the experimental data. However, in a few regions we observe discrepancies between the predicted and measured values. Further improvement of the predictions is required to describe the ZZ+jets production in the whole phase space.
Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system, in the on-shell ZZ region
Differential cross sections normalized to the fiducial cross section as a function of the number of jets with $p_T > 30$ GeV
Differential cross sections normalized to the fiducial cross section as a function of the $p_T$ of the highest-$p_T$ jet
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.
$\mathrm{J}\mspace{-2mu}/\mspace{-2mu}\psi \to \mu\mu\mu\mu$ branching fraction
$\mathcal{B}(\mathrm{J}\mspace{-2mu}/\mspace{-2mu}\psi \to \mu\mu\mu\mu)$ / $\mathcal{B}(\mathrm{J}\mspace{-2mu}/\mspace{-2mu}\psi \to \mu\mu)$ ratio
A study of the anomalous couplings of the Higgs boson to vector bosons, including $CP$-violation effects, has been conducted using its production and decay in the WW channel. This analysis is performed on proton-proton collision data collected with the CMS detector at the CERN LHC during 2016-2018 at a center-of-mass energy of 13 TeV, and corresponds to an integrated luminosity of 138 fb$^{-1}$. The different-flavor dilepton (e$\mu$) final state is analyzed, with dedicated categories targeting gluon fusion, electroweak vector boson fusion, and associated production with a W or Z boson. Kinematic information from associated jets is combined using matrix element techniques to increase the sensitivity to anomalous effects at the production vertex. A simultaneous measurement of four Higgs boson couplings to electroweak vector bosons is performed in the framework of a standard model effective field theory. All measurements are consistent with the expectations for the standard model Higgs boson and constraints are set on the fractional contribution of the anomalous couplings to the Higgs boson production cross section.
Expected profiled likelihood on $f_{a2}$ using Approach 1. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.
Observed profiled likelihood on $f_{a2}$ using Approach 1. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.
Expected profiled likelihood on $f_{\Lambda1}$ using Approach 1. The signal strength modifiers are treated as free parameters. Axis scales are varied to improve the visibility of important features.
A search is presented for baryon number violating interactions in top quark production and decay. The analysis uses data from proton-proton collisions at a center-of-mass energy of 13 TeV, collected with the CMS detector at the LHC with an integrated luminosity of 138 fb$^{-1}$. Candidate events are selected by requiring two oppositely-charged leptons (electrons or muons) and exactly one jet identified as originating from a bottom quark. Multivariate discriminants are used to separate the signal from the background. No significant deviation from the standard model prediction is observed. Upper limits are placed on the strength of baryon number violating couplings. For the first time the production of single top quarks via baryon number violating interactions is studied. This allows the search to set the most stringent constraints to date on the branching fraction of the top quark decay to a lepton, an up-type quark (u or c), and a down-type quark (d, s, or b). The results improve the previous bounds by three to six orders of magnitude based on the fermion flavor combination of the baryon number violating interactions.
The observed upper limits on the branching fractions of the top quark BNV decays are shown with circle and triangle shapes for electron and muon couplings, respectively. The observed limits corresponding to the $C_t$ and $C_s$ coefficients are shown with filled and open markers, respectively. The yellow light (green dark) bands indicate the range within plus or minus one (two) standard deviations bands around the expected limits.
A search for the production of long-lived particles in proton-proton collisions at a center-of-mass energy of 13 TeV at the CERN LHC is presented. The search is based on data collected by the CMS experiment in 2016-2018, corresponding to a total integrated luminosity of 137 fb$^{-1}$. This search is designed to be sensitive to long-lived particles with mean proper decay lengths between 0.1 and 1000 $\mu$m, whose decay products produce a final state with at least one displaced vertex and missing transverse momentum. A machine learning algorithm, which improves the background rejection power by more than an order of magnitude, is applied to improve the sensitivity. The observation is consistent with the standard model background prediction, and the results are used to constrain split supersymmetry (SUSY) and gauge-mediated SUSY breaking models with different gluino mean proper decay lengths and masses. This search is the first CMS search that shows sensitivity to hadronically decaying long-lived particles from signals with mass differences between the gluino and neutralino below 100 GeV. It sets the most stringent limits to date for split-SUSY models and gauge-mediated SUSY breaking models with gluino proper decay length less than 6 $\mu$m.
Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of 3. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1900 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.
Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of 3. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1800 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.
Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of 4. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1900 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.
Using proton-proton collision data corresponding to an integrated luminosity of 140 fb$^{-1}$ collected by the CMS experiment at $\sqrt{s}$ = 13 TeV, the $\Lambda_\text{b}^0$$\to$ J/$\psi\Xi^-$K$^+$ decay is observed for the first time, with a statistical significance exceeding 5 standard deviations. The relative branching fraction, with respect to the $\Lambda_\text{b}^0$$\to$$\psi$(2S)$\Lambda$ decay, is measured to be $\mathcal{B}$($\Lambda_\text{b}^0$$\to$ J/$\psi\Xi^-$K$^+$)/$\mathcal{B}$( $\Lambda_\text{b}^0$$\to$$\psi$(2S)$\Lambda$) = [3.38 $\pm$ 1.02 $\pm$ 0.61 $\pm$ 0.03]%, where the first uncertainty is statistical, the second is systematic, and the third is related to the uncertainties in $\mathcal{B}$($\psi$(2S) $\to$ J/$\psi\pi^+\pi^-$) and $\mathcal{B}$($\Xi^-$ $\to$ $\Lambda\pi^-$).
The measured branching fraction ratio
A test of lepton flavor universality in B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$ and B$^{\pm}$$\to$ K$^{\pm}$e$^+$e$^-$ decays, as well as a measurement of differential and integrated branching fractions of a nonresonant B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$ decay are presented. The analysis is made possible by a dedicated data set of proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded in 2018, by the CMS experiment at the LHC, using a special high-rate data stream designed for collecting about 10 billion unbiased b hadron decays. The ratio of the branching fractions $\mathcal{B}$(B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$) to $\mathcal{B}$(B$^{\pm}$$\to$ K$^{\pm}$e$^+$e$^-$) is determined from the measured double ratio $R$(K) of these decays to the respective branching fractions of the B$^\pm$$\to$ J/$\psi$K$^\pm$ with J/$\psi$$\to$$\mu^+\mu^-$ and e$^+$e$^-$ decays, which allow for significant cancellation of systematic uncertainties. The ratio $R$(K) is measured in the range 1.1 $\lt q^2 \lt$ 6.0 GeV$^2$, where $q$ is the invariant mass of the lepton pair, and is found to be $R$(K) = 0.78$^{+0.47}_{-0.23}$, in agreement with the standard model expectation $R$(K) $\approx$ 1. This measurement is limited by the statistical precision of the electron channel. The integrated branching fraction in the same $q^2$ range, $\mathcal{B}$(B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$) = (12.42 $\pm$ 0.68) $\times$ 10$^{-8}$, is consistent with the present world-average value and has a comparable precision.
The differential $\text{B}^+ \to \text{K}^+\mu^+\mu^-$ branching fraction measured in the individual $q^2$ bins. The uncertainties in the yields are statistical uncertainties from the fit, while the branching fraction uncertainties include both the statistical and systematic components.
Differential branching fraction $d\mathcal{B}/dq^2$, with theoretical predictions obtained with the HEPFiT, SuperIso, Flavio, and EOS packages. The HEPFiT predictions are available only for $q^2 < 8\ \mathrm{GeV}^2$.
Relative uncertainties in the differential branching fraction measurement of $\mathrm{B}^+\to\mathrm{K}^+\mu^+\mu^-$ per $q^2$ bin.
A measurement is presented of the primary Lund jet plane (LJP) density in inclusive jet production in proton-proton collisions. The analysis uses 138 fb$^{-1}$ of data collected by the CMS experiment at $\sqrt{s}$ = 13 TeV. The LJP, a representation of the phase space of emissions inside jets, is constructed using iterative jet declustering. The transverse momentum $k_\mathrm{T}$ and the splitting angle $\Delta R$ of an emission relative to its emitter are measured at each step of the jet declustering process. The average density of emissions as function of $\ln(k_\mathrm{T}$/GeV) and $\ln(R/\Delta R)$ is measured for jets with distance parameters $R$ = 0.4 or 0.8, transverse momentum $p_\mathrm{T} \gt$ 700 GeV, and rapidity $\vert y\vert \lt $ 1.7. The jet substructure is measured using the charged-particle tracks of the jet. The measured distributions, unfolded to the level of stable particles, are compared with theoretical predictions from simulations and with perturbative quantum chromodynamics calculations. Due to the ability of the LJP to factorize physical effects, these measurements can be used to improve different aspects of the physics modeling in event generators.
Primary Lund jet plane density for AK4 jets in a one-dimensional representation with bin indices for MC tuning purposes. The mapping between the bin indices and the physical binning can be imported from the XML file attached to this HepData record using the TUnfoldBinningXML class of ROOT (qualitatively, it corresponds to slicing the Lund plane horizontally from low kT to high kT). All systematic uncertainties are bin-to-bin fully correlated (allowing for sign-changes bin-to-bin), with the exception of the statistical uncertainties from data and MC, for which a separate correlation matrix is provided in this HepData record.
Primary Lund jet plane density for AK8 jets in a one-dimensional representation with bin indices for MC tuning purposes. The mapping between the bin indices and the physical binning can be imported from the XML file attached to this HepData record using the TUnfoldBinningXML class of ROOT (qualitatively, it corresponds to slicing the Lund plane horizontally from low kT to high kT). All systematic uncertainties are bin-to-bin fully correlated (allowing for sign-changes bin-to-bin), with the exception of the statistical uncertainties from data and MC, for which a separate correlation matrix is provided in this HepData record.
A search is performed for charged-lepton flavor violating processes in top quark (t) production and decay. The data were collected by the CMS experiment from proton-proton collisions at a center-of-mass energy of 13 TeV and correspond to an integrated luminosity of 138 fb$^{-1}$. The selected events are required to contain one opposite-sign electron-muon pair, a third charged lepton (electron or muon), and at least one jet of which no more than one is associated with a bottom quark. Boosted decision trees are used to distinguish signal from background, exploiting differences in the kinematics of the final states particles. The data are consistent with the standard model expectation. Upper limits at 95% confidence level are placed in the context of effective field theory on the Wilson coefficients, which range between 0.024-0.424 TeV$^{-2}$ depending on the flavor of the associated light quark and the Lorentz structure of the interaction. These limits are converted to upper limits on branching fractions involving up (charm) quarks, t$\to$e$\mu$u (t$\to$e$\mu$c), of 0.032 (0.498)$\times$10$^{-6}$, 0.022 (0.369)$\times$10$^{-6}$, and 0.012 (0.216)$\times$10$^{-6}$ for tensor-like, vector-like, and scalar-like interactions, respectively.
The expected and observed upper limits on CLFV Wilson coefficients. The Limits on the Wilson coefficients are extracted from the upper limits on the cross sections.
The expected and observed upper limits on top quark CLFV branching fractions. The Limits on the top quark CLFV branching fractions are extracted from the upper limits on the Wilson coefficients.
The first search for narrow resonances decaying to three well-separated hadronic jets is presented. The search uses proton-proton collision data corresponding to an integrated luminosity of 138 fb$^{-1}$ at $\sqrt{s}$ = 13 TeV, collected at the CERN LHC. No significant deviations from the background predictions are observed between 1.75-9.00 TeV. The results provide the first mass limits on a right-handed boson Z$_{\mathrm{R}}$ decaying to three gluons, an excited quark decaying via a vector boson to three quarks, as well as updated limits on a Kaluza-Klein gluon decaying via a radion to three gluons.
Observed and expected (background-only fitted) invariant mass spectra of trijet events. Data spectra from three years are fitted separately and the sum is shown in the figure. The fitting function used is ${ d N}/{ d m} = p_{0}(1-x)^{p_{1}}/x^{\sum_{i=2}^{3} p_{i}\log^{i-2}(x)}$. The fitted parameters are $p_{1} = 7.350, p_{2} = 6.926, p_{3} = 0.388$ for 2016, $p_{1} = 8.308, p_{2} = 5.931, p_{3} = 0.167$ for 2017 and $p_{1} = 8.770, p_{2} = 5.617, p_{3} = 0.106$ for 2018. $p_{0}$ is the normalization parameter and its exact value is irrelevant.
Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to ggg) \mathcal{A}$ for a 3-body decay trijet resonance with $\Gamma_{X}\sim 3\% m_{X}$. The acceptance $\mathcal{A}$ is defined as $\mathcal{A} = N$(events with $m_{X}^{GEN} > 85\% m_{X}^{input}$) / $N$(events generated in the full phase space defined by the CMS default generator settings). Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions assuming SM-like couplings are depicted with the red curve.
Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to ggg) \mathcal{A}$ for a 3-body decay trijet resonance with $\Gamma_{X}\sim 0.01\% m_{X}$. The acceptance $\mathcal{A}$ is defined as $\mathcal{A} = N$(events with $m_{X}^{GEN} > 85\% m_{X}^{input}$) / $N$(events generated in the full phase space defined by the CMS default generator settings). Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions are depicted with the red curve.
Measurements of inclusive and normalized differential cross sections of the associated production of top quark-antiquark and bottom quark-antiquark pairs, ttbb, are presented. The results are based on data from proton-proton collisions collected by the CMS detector at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The cross sections are measured in the lepton+jets decay channel of the top quark pair, using events containing exactly one isolated electron or muon and at least five jets. Measurements are made in four fiducial phase space regions, targeting different aspects of the ttbb process. Distributions are unfolded to the particle level through maximum likelihood fits, and compared with predictions from several event generators. The inclusive cross section measurements of this process in the fiducial phase space regions are the most precise to date. In most cases, the measured inclusive cross sections exceed the predictions with the chosen generator settings. The only exception is when using a particular choice of dynamic renormalization scale, $\mu_\mathrm{R}=\frac{1}{2} \prod_{i=\mathrm{t, \bar{t}, b, \bar{b}}} m_{\mathrm{T},i}^{1/4}$, where $m_{\mathrm{T},i}^2=m_i^2+p^2_{\mathrm{T},i}$ are the transverse masses of top and bottom quarks. The differential cross sections show varying degrees of compatibility with the theoretical predictions, and none of the tested generators with the chosen settings simultaneously describe all the measured distributions.
Fiducial cross sections from the measurements of all observables, compared to predictions from different ttbb simulation approaches. For each of the normalized differential measurements the fiducial cross section in the respective phase space is also determined. In the paper only one representative observable is quoted for each fiducial phase space, while here the measured cross section with the uncertainties from the fit to the respective observable is summarized.
Compatibility of normalized differential cross section measurements with modeling predictions. The compatibility is quantified with z scores for each of the theoretical predictions, given the unfolded normalized differential cross sections and their covariances. A lower value indicates a better agreement between prediction and measurement. A value of z = 2 indicates a p-value of 5%. In the calculation of the z score only the measurement uncertainties and the statistical uncertainties of the modeling predictions are taken into account
Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.
Measurements of the charge-dependent two-particle angular correlation function in proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV and lead-lead (PbPb) collisions at$\sqrt{s_\mathrm{NN}}$ = 5.02 TeV are reported. The pPb and PbPb datasets correspond to integrated luminosities of 186\nbinv and 0.607 nb$^{-1}$, respectively, and were collected using the CMS detector at the CERN LHC. The charge-dependent correlations are characterized by balance functions of same- and opposite-sign particle pairs. The balance functions, which contain information about the creation time of charged particle pairs and the development of collectivity, are studied as functions of relative pseudorapidity ($\Delta \eta$) and relative azimuthal angle ($\Delta \phi$), for various multiplicity and transverse momentum ($p_\mathrm{T}$) intervals. A multiplicity dependence of the balance function is observed in $\Delta \eta$ and $\Delta \phi$ for both systems. The width of the balance functions decreases towards high-multiplicity collisions in the momentum region $\lt$2 GeV, for pPb and PbPb results. No multiplicity dependence is observed at higher transverse momentum. The data are compared with HYDJET, HIJING and AMPT generator predictions, none of which capture completely the multiplicity dependence seen in the data.
Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 30-40% centrality
Quasireal photons exchanged in relativistic heavy ion interactions are powerful probes of the gluonic structure of nuclei. The coherent J/$\psi$ photoproduction cross section in ultraperipheral lead-lead collisions is measured as a function of photon-nucleus center-of-mass energies per nucleon (W$^\text{Pb}_{\gamma\text{N}}$), over a wide range of 40 $\lt$ W$^\text{Pb}_{\gamma\text{N}}$$\lt$ 400 GeV. Results are obtained using data at the nucleon-nucleon center-of-mass energy of 5.02 TeV collected by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of 1.52 nb$^{-1}$. The cross section is observed to rise rapidly at low W$^\text{Pb}_{\gamma\text{N}}$, and plateau above W$^\text{Pb}_{\gamma\text{N}}$$\approx$ 40 GeV, up to 400 GeV, a new regime of small Bjorken-$x$ ($\approx$ 6 $\times$ 10$^{-5}$) gluons being probed in a heavy nucleus. The observed energy dependence is not predicted by current quantum chromodynamic models.
The total coherent $\mathrm{J}/\psi$ photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$, measured in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$ values used correspond to the center of each rapidity range. The theoretical uncertainties is due to the uncertainties in the photon flux.
The total coherent $\mathrm{J}/\psi$ photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$, measured in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$ values used correspond to the center of each rapidity range. The theoretical uncertainties is due to the uncertainties in the photon flux.
The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.
For the first time at LHC energies, the forward rapidity gap spectra from proton-lead collisions for both proton and lead dissociation processes are presented. The analysis is performed over 10.4 units of pseudorapidity at a center-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV, almost 300 times higher than in previous measurements of diffractive production in proton-nucleus collisions. For lead dissociation processes, which correspond to the pomeron-lead event topology, the EPOS-LHC generator predictions are a factor of two below the data, but the model gives a reasonable description of the rapidity gap spectrum shape. For the pomeron-proton topology, the EPOS-LHC, QGSJET II, and HIJING predictions are all at least a factor of five lower than the data. The latter effect might be explained by a significant contribution of ultra-peripheral photoproduction events mimicking the signature of diffractive processes. These data may be of significant help in understanding the high energy limit of quantum chromodynamics and for modeling cosmic ray air showers.
The total energy of all Particle Flow candidatesin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV
Unfolded diffraction enhanced differential cross section for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, defined on the region $|\eta| < 5.2$, compared to hadron level predictions of the MC generators. The data are corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap. The corrections are obtained using the EPOS-LHC MC sample. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles
Unfolded diffraction enhanced differential cross section for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, defined on the region $|\eta| < 5.2$, compared to hadron level predictions of the MC generators. The data are corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap. The corrections are obtained using the EPOS-LHC MC sample. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles
A search for long-lived particles decaying in the outer regions of the CMS silicon tracker or in the calorimeters is presented. The search is based on a data sample of proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded with the CMS detector at the LHC in 2016-2018, corresponding to an integrated luminosity of 138 fb$^{-1}$. A novel technique, using trackless and out-of-time jet information combined in a deep neural network discriminator, is employed to identify decays of long-lived particles. The results are interpreted in a simplified model of chargino-neutralino production, where the neutralino is the next-to-lightest supersymmetric particle, is long-lived, and decays to a gravitino and either a Higgs or Z boson. This search is most sensitive to neutralino proper decay lengths of approximately 0.5 m, for which masses up to 1.18 TeV are excluded at 95% confidence level. The current search is the best result to date in the mass range from the kinematic limit imposed by the Higgs mass up to 1.8 TeV.
The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).
The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).
The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).
We report the triton ($t$) production in mid-rapidity ($|y| <$ 0.5) Au+Au collisions at $\sqrt{s_\mathrm{NN}}$= 7.7--200 GeV measured by the STAR experiment from the first phase of the beam energy scan at the Relativistic Heavy Ion Collider (RHIC). The nuclear compound yield ratio ($\mathrm{N}_t \times \mathrm{N}_p/\mathrm{N}_d^2$), which is predicted to be sensitive to the fluctuation of local neutron density, is observed to decrease monotonically with increasing charged-particle multiplicity ($dN_{ch}/d\eta$) and follows a scaling behavior. The $dN_{ch}/d\eta$ dependence of the yield ratio is compared to calculations from coalescence and thermal models. Enhancements in the yield ratios relative to the coalescence baseline are observed in the 0%-10% most central collisions at 19.6 and 27 GeV, with a significance of 2.3$\sigma$ and 3.4$\sigma$, respectively, giving a combined significance of 4.1$\sigma$. The enhancements are not observed in peripheral collisions or model calculations without critical fluctuation, and decreases with a smaller $p_{T}$ acceptance. The physics implications of these results on the QCD phase structure and the production mechanism of light nuclei in heavy-ion collisions are discussed.
Collision energy, centrality, and p_{T} dependence of light nuclei yield, 0%-10% and 40%-80% centrality
Invariant p_{T} spectra of primordial antiprotons in Au+Au collisions at SQRT(s_NN) = 11.5 GeV at 40-80% centrality
Three searches are presented for signatures of physics beyond the standard model (SM) in $\tau\tau$ final states in proton-proton collisions at the LHC, using a data sample collected with the CMS detector at $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. Upper limits at 95% confidence level (CL) are set on the products of the branching fraction for the decay into $\tau$ leptons and the cross sections for the production of a new boson $\phi$, in addition to the H(125) boson, via gluon fusion (gg$\phi$) or in association with b quarks, ranging from $\mathcal{O}$(10 pb) for a mass of 60 GeV to 0.3 fb for a mass of 3.5 TeV each. The data reveal two excesses for gg$\phi$ production with local $p$-values equivalent to about three standard deviations at $m_\phi$ = 0.1 and 1.2 TeV. In a search for $t$-channel exchange of a vector leptoquark U$_1$, 95% CL upper limits are set on the dimensionless U$_1$ leptoquark coupling to quarks and $\tau$ leptons ranging from 1 for a mass of 1 TeV to 6 for a mass of 5 TeV, depending on the scenario. In the interpretations of the $M_\mathrm{h}^{125}$ and $M_\mathrm{h, EFT}^{125}$ minimal supersymmetric SM benchmark scenarios, additional Higgs bosons with masses below 350 GeV are excluded at 95% CL.
Observed and expected distributions of the variable chosen for statistical inference in the low-mass analysis $m_{\tau\tau}$. Numerical values provided in this table correspond to Figure 22 of the auxilliary material of the publication, but restricted to 2016 data-taking year. All distributions are considered after a fit to data is performed using a background-only model, which includes the $\text{H}(125)$ boson. Some details on how the distributions should be used: 1) All given uncertainties correspond to systematic variations of $\pm1\sigma$. 2) Upper values ('plus' in the yaml file) correspond to an upward systematic variation of the parameter ($+1\sigma$). 3) Lower values ('minus' in the yaml file) correspond to a downward systematic variation of the parameter ($-1\sigma$). 4) These variations can have both positive and negative values, depending on the modelled effect. 5) Uncertainties with the same name should be treated as correlated, consistently across the upper and lower variations. 6) Systematic uncertainties with 'prop_' in the name treat limited background statistics per histogram bin, and are deployed with 'Barlow-Beeston-lite' approach. Details in https://arxiv.org/abs/1103.0354 section 5 7) Remaining systematic uncertainties alter the normalization, the shape, or both for a distribution. The nuisance parameter for such an uncertainty is mapped separately on the normalization and the shape variation components of the uncertainty. For normalization, $\ln$ mapping is used, for shape a spline. Details in https://cms-analysis.github.io/HiggsAnalysis-CombinedLimit/part2/settinguptheanalysis/#binned-shape-analysis 8) All nuisance parameters for the systematic uncertainties are modelled with a Gaussian pdf. 9) Gluon fusion contributions are all scaled to 1 pb. Please combine them using either the scale factors from 'Table SM Gluon Fusion Fractions', or using your own composition.
Observed and expected distributions of the variable chosen for statistical inference in the low-mass analysis $m_{\tau\tau}$. Numerical values provided in this table correspond to Figure 22 of the auxilliary material of the publication, but restricted to 2017 data-taking year. All distributions are considered after a fit to data is performed using a background-only model, which includes the $\text{H}(125)$ boson. Some details on how the distributions should be used: 1) All given uncertainties correspond to systematic variations of $\pm1\sigma$. 2) Upper values ('plus' in the yaml file) correspond to an upward systematic variation of the parameter ($+1\sigma$). 3) Lower values ('minus' in the yaml file) correspond to a downward systematic variation of the parameter ($-1\sigma$). 4) These variations can have both positive and negative values, depending on the modelled effect. 5) Uncertainties with the same name should be treated as correlated, consistently across the upper and lower variations. 6) Systematic uncertainties with 'prop_' in the name treat limited background statistics per histogram bin, and are deployed with 'Barlow-Beeston-lite' approach. Details in https://arxiv.org/abs/1103.0354 section 5 7) Remaining systematic uncertainties alter the normalization, the shape, or both for a distribution. The nuisance parameter for such an uncertainty is mapped separately on the normalization and the shape variation components of the uncertainty. For normalization, $\ln$ mapping is used, for shape a spline. Details in https://cms-analysis.github.io/HiggsAnalysis-CombinedLimit/part2/settinguptheanalysis/#binned-shape-analysis 8) All nuisance parameters for the systematic uncertainties are modelled with a Gaussian pdf. 9) Gluon fusion contributions are all scaled to 1 pb. Please combine them using either the scale factors from 'Table SM Gluon Fusion Fractions', or using your own composition.
Observed and expected distributions of the variable chosen for statistical inference in the low-mass analysis $m_{\tau\tau}$. Numerical values provided in this table correspond to Figure 22 of the auxilliary material of the publication, but restricted to 2018 data-taking year. All distributions are considered after a fit to data is performed using a background-only model, which includes the $\text{H}(125)$ boson. Some details on how the distributions should be used: 1) All given uncertainties correspond to systematic variations of $\pm1\sigma$. 2) Upper values ('plus' in the yaml file) correspond to an upward systematic variation of the parameter ($+1\sigma$). 3) Lower values ('minus' in the yaml file) correspond to a downward systematic variation of the parameter ($-1\sigma$). 4) These variations can have both positive and negative values, depending on the modelled effect. 5) Uncertainties with the same name should be treated as correlated, consistently across the upper and lower variations. 6) Systematic uncertainties with 'prop_' in the name treat limited background statistics per histogram bin, and are deployed with 'Barlow-Beeston-lite' approach. Details in https://arxiv.org/abs/1103.0354 section 5 7) Remaining systematic uncertainties alter the normalization, the shape, or both for a distribution. The nuisance parameter for such an uncertainty is mapped separately on the normalization and the shape variation components of the uncertainty. For normalization, $\ln$ mapping is used, for shape a spline. Details in https://cms-analysis.github.io/HiggsAnalysis-CombinedLimit/part2/settinguptheanalysis/#binned-shape-analysis 8) All nuisance parameters for the systematic uncertainties are modelled with a Gaussian pdf. 9) Gluon fusion contributions are all scaled to 1 pb. Please combine them using either the scale factors from 'Table SM Gluon Fusion Fractions', or using your own composition.
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.
Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.
Results are presented from a search for physics beyond the standard model in proton-proton collisions at $\sqrt{s} =$ 13 TeV in channels with two Higgs bosons, each decaying via the process H $\to$$\mathrm{b\bar{b}}$, and large missing transverse momentum. The search uses a data sample corresponding to an integrated luminosity of 137 fb$^{-1}$ collected by the CMS experiment at the CERN LHC. The search is motivated by models of supersymmetry that predict the production of neutralinos, the neutral partners of the electroweak gauge and Higgs bosons. The observed event yields in the signal regions are found to be consistent with the standard model background expectations. The results are interpreted using simplified models of supersymmetry. For the electroweak production of nearly mass-degenerate higgsinos, each of whose decay chains yields a neutralino ($\tilde{\chi}^0_1$) that in turn decays to a massless goldstino and a Higgs boson, $\tilde{\chi}^0_1$ masses in the range 175 to 1025 GeV are excluded at 95% confidence level. For the strong production of gluino pairs decaying via a slightly lighter $\tilde{\chi}^0_2$ to H and a light $\tilde{\chi}^0_1$, gluino masses below 2330 GeV are excluded.
Pre-fit background covariance matrix $\sigma_{xy}$ for the 22 analysis bins, ordered as in Fig. 10.
Pre-fit background correlation matrix $\rho_{xy}$ for the 22 analysis bins, ordered as in Fig. 10.
Efficiency vs $m(\widetilde{\chi}^0_1)$ for SMS model TChiHH-G. The denominator includes all 22 signal regions, and assumes $\mathcal{B}(\mathrm{H}$-->$\mathrm{b}\overline{\mathrm{b}})=100\%$.
A search is presented for single production of a vector-like T quark with charge 2/3 $e$, in the decay channel featuring a top quark and a Z boson, with the top quark decaying hadronically and the Z boson decaying to neutrinos. The search uses data collected by the CMS experiment in proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 137 fb$^{-1}$ recorded at the CERN LHC in 2016-2018. The search is sensitive to a T quark mass between 0.6 and 1.8 TeV with decay widths ranging from negligibly small up to 30% of the T quark mass. Reconstruction strategies for the top quark are based on the degree of Lorentz boosting of its final state. At 95% confidence level, the upper limit on the product of the cross section and branching fraction for a T quark of small decay width varies between 15 and 602 fb, depending on its mass. For a T quark with decay widths between 10 and 30% of its mass, this upper limit ranges between 16 and 836 fb. For most of the studied range, the results provide the best limits to date. This is the first search for single T quark production based on the full Run 2 data set of the LHC.
Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged validation region with 0 forward jets. Data are compared to backgrounds are after the data-driven extraction is performed.
Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged validation region with at least 1 forward jet. Data are compared to backgrounds are after the data-driven extraction is performed.
Distribution of the transverse mass $M_T$ reconstructed from the top quark and the missing transverse energy in the merged validation region with 0 forward jets. Data are compared to backgrounds are after the data-driven extraction is performed.
The chiral magnetic effect (CME) is predicted to occur as a consequence of a local violation of $\cal P$ and $\cal CP$ symmetries of the strong interaction amidst a strong electro-magnetic field generated in relativistic heavy-ion collisions. Experimental manifestation of the CME involves a separation of positively and negatively charged hadrons along the direction of the magnetic field. Previous measurements of the CME-sensitive charge-separation observables remain inconclusive because of large background contributions. In order to better control the influence of signal and backgrounds, the STAR Collaboration performed a blind analysis of a large data sample of approximately 3.8 billion isobar collisions of $^{96}_{44}$Ru+$^{96}_{44}$Ru and $^{96}_{40}$Zr+$^{96}_{40}$Zr at $\sqrt{s_{\rm NN}}=200$ GeV. Prior to the blind analysis, the CME signatures are predefined as a significant excess of the CME-sensitive observables in Ru+Ru collisions over those in Zr+Zr collisions, owing to a larger magnetic field in the former. A precision down to 0.4% is achieved, as anticipated, in the relative magnitudes of the pertinent observables between the two isobar systems. Observed differences in the multiplicity and flow harmonics at the matching centrality indicate that the magnitude of the CME background is different between the two species. No CME signature that satisfies the predefined criteria has been observed in isobar collisions in this blind analysis.
fig4_left_top_isobarpaper_star_purple_solid_v2_subEv_ruru.
A search for long-lived particles (LLPs) produced in decays of standard model (SM) Higgs bosons is presented. The data sample consists of 137 fb$^{-1}$ of proton-proton collisions at $\sqrt{s} =$ 13 TeV, recorded at the LHC in 2016-2018. A novel technique is employed to reconstruct decays of LLPs in the endcap muon detectors. The search is sensitive to a broad range of LLP decay modes and to masses as low as a few GeV. No excess of events above the SM background is observed. The most stringent limits to date on the branching fraction of the Higgs boson to LLPs subsequently decaying to quarks and $\tau^+\tau^-$ are found for proper decay lengths greater than 6, 20, and 40 m, for LLP masses of 7, 15, and 40 GeV, respectively.
The cluster efficiency in bins of hadronic and EM energy in region A. Region A is defined as 391 cm $< r <$ 695.5 cm and 400 cm $< |z| <$ 671 cm. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.
The cluster efficiency in bins of hadronic and EM energy in region A. Region A is defined as 391 cm $< r <$ 695.5 cm and 400 cm $< |z| <$ 671 cm. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.
The cluster efficiency in bins of hadronic and EM energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.
The production cross section of a top quark pair in association with a photon is measured in proton-proton collisions at a center-of-mass energy of 13 TeV. The data set, corresponding to an integrated luminosity of 137 fb$^{-1}$, was recorded by the CMS experiment during the 2016-2018 data taking of the LHC. The measurements are performed in a fiducial volume defined at the particle level. Events with an isolated, highly energetic lepton, at least three jets from the hadronization of quarks, among which at least one is b tagged, and one isolated photon are selected. The inclusive fiducial $\mathrm{t\overline{t}}\gamma$ cross section, for a photon with transverse momentum greater than 20 GeV and pseudorapidity $\lvert \eta\rvert$$\lt$ 1.4442, is measured to be 798 $\pm$ 7 (stat) $\pm$ 48 (syst) fb, in good agreement with the prediction from the standard model at next-to-leading order in quantum chromodynamics. The differential cross sections are also measured as a function of several kinematic observables and interpreted in the framework of the standard model effective field theory (EFT), leading to the most stringent direct limits to date on anomalous electromagnetic dipole moment interactions of the top quark and the photon.
The measured inclusive ttgamma cross section in the fiducial phase space compared to the prediction from simulation using Madgraph_aMC@NLO at a center-of-mass energy of 13 TeV.
The measured inclusive ttgamma cross section in the fiducial phase space compared to the prediction from simulation using Madgraph_aMC@NLO at a center-of-mass energy of 13 TeV.
The first observation of the electroweak (EW) production of a Z boson, a photon, and two forward jets (Z$\gamma$jj) in proton-proton collisions at a center-of-mass energy of 13 TeV is presented. A data set corresponding to an integrated luminosity of 137 fb$^{-1}$, collected by the CMS experiment at the LHC in 2016-2018 is used. The measured fiducial cross section for EW Z$\gamma$jj is $\sigma_{\mathrm{EW}}$ = 5.21 $\pm$ 0.52 (stat) $\pm$ 0.56 (syst) fb = 5.21 $\pm$ 0.76 fb. Single-differential cross sections in photon, leading lepton, and leading jet transverse momenta, and double-differential cross sections in $m_{\mathrm{jj}}$ and $\lvert\Delta\eta_{\mathrm{jj}}\rvert$ are also measured. Exclusion limits on anomalous quartic gauge couplings are derived at 95% confidence level in terms of the effective field theory operators $\mathrm{M}_{0}$ to $\mathrm{M}_{5}$, $\mathrm{M}_{7}$, $\mathrm{T}_{0}$ to $\mathrm{T}_{2}$, and $\mathrm{T}_{5}$ to $\mathrm{T}_{9}$.
The expected and observed limits on the aQGC parameters at 95% confidence level. The last column presents the scattering energy values for which the amplitude would violate unitarity for the observed value of the aQGC parameter. All coupling parameter limits are set in TeV$^{-4}$, whereas the unitarity bounds are in TeV.
According to first-principle lattice QCD calculations, the transition from quark-gluon plasma to hadronic matter is a smooth crossover in the region μB ≤ T c. In this range the ratio, C6=C2, of net-baryon distributions are predicted to be negative. In this Letter, we report the first measurement of the midrapidity net-proton C6=C2 from 27, 54.4, and 200 GeV Au þ Au collisions at the Relativistic Heavy Ion Collider (RHIC). The dependence on collision centrality and kinematic acceptance in (p T , y) are analyzed. While for 27 and 54.4 GeV collisions the C6=C2 values are close to zero within uncertainties, it is observed that for 200 GeV collisions, the C6=C2 ratio becomes progressively negative from peripheral to central collisions. Transport model calculations without critical dynamics predict mostly positive values except for the most central collisions within uncertainties. These observations seem to favor a smooth crossover in the high-energy nuclear collisions at top RHIC energy.
Event by event net-proton multiplicity, $\Delta N_{p}$, distributions for Au+Au collisions at √sNN = 27, 54.4, and 200 GeV in 0-10% and 30-40% centralities at midrapidity (|y| < 0.5) for the transverse momentum range of 0.4 < $p_{T}$ (GeV/c) < 2.0. These distributions are normalized by the corresponding numbers of events and are not corrected for detector efficiencies. Statistical uncertainties are shown as vertical lines. The dashed lines show the Skellam distributions for each collision energy and centrality. The bottom panel shows the ratio of the data to the Skellam expectations.
Event by event net-proton multiplicity, $\Delta N_{p}$, distributions for Au+Au collisions at √sNN = 27, 54.4, and 200 GeV in 0-10% and 30-40% centralities at midrapidity (|y| < 0.5) for the transverse momentum range of 0.4 < $p_{T}$ (GeV/c) < 2.0. These distributions are normalized by the corresponding numbers of events and are not corrected for detector efficiencies. Statistical uncertainties are shown as vertical lines. The dashed lines show the Skellam distributions for each collision energy and centrality. The bottom panel shows the ratio of the data to the Skellam expectations.
Event by event net-proton multiplicity, $\Delta N_{p}$, distributions for Au+Au collisions at √sNN = 27, 54.4, and 200 GeV in 0-10% and 30-40% centralities at midrapidity (|y| < 0.5) for the transverse momentum range of 0.4 < $p_{T}$ (GeV/c) < 2.0. These distributions are normalized by the corresponding numbers of events and are not corrected for detector efficiencies. Statistical uncertainties are shown as vertical lines. The dashed lines show the Skellam distributions for each collision energy and centrality. The bottom panel shows the ratio of the data to the Skellam expectations.
Forum