Past, Present and Future

SFT group meeting CERN 15 April 2005

PYTHIA:

Past, Present and Future

Torbj ¨ orn Sj ¨ ostrand

CERN/PH and

Department of Theoretical Physics, Lund University

The Physics of Event Generators The Event Generator Landscape

The move to C++

PYTHIA8 sneak preview

Who was Pythia?

The Oracle of Delphi: ca. 1000 B.C. — 390 A.D.

What is PYTHIA?

“real life”

Machine ⇒ events produce events

“virtual reality”

Event Generator

observe & store events

Detector,

Detector Simulation

what is

knowable? Event Reconstruction

compare real and

simulated data Physics Analysis

conclusions, articles, talks, . . .

“quick and dirty”

where and why?

• detector requirements

• analysis strategies

• acceptance corrections

physics is complex

What is PYTHIA?

“real life”

Machine ⇒ events LHC

produce events

“virtual reality”

Event Generator PYTHIA

observe & store events Detector,

ATLAS,CMS,LHC-B,ALICE

Detector Simulation Geant4, LCG what is

knowable? Event Reconstruction ORCA, ATHENA compare real and

simulated data Physics Analysis ROOT, JetClu

conclusions, articles, talks, . . .

“quick and dirty”

where and why?

• detector requirements

• analysis strategies

• acceptance corrections

physics is complex

A tour to Monte Carlo

. . . because Einstein was wrong: God does throw dice!

Quantum mechanics: amplitudes =⇒ probabilities

Anything that possibly can happen, will! (but more or less often)

The structure of an event

Warning: schematic only, everything simplified, nothing to scale, . . .

p

p/p

Incoming beams: parton densities

p

p/p

u g

W

d

Hard subprocess: described by matrix elements

p

p/p

u g

W

d

c s

Resonance decays: correlated with hard subprocess

p

p/p

u g

W

d

c s

Initial-state radiation: spacelike parton showers

p

p/p

u g

W

d

c s

Final-state radiation: timelike parton showers

p

p/p

u g

W

d

c s

Multiple parton–parton interactions . . .

p

p/p

u g

W

d

c s

. . . with its initial- and final-state radiation

Beam remnants and other outgoing partons

Everything is connected by colour confinement strings

Recall! Not to scale: strings are of hadronic widths

The strings fragment to produce primary hadrons

Many hadrons are unstable and decay further

Detector.gif (GIF Image, 460x434 pixels) http://atlas.web.cern.ch/Atlas/Detector.gif

1 of 1 02/06/2005 01:49 PM

These are the particles that hit the detector

The Monte Carlo method

Want to generate events in as much detail as Mother Nature

=⇒ get average and fluctutations right

=⇒ make random choices, ∼ as in nature

σ

= σ

P

(appropriately summed & integrated over non-distinguished final states) where P

= P

P

P

P

P

P

P

with P

=

P

=

P

= . . . in its turn

=⇒ divide and conquer

an event with n particles involves O (10n) random choices, (flavour, mass, momentum, spin, production vertex, lifetime, . . . ) LHC: ∼ 100 charged and ∼ 200 neutral (+ intermediate stages)

=⇒ several thousand choices

(of O (100) different kinds)

Generator Landscape

Hard Processes

Resonance Decays

Parton Showers Underlying Event

Hadronization

Ordinary Decays

General-Purpose

HERWIG

PYTHIA

ISAJET

SHERPA

Specialized a lot

HDECAY, . . .

Ariadne/LDC, NLLjet

DPMJET

none (?)

TAUOLA, EvtGen

specialized often best at given task, but need General-Purpose core

PYTHIA Process Library

No. Subprocess Hard QCD processes:

11 fifj → fifj

12 fifi→ fkfk

13 fifi→ gg 28 fig → fⁱg 53 gg → f^kfk

68 gg → gg Soft QCD processes:

91 elastic scattering 92 single diffraction (XB) 93 single diffraction (AX) 94 double diffraction 95 low-p⊥production Open heavy flavour:

(also fourth generation) 81 fifi→ QkQ_k 82 gg → Q^kQ_k 83 qifj → Qkfl

84 gγ → Q^kQ_k 85 γγ → F^kFk

Closed heavy flavour:

86 gg → J/ψg 87 gg → χ^0cg 88 gg → χ1cg 89 gg → χ^2cg 104 gg → χ^0c 105 gg → χ^2c 106 gg → J/ψγ 107 gγ → J/ψg 108 γγ → J/ψγ W/Z production:

1 fifi→ γ^∗/Z⁰ 2 fifj → W^± 22 fifi→ Z⁰Z⁰ 23 fifj → Z⁰W^± 25 fifi→ W⁺W⁻ 15 fifi→ gZ⁰ 16 fifj → gW^± 30 fig → fⁱZ⁰ 31 fig → f^kW^± 19 fifi→ γZ⁰ 20 fifj → γW^± 35 fiγ → fⁱZ⁰

No. Subprocess 36 fiγ → f^kW^± 69 γγ → W⁺W⁻ 70 γW^±→ Z⁰W^± Prompt photons:

14 fifi→ gγ 18 fifi→ γγ 29 fig → fⁱγ 114 gg → γγ 115 gg → gγ Deeply Inel. Scatt.:

10 fifj → fkfl

99 γ^∗q → q Photon-induced:

33 fiγ → fⁱg 34 fiγ → fⁱγ 54 gγ → f^kfk

58 γγ → f^kfk

131 fiγ_T^∗→ fⁱg 132 fiγ_L^∗→ fig 133 fiγT^∗→ fiγ 134 fiγL^∗→ fiγ 135 gγT^∗→ fifi

136 gγL^∗→ fifi

137 γT^∗γ^∗T→ fifi

138 γT^∗γ^∗L→ fifi

139 γL^∗γ^∗T→ fifi

140 γL^∗γ^∗L→ fifi

80 qiγ → qkπ^± Light SM Higgs:

3 fifi→ h⁰ 24 fifi→ Z⁰h⁰ 26 fifj → W^±h⁰ 32 fig → fⁱh⁰ 102 gg → h⁰ 103 γγ → h⁰ 110 fifi→ γh⁰ 111 fifi→ gh⁰ 112 fig → fⁱh⁰ 113 gg → gh⁰ 121 gg → Q^kQ_kh⁰ 122 qiq_i→ QkQ_kh⁰ 123 fifj → fifjh⁰ 124 fifj → fkflh⁰

No. Subprocess New gauge bosons:

141 fifi→ γ/Z⁰/Z^{0 0} 142 fifj → W^{0 +} 144 fifj → R Heavy SM Higgs:

5 Z⁰Z⁰→ h⁰ 8 W⁺W⁻→ h⁰ 71 Z⁰_LZ⁰_L→ Z⁰LZ⁰_L 72 Z⁰_LZ⁰_L→ W⁺LW⁻_L 73 Z⁰LW^±_L → Z⁰LW_L^± 76 W⁺_LW_L⁻→ Z⁰LZ⁰L

77 W^±_LW^±_L → WL^±W^±_L BSM Neutral Higgs:

151 fifi→ H⁰ 152 gg → H⁰ 153 γγ → H⁰ 171 fifi→ Z⁰H⁰ 172 fifj → W^±H⁰ 173 fifj → fifjH⁰ 174 fifj → fkflH⁰ 181 gg → Q^kQ_kH⁰ 182 qiq_i→ QkQ_kH⁰ 183 fifi→ gH⁰ 184 fig → fⁱH⁰ 185 gg → gH⁰ 156 fifi→ A⁰ 157 gg → A⁰ 158 γγ → A⁰ 176 fifi→ Z⁰A⁰ 177 fifj → W^±A⁰ 178 fifj → fifjA⁰ 179 fifj → fkflA⁰ 186 gg → Q^kQ_kA⁰ 187 qiq_i→ QkQ_kA⁰ 188 fifi→ gA⁰ 189 fig → fⁱA⁰ 190 gg → gA⁰ Charged Higgs:

143 fifj → H⁺ 161 fig → f^kH⁺ 401 gg → tbH⁺ 402 qq → tbH⁺

No. Subprocess Higgs pairs:

297 fifj→ H^±h⁰ 298 fifj→ H^±H⁰ 299 fifi→ A⁰h⁰ 300 fifi→ A⁰H⁰ 301 fifi→ H⁺H⁻ Leptoquarks:

145 qi`j → LQ

162 qg → `L^Q 163 gg → L^QLQ

164 qiq_i→ LQLQ

Technicolor:

149 gg → η^tc 191 fifi→ ρ⁰tc

192 fifj→ ρ⁺tc

193 fifi→ ωtc⁰

194 fifi→ fkfk

195 fifj→ fkfl

361 fifi→ W⁺LW⁻_L 362 fifi→ W^±Lπ_tc^∓ 363 fifi→ π⁺tcπ_tc⁻ 364 fifi→ γπtc⁰

365 fifi→ γπ⁰⁰tc

366 fifi→ Z⁰πtc⁰

367 fifi→ Z⁰π⁰⁰tc

368 fifi→ W^±π_tc^∓ 370 fifj→ W^±LZ⁰L

371 fifj→ W^±Lπtc⁰

372 fifj→ π^±tcZ⁰L

373 fifj→ π^±tcπtc⁰

374 fifj→ γπtc^±

375 fifj→ Z⁰π_tc^± 376 fifj→ W^±πtc⁰

377 fifj→ W^±π⁰⁰tc

381 qiqj→ qiqj

382 qiq_i→ qkq_k 383 qiq_i→ gg 384 fig → fⁱg 385 gg → q^kqk

386 gg → gg 387 fifi→ QkQ_k 388 gg → Q^kQ_k

No. Subprocess Compositeness:

146 eγ → e^∗ 147 dg → d^∗ 148 ug → u^∗ 167 qiqj → d^∗qk

168 qiqj → u^∗qk

169 qiq_i→ e^±e^∗∓

165 fifi(→ γ^∗/Z⁰) → f^kfk

166 fifj(→ W^±) → f^kfl

Extra Dimensions:

391 ff → G^∗ 392 gg → G^∗ 393 qq → gG^∗ 394 qg → qG^∗ 395 gg → gG^∗ Left–right symmetry:

341 `i`j → H^±±L

342 `i`j → H^±±_R 343 `<sup>±</sup><sub>i</sub>γ → H<sup>±±</sup>L e<sup>∓</sup> 344 `^±_iγ → H^±±R e^∓ 345 `<sup>±</sup><sub>i</sub>γ → H<sup>±±</sup>L µ<sup>∓</sup> 346 `^±_iγ → H^±±R µ^∓ 347 `<sup>±</sup><sub>i</sub>γ → H<sup>±±</sup>L τ<sup>∓</sup> 348 `^±_iγ → H^±±R τ^∓ 349 fifi→ H⁺⁺L H⁻⁻_L 350 fifi→ H⁺⁺_R H⁻⁻_R 351 fifj → fkflH^±±_L 352 fifj → fkflH^±±_R 353 fifi→ Z⁰R

354 fifj → W^±R

SUSY:

201 fifi→ ˜eL˜e^∗_L 202 fifi→ ˜eR˜e^∗_R 203 fifi→ ˜eL˜e^∗_R+ 204 fifi→ ˜µLµ˜^∗_L 205 fifi→ ˜µRµ˜^∗_R 206 fifi→ ˜µLµ˜^∗_R+ 207 fifi→ ˜τ1τ˜1^∗

208 fifi→ ˜τ2τ˜₂^∗ 209 fifi→ ˜τ1τ˜₂^∗+

No. Subprocess 210 fifj → ˜`Lν˜`^∗+ 211 fifj → ˜τ1ν˜τ^∗+ 212 fifj → ˜τ2ν˜τ^∗+ 213 fifi→ ˜ν`ν˜`∗

214 fifi→ ˜ντν˜^∗τ

216 fifi→ ˜χ1χ˜1

217 fifi→ ˜χ2χ˜2

218 fifi→ ˜χ3χ˜3

219 fifi→ ˜χ4χ˜4

220 fifi→ ˜χ1χ˜2

221 fifi→ ˜χ1χ˜3

222 fifi→ ˜χ1χ˜4

223 fifi→ ˜χ2χ˜3

224 fifi→ ˜χ2χ˜4

225 fifi→ ˜χ3χ˜4

226 fifi→ ˜χ^±₁χ˜^∓₁ 227 fifi→ ˜χ^±₂χ˜^∓₂ 228 fifi→ ˜χ^±₁χ˜^∓₂ 229 fifj → ˜χ1χ˜^±₁ 230 fifj → ˜χ2χ˜^±₁ 231 fifj → ˜χ3χ˜^±₁ 232 fifj → ˜χ4χ˜^±₁ 233 fifj → ˜χ1χ˜^±₂ 234 fifj → ˜χ2χ˜^±₂ 235 fifj → ˜χ3χ˜^±₂ 236 fifj → ˜χ4χ˜^±₂ 237 fifi→ ˜g ˜χ1

238 fifi→ ˜g ˜χ2

239 fifi→ ˜g ˜χ3

240 fifi→ ˜g ˜χ4

241 fifj → ˜g ˜χ^±₁ 242 fifj → ˜g ˜χ^±₂ 243 fifi→ ˜g˜g 244 gg → ˜g˜g 246 fig → ˜q^iLχ˜1

247 fig → ˜q^iRχ˜1

248 fig → ˜q^iLχ˜2

249 fig → ˜q^iRχ˜2

No. Subprocess 250 fig → ˜q^iLχ˜3

251 fig → ˜q^iRχ˜3

252 fig → ˜q^iLχ˜4

253 fig → ˜q^iRχ˜4

254 fig → ˜q^{j L}χ˜^±₁ 256 fig → ˜q^{j L}χ˜^±₂ 258 fig → ˜q^iL˜g 259 fig → ˜qiR˜g 261 fifi→ ˜t1˜t^∗1

262 fifi→ ˜t2˜t^∗2

263 fifi→ ˜t1˜t^∗2+ 264 gg → ˜t¹˜t^∗1

265 gg → ˜t²˜t^∗2

271 fifj → ˜qiLq˜j L

272 fifj → ˜qiR˜qj R

273 fifj → ˜qiLq˜j R+ 274 fifj → ˜qiLq˜^∗j L

275 fifj → ˜qiR˜q^∗j R

276 fifj → ˜qiLq˜^∗_{j R}+ 277 fifi→ ˜qj Lq˜^∗_{j L} 278 fifi→ ˜qj R˜q^∗_{j R} 279 gg → ˜q^iLq˜^∗i L

280 gg → ˜q^iR˜q^∗i R

281 bqi→ ˜b1˜qiL

282 bqi→ ˜b2˜qiR

283 bqi→ ˜b1˜qiR+ 284 bq_i→ ˜b1˜q^∗i L

285 bq_i→ ˜b2˜q^∗i R

286 bq_i→ ˜b1˜q^∗_{i R}+ 287 fifi→ ˜b1˜b^∗1

288 fifi→ ˜b2˜b^∗2

289 gg → ˜b¹˜b^∗1

290 gg → ˜b²˜b^∗2

291 bb → ˜b1b˜1

292 bb → ˜b²b˜2

293 bb → ˜b¹b˜2

294 bg → ˜b¹˜g 295 bg → ˜b²˜g 296 bb → ˜b¹b˜^∗2+

The Les Houches Accord

Specialized Generator

=⇒ Hard Process

Les Houches Interface

HERWIG or PYTHIA (Resonance Decays) Parton Showers

Underlying Event Hadronization Ordinary Decays

Some Specialized Generators:

• AcerMC: ttbb, . . .

• ALPGEN: W/Z+ ≤ 6j,

nW + mZ + kH+ ≤ 3j, . . .

• AMEGIC++: generic LO

• CompHEP: generic LO

• GRACE+Bases/Spring:

generic LO+ some NLO loops

• GR@PPA: bbbb

• MadCUP: W/Z+ ≤ 3j, ttbb

• MadGraph+HELAS: generic LO

• MCFM: NLO W/Z+ ≤ 2j, WZ, WH, H+ ≤ 1j

• O’Mega+WHIZARD: generic LO

• VECBOS: W/Z+ ≤ 4j

Apologies for all unlisted programs

The Bigger Picture

Process Selection Resonance Decays

Parton Showers Multiple Interactions

Beam Remnants

Hadronization Ordinary Decays

Detector Simulation ME Generator

ME Expression

SUSY/. . . spectrum calculation

Phase Space Generation

PDF Library

τ Decays

B Decays

=⇒ need standardized interfaces (LHAPDF, SUSY LHA, . . . )

PYTHIA history

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

   

  

  

  

  

  

  

  

  

  

    

ll

   

   

   

   

   

   

   

   

   

   



 

 

 

 

 

 

 

 

    

    

    

    

    

    

PYTHIA

lh hh

1978 JETSET versions 1–7 string frag.

e

e

, FSR

1982 PYTHIA

versions 1–5 pp, ISR, MI 1997

PYTHIA 6.1

2001

PYTHIA 6.2 2003

PYTHIA 6.3

All in Fortran 77

2005

Time axis not to scale

PYTHIA standalone, but other programs rely on PYTHIA:

LEPTO

ARIADNE/LDC

RAPGAP/CASCADE POMPYT

HIJING

SHERPA

EVTGEN

. . .

On To C++

Currently HERWIG and PYTHIA are successfully being used, also in new LHC environments, using C++ wrappers

Q: Why rewrite?

A1: Need to clean up!

A2: Fortran 77 is limiting Q: Why C++?

A1: All the reasons for ROOT, Geant4, . . . (“a better language”, industrial standard, . . . )

A2: Young experimentalists will expect C++

(educational and professional continuity) A3: Only game in town! Fortran 90

So far mixed experience:

• Conversion effort: everything takes longer and costs more (as for LHC machine, detectors and software)

• The physics hurdle is as steep as the C++ learning curve

C++ Players

PYTHIA7 project =⇒ ThePEG

Toolkit for High Energy Physics Event Generation (L. L ¨onnblad; S. Gieseke, A. Ribon, P. Richardson)

HERWIG++: complete reimplementation

(B.R. Webber; S. Gieseke, A. Ribon, P. Richardson, M. Seymour, P. Stephens, 3 new)

ARIADNE/LDC: to do ISR/FSR showers, multiple interactions (L. L ¨onnblad; N. Lavesson)

SHERPA: in C++ from start, partly wrappers to PYTHIA Fortran (F. Krauss; T. Gleisberg, S. Hoeche, A. Schaelicke,

S. Schumann, J. Winter)

PYTHIA8: restart to write complete event generator

(T. Sj ¨ostrand, (S. Mrenna?, P. Skands?))

What is ThePEG?

Toolkit for High Energy Physics Event Generation CLHEP

utilities

ThePEG

basic structure

HERWIG++

physics modules

PYTHIA7

physics modules

Ariadne/LDC

physics modules

?

PYTHIA8 · · ·

not SHERPA

The new generator Herwig++

A completely new event generator in C++

• Aiming at full multi–purpose generator for LHC and future colliders.

• Preserving main features of HERWIG such as – angular ordered parton shower

– cluster hadronization

• New features and improvements – covariant shower formulation

– improved parton shower evolution for heavy quarks

– consistent radiation from unstable particles (multiscale evolution)

Growth of Fortran HERWIG

Bryan Webber, QCD Simulation for LHC and Herwig++, KEK, 6 April 2004 17

What’s next?

Near Future. . .

H Initial state shower:

• Complete implementation and tests.

H Refine e⁺e⁻:

• Full CKKW ME+PS matching.

• Precision tune to LEP data should be possible.

H with IS and FS showers running:

• we can start to test Drell–Yan and jets in pp collisions.

• cross check with Tevatron data and finally make predictions for the LHC.

H Underlying Event.

H Hadronic Decays: NEW! many new decayers, τ –decays, Spin correlations (P Richardson).

H New Ideas: soft gluons, improved shower algorithm, NLO, . . .

Schedule?

• Ready for LHC!

Mike Seymour, Moriond 2005 33

SHERPA

MC for LHC 4 Mike Seymour

PYTHIA8: A fresh start

Problem: Leif not committed to PYTHIA, no other manpower Solution?: take a sabbatical and work “full-time”!

Tentative schedule:

time date processes final states

0 = 1 Sept. 2004 — —

1 = 1 Sept. 2005 LHA-style input incomplete draft 2 = 1 Sept. 2006 a few processes complete, buggy(?) 3 = 1 Sept. 2007 more processes stable, debugged

. . . but don’t forget Murphy’s law Objectives:

• concentrate on physics, not administration

• clean up, keep the most recent models

• pure standard C++, no fancy programming tricks

• independent of ThePEG (or anything else), but

• Les Houches Accord style input central

• interface to ThePEG later written by Leif (?)