Simulation in Particle Physics David Grellscheid What are the - - PowerPoint PPT Presentation

▶
simulation in particle physics
SMART_READER_LITE
LIVE PREVIEW

Simulation in Particle Physics David Grellscheid What are the - - PowerPoint PPT Presentation

Oct 18, 2022 579 likes •1.1k views

Simulation in Particle Physics David Grellscheid What are the fundamental building blocks of Nature? T-Shirt, explain rules here T-Shirt, explain rules here 4e candidate with m 4e = 124.6 GeV p T (electrons)= 24.9, 53.9, 61.9, 17.8 GeV m

slide-1
SLIDE 1

David Grellscheid

Simulation in Particle Physics

slide-2
SLIDE 2

What are the fundamental building blocks of Nature?

slide-3
SLIDE 3
slide-4
SLIDE 4
slide-5
SLIDE 5
slide-6
SLIDE 6
slide-7
SLIDE 7
slide-8
SLIDE 8
slide-9
SLIDE 9
slide-10
SLIDE 10

T-Shirt, explain rules here

slide-11
SLIDE 11

T-Shirt, explain rules here

slide-12
SLIDE 12
slide-13
SLIDE 13
slide-14
SLIDE 14
slide-15
SLIDE 15
slide-16
SLIDE 16
slide-17
SLIDE 17

ATLAS: Status of SM Higgs searches, 4/7/2012 37

4e candidate with m4e= 124.6 GeV

pT (electrons)= 24.9, 53.9, 61.9, 17.8 GeV m12= 70.6 GeV, m34= 44.7 GeV 12 reconstructed vertices

slide-18
SLIDE 18

Theory Experiment

slide-19
SLIDE 19

Theory Experiment ?

slide-20
SLIDE 20

Proton bunches 25 ns apart 10 protons per bunch

11

~25 collisions per crossing

slide-21
SLIDE 21

~100 million readout channels, every 25 ns

slide-22
SLIDE 22

CMS

slide-23
SLIDE 23

CMS After zeroes removed, 1.6 MB / event


slide-24
SLIDE 24

CMS After zeroes removed, 1.6 MB / event
 * 40 M events / s = 64 TB / s ?

slide-25
SLIDE 25

CMS After zeroes removed, 1.6 MB / event
 * 40 M events / s = 64 TB / s ?

Can’t save everything

slide-26
SLIDE 26

Incoming event rate per second Outgoing event rate per second Reduction factor Level 1 40 000 000 100 000 400 Level 2 100 000 3 000 30 Level 3 3 000 200 15

Trigger system to keep only interesting events

slide-27
SLIDE 27

Incoming event rate per second Outgoing event rate per second Reduction factor Level 1 40 000 000 100 000 400 Level 2 100 000 3 000 30 Level 3 3 000 200 15

Trigger system to keep only interesting events 200 events / s * 1.6 MB / event = 320 MB / s

slide-28
SLIDE 28

Incoming event rate per second Outgoing event rate per second Reduction factor Level 1 40 000 000 100 000 400 Level 2 100 000 3 000 30 Level 3 3 000 200 15

Trigger system to keep only interesting events 200 events / s * 1.6 MB / event = 320 MB / s

= ~ 3200 TB / year raw data

slide-29
SLIDE 29

Incoming event rate per second Outgoing event rate per second Reduction factor Level 1 40 000 000 100 000 400 Level 2 100 000 3 000 30 Level 3 3 000 200 15

Trigger system to keep only interesting events 200 events / s * 1.6 MB / event = 320 MB / s

= ~ 3200 TB / year raw data

Analysis is done offline,
 ~3000 collaboration members should have equal access to data worldwide

slide-30
SLIDE 30

[GeV]

4l

m

Events / 3 GeV

2 4 6 8 10 12 [GeV]

4l

m

Events / 3 GeV

2 4 6 8 10 12

Data Z+X *,ZZ

  • Z

=126 GeV

H

m

µ , 2e2 µ 7 TeV 4e, 4 µ , 2e2 µ 8 TeV 4e, 4

CMS Preliminary

  • 1

= 8 TeV, L = 5.26 fb s ;

  • 1

= 7 TeV, L = 5.05 fb s

[GeV]

4l

m

80 100 120 140 160 180

2011+2012%

ATLAS: Status of SM Higgs searches, 4/7/2012

μ

μ
slide-31
SLIDE 31

[GeV]

4l

m

Events / 3 GeV

2 4 6 8 10 12 [GeV]

4l

m

Events / 3 GeV

2 4 6 8 10 12

Data Z+X *,ZZ

  • Z

=126 GeV

H

m

µ , 2e2 µ 7 TeV 4e, 4 µ , 2e2 µ 8 TeV 4e, 4

CMS Preliminary

  • 1

= 8 TeV, L = 5.26 fb s ;

  • 1

= 7 TeV, L = 5.05 fb s

[GeV]

4l

m

80 100 120 140 160 180

2011+2012%

Re

  • ATLAS: Status of SM Higgs searches, 4/7/2012

μ

μ
slide-32
SLIDE 32

[GeV]

4l

m

Events / 3 GeV

2 4 6 8 10 12 [GeV]

4l

m

Events / 3 GeV

2 4 6 8 10 12

Data Z+X *,ZZ

  • Z

=126 GeV

H

m

µ , 2e2 µ 7 TeV 4e, 4 µ , 2e2 µ 8 TeV 4e, 4

CMS Preliminary

  • 1

= 8 TeV, L = 5.26 fb s ;

  • 1

= 7 TeV, L = 5.05 fb s

[GeV]

4l

m

80 100 120 140 160 180

2011+2012%

Re

  • Need to get theory

predictions.

ATLAS: Status of SM Higgs searches, 4/7/2012

μ

μ
slide-33
SLIDE 33
slide-34
SLIDE 34
slide-35
SLIDE 35

Monte Carlo event generators

Re

slide-36
SLIDE 36
  • General purpose MC event generator
  • 30-year history in its F77 implementation; 


Herwig 7 is a complete redesign from ground up in 
 C++, started 2003. (more and more python creeping in)

  • currently ~15 collaboration members


in Durham, Karlsruhe, Manchester, Zurich

  • Coordination is hard. Technically: Trello / Trac / Jenkins.


Don’t have a good model sociologically.

  • Main reference: arXiv:0803.0883, 1101.2599


https://herwig.hepforge.org/

Herwig details

slide-37
SLIDE 37

ThePEG Toolkit for high energy physics event generation Herwig++ Box of physics implementations Each building block is
 a compiled C++ class

[Leif Lönnblad]

slide-38
SLIDE 38

ThePEG Repository plaintext setup files no more compilation needed here

slide-39
SLIDE 39

Default Setup

EventGenerator RNG Physics Model EventHandler

Decays Hadron-
 ization Parton 
 Shower

Hard
 Subprocess

slide-40
SLIDE 40

Default Setup

EventGenerator RNG Physics Model EventHandler

Decays Hadron-
 ization Parton 
 Shower

Hard
 Subprocess

create ThePEG::StandardEventHandler /Herwig/LHCHandler set LHCHandler:LuminosityFunction FixedLHCLuminosity insert LHCHandler:SubProcessHandlers[0] /Herwig/SimpleQCD set LHCHandler:CascadeHandler /Herwig/ShowerHandler set LHCHandler:HadronizationHandler /Herwig/ClusterHadHandler set LHCHandler:DecayHandler /Herwig/DecayHandler [...] create ThePEG::EventGenerator /Herwig/LHCGenerator ThePEG.so
 set LHCGenerator:EventHandler /Herwig/LHCHandler
 [...] set LHCHandler:BeamA /Herwig/Particles/p+ set LHCHandler:BeamB /Herwig/Particles/p+ set FixedLHCLuminosity:Energy 14000.0 [...]

slide-41
SLIDE 41

Default Setup

EventGenerator RNG Physics Model EventHandler

Decays Hadron-
 ization Parton 
 Shower

Hard
 Subprocess

create ThePEG::StandardEventHandler /Herwig/LHCHandler set LHCHandler:LuminosityFunction FixedLHCLuminosity insert LHCHandler:SubProcessHandlers[0] /Herwig/SimpleQCD set LHCHandler:CascadeHandler /Herwig/ShowerHandler set LHCHandler:HadronizationHandler /Herwig/ClusterHadHandler set LHCHandler:DecayHandler /Herwig/DecayHandler [...] create ThePEG::EventGenerator /Herwig/LHCGenerator ThePEG.so
 set LHCGenerator:EventHandler /Herwig/LHCHandler
 [...] set LHCHandler:BeamA /Herwig/Particles/p+ set LHCHandler:BeamB /Herwig/Particles/p+ set FixedLHCLuminosity:Energy 14000.0 [...]

Arbitrary user extensions use dlopen():

create DGrell::Myclass /DGrell/Myclass DGrellHwPlugin.so

Main code never needs recompilation.

slide-42
SLIDE 42

Default Setup

EventGenerator RNG Physics Model EventHandler

Decays Hadron-
 ization Parton 
 Shower

Hard
 Subprocess

create ThePEG::StandardEventHandler /Herwig/LHCHandler set LHCHandler:LuminosityFunction FixedLHCLuminosity insert LHCHandler:SubProcessHandlers[0] /Herwig/SimpleQCD set LHCHandler:CascadeHandler /Herwig/ShowerHandler set LHCHandler:HadronizationHandler /Herwig/ClusterHadHandler set LHCHandler:DecayHandler /Herwig/DecayHandler [...] create ThePEG::EventGenerator /Herwig/LHCGenerator ThePEG.so
 set LHCGenerator:EventHandler /Herwig/LHCHandler
 [...] set LHCHandler:BeamA /Herwig/Particles/p+ set LHCHandler:BeamB /Herwig/Particles/p+ set FixedLHCLuminosity:Energy 14000.0 [...]

slide-43
SLIDE 43

Default Setup

EventGenerator RNG Physics Model EventHandler

Decays Hadron-
 ization Parton 
 Shower

Hard
 Subprocess

slide-44
SLIDE 44

Event record flow

Decays Hadron-
 ization Parton 
 Shower

Hard
 Subprocess

Herwig
 Analysis

∅

slide-45
SLIDE 45

Event record flow

Decays Hadron-
 ization Parton 
 Shower

Hard
 Subprocess

Herwig
 Analysis Rivet LHC
 expts

  • ther


users

HepMC

∅

slide-46
SLIDE 46

Event record flow

Hard
 Subprocess

Decays Hadron-
 ization Parton 
 Shower Herwig
 Analysis Rivet LHC
 expts

  • ther


users

HepMC

∅

slide-47
SLIDE 47

Event record flow

Hard
 Subprocess

Decays Hadron-
 ization Parton 
 Shower Herwig
 Analysis Rivet LHC
 expts

  • ther


users

HepMC

2→2 ME
 construction

LHE
 file reader

hand-coded
 MEs


(hel. amplitudes)

∅

slide-48
SLIDE 48

Event record flow

Hard
 Subprocess

Decays Hadron-
 ization Parton 
 Shower Herwig
 Analysis Rivet LHC
 expts

  • ther


users

HepMC

2→2 ME
 construction

LHE
 file reader

hand-coded
 MEs


(hel. amplitudes)

∅

UFO

slide-49
SLIDE 49

Simulated data sets of millions of events

slide-50
SLIDE 50

August 1991

World Wide Web

The WorldWideWeb (W3) is a wide-area hypermedia information retrieval initiative aiming to give universal access to a large universe of documents. Everything there is online about W3 is linked directly or indirectly to this document, including an executive summary of the project, Mailing lists , Policy , November's W3 news , Frequently Asked Questions . What's out there? Pointers to the world's online information, subjects , W3 servers, etc. Help

  • n the browser you are using

Software Products A list of W3 project components and their current state. (e.g. Line Mode ,X11 Viola , NeXTStep , Servers , Tools , Mail robot , Library ) Technical Details of protocols, formats, program internals etc Bibliography Paper documentation on W3 and references. People A list of some people involved in the project. History A summary of the history of the project. How can I help ? If you would like to support the web.. Getting code Getting the code by anonymous FTP , etc.