5 | | The tutorial consists in implementing the so-called ''Hill Model'', described in [http://arxiv.org/abs/hep-ph/0603082 hep-ph/0603082], into !MadGraph, and to simulate the full chain down to PGS level events for a sample of ''tt~h'' events at the !LHC. A very brief description of the model can be found [attachment:MC4BSM_FeynRules.pdf here]. |
6 | | |
7 | | * The !MC4BSM [http://particle.physics.ucdavis.edu/workshops/doku.php?id=2009:mc4bsm_wishlist/ website]. |
8 | | |
9 | | === Implementing the model into !FeynRules === |
10 | | |
11 | | [http://feynrules.phys.ucl.ac.be/ FeynRules] is a Mathematica® package that allows the calculation of Feynman rules in momentum space for any !QFT physics model. The Feynman rules calculated by the code can then be used to implement the new physics model into other existing tools, such as MC generators. |
12 | | |
13 | | Usefull material for the tutorial: |
14 | | |
15 | | * A ready-to-go !FeynRules model file for the Hill model: [attachment:HillModel.fr.tar.gz HillModel.fr.tar.gz]. |
16 | | * An example notebook of how to run !FeynRules: [attachment:MC4BSM.nb.tar.gz MC4BSM.nb]. |
17 | | * A short introduction to implementing new models into !FeynRules, and how to run the !MadGraph interface: [attachment:MC4BSM_FeynRules.pdf MC4BSM_FR.pdf]. |
18 | | |
19 | | === Generating parton level events with MG/ME === |
20 | | |
21 | | After the model has been implemented into !FeynRules, it can easily be exported to !MadGraph via the corresponding !FeynRules interface. The model can then be used just like any other built-in !MadGraph model. |
22 | | |
23 | | In this tutorial we use the !FeynRules implementation of the Hill model to generate a sample of ''pp>tt~h1'' events at the !LHC @ 14TeV. |
24 | | |
25 | | * A short introduction of how to import !FeynRules models into MG/ME: [attachment:MG_MC4BSM.pdf MC4BSM_MG.pdf]. |
26 | | |
27 | | === Decaying the BSM particles with BRIDGE === |
28 | | |
29 | | We can use [http://www.lepp.cornell.edu/public/theory/BRIDGE/ BRIDGE] to |
30 | | |
31 | | * Compute the branching ratios and the decay tables for the BSM particles contained in the Hill model. |
32 | | * Decay the ''tt~h1'' final state obtained in the previous step. |
33 | | |
34 | | === Pythis and PGS === |
35 | | |
36 | | Finally, we can run [http://home.thep.lu.se/~torbjorn/Pythia.html Pythia] and [http://www.physics.ucdavis.edu/~conway/research/software/pgs/pgs.html PGS] to obtained hadron-level and detector level events. |
| 10 | === Install on Windows === |
| 11 | The different program are foreseen to run on !Linux/Mac. If you don't have Linux or Mac, the best solution is to install Linux (I would advice to install Ubuntu) in dual boot and follow the instructions for Linux. |
| 12 | This video might help you to do so: http://video.google.com/videoplay?docid=-2369893842637434537 |
| 13 | virtual machine are fine, but quite often the library readline is not present on the system. This package is not mandatory but allozs mg5 to have history and nice auto-completions. |
| 14 | So you might consider to install that library first and the compile python2.6. |
| 17 | === Install on Linux === |
| 18 | |
| 19 | * '''python 2.6/2.7'''[[BR]] |
| 20 | !MadGraph 5 requires python 2.6 or python 2.7. In order to check your version of python, you can type: |
| 21 | python --version. Some of the repositories are still on python 2.5. So if updating python via your repository manager didn't work. |
| 22 | you can download directly python from the following link: |
| 23 | http://www.python.org/download/ |
| 24 | and follow instructions. |
| 25 | * '''Madgraph5'''[[BR]] |
| 26 | You will find madgraph5 package on the following page: https://launchpad.net/madgraph5. |
| 27 | For this program, you just need to untar it. |
| 28 | To check if mg5 is correctly install you directly try to run it by doing: |
| 29 | ./bin/mg5 |
| 30 | If you don't have a valid python version. It will fail directly. |
| 31 | * '''MadAnalysis'''[[BR]] |
| 32 | This programs is not mandatory but allows to have partonic distributions done automatically. |
| 33 | In order to install this program, you need to launch mg5 (./bin/mg5) and then type |
| 34 | {{{ |
| 35 | install MadAnalysis |
| 36 | }}} |
| 37 | |
| 38 | === Install on !MacOs === |
| 39 | |
| 40 | 1. '''Basic program''' |
| 41 | * perl |
| 42 | * bash [[BR]] |
| 43 | Those programs are by default present in any Mac distribution, |
| 44 | but just in case you should check that they are all installed in your computer (by doing perl --version). |
| 45 | * gmake [[BR]] |
| 46 | By default, makefile are not recognize on mac. In order to check, you can do 'make --version'. If it says: |
| 47 | -bash: make: command not found |
| 48 | then you should install it. The easiest is to install the xcode. |
| 49 | 1. For !MacOs 10.5: https://connect.apple.com/cgi-bin/WebObjects/MemberSite.woa/wa/getSoftware?bundleID=20414 |
| 50 | 2. For !MacOs 10.6: http://connect.apple.com/cgi-bin/WebObjects/MemberSite.woa/wa/getSoftware?bundleID=20792 |
| 51 | 3. For !MacOs 10.7: http://itunes.apple.com/us/app/xcode/id448457090?mt=12 |
| 52 | In order to be able to download it, you will need to a apple developer account (which is free) |
| 53 | * gfortran 4.x (To check the version do gfortran do: gfortran --version) |
| 54 | This program is often not included in standard mac version. In order to download it, you can download it from: |
| 55 | a. http://prdownloads.sourceforge.net/hpc/gcc-lion.tar.gz?download for !MacOs 10.7 |
| 56 | b. http://prdownloads.sourceforge.net/hpc/gcc-snwleo-intel-bin.tar.gz?download for !MacOs 10.6 |
| 57 | c. http://sourceforge.net/projects/hpc/files/hpc/gcc/gcc-leopard-intel-bin.tar.gz/download for !MacOs 10.5 |
| 58 | |
| 59 | 2. '''!MadGraph5''' |
| 60 | * '''python''' |
| 61 | a. The default python version of !MacOs 10.6 is python 2.6 and is (in 99% of the case) fully working with Madgraph5, so nothing to install for you. |
| 62 | b. The python version of !MacOs 10.7 is python 2.7. And this one works perfectly, so nothing to install for you. |
| 63 | c. For other version of MacOS, you will need to install python2.7 from the python web site: http://www.python.org/download/ |
| 64 | * '''Madgraph5'''[[BR]] |
| 65 | You will find madgraph5 package on the following page: https://launchpad.net/madgraph5. |
| 66 | For this program, you just need to untar it. |
| 67 | To check if mg5 is correctly install you directly try to run it by doing: |
| 68 | ./bin/mg5 |
| 69 | If you don't have a valid python version. It will fail directly. |
| 70 | * '''MadAnalysis'''[[BR]] |
| 71 | This programs is not mandatory but allows to have partonic distributions done automatically. |
| 72 | In order to install this program, you need to launch mg5 (./bin/mg5) and then type |
| 73 | {{{ |
| 74 | install MadAnalysis |
| 75 | }}} |
| 76 | |
| 77 | === pre-workshop exercises === |
| 78 | As a pre-workshop exercises and a proof that the full instalation is succesfull. |
| 79 | Please launch mg5 |
| 80 | ./bin/mg5 |
| 81 | and type the following command: |
| 82 | {{{ |
| 83 | tutorial |
| 84 | }}} |
| 85 | and follow the tutorial instructions on the screen. |
42 | | -- Main.ClaudeDuhr - 14 Apr 2009 |
| 90 | * If you don't have run the tutorial before the lectures please do it now. |
| 91 | Please launch mg5 |
| 92 | ./bin/mg5 |
| 93 | and type the following command: |
| 94 | {{{ |
| 95 | tutorial |
| 96 | }}} |
| 97 | and follow the tutorial instructions on the screen. |
| 98 | |
| 99 | * The exercice of this tutorial is to |
| 100 | 1. Use the model generatedby Feynrules |
| 101 | 2. generate $ p p \to U \bar U $ |
| 102 | 3. generate some events |
| 103 | 4. pass the events to Pythia |
| 104 | 5. Restart the chain but this time by including all the decay chain in MG5. |
| 105 | |
| 106 | Note that in order to have the correct cross-sections in the second cases or to have the correct branching ratio computed by Pythia in the first case, we need first to compute the width and the branching ratio of the non standard model particles. |
| 107 | |
| 108 | * '''The model'''[[BR]] |
| 109 | The simple way to have access to a model in MG5 is to put in the directory: MG5_DIR/models |
| 110 | after that you can simply import it in MG5 by doing |
| 111 | {{{ |
| 112 | import model MODELNAME |
| 113 | }}} |
| 114 | or |
| 115 | {{{ |
| 116 | import model MODELNAME --modelname |
| 117 | }}} |
| 118 | The option --modelname tells MG5 to use the name of the particles defines in the UFO model, and not the usual MG5 conventions for the particles of the SM and the MSSM. In principle this shouldn't change anything in the model that you load. |
| 119 | For the first time that you import a new model, it's allow a good points to test it. For this |
| 120 | 1. It's advice to run MG5 in debug mode (./bin/mg5 --debug) |
| 121 | 2. It's useful to test the Lorentz invariance/Gauge invariance/... of at least a couple of processes. This is one examples: |
| 122 | {{{ |
| 123 | import model MODELNAME |
| 124 | check p p > urv urv~ |
| 125 | }}} |
| 126 | |
| 127 | * ''' Computation of the width and branching ratio''' |
| 128 | The first element in order to have a valid generations is to have valid input parameters. This includes in particular to have the correct width (and the correct branching ratio for Pythia). |
| 129 | If you tell all the channel of desintegration, MG5 is able to compute all those parameters. |
| 130 | {{{ |
| 131 | import model MODELNAME |
| 132 | generate urv > u p1 |
| 133 | add process urv > u p2 |
| 134 | add process [ADD all 1->2 decay and relevant 1->3 decay, for both interesting particles] |
| 135 | output |
| 136 | launch |
| 137 | }}} |
| 138 | This will creates to you one valid param_card with the computed width and the computed branching ratio. |
| 139 | Note that the solution (and the valid files to produces it are in attachment of this page) |
| 140 | Note that particles which doesn't have a specify decays keeps the width defines by the UFO Model. This might be problematic for p1 (If you don't set the width of p1 to zero in the model) |
| 141 | |
| 142 | * '''Generation of events (Decay done in pythia)''' |
| 143 | {{{ |
| 144 | import model MODELNAME |
| 145 | generate p p > urv urv~ |
| 146 | output |
| 147 | launch |
| 148 | }}} |
| 149 | One of the questions will be to ask you if you want to edit the param_card. At this stage, you can enter the path of the previously created param_card.dat. |
| 150 | |
| 151 | * '''Generation of events (Decay done in MG5)''' |
| 152 | Add the decay in MG5 (such that the full spin correlation is taken into-account) and generate events. |
| 153 | In order to learn the syntax of the decay chain. Please type |
| 154 | {{{ |
| 155 | help generate |
| 156 | }}} |
| 157 | Solution in attachments. |