Ticket #668: vbs.v2.2.2.sin

# Vector boson scattering with WHIZARD
########################################################
# ulrike.schnoor@physik.tu-dresden.de
# anomalous couplings a4, a5 in SM_km.mdl
# alpha_s = 0: EW
########################################################
model = SM #_km

#set alpha_s EW
alphas = 0

# Set up the process
alias nu = n1:n2:n3
alias NU = N1:N2:N3
alias neutrino = n1:n2:n3:N1:N2:N3
alias lep = e1:e2#:e3
alias LEP = E1:E2#:E3
alias lepton = e1:E1:e2:E2:e3:E3
alias qf = u:d:s:c:U:D:S:C:g
alias allleptons = lepton:neutrino

process vbslplp = qf, qf => qf, qf, LEP, LEP, nu, nu
process vbslmlm = qf, qf => qf, qf, lep, lep, NU, NU
#process vbstplp = qf, qf => qf, qf, E3, LEP, n3, nu
#process vbstmlm = qf, qf => qf, qf, e3, lep, N3, NU
#process vbstauptaup = qf, qf => qf, qf, E3, E3, n3, n3
#process vbstaumtaum = qf, qf => qf, qf, e3, e3, N3, N3

#===============================================================================

#ATLAS settings
mtop = 172.5
mW = 80.399
mZ = 91.1876
wZ = 2.4952
wW = 2.085

# Set up parameters
scale = eval M [collect [allleptons]]
ms = 0   mtau = 1.77705
mH = 126 GeV
wH = 0.00418 GeV
me = 0  mmu = 0
mc = 0

?phs_keep_nonresonant = true

# Model parameters  =============================================================
#a4 = 0
#a5 = 0

#VBS cuts fiducial volume ======================================================
cuts = all Dist > 0.4 [qf,qf] and all -5.0 < Eta < 5.0 [qf] and all Pt > 10
 GeV [qf] and all -5.0 < Eta < 5.0 [lepton] and all Pt > 8 GeV [lepton]

compile

# LHC energy===================================================================
sqrts = 14 TeV        
$lhapdf_file = "cteq6ll.LHpdf"                                
beams = p, p => lhapdf
#================================================================================                     
integrate (vbslplp,vbslmlm)  {iterations= 24:80000, 20:100000}                 
show(results)                           

#================================================================================                     
n_events = 50000
#sample_format = hepmc
sample_format = lhef
?unweighted = true
simulate (vbslplp,vbslmlm) {checkpoint = 500 ?keep_beams = false}