MSEL = 1, 2
No higher-order processes are explicitly included, nor any higher-order loop corrections to the processes. However, by initial- and final-state QCD radiation, multijet events are being generated, starting from the above processes. The shower rate of multijet production is clearly uncertain by some amount, especially for well-separated jets.
A string-based fragmentation scheme such as the Lund model needs cross sections for the different colour flows; these have been calculated in [Ben84] and differ from the usual calculations by interference terms of the order . By default, the standard colour-summed QCD expressions for the differential cross sections are used. In this case, the interference terms are distributed among the various colour flows according to the pole structure of the terms. However, the interference terms can be excluded, by changing MSTP(34)
As an example, consider subprocess 28,
. The total
cross section for this process, obtained by summing and squaring the
Feynman -, -, and -channel graphs, is
On the other hand, the cross sections for the two possible colour
flows of this subprocess are [Ben84]
The sum of the and contributions is
All the matrix elements in this group are for massless quarks (although final-state quarks are of course put on the mass shell). As a consequence, cross sections are divergent for , and some kind of regularization is required. Normally you are expected to set the desired value in CKIN(3).
The new flavour produced in the annihilation processes (ISUB = 12 and 53) is determined by the flavours allowed for gluon splitting into quark-antiquark; see switch MDME.
Subprocess 96 is special among all the ones in the program. In terms of the basic cross section, it is equivalent to the sum of the other ones, i.e. 11, 12, 13, 28, 53 and 68. The phase space is mapped differently, however, and allows as input variable. This is especially useful in the context of the multiple interactions machinery, see subsection , where potential scatterings are considered in order of decreasing , with a form factor related to the probability of not having another scattering with a larger than the considered one. You are not expected to access process 96 yourself. Instead it is automatically initialized and used either if process 95 is included or if multiple interactions are switched on. The process will then appear in the maximization information output, but not in the cross section table at the end of a run. Instead, the hardest scattering generated within the context of process 95 is reclassified as an event of the 11, 12, 13, 28, 53 or 68 kinds, based on the relative cross section for these in the point chosen. Further multiple interactions, subsequent to the hardest one, also do not show up in cross section tables.