Commit 97808bd4 authored by Simon Brass's avatar Simon Brass

benchmark processes

parent 4c37a744
model = SM
process ee6q = e1, E1 => bottom, B, u, U, g, g
compile
alias lepton = e1:E1:e2:E2:e3:E3
alias jet = u:U:d:D:s:S:c:C:bottom:B:g
me = 0.0
mW = 80.419
mZ = 91.1882
mH = 130
wW = 2.12
wZ = 2.4952
wH = 0.00429
GF = 1.16639E-5
alphas = 0.0891
mmu = 0.1056583
mtau = 1.777
ms = 0.2
mc = 1.3
mb = 4.8
mtop = 174.3
wtop = 1.6
seed=0
sqrts = 500 GeV
beams = e1, E1
show(model)
# scale = 500 GeV
# ?alpha_s_is_fixed = false
cuts = all E > 10 GeV [lepton]
and all E > 10 GeV [jet]
and all M > 10 GeV [combine [jet, jet]]
and all 5 degree < Theta < 175 degree [lepton, lepton]
and all 5 degree < Theta < 175 degree [jet, lepton]
and all 5 degree < Theta < 175 degree [lepton]
and all 5 degree < Theta < 175 degree [jet]
integrate(ee6q){iterations = 5:10000:"gw", 5:10000}
show(results)
# arXiv:hep-ph/0206070
#
# Table 3, Row 2
include ("/home/users/brass/workspace/sindarin/benchmark/ee6f.sin")
?vis_history = false
process ee3munutaunu = e1, E1 => e2, E2, e2, N2, E3, n3
integrate (ee3munutaunu) { iterations = 5:10000:"gw", 5:10000 }
show(results)
# arXiv:hep-ph/0206070
include ("/home/users/brass/workspace/sindarin/benchmark/ee6f.sin")
process ee4mu2nu = e1, E1 => e2, E2, e2, E2, n2, N2
integrate (ee4mu2nu) { iterations = 5:10000:"gw", 5:10000 }
show(results)
# arXiv:hep-ph/0206070
#
# THIS IS A INCLUDE FILE
model = SM
seed = 12345
me = 0.51099907E-3
mW = 80.419
mZ = 91.1882
mH = 170
wW = 2.12
wZ = 2.4952
wH = 0.3835
GF = 1.16639E-5
alphas = 0.1181
mtop = 174.3
wtop = 1.6
alias quark = d:u:s:c:b:t
alias antiquark = D:U:S:C:B:T
alias lepton = e1:E1:e2:E2:e3:E3
alias neutrinos = n1:N1:n2:N2:n3:N3
cuts = all M > 10 [quark,antiquark]
and all E > 10 [quark:antiquark]
and all E > 10 [lepton:neutrinos]
and all Theta > 5 degree [lepton,lepton]
and all Theta > 5 degree [lepton,quark]
and all Theta > 5 degree [lepton]
and all Theta > 5 degree [quark]
sqrts = 500
beams = e1, E1
?use_vamp_equivalences = false
?vis_history = false
\ No newline at end of file
model = SM
process ee = e1, E1 => e2, E2
sqrts = 250 GeV
?vis_history = false
plot scan_iterations
scan real iteration = ((100 => 100000 /*/ 100)) {
integrate (ee) { iterations = 5:1000:"gw", 5:iteration }
record scan_iterations (iteration, integral (ee), error(ee))
}
write_analysis { $out_file = "scan-iterations.dat" }
model = SM
process ee = e1, E1 => e2, E2
sqrts = 50 GeV
! n_events = 10
! sample_format = lhef
! iterations = 5:100:"gw", 3:10000
! integrate (ee)
plot scan_seed
scan seed = ((1 => 100 /+ 1)) {
integrate (ee) { iterations = 5:1000:"gw", 5:10000 }
record scan_seed (seed, integral (ee), error(ee))
}
compile_analysis { $out_file = "scan-seed.dat" }
model = SM
process ee = e1, E1 => e2, E2
sqrts = 50 GeV
! n_events = 10
! sample_format = lhef
! iterations = 5:100:"gw", 3:10000
! integrate (ee)
plot scan_seed
scan seed = ((1 => 100 /+ 1)) {
integrate (ee) { iterations = 2:1000:"gw", 1:5000 }
record scan_seed (seed, integral (ee), error(ee))
}
compile_analysis { $out_file = "scan-seed.dat" }
model = SM
process ee = e1, E1 => e2, E2
seed = 12345
compile
plot scan_sqrts
?vis_history = false
scan sqrts = ((75.0 GeV => 100.0 GeV /+/ 250)) {
integrate (ee) { iterations = 2:1000:"gw", 1:5000 }
record scan_sqrts (sqrts, integral (ee) / 1000, error(ee) / 1000)
}
write_analysis { $out_file = "scan-sqrts.dat" }
# $Id: W-endpoint.sin 7444 2016-02-17 15:37:20Z jr_reuter $
# Simple complete physics example: W endpoint in the Standard Model
########################################################################
#
# Copyright (C) 1999-2016 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. Cuts $p_T($visible$)\ge10 \textrm{ GeV}$."
$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 = 2000 }
compile_analysis { $out_file = "W-endpoint.dat" }
model = SM
process ggqqW = gg, gg => quark, antiquark, e1, N1
process ggqqWg = gg, gg => quark, antiquark, e1, N1, g
process ggqqWgg = gg, gg => quark, antiquark, e1, N1, g, g
process ggqqWggg = gg, gg => quark, antiquark, e1, N1, g, g, g
alias quark = d:u:s:c
alias antiquark = D:U:S:C
# Aliases for cuts' definition
alias lepton = e1:N1
alias jet = u:U:d:D:s:S:c:D:g
me = 0.0
mW = 80.419
mZ = 91.1882
mH = 130
wW = 2.049
wZ = 2.443
GF = 1.16639E-5
alphas = 0.1178
ms = 0
mc = 0
sqrts = 14 TeV
beams = p, p => lhapdf
$lhapdf_file = "CTEQ6L"
cuts = all Pt > 20 GeV [lepton, jet]
and all -2.5 < Eta < 2.5 [lepton, jet]
and all Dist > 0.7 [jet, jet]
and all Dist > 0.4 [lepton, jet]
integrate (ggqqW) {iterations = 5:10000:"gw", 5:10000}
\ No newline at end of file
model = SM
process ggqqW = gg, gg => quark, antiquark, e1, N1
process ggqqWg = gg, gg => quark, antiquark, e1, N1, g
process ggqqWgg = gg, gg => quark, antiquark, e1, N1, g, g
process ggqqWggg = gg, gg => quark, antiquark, e1, N1, g, g, g
alias quark = d:u:s:c
alias antiquark = D:U:S:C
# Aliases for cuts' definition
alias lepton = e1:N1
alias jet = u:U:d:D:s:S:c:D:g
me = 0.0
mW = 80.419
mZ = 91.1882
mH = 130
wW = 2.049
wZ = 2.443
GF = 1.16639E-5
alphas = 0.1178
mu = md = ms = mc = 0
sqrts = 14 TeV
beams = p, p => lhapdf
$lhapdf_file = "CTEQ6L"
cuts = all Pt > 20 GeV [lepton, jet]
and all -2.5 < Eta < 2.5 [lepton, jet]
and all Dist > 0.7 [jet, jet]
and all Dist > 0.4 [lepton, jet]
integrate (ggqqWg) {iterations = 5:1000000:"gw", 5:1000000}
\ No newline at end of file
model = SM
process ggqqW = gg, gg => quark, antiquark, e1, N1
process ggqqWg = gg, gg => quark, antiquark, e1, N1, g
process ggqqWgg = gg, gg => quark, antiquark, e1, N1, g, g
process ggqqWggg = gg, gg => quark, antiquark, e1, N1, g, g, g
alias quark = d:u:s:c
alias antiquark = D:U:S:C
# Aliases for cuts' definition
alias lepton = e1:N1
alias jet = u:U:d:D:s:S:c:D:g
me = 0.0
mW = 80.419
mZ = 91.1882
mH = 130
wW = 2.049
wZ = 2.443
GF = 1.16639E-5
alphas = 0.1178
mu = md = ms = mc = 0
sqrts = 14 TeV
beams = p, p => lhapdf
$lhapdf_file = "CTEQ6L"
cuts = all Pt > 20 GeV [lepton, jet]
and all -2.5 < Eta < 2.5 [lepton, jet]
and all Dist > 0.7 [jet, jet]
and all Dist > 0.4 [lepton, jet]
integrate (ggqqWgg) {iterations = 5:1250000:"gw", 5:1000000}
\ No newline at end of file
model = SM
process ggqqW = gg, gg => quark, antiquark, e1, N1
process ggqqWg = gg, gg => quark, antiquark, e1, N1, g
process ggqqWgg = gg, gg => quark, antiquark, e1, N1, g, g
process ggqqWggg = gg, gg => quark, antiquark, e1, N1, g, g, g
alias quark = d:u:s:c
alias antiquark = D:U:S:C
# Aliases for cuts' definition
alias lepton = e1:N1
alias jet = u:U:d:D:s:S:c:D:g
me = 0.0
mW = 80.419
mZ = 91.1882
mH = 130
wW = 2.049
wZ = 2.443
GF = 1.16639E-5
alphas = 0.1178
mu = md = ms = mc = 0
sqrts = 14 TeV
beams = p, p => lhapdf
$lhapdf_file = "CTEQ6L"
cuts = all Pt > 20 GeV [lepton, jet]
and all -2.5 < Eta < 2.5 [lepton, jet]
and all Dist > 0.7 [jet, jet]
and all Dist > 0.4 [lepton, jet]
integrate (ggqqWggg) {iterations = 5:2000000:"gw", 5:1000000}
\ No newline at end of file
model = SM
process h_decay_WWstar = H => Wp, e1, N1
sqrts = 500 GeV
me = 0
mH = 125.36 GeV
plot integral_error
?vis_history = false
integrate (h_decay_WWstar) { iterations = 5:5000:"gw", 5:10000 }
record integral_error (sqrts, integral(h_decay_WWstar), error(h_decay_WWstar))
write_analysis { $out_file = "result.dat" }
show (results)
\ No newline at end of file
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment