Ticket #270: lepjets.sin

model = SM

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


process jets2 = e1, E1 -> q, Q
process jets3 = e1, E1 -> q, Q, g
process jets4 = e1, E1 -> q, Q, g, g 
!process jets5 = e1, E1 -> q, Q, g, g, g 

compile

alias lepton = e1:e2:e3:E1:E2:E3
cuts = all Dist > 0.7 [j, j]
        and
        all Pt > 20 GeV [j]

sqrts = 500
iterations = 5:20000,5:20000

?read_color_factors = false

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

?write_lhef = true

n_events = 200
simulate (jets2, jets3, jets4)!, jets5)


exec ("../bin/interface")
exec ("rm *.phs")