1 | # $Id: W-endpoint.sin 3132 2011-04-07 15:12:02Z jr_reuter $ |
---|
2 | # Simple complete physics example: W endpoint in the Standard Model |
---|
3 | ######################################################################## |
---|
4 | # |
---|
5 | # Copyright (C) 1999-2011 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 | # Christian Speckner <christian.speckner@physik.uni-freiburg.de> |
---|
10 | # |
---|
11 | # WHIZARD is free software; you can redistribute it and/or modify it |
---|
12 | # under the terms of the GNU General Public License as published by |
---|
13 | # the Free Software Foundation; either version 2, or (at your option) |
---|
14 | # any later version. |
---|
15 | # |
---|
16 | # WHIZARD is distributed in the hope that it will be useful, but |
---|
17 | # WITHOUT ANY WARRANTY; without even the implied warranty of |
---|
18 | # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
---|
19 | # GNU General Public License for more details. |
---|
20 | # |
---|
21 | # You should have received a copy of the GNU General Public License |
---|
22 | # along with this program; if not, write to the Free Software |
---|
23 | # Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. |
---|
24 | # |
---|
25 | ######################################################################## |
---|
26 | |
---|
27 | model = SM |
---|
28 | |
---|
29 | # Set up the process |
---|
30 | alias parton = u:U:d:D:g |
---|
31 | alias jet = parton |
---|
32 | ### We could have added strange and charm jets: |
---|
33 | # alias jet = parton:s:S:c:C |
---|
34 | alias lepton = e1:e2 |
---|
35 | alias neutrino = n1:N1:n2:N2 |
---|
36 | |
---|
37 | process enj = parton, parton => lepton, neutrino, jet |
---|
38 | compile |
---|
39 | |
---|
40 | ms = 0 |
---|
41 | mc = 0 |
---|
42 | me = 0 |
---|
43 | mmu = 0 |
---|
44 | |
---|
45 | # LHC design energy |
---|
46 | sqrts = 14 TeV |
---|
47 | beams = p, p => lhapdf { $lhapdf_file = "cteq5l.LHgrid" } |
---|
48 | |
---|
49 | # Define reasonable cuts and integrate the cross section(s) |
---|
50 | # in order to initialize the phase space grids for similation |
---|
51 | cuts = all Pt >= 10 GeV [jet:lepton] |
---|
52 | integrate (enj) { iterations = 5:20000 } |
---|
53 | |
---|
54 | # Define title etc. as global variables, that will be used by PLOT |
---|
55 | $description = |
---|
56 | "A WHIZARD 2.0 Example." |
---|
57 | $y_label = "$N_{\textrm{events}}$" |
---|
58 | |
---|
59 | # Allocate plots |
---|
60 | $title = "$W$ Endpoint (``Jacobian Peak'') in $pp\to \ell\bar\nu j$" |
---|
61 | $x_label = "$p_T^\ell$/GeV" |
---|
62 | histogram pt_lepton (0 GeV, 80 GeV, 2 GeV) |
---|
63 | |
---|
64 | $title = "Jet-$p_T$ in $pp\to \ell\bar\nu j$" |
---|
65 | $x_label = "$p_T^j$/GeV" |
---|
66 | histogram pt_jet (0 GeV, 80 GeV, 2 GeV) |
---|
67 | |
---|
68 | $title = "Lepton Energy in $pp\to \ell\bar\nu j$" |
---|
69 | $x_label = "$E^\ell$/GeV" |
---|
70 | histogram e_lepton (0 GeV, 160 GeV, 4 GeV) |
---|
71 | |
---|
72 | $title = "Jet Energy in $pp\to \ell\bar\nu j$" |
---|
73 | $x_label = "$E^j$/GeV" |
---|
74 | histogram e_jet (0 GeV, 80 GeV, 2 GeV) |
---|
75 | |
---|
76 | analysis = record pt_lepton (eval Pt [extract index 1 [sort by Pt [lepton]]]); |
---|
77 | record pt_jet (eval Pt [extract index 1 [sort by Pt [jet]]]); |
---|
78 | record e_lepton (eval E [extract index 1 [sort by Pt [lepton]]]); |
---|
79 | record e_jet (eval E [extract index 1 [sort by Pt [jet]]]) |
---|
80 | |
---|
81 | simulate (enj) { n_events = 1000 } |
---|
82 | |
---|
83 | compile_analysis { $out_file = "W-endpoint.dat" } |
---|