Date

Combination and interpretation of differential Higgs boson production cross sections in proton-proton collisions at $\sqrt{s}$ = 13 TeV

The CMS collaboration Chekhovsky, Vladimir ; Hayrapetyan, Aram ; Makarenko, Vladimir ; et al.
CMS-HIG-23-013, 2025.
Inspire Record 2913615 DOI 10.17182/hepdata.156816

Precision measurements of Higgs boson differential production cross sections are a key tool to probe the properties of the Higgs boson and test the standard model. New physics can affect both Higgs boson production and decay, leading to deviations from the distributions that are expected in the standard model. In this paper, combined measurements of differential spectra in a fiducial region matching the experimental selections are performed, based on analyses of four Higgs boson decay channels ($\gamma\gamma$, ZZ$^{(*)}$, WW$^{(*)}$, and $\tau\tau$) using proton-proton collision data recorded with the CMS detector at $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The differential measurements are extrapolated to the full phase space and combined to provide the differential spectra. A measurement of the total Higgs boson production cross section is also performed using the $\gamma\gamma$ and ZZ decay channels, with a result of 53.4$^{+2.9}_{-2.9}$ (stat)$^{+1.9}_{-1.8}$ (syst) pb, consistent with the standard model prediction of 55.6 $\pm$ 2.5 pb. The fiducial measurements are used to compute limits on Higgs boson couplings using the $\kappa$-framework and the SM effective field theory.

17 data tables

Observed best fit differential cross section for the $p_{T}^{H}$ observable

Observed best fit differential cross section for the $N_{jets}$ observable

Observed best fit differential cross section for the $p_{T}^{j1}$ (GeV) observable

More…

Energy-scaling behavior of intrinsic transverse momentum parameters in Drell-Yan simulation

The CMS collaboration Hayrapetyan, Aram ; Tumasyan, Armen ; Adam, Wolfgang ; et al.
Phys.Rev.D 111 (2025) 072003, 2025.
Inspire Record 2839223 DOI 10.17182/hepdata.154142

An analysis is presented based on models of the intrinsic transverse momentum (intrinsic $k_\mathrm{T}$) of partons in nucleons by studying the dilepton transverse momentum in Drell-Yan events. Using parameter tuning in event generators and existing data from fixed-target experiments and from hadron colliders, our investigation spans three orders of magnitude in center-of-mass energy and two orders of magnitude in dilepton invariant mass. The results show an energy-scaling behavior of the intrinsic $k_\mathrm{T}$ parameters, independent of the dilepton invariant mass at a given center-of-mass energy.

45 data tables

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP4 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP3 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

More…

ATLAS Run 2 searches for electroweak production of supersymmetric particles interpreted within the pMSSM

The ATLAS collaboration Aad, Georges ; Abbott, Braden Keim ; Abeling, Kira ; et al.
JHEP 05 (2024) 106, 2024.
Inspire Record 2755168 DOI 10.17182/hepdata.149493

A summary of the constraints from searches performed by the ATLAS Collaboration for the electroweak production of charginos and neutralinos is presented. Results from eight separate ATLAS searches are considered, each using 140 fb$^{-1}$ of proton-proton data at a centre-of-mass energy of $\sqrt{s}$=13 TeV collected at the Large Hadron Collider during its second data-taking run. The results are interpreted in the context of the 19-parameter phenomenological minimal supersymmetric standard model, where R-parity conservation is assumed and the lightest supersymmetric particle is assumed to be the lightest neutralino. Constraints from previous electroweak, flavour and dark matter related measurements are also considered. The results are presented in terms of constraints on supersymmetric particle masses and are compared with limits from simplified models. Also shown is the impact of ATLAS searches on parameters such as the dark matter relic density and the spin-dependent and spin-independent scattering cross-sections targeted by direct dark matter detection experiments. The Higgs boson and Z boson `funnel regions', where a low-mass neutralino would not oversaturate the dark matter relic abundance, are almost completely excluded by the considered constraints. Example spectra for non-excluded supersymmetric models with light charginos and neutralinos are also presented.

2 data tables

SLHA files and exclusion information (in CSV format) are available to download for the pMSSM models in this paper. Please refer to <a href="https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2020-15/inputs/ATLAS_EW_pMSSM_Run2.html">this web page</a> for download links along with a description of the contents.

SLHA files and exclusion information (in CSV format) are available to download for the pMSSM models in this paper. Please refer to <a href="https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2020-15/inputs/ATLAS_EW_pMSSM_Run2.html">this web page</a> for download links along with a description of the contents.


Constraining the Higgs boson self-coupling from single- and double-Higgs production with the ATLAS detector using $pp$ collisions at $\sqrt{s}=13$ TeV

The ATLAS collaboration Aad, Georges ; Abbott, Braden Keim ; Abbott, Dale ; et al.
Phys.Lett.B 843 (2023) 137745, 2023.
Inspire Record 2175556 DOI 10.17182/hepdata.135471

Constraints on the Higgs boson self-coupling are set by combining double-Higgs boson analyses in the $b\bar{b}b\bar{b}$, $b\bar{b}\tau^+\tau^-$ and $b\bar{b} \gamma \gamma$ decay channels with single-Higgs boson analyses targeting the $\gamma \gamma$, $ZZ^*$, $WW^*$, $\tau^+ \tau^-$ and $b\bar{b}$ decay channels. The data used in these analyses were recorded by the ATLAS detector at the LHC in proton$-$proton collisions at $\sqrt{s}=13$ TeV and correspond to an integrated luminosity of 126$-$139 fb$^{-1}$. The combination of the double-Higgs analyses sets an upper limit of $\mu_{HH} < 2.4$ at 95% confidence level on the double-Higgs production cross-section normalised to its Standard Model prediction. Combining the single-Higgs and double-Higgs analyses, with the assumption that new physics affects only the Higgs boson self-coupling ($\lambda_{HHH}$), values outside the interval $-0.4< \kappa_{\lambda}=(\lambda_{HHH}/\lambda_{HHH}^{\textrm{SM}})< 6.3$ are excluded at 95% confidence level. The combined single-Higgs and double-Higgs analyses provide results with fewer assumptions, by adding in the fit more coupling modifiers introduced to account for the Higgs boson interactions with the other Standard Model particles. In this relaxed scenario, the constraint becomes $-1.4 < \kappa_{\lambda} < 6.1$ at 95% CL.

44 data tables

Observed and expected 95% CL upper limits on the signal strength for double-Higgs production from the bbbb, bb$\tau\tau$ and bb$\gamma\gamma$ decay channels, and their statistical combination. The value $m_H$ = 125.09 GeV is assumed when deriving the predicted SM cross-section. The expected limit and the corresponding error bands are derived assuming the absence of the HH process and with all nuisance parameters profiled to the observed data.

Observed and expected 95% CL exclusion limits on the production cross-sections of the combined ggF HH and VBF HH processes as a function of $\kappa_\lambda$, for the three double-Higgs search channels and their combination. The expected limits assume no HH production. The red line shows the theory prediction for the combined ggF HH and VBF HH cross-section as a function of $\kappa_\lambda$ where all parameters and couplings are set to their SM values except for $\kappa_\lambda$. The band surrounding the red cross-section lines indicate the theoretical uncertainty of the predicted cross-section.

Observed and expected 95% CL exclusion limits on the production cross-sections of the VBF HH process as a function of $\kappa_{2V}$, for the three double-Higgs search channels and their combination. The expected limits assume no VBF HH production. The red line shows the predicted VBF HH cross-section as a function of $\kappa_{2V}$. The bands surrounding the red cross-section lines indicate the theoretical uncertainty of the predicted cross-section. The uncertainty band is smaller than the width of the plotted line.

More…

Direct observation of the dead-cone effect in QCD

The ALICE collaboration Acharya, S. ; Acharya, S. ; Adamova, D. ; et al.
Nature 605 (2022) 440-446, 2022.
Inspire Record 1867966 DOI 10.17182/hepdata.130725

In particle collider experiments, elementary particle interactions with large momentum transfer produce quarks and gluons (known as partons) whose evolution is governed by the strong force, as described by the theory of quantum chromodynamics (QCD). These partons subsequently emit further partons in a process that can be described as a parton shower which culminates in the formation of detectable hadrons. Studying the pattern of the parton shower is one of the key experimental tools for testing QCD. This pattern is expected to depend on the mass of the initiating parton, through a phenomenon known as the dead-cone effect, which predicts a suppression of the gluon spectrum emitted by a heavy quark of mass $m_{\rm{Q}}$ and energy $E$, within a cone of angular size $m_{\rm{Q}}$/$E$ around the emitter. Previously, a direct observation of the dead-cone effect in QCD had not been possible, owing to the challenge of reconstructing the cascading quarks and gluons from the experimentally accessible hadrons. We report the direct observation of the QCD dead cone by using new iterative declustering techniques to reconstruct the parton shower of charm quarks. This result confirms a fundamental feature of QCD. Furthermore, the measurement of a dead-cone angle constitutes a direct experimental observation of the non-zero mass of the charm quark, which is a fundamental constant in the standard model of particle physics.

1 data table

The $R(\theta)$ variable for charm/inclusive emissions in three bins of $E_{Rad}$: 5-10, 10-20 and 20-35 GeV.


Updated MiniBooNE Neutrino Oscillation Results with Increased Data and New Background Studies

The MiniBooNE collaboration Aguilar-Arevalo, A.A. ; Brown, B.C. ; Conrad, J.M. ; et al.
Phys.Rev.D 103 (2021) 052002, 2021.
Inspire Record 1804293 DOI 10.17182/hepdata.114365

The MiniBooNE experiment at Fermilab reports a total excess of $638.0 \pm 132.8$ electron-like events ($4.8 \sigma$) from a data sample corresponding to $18.75 \times 10^{20}$ protons-on-target in neutrino mode, which is a 46\% increase in the data sample with respect to previously published results, and $11.27 \times 10^{20}$ protons-on-target in antineutrino mode. The additional statistics allow several studies to address questions on the source of the excess. First, we provide two-dimensional plots in visible energy and cosine of the angle of the outgoing lepton, which can provide valuable input to models for the event excess. Second, we test whether the excess may arise from photons that enter the detector from external events or photons exiting the detector from $\pi^0$ decays in two model independent ways. Beam timing information shows that almost all of the excess is in time with neutrinos that interact in the detector. The radius distribution shows that the excess is distributed throughout the volume, while tighter cuts on the fiducal volume increase the significance of the excess. We conclude that models of the event excess based on entering and exiting photons are disfavored.

15 data tables

The frequentist $1\sigma$ confidence region in $\sin^2(2\theta)$ $\Delta m^2$ for a 2-neutrino muon-to-electron oscillation fit.

The frequentist $90\%$ confidence region in $\sin^2(2\theta)$ $\Delta m^2$ for a 2-neutrino muon-to-electron oscillation fit.

The frequentist $99\%$ confidence region in $\sin^2(2\theta)$ $\Delta m^2$ for a 2-neutrino muon-to-electron oscillation fit.

More…

Beam energy dependence of net-$\Lambda$ fluctuations measured by the STAR experiment at RHIC

The STAR collaboration Adam, Jaroslav ; Adamczyk, Leszek ; Adams, Joseph ; et al.
Phys.Rev.C 102 (2020) 024903, 2020.
Inspire Record 1776194 DOI 10.17182/hepdata.113523

The measurements of particle multiplicity distributions have generated considerable interest in understanding the fluctuations of conserved quantum numbers in the Quantum Chromodynamics (QCD) hadronization regime, in particular near a possible critical point and near the chemical freeze-out. We report the measurement of efficiency and centrality bin width corrected cumulant ratios ($C_{2}/C_{1}$, $C_{3}/C_{2}$) of net-$\Lambda$ distributions, in the context of both strangeness and baryon number conservation, as a function of collision energy, centrality and rapidity. The results are for Au + Au collisions at five beam energies ($\sqrt{s_{NN}}$ = 19.6, 27, 39, 62.4 and 200 GeV) recorded with the Solenoidal Tracker at RHIC (STAR). We compare our results to the Poisson and negative binomial (NBD) expectations, as well as to Ultra-relativistic Quantum Molecular Dynamics (UrQMD) and Hadron Resonance Gas (HRG) model predictions. Both NBD and Poisson baselines agree with data within the statistical and systematic uncertainties. The ratios of the measured cumulants show no features of critical fluctuations. The chemical freeze-out temperatures extracted from a recent HRG calculation, which was successfully used to describe the net-proton, net-kaon and net-charge data, indicate $\Lambda$ freeze-out conditions similar to those of kaons. However, large deviations are found when comparing to temperatures obtained from net-proton fluctuations. The net-$\Lambda$ cumulants show a weak, but finite, dependence on the rapidity coverage in the acceptance of the detector, which can be attributed to quantum number conservation.

35 data tables

Centrality dependence of single cumulants C1, of net-lambda multiplicity distributions at Au + Au collision 19.6 GeV. Values are shown with NBD, Poisson and UrQMD predictions. Npart values are from Phys. Rev. C 104, 024902 (2021) and they are little different than the values shown in the original paper.

Centrality dependence of single cumulants C1, of net-lambda multiplicity distributions at Au + Au collision 27 GeV. Values are shown with NBD, Poisson and UrQMD predictions.Npart values are from Phys. Rev. C 104, 024902 (2021) and they are little different than the values shown in the original paper.

Centrality dependence of single cumulants C1, of net-lambda multiplicity distributions at Au + Au collision 39 GeV. Values are shown with NBD, Poisson and UrQMD predictions.Npart values are from Phys. Rev. C 104, 024902 (2021) and they are little different than the values shown in the original paper.

More…

Bulk Properties of the System Formed in Au+Au Collisions at $\sqrt{s_{\mathrm{NN}}}$ = 14.5 GeV

The STAR collaboration Adam, Jaroslav ; Adamczyk, Leszek ; Adams, Joseph ; et al.
Phys.Rev.C 101 (2020) 024905, 2020.
Inspire Record 1748776 DOI 10.17182/hepdata.103857

We report systematic measurements of bulk properties of the system created in Au+Au collisions at $\sqrt{s_{\mathrm{NN}}}$ = 14.5 GeV recorded by the STAR detector at the Relativistic Heavy Ion Collider (RHIC).The transverse momentum spectra of $\pi^{\pm}$, $K^{\pm}$ and $p(\bar{p})$ are studied at mid-rapidity ($|y| < 0.1$) for nine centrality intervals. The centrality, transverse momentum ($p_T$),and pseudorapidity ($\eta$) dependence of inclusive charged particle elliptic flow ($v_2$), and rapidity-odd charged particles directed flow ($v_{1}$) results near mid-rapidity are also presented. These measurements are compared with the published results from Au+Au collisions at other energies, and from Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ = 2.76 TeV. The results at $\sqrt{s_{\mathrm{NN}}}$ = 14.5 GeV show similar behavior as established at other energies and fit well in the energy dependence trend. These results are important as the 14.5 GeV energy fills the gap in $\mu_B$, which is of the order of 100 MeV,between $\sqrt{s_{\mathrm{NN}}}$ =11.5 and 19.6 GeV. Comparisons of the data with UrQMD and AMPT models show poor agreement in general.

42 data tables

The $p_{T}$ spectra of proton measured at midrapidity (|y|<0.1) in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV. Spectra are plotted for nine centrality classes, with some spectra multiplied by a scale factor to improve clarity, as indicated in the legend

The $p_{T}$ spectra of antiproton measured at midrapidity (|y|<0.1) in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV. Spectra are plotted for nine centrality classes, with some spectra multiplied by a scale factor to improve clarity, as indicatedin the legend

The $p_{T}$ spectra of $\pi^{+}$ measured at midrapidity (|y|<0.1) in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV. Spectra are plotted for nine centrality classes, with some spectra multiplied by a scale factor to improve clarity, as indicatedin the legend

More…

Version 2
Collision energy dependence of second-order off-diagonal and diagonal cumulants of net-charge, net-proton and net-kaon multiplicity distributions in Au+Au collisions

The STAR collaboration Adam, Jaroslav ; Adamczyk, Leszek ; Adams, Joseph ; et al.
Phys.Rev.C 100 (2019) 014902, 2019.
Inspire Record 1724809 DOI 10.17182/hepdata.105908

We report the first measurements of a complete second-order cumulant matrix of net-charge, net-proton, and net-kaon multiplicity distributions for the first phase of the beam energy scan program at RHIC. This includes the centrality and, for the first time, the pseudorapidity window dependence of both diagonal and off-diagonal cumulants in Au+Au collisions at \sNN~= 7.7-200 GeV. Within the available acceptance of $|\eta|<0.5$, the cumulants grow linearly with the pseudorapidity window. Relative to the corresponding measurements in peripheral collisions, the ratio of off-diagonal over diagonal cumulants in central collisions indicates an excess correlation between net-charge and net-kaon, as well as between net-charge and net-proton. The strength of such excess correlation increases with the collision energy. The correlation between net-proton and net-kaon multiplicity distributions is observed to be negative at \sNN~= 200 GeV and change to positive at the lowest collision energy. Model calculations based on non-thermal (UrQMD) and thermal (HRG) production of hadrons cannot explain the data. These measurements will help map the QCD phase diagram, constrain hadron resonance gas model calculations, and provide new insights on the energy dependence of baryon-strangeness correlations. An erratum has been added to address the issue of self-correlation in the previously considered efficiency correction for off-diagonal cumulant measurement. Previously considered unidentified (net-)charge correlation results ($\sigma^{11}_{Q,p}$ and $\sigma^{11}_{Q,k})$ are now replaced with identified (net-)charge correlation ($\sigma^{11}_{Q^{PID},p}$ and $\sigma^{11}_{Q^{PID},k}$)

31 data tables

The dependence of efficiency corrected second-order diagonal and off-diagonal cumulants on the width of the η-window. The filled and open circles represent 0-5% and 70-80% central collisions respectively. The shaded band represents the systematic uncertainty. The statistical uncertainties are within the marker size and solid lines are UrQMD calculations.

The dependence of efficiency corrected second-order diagonal and off-diagonal cumulants on the width of the η-window. The filled and open circles represent 0-5% and 70-80% central collisions respectively. The shaded band represents the systematic uncertainty. The statistical uncertainties are within the marker size and solid lines are UrQMD calculations.

Centrality dependence of efficiency corrected second-order diagonal cumulants of net-proton, net-kaon and net-pion (top to bottom) of the multiplicity distributions for Au+Au collisions at GeV (left to right) within kinematic range of |η| < 0.5 and 0.4 < pT < 1.6 GeV/c. The boxes represent the systematic error. The statistical error bars are within the marker size. The dashed lines represent scaling predicted by central limit theorem and the solid lines are UrQMD calculations.

More…

Light isovector resonances in $\pi^- p \to \pi^-\pi^-\pi^+ p$ at 190 GeV/${\it c}$

The COMPASS collaboration Aghasyan, M. ; Alexeev, M.G. ; Alexeev, G.D. ; et al.
Phys.Rev.D 98 (2018) 092003, 2018.
Inspire Record 1655631 DOI 10.17182/hepdata.82958

We have performed the most comprehensive resonance-model fit of $\pi^-\pi^-\pi^+$ states using the results of our previously published partial-wave analysis (PWA) of a large data set of diffractive-dissociation events from the reaction $\pi^- + p \to \pi^-\pi^-\pi^+ + p_\text{recoil}$ with a 190 GeV/$c$ pion beam. The PWA results, which were obtained in 100 bins of three-pion mass, $0.5 &lt; m_{3\pi} &lt; 2.5$ GeV/$c^2$, and simultaneously in 11 bins of the reduced four-momentum transfer squared, $0.1 &lt; t' &lt; 1.0$ $($GeV$/c)^2$, are subjected to a resonance-model fit using Breit-Wigner amplitudes to simultaneously describe a subset of 14 selected waves using 11 isovector light-meson states with $J^{PC} = 0^{-+}$, $1^{++}$, $2^{++}$, $2^{-+}$, $4^{++}$, and spin-exotic $1^{-+}$ quantum numbers. The model contains the well-known resonances $\pi(1800)$, $a_1(1260)$, $a_2(1320)$, $\pi_2(1670)$, $\pi_2(1880)$, and $a_4(2040)$. In addition, it includes the disputed $\pi_1(1600)$, the excited states $a_1(1640)$, $a_2(1700)$, and $\pi_2(2005)$, as well as the resonancelike $a_1(1420)$. We measure the resonance parameters mass and width of these objects by combining the information from the PWA results obtained in the 11 $t'$ bins. We extract the relative branching fractions of the $\rho(770) \pi$ and $f_2(1270) \pi$ decays of $a_2(1320)$ and $a_4(2040)$, where the former one is measured for the first time. In a novel approach, we extract the $t'$ dependence of the intensity of the resonances and of their phases. The $t'$ dependence of the intensities of most resonances differs distinctly from the $t'$ dependence of the nonresonant components. For the first time, we determine the $t'$ dependence of the phases of the production amplitudes and confirm that the production mechanism of the Pomeron exchange is common to all resonances.

2 data tables

Real and imaginary parts of the normalized transition amplitudes $\mathcal{T}_a$ of the 14 selected partial waves in the 1100 $(m_{3\pi}, t')$ cells (see Eq. (12) in the paper). The wave index $a$ represents the quantum numbers that uniquely define the partial wave. The quantum numbers are given by the shorthand notation $J^{PC} M^\varepsilon [$isobar$] \pi L$. We use this notation to label the transition amplitudes in the column headers. The $m_{3\pi}$ values that are given in the first column correspond to the bin centers. Each of the 100 $m_{3\pi}$ bins is 20 MeV/$c^2$ wide. Since the 11 $t'$ bins are non-equidistant, the lower and upper bounds of each $t'$ bin are given in the column headers. The transition amplitudes define the spin-density matrix elements $\varrho_{ab}$ for waves $a$ and $b$ according to Eq. (18). The spin-density matrix enters the resonance-model fit via Eqs. (33) and (34). The transition amplitudes are normalized via Eqs. (9), (16), and (17) such that the partial-wave intensities $\varrho_{aa} = |\mathcal{T}_a|^2$ are given in units of acceptance-corrected number of events. The relative phase $\Delta\phi_{ab}$ between two waves $a$ and $b$ is given by $\arg(\varrho_{ab}) = \arg(\mathcal{T}_a) - \arg(\mathcal{T}_b)$. Note that only relative phases are well-defined. The phase of the $1^{++}0^+ \rho(770) \pi S$ wave was set to $0^\circ$ so that the corresponding transition amplitudes are real-valued. In the PWA model, some waves are excluded in the region of low $m_{3\pi}$ (see paper and [Phys. Rev. D 95, 032004 (2017)] for a detailed description of the PWA model). For these waves, the transition amplitudes are set to zero. The tables with the covariance matrices of the transition amplitudes for all 1100 $(m_{3\pi}, t')$ cells can be downloaded via the 'Additional Resources' for this table.

Decay phase-space volume $I_{aa}$ for the 14 selected partial waves as a function of $m_{3\pi}$, normalized such that $I_{aa}(m_{3\pi} = 2.5~\text{GeV}/c^2) = 1$. The wave index $a$ represents the quantum numbers that uniquely define the partial wave. The quantum numbers are given by the shorthand notation $J^{PC} M^\varepsilon [$isobar$] \pi L$. We use this notation to label the decay phase-space volume in the column headers. The labels are identical to the ones used in the column headers of the table of the transition amplitudes. $I_{aa}$ is calculated using Monte Carlo integration techniques for fixed $m_{3\pi}$ values, which are given in the first column, in the range from 0.5 to 2.5 GeV/$c^2$ in steps of 10 MeV/$c^2$. The statistical uncertainties given for $I_{aa}$ are due to the finite number of Monte Carlo events. $I_{aa}(m_{3\pi})$ is defined in Eq. (6) in the paper and appears in the resonance model in Eqs. (19) and (20).