For protons, many sets exist on the market. These are obtained by fits to experimental data, constrained so that the dependence is in accordance with the standard QCD evolution equations. The current default in PYTHIA is GRV 94L [Glü95], a simple leading-order fit. Several other sets are found in PYTHIA. The complete list is:
There is a steady flow of new parton-distribution sets on the market. To keep track of all of them is a major work on its own. Therefore PYTHIA contains an interface to an external library of parton distribution functions, PDFLIB [Plo93]. This is a truly encyclopedic collection of almost all proton, pion and photon parton distributions proposed since the late 70's. Three dummy routines come with the PYTHIA package, so as to avoid problems with unresolved external references if PDFLIB is not linked. One should also note that PYTHIA does not check the results, but assumes that sensible answers will be returned, also outside the nominal range of a set. Only the sets that come with PYTHIA have been suitably modified to provide reasonable answers outside their nominal domain of validity.
From the proton parton distributions, those of the neutron are obtained by isospin conjugation, i.e. and .
The program does allow for incoming beams of a number of hyperons: , , and . Here one has essentially no experimental information. One could imagine to construct models in which valence quarks are found at larger average values than valence and ones, because of the larger -quark mass. However, hyperon beams is a little-used part of the program, included only for a few specific studies. Therefore a simple approach has been taken, in which an average valence quark distribution is constructed as , according to which each valence quark in a hyperon is assumed to be distributed. Sea-quark and gluon distributions are taken as in the proton. Any proton parton distribution set may be used with this procedure.