# $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_CKM

# Alias definition
include ( "/group/ilc/soft/gcc481/whizard2/share/include/aliases.sin" )

#
#      e+e- --> mu+mu- and tau+tau-
#
# process P2f_z_l = e1, E1 => ( e2, E2 ) + (e3, E3 ) { process_num_id = 250106 }
polarized e3, E3, "W+", "W-", Z

process P2f_tautau = e1, E1 => (e3, E3 ) { process_num_id = 250106 }

compile

#
logical ?is_higgs_process = false

include ("/home/ilc/miyamoto/work/160106-tauola/checkrun/sindarin/pythia6-parameters.sin" )

sqrts = 250 GeV

beams = e1, E1 => beam_events => isr
$beam_events_file = "/group/ilc/soft/gcc481/GP/beamdata/lumiee250/lumiee250-run1.dat"
beams_pol_density = @(-1), @(+1)

integrate ( P2f_tautau ) { ?polarized_events = true }

printf "#Total cross section : P2f_tautau %13e fb" ( integral( P2f_tautau ) )

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

n_events = 100
simulate ( P2f_tautau ) { ?polarized_events = true }