Ticket #269: test.sin

model = Threeshl

alias pr = u:U:d:D:g

process born = pr, pr -> e2, N2, u, D

compile

mhw = 500

cuts = all Pt > 50 GeV [e2]
       and all Pt > 50 GeV [N2]
       and all Pt > 50 GeV [u]
       and all Pt > 50 GeV [D]
       and all 70 GeV < M < 90 GeV [e2, N2]
       and all 70 GeV < M < 90 GeV [u, D]
 
sqrts = 7 TeV

beams = p, p -> lhapdf

integrate (born) { iterations = 5:5000, 10:5000 }

histogram m_inv_mmqq (450,550,5)

n_events = 10000

analysis = record m_inv_mmqq (eval M / 1 GeV [combine[combine[e2,N2],combine[u,D]]])

?write_lhef = true

$file_lhef = "test.lhe"

simulate (born)

$analysis_filename = "test"
write_analysis