The Monte Carlo event generator DPMJET-III 8th International - - PowerPoint PPT Presentation

the monte carlo event generator dpmjet iii
SMART_READER_LITE
LIVE PREVIEW

The Monte Carlo event generator DPMJET-III 8th International - - PowerPoint PPT Presentation

The Monte Carlo event generator DPMJET-III 8th International Workshop on Multiple Partonic Interactions at the LHC Anatoli Fedynitch (DESY Zeuthen) with Ralph Engel (KIT) DPMJET-III Tevatron SPS LHC years, max. E per nucleon DPMJET II.5 EOL!


slide-1
SLIDE 1

8th International Workshop on Multiple Partonic Interactions at the LHC

The Monte Carlo event generator DPMJET-III

Anatoli Fedynitch (DESY Zeuthen) with Ralph Engel (KIT)

slide-2
SLIDE 2

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 2

DPMJET-III

EOL! DTUJET DTUNUC DPMJET II DPMJET II.5 DTUNUC 2 PHOJET DPMJET III (3.0.6) h-h, e/ɣ-h, e/ɣ-e/ɣ h/A-A/A, ɣ-h/A, ɣ-ɣ; h-A, A-A @ cosmic ray energies p-p collisions,

  • nly gluon jets

h-A, A-A, no minijets

GVDM

years, max. E per nucleon SPS Tevatron LHC DPMJET III 2016+ (no separate PHOJET distribution)

Hahn & Ranft (1990) Möhring & Ranft (1990) Rösler & Ranft (1996) Ranft (1997) Engel & Ranft (1995) Rösler, Engel & Ranft(2000) AF & Engel (2016)

slide-3
SLIDE 3

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 3

What’s new in DPMJET 2016 > DPMJET III at cosmic ray energies

§ Extrapolation up to PeV (center-of-mass) § Simultaneous handling of multiple projectile-target combinations (cascade codes)

> Minimum bias

§ New PDF (GRV94 -> CT14 LO) § Re-tune § More accurate MPI behavior § Cross section fits

> Untouched:

§ Nuclear phenomenology § Heavy-Ions

> Distributed in

§ FLUKA 2017 § CORSIKA 2017 § CRMC, standalone in progress § Will become available at https://dpmjetiii.hepforge.org

slide-4
SLIDE 4

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 4

Hadron – hadron amplitudes in DPMJET (identical to PHOJET)

> Unitarization of soft, hard and diffractive amplitudes (together) via Eikonal scheme (AGK cutting rules)

Lowest order enhanced graphs

slide-5
SLIDE 5

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 5

Cross sections from optical theorem

> Free Regge parameters are fitted to data using total, elastic, inelastic and t-slope

soft: hard: diffractive (e.g.): not free! “free”?

> From we identify following relationships

CT14 LO

slide-6
SLIDE 6

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 6

Accurate fits for hadron-nucleon cross sections

> In an air-shower cascade all kinds of hadrons can re- interact > DPMJET is the only model working from particle production threshold up to highest cosmic ray energies (1021 GeV)

slide-7
SLIDE 7

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 7

DPMJET III in hA/AA/gA

Glauber approximation

> Each individual NN interaction is simulated with PHOJET (for A = 1 identical to PHOJET) > DPMJET tracks if valence quarks were involved (re-interaction of remnants) > Bookkeeping of used momentum fractions x for subsequent interactions > Partonic state from PHOJET is color connected but not fragmented > Strings localized between interaction partners in b-space

slide-8
SLIDE 8

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 8

Nuclear phenomenology and construction of (nuclear) remnants

Localized partonic strings String fusion + percolation (Nch↓ , <pT>↑) “fused” partonic strings Fragmentation (LUND/PYTHIA6) Fast and outbound directed particles escape “slow” particles re-interact with remnant “intranuclear cascade”

+Eex,t +Eex,p

Remaining excitation energy dissipated via Evaporation, Pre-equ… (link with FLUKA)

g De-

excitation

hep-ph/0403084

slide-9
SLIDE 9

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 9

Some AA results

Au – Au @ 130 GeV

F.W.Bopp, J. Ranft, R.Engel, S.Roesler arXiv:hep-ph/0403084v1

d – Au @ 200 GeV string fusion was needed, otherwise multiplicity 30-40% too high Best description of first p-Pb data

ALICE Coll.,PRL 110, 032301 (2013)

slide-10
SLIDE 10

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 10

MPI

> MPI (uncorrelated) from > Pomeron has its own partonic structure, i.e. pomeron-hadron, pomeron-pomeron scattering can have multiple parton interactions (recursion)

Arbitrary complex topologies

slide-11
SLIDE 11

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 11

New energy dependence of pT-cutoff

> Somewhat artificial separation of soft and hard interactions > Inspired by saturation in dipole models, e.g. Golec-Biernat and Wusthoff, PRD 59(1) 014017 (1999)

  • ld (empirical)

new GBW-type Is there a choice?

AF and R. Engel, CERN-Proceedings-2015-001

slide-12
SLIDE 12

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 12

Fits to cross sections

> normal: one good choice of pT-cut, σtot ~ sΔ > high: means low σQCD, σtot ~ sΔ > low: means high σQCD, σtot ~ ln2(s) bound by unitarity > All choices fit these observables “satisfactory” > Hard and diffractive cross section are related

slide-13
SLIDE 13

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 13

Tuned to minimum-bias

Fragmentation parameters are universal for hard and soft strings

slide-14
SLIDE 14

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 14

But…

Some models show characteristic flattening of the high Nch tail

slide-15
SLIDE 15

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 15

Deformation of multiplicity distribution from approaching the black disk limit

> Simplified model:

§ Only one pomeron graph with soft and hard part § Eikonal is only § Impact parameter shape is gaussian

> Width of b-shape scales with ln(s) for soft (constant for hard part) > Hard intercept fitted to minijet cross section ~0.2 > Without some mechanism to suppress nhard, the shape of multiplicty distribution will flatten at some point

slide-16
SLIDE 16

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 16

Cosmic ray applications

> For cosmic ray applications extrapolation to very high energies and p-Air interactions are crucial > Listed models are adjusted using only 7 TeV data > Model predictions diverge stronger for p+O ➡ input from LHC needed

  • T. Pierog, UHECR 2016
slide-17
SLIDE 17

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 17

Summary and conclusion

> DPMJET III 2016 is an improved general purpose and cosmic ray Monte Carlo > Re-tuned using minimum-bias data from all LHC experiments at 7 TeV > Air-shower cascade/cosmic ray extension: Ultra high energies, hadron-nucleon/nucleus cross sections > Problems with multiplicity distributions remain (as in many models) ➡ new phenomenology needed or impact parameter dependent cross sections > Next steps:

§ Finalize public versions (clean up, small technical things..) § Connect with LHC MC tools (Rivet/Professor, CRMC,..) § Air-shower simulations § Look at 13 TeV results!

slide-18
SLIDE 18

Anatoli Fedynitch | 8th MPI @ LHC, San Cristóbal de las Casas, Chiapas, Mexico | 2016/11/29 | Page 18

Where nuclear remnants matter

> Radiation protection! > Machine design

§ LHC accelerator elements § Magnets § Future accelerators § …

> Beam halo > Main application fields of FLUKA

Braun et al., Accelerators and Beams 17, 021006 (2014) NA49 @ 158 A GeV