Tunes
Since some physics aspects cannot be derived from first principles,
this program contains many parameters that represent a true
uncertainty in our understanding of nature. Particularly afflicted
are the areas of hadronization and multiple interactions, which both
involve nonperturbative QCD physics.
Technically, PYTHIA parameters can be varied independently of each
other, but the physical requirement of a sensible description of a set
of data leads to correlations and anticorrelations between the
parameters. Hence the need to produce tunes, not of one parameter at
a time, but simultaneously for a group of them. A well-known such
example is parton densities, where combined tunes to a wide range of data
have been produced, that can then be obtained prepackaged.
Given the many PYTHIA parameters to be tuned, it is convenient to
divide the task into subtasks. Firstly, if we assume jet universality,
hadronization and final-state parton showers should be tuned to
e^+e^- annihilation data, notably from LEP1, since this
offers the cleanest environment. Secondly, with such parameters fixed,
hadron collider data should be studied to pin down multiple interactions
and other further aspects, such as initial-state radiation. Ideally this
would be done separately for diffractive and non-diffractive events,
although it is not possible to have a clean separation. (Thirdly
would come anything else, such as physics with photon beams, which
involve further parameters, but that is beyond the current scope.)
The first step in this program has now been taken, with a tune to LEP1
data by Hendrik Hoeth, using the Rivet + Professor framework. Starting
with version 8.125 it defines the default values for hadronization
parameters and timelike showers.
The situation is worse for multiple interactions, where PYTHIA 8 is more
different from PYTHIA 6. Nevertheless, a first simple tune is now
available, appropriately called "Tune 1", and became default starting with
version 8.127.
It was noted, in particular by Hendrik Hoeth, that the
program had a tension between parameters needed to describe minimum-bias
and underlying-event activity. Therefore some further physics features
have been introduced in the code itself, which are made default as of
8.140. This version also includes two new draft tunes, 2C and 2M, based
on the CTEQ 6L1 and the MRST LO** PDF sets, respectively. These have been
made by hand, as a prequel to complete Professor-style tunings.
The very first data to come out of the LHC shows a higher rapidity plateau
than predicted for current PYTHIA 6 tunes, also for the lower energies.
This may suggest some tension in the data. Therefore two alternatives,
3C and 3M, have been produced by a few brute-force changes of 2C and 2M.
This includes a reduced admixture of diffractive topologies and a steeper
multiplicity variation with energy.
In the future we hope to see further PYTHIA 8 tunes appear. Like with
parton distributions, there is likely to be several tunes, because
different sets of data will pull in different directions, by imperfections
in the model or in the data, and by differences in the chosen
tuning strategies. We therefore propose to collect some of these tunes
here, in a prepackaged form. Of course, in all cases it is a matter
of setting values for parameters already defined elsewhere, so the
tunes offer no new functionality, only a more convenient setup.
You should be aware that the evolution of the program will not guarantee
complete backwards compatibility between versions. Most obviously this
concerns bug fixes. But also for some other major changes, like the
introduction of the new diffractive machinery, the default behaviour
of old tunes has been changed retroactively. (Which should be fine for
diffraction, since previous tunes were not based on data strongly
influenced by diffraction.)
If you set either the Tune:ee
and Tune:pp
modes below non-zero then all parameters used in the respective tune
will be set accordingly when pythia.init(...)
is called.
You can check this by calling pythia.settings.listChanged()
before and after initialization; before only the tune modes are
nondefault, afterwards all the non-default-valued parameters in the
tune appear. Therefore you cannot directly combine a tune option with
your own choices for some of the parameters used in the tune, since
the values you set before pythia.init(...)
would be
overwritten at that point.
Instead you can yourself call the method that sets up the tunes,
pythia.initTunes(int eeTune, int ppTune)
.
Here eeTune
and ppTune
encode the same
options as the already-introduced Tune:ee
and
Tune:pp
. So the procedure would be first to call
pythia.initTunes(...)
, then do the further changes
you wish, and finally initialize with Tune:ee
and
Tune:pp
both zero.
mode
Tune:ee
(default = 0
; minimum = 0
; maximum = 3
)
Choice of tune to e^+e^- data, mainly for the hadronization
and timelike-showering aspects of PYTHIA.
option
0 : no values are overwritten at initialization,
so you can set the individual parameters as you wish.
option
1 : the original PYTHIA 8 parameter set, based on some
very old flavour studies (with JETSET around 1990) and a simple tune
of alpha_strong to three-jet shapes to the new
pT-ordered shower. These were the default values before
version 8.125.
option
2 : a tune by Marc Montull to the LEP 1 particle
composition, as published in the RPP (August 2007). No related (re)tune
to event shapes has been performed, however.
option
3 : a tune to a wide selection of LEP1 data by Hendrik
Hoeth within the Rivet + Professor framework, both to hadronization and
timelike-shower parameters (June 2009). These are the default values
starting from version 8.125, so currently there is no need for this
option.
mode
Tune:pp
(default = 0
; minimum = 0
; maximum = 6
)
Choice of tune to pp / ppbar data, mainly for the
initial-state-radiation, multiple-interactions and beam-remnants
aspects of PYTHIA.
option
0 : no values are overwritten at initialization,
so you can set the individual parameters as you wish. Most default
values are based on "Tune 1", option 2 below, but some new options
introduced in 8.140 means that the two no longer agree.
option
1 : default used up to version 8.126, based on
some early and primitive comparisons with data.
option
2 : "Tune 1", default in 8.127 - 8.139, based on some
data comparisons by Peter Skands. Largely but not wholly overlaps
with the default option 0.
option
3 : "Tune 2C", new draft tune, introduced with 8.140.
It uses the CTEQ 6L1 PDF, and is intended to give good agreement with
much of the published CDF data.
option
4 : "Tune 2M", new draft tune, introduced with 8.140.
It is uses the MRST LO** PDF, which has a momentum sum somewhat above
unity, which is compensated by a smaller alpha_s than in the
previous tune. Again it is intended to give good agreement with much of
the published CDF data.
option
5 : "Tune 3C", new draft tune, introduced with 8.140.
Starts out from tune 2C, but with a reduced cross section for
diffraction, plus modified multiple interactions parameters to give
a higher and more rapidly increasing charged pseudorapidity plateau,
for better agreement with some early key LHC numbers.
option
6 : "Tune 3M" tune, introduced with 8.140.
Starts out from tune 2M and then does the same kind of modifications
as for 3C.