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 | |
---|
28 | model = SM |
---|
29 | |
---|
30 | # This is the CC10 charged current process important |
---|
31 | # for LEP physics. No flavor summation for jets here. |
---|
32 | process cc10 = e1, E1 => e2, N2, u, D |
---|
33 | |
---|
34 | compile |
---|
35 | |
---|
36 | # We set the muon mass to zero. |
---|
37 | mmu = 0 |
---|
38 | |
---|
39 | sqrts = 209 GeV |
---|
40 | |
---|
41 | integrate (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 | |
---|
46 | luminosity = 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" |
---|
57 | histogram 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" |
---|
61 | histogram e_muon (0 GeV, 209 GeV, 4) |
---|
62 | |
---|
63 | analysis = record m_jets (eval M [u,D]); |
---|
64 | record e_muon (eval E [e2]) |
---|
65 | |
---|
66 | simulate (cc10) |
---|
67 | |
---|
68 | compile_analysis { $out_file = "cc10.dat" } |
---|