Ticket #464: hw.sin

########################################################
# Vector boson scattering 
# higgs width dependence study
# ulrike.schnoor@physik.tu-dresden.de
# alpha_s = 0: EW, alpha_s = a(mZ) from pdf: EWQCD
#########################################################
model = SM

# Set variables
real HiggsMass = 120 GeV
real HiggsWidth = 3 MeV

#set alpha_s (EW: 0; EWQCD: alpha_s from pdf @ 91.188)
scale = 91.188
alphas = 0

#?vis_channels = true
?phs_keep_nonresonant = true

# Set up the process
alias lepton = e1:E1:e2:E2
alias qf = u:d
process vbs1 = u:d, u:d => u, d, E1, e2, nue, numubar 
process vbs2 = u:d, u:d => u, d, E1, e2, nue, numubar 
process vbs3 = u:d, u:d => u, d, E1, e2, nue, numubar 

# Set up parameters
ms = 0   mtau = 0
mH = HiggsMass wH = HiggsWidth
me = 0  mmu = 0
mc = 0

# specify cuts
cuts = all Pt > 15 GeV [lepton] and all Pt > 15 GeV [qf] and all Dist > 1.0 [qf,qf] and all M > 20 GeV [lepton,lepton]

compile

# LHC energy
sqrts = 14 TeV
beams = p, p => pdf_builtin

#?phs_only = true
?phs_step_mapping = false
integrate (vbs1) { iterations = 10:1000000, 5:1000000  }

?phs_step_mapping = true
?phs_step_mapping_exp = true
integrate (vbs2) { iterations = 10:1000000, 5:1000000  }

?phs_step_mapping = true
?phs_step_mapping_exp = false
integrate (vbs3) { iterations = 10:1000000, 5:1000000  }


#n_events = 1000
#sample_format = hepmc
#?unweighted = false
#simulate (vbs) { checkpoint = 100 }
show(results)