whizard is hosted by Hepforge, IPPP Durham

ExampleCC10: LEP_cc10.sin

File LEP_cc10.sin, 2.2 KB (added by Juergen Reuter, 10 years ago)

SINDARIN script for this example.

Line 
1# $Id: LEP_cc10.sin 2227 2010-04-05 04:46:08Z reuter $
2# Simple complete physics example: W endpoint in the Standard Model
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
28model = SM
29
30# This is the CC10 charged current process important
31# for LEP physics. No flavor summation for jets here.
32process cc10 = e1, E1 => e2, N2, u, D
33
34compile
35
36# We set the muon mass to zero.
37mmu = 0
38
39sqrts = 209 GeV
40
41integrate (cc10) { iterations = 12:20000 }
42
43# We set a luminosity of 10 inverse femtobarn. Note that this is roughly
44# two orders of magnitude higher than the final LEP2 integrated luminosity.
45
46luminosity = 10
47
48# Define title etc. as global variables, that will be used by PLOT
49$description =
50  "A WHIZARD 2.2 Example.
51   Charged current CC10 process from LEP 2."
52$y_label = "$N_{\textrm{events}}$"
53
54# Allocate plots
55$title = "Di-jet invariant mass $M_{jj}$ in $e^+e^- \to \mu^- \bar\nu_\mu u \bar d$"
56$x_label = "$M_{jj}$/GeV"
57histogram m_jets (70 GeV, 90 GeV, 0.5 GeV)
58
59$title = "Muon energy $E_\mu$ in $e^+e^- \to \mu^- \bar\nu_\mu u \bar d$"
60$x_label = "$E_\mu$/GeV"
61histogram e_muon (0 GeV, 209 GeV, 4)
62
63analysis = record m_jets (eval M [u,D]);
64           record e_muon (eval E [e2])
65
66simulate (cc10)
67
68compile_analysis { $out_file = "cc10.dat" }