MadGraph5. Olivier Mattelaer CP3/FNRS

Transcription

1 MadGraph5 Olivier Mattelaer CP3/FNRS MG5: J. Alwall / M. herquet / F. Maltoni / OM /T. Stelzer ALOHA: P. Aquino / W. Link / F. Maltoni / OM / T. Stelzer UFO: C. Degrande / C. Duhr / B. Fuks / D. Grellsheid/ OM T. Reiter AMC@NLO: J. Alwall, R.Frederix, S. Frixione, V. Hirschi, F.Maltoni, OM, P. Torrielli, M. Zaro and a lot of external collaborators

2 Plan BSM Core QCD v5 Tools 2

3 Plan BSM Core QCD v5 Tools 3

4 What is MadGraph5? Diagram generator in python. returns code to compute matrix-element based on the Helicity Amplitude formalism Various output format MadEvent: Leading-order cross-section and event generation NLO cross-section and event generation (matched to the shower) Pythia8: export the matrix element inside the pythia8 framework Tools: MadSpin, MadWeight, MadDM,... 4

5 Recap: What is MG5 Can handle ANY model (that can be written as a Lagrangian), conveniently output by FeynRules Any Lorentz structure for any spin (up to 2) and color (6tets, ε ijk ) Multiparticle vertices for any multiplicity, multifermion vertices Super fast process generation Unlimited-length decay chains with full BW and spin effects Event generation speedups by orders of magnitude Super-user-friendly command line interfaces 5

6 Core news Lots of speedups and improvements, including Huge speedup of gridpacks Vast speedup for long decay chains with multiparticle decays Huge improvements in user interface Multi cluster support New 4-flavor matching and VBF-type matching 4 fermion vertices in FR+MG5 (except Majorana) Spin 3/2 particles in FR+MG5 Complex mass scheme Feynman gauge Handling of negative weights On-the-fly 2-body decay width calculations ( Auto width ) 6

7 Plan UFO / ALOHA BSM Core QCD v5 Tools 7

8 UFO = Universal Feynrules Output New Model Format Gosam/ Herwig++/ MG5 Fully generic color/lorentz/... Automatic Creation of HELAS routine for ANY BSM theory Fortran / C++/Python 8

9 ALOHA ALOHA UFO Helicity Brussels October Tim Stelzer

10 ALOHA ALOHA UFO Helicity Options: Standard (HELAS) Feynman gauge Complex-mass scheme Numerator (for Loop) Open-Loop Brussels October Tim Stelzer

11 Any BSM should be possible in a fully automatic and efficient way! Some restriction applies: Only local theory Theory should respect CPT and lorentz invariance (all indices should be contracted) Color supported up to dimension 8 (including sextet and epsilon structure) Spin supported up to spin 2 (including spin3/2) No four fermion interaction with fermion-flow violation / majorana in the same model 10

12 Plan BSM Core QCD v5 Tools 11

13 ) Tools Tools Utility Progress MadAnalysis5 Plotting distributions -1 = 10 fb int # mu ( L Released PT [ mu ] (GeV/c) [E.Conte, B. Fuks: CPC 184 (2013) ] 12

14 ) Tools Tools Utility Progress MadAnalysis5 Plotting distributions -1 = 10 fb int # mu ( L Released PT [ mu ] (GeV/c) MadDecay Automatic width computation 95% Decay with full (LO) spin- [J.Alwall, C.Duhr, B.Fuks, OM, D.Ozturk, CH Shen] 12

15 ) Tools Tools Utility Progress MadAnalysis5 Plotting distributions -1 = 10 fb int # mu ( L Released PT [ mu ] (GeV/c) MadDecay Automatic width computation 95% MadSpin Decay with full (LO) spincorrelation Released [P. Artoisenet, R. Frederix, OM, R. Rietkerk: ] 12

16 ) Tools Tools Utility Progress MadAnalysis5 Plotting distributions -1 = 10 fb int # mu ( L Released PT [ mu ] (GeV/c) MadDecay Automatic width computation 95% MadSpin re-weight module Decay with full (LO) spincorrelation add weight to event according to a new theoretical hypothesis Specific integrator for the Matrix Released golden master [OM] 12

17 ) Tools Tools Utility Progress MadAnalysis5 Plotting distributions -1 = 10 fb int # mu ( L Released PT [ mu ] (GeV/c) MadDecay Automatic width computation 95% MadSpin re-weight module MadWeight5 Decay with full (LO) spincorrelation add weight to event according to a new theoretical hypothesis Specific integrator for the Matrix Element Method Released golden master in beta [P.Artoisenet, OM] 12

18 ) Tools Tools Utility Progress MadAnalysis5 Plotting distributions -1 = 10 fb int # mu ( L Released PT [ mu ] (GeV/c) MadDecay Automatic width computation 95% MadSpin re-weight module MadWeight5 MadDM Decay with full (LO) spincorrelation add weight to event according to a new theoretical hypothesis Specific integrator for the Matrix Element Method Relic density and direct detection limit for any BSM model [M. Backovic, K.C Kong, M. McCaskey] Released golden master in beta Released 12

19 Plan BSM Core QCD v5 Tools Matching/Merging 13

20 Matching in MG + Pythia6 Alwall et al. [arxiv: ], Alwall, de Visscher, Maltoni [arxiv: ], MadGraph wiki MLM-style matching Classic cone-jet matching (a la AlpGen), k T -jet matching, and shower-k T matching Official CMS SM background simulation Excellent agreement with data across the line Fully supports matching for any NP process Restriction: Matching of jets from core process, not from decay products (let Pythia ME corr. take care of that) 14

21 MLM algorithm 1. Generate ME events (with different parton multiplicities) using parton-level cuts (p T ME /ΔR or k T ME ) 2. Cluster each event and reweight α s and PDFs based on the scales in the corresponding clustering vertices 3. Run the parton shower with starting scale t 0 = m T. 4. Check that the number of jets after parton shower is the same as ME partons, and that all jets after parton shower are matched to the ME partons at a scale Q match. If yes, keep the event. If no, reject the event. Q match is called the matching scale. 5. For highest multiplicity, allow radiation < lowest ME scale 15

22 Matching with b-quarks p MT (Z)M T (b) for factorization scale Gives smooth matching to Pythia PS Shower ME 16

23 MadGraph FKS FKS CutTools 17

24 18

25 Why automation? Time: Less tools, means more time for physics Robust: Easier to test, to trust Easy: One framework/tool to learn 18

26 Why automation? Time: Less tools, means more time for physics Robust: Easier to test, to trust Easy: One framework/tool to learn Why NLO? Reliable prediction of the total rate Reduction of the theoretical uncertainty 18

27 Why automation? Time: Less tools, means more time for physics Robust: Easier to test, to trust Easy: One framework/tool to learn Why NLO? Reliable prediction of the total rate Reduction of the theoretical uncertainty Why matched to the PS? Parton are not an detector observables Matching cure some fix-order ill behaved observables 18

28 DEMO Is it really automatic? 19

29 DEMO 1) Download the code 20

30 launch the code [./bin/mg5] Exactly like MG5!!! 21

31 You can enter ANY process! add [QCD] for NLO functionalities generate p p > t t~ [QCD] generate p p > e+ e- mu+ mu- [QCD] generate p p > w+ w- j j [QCD] 22

32 Create your code output PATH Run it: launch [PATH] 23

33 Create your code output PATH Run it: launch [PATH] 23

34 Create your code output PATH Run it: launch [PATH] First Question: 23

35 Create your code output PATH Run it: launch [PATH] Second Question: 23

36 Compilation Check Poles cancelation 24

37 Compilation Check Poles cancelation 24

38 Compilation Check Poles cancelation 24

39 Integration Events Generation 25

40 Unweight Events Main Results The Shower Next release will be much faster 4 times faster for this simple process (so around 4 min) 26

41 DEMO Is it really automatic? 27

42 DEMO Is it really automatic? As much as LO! 27

43 Conclusion BSM Core QCD v5 Tools 28

44 Conclusion BSM Core MG5 is available on the web Big improvement compare to MG4 QCD v5 Tools 28

45 Conclusion BSM Core MG5 is available on the web Big improvement compare to MG4 QCD v5 Tools Any BSM theory can be runned in MG5 in a fully automatic way 28

46 Conclusion BSM Core MG5 is available on the web Big improvement compare to MG4 QCD v5 Tools Any BSM theory can be runned in MG5 in a fully automatic way Various tools available to help the community 28

47 Conclusion BSM Core MG5 is available on the web Big improvement compare to MG4 QCD v5 Tools Any BSM theory can be runned in MG5 in a fully automatic way Various tools available to help the community Fully automatic NLO generator 28

48 Plan BSM Core QCD v5 Tools DECAY MADWEIGHT MadAnalysis5 29

49 MadWeight Tools for the Matrix Element Method Reweighting based on the matrix-element Technique used for the top mass measurment Now Included in MG5 (still in beta) Include ISR corrections Better multi-process treatment 30

50 MadWeight Tools for the Matrix Element Method Reweighting based on the matrix-element P(p vis )= 1 vis d dx 1 dx 2 M (p) 2 W (p, p vis ) Technique used for the top mass measurment Now Included in MG5 (still in beta) Include ISR corrections Better multi-process treatment 30

51 MadDECAY MG is now able to create the decay table for pythia. But the User should provide all the possible decay This project will allow to find the relevant channel of integration automatically by checking the relevance of one additional decay level 3 3 ve~ 4 b~ w- w- e- 1 h 1 h 1 h 3 w+ b w+ 2 2 No Phase-space 2 Possibility to add a module for decaying the final state particles with either a part OR the full spin-correlation. For ANY BSM theories 31

52 Computing the widths mg5> compute_widths Z First evaluate 2>2 and 2>3 contribution Compute ONLY the relevant contribution Write the new param_card.dat 32

53 MadWeight Associate to each experimental event characterised by p vis, the probability P(p vis ) to be produced and observed following a theoretical assumption 33

54 MadWeight Associate to each experimental event characterised by p vis, the probability P(p vis ) to be produced and observed following a theoretical assumption 33

55 MadWeight Associate to each experimental event characterised by p vis, the probability P(p vis ) to be produced and observed following a theoretical assumption is the squared matrix element 33

56 MadWeight Associate to each experimental event characterised by p vis, the probability P(p vis ) to be produced and observed following a theoretical assumption is the squared matrix element is the transfer function 33

57 MadWeight Associate to each experimental event characterised by p vis, the probability P(p vis ) to be produced and observed following a theoretical assumption is the squared matrix element is the transfer function is the phase-space integral 33

58 MadWeight Associate to each experimental event characterised by p vis, the probability P(p vis ) to be produced and observed following a theoretical assumption P(p vis )= 1 vis d dx 1 dx 2 M (p) 2 W (p, p vis ) M (p) 2 d dx 1 dx 2 vis is the squared matrix element W (p, p vis ) is the transfer function is the phase-space integral is the cross-section (after cuts) 33

59 Introduction My first analysis Advanced functionalities Developer-friendly mode Summary MadAnalysis overview Scope Analysis of event files produced by Monte Carlo tools at parton level, hadron level or after detector simulation. Definition of various selection cuts on the input samples. Production of histograms for di erent distributions. Results of the analysis summed up by a S/B-like ratio table. Computing details Interface written in Python and ROOT ; kernel in C++. Possible output in ROOT, HTML, L A TEX. Website https : //server 06.fynu.ucl.ac.be/projects/madanalysis Please send us your comments and suggestions (tickets on the wiki). 34

60 Introduction My first analysis Advanced functionalities Developer-friendly mode Summary Step 5 : diplaying results # 35