Onia Processes

Production of any 3S1, 3PJ, and 3DJ charmonium and bottomonium states via the colour-singlet and colour-octet mechanisms. This includes by default, but is not limited to, production of the 3S1 J/psi and Upsilon and their radially excited states, as well as the 3PJ chi states and the 3D1 psi(3770). In each process the heavy quark content, either ccbar or bbbar, is followed by a round-bracketed expression which specifies the physical state in spectroscopic notation, (2S+1) L J. Proceding this is a square-bracketed expression, also in spectroscopic notation, which specifies the Fock state through which the process occurs, where (1) indicates a colour-singlet state and (8) a colour-octet state.

The unphysical colour-octet states follow the id scheme of 99 n_q n_s n_r n_L n_J where n_q is the quark flavour of the state and n_s is the colour-octet state type. Here 0 is 3S1, 1 is 1S0, and 2 is 3PJ. All remaining numbers follow the standard PDG numbering scheme. If a physical state is requested without a corresponding colour-octet state, a colour-octet state is automatically added to the ParticleData when a colour-octet process is selected. The colour-octet state is created with a mass given by the mass of the physical state plus the singlet-octet mass splitting parameter Onia:massSplit, which is by default set at 200 MeV, and decays exclusively to a gluon and the physical state. If the user wishes to manually set the mass splitting for each colour-octet state individually then Onia:forceMassSplit can be set to off. By default the widths of the octet states are set to vanish. This is not realistic, given their presumably rather rapid decay, but a nonvanishing width is not likely to have any measurable consequences that go beyond what comes from viewing the singlet-octet mass splitting as an effective parameter.

The original Fortran code for these processes has been contributed by Stefan Wolf [unpublished]. For the C++ version only the unpolarized expressions are retained, since the theoretical predictions of the colour-octet model anyway do not agree with the experimental observations. Furthermore, the polarization effects are modest, so isotropic decay is not a bad starting point. Such an event sample can afterwards be reweighted at will by the user, to test various assumptions. The expressions for the colour-singlet production of the 3S1 and 3PJ states can be found in [Bai83] and [Gas87]. Colour-octet expressions can be found in [Cho96] for the 1S0, 3S1, and 3PJ states, and the matrix elements for the 3DJ states are taken from [Yua98].

The implementation of charmonium and bottomonium production, including the colour-octet production mechanism, requires information on long-distance NRQCD matrix elements for the various wavefunctions involved. Default values for these are encoded in the O parameters and are taken from [Nas00]; see also [Bar07]. The 3DJ long-distance matrix elements are extracted from [Yua98].

Note that states that differ only by the radial excitation number n_r share the same short-distance matrix elements. The program has therefore been written such that further radial excitations can be easily added by editing this file, without requiring a recompilation of the code. All related arrays must be expanded in exactly the same way, however, i.e. the code of the colour singlet state, the long-distance matrix elements and the individual process on/off switches.

Double production of charmonium and bottomonium 3S1 states is also available, but with only the colour-singlet processes included. The short-distance matrix elements for these processes can be found in [Hum83] and [Qia02]. Arrays seperate to those used for single charmonium and bottomonium production set the long-distance NRQCD matrix elements. For consistency these values should be the same as those provided for single charmonium and bottomonium production, but are seperated to allow further control by the user if needed.

The description of final-state radiation is in this case based on some further model assumptions.

Most of the processes below are divergent in the limit pT → 0, and therefore a pTmin scale should be set. Comparisons with data indicate that this divergence can be tamed the same way as for the normal QCD 2 → 2 cross sections [Bar07,Kra08], which makes sense, since they are all dominated by the same kind of t-channel gluon exchange. It is therefore possible to use the SuppressSmallPT user hook to impose a reweighting that cancels the low-pT divergence.

An eikonalized description of these processes, excluding double onia production, is included in the multiparton-interactions framework. Here the low-pT damping is automatic, and additionally the framework is more consistent (e.g. with respect to energy-momentum constraints and the impact-parameter description) for events where the onium production is not the hardest subprocess, as would often be the case in the low-pT limit.

flag  Onia:forceMassSplit   (default = on)
Force the mass splitting between the colour-singlet states and their corresponding colour-octet state to be Onia:massSplit.

parm  Onia:massSplit   (default = 0.2; minimum = 0.0; maximum = 1.0)
Mass splitting in GeV between the physical colour-singlet states and their corresponding colour-octet state.

flag  Onia:all   (default = off)
Common switch for the group of onia production.

flag  Onia:all(3S1)   (default = off)
Common switch for the group of 3S1 onia production, e.g. J/psi and Upsilon.

flag  Onia:all(3PJ)   (default = off)
Common switch for the group of 3PJ onia production, e.g. chi_c and chi_b.

flag  Onia:all(3DJ)   (default = off)
Common switch for the group of 3DJ onia production, e.g. psi(3770).

flag  Charmonium:all   (default = off)
Common switch for the group of charmonium production, e.g. J/psi and chi_c.

flag  Bottomonium:all   (default = off)
Common switch for the group of bottomonium production, e.g. Upsilon and chi_b.

Charmonium 3S1 States

Warning: changed fvec, mvec or pvec values must be provided as a comma-separated list with the right number of elements, without any blanks inside the list.

mvec  Charmonium:states(3S1)   (default = 443,100443; minimum = 0)
The 3S1 charmonium states that can be produced from the following processes. Note that all vectors within this section, either of flags or parameters, must be the same length as this vector.

pvec  Charmonium:O(3S1)[3S1(1)]   (default = 1.16,0.76; minimum = 0.0)
The colour-singlet long-distance matrix elements <O[3S1(1)]> for the 3S1 charmonium states.

pvec  Charmonium:O(3S1)[3S1(8)]   (default = 0.0119,0.0050; minimum = 0.0)
The colour-octet long-distance matrix elements <O[3S1(8)]> for the 3S1 charmonium states.

pvec  Charmonium:O(3S1)[1S0(8)]   (default = 0.01,0.004; minimum = 0.0)
The colour-octet long-distance matrix elements <O[1S0(8)]> for the 3S1 charmonium states.

pvec  Charmonium:O(3S1)[3P0(8)]   (default = 0.01,0.004; minimum = 0.0)
The colour-octet long-distance matrix elements <O[3P0(8)]>/m_Q^2 for the 3S1 charmonium states. The remaining <O[3PJ(8)]>/m_Q^2 are calculated from these long-distance matrix elements.

fvec  Charmonium:gg2ccbar(3S1)[3S1(1)]g   (default = off,off)
Colour-singlet production of 3S1 charmonium states via g g → ccbar[3S1(1)] g. Code 401.

fvec  Charmonium:gg2ccbar(3S1)[3S1(8)]g   (default = off,off)
Colour-octet production of 3S1 charmonium states via g g → ccbar[3S1(8)] g. Code 402.

fvec  Charmonium:qg2ccbar(3S1)[3S1(8)]q   (default = off,off)
Colour-octet production of 3S1 charmonium states via q g → ccbar[3S1(8)] q. Code 403.

fvec  Charmonium:qqbar2ccbar(3S1)[3S1(8)]g   (default = off,off)
Colour-octet production of 3S1 charmonium states via q qbar → ccbar[3S1(8)] g. Code 404.

fvec  Charmonium:gg2ccbar(3S1)[1S0(8)]g   (default = off,off)
Colour-octet production of 3S1 charmonium states via g g → ccbar[1S0(8)] g. Code 405.

fvec  Charmonium:qg2ccbar(3S1)[1S0(8)]q   (default = off,off)
Colour-octet production of 3S1 charmonium states via q g → ccbar[1S0(8)] q. Code 406.

fvec  Charmonium:qqbar2ccbar(3S1)[1S0(8)]g   (default = off,off)
Colour-octet production of 3S1 charmonium states via q qbar → ccbar[1S0(8)] g. Code 407.

fvec  Charmonium:gg2ccbar(3S1)[3PJ(8)]g   (default = off,off)
Colour-octet production of 3S1 charmonium states via g g → ccbar[3PJ(8)] g. Code 408.

fvec  Charmonium:qg2ccbar(3S1)[3PJ(8)]q   (default = off,off)
Colour-octet production of 3S1 charmonium states via q g → ccbar[3PJ(8)] q. Code 409.

fvec  Charmonium:qqbar2ccbar(3S1)[3PJ(8)]g   (default = off,off)
Colour-octet production of 3S1 charmonium states via q qbar → ccbar[3SJ(8)] g. Code 410.

Charmonium 3PJ States

Warning: changed fvec, mvec or pvec values must be provided as a comma-separated list with the right number of elements, without any blanks inside the list.

mvec  Charmonium:states(3PJ)   (default = 10441,20443,445)
The 3PJ charmonium states that can be produced from the following processes. Note that all vectors within this section, either of flags or parameters, must be the same length as this vector.

pvec  Charmonium:O(3PJ)[3P0(1)]   (default = 0.05,0.05,0.05; minimum = 0.0)
The color-singlet long-distance matrix elements <O[3P0(1)]>/m_Q^2 for the 3PJ charmonium states. The remaining <O[3PJ(1)]>/m_Q^2 are calculated from these long-distance matrix elements.

pvec  Charmonium:O(3PJ)[3S1(8)]   (default = 0.0031,0.0031,0.0031; minimum = 0.0)
The color-singlet long-distance matrix elements <O[3S1(8)]> for the 3PJ charmonium states.

fvec  Charmonium:gg2ccbar(3PJ)[3PJ(1)]g   (default = off,off,off)
Colour-singlet production of 3PJ charmonium states via g g → ccbar[3PJ(1)] g. Code 411.

fvec  Charmonium:qg2ccbar(3PJ)[3PJ(1)]q   (default = off,off,off)
Colour-singlet production of 3PJ charmonium states via q g → ccbar[3PJ(1)] q. Code 412.

fvec  Charmonium:qqbar2ccbar(3PJ)[3PJ(1)]g   (default = off,off,off)
Colour-singlet production of 3PJ charmonium states via q qbar → ccbar[3PJ(1)] g. Code 413.

fvec  Charmonium:gg2ccbar(3PJ)[3S1(8)]g   (default = off,off,off)
Colour-octet production of 3PJ charmonium states via g g → ccbar[3S1(8)] g. Code 414.

fvec  Charmonium:qg2ccbar(3PJ)[3S1(8)]q   (default = off,off,off)
Colour-octet production of 3PJ charmonium states via q g → ccbar[3S1(8)] q. Code 415.

fvec  Charmonium:qqbar2ccbar(3PJ)[3S1(8)]g   (default = off,off,off)
Colour-octet production of 3PJ charmonium states via q qbar → ccbar[3S1(8)] g. Code 416.

Charmonium 3DJ States

Warning: changed fvec, mvec or pvec values must be provided as a comma-separated list with the right number of elements, without any blanks inside the list.

mvec  Charmonium:states(3DJ)   (default = 30443)
The 3DJ charmonium states that can be produced from the following processes. Note that all vectors within this section, either of flags or parameters, must be the same length as this vector.

pvec  Charmonium:O(3DJ)[3D1(1)]   (default = 0.161; minimum = 0.0)
The color-singlet long-distance matrix elements <O[3D1(1)]> for the 3PJ charmonium states. For a 3DJ charmonium state where J is not 1 the long distance matrix element <O[3DJ(1)]> is calculated by (2J+1)<O[3D1(1)]/3> using leading order spin symmetry relations.

pvec  Charmonium:O(3DJ)[3P0(8)]   (default = 0.01; minimum = 0.0)
The colour-octet long-distance matrix elements <O[3P0(8)]>/m_Q^2 for the 3DJ charmonium states. The remaining <O[3PJ(8)]>/m_Q^2 are calculated from these long-distance matrix elements.

fvec  Charmonium:gg2ccbar(3DJ)[3DJ(1)]g   (default = off)
Colour-singlet production of 3PJ charmonium states via g g → ccbar[3DJ(1)] g. Code 417.

fvec  Charmonium:gg2ccbar(3DJ)[3PJ(8)]g   (default = off)
Colour-octet production of 3DJ charmonium states via g g → ccbar[3PJ(8)] g. Code 418.

fvec  Charmonium:qg2ccbar(3DJ)[3PJ(8)]q   (default = off)
Colour-octet production of 3DJ charmonium states via q g → ccbar[3PJ(8)] q. Code 419.

fvec  Charmonium:qqbar2ccbar(3DJ)[3PJ(8)]g   (default = off)
Colour-octet production of 3DJ charmonium states via q qbar → ccbar[3PJ(8)] g. Code 420.

Double Charmonium 3S1 States

Warning: changed fvec, mvec or pvec values must be provided as a comma-separated list with the right number of elements, without any blanks inside the list.

mvec  Charmonium:states(3S1)1   (default = 443,443,100443)
The first 3S1 charmonium states of the charmonium pairs that can be produced from the following processes. Note that all vectors within this section, either of modes, flags, or parameters, must be the same length as this vector.

mvec  Charmonium:states(3S1)2   (default = 443,100443,100443)
The second 3S1 charmonium states of the charmonium pairs that can be produced from the following processes.

pvec  Charmonium:O(3S1)[3S1(1)]1   (default = 1.16,1.16,0.76; minimum = 0.0)
The color-singlet long-distance matrix elements <O[3S1(1)]> for the first 3S1 charmonium states.

pvec  Charmonium:O(3S1)[3S1(1)]2   (default = 1.16,0.76,0.76; minimum = 0.0)
The color-singlet long-distance matrix elements <O[3S1(1)]> for the second 3S1 charmonium states.

fvec  Charmonium:gg2doubleccbar(3S1)[3S1(1)]   (default = off,off,off)
Colour-singlet production of 3S1 charmonium state pairs via g g → ccbar[3S1(1)] ccbar[3S1(1)] . Code 421.

fvec  Charmonium:qqbar2doubleccbar(3S1)[3S1(1)]   (default = off,off,off)
Colour-singlet production of 3S1 charmonium state pairs via q qbar → ccbar[3S1(1)] ccbar[3S1(1)] . Code 422.

Bottomonium 3S1 States

Warning: changed fvec, mvec or pvec values must be provided as a comma-separated list with the right number of elements, without any blanks inside the list.

mvec  Bottomonium:states(3S1)   (default = 553,100553,200553; minimum = 0)
The 3S1 bottomonium states that can be produced from the following processes. Note that all vectors within this section, either of flags or parameters, must be the same length as this vector.

pvec  Bottomonium:O(3S1)[3S1(1)]   (default = 9.28,4.63,3.54; minimum = 0.0)
The colour-singlet long-distance matrix elements <O[3S1(1)]> for the 3S1 bottomonium states.

pvec  Bottomonium:O(3S1)[3S1(8)]   (default = 0.15,0.045,0.075; minimum = 0.0)
The colour-octet long-distance matrix elements <O[3S1(8)]> for the 3S1 bottomonium states.

pvec  Bottomonium:O(3S1)[1S0(8)]   (default = 0.02,0.006,0.01; minimum = 0.0)
The colour-octet long-distance matrix elements <O[1S0(8)]> for the 3S1 bottomonium states.

pvec  Bottomonium:O(3S1)[3P0(8)]   (default = 0.02,0.006,0.01; minimum = 0.0)
The colour-octet long-distance matrix elements <O[3P0(8)]>/m_Q^2 for the 3S1 bottomonium states. The remaining <O[3PJ(8)]>/m_Q^2 are calculated from these long-distance matrix elements.

fvec  Bottomonium:gg2bbbar(3S1)[3S1(1)]g   (default = off,off,off)
Colour-singlet production of 3S1 bottomonium states via g g → bbbar[3S1(1)] g. Code 501.

fvec  Bottomonium:gg2bbbar(3S1)[3S1(8)]g   (default = off,off,off)
Colour-octet production of 3S1 bottomonium states via g g → bbbar[3S1(8)] g. Code 502.

fvec  Bottomonium:qg2bbbar(3S1)[3S1(8)]q   (default = off,off,off)
Colour-octet production of 3S1 bottomonium states via q g → bbbar[3S1(8)] q. Code 503.

fvec  Bottomonium:qqbar2bbbar(3S1)[3S1(8)]g   (default = off,off,off)
Colour-octet production of 3S1 bottomonium states via q qbar → bbbar[3S1(8)] g. Code 504.

fvec  Bottomonium:gg2bbbar(3S1)[1S0(8)]g   (default = off,off,off)
Colour-octet production of 3S1 bottomonium states via g g → bbbar[1S0(8)] g. Code 505.

fvec  Bottomonium:qg2bbbar(3S1)[1S0(8)]q   (default = off,off,off)
Colour-octet production of 3S1 bottomonium states via q g → bbbar[1S0(8)] q. Code 506.

fvec  Bottomonium:qqbar2bbbar(3S1)[1S0(8)]g   (default = off,off,off)
Colour-octet production of 3S1 bottomonium states via q qbar → bbbar[1S0(8)] g. Code 507.

fvec  Bottomonium:gg2bbbar(3S1)[3PJ(8)]g   (default = off,off,off)
Colour-octet production of 3S1 bottomonium states via g g → bbbar[3PJ(8)] g. Code 508.

fvec  Bottomonium:qg2bbbar(3S1)[3PJ(8)]q   (default = off,off,off)
Colour-octet production of 3S1 bottomonium states via q g → bbbar[3PJ(8)] q. Code 509.

fvec  Bottomonium:qqbar2bbbar(3S1)[3PJ(8)]g   (default = off,off,off)
Colour-octet production of 3S1 bottomonium states via q qbar → bbbar[3SJ(8)] g. Code 510.

Bottomonium 3PJ States

Warning: changed fvec, mvec or pvec values must be provided as a comma-separated list with the right number of elements, without any blanks inside the list.

mvec  Bottomonium:states(3PJ)   (default = 10551,20553,555)
The 3PJ bottomonium states that can be produced from the following processes. Note that all vectors within this section, either of flags or parameters, must be the same length as this vector.

pvec  Bottomonium:O(3PJ)[3P0(1)]   (default = 0.085,0.085,0.085; minimum = 0.0)
The color-singlet long-distance matrix elements <O[3P0(1)]>/m_Q^2 for the 3PJ bottomonium states. The remaining <O[3PJ(1)]>/m_Q^2 are calculated from these long-distance matrix elements.

pvec  Bottomonium:O(3PJ)[3S1(8)]   (default = 0.04,0.04,0.04; minimum = 0.0)
The color-singlet long-distance matrix elements <O[3S1(8)]> for the 3PJ bottomonium states.

fvec  Bottomonium:gg2bbbar(3PJ)[3PJ(1)]g   (default = off,off,off)
Colour-singlet production of 3PJ bottomonium states via g g → bbbar[3PJ(1)] g. Code 511.

fvec  Bottomonium:qg2bbbar(3PJ)[3PJ(1)]q   (default = off,off,off)
Colour-singlet production of 3PJ bottomonium states via q g → bbbar[3PJ(1)] q. Code 512.

fvec  Bottomonium:qqbar2bbbar(3PJ)[3PJ(1)]g   (default = off,off,off)
Colour-singlet production of 3PJ bottomonium states via q qbar → bbbar[3PJ(1)] g. Code 513.

fvec  Bottomonium:gg2bbbar(3PJ)[3S1(8)]g   (default = off,off,off)
Colour-octet production of 3PJ bottomonium states via g g → bbbar[3S1(8)] g. Code 514.

fvec  Bottomonium:qg2bbbar(3PJ)[3S1(8)]q   (default = off,off,off)
Colour-octet production of 3PJ bottomonium states via q g → bbbar[3S1(8)] q. Code 515.

fvec  Bottomonium:qqbar2bbbar(3PJ)[3S1(8)]g   (default = off,off,off)
Colour-octet production of 3PJ bottomonium states via q qbar → bbbar[3S1(8)] g. Code 516.

Bottomonium 3DJ States

Warning: changed fvec, mvec or pvec values must be provided as a comma-separated list with the right number of elements, without any blanks inside the list.

mvec  Bottomonium:states(3DJ)  
The 3DJ bottomonium states that can be produced from the following processes. Currently, no 3DJ states are included in the default ParticleData and so none are included here. Note that all vectors within this section, either of flags or parameters, must be the same length as this vector.

pvec  Bottomonium:O(3DJ)[3D1(1)]   (; minimum = 0.0)
The color-singlet long-distance matrix elements <O[3D1(1)]> for the 3PJ bottomonium states. For a 3DJ bottomonium state where J is not 1 the long distance matrix element <O[3DJ(1)]> is calculated by (2J+1)<O[3D1(1)]/3> using leading order spin symmetry relations.

pvec  Bottomonium:O(3DJ)[3P0(8)]   (; minimum = 0.0)
The colour-octet long-distance matrix elements <O[3P0(8)]>/m_Q^2 for the 3DJ bottomonium states. The remaining <O[3PJ(8)]>/m_Q^2 are calculated from these long-distance matrix elements.

fvec  Bottomonium:gg2bbbar(3DJ)[3DJ(1)]g  
Colour-singlet production of 3PJ bottomonium states via g g → bbbar[3DJ(1)] g. Code 517.

fvec  Bottomonium:gg2bbbar(3DJ)[3PJ(8)]g  
Colour-octet production of 3DJ bottomonium states via g g → bbbar[3PJ(8)] g. Code 518.

fvec  Bottomonium:qg2bbbar(3DJ)[3PJ(8)]q  
Colour-octet production of 3DJ bottomonium states via q g → bbbar[3PJ(8)] q. Code 519.

fvec  Bottomonium:qqbar2bbbar(3DJ)[3PJ(8)]g  
Colour-octet production of 3DJ bottomonium states via q qbar → bbbar[3PJ(8)] g. Code 520.

Double Bottomonium 3S1 States

Warning: changed fvec, mvec or pvec values must be provided as a comma-separated list with the right number of elements, without any blanks inside the list.

mvec  Bottomonium:states(3S1)1   (default = 553,553,553,100553,100553,200553)
The first 3S1 bottomonium states of the bottomonium pairs that can be produced from the following processes. Note that all vectors within this section, either of modes, flags, or parameters, must be the same length as this vector.

mvec  Bottomonium:states(3S1)2   (default = 553,100553,200553,100553,200553,200553)
The second 3S1 bottomonium states of the bottomonium pairs that can be produced from the following processes.

pvec  Bottomonium:O(3S1)[3S1(1)]1   (default = 9.28,9.28,9.28,4.63,4.63,3.54; minimum = 0.0)
The color-singlet long-distance matrix elements <O[3S1(1)]> for the first 3S1 bottomonium states.

pvec  Bottomonium:O(3S1)[3S1(1)]2   (default = 9.28,4.63,3.54,4.63,3.54,3.54; minimum = 0.0)
The color-singlet long-distance matrix elements <O[3S1(1)]> for the second 3S1 bottomonium states.

fvec  Bottomonium:gg2doublebbbar(3S1)[3S1(1)]   (default = off,off,off,off,off,off)
Colour-singlet production of 3S1 bottomonium state pairs via g g → bbbar[3S1(1)] bbbar[3S1(1)] . Code 521.

fvec  Bottomonium:qqbar2doublebbbar(3S1)[3S1(1)]   (default = off,off,off,off,off,off)
Colour-singlet production of 3S1 bottomonium state pairs via q qbar → bbbar[3S1(1)] bbbar[3S1(1)] . Code 522.