Commit ec8eb348 authored by Simon Brass's avatar Simon Brass

remove old sindarin files from benchmark directory

parent 55404061
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 3
include ("/home/users/brass/workspace/sindarin/benchmark/ee6f.sin")
?vis_history = false
process ee3munutaunu = e1, E1 => e2, E2, n2, N2, e3, E3
integrate (ee3munutaunu) { iterations = 5:10000:"gw", 5:10000 }
show(results)
# arXiv:hep-ph/0206070
#
# Table 3, Row 4
include ("/home/users/brass/workspace/sindarin/benchmark/ee6f.sin")
?vis_history = false
process ee3mu2nue = e1, E1 => e1, N1, e2, E2, n2, E2
integrate (ee3mu2nue) { 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/whizard-sindarin/benchmark/ee6f.sin")
process ee4mu2nu = e1, E1 => e2, E2, e2, E2, n2, N2
integrate (ee4mu2nu) { iterations = 10:100000:"gw", 10:100000 }
show(results)
model = SM
process ee = e1, E1 => e2, E2
seed = 12345
sqrts = 150 GeV
?vis_history = false
integrate (ee) { iterations = 5:{{ N_CALL }}:"gw", 5:{{ N_CALL }} }
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
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" }
# ---
# file benchmark/eemumu/task.yml
- hosts: notebook cluster
tasks:
- name: deploy e1, E1 => e2, E2
template:
src: eemumu.sin.j2
dest: "{{ SIN_LOCAL }}/eemumu-{{ item.name }}.sin"
with_items:
- { name: 10e5, n_calls: 10000 }
- { name: 10e6, n_calls: 100000 }
- { name: 10e7, n_calls: 1000000 }
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
seed = 82618948
?vis_history = false
mZ = 91.1876
mW = 80.385
mb = 4.2
mH = 1.250000000000E+02
me = 0.0
ms = 0
mc = 0
mmu = 0
mtop = 172.0
GF = 1.1663787E-005
wZ = 2.443
wW = 2.049000000000E+00
wH = 4.143000000000E-03
?alpha_s_is_fixed = false
?alpha_s_from_mz = true
?alpha_s_from_lambda_qcd = false
alphas = 0.118 ! (Z pole, NLL running to mu_h)
alpha_s_order = 2
alpha_s_nf = 5
scale = mtop
$method = "threshold"
$born_me_method = "threshold"
$loop_me_method = "threshold"
$correlation_me_method = "threshold"
$real_tree_me_method = "threshold"
$soft_mismatch_me_method = "threshold"
?openloops_use_cms = true
alpha_power = 4
alphas_power = 0
sqrts = 500 GeV
sample_format = hepmc
?write_raw = false
?unweighted = false
?hepmc_output_cross_section = true
?negative_weights = true
$method = "openloops"
$born_me_method = "openloops"
$loop_me_method = "openloops"
$correlation_me_method = "openloops"
$real_tree_me_method = "openloops"
$soft_mismatch_me_method = "openloops"
mult_call_real = 20
mult_call_virt = 0.75
?openloops_use_cms = true
?negative_weights = true
sqrts = 350
library = nlo_WbWb_Real
real scale_var = 0.0
scale_var = 0.125
sqrts = 350
scale = scale_var * mtop
! wtop = 0
! alpha_power = 1
! iterations = 5:400:"gw",5:800
! process top_width_nlo = t => Wp, b { nlo_calculation = "Full" } #KEEP
! integrate (top_width_nlo)
! wtop = integral(top_width_nlo)
wtop = 1.5
alpha_power = 4
iterations = 10:40000:"gw",5:20000
process nlo_WbWb_Real = e1, E1 => Wp, Wm, b, B { nlo_calculation = "Real" }
integrate (nlo_WbWb_Real)
printf "RES %8e %8e %8e" (scale_var, integral(nlo_WbWb_Real), error(nlo_WbWb_Real))
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