Ticket #397: foo.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 ubdb    = u, b => E1, n1, d, b
process ubdbx   = u, B => E1, n1, d, B
process cbsbx   = c, B => E1, n1, s, B
process bubd    = b, u => E1, n1, b, d
process bubdx   = B, u => E1, n1, B, d
process bcbsx   = B, c => E1, n1, B, s
process ubbdbb1 = B:u, b:u => E1, n1, b:B, d
process ubbdbb2 = b:u, B:u => E1, n1, b:B, d
process ubbdbb3 = b:B, u => E1, n1, b:B, d
process ubbdbb4 = u, b:B => E1, n1, b:B, d
process cbbsbb1 = B:u, b:c => E1, n1, b:B, d:s
process cbbsbb2 = b:c, B:u => E1, n1, b:B, d:s

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 (ubdb,bubd,ubdbx,bubdx,ubbdbb1,ubbdbb2,ubbdbb3,ubbdbb4) { iterations = 20:50000 }