Ticket #454: Wprod2.sin

model = SM

ms = 0
mc = 0
alias proton = u:d:s:c:U:D:S:C:g

process Wprod = proton, proton => e1:E1, n1:N1

sqrts = 1960 GeV
beams = p, pbar => lhapdf

integrate(Wprod) { iterations = 2:20000, 2:20000 }