The full machinery needed for a multiparton system is very complicated, and is described in detail in [Sjö84]. The following outline is far from complete, and is complicated nonetheless. The main message to be conveyed is that a Lorentz covariant algorithm exists for handling an arbitrary parton configuration, but that the necessary machinery is more complex than in either cluster or independent fragmentation.

Assume partons, with ordering along the string, and related four-momenta, given by . The initial string then contains separate pieces. The string piece between the quark and its neighbouring gluon is, in four-momentum space, spanned by one side with four-momentum and another with . The factor of 1/2 in the second expression comes from the fact that the gluon shares its energy between two string pieces. The indices `' and `' denotes direction towards the and end, respectively. The next string piece, counted from the quark end, is spanned by and , and so on, with the last one being and .

For the algorithm to work, it is important that all be light-cone-like, i.e. . Since gluons are massless, it is only the two endpoint quarks which can cause problems. The procedure here is to create new vectors for each of the two endpoint regions, defined to be linear combinations of the old ones for the same region, with coefficients determined so that the new vectors are light-cone-like. De facto, this corresponds to replacing a massive quark at the end of a string piece with a massless quark at the end of a somewhat longer string piece. With the exception of the added fictitious piece, which anyway ends up entirely within the heavy hadron produced from the heavy quark, the string motion remains unchanged by this.

In the continued string motion, when new string regions appear as time goes by, the first such string regions that appear can be represented as being spanned by one and another four-vector, with and not necessarily adjacent. For instance, in the case, the `third' string region is spanned by and . Later on in the string evolution history, it is also possible to have regions made up of two or two momenta. These appear when an endpoint quark has lost all its original momentum, has accreted the momentum of a gluon, and is now re-emitting this momentum. In practice, these regions may be neglected. Therefore only pieces made up by a pair of momenta are considered in the program.

The allowed string regions may be ordered in an abstract parameter
plane, where the indices of the four-momentum pairs define
the position of each region along the two (parameter plane)
coordinate axes. In this plane the fragmentation procedure can be
described as a sequence of steps, starting at the quark end,
where , and ending at the antiquark one,
. Each step is taken from an `old'
pair production vertex, to the
production vertex of a `new'
pair, and the string piece between these two string
breaks represent a hadron. Some steps may be taken within one and
the same region, while others may have one vertex in one region
and the other vertex in another region. Consistency requirements,
like energy-momentum conservation, dictates that vertex and
region values be ordered in a monotonic sequence, and that
the vertex positions are monotonically ordered inside each region.
The four-momentum of each hadron can be read off, for ()
momenta, by projecting the separation between the old and the new
vertex on to the () axis. If the four-momentum fraction of
taken by a hadron is denoted ,
then the total hadron four-momentum is given by

The pairs are the transverse momenta produced at
the two string breaks, and the
pairs
four-vectors transverse to the string directions in the regions
of the respective string breaks:

(255) |

The fact that the hadron should be on mass shell, , puts
one constraint on where a new breakup may be, given that the old
one is already known, just as eq. () did in the
simple 2-jet case. The remaining degree of freedom is, as before,
to be given by
the fragmentation function . The interpretation of the
is only well-defined for a step entirely constrained to one of the
initial string regions, however, which is not enough. In the
2-jet case, the values can be related to the proper times
of string breaks, as follows. The variable
, with the string tension and
the proper time between the production vertex of the
partons and the breakup point, obeys an iterative relation of the
kind

(256) |

Here represents the value at the and endpoints, and and the values at the old and new breakup vertices needed to produce a hadron with transverse mass , and with the of the step chosen according to . The proper time can be defined in an unambiguous way, also over boundaries between the different string regions, so for multijet events the variable may be interpreted just as an auxiliary variable needed to determine the next value. (In the Lund symmetric fragmentation function derivation, the variable actually does appear naturally, so the choice is not as arbitrary as it may seem here.) The mass and constraints together are sufficient to determine where the next string breakup is to be chosen, given the preceding one in the iteration scheme. In reality, several ambiguities remain, however none of these are of importance to the overall picture.

The algorithm for finding the next breakup then works something
like follows. Pick a hadron, , and , and calculate the next
. If the old breakup is in the region , and if the
new breakup is also assumed to be in the same region, then the
and constraints can be reformulated in terms of
the fractions and the hadron must take of the
total four-vectors and :

Here the coefficients are fairly simple expressions, obtainable by squaring eq. (), while are slightly more complicated in that they depend on the position of the old string break, but both the and the are explicitly calculable. What remains is an equation system with two unknowns, and . The absence of any quadratic terms is due to the fact that all , i.e. to the choice of a formulation based on light-cone-like longitudinal vectors. Of the two possible solutions to the equation system (elimination of one variable gives a second degree equation in the other), one is unphysical and can be discarded outright. The other solution is checked for whether the values are actually inside the physically allowed region, i.e. whether the values of the current step, plus whatever has already been used up in previous steps, are less than unity. If yes, a solution has been found. If no, it is because the breakup could not take place inside the region studied, i.e. because the equation system was solved for the wrong region. One therefore has to change either index or index above by one step, i.e. go to the next nearest string region. In this new region, a new equation system of the type in eq. () may be written down, with new coefficients. A new solution is found and tested, and so on until a physically acceptable solution is found. The hadron four-momentum is now given by an expression of the type (). The breakup found forms the starting point for the new step in the fragmentation chain, and so on. The final joining in the middle is done as in the 2-jet case, with minor extensions.