Commit 3ab24c33 authored by Simon Brass's avatar Simon Brass

add arXiv:0708.4233 validation

parent 5338f404
......@@ -3,15 +3,18 @@ localhost ansible_connection=local
[notebook:vars]
SIN_LOCAL=/home/simon/workspace/run/sindarin
WHIZARD_SCRIPTS=/home/simon/whizard-scripts
[caesium]
caesium ansible_user=brass
[caesium:vars]
SIN_LOCAL=/home/users/brass/workspace/run/sindarin
WHIZARD_SCRIPTS=/home/users/brass/whizard-scripts
[cluster]
slc ansible_user=gg662
[cluster:vars]
SIN_LOCAL=/home/gg662/workspace/run/sindarin
\ No newline at end of file
SIN_LOCAL=/home/gg662/workspace/run/sindarin
WHIZARD_SCRIPTS=/home/gg662/whizard-scripts
\ No newline at end of file
......@@ -11,7 +11,11 @@ alias quark = d:u:s:c
alias antiquark = D:U:S:C
# Aliases for cuts' definition
alias lepton = e1:e2:e3:E1:E2:E3
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
me = 0
......@@ -19,16 +23,18 @@ mmu = 0
mtau = 0
mW = 80.419
mZ = 91.1882
mH = 130
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
mtop = 174
mtop = 174.0 GeV
wtop = 1.523 GeV
sqrts = 14 TeV
beams = p, p => pdf_builtin
......
......@@ -3,7 +3,6 @@
# --------------------------------------------------
- name: copy configuraion file for e1, E1 => 6f
copy:
src: ee6f.sin
......@@ -29,7 +28,27 @@
# --------------------------------------------------
- name: deploy e1, E1 => e2, E2
- name: deploy e1, E1 => b, B, u, D, d, U, b, B (w/o beampol. and isr) (500 GeV)
template:
src: eebbuddubb-simple.sin.j2
dest: "{{ SIN_LOCAL }}/eebbuddubb-simple-{{ item.name }}.sin"
with_items:
- { name: 5e6, n_calls: 500000 }
- { name: 1e7, n_calls: 1000000 }
- { name: 2e7, n_calls: 2000000 }
- name: deploy e1, E1 => b, B, u, D, d, U, b, B (500 GeV)
template:
src: eebbuddubb.sin.j2
dest: "{{ SIN_LOCAL }}/eebbuddubb-{{ item.name }}.sin"
with_items:
- { name: 5e6, n_calls: 500000 }
- { name: 1e7, n_calls: 1000000 }
- { name: 2e7, n_calls: 2000000 }
# --------------------------------------------------
- name: deploy e1, E1 => e2, E2 (150 GeV)
template:
src: eemumu.sin.j2
dest: "{{ SIN_LOCAL }}/eemumu-{{ item.name }}.sin"
......@@ -38,32 +57,85 @@
- { name: 10e6, n_calls: 100000 }
- { name: 10e7, n_calls: 1000000 }
- name: deploy e1, E1 => e2, E2 (50 GeV)
template:
src: ee2mu.sin.j2
dest: "{{ SIN_LOCAL }}/ee2mu-{{ item.name }}.sin"
with_items:
- { name: 10e5, n_calls: 10000 }
- { name: 10e6, n_calls: 100000 }
- { name: 10e7, n_calls: 1000000 }
# --------------------------------------------------
- name: copy basic configuration for gl, gl => WW, qq
copy:
src: ggW.sin
dest: "{{ SIN_LOCAL }}/ggW.sin"
- block:
- name: copy basic configuration for gl, gl => WW, qq
copy:
src: ggW.sin
dest: "{{ SIN_LOCAL }}/ggW.sin"
- name: deploy gl, gl => W, W, q, q
template:
src: ggqqW.sin.j2
dest: "{{ SIN_LOCAL }}/ggqqW.sin"
- name: deploy gl, gl => W, q, q, gl {0, 3}
template:
src: "{{ item }}"
dest: "{{SIN_LOCAL }}/{{ item |replace('.j2', '') }}"
with_items:
- ggqqW.sin.j2
- ggqqWg.sin.j2
- ggqqWgg.sin.j2
- ggqqWggg.sin.j2
- name: deploy gl, gl => W, W, q, q, gl
template:
src: ggqqW.sin.j2
dest: "{{ SIN_LOCAL }}/ggqqWg.sin"
- name: collect gl, gl => W, q, q, gl {0, 3}
shell: "{{ WHIZARD_SCRIPTS }}/slurm/collect-result.sh -f results/{{ item }}.dat -h '# {{ item }}' benchmark-mpi/{{ item }}"
args:
chdir: ~/workspace/run
with_items:
- ggqqW
- ggqqWg
- ggqqWgg
- ggqqWggg
- name: deploy gl, gl => W, W, q, q, gl, gl
template:
src: ggqqWgg.sin.j2
dest: "{{ SIN_LOCAL }}/ggqqWgg.sin"
- name: deploy jet, jet => W, jet(1-4)
template:
src: "{{ item }}"
dest: "{{SIN_LOCAL }}/{{ item |replace('.j2', '') }}"
with_items:
- jjWj.sin.j2
- jjWjj.sin.j2
- jjWjjj.sin.j2
- jjWjjjj.sin.j2
- name: deploy gl, gl => W, W, q, q, gl, gl, gl
template:
src: ggqqWggg.sin.j2
dest: "{{ SIN_LOCAL }}/ggqqWggg.sin"
- name: collect jet, jet => W, jet {1,4}
shell: "{{ WHIZARD_SCRIPTS }}/slurm/collect-result.sh -f results/{{ item }}.dat -h '# {{ item }}' benchmark-mpi/{{ item }}"
args:
chdir: ~/workspace/run
with_items:
- jjWj
- jjWjj
- jjWjjj
- jjWjjjj
- name: deploy jet, jet => W, l, l, *
template:
src: "{{ item }}"
dest: "{{SIN_LOCAL }}/{{ item |replace('.j2', '') }}"
with_items:
- jjWWbb.sin.j2
- jjWWbbj.sin.j2
- jjWWbbjj.sin.j2
- jjWWbbbb.sin.j2
- name: collect jet, jet => W, jet {1,4}
shell: "{{ WHIZARD_SCRIPTS }}/slurm/collect-result.sh -f results/{{ item }}.dat -h '# {{ item }}' benchmark-mpi/{{ item }}"
args:
chdir: ~/workspace/run
with_items:
- jjWWbb
- jjWWbbj
- jjWWbbjj
- jjWWbbbb
tags:
- whizard
# --------------------------------------------------
......@@ -72,6 +144,9 @@
args:
chdir: ~/workspace/run
with_items:
- { name: ee2mu-10e5 }
- { name: ee2mu-10e6 }
- { name: ee2mu-10e7 }
- { name: ee4mu-10e5 }
- { name: ee4mu-10e6 }
- { name: ee4mu-10e7 }
......@@ -85,5 +160,8 @@
- { name: ee4mu2nu-10e5 }
- { name: ee4mu2nu-10e6 }
- { name: ee4mu2nu-10e7 }
- { name: eebbuddubb-5e6 }
- { name: eebbuddubb-1e7 }
- { name: eebbuddubb-2e7 }
tags:
collect
# arXiv:0708.4233v2
#
# Section 3.1.2
~~~~#
#
# Table 1, Row 1
include("{{ SIN_LOCAL }}/ggW.sin")
process ggqqW = gl, gl~/`~ => quark, antiquark, e1, N1
process ggqqW = gl, gl => quark, antiquark, e1, N1
# calls(adapt) = 5e7, calls(integ) = 5e7
integrate (ggqqW) {iterations = 10:500000:"gw", 5:1000000}
show (results)
printf "RES %8e %8e" (integral(ggqqW), error(ggqqW))
# arXiv:0708.4233v2
#
# Section 3.1.2
#
# Table 1, Row 2
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))
# arXiv:0708.4233v2
#
# Section 3.1.2
#
# Table 1, Row 2
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))
# arXiv:0708.4233v2
#
# Section 3.1.2
#
# Table 1, Row 2
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))
# arXiv:0708.4233v2
#
# Section 3.1.2
#
# Table 1, Row 2
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))
# arXiv:0708.4233v2
#
# Section 3.1.2
#
# Table 1, Row 2
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}
printf "RES %8e %8e" (integral(j1), error(j1))
# arXiv:0708.4233v2
#
# Section 3.1.2
#
# Table 1, Row 2
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}
printf "RES %8e %8e" (integral(j2), error(j2))
# arXiv:0708.4233v2
#
# Section 3.1.2
#
# Table 1, Row 2
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}
printf "RES %8e %8e" (integral(j3), error(j3))
# arXiv:0708.4233v2
#
# Section 3.1.2
#
# Table 1, Row 2
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}
printf "RES %8e %8e" (integral(j4), error(j4))
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