| 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 |