Resonance Production

The simplest way to produce a resonance is by a process.
If the decay of the resonance is not considered, the cross-section
formula does not depend on , but takes the form

In published formulae, cross sections are often given in the zero-width approximation, i.e. , where is the mass of the resonance. Introducing the scaled mass , this corresponds to a delta function , which can be used to eliminate the integral over .

However, what we normally want to do is replace the
function by the appropriate Breit-Wigner shape. For a resonance
width this is achieved by the replacement

An improved description of resonance shapes is obtained if the width is made -dependent (occasionally also referred to as mass-dependent width, since is not always the resonance mass), see e.g. [Ber89]. To first approximation, this means that the expression is to be replaced by , both in the numerator and the denominator. An intermediate step is to perform this replacement only in the numerator. This is convenient when not only -channel resonance production is simulated but also non-resonance - or -channel graphs are involved, since mass-dependent widths in the denominator here may give an imperfect cancellation of divergences. (More about this below.)

To be more precise, in the program the quantity
is introduced, and the Breit-Wigner is written as

A decay to a fermion pair, , gives no contribution below threshold, i.e. for . Above threshold, is proportional to , multiplied by a threshold factor for the vector part of a spin 1 resonance, by for the axial vector part, by for a scalar resonance and by for a pseudoscalar one. Here . For the decay into unequal masses, e.g. of the , corresponding but more complicated expressions are used.

For decays into a quark pair, a first-order strong correction factor is included in . This is the correct choice for all spin 1 colourless resonances, but is here used for all resonances where no better knowledge is available. Currently the major exception is top decay, where the factor is used to approximate loop corrections [Jez89]. The second-order corrections are often known, but then are specific to each resonance, and are not included. An option exists for the resonances, where threshold effects due to bound-state formation are taken into account in a smeared-out, average sense, see eq. ().

For other decay channels, not into fermion pairs, the dependence is typically more complicated. An example would be the decay , with a nontrivial threshold and a subtle energy dependence above that [Sey95a]. Since a Higgs with could still decay in this channel, it is in fact necessary to perform a two-dimensional integral over the Breit-Wigner mass distributions to obtain the correct result (and this has to be done numerically, at least in part). Fortunately, a Higgs particle lighter than is sufficiently narrow that the integral only needs to be performed once and for all at initialization (whereas most other partial widths are recalculated whenever needed). Channels that proceed via loops, such as , also display complicated threshold behaviours.

The coupling structure within the electroweak sector is usually
(re)expressed in terms of gauge boson masses,
and
, i.e. factors of
are replaced according to

(87) |

For Higgs particles and technipions, fermion masses enter not only
in the kinematics but also as couplings. The latter kind of quark masses
(but not the former, at least not in the program) are running with the
scale of the process, i.e. normally the resonance mass. The expression
used is [Car96]

In summary, we see that an dependence may enter several different ways into the expressions from which the total is built up.

When only decays to a specific final state are considered, the in the denominator remains the sum over all allowed decay channels, but the numerator only contains the term of the final state considered.

If the combined production and decay process is
considered, the same dependence is implicit in the
coupling structure of as one would have had in
, i.e. to first approximation there is a symmetry between
couplings of a resonance to the initial and to the final state.
The cross section
is therefore, in the program, written in the form

(90) |

(91) |

The analogy between the and cannot be pushed too far, however. The two differ in several important aspects. Firstly, colour factors appear reversed: the decay contains a colour factor enhancement, while is instead suppressed by a factor . Secondly, the first-order correction factor for the final state has to be replaced by a more complicated factor for the initial state. This factor is not known usually, or it is known (to first non-trivial order) but too lengthy to be included in the program. Thirdly, incoming partons as a rule are space-like. All the threshold suppression factors of the final-state expressions are therefore irrelevant when production is considered. In sum, the analogy between and is mainly useful as a consistency cross-check, while the two usually are calculated separately. Exceptions include the rather messy loop structure involved in and , which is only coded once.

It is of some interest to consider the observable resonance shape when the effects of parton distributions are included. In a hadron collider, to first approximation, parton distributions tend to have a behaviour roughly like for small -- this is why is replaced by in eq. (). Instead, the basic parton-distribution behaviour is shifted into the factor of in the integration phase space , cf. eq. (). When convoluted with the Breit-Wigner shape, two effects appear. One is that the overall resonance is tilted: the low-mass tail is enhanced and the high-mass one suppressed. The other is that an extremely long tail develops on the low-mass side of the resonance: when , eq. () with gives a , which exactly cancels the factor mentioned above. Naïvely, the integral over , , therefore gives a net logarithmic divergence of the resonance shape when . Clearly, it is then necessary to consider the shape of the parton distributions in more detail. At not-too-small , the evolution equations in fact lead to parton distributions more strongly peaked than , typically with , and therefore a divergence like in the cross-section expression. Eventually this divergence is regularized by a closing of the phase space, i.e. that vanishes faster than , and by a less drastic small- parton-distribution behaviour when .

The secondary peak at small may give a rather high cross section, which can even rival that of the ordinary peak around the nominal mass. This is the case, for instance, with production. Such a peak has never been observed experimentally, but this is not surprising, since the background from other processes is overwhelming at low . Thus a lepton of one or a few GeV of transverse momentum is far more likely to come from the decay of a charm or bottom hadron than from an extremely off-shell of a mass of a few GeV. When resonance production is studied, it is therefore important to set limits on the mass of the resonance, so as to agree with the experimental definition, at least to first approximation. If not, cross-section information given by the program may be very confusing.

Another problem is that often the matrix elements really are valid only in the resonance region. The reason is that one usually includes only the simplest -channel graph in the calculation. It is this `signal' graph that has a peak at the position of the resonance, where it (usually) gives much larger cross sections than the other `background' graphs. Away from the resonance position, `signal' and `background' may be of comparable order, or the `background' may even dominate. There is a quantum mechanical interference when some of the `signal' and `background' graphs have the same initial and final state, and this interference may be destructive or constructive. When the interference is non-negligible, it is no longer meaningful to speak of a `signal' cross section. As an example, consider the scattering of longitudinal 's, , where the `signal' process is -channel exchange of a Higgs. This graph by itself is ill-behaved away from the resonance region. Destructive interference with `background' graphs such as -channel exchange of a Higgs and - and -channel exchange of a is required to save unitarity at large energies.

In colliders, the parton distribution is peaked at rather than at . The situation therefore is the opposite, if one considers e.g. production in a machine running at energies above : the resonance-peak tail towards lower masses is suppressed and the one towards higher masses enhanced, with a sharp secondary peak at around the nominal energy of the machine. Also in this case, an appropriate definition of cross sections therefore is necessary -- with additional complications due to the interference between and . When other processes are considered, problems of interference with background appears also here. Numerically the problems may be less pressing, however, since the secondary peak is occurring in a high-mass region, rather than in a more complicated low-mass one. Further, in there is little uncertainty from the shape of the parton distributions.

In processes where a pair of resonances are produced, e.g.
, cross section are almost always given
in the zero-width approximation for the resonances. Here
two substitutions of the type

There is no particular reason why the full mass-dependence could not be introduced, except for the extra work and time consumption needed for each process. In fact, the matrix elements for several and production processes do contain the full expressions. On the other hand, the matrix elements given in the literature are often valid only when the resonances are almost on the mass shell, since some graphs have been omitted. As an example, the process is dominated by when each of the two lepton pairs is close to in mass, but in general also receives contributions e.g. from , followed by and . The latter contributions are neglected in cross sections given in the zero-width approximation.

Widths may induce gauge invariance problems, in particular when the -channel graph interferes with - or -channel ones. Then there may be an imperfect cancellation of contributions at high energies, leading to an incorrect cross section behaviour. The underlying reason is that a Breit-Wigner corresponds to a resummation of terms of different orders in coupling constants, and that therefore effectively the -channel contributions are calculated to higher orders than the - or -channel ones, including interference contributions. A specific example is , where -channel exchange interferes with -channel exchange. In such cases, a fixed width is used in the denominator. One could also introduce procedures whereby the width is made to vanish completely at high energies, and theoretically this is the cleanest, but the fixed-width approach appears good enough in practice.

Another gauge invariance issue is when two particles of the same kind
are produced in a pair, e.g.
. Matrix elements are
then often calculated for one common mass, even though in real
life the masses . The proper gauge invariant procedure
to handle this would be to study the full six-fermion state obtained
after the two
decays, but that may
be overkill if indeed the 's are close to mass shell. Even when only
equal-mass matrix elements are available, Breit-Wigners are therefore
used to select two separate masses and . From these two
masses, an average mass is constructed so that the
velocity factor of eq. () is retained,

(93) |

(94) |

Processes with one final-state resonance and another ordinary final-state product, e.g. , are treated in the same spirit as the processes with two resonances, except that only one mass need be selected according to a Breit-Wigner.