| 1 | # $Id: LEP_higgs.sin 2293 2010-04-11 23:57:50Z jr_reuter $ |
|---|
| 2 | # Simple complete physics example: Higgs search at LEP |
|---|
| 3 | ######################################################################## |
|---|
| 4 | # |
|---|
| 5 | # Copyright (C) 1999-2014 by |
|---|
| 6 | # Wolfgang Kilian <kilian@physik.uni-siegen.de> |
|---|
| 7 | # Thorsten Ohl <ohl@physik.uni-wuerzburg.de> |
|---|
| 8 | # Juergen Reuter <juergen.reuter@desy.de> |
|---|
| 9 | # with contributions from |
|---|
| 10 | # Christian Speckner <cnspeckn@googlemail.com> |
|---|
| 11 | # |
|---|
| 12 | # WHIZARD is free software; you can redistribute it and/or modify it |
|---|
| 13 | # under the terms of the GNU General Public License as published by |
|---|
| 14 | # the Free Software Foundation; either version 2, or (at your option) |
|---|
| 15 | # any later version. |
|---|
| 16 | # |
|---|
| 17 | # WHIZARD is distributed in the hope that it will be useful, but |
|---|
| 18 | # WITHOUT ANY WARRANTY; without even the implied warranty of |
|---|
| 19 | # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
|---|
| 20 | # GNU General Public License for more details. |
|---|
| 21 | # |
|---|
| 22 | # You should have received a copy of the GNU General Public License |
|---|
| 23 | # along with this program; if not, write to the Free Software |
|---|
| 24 | # Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. |
|---|
| 25 | # |
|---|
| 26 | ######################################################################## |
|---|
| 27 | |
|---|
| 28 | model = NMSSM |
|---|
| 29 | |
|---|
| 30 | alias n = n1:n2:n3:N1:N2:N3 |
|---|
| 31 | alias q = u:d:s:c |
|---|
| 32 | alias Q = U:D:S:C |
|---|
| 33 | alias l = e1:e2:E1:E2 |
|---|
| 34 | |
|---|
| 35 | # The missing-energy channel |
|---|
| 36 | process zh = e1, E1 => Z,h1 |
|---|
| 37 | |
|---|
| 38 | ?vis_history = false |
|---|
| 39 | |
|---|
| 40 | sqrts = 250 GeV |
|---|
| 41 | |
|---|
| 42 | # Set the input parameters |
|---|
| 43 | # Running b mass |
|---|
| 44 | mb = 0 |
|---|
| 45 | ms = 0 |
|---|
| 46 | mc = 0 |
|---|
| 47 | mneu1 = 15 |
|---|
| 48 | mneu2 = 85 |
|---|
| 49 | # Lepton mass |
|---|
| 50 | me = 0 |
|---|
| 51 | mmu = 0 |
|---|
| 52 | mtau = 0 |
|---|
| 53 | # higgs parameters |
|---|
| 54 | mh0_2 = 60 |
|---|
| 55 | |
|---|
| 56 | |
|---|
| 57 | process zdecl = Z => l,l |
|---|
| 58 | process hdec = h1 => neu1,neu2 |
|---|
| 59 | process k2dec = neu2 => h2,neu1 |
|---|
| 60 | process h2dec = h2 => b,B |
|---|
| 61 | |
|---|
| 62 | ?ps_isr_active = true |
|---|
| 63 | ?ps_fsr_active = true |
|---|
| 64 | ?ps_use_PYTHIA_shower = true |
|---|
| 65 | ?hadronization_active = true |
|---|
| 66 | #?ps_PYTHIA_verbose = true |
|---|
| 67 | |
|---|
| 68 | ps_max_n_flavors = 5 |
|---|
| 69 | ps_mass_cutoff = 1 |
|---|
| 70 | |
|---|
| 71 | $ps_PYTHIA_PYGIVE = "MSTP(125)=2; MSTU(21)=1; IMSS(1)=11; IMSS(13)=1" |
|---|
| 72 | |
|---|
| 73 | real ycut = mlm_ETclusminE**2 |
|---|
| 74 | cuts = all kT > ycut [j, j] |
|---|
| 75 | cuts = all M >= 10 GeV [q,Q] |
|---|
| 76 | |
|---|
| 77 | integrate(zh,zdecl) |
|---|
| 78 | ?slha_read_decays = true |
|---|
| 79 | read_slha("nmssm.slha") |
|---|
| 80 | mneu1 = 15 |
|---|
| 81 | mneu2 = 85 |
|---|
| 82 | mh0_2 = 60 |
|---|
| 83 | integrate(hdec,k2dec,h2dec) |
|---|
| 84 | unstable Z(zdecl) |
|---|
| 85 | unstable h1(hdec) |
|---|
| 86 | unstable neu2(k2dec) |
|---|
| 87 | unstable h2(h2dec) |
|---|
| 88 | |
|---|
| 89 | # luminosity = 10 |
|---|
| 90 | # ?write_raw = false |
|---|
| 91 | simulate (zh) {n_events = 10 sample_format = hepmc } |
|---|
| 92 | |
|---|
| 93 | show ( results ) |
|---|
| 94 | compile |
|---|
