whizard is hosted by Hepforge, IPPP Durham

ExampleWendpoint: W-endpoint.sin

File W-endpoint.sin, 2.7 KB (added by Juergen Reuter, 13 years ago)
Line 
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
27model = SM
28
29# Set up the process
30alias parton = u:U:d:D:g
31alias jet = parton
32### We could have added strange and charm jets:
33# alias jet = parton:s:S:c:C
34alias lepton = e1:e2
35alias neutrino = n1:N1:n2:N2
36
37process enj = parton, parton => lepton, neutrino, jet
38compile
39
40ms = 0
41mc = 0
42me = 0
43mmu = 0
44
45# LHC design energy
46sqrts = 14 TeV
47beams = 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
51cuts = all Pt >= 10 GeV [jet:lepton]
52integrate (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"
62histogram 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"
66histogram 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"
70histogram 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"
74histogram e_jet (0 GeV, 80 GeV, 2 GeV)
75
76analysis = 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
81simulate (enj) { n_events = 1000 }
82
83compile_analysis { $out_file = "W-endpoint.dat" }