- Page 1 and 2: hep-ph/0603175LU TP 06-13FERMILAB-P
- Page 4 and 5: Contents1 Introduction 11.1 The Com
- Page 6 and 7: 8.6.3 Left-Right Symmetry and Doubl
- Page 8 and 9: 13 Particles and Their Decays 38413
- Page 10 and 11: another sense, the overall energy f
- Page 13 and 14: p ⊥ -ordered parton showers and a
- Page 15 and 16: you expect to happen if you do not
- Page 17 and 18: The history of the Pythia program i
- Page 19 and 20: andom, according to the relative co
- Page 21 and 22: Of course, the above ‘resonance
- Page 23 and 24: 2.2.1 Matrix elementsMatrix element
- Page 25 and 26: analysis strongly suggests the scal
- Page 27 and 28: y calling PYEVNW instead of PYEVNT,
- Page 29 and 30: charge Casimir operators, N C /C F
- Page 31 and 32: 3 Program OverviewThis section cont
- Page 33 and 34: Table 2: The main versions of Pythi
- Page 35 and 36: • A new master event-generation r
- Page 37 and 38: A test program, PYTEST, is included
- Page 39 and 40: argument is used with a value that
- Page 41 and 42: &XMI(2,240),Q2MI(240),IMISEP(0:240)
- Page 43 and 44: original u and u. The parentheses e
- Page 45: C...End of particle and event loops
- Page 49 and 50: third the particle-flavour code (wh
- Page 51 and 52: 4 Monte Carlo TechniquesQuantum mec
- Page 53 and 54: different g i ’s, it is possible
- Page 55 and 56: The last equality is most easily se
- Page 57 and 58: Purpose: to generate a (pseudo)rand
- Page 59 and 60: 5 The Event RecordThe event record
- Page 61 and 62: Table 4: Gauge boson and other fund
- Page 63 and 64: Table 7: Meson codes, part 1.KF Nam
- Page 65: Table 10: QCD effective states.KF P
- Page 68 and 69: is then required. Normally this mea
- Page 70 and 71: quarks. Note that the positions nee
- Page 72 and 73: original entries appear with pointe
- Page 74 and 75: consistent summary, but then in div
- Page 76 and 77: = 3 : a documentation line, defined
- Page 78 and 79: 6 The Old Electron-Positron Annihil
- Page 80 and 81: 6.1.2 First-order QCD matrix elemen
- Page 82 and 83: 6.1.4 Second-order three-jet matrix
- Page 84 and 85: evaluated and compared with a param
- Page 86 and 87: The angular orientation of a 3- or
- Page 88 and 89: these offer unique possibilities to
- Page 90 and 91: 6.3.3 Common-block variablesThe sta
- Page 92 and 93: 3-jet cross section, so large that
- Page 94 and 95: fraction of events containing a rad
- Page 96 and 97:
! PYSPHE to work on it (unusual)CAL
- Page 98 and 99:
Pdflib has been frozen in recent ye
- Page 100 and 101:
modelling of γp and γγ events, s
- Page 102 and 103:
One can obtain the positron distrib
- Page 104 and 105:
7.2 Kinematics and Cross Section fo
- Page 106 and 107:
7.3 Resonance ProductionThe simples
- Page 108 and 109:
y a factor 1/N C = 1/3. Secondly, t
- Page 110 and 111:
This choice certainly is not unique
- Page 112 and 113:
The phase-space points tested at in
- Page 114 and 115:
spin-1 resonance the expression is,
- Page 116 and 117:
probabilities.In some processes, su
- Page 118 and 119:
where decay channels are chosen acc
- Page 120 and 121:
and for a resonance-antiresonance p
- Page 122 and 123:
The Regge formulae above for single
- Page 124 and 125:
emnant, but this remnant has a larg
- Page 126 and 127:
called direct, our direct×VMD and
- Page 128 and 129:
are also covered in other places in
- Page 130 and 131:
Table 17: Subprocess codes, part 2.
- Page 132 and 133:
Table 19: Subprocess codes, part 4.
- Page 134 and 135:
Table 21: Subprocess codes, part 6.
- Page 136 and 137:
Table 23: Subprocess codes, part 8.
- Page 138 and 139:
Table 25: Subprocess codes, part 10
- Page 140 and 141:
8.2 QCD ProcessesObviously most pro
- Page 142 and 143:
1. if MSEL = 4 - 8 then the flavour
- Page 144 and 145:
into the cross section. However, cu
- Page 146 and 147:
Uncertainties come from a number of
- Page 148 and 149:
If you are only concerned with stan
- Page 150 and 151:
for MSTP(14) = 1 and the VMD part o
- Page 152 and 153:
vertex. MSTP(66) offer some alterna
- Page 154 and 155:
calculation is preferred. Some caut
- Page 156 and 157:
contains a full two-Higgs-multiplet
- Page 158 and 159:
8.5.2 Heavy Standard Model HiggsISU
- Page 160 and 161:
parton-distribution-function approa
- Page 162 and 163:
The advantage of the explicit pair
- Page 164 and 165:
handed neutrinos. The Higgs fields
- Page 166 and 167:
former and ŝ for the latter. To ma
- Page 168 and 169:
mode is therefore tt, if kinematica
- Page 170 and 171:
through a vector-dominance mechanis
- Page 172 and 173:
where α a = ga/4π.2The phenomenol
- Page 174 and 175:
has two scalar partners, one associ
- Page 176 and 177:
above. To obtain proper linking wit
- Page 178 and 179:
or (ii) the masses of the first- an
- Page 180 and 181:
mSUGRA model input parameters shoul
- Page 182 and 183:
stable. If this is not the case, MW
- Page 184 and 185:
is no a priori generic prediction f
- Page 186 and 187:
tools exist. Therefore the producti
- Page 188 and 189:
’gamma/lepton’ machinery. There
- Page 190 and 191:
numerical precision may suffer; if
- Page 192 and 193:
call, e.g. at the end of the run, o
- Page 194 and 195:
on so as to give the full (paramete
- Page 196 and 197:
MSEL = 10) (ISUB = 19, 20, 22, 23,
- Page 198 and 199:
CKIN(7), CKIN(8) : (D = −10., 10.
- Page 200 and 201:
For two resonances that are identic
- Page 202 and 203:
= 1 : you should set couplings for
- Page 204 and 205:
or GVMD; we will use both interchan
- Page 206 and 207:
MSTP(17) : (D = 4) possibility of a
- Page 208 and 209:
h 0 and H 0 . It is mainly intended
- Page 210 and 211:
correct only in the Standard Model
- Page 212 and 213:
for three technicolors, based on sc
- Page 214 and 215:
esponsibility on you.Note 2: when P
- Page 216 and 217:
= 1 : the recommended standard for
- Page 218 and 219:
where a reconnection is actually ma
- Page 220 and 221:
MSTP(145) : (D = 0) choice of polar
- Page 222 and 223:
PARP(32) : (D = 2.0) special K fact
- Page 224 and 225:
PARP(121) : (D = 1.) the maxima obt
- Page 226 and 227:
Q2 : the momentum transfer scale Q
- Page 228 and 229:
1-PARU(136) corresponding to the sa
- Page 230 and 231:
cot θ 3 .ITCM(5) : (D = 0) presenc
- Page 232 and 233:
are SUGRA, SSMSSM and VISAJE, locat
- Page 234 and 235:
IMSS(22) : (D = 0) Read-in of SLHA
- Page 236 and 237:
it is the constant of proportionali
- Page 238 and 239:
k = 1 : only possibility.ISUB = 13,
- Page 240 and 241:
number is obtained as a by-product
- Page 242 and 243:
vanishing if no splitting is needed
- Page 244 and 245:
CALL PYINIT(’CMS’,’p’,’pb
- Page 246 and 247:
C...Histograms.CALL PYBOOK(1,’dn_
- Page 248 and 249:
Of course, the separation of which
- Page 250 and 251:
does imply an infinite total cross
- Page 252 and 253:
Pythia. The solution adopted here
- Page 254 and 255:
est strategy would vary, depending
- Page 256 and 257:
tered during the course of the run,
- Page 258 and 259:
NPRUP entries of the XSECUP, XERRUP
- Page 260 and 261:
should be denoted by using the valu
- Page 262 and 263:
more information than is provided b
- Page 264 and 265:
9.9.3 An exampleTo exemplify the ab
- Page 266 and 267:
DOUBLE PRECISION EBMUP,XSECUP,XERRU
- Page 268 and 269:
momenta and masses can be reconstru
- Page 270 and 271:
quantity, the effect of such events
- Page 272 and 273:
9.10 Interfaces to Other Generators
- Page 274 and 275:
= 2 : assign the interference contr
- Page 276 and 277:
value before the routine is called.
- Page 278 and 279:
• PYSGQC : normal QCD processes,
- Page 280 and 281:
afterwards.)MINT(17), MINT(18) : fl
- Page 282 and 283:
= 0 : normal initialization.= 1 : i
- Page 284 and 285:
agrees with the Bjorken definition,
- Page 286 and 287:
action; used to find what is left f
- Page 288 and 289:
VINT(319) : photon flux factor in P
- Page 290 and 291:
COMMON/PYINT4/MWID(500),WIDS(500,5)
- Page 292 and 293:
Purpose: to store information on to
- Page 294 and 295:
10 Initial- and Final-State Radiati
- Page 296 and 297:
for the original parton a. On the o
- Page 298 and 299:
180 ◦ and therefore the emission
- Page 300 and 301:
anching activity of the gluon. (Wit
- Page 302 and 303:
anching of both the q and the q, in
- Page 304 and 305:
10.2.5 Other final-state shower asp
- Page 306 and 307:
wide selection of generic colour an
- Page 308 and 309:
e a function of the transverse mome
- Page 310 and 311:
That way we hope to achieve the mos
- Page 312 and 313:
step by step one moves ‘backwards
- Page 314 and 315:
in hadronic interactions, with the
- Page 316 and 317:
happens it is almost always for one
- Page 318 and 319:
need be.The scheme presented above
- Page 320 and 321:
with the new PYPTIS one. The advant
- Page 322 and 323:
of the system.)CALL PYPTFS(NPART,IP
- Page 324 and 325:
Note:a few non-standard options hav
- Page 326 and 327:
g → gg and reduced angle for g
- Page 328 and 329:
(MSTP(42) = 2) and MSTJ(46) = 3).PA
- Page 330 and 331:
= 3 : the k ⊥ of the anomalous/GV
- Page 332 and 333:
to the cut implied by PARP(65).PARP
- Page 334 and 335:
outlines where the intermediate and
- Page 336 and 337:
11.1.4 Primordial k ⊥It is custom
- Page 338 and 339:
mended first to read this section,
- Page 340 and 341:
In general, there is quite a strong
- Page 342 and 343:
• Scatterings of the gg → gg ty
- Page 344 and 345:
where the denominator comes from th
- Page 346 and 347:
hard processes in pp or pp collisio
- Page 348 and 349:
proton net flavour content, and ene
- Page 350 and 351:
can be used here, from sharing the
- Page 352 and 353:
There is also the matter of a lower
- Page 354 and 355:
The program needs to know the assum
- Page 356 and 357:
Λ FSR scale: PARJ(81) = 0.14D0;Bea
- Page 358 and 359:
‘correct’ total cross section.M
- Page 360 and 361:
generated with the new model (in PY
- Page 362 and 363:
e attached to the remnant, PARP(80)
- Page 364 and 365:
Determine physical PT2MAX and PT2MI
- Page 366 and 367:
COMMON/PYCTAG/NCT, MCT(4000,2)Purpo
- Page 368 and 369:
string is assumed to have no transv
- Page 370 and 371:
Only two baryon multiplets are incl
- Page 372 and 373:
Additional parameters include the r
- Page 374 and 375:
or the q one, ‘left-right symmetr
- Page 376 and 377:
are changed, some retuning should b
- Page 378 and 379:
end in the initial stages of the st
- Page 380 and 381:
The fact that the hadron should be
- Page 382 and 383:
(including the one between the two
- Page 384 and 385:
the ratio of rescaled jet momentum
- Page 386 and 387:
duly randomized, but properly it wo
- Page 388 and 389:
Other models include a simplified i
- Page 390 and 391:
The negative sign is exactly what w
- Page 392 and 393:
13 Particles and Their DecaysPartic
- Page 394 and 395:
multiplet m 0 kpseudoscalars and ve
- Page 396 and 397:
toIn Dalitz decays, π 0 or η →
- Page 398 and 399:
fragmentation description. This doe
- Page 400 and 401:
14 The Fragmentation and Decay Prog
- Page 402 and 403:
Initially the partons must be given
- Page 404 and 405:
SUBROUTINE PYZDIS(KFL1,KFL3,PR,Z) :
- Page 406 and 407:
MSTU(4) : (D = 4000) number of line
- Page 408 and 409:
full. Is reset at each PYEXEC call.
- Page 410 and 411:
= 4 : transverse momenta are compen
- Page 412 and 413:
included in the event. If the showe
- Page 414 and 415:
allowed to shower, is 0 if no pair
- Page 416 and 417:
an independent gluon jet generated
- Page 418 and 419:
• Increase PARJ(1) by about a fac
- Page 420 and 421:
With five possible flavours for q 1
- Page 422 and 423:
LFN :is intended for the program au
- Page 424 and 425:
where both KTAB1 and KTAB3 are know
- Page 426 and 427:
a fourth generation is small.Remark
- Page 428 and 429:
code it is possible to define chann
- Page 430 and 431:
to be borne out by comparisons with
- Page 432 and 433:
CALL PYEXEC! particles decay, excep
- Page 434 and 435:
15 Event Study and Analysis Routine
- Page 436 and 437:
it can be further specified to +z a
- Page 438 and 439:
= 8 : KF code (K(I,2)) for a remain
- Page 440 and 441:
15.2 Event ShapesIn this section we
- Page 442 and 443:
To find the thrust value and axis t
- Page 444 and 445:
Hence rapidity, azimuthal angle and
- Page 446 and 447:
as a starting point for the iterati
- Page 448 and 449:
may be performed in E rather than i
- Page 450 and 451:
Here n = ∑ j n j is the total mul
- Page 452 and 453:
MSTU(63) and PARU(61) - PARU(63). I
- Page 454 and 455:
H40 : H 4 /H 0 .Remark: the number
- Page 456 and 457:
V(I,2) - V(I,5) = statistical error
- Page 458 and 459:
PYTHRU, PYCLUS and PYCELL.= 1 : sto
- Page 460 and 461:
MSTU(53)).PARU(57) : (D = 3.2) maxi
- Page 462 and 463:
Purpose: to dump the contents of ex
- Page 464 and 465:
References[Abb87]A. Abbasabadi and
- Page 466 and 467:
[Bai83] R. Baier and R. Rückl, Z.
- Page 468 and 469:
[Bra64] S. Brandt, Ch. Peyrou, R. S
- Page 470 and 471:
[Djo97] A. Djouadi, J. Kalinowski a
- Page 472 and 473:
[Gie04]S. Gieseke, A. Ribon, M.H. S
- Page 474 and 475:
[Ing80] G. Ingelman and T. Sjöstra
- Page 476 and 477:
[Lan00] K. Lane, T. Rador and E. Ei
- Page 478 and 479:
[OPA91] OPAL Collaboration, M.Z. Ak
- Page 480 and 481:
[Sjö92][Sjö92a][Sjö92b]T. Sjöst
- Page 482 and 483:
Appendix A: Subprocess Summary Tabl
- Page 484 and 485:
No. Subprocess No. SubprocessNo. Su
- Page 486 and 487:
Appendix B: Index of Subprograms an
- Page 488 and 489:
PYK function 429PYKCUT subroutine 1