Commit dc28166b authored by Simon Brass's avatar Simon Brass

[benchmark] reuse process file from WK

parent 2ed13e82
......@@ -7,41 +7,56 @@ model = SM
seed = 12345
alias quark = d:u:s:c
alias antiquark = D:U:S:C
# Aliases for cuts' definition
alias q = d:u:s:c
alias Q = D:U:S:C
alias l = e1:e2:e3
alias L = E1:E2:E3
alias n = n1:n2:n3
alias N = N1:N2:N3
alias lepton = l:L
alias jet = quark:antiquark:gl
alias j = q:Q:g
process g2 = g, g => e1, N1, q, Q
process g3 = g, g => e1, N1, q, Q, g
process g4 = g, g => e1, N1, q, Q, g, g
process g5 = g, g => e1, N1, q, Q, g, g, g
process j1 = j, j => e1, N1, j
process j2 = j, j => e1, N1, j, j
process j3 = j, j => e1, N1, j, j, j
process j4 = j, j => e1, N1, j, j, j, j
process b0 = j, j => l, L, n, N, b, B
process b1 = j, j => l, L, n, N, b, B, j
process b2 = j, j => l, L, n, N, b, B, j, j
process bb = j, j => l, L, n, N, b, B, b, B
compile
me = 0
mmu = 0
mtau = 0
mW = 80.419
mZ = 91.1882
mH = 120
wW = 2.049
wZ = 2.443
wH = 0.0036 GeV
GF = 1.16639E-5 / (1 GeV)^2
alphas = 0.1178
ms = 0
mc = 0
mb = 4.2
mb = 4.2 GeV
mtop = 174.0 GeV
wtop = 1.523 GeV
GF = 1.16639e-5 / (1 GeV)^2
mZ = 91.1882 GeV
mW = 80.419 GeV
mH = 120.0 GeV
alphas = 0.1178
wZ = 2.443 GeV
wW = 2.049 GeV
wH = 0.0036 GeV
cuts =
all Pt > 20 GeV [j:lepton] and
all -2.5 < Eta < 2.5 [j:lepton] and
all Dist > 0.7 [j,j] and
all Dist > 0.4 [j,lepton]
sqrts = 14 TeV
beams = p, p => pdf_builtin
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]
?vis_history = false
\ No newline at end of file
......@@ -28,14 +28,14 @@ logical ?is_higgs_process = true
sqrts = 500 GeV
beams = e1, E1 => isr
beams_pol_density = @(-1), @(+1)
! beams = e1, E1 => isr
! beams_pol_density = @(-1), @(+1)
?isr_recoil = true
real default_M_cut = 10 GeV
cuts = all M > default_M_cut [ quarks, quarks ]
! ?isr_recoil = false
real default_M_cut = 50 GeV ! 10 GeV
cuts = all M > default_M_cut [ quarks, quarks ] ! and sqrts_hat > 400 GeV
iterations = 5:{{ item.n_calls }}:"gw", 1:{{ 3 * item.n_calls }}
iterations = 10:{{ item.n_calls }}:"gw", 2:{{ 3 * item.n_calls }}
relative_error_goal = 0.005
......
......@@ -6,9 +6,7 @@
include("{{ SIN_LOCAL }}/ggW.sin")
process ggqqW = gl, gl => quark, antiquark, e1, N1
# calls(adapt) = 5e7, calls(integ) = 5e7
integrate (ggqqW) {iterations = 10:500000:"gw", 5:1000000}
integrate (g2) {iterations = 10:500000:"gw", 5:1000000}
printf "RES %8e %8e" (integral(ggqqW), error(ggqqW))
printf "RES %8e %8e" (integral(g2), error(g2))
......@@ -6,9 +6,7 @@
include("{{ SIN_LOCAL }}/ggW.sin")
process ggqqWg = gl, gl => quark, antiquark, e1, N1, g
# calls(adapt) = 5e7, calls(integ) = 5e7
integrate (ggqqWg) {iterations = 10:500000:"gw", 5:1000000}
integrate (g3) {iterations = 10:500000:"gw", 5:1000000}
show (results)
printf "RES %8e %8e" (integral(g3), error(g3))
......@@ -6,9 +6,7 @@
include("{{ SIN_LOCAL }}/ggW.sin")
process ggqqWgg = gl, gl => quark, antiquark, e1, N1, g, g
# calls(adapt) = 7.5e7, calls(integ) = 5e7
integrate (ggqqWgg) {iterations = 15:500000:"gw", 5:1000000}
integrate (g4) {iterations = 15:500000:"gw", 5:1000000}
show (results)
printf "RES %8e %8e" (integral(g4), error(g4))
......@@ -6,9 +6,7 @@
include("{{ SIN_LOCAL }}/ggW.sin")
process ggqqWggg = gl, gl => quark, antiquark, e1, N1, g, g, g
# calls(adapt) = 20e7, calls(integ) = 5e7
integrate (ggqqWggg) {iterations = 20:1000000:"gw", 5:1000000}
integrate (g5) {iterations = 20:1000000:"gw", 5:1000000}
show (results)
printf "RES %8e %8e" (integral(g5), error(g5))
......@@ -6,8 +6,6 @@
include("{{ SIN_LOCAL }}/ggW.sin")
process b0 = jet, jet => l, L, n, N, b, B
integrate (b0) {iterations = 10:100000:"gw", 5:1000000}
printf "RES %8e %8e" (integral(b0), error(b0))
......@@ -6,8 +6,6 @@
include("{{ SIN_LOCAL }}/ggW.sin")
process b3 = jet, jet => l, L, n, N, b, B, b, B
integrate (b3) {iterations = 10:100000:"gw", 5:1000000}
printf "RES %8e %8e" (integral(b3), error(b3))
......@@ -6,8 +6,6 @@
include("{{ SIN_LOCAL }}/ggW.sin")
process b1 = jet, jet => l, L, n, N, b, B, jet
integrate (b1) {iterations = 10:100000:"gw", 5:1000000}
printf "RES %8e %8e" (integral(b1), error(b1))
......@@ -6,8 +6,6 @@
include("{{ SIN_LOCAL }}/ggW.sin")
process b2 = jet, jet => l, L, n, N, b, B, jet, jet
integrate (b2) {iterations = 10:100000:"gw", 5:1000000}
printf "RES %8e %8e" (integral(b2), error(b2))
......@@ -6,8 +6,6 @@
include("{{ SIN_LOCAL }}/ggW.sin")
process j1 = jet, jet => e1, N1, jet
# calls(adapt) = 5e7, calls(integ) = 5e7
integrate (j1) {iterations = 10:500000:"gw", 5:1000000}
......
......@@ -6,8 +6,6 @@
include("{{ SIN_LOCAL }}/ggW.sin")
process j2 = jet, jet => e1, N1, jet, jet
# calls(adapt) = 5e7, calls(integ) = 5e7
integrate (j2) {iterations = 10:500000:"gw", 5:1000000}
......
......@@ -6,8 +6,6 @@
include("{{ SIN_LOCAL }}/ggW.sin")
process j3 = jet, jet => e1, N1, jet, jet, jet
# calls(adapt) = 5e7, calls(integ) = 5e7
integrate (j3) {iterations = 10:500000:"gw", 5:1000000}
......
......@@ -6,8 +6,6 @@
include("{{ SIN_LOCAL }}/ggW.sin")
process j4 = jet, jet => e1, N1, jet, jet, jet, jet
# calls(adapt) = 5e7, calls(integ) = 5e7
integrate (j4) {iterations = 10:500000:"gw", 5:1000000}
......
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