This file is indexed.

/usr/share/doc/pythia8-doc/html/MultipartonInteractions.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
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
<html>
<head>
<title>Multiparton Interactions</title>
<link rel="stylesheet" type="text/css" href="pythia.css"/>
<link rel="shortcut icon" href="pythia32.gif"/>
</head>
<body>

<h2>Multiparton Interactions</h2>

The starting point for the multiparton interactions physics scenario in
PYTHIA is provided by [<a href="Bibliography.html" target="page">Sjo87</a>]. Recent developments have 
included a more careful study of flavour and colour correlations, 
junction topologies and the relationship to beam remnants 
[<a href="Bibliography.html" target="page">Sjo04</a>], interleaving with initial-state radiation 
[<a href="Bibliography.html" target="page">Sjo05</a>], making use of transverse-momentum-ordered
initial- and final-state showers, with the extension to fully 
interleaved evolution covered in [<a href="Bibliography.html" target="page">Cor10a</a>]. A framework to 
handle rescattering is described in [<a href="Bibliography.html" target="page">Cor09</a>].

<p/>
A big unsolved issue is how the colour of all these subsystems is 
correlated. For sure there is a correlation coming from the colour
singlet nature of the incoming beams, but in addition final-state 
colour rearrangements may change the picture. Indeed such extra
effects appear necessary to describe data, e.g. on 
<i>&lt;pT&gt;(n_ch)</i>. A simple implementation of colour 
rearrangement is found as part of the
<a href="BeamRemnants.html" target="page">beam remnants</a> description.

<h3>Main variables</h3>

<h4>Matching to hard process</h4>

The maximum <i>pT</i> to be allowed for multiparton interactions is
related to the nature of the hard process itself. It involves a
delicate balance between not double-counting and not leaving any
gaps in the coverage. The best procedure may depend on information 
only the user has: how the events were generated and mixed (e.g. with 
Les Houches Accord external input), and how they are intended to be 
used. Therefore a few options are available, with a sensible default 
behaviour.
<p/><code>mode&nbsp; </code><strong> MultipartonInteractions:pTmaxMatch &nbsp;</strong> 
 (<code>default = <strong>0</strong></code>; <code>minimum = 0</code>; <code>maximum = 2</code>)<br/>
Way in which the maximum scale for multiparton interactions is set
to match the scale of the hard process itself.
<br/><code>option </code><strong> 0</strong> : <b>(i)</b> if the final state of the hard process 
(not counting subsequent resonance decays) contains only quarks 
(<i>u, d, s, c, b</i>), gluons and photons then <i>pT_max</i> 
is chosen to be the factorization scale for internal processes 
and the <code>scale</code> value for Les Houches input; 
<b>(ii)</b> if not, interactions are allowed to go all the way up 
to the kinematical limit. 
The reasoning is that the former kind of processes are generated by
the multiparton-interactions machinery and so would double-count hard
processes if allowed to overlap the same <i>pT</i> range, 
while no such danger exists in the latter case.
  
<br/><code>option </code><strong> 1</strong> : always use the factorization scale for an internal
process and the <code>scale</code> value for Les Houches input, 
i.e. the lower value. This should avoid double-counting, but
may leave out some interactions that ought to have been simulated.
  
<br/><code>option </code><strong> 2</strong> : always allow multiparton interactions up to the 
kinematical limit. This will simulate all possible event topologies, 
but may lead to double-counting.
  
<br/><b>Note:</b> If a "second hard" process is present, the two 
are analyzed separately for the default 0 option. It is enough that 
one of them only consists of quarks, gluons and photons to restrict
the <i>pT</i> range. The maximum for MPI is then set by the hard
interaction with lowest scale.
  

<h4>Cross-section parameters</h4>

The rate of interactions is determined by 
<p/><code>parm&nbsp; </code><strong> MultipartonInteractions:alphaSvalue &nbsp;</strong> 
 (<code>default = <strong>0.127</strong></code>; <code>minimum = 0.06</code>; <code>maximum = 0.25</code>)<br/>
The value of <i>alpha_strong</i> at <i>m_Z</i>. Default value is 
picked equal to the one used in CTEQ 5L. 
   

<p/>
The actual value is then regulated by the running to the scale 
<i>pT^2</i>, at which it is evaluated
<p/><code>mode&nbsp; </code><strong> MultipartonInteractions:alphaSorder &nbsp;</strong> 
 (<code>default = <strong>1</strong></code>; <code>minimum = 0</code>; <code>maximum = 2</code>)<br/>
The order at which <i>alpha_strong</i> runs at scales away from 
<i>m_Z</i>.
<br/><code>option </code><strong> 0</strong> : zeroth order, i.e. <i>alpha_strong</i> is kept 
fixed.  
<br/><code>option </code><strong> 1</strong> : first order, which is the normal value.  
<br/><code>option </code><strong> 2</strong> : second order. Since other parts of the code do 
not go to second order there is no strong reason to use this option, 
but there is also nothing wrong with it.  
  

<p/>
QED interactions are regulated by the <i>alpha_electromagnetic</i>
value at the <i>pT^2</i> scale of an interaction.
 
<p/><code>mode&nbsp; </code><strong> MultipartonInteractions:alphaEMorder &nbsp;</strong> 
 (<code>default = <strong>1</strong></code>; <code>minimum = -1</code>; <code>maximum = 1</code>)<br/>
The running of <i>alpha_em</i> used in hard processes.
<br/><code>option </code><strong> 1</strong> : first-order running, constrained to agree with
<code>StandardModel:alphaEMmZ</code> at the <i>Z^0</i> mass.
  
<br/><code>option </code><strong> 0</strong> : zeroth order, i.e. <i>alpha_em</i> is kept 
fixed at its value at vanishing momentum transfer.  
<br/><code>option </code><strong> -1</strong> : zeroth order, i.e. <i>alpha_em</i> is kept 
fixed, but at <code>StandardModel:alphaEMmZ</code>, i.e. its value
at the <i>Z^0</i> mass.
   
  

<p/>
Note that the choices of <i>alpha_strong</i> and <i>alpha_em</i> 
made here override the ones implemented in the normal process machinery, 
but only for the interactions generated by the 
<code>MultipartonInteractions</code> class.

<p/>
In addition there is the possibility of a global rescaling of 
cross sections (which could not easily be accommodated by a 
changed <i>alpha_strong</i>, since <i>alpha_strong</i> runs)
<p/><code>parm&nbsp; </code><strong> MultipartonInteractions:Kfactor &nbsp;</strong> 
 (<code>default = <strong>1.0</strong></code>; <code>minimum = 0.5</code>; <code>maximum = 4.0</code>)<br/>
Multiply all cross sections by this fix factor.
  

<p/> 
The processes used to generate multiparton interactions form a subset
of the standard library of hard processes. The input is slightly
different from the standard hard-process machinery, however, 
since incoming flavours, the <i>alpha_strong</i> value and most
of the kinematics are already fixed when the process is called.
It is possible to regulate the set of processes that are included in the
multiparton-interactions framework.

<p/><code>mode&nbsp; </code><strong> MultipartonInteractions:processLevel &nbsp;</strong> 
 (<code>default = <strong>3</strong></code>; <code>minimum = 0</code>; <code>maximum = 3</code>)<br/>
Set of processes included in the machinery.
<br/><code>option </code><strong> 0</strong> : only the simplest <i>2 -> 2</i> QCD processes
between quarks and gluons, giving no new flavours, i.e. dominated by
<i>t</i>-channel gluon exchange.  
<br/><code>option </code><strong> 1</strong> : also <i>2 -> 2</i> QCD processes giving new flavours
(including charm and bottom), i.e. proceeding through <i>s</i>-channel 
gluon exchange.  
<br/><code>option </code><strong> 2</strong> : also <i>2 -> 2</i> processes involving one or two
photons in the final state, <i>s</i>-channel <i>gamma</i>
boson exchange and <i>t</i>-channel <i>gamma/Z^0/W^+-</i>
boson exchange.  
<br/><code>option </code><strong> 3</strong> : also charmonium and bottomonium production, via
colour singlet and colour octet channels.   
  

<h4>Cross-section regularization</h4>

There are two complementary ways of regularizing the small-<i>pT</i> 
divergence, a sharp cutoff and a smooth dampening. These can be 
combined as desired, but it makes sense to coordinate with how the 
same issue is handled in <a href="SpacelikeShowers.html" target="page">spacelike 
showers</a>. Actually, by default, the parameters defined here are 
used also for the spacelike showers, but this can be overridden.

<p/>
Regularization of the divergence of the QCD cross section for 
<i>pT -> 0</i> is obtained by a factor <i>pT^4 / (pT0^2 + pT^2)^2</i>, 
and by using an <i>alpha_s(pT0^2 + pT^2)</i>. An energy dependence 
of the <i>pT0</i> choice is introduced by two further parameters, 
so that <i>pT0Ref</i> is the <i>pT0</i> value for the reference 
CM energy, <i>pT0Ref = pT0(ecmRef)</i>.   
<br/><b>Warning:</b> if a large <i>pT0</i> is picked for multiparton 
interactions, such that the integrated interaction cross section is 
below the nondiffractive inelastic one, this <i>pT0</i> will 
automatically be scaled down to cope.

<p/>
The actual <i>pT0</i> parameter used at a given CM energy scale, 
<i>ecmNow</i>, is obtained as
<br/><i>
     pT0 = pT0(ecmNow) = pT0Ref * (ecmNow / ecmRef)^ecmPow 
</i><br/>
where <i>pT0Ref</i>, <i>ecmRef</i> and <i>ecmPow</i> are the 
three parameters below.

<p/><code>parm&nbsp; </code><strong> MultipartonInteractions:pT0Ref &nbsp;</strong> 
 (<code>default = <strong>2.15</strong></code>; <code>minimum = 0.5</code>; <code>maximum = 10.0</code>)<br/>
The <i>pT0Ref</i> scale in the above formula.
<br/><b>Note:</b> <i>pT0Ref</i> is one of the key parameters in a
complete PYTHIA tune. Its value is intimately tied to a number of other
choices, such as that of colour flow description, so unfortunately it is
difficult to give an independent meaning to <i>pT0Ref</i>.
  

<p/><code>parm&nbsp; </code><strong> MultipartonInteractions:ecmRef &nbsp;</strong> 
 (<code>default = <strong>1800.0</strong></code>; <code>minimum = 1.</code>)<br/>
The <i>ecmRef</i> reference energy scale introduced above.
  

<p/><code>parm&nbsp; </code><strong> MultipartonInteractions:ecmPow &nbsp;</strong> 
 (<code>default = <strong>0.24</strong></code>; <code>minimum = 0.0</code>; <code>maximum = 0.5</code>)<br/>
The <i>ecmPow</i> energy rescaling pace introduced above.
  

<p/>
Alternatively, or in combination, a sharp cut can be used.
<p/><code>parm&nbsp; </code><strong> MultipartonInteractions:pTmin &nbsp;</strong> 
 (<code>default = <strong>0.2</strong></code>; <code>minimum = 0.1</code>; <code>maximum = 10.0</code>)<br/>
Lower cutoff in <i>pT</i>, below which no further interactions 
are allowed. Normally <i>pT0</i> above would be used to provide 
the main regularization of the cross section for <i>pT -> 0</i>, 
in which case <i>pTmin</i> is used  mainly for technical reasons. 
It is possible, however, to set <i>pT0Ref = 0</i> and use 
<i>pTmin</i> to provide a step-function regularization, or to 
combine them in intermediate approaches. Currently <i>pTmin</i> 
is taken to be energy-independent.  
  

<p/>
G&ouml;sta Gustafson has proposed (private communication, unpublished)
that the amount of screening, as encapsulated in the <i>pT0</i> 
parameter, fluctuates from one event to the next. Specifically, 
high-activity event are more likely to lead to interactions at large 
<i>pT</i> scales, but the high activity simultaneously leads to a 
larger screening of interactions at smaller <i>pT</i>. Such a scenario 
can approximately be simulated by scaling up the <i>pT0</i> by a 
factor <i>sqrt(n)</i>, where <i>n</i> is the number of interactions 
considered so far, including the current one. That is, for the first 
interaction the dampening factor is <i>pT^4 / (pT0^2 + pT^2)^2</i>, 
for the second <i>pT^4 / (2 pT0^2 + pT^2)^2</i>, for the third
<i>pT^4 / (3 pT0^2 + pT^2)^2</i>, and so on. Optionally the scheme
may also be applied to ISR emissions. For simplicity the same 
<i>alpha_s(pT0^2 + pT^2)</i> is used throughout. Note that, in this 
scenario the <i>pT0</i> scale must be lower than in the normal case 
to begin with, since it later is increased back up. Also note that the 
idea with this scenario is to propose an alternative to colour 
reconnection to understand the rise of <i>&lt;pT&gt;(n_ch)</i>, 
so that the amount of colour reconnection should be reduced.
<p/><code>mode&nbsp; </code><strong> MultipartonInteractions:enhanceScreening &nbsp;</strong> 
 (<code>default = <strong>0</strong></code>; <code>minimum = 0</code>; <code>maximum = 2</code>)<br/>
Choice to activate the above screening scenario, i.e. an increasing
effective <i>pT0</i> for consecutive interactions.
<br/><code>option </code><strong> 0</strong> : No activity-dependent screening, i.e. <i>pT0</i>
is fixed.  
<br/><code>option </code><strong> 1</strong> : The <i>pT0</i> scale is increased as a function
of the number of MPI's, as explained above. ISR is not affected,
but note that, if <code>SpaceShower:samePTasMPI</code> is on,
then <code>MultipartonInteractions:pT0Ref</code> is used also for ISR, 
which may or may not be desirable.
   
<br/><code>option </code><strong> 2</strong> : Both MPI and ISR influence and are influenced by the
screening. That is, the dampening is reduced based on the total number 
of MPI and ISR steps considered so far, including the current one.
This dampening is implemented both for MPI and for ISR emissions, 
for the latter provided that <code>SpaceShower:samePTasMPI</code> is on 
(default).  
  
  

<h4>Impact-parameter dependence</h4>
 
The choice of impact-parameter dependence is regulated by several
parameters. The ones listed here refer to nondiffractive topologies
only, while their equivalents for diffractive events are put in the 
<a href="Diffraction.html" target="page">Diffraction</a> description. Note that
there is currently no <code>bProfile = 4</code> option for diffraction.
Other parameters are assumed to agree between diffractive and 
nondiffractive topologies.

<p/><code>mode&nbsp; </code><strong> MultipartonInteractions:bProfile &nbsp;</strong> 
 (<code>default = <strong>1</strong></code>; <code>minimum = 0</code>; <code>maximum = 4</code>)<br/>
Choice of impact parameter profile for the incoming hadron beams.
<br/><code>option </code><strong> 0</strong> : no impact parameter dependence at all.  
<br/><code>option </code><strong> 1</strong> : a simple Gaussian matter distribution; 
no free parameters.  
<br/><code>option </code><strong> 2</strong> : a double Gaussian matter distribution, 
with the two free parameters <i>coreRadius</i> and 
<i>coreFraction</i>.  
<br/><code>option </code><strong> 3</strong> : an overlap function, i.e. the convolution of 
the matter distributions of the two incoming hadrons, of the form
<i>exp(- b^expPow)</i>, where <i>expPow</i> is a free 
parameter.   
<br/><code>option </code><strong> 4</strong> : a Gaussian matter distribution with a width
that varies according to the selected <i>x</i> value of an interaction,
<i>1. + a1 log (1 / x)</i>, where <i>a1</i> is a free parameter.
Note that once <i>b</i> has been selected for the hard process,
it remains fixed for the remainder of the evolution.
  
  

<p/><code>parm&nbsp; </code><strong> MultipartonInteractions:coreRadius &nbsp;</strong> 
 (<code>default = <strong>0.4</strong></code>; <code>minimum = 0.1</code>; <code>maximum = 1.</code>)<br/>
When assuming a double Gaussian matter profile, <i>bProfile = 2</i>,
the inner core is assumed to have a radius that is a factor
<i>coreRadius</i> smaller than the rest.
   

<p/><code>parm&nbsp; </code><strong> MultipartonInteractions:coreFraction &nbsp;</strong> 
 (<code>default = <strong>0.5</strong></code>; <code>minimum = 0.</code>; <code>maximum = 1.</code>)<br/>
When assuming a double Gaussian matter profile, <i>bProfile = 2</i>,
the inner core is assumed to have a fraction <i>coreFraction</i> 
of the matter content of the hadron.
   

<p/><code>parm&nbsp; </code><strong> MultipartonInteractions:expPow &nbsp;</strong> 
 (<code>default = <strong>1.</strong></code>; <code>minimum = 0.4</code>; <code>maximum = 10.</code>)<br/>
When <i>bProfile = 3</i> it gives the power of the assumed overlap 
shape <i>exp(- b^expPow)</i>. Default corresponds to a simple 
exponential drop, which is not too dissimilar from the overlap 
obtained with the standard double Gaussian parameters. For 
<i>expPow = 2</i> we reduce to the simple Gaussian, <i>bProfile = 1</i>, 
and for <i>expPow -> infinity</i> to no impact parameter dependence 
at all, <i>bProfile = 0</i>. For small <i>expPow</i> the program 
becomes slow and unstable, so the min limit must be respected.
   

<p/><code>parm&nbsp; </code><strong> MultipartonInteractions:a1 &nbsp;</strong> 
 (<code>default = <strong>0.15</strong></code>; <code>minimum = 0.</code>; <code>maximum = 2.</code>)<br/>
When <i>bProfile = 4</i>, this gives the <i>a1</i> constant in the
Gaussian width. When <i>a1 = 0.</i>, this reduces back to the single
Gaussian case.
  

<p/><code>mode&nbsp; </code><strong> MultipartonInteractions:bSelScale &nbsp;</strong> 
 (<code>default = <strong>1</strong></code>; <code>minimum = 1</code>; <code>maximum = 3</code>)<br/>
The selection of impact parameter is related to the scale of the
hard process: the harder this scale is, the more central the collision.
In practice this centrality saturates quickly, however, and beyond 
a scale of roughly 20 GeV very little changes. (The relevant 
quantity is that the QCD jet cross section above the scale should be
a tiny fraction of the total cross section.) In <i>2 -> 1</i> and
<i>2 -> 2</i> processes traditional scale choices work fine, but 
ambiguities arise for higher multiplicities, in particular when the
scale is used for matching between the multiparton matrix elements 
and parton showers. Then the event scale may 
be chosen as that of a very low-<i>pT</i> parton, i.e. suggesting a
peripheral collision, while the much harder other partons instead 
would favour a central collision. Therefore the default here is to  
override whatever scale value have been read in from an LHEF, say.
Notice that the scale used here is decoupled from the maximum scale 
for MPIs (<code>MultipartonInteractions:pTmaxMatch</code>).
<br/><code>option </code><strong> 1</strong> : 
Use the mass for a <i>2 -> 1</i> process. 
For <i>2 -> n, n > 1</i> processes order the particles in 
falling <i>mmT = m + mT</i> and then let the scale be 
<i>(mmT_1 + mmT_2)/2 + mmT_3/3 + mmT_4/4 + ... + mmT_n/n</i>.
This is constructed always to be above <i>m1</i>, and to assign 
decreasing importance to softer particles that are less likely to 
be associated with the hard process.  
<br/><code>option </code><strong> 2</strong> : Use the <code>scale</code> parameter of the event. 
  
<br/><code>option </code><strong> 3</strong> : use the same scale as chosen by the rules for
<code>MultipartonInteractions:pTmaxMatch</code>.  
  

<h4>Rescattering</h4>
 
It is possible that a parton may rescatter, i.e. undergo a further
interaction subsequent to the first one. The machinery to model this
kind of physics has only recently become fully operational
[<a href="Bibliography.html" target="page">Cor09</a>], and is therefore not yet so well explored. 

<p/>
The rescattering framework has ties with other parts of the program,
notably with the <a href="BeamRemnants.html" target="page">beam remnants</a>.

<p/><code>flag&nbsp; </code><strong> MultipartonInteractions:allowRescatter &nbsp;</strong> 
 (<code>default = <strong>off</strong></code>)<br/>
Switch to allow rescattering of partons; on/off = true/false.<br/>
<b>Note:</b> the rescattering framework has not yet been implemented
for the <code>MultipartonInteractions:bProfile = 4</code> option,
and can therefore not be switched on in that case.  
<b>Warning:</b> use with caution since machinery is still not 
so well tested.
  

<p/><code>flag&nbsp; </code><strong> MultipartonInteractions:allowDoubleRescatter &nbsp;</strong> 
 (<code>default = <strong>off</strong></code>)<br/>
Switch to allow rescattering of partons, where both incoming partons 
have already rescattered; on/off = true/false. Is only used if 
<code>MultipartonInteractions:allowRescatter</code> is switched on.<br/>
<b>Warning:</b> currently there is no complete implementation that 
combines it with shower evolution, so you must use 
<code>PartonLevel:ISR = off</code> and <code>PartonLevel:FSR = off</code>.
If not, a warning will be issued and double rescattering will not be 
simulated. The rate also comes out to be much lower than for single 
rescattering, so to first approximation it can be neglected.
  

<p/><code>mode&nbsp; </code><strong> MultipartonInteractions:rescatterMode &nbsp;</strong> 
 (<code>default = <strong>0</strong></code>; <code>minimum = 0</code>; <code>maximum = 4</code>)<br/>
Selection of which partons rescatter against unscattered partons
from the incoming beams A and B, based on their rapidity value 
<i>y</i> in the collision rest frame. Here <i>ySep</i> is
shorthand for <code>MultipartonInteractions:ySepRescatter</code> and
<i>deltaY</i> for <code>MultipartonInteractions:deltaYRescatter</code>,
defined below. The description is symmetric between the two beams,
so only one case is described below. 
<br/><code>option </code><strong> 0</strong> : only scattered partons with <i>y > 0</i> 
can collide with unscattered partons from beam B.   
<br/><code>option </code><strong> 1</strong> : only scattered partons with <i>y > ySep</i> 
can collide with unscattered partons from beam B.  
<br/><code>option </code><strong> 2</strong> : the probability for a scattered parton to be considered 
as a potential rescatterer against unscattered partons in beam B increases 
linearly from zero at <i>y = ySep - deltaY</i> to unity at 
<i>y = ySep + deltaY</i>.  
<br/><code>option </code><strong> 3</strong> : the probability for a scattered parton to be considered 
as a potential rescatterer against unscattered partons in beam B increases 
with <i>y</i> according to 
<i>(1/2) * (1 + tanh( (y - ySep) / deltaY))</i>.  
<br/><code>option </code><strong> 4</strong> : all partons are potential rescatterers against both 
beams.  
  

<p/><code>parm&nbsp; </code><strong> MultipartonInteractions:ySepRescatter &nbsp;</strong> 
 (<code>default = <strong>0.</strong></code>)<br/>
used for some of the <code>MultipartonInteractions:rescatterMode</code>
options above, as the rapidity for which a scattered parton has a 50%
probability to be considered as a potential rescatterer. 
A <i>ySep > 0</i> generally implies that some central partons cannot
rescatter at all, while a <i>ySep < 0</i>  instead allows central 
partons to scatter against either beam.
   

<p/><code>parm&nbsp; </code><strong> MultipartonInteractions:deltaYRescatter &nbsp;</strong> 
 (<code>default = <strong>1.</strong></code>; <code>minimum = 0.1</code>)<br/>
used for some of the <code>MultipartonInteractions:rescatterMode</code>
options above, as the width of the rapidity transition region, where the 
probability rises from zero to unity that a scattered parton is considered 
as a potential rescatterer.
   


<h3>Further variables</h3>

These should normally not be touched. Their only function is for
cross-checks.

<p/><code>mode&nbsp; </code><strong> MultipartonInteractions:nQuarkIn &nbsp;</strong> 
 (<code>default = <strong>5</strong></code>; <code>minimum = 0</code>; <code>maximum = 5</code>)<br/>
Number of allowed incoming quark flavours in the beams; a change 
to 4 would thus exclude <i>b</i> and <i>bbar</i> as incoming 
partons, etc.
  

<p/><code>mode&nbsp; </code><strong> MultipartonInteractions:nSample &nbsp;</strong> 
 (<code>default = <strong>1000</strong></code>; <code>minimum = 100</code>)<br/>
The allowed <i>pT</i> range is split (unevenly) into 100 bins, 
and in each of these the interaction cross section is evaluated in 
<i>nSample</i> random phase space points. The full integral is used 
at initialization, and the differential one during the run as a
"Sudakov form factor" for the choice of the hardest interaction.
A larger number implies increased accuracy of the calculations.
  

<h3>Technical notes</h3>

Relative to the articles mentioned above, not much has happened.
The main news is a technical one, that the phase space of the 
<i>2 -> 2</i> (massless) QCD processes is now sampled in 
<i>dy_3 dy_4 dpT^2</i>, where <i>y_3</i> and <i>y_4</i> are 
the rapidities of the two produced partons. One can show that
<br/><i>
    (dx_1 / x_1) * (dx_2 / x_2) * d(tHat) = dy_3 * dy_4 * dpT^2
</i><br/>
Furthermore, since cross sections are dominated by the "Rutherford"
one of <i>t</i>-channel gluon exchange, which is enhanced by a 
factor of 9/4 for each incoming gluon, effective structure functions 
are defined as
<br/><i>
    F(x, pT2) = (9/4) * xg(x, pT2) + sum_i xq_i(x, pT2) 
</i><br/>
With this technical shift of factors 9/4 from cross sections to parton 
densities, a common upper estimate of 
<br/><i>
    d(sigmaHat)/d(pT2) &lt; pi * alpha_strong^2 / pT^4   
</i><br/>
is obtained. 

<p/>
In fact this estimate can be reduced by a factor of 1/2 for the 
following reason: for any configuration <i>(y_3, y_4, pT2)</i> also 
one with <i>(y_4, y_3, pT2)</i> lies in the phase space. Not both 
of those can enjoy being enhanced by the <i>tHat -> 0</i> 
singularity of 
<br/><i>
    d(sigmaHat) propto 1/tHat^2. 
</i><br/>
Or if they are, which is possible with identical partons like 
<i>q q -> q q</i> and <i>g g -> g g</i>, each singularity comes 
with half the strength. So, when integrating/averaging over the two 
configurations, the estimated <i>d(sigmaHat)/d(pT2)</i> drops. 
Actually, it drops even further, since the naive estimate above is 
based on
<br/><i>
    (4 /9) * (1 + (uHat/sHat)^2) &lt; 8/9 &lt; 1
</i><br/>
The 8/9 value would be approached for <i>tHat -> 0</i>, which 
implies <i>sHat >> pT2</i> and thus a heavy parton-distribution 
penalty, while parton distributions are largest for 
<i>tHat = uHat = -sHat/2</i>, where the above expression 
evaluates to 5/9. A fudge factor is therefore introduced to go the 
final step, so it can easily be modified when further non-Rutherford 
processes are added, or should parton distributions change significantly.

<p/>
At initialization, it is assumed that  
<br/><i>
    d(sigma)/d(pT2) &lt; d(sigmaHat)/d(pT2) * F(x_T, pT2) * F(x_T, pT2)
       * (2 y_max(pT))^2
</i><br/>
where the first factor is the upper estimate as above, the second two
the parton density sum evaluated at <i>y_3 = y_ 4 = 0</i> so that 
<i>x_1 = x_2 = x_T = 2 pT / E_cm</i>, where the product is expected 
to be maximal, and the final is the phase space for
<i>-y_max &lt; y_{3,4} &lt; y_max</i>.
The right-hand side expression is scanned logarithmically in <i>y</i>, 
and a <i>N</i> is determined such that it always is below 
<i>N/pT^4</i>.

<p/>
To describe the dampening of the cross section at <i>pT -> 0</i> by
colour screening, the actual cross section is multiplied by a 
regularization factor <i>(pT^2 / (pT^2 + pT0^2))^2</i>, and the 
<i>alpha_s</i> is evaluated at a scale <i>pT^2 + pT0^2</i>, 
where <i>pT0</i> is a free parameter of the order of 2 - 4 GeV. 
Since <i>pT0</i> can be energy-dependent,  an ansatz
<br/><i>
    pT0(ecm) = pT0Ref * (ecm/ecmRef)^ecmPow
</i><br/>
is used, where <i>ecm</i> is the current CM frame energy, 
<i>ecmRef</i> is an arbitrary reference energy where <i>pT0Ref</i> 
is defined, and <i>ecmPow</i> gives the energy rescaling pace. For 
technical reasons, also an absolute lower <i>pT</i> scale <i>pTmin</i>, 
by default 0.2 GeV, is introduced. In principle, it is possible to 
recover older scenarios with a sharp <i>pT</i> cutoff by setting 
<i>pT0 = 0</i> and letting <i>pTmin</i> be a larger number. 

<p/>
The above scanning strategy is then slightly modified: instead of
an upper estimate <i>N/pT^4</i> one of the form 
<i>N/(pT^2 + r * pT0^2)^2</i> is used. At first glance, <i>r = 1</i> 
would seem to be fixed by the form of the regularization procedure, 
but this does not take into account the nontrivial dependence on 
<i>alpha_s</i>, parton distributions and phase space. A better 
Monte Carlo efficiency is obtained for <i>r</i> somewhat below unity, 
and currently <i>r = 0.25</i> is hardcoded.

In the generation a trial <i>pT2</i> is then selected according to
<br/><i>
    d(Prob)/d(pT2) = (1/sigma_ND) * N/(pT^2 + r * pT0^2)^2 * ("Sudakov")
</i><br/>
For the trial <i>pT2</i>, a <i>y_3</i> and a <i>y_4</i> are then 
selected, and incoming flavours according to the respective 
<i>F(x_i, pT2)</i>, and then the cross section is evaluated for this 
flavour combination. The ratio of trial/upper estimate gives the 
probability of survival.

<p/>
Actually, to profit from the factor 1/2 mentioned above, the cross
section for the combination with <i>y_3</i> and <i>y_4</i> 
interchanged is also tried, which corresponds to exchanging <i>tHat</i>
and <i>uHat</i>, and the average formed, while the final kinematics 
is given by the relative importance of the two.

<p/>
Furthermore, since large <i>y</i> values are disfavoured by dropping 
PDF's, a factor 
<br/><i>
   WT_y = (1 - (y_3/y_max)^2) * (1 - (y_4/y_max)^2) 
</i><br/>
is evaluated, and used as a survival probability before the more
time-consuming PDF+ME evaluation, with surviving events given a 
compensating weight <i>1/WT_y</i>. 

<p/>
An impact-parameter dependence is also allowed. Based on the hard 
<i>pT</i> scale of the first interaction, and enhancement/depletion 
factor is picked, which multiplies the rate of subsequent interactions.

<p/>
Parton densities are rescaled and modified to take into account the 
energy-momentum and flavours kicked out by already-considered 
interactions.

</body>
</html>

<!-- Copyright (C) 2013 Torbjorn Sjostrand -->