| | 1 | |
| | 2 | |
| | 3 | ==== How do I change the mass of a particle that is assumed massless in MadGraph/MadEvent ? ==== |
| | 4 | |
| | 5 | Several particles (light quarks, charm, he electron and the muon,..) are assumed to be massless in the SM implementation. |
| | 6 | |
| | 7 | ===== '''BEFORE CODE GENERATION''' ===== |
| | 8 | |
| | 9 | One can change this by dowloading the full package and defining a new model. |
| | 10 | |
| | 11 | {{{ |
| | 12 | > cd Models |
| | 13 | > cp sm sm_mc |
| | 14 | > cd sm_mc |
| | 15 | }}} |
| | 16 | |
| | 17 | and then edit the [http://cp3wks05.fynu.ucl.ac.be/cgi-bin/ViewVC/viewvc.cgi/MG_ME/Models/sm/particles.dat?view=markup particles.dat] by changing |
| | 18 | |
| | 19 | {{{ |
| | 20 | c c~ F S ZERO ZERO T c 4. }}} |
| | 21 | |
| | 22 | into |
| | 23 | |
| | 24 | {{{ |
| | 25 | c c~ F S CMASS ZERO T c 4. }}} |
| | 26 | |
| | 27 | You then need to make sure this parameter is read in the lh_readin routine in [http://cp3wks05.fynu.ucl.ac.be/cgi-bin/ViewVC/viewvc.cgi/MG_ME/Models/sm/couplings.f?view=markup couplings.f]. You have to add it also in the [http://cp3wks05.fynu.ucl.ac.be/cgi-bin/ViewVC/viewvc.cgi/MG_ME/Models/sm/param_card.dat?view=markup param_card.dat] if it is not |
| | 28 | already there. Then you're ready to generate processes and events. |
| | 29 | |
| | 30 | |
| | 31 | ===== '''AFTER EVENT GENERATION''' ===== |
| | 32 | |
| | 33 | One can change this by acting directly on the event file with suitable script that change the mass and the four momenta of the interested particle |
| | 34 | in the final state '''by also reshuffling''' all other momenta to conserve overall momentum. Pythia will infact reject events where the energy momentum checkcount is not satisfied. |
| | 35 | |
| | 36 | |
| | 37 | -- Main.FabioMaltoni - 02 Mar 2009 |