Page MenuHomeHEPForge
Files F7879108

No OneTemporary
Actions

View Options

Index: trunk/share/examples/Z-lineshape.sin
===================================================================
--- trunk/share/examples/Z-lineshape.sin (revision 5322)
+++ trunk/share/examples/Z-lineshape.sin (revision 5323)
@@ -1,63 +1,64 @@
# $Id$
# Simple complete physics example: Z lineshape in the Standard Model
########################################################################
#
# Copyright (C) 1999-2014 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
alias lep = e1:E1:e2:E2
alias prt = lep:A
process born = e1, E1 => e2, E2
process rc = e1, E1 => e2, E2, A
compile
cuts = all E >= 100 MeV [prt]
and all abs (cos(Theta)) <= 0.99 [prt]
and all M2 >= (1 GeV)^2 [prt, prt]
# Define title etc. as global variables, that will be used by PLOT:
-$description = "A WHIZARD 2.0 Example"
+$description = "A WHIZARD 2.2 Example"
$x_label = "$\sqrt{s}$/GeV"
$y_label = "$\sigma(s)$/pb"
# Allocate one plot
$title = "The Z Lineshape in $e^+e^-\to\mu^+\mu^-$"
plot lineshape_born { x_min = 88 GeV x_max = 95 GeV }
# Allocate another plot
$title = "The Z Lineshape in $e^+e^-\to\mu^+\mu^-\gamma$"
plot lineshape_rc { x_min = 88 GeV x_max = 95 GeV }
scan sqrts = ((88.0 GeV => 90.0 GeV /+ 0.5 GeV),
(90.1 GeV => 91.9 GeV /+ 0.1 GeV),
(92.0 GeV => 95.0 GeV /+ 0.5 GeV)) {
beams = e1, E1
- integrate (born) { iterations = 1:2000 }
+ integrate (born) { iterations = 2:1000:"gw", 1:2000 }
record lineshape_born (sqrts, integral (born) / 1000)
- integrate (rc) { iterations = 5:5000 }
+ integrate (rc) { iterations = 5:3000:"gw", 2:5000 }
record lineshape_rc (sqrts, integral (rc) / 1000)
}
compile_analysis { $out_file = "Z-lineshape.dat" }
Index: trunk/share/examples/W-endpoint.sin
===================================================================
--- trunk/share/examples/W-endpoint.sin (revision 5322)
+++ trunk/share/examples/W-endpoint.sin (revision 5323)
@@ -1,84 +1,85 @@
# $Id$
# Simple complete physics example: W endpoint in the Standard Model
########################################################################
#
# Copyright (C) 1999-2014 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
# Set up the process
alias parton = u:U:d:D:g
alias jet = parton
### We could have added strange and charm jets:
# alias jet = parton:s:S:c:C
alias lepton = e1:e2
alias neutrino = n1:N1:n2:N2
process enj = parton, parton => lepton, neutrino, jet
compile
ms = 0
mc = 0
me = 0
mmu = 0
# LHC design energy
sqrts = 14 TeV
beams = p, p => pdf_builtin
$pdf_builtin_set = "cteq6l"
# Define reasonable cuts and integrate the cross section(s)
# in order to initialize the phase space grids for similation
cuts = all Pt >= 10 GeV [jet:lepton]
integrate (enj) { iterations = 5:20000:"gw", 3:10000 }
# Define title etc. as global variables, that will be used by PLOT
$description =
"A WHIZARD 2.2 Example."
$y_label = "$N_{\textrm{events}}$"
# Allocate plots
$title = "$W$ Endpoint (``Jacobian Peak'') in $pp\to \ell\bar\nu j$"
$x_label = "$p_T^\ell$/GeV"
histogram pt_lepton (0 GeV, 80 GeV, 2 GeV)
$title = "Jet-$p_T$ in $pp\to \ell\bar\nu j$"
$x_label = "$p_T^j$/GeV"
histogram pt_jet (0 GeV, 80 GeV, 2 GeV)
$title = "Lepton Energy in $pp\to \ell\bar\nu j$"
$x_label = "$E^\ell$/GeV"
histogram e_lepton (0 GeV, 160 GeV, 4 GeV)
$title = "Jet Energy in $pp\to \ell\bar\nu j$"
$x_label = "$E^j$/GeV"
histogram e_jet (0 GeV, 80 GeV, 2 GeV)
analysis = record pt_lepton (eval Pt [extract index 1 [sort by Pt [lepton]]]);
record pt_jet (eval Pt [extract index 1 [sort by Pt [jet]]]);
record e_lepton (eval E [extract index 1 [sort by Pt [lepton]]]);
record e_jet (eval E [extract index 1 [sort by Pt [jet]]])
simulate (enj) { n_events = 1000 }
compile_analysis { $out_file = "W-endpoint.dat" }

File Metadata

Mime Type
text/x-diff
Expires
Tue, Nov 19, 7:36 PM (1 d, 7 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3805854
Default Alt Text
(5 KB)

Event Timeline