Ticket #175: comp.sin

alias q = u:d
alias Q = U:D:S
alias j = u:d:U:D:g


process jets2 = e1, E1 -> u, U
process jets3 = e1, E1 -> u, U, g
process jets4 = e1, E1 -> u, U, g, g

compile

alias lepton = e1:e2:e3:E1:E2:E3
cuts = all Dist > 0.2 [j, j]
        and
        all Pt > 5 GeV [j]
!       and
!       all sqrt(Eta^2) < 5 [j]

sqrts = 500
iterations = 5:10000,3:10000

!?alpha_s_is_fixed = true
alphas = 0.118  
!alpha_s_order = 1

beams = e1, E1
integrate (jets2, jets3, jets4)

?write_lhef = true
$file_lhef = "mecomparision5.lhef"

n_events = 10
simulate (jets2, jets3, jets4)