Ticket #397: foo2.sin

model = SM

alias q   = u:d:s:c
alias Q   = U:D:S:C
alias jet = q:Q:gl
alias bot = b:B
alias pr  = q:Q:b:B:gl

process dbub    = d, B => e1, N1, u, B
process dbubx   = d, b => e1, n1, u, b
process sbcbx   = s, b => e1, N1, c, b
process bdbu    = B, d => e1, N1, B, u
process bdbux   = b, d => e1, N1, b, u
process bsbcx   = b, s => e1, N1, b, c
process dbbubb1 = b:d, B:d => e1, N1, B:b, u
process dbbubb2 = B:d, b:d => e1, N1, B:b, u
process dbbubb3 = B:b, d => e1, N1, B:b, u
process dbbubb4 = d, B:b => e1, N1, B:b, u
process sbbcbb1 = b:d, B:s => e1, N1, B:b, u:c
process sbbcbb2 = B:s, b:d => e1, N1, B:b, u:c

compile

# Cuts
cuts = all Pt > 25 GeV [e1:N1]
   and all Pt > 30 GeV [bot]
   and all M2 < (225 GeV)^2 [bot,combine[e1,N1]]
   and all M2 > (150 GeV)^2 [bot,combine[e1,N1]]
   and all Pt > 50 GeV [jet]
   and (   all  2.5 < Eta <  5.0 [jet]
        or all -2.5 > Eta > -5.0 [jet])

# Beam specification and integration
sqrts = 14 TeV
beams = p,p => lhapdf 

integrate (dbub,bdbu,dbubx,bdbux,dbbubb1,dbbubb2,dbbubb3,dbbubb4) { iterations = 20:50000 }