Ticket #769: NMSSM_zh_h0k12.sin

# $Id: LEP_higgs.sin 2293 2010-04-11 23:57:50Z jr_reuter $
# Simple complete physics example: Higgs search at LEP
########################################################################
#
# Copyright (C) 1999-2014 by
#     Wolfgang Kilian <kilian@physik.uni-siegen.de>
#     Thorsten Ohl <ohl@physik.uni-wuerzburg.de>
#     Juergen Reuter <juergen.reuter@desy.de>
#     with contributions from
#     Christian Speckner <cnspeckn@googlemail.com>
#
# WHIZARD is free software; you can redistribute it and/or modify it
# under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2, or (at your option)
# any later version.
#
# WHIZARD is distributed in the hope that it will be useful, but
# WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
#
########################################################################

model = SM_CKM

alias n   = n1:n2:n3:N1:N2:N3
alias q   = u:d:s:c
alias Q   = U:D:S:C
alias l  = e1:e2:E1:E2

# The missing-energy channel
process zh = e1, E1 => l,l,H

?vis_history = true

sqrts = 250 GeV

# Set the input parameters
# Running b mass
mb = 2.9 GeV
ms = 0
mc = 0
#mneu1 = 15
#mneu2 = 85
# Lepton mass
me = 0
mmu = 0
mtau = 0
# higgs parameters
#mh0_2 = 60

#process hdec  = h1 => neu1,neu2
#process k2dec = neu2 => h2,neu1
#process h2dec = h2 => b,B
process hdec = H => b,B

?ps_fsr_active = true
#?ps_isr_active = false
?hadronization_active = true
$shower_method = "PYTHIA6"

ps_max_n_flavors = 5
ps_mass_cutoff = 1

#$ps_PYTHIA_PYGIVE = "MSTP(125)=2; MSTU(21)=1"

real ycut = mlm_ETclusminE**2
cuts = all kT > ycut [j, j]
cuts = all M >= 10 GeV [q,Q]

integrate(zh) {beams = e1, E1 => isr isr_mass = 0.000511}

#?slha_read_decays = true
#read_slha("nmssm.slha")
#mneu1 = 15
#mneu2 = 85
#mh0_2 = 60

integrate(hdec)
unstable H(hdec)
#integrate(k2dec)
#unstable neu2(k2dec)
#integrate(h2dec)
#unstable h2(h2dec)

luminosity = 1 
?write_raw = false
?rebuild_events = true
#simulate (zh) {n_events = 100 sample_format = stdhep $extension_stdhep = "stdhep"}
#simulate (zh) {n_events = 100 sample_format = hepevt}
simulate (zh) {n_events = 100 sample_format = athena}

show ( results )