/usr/share/doc/pythia8-doc/html/BeamRemnants.html is in pythia8-doc-html 8.1.80-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 | <html>
<head>
<title>Beam Remnants</title>
<link rel="stylesheet" type="text/css" href="pythia.css"/>
<link rel="shortcut icon" href="pythia32.gif"/>
</head>
<body>
<h2>Beam Remnants</h2>
<h3>Introduction</h3>
The <code>BeamParticle</code> class contains information on all partons
extracted from a beam (so far). As each consecutive multiparton interaction
defines its respective incoming parton to the hard scattering a
new slot is added to the list. This information is modified when
the backwards evolution of the spacelike shower defines a new
initiator parton. It is used, both for the multiparton interactions
and the spacelike showers, to define rescaled parton densities based
on the <i>x</i> and flavours already extracted, and to distinguish
between valence, sea and companion quarks. Once the perturbative
evolution is finished, further beam remnants are added to obtain a
consistent set of flavours. The current physics framework is further
described in [<a href="Bibliography.html" target="page">Sjo04</a>].
<p/>
The introduction of <a href="MultipartonInteractions.html" target="page">rescattering</a>
in the multiparton interactions framework further complicates the
processing of events. Specifically, when combined with showers,
the momentum of an individual parton is no longer uniquely associated
with one single subcollision. Nevertheless the parton is classified
with one system, owing to the technical and administrative complications
of more complete classifications. Therefore the addition of primordial
<i>kT</i> to the subsystem initiator partons does not automatically
guarantee overall <i>pT</i> conservation. Various tricks are used to
minimize the mismatch, with a brute force shift of all parton
<i>pT</i>'s as a final step.
<p/>
Much of the above information is stored in a vector of
<code>ResolvedParton</code> objects, which each contains flavour and
momentum information, as well as valence/companion information and more.
The <code>BeamParticle</code> method <code>list()</code> shows the
contents of this vector, mainly for debug purposes.
<p/>
The <code>BeamRemnants</code> class takes over for the final step
of adding primordial <i>kT</i> to the initiators and remnants,
assigning the relative longitudinal momentum sharing among the
remnants, and constructing the overall kinematics and colour flow.
This step couples the two sides of an event, and could therefore
not be covered in the <code>BeamParticle</code> class, which only
considers one beam at a time.
<p/>
The methods of these classes are not intended for general use,
and so are not described here.
<p/>
In addition to the parameters described on this page, note that the
choice of <a href="PDFSelection.html" target="page">parton densities</a> is made
in the <code>Pythia</code> class. Then pointers to the pdf's are handed
on to <code>BeamParticle</code> at initialization, for all subsequent
usage.
<h3>Primordial <i>kT</i></h3>
The primordial <i>kT</i> of initiators of hard-scattering subsystems
are selected according to Gaussian distributions in <i>p_x</i> and
<i>p_y</i> separately. The widths of these distributions are chosen
to be dependent on the hard scale of the central process and on the mass
of the whole subsystem defined by the two initiators:
<br/><i>
sigma = (sigma_soft * Q_half + sigma_hard * Q) / (Q_half + Q)
* m / (m_half + m)
</i><br/>
Here <i>Q</i> is the hard-process renormalization scale for the
hardest process and the <i>pT</i> scale for subsequent multiparton
interactions, <i>m</i> the mass of the system, and
<i>sigma_soft</i>, <i>sigma_hard</i>, <i>Q_half</i> and
<i>m_half</i> parameters defined below. Furthermore each separately
defined beam remnant has a distribution of width <i>sigma_remn</i>,
independently of kinematical variables.
<p/><code>flag </code><strong> BeamRemnants:primordialKT </strong>
(<code>default = <strong>on</strong></code>)<br/>
Allow or not selection of primordial <i>kT</i> according to the
parameter values below.
<p/><code>parm </code><strong> BeamRemnants:primordialKTsoft </strong>
(<code>default = <strong>0.5</strong></code>; <code>minimum = 0.</code>)<br/>
The width <i>sigma_soft</i> in the above equation, assigned as a
primordial <i>kT</i> to initiators in the soft-interaction limit.
<p/><code>parm </code><strong> BeamRemnants:primordialKThard </strong>
(<code>default = <strong>2.0</strong></code>; <code>minimum = 0.</code>)<br/>
The width <i>sigma_hard</i> in the above equation, assigned as a
primordial <i>kT</i> to initiators in the hard-interaction limit.
<p/><code>parm </code><strong> BeamRemnants:halfScaleForKT </strong>
(<code>default = <strong>1.</strong></code>; <code>minimum = 0.</code>)<br/>
The scale <i>Q_half</i> in the equation above, defining the
half-way point between hard and soft interactions.
<p/><code>parm </code><strong> BeamRemnants:halfMassForKT </strong>
(<code>default = <strong>1.</strong></code>; <code>minimum = 0.</code>)<br/>
The scale <i>m_half</i> in the equation above, defining the
half-way point between low-mass and high-mass subsystems.
(Kinematics construction can easily fail if a system is assigned
a primordial <i>kT</i> value higher than its mass, so the
mass-dampening is intended to reduce some troubles later on.)
<p/><code>parm </code><strong> BeamRemnants:primordialKTremnant </strong>
(<code>default = <strong>0.4</strong></code>; <code>minimum = 0.</code>)<br/>
The width <i>sigma_remn</i>, assigned as a primordial <i>kT</i>
to beam-remnant partons.
<p/>
A net <i>kT</i> imbalance is obtained from the vector sum of the
primordial <i>kT</i> values of all initiators and all beam remnants.
This quantity is compensated by a shift shared equally between
all partons, except that the dampening factor <i>m / (m_half + m)</i>
is again used to suppress the role of small-mass systems.
<p/>
Note that the current <i>sigma</i> definition implies that
<i><pT^2> = <p_x^2>+ <p_y^2> = 2 sigma^2</i>.
It thus cannot be compared directly with the <i>sigma</i>
of nonperturbative hadronization, where each quark-antiquark
breakup corresponds to <i><pT^2> = sigma^2</i> and only
for hadrons it holds that <i><pT^2> = 2 sigma^2</i>.
The comparison is further complicated by the reduction of
primordial <i>kT</i> values by the overall compensation mechanism.
<p/><code>flag </code><strong> BeamRemnants:rescatterRestoreY </strong>
(<code>default = <strong>off</strong></code>)<br/>
Is only relevant when <a href="MultipartonInteractions.html" target="page">rescattering</a>
is switched on in the multiparton interactions scenario. For a normal
interaction the rapidity and mass of a system is preserved when
primordial <i>kT</i> is introduced, by appropriate modification of the
incoming parton momenta. Kinematics construction is more complicated for
a rescattering, and two options are offered. Differences between these
can be used to explore systematic uncertainties in the rescattering
framework.<br/>
The default behaviour is to keep the incoming rescattered parton as is,
but to modify the unrescattered incoming parton so as to preserve the
invariant mass of the system. Thereby the rapidity of the rescattering
is modified.<br/>
The alternative is to retain the rapidity (and mass) of the rescattered
system when primordial <i>kT</i> is introduced. This is made at the
expense of a modified longitudinal momentum of the incoming rescattered
parton, so that it does not agree with the momentum it ought to have had
by the kinematics of the previous interaction.<br/>
For a double rescattering, when both incoming partons have already scattered,
there is no obvious way to retain the invariant mass of the system in the
first approach, so the second is always used.
<h3>Colour flow</h3>
The colour flows in the separate subprocesses defined in the
multiparton-interactions scenario are tied together via the assignment
of colour flow in the beam remnant. This is not an unambiguous
procedure, but currently no parameters are directly associated with it.
However, a simple "minimal" procedure of colour flow only via the beam
remnants does not result in a scenario in
agreement with data, notably not a sufficiently steep rise of
<i><pT>(n_ch)</i>. The true origin of this behaviour and the
correct mechanism to reproduce it remains one of the big unsolved issues
at the borderline between perturbative and nonperturbative QCD.
As a simple attempt, an additional step is introduced, wherein the gluons
of a lower-<i>pT</i> system are merged with the ones in a higher-pT one.
<p/><code>flag </code><strong> BeamRemnants:reconnectColours </strong>
(<code>default = <strong>on</strong></code>)<br/>
Allow or not a system to be merged with another one.
<p/><code>parm </code><strong> BeamRemnants:reconnectRange </strong>
(<code>default = <strong>10.0</strong></code>; <code>minimum = 0.</code>; <code>maximum = 10.</code>)<br/>
A system with a hard scale <i>pT</i> can be merged with one of a
harder scale with a probability that is
<i>pT0_Rec^2 / (pT0_Rec^2 + pT^2)</i>, where
<i>pT0_Rec</i> is <code>reconnectRange</code> times <i>pT0</i>,
the latter being the same energy-dependent dampening parameter as
used for multiparton interactions.
Thus it is easy to merge a low-<i>pT</i> system with any other,
but difficult to merge two high-<i>pT</i> ones with each other.
<p/>
The procedure is used iteratively. Thus first the reconnection probability
<i>P = pT0_Rec^2 / (pT0_Rec^2 + pT^2)</i> of the lowest-<i>pT</i>
system is found, and gives the probability for merger with the
second-lowest one. If not merged, it is tested with the third-lowest one,
and so on. For the <i>m</i>'th higher system the reconnection
probability thus becomes <i>(1 - P)^(m-1) P</i>. That is, there is
no explicit dependence on the higher <i>pT</i> scale, but implicitly
there is via the survival probability of not already having been merged
with a lower-<i>pT</i> system. Also note that the total reconnection
probability for the lowest-<i>pT</i> system in an event with <i>n</i>
systems becomes <i>1 - (1 - P)^(n-1)</i>. Once the fate of the
lowest-<i>pT</i> system has been decided, the second-lowest is considered
with respect to the ones above it, then the third-lowest, and so on.
<p/>
Once it has been decided which systems should be joined, the actual merging
is carried out in the opposite direction. That is, first the hardest
system is studied, and all colour dipoles in it are found (including to
the beam remnants, as defined by the holes of the incoming partons).
Next each softer system to be merged is studied in turn. Its gluons are,
in decreasing <i>pT</i> order, inserted on the colour dipole <i>i,j</i>
that gives the smallest <i>(p_g p_i)(p_g p_j)/(p_i p_j)</i>, i.e.
minimizes the "disturbance" on the existing dipole, in terms of
<i>pT^2</i> or <i>Lambda</i> measure (string length). The insertion
of the gluon means that the old dipole is replaced by two new ones.
Also the (rather few) quark-antiquark pairs that can be traced back to
a gluon splitting are treated in close analogy with the gluon case.
Quark lines that attach directly to the beam remnants cannot be merged
but are left behind.
<p/>
The joining procedure can be viewed as a more sophisticated variant of
the one introduced already in [<a href="Bibliography.html" target="page">Sjo87</a>]. Clearly it is ad hoc.
It hopefully captures some elements of truth. The lower <i>pT</i> scale
a system has the larger its spatial extent and therefore the larger its
overlap with other systems. It could be argued that one should classify
individual initial-state partons by <i>pT</i> rather than the system
as a whole. However, for final-state radiation, a soft gluon radiated off
a hard parton is actually produced at late times and therefore probably
less likely to reconnect. In the balance, a classification by system
<i>pT</i> scale appears sensible as a first try.
<p/>
Note that the reconnection is carried out before resonance decays are
considered. Colour inside a resonance therefore is not reconnected.
This is a deliberate choice, but certainly open to discussion and
extensions at a later stage, as is the rest of this procedure.
<h3>Further variables</h3>
<p/><code>mode </code><strong> BeamRemnants:maxValQuark </strong>
(<code>default = <strong>3</strong></code>; <code>minimum = 0</code>; <code>maximum = 5</code>)<br/>
The maximum valence quark kind allowed in acceptable incoming beams,
for which multiparton interactions are simulated. Default is that hadrons
may contain <i>u</i>, <i>d</i> and <i>s</i> quarks,
but not <i>c</i> and <i>b</i> ones, since sensible
kinematics has not really been worked out for the latter.
<p/><code>mode </code><strong> BeamRemnants:companionPower </strong>
(<code>default = <strong>4</strong></code>; <code>minimum = 0</code>; <code>maximum = 4</code>)<br/>
When a sea quark has been found, a companion antisea quark ought to be
nearby in <i>x</i>. The shape of this distribution can be derived
from the gluon mother distribution convoluted with the
<i>g -> q qbar</i> splitting kernel. In practice, simple solutions
are only feasible if the gluon shape is assumed to be of the form
<i>g(x) ~ (1 - x)^p / x</i>, where <i>p</i> is an integer power,
the parameter above. Allowed values correspond to the cases programmed.
<br/>
Since the whole framework is approximate anyway, this should be good
enough. Note that companions typically are found at small <i>Q^2</i>,
if at all, so the form is supposed to represent <i>g(x)</i> at small
<i>Q^2</i> scales, close to the lower cutoff for multiparton interactions.
<p/>
When assigning relative momentum fractions to beam-remnant partons,
valence quarks are chosen according to a distribution like
<i>(1 - x)^power / sqrt(x)</i>. This <i>power</i> is given below
for quarks in mesons, and separately for <i>u</i> and <i>d</i>
quarks in the proton, based on the approximate shape of low-<i>Q^2</i>
parton densities. The power for other baryons is derived from the
proton ones, by an appropriate mixing. The <i>x</i> of a diquark
is chosen as the sum of its two constituent <i>x</i> values, and can
thus be above unity. (A common rescaling of all remnant partons and
particles will fix that.) An additional enhancement of the diquark
momentum is obtained by its <i>x</i> value being rescaled by the
<code>valenceDiqEnhance</code> factor.
<p/><code>parm </code><strong> BeamRemnants:valencePowerMeson </strong>
(<code>default = <strong>0.8</strong></code>; <code>minimum = 0.</code>)<br/>
The abovementioned power for valence quarks in mesons.
<p/><code>parm </code><strong> BeamRemnants:valencePowerUinP </strong>
(<code>default = <strong>3.5</strong></code>; <code>minimum = 0.</code>)<br/>
The abovementioned power for valence <i>u</i> quarks in protons.
<p/><code>parm </code><strong> BeamRemnants:valencePowerDinP </strong>
(<code>default = <strong>2.0</strong></code>; <code>minimum = 0.</code>)<br/>
The abovementioned power for valence <i>d</i> quarks in protons.
<p/><code>parm </code><strong> BeamRemnants:valenceDiqEnhance </strong>
(<code>default = <strong>2.0</strong></code>; <code>minimum = 0.5</code>; <code>maximum = 10.</code>)<br/>
Enhancement factor for valence diquarks in baryons, relative to the
simple sum of the two constituent quarks.
<p/><code>flag </code><strong> BeamRemnants:allowJunction </strong>
(<code>default = <strong>on</strong></code>)<br/>
The <code>off</code> option is intended for debug purposes only, as
follows. When more than one valence quark is kicked out of a baryon
beam, as part of the multiparton interactions scenario, the subsequent
hadronization is described in terms of a junction string topology.
This description involves a number of technical complications that
may make the program more unstable. As an alternative, by switching
this option off, junction configurations are rejected (which gives
an error message that the remnant flavour setup failed), and the
multiparton interactions and showers are redone until a
junction-free topology is found.
</body>
</html>
<!-- Copyright (C) 2013 Torbjorn Sjostrand -->
|