# $Id: LEP_higgs.sin 2293 2010-04-11 23:57:50Z jr_reuter $
# Simple complete physics example: Higgs search at LEP
########################################################################
#
# Copyright (C) 1999-2015 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
# model = SM_Higgs
model = SM_CKM

# The Higgsstrahlung process
process e2e2h = e1, E1 => e2, E2, h
polarized e1, E1, e2, E2, e3, E3, "W+", "W-", Z

compile

# Set the input parameters
mH = 125 GeV
wH = 4.3 MeV

# ##################################################################################
# Parameters for Pythia hadronization
# ?ps_isr_active = true 
# Do we need to set ps_isr_active = true, when isr is set by beams command

?ps_fsr_active = true
?hadronization_active = true
$shower_method = "PYTHIA6"

# h--> tau+ tau-
$ps_PYTHIA_PYGIVE = "MSTJ(28)=2;PMAS(25,1)=125.0;PMAS(25,2)=0.0043;MSTJ(41)=2; MSTU(22)=20; PARJ(21)=0.40000; PARJ(41)=0.11000; PARJ(42)=0.52000; PARJ(81)=0.25000;PARJ(82)=1.90000; MSTJ(11)=3; PARJ(54)=-0.03100; PARJ(55)=-0.00200;PARJ(1)=0.08500; PARJ(3)=0.45000; PARJ(4)=0.02500; PARJ(2)=0.31000;PARJ(11)=0.60000; PARJ(12)=0.40000; PARJ(13)=0.72000; PARJ(14)=0.43000;PARJ(15)=0.08000; PARJ(16)=0.08000; PARJ(17)=0.17000; MSTP(3)=1;MWID(25)=2;BRAT(212)=0.00044;BRAT(213)=0.0268;BRAT(214)=0.578;BRAT(219)=0.000221;BRAT(220)=0.0637;BRAT(222)=0.0856;BRAT(223)=0.0023;BRAT(224)=0.00155;BRAT(225)=0.0267;BRAT(226)=0.216;MDME(210,1)=0;MDME(211,1)=0;MDME(212,1)=0;MDME(213,1)=0;MDME(214,1)=0;MDME(215,1)=0;MDME(218,1)=0;MDME(219,1)=0;MDME(220,1)=1;MDME(222,1)=0;MDME(223,1)=0;MDME(224,1)=0;MDME(225,1)=0;MDME(226,1)=0; MDME(174,1)=0;MDME(175,1)=0;MDME(176,1)=0;MDME(177,1)=0;MDME(178,1)=0;MDME(179,1)=0;MDME(182,1)=0;MDME(183,1)=0;MDME(184,1)=0;MDME(185,1)=0;MDME(186,1)=1;MDME(187,1)=0; MDME(190,1)=0;MDME(191,1)=0;MDME(192,1)=0;MDME(194,1)=0;MDME(195,1)=0;MDME(196,1)=0;MDME(198,1)=0;MDME(199,1)=0;MDME(200,1)=0;MDME(206,1)=0;MDME(207,1)=0;MDME(208,1)=1;"
# MSTJ(28)=2 to call external Tau decay routine , MSTJ(28)=0 is not to call
# ! MDME(219,1)=0;  To suppress h-> mu mu decay
# !!!Event_generation_branching_ratio=0.578
# ! Hdecay mode
# !   dd             uu             ss              cc       bb              tt           e1e1
# !  MDME(210,1)=0;MDME(211,1)=0;MDME(212,1)=0;MDME(213,1)=0;MDME(214,1)=0;MDME(215,1)=0;MEME(218,1)=0;
# !  e2e2           e3e3            gg          gammagamma    gammaZ0         ZZ           WW
# !  MDME(219,1)=0;MDME(220,1)=0;MDME(222,1)=0;MDME(223,1)=0;MDME(224,1)=0;MDME(225,1)=1;MDME(226,1)=0;
# ! Z decay mode
# !  dd            uu             ss             cc            bb            tt           e1e1
# !  MDME(174,1)=0;MDME(175,1)=0;MDME(176,1)=0;MDME(177,1)=0;MDME(178,1)=0;MDME(179,1)=0;MDME(182,1)=0;
# !  n1n1          e2e2           n2n2           e3e3         n3n3
# !  MDME(183,1)=1;MDME(184,1)=0;MDME(185,1)=1;MDME(186,1)=0;MDME(187,1)=1;
# !
# ! W decay mode
# !  dbar-u         dbar-c         dbar-t        sbar-u       sbar-c         sbar-t
# !  MDME(190,1)=0;MDME(191,1)=0;MDME(192,1)=0;MDME(194,1)=0;MDME(195,1)=0;MDME(196,1)=0;
# !   bbar-u        bbar-c         bbar-t       e+-nu_e       mu+-nu_mu      tau+-nu_tau
# !  MDME(198,1)=0;MDME(199,1)=0;MDME(200,1)=0;MDME(206,1)=0;MDME(207,1)=0;MDME(208,1)=0;
# !
# ?ps_PYTHIA_verbose = true

sqrts = 250 GeV

beams = e1, E1 => beam_events => isr
$beam_events_file = "%%BEAM_EVENTS_FILE%%"

beams_pol_density = @(-1), @(+1)

integrate( e2e2h ) {  
  iterations = 1:100000, 10:100000 
  ?polarized_events = true  
}

show( results )
#

# sample_format = lcio
sample_format = stdhep_ev4
$sample = "E250.Pe2e2h_h2tautau.Gwhizard2-beam.eL.pR"
?write_raw = false

n_events = %%N_EVENTS%%
simulate ( e2e2h ) {  ?polarized_events = true  }