Commit 2ed13e82 authored by Simon Brass's avatar Simon Brass

add LC generator group sin file

parent f9a496ba
model = SM
alias x = u:c
alias y = d:s:b
alias X = U:C
alias Y = D:S:B
alias uq = u:U
alias dq = d:D
alias cq = c:C
alias sq = s:S
alias bq = b:B
alias up_type_q = u:c:U:C
alias down_type_q = d:s:b:D:S:B
alias quarks = u:d:s:c:b:U:D:S:C:B
ms = 0 GeV
mc = 0 GeV
mb = 0 GeV
alphas = 0
process Pbbuddubb = e1, E1 => b,B,u,D,d,U,b,B
compile
logical ?is_higgs_process = true
sqrts = 500 GeV
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 ]
iterations = 5:{{ item.n_calls }}:"gw", 1:{{ 3 * item.n_calls }}
relative_error_goal = 0.005
integrate ( Pbbuddubb )
printf "RES %8i %8e %8e" ({{ item.n_calls }}, integral(Pbbuddubb), error(Pbbuddubb))
model = SM
alias x = u:c
alias y = d:s:b
alias X = U:C
alias Y = D:S:B
alias uq = u:U
alias dq = d:D
alias cq = c:C
alias sq = s:S
alias bq = b:B
alias up_type_q = u:c:U:C
alias down_type_q = d:s:b:D:S:B
alias quarks = u:d:s:c:b:U:D:S:C:B
ms = 0 GeV
mc = 0 GeV
mb = 0 GeV
alphas = 0
process Pbbuddubb = e1, E1 => b,B,u,D,d,U,b,B
compile
logical ?is_higgs_process = true
sqrts = 500 GeV
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 ]
iterations = 5:{{ item.n_calls }}:"gw", 5:{{ item.n_calls }}
relative_error_goal = 0.005
integrate ( Pbbuddubb )
printf "RES %8i %8e %8e" ({{ item.n_calls }}, integral(Pbbuddubb), error(Pbbuddubb))
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