herwig is hosted by Hepforge, IPPP Durham
Herwig  7.1.5
Todo List
Member Herwig::ClusterFissioner::cut (stack< ClusterPtr > &, ClusterVector &, tPVector &finalhadrons, bool softUEisOn)
is the following still true? For normal, non-beam clusters, a power-like mass distribution is used, whereas for beam clusters a fast-decreasing exponential mass distribution is used instead. This avoids many iterative splitting which could produce an unphysical large transverse energy from a supposed soft beam remnant process.
Member Herwig::HadronSelector::lightestHadronPair (tcPDPtr ptr1, tcPDPtr ptr2, tcPDPtr ptr3=PDPtr()) const

At the moment it does nothing in the case that also par3 is present.

The method assumes par3 == PDPtr() (otherwise we don't know how to proceed: a possible, trivial way would be to randomly select two of the three (anti-)quarks and treat them as a (anti-)diquark, reducing the problem to two components as treated below. In the normal (two components) situation, the strategy is the following: treat in the same way the two possibilities: (d dbar) (i=0) and (u ubar) (i=1) as the pair quark-antiquark necessary to form a pair of hadrons containing the input flavour par1 and par2; finally, select the one that produces the lightest pair of hadrons, compatible with the charge conservation constraint.