| 1 | |
| 2 | |
| 3 | == Quarkonium MadGraph == |
| 4 | |
| 5 | |
| 6 | |
| 7 | === People === |
| 8 | 1. Pierre Artoisenet |
| 9 | 1. Fabio Maltoni |
| 10 | 1. Tim Stelzer |
| 11 | |
| 12 | === Aim of the project === |
| 13 | Build an event generator for quarkonium physics. |
| 14 | |
| 15 | === Tools === |
| 16 | * Madgraph/MadEvent v4 and StandAlone |
| 17 | |
| 18 | === Work === |
| 19 | We have identified two phases: |
| 20 | |
| 21 | 1. Add to MadGraph the possibility of generating pairs of heavy quarks (same or different flavors) with given quantum numbers, of spin, total angular momentum, color: %$^{2S+1}L_J[1,8]$%, for %$L=0,1$%. |
| 22 | 1. Interface it to MadEvent for event generation. |
| 23 | |
| 24 | ==== Physics results and papers writeup ==== |
| 25 | |
| 26 | The paper phase 1 is out! |
| 27 | [http://arxiv.org/abs/0712.2770 MadOnia paper] |
| 28 | |
| 29 | |
| 30 | === Preliminary code === |
| 31 | The code phase 1 is done. The quarkonium amplitude generator has been interfaced with a basic phase space generator to produce cross sections and distributions. |
| 32 | The code is available under request. |
| 33 | In order to generate heavy quarkonium cross sections, you must: |
| 34 | |
| 35 | '''1. Copy the directory Template into my_proc''' |
| 36 | |
| 37 | '''2. Edit the proc_card.dat in my_proc/Card/proc_card.dat.''' |
| 38 | |
| 39 | You have to enter the parton-level process |
| 40 | |
| 41 | IMPORTANT: THE HEAVY QQ bar PAIR MUST BE PLACED AT THE END OF THE LINE DEFINING THE PROCESS !! |
| 42 | |
| 43 | The quantum numbers of the QQ bar pair are written in squared brackets using the spectroscopic notation |
| 44 | |
| 45 | [(2S+1)LJC] with S=spin, L=orbital momentum, J= total momentum, C=color |
| 46 | |
| 47 | For example, type gg>gcc~[3S11] if you want to consider gg > g + charmonium with S=1, L=0 J=1, c=1 |
| 48 | (J/psi production via a color-singlet transition) |
| 49 | |
| 50 | see also an example of proc_card.dat in the Template. |
| 51 | |
| 52 | '''3. go into my_proc/bin and type "./newprocess"''' |
| 53 | |
| 54 | '''4. Fill in the two input files my_proc/Cards/param_card.dat and my_proc/Cards/run_card.dat''' |
| 55 | |
| 56 | '''5. go to my_proc/SubProcesses/P???''' |
| 57 | |
| 58 | P??? is a specific process that you have requested in the proc_card.dat. |
| 59 | |
| 60 | '''6. type "make" and then ./cross_section''' |
| 61 | |
| 62 | the PT and rapidity distributions of the quarkonium state are written in the file histogram.dat |
| 63 | |
| 64 | |
| 65 | Remarks: |
| 66 | |
| 67 | * We use the BBL normalization for the long distance matrix elements. |
| 68 | |
| 69 | |
| 70 | === MadOnia version II === |
| 71 | |
| 72 | |
| 73 | 11/2008: Phase 2 is done! MG_ME is been upgraded with the MadOnia implementation. Production of quarkonium unweighted events has been successfully tested. |
| 74 | |
| 75 | |
| 76 | |
| 77 | -- Main.PierreArtoisenet - 25 Nov 2008 |
| 78 | |
| 79 | |
| 80 | |
| 81 | |
| 82 | |
| 83 | |