1 | <html> |
---|
2 | <head> |
---|
3 | <title>MadGraph 5 Processes</title> |
---|
4 | <link rel="stylesheet" type="text/css" href="pythia.css"/> |
---|
5 | <link rel="shortcut icon" href="pythia32.gif"/> |
---|
6 | </head> |
---|
7 | <body> |
---|
8 | |
---|
9 | <h2>MadGraph 5 Processes</h2> |
---|
10 | |
---|
11 | By far the easiest way to implement new processes into PYTHIA 8 is |
---|
12 | by using the matrix-element generator MadGraph 5. This program has |
---|
13 | an option to output the results of a matrix-element calculation |
---|
14 | as a set of PYTHIA 8 C++ classes (plus further auxiliary code), |
---|
15 | that can then be linked and used as |
---|
16 | <a href="SemiInternalProcesses.html" target="page">semi-internal</a> processes, |
---|
17 | meaning they are handled identically with normal internal ones. |
---|
18 | This way, MadGraph 5 can be used to implement processes from |
---|
19 | any model that can be written in terms of a Lagrangian. Any |
---|
20 | <i>2 -> 1</i>, <i>2 -> 2</i> and <i>2 -> 3</i> processes |
---|
21 | can be implemented, the limit being set by the absence of efficient |
---|
22 | phase space generator algorithms for higher multiplicities in |
---|
23 | PYTHIA. Features such as <i>s</i>-channel resonances are |
---|
24 | automatically implemented in the process classes. Besides the process |
---|
25 | library and necessary model files, also an example main program is |
---|
26 | generated for each set of processes, which can be easily modified to |
---|
27 | perform the desired analyses. |
---|
28 | |
---|
29 | <p/> |
---|
30 | In order to create a PYTHIA 8 process library with MadGraph 5, first |
---|
31 | download the MadGraph 5 package from |
---|
32 | <a href="https://launchpad.net/madgraph5" target="page"> |
---|
33 | https://launchpad.net/madgraph5</a>, and untar the package. You can |
---|
34 | then specify the location of your <code>pythia81xx</code> directory |
---|
35 | in the file <code>input/mg5_configuration.txt</code>: |
---|
36 | <br/><code>pythia8_path = ./pythia81xx</code> |
---|
37 | <br/>The location can be either relative (to the directory |
---|
38 | <code>MadGraph5_v_x_x_x/.</code>) or absolute. |
---|
39 | |
---|
40 | <p/> |
---|
41 | For any model that is already implemented in the MadGraph 5 package, |
---|
42 | you can directly use the model. Start the MadGraph 5 interface |
---|
43 | <code>bin/mg5</code>, and do: |
---|
44 | <pre> |
---|
45 | import model model_name |
---|
46 | generate your_process_in_mg5_syntax |
---|
47 | add process your_next_process_in_mg5_syntax |
---|
48 | ... |
---|
49 | output pythia8 [path_to_pythia81xx_directory] |
---|
50 | </pre> |
---|
51 | |
---|
52 | <p/> |
---|
53 | For examples of MG5 process syntax, please see |
---|
54 | <a href="http://madgraph.phys.ucl.ac.be/EXAMPLES/example_mg5.html" |
---|
55 | target="page">http://madgraph.phys.ucl.ac.be/EXAMPLES/example_mg5.html</a> |
---|
56 | or type <code>help generate</code>. If you specified the path to the |
---|
57 | <code>pythia81xx</code> directory in the <code>mg5_configuration</code> |
---|
58 | file, you do not need to enter it in the <code>output</code> command. |
---|
59 | |
---|
60 | <p/> |
---|
61 | If your preferred model is found on the FeynRules model wiki page, |
---|
62 | <a href="http://feynrules.irmp.ucl.ac.be/wiki/ModelDatabaseMainPage" |
---|
63 | target="page">http://feynrules.irmp.ucl.ac.be/wiki/ModelDatabaseMainPage</a>, |
---|
64 | download the UFO (Universal FeynRules Output) tar file for the model, |
---|
65 | untar in the <code>models/</code> directory, and use as above. |
---|
66 | |
---|
67 | <p/> |
---|
68 | If you want to implement a new model which has not yet been implemented, |
---|
69 | you can do this either using the Mathematica package FeynRules (see |
---|
70 | <a href="http://feynrules.irmp.ucl.ac.be/" target="page"> |
---|
71 | http://feynrules.irmp.ucl.ac.be/</a>) or directly edit the UFO model |
---|
72 | files of the most similar model in the <code>models/</code> directory. |
---|
73 | |
---|
74 | <p/> |
---|
75 | The resulting output from the <code>output pythia8</code> command is: |
---|
76 | <ul> |
---|
77 | <li>A process directory <code>Processes_modelname</code> with the |
---|
78 | model information and the files needed for all processes defined for |
---|
79 | this model, placed in the <code>pythia81xx</code> main directory. |
---|
80 | The model files are <code>Parameters_modelname.h/cc</code> and |
---|
81 | <code>HelAmps_modelname.h/cc</code>, and the process files for each |
---|
82 | process class (with the same mass, spin and color of the initial/final |
---|
83 | state particles) are called <code>Sigma_modelname_processname.h/cc</code>. |
---|
84 | The directory also contains a <code>makefile</code> and a model parameter |
---|
85 | file <code>param_card_modelname.dat</code>.</li> |
---|
86 | <li>An example main program in the directory <code>examples/</code> |
---|
87 | (in the <code>pythia81xx</code> main directory) called |
---|
88 | <code>main_modelname_N.cc</code> and a corresponding makefile |
---|
89 | <code>Makefile_modelname_N</code>. This main program links in the |
---|
90 | process classes in the process directory described above. To run the |
---|
91 | example main program, just go to the <code>examples/</code> |
---|
92 | directory and run |
---|
93 | <br/><code>make -f Makefile_modelname_N</code> |
---|
94 | <br/>or run <code>launch</code> directly inside the MadGraph 5 |
---|
95 | command line interface.</li> |
---|
96 | </ul> |
---|
97 | |
---|
98 | <p/> |
---|
99 | Note that in order for PYTHIA 8 to be able to automatically decay any |
---|
100 | new particles, it is necessary to specify the branching ratios of the |
---|
101 | particles in the <code>param_card</code> file, see |
---|
102 | [<a href="Bibliography.html" target="page">Ska04,Alw07</a>] for details. |
---|
103 | |
---|
104 | <p/> |
---|
105 | For further technical details, please see the MadGraph 5 release paper |
---|
106 | [<a href="Bibliography.html" target="page">Alw11</a>] and the |
---|
107 | <a href="SemiInternalProcesses.html" target="page">semi-internal</a> processes page. |
---|
108 | |
---|
109 | <p/> |
---|
110 | Of course, as with MadGraph 4, MadGraph 5 can also output |
---|
111 | files of parton-level events according to the |
---|
112 | <a href="LesHouchesAccord.html" target="page">LHEF</a> standard, |
---|
113 | that can be read in and processed further by PYTHIA 8. |
---|
114 | The advantage is that then the MadGraph 5 phase space generator |
---|
115 | can be used, which opens up for processes with more than three |
---|
116 | particles in the final state. The disadvantages are that it is less |
---|
117 | easy to mix and match with existing PYTHIA processes, and that one |
---|
118 | needs to regenerate and store large LHEF files for different |
---|
119 | kinematics cuts or parameter values. |
---|
120 | |
---|
121 | <p/> |
---|
122 | Please cite the MadGraph 5 release paper [<a href="Bibliography.html" target="page">Alw11</a>] if you use |
---|
123 | MadGraph 5 to generate process libraries for PYTHIA 8. |
---|
124 | |
---|
125 | </body> |
---|
126 | </html> |
---|
127 | |
---|
128 | <!-- Copyright (C) 2012 Torbjorn Sjostrand --> |
---|