eebbuddubb.sin.j2 749 Bytes
Newer Older
Simon Brass's avatar
Simon Brass committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
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))