Next: Cross-section Calculations Up: Process Generation Previous: Kinematics and Cross Section   Contents

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

 (84)

Here the physics is contained in the cross section . The scale is usually taken to be .

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

 (85)

In this formula the resonance width is a constant.

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

 (86)

The factor is evaluated as a sum over all possible final-state channels, . Each decay channel may have its own dependence, as follows.

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)

Having done that, is allowed to run [Kle89], and is evaluated at the scale. Thereby the relevant electroweak loop correction factors are recovered at the scale. However, the option exists to go the other way and eliminate in favour of . Currently is not allowed to run.

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]

 (88)

Here is the input mass at a reference scale , defined in the scheme. Typical choices are either or ; the latter would be relevant if the reference scale is chosen at the threshold. Both and are as given in .

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

 (89)

As a simple example, the cross section for the process can be written as
 (90)

where
 (91)

If the effects of several initial and/or final states are studied, it is straightforward to introduce an appropriate summation in the numerator.

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

 (92)

are used to introduce mass distributions for the two resonance masses, i.e. and . In the formula, is the nominal mass and the actually selected one. The phase-space integral over , and in eq. () is then extended to involve also and . The effects of the mass-dependent width is only partly taken into account, by replacing the nominal masses and in the expression by the actually generated ones (also e.g. in the relation between and ), while the widths are evaluated at the nominal masses. This is the equivalent of a simple replacement of by in the numerator of eq. (), but not in the denominator. In addition, the full threshold dependence of the widths, i.e. the velocity-dependent factors, is not reproduced.

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)

This choice certainly is not unique, but normally should provide a sensible behaviour, also around threshold. Of course, the differential cross section is no longer guaranteed to be gauge invariant when gauge bosons are involved, or positive definite. The program automatically flags the latter situation as unphysical. The approach may well break down when either or both particles are far away from mass shell. Furthermore, the preliminary choice of scattering angle is also retained. Instead of the correct and of eq. (), modified
 (94)

can then be obtained. The , and are now used in the matrix elements to decide whether to retain the event or not.

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.

Next: Cross-section Calculations Up: Process Generation Previous: Kinematics and Cross Section   Contents
Stephen Mrenna 2007-10-30