Ticket #436: tty_2step_2_1.sin

model = SM

mc=0 GeV
ms=0 GeV
me=0 GeV
mmu=0 GeV

alias ini1 = gl:u:U:d:D:s:S:c:C
alias ini2 = gl:u:U:d:D:s:S:c:C
alias f = u:d:c:s:e1:e2:n1:n2
alias af = U:D:C:S:E1:E2:N1:N2
alias l = e1:e2
alias n = n1:n2
alias al = E1:E2
alias an = N1:N2
alias q = u:d:c:s
alias aq = U:D:C:S
alias taus = e3:E3

process ttbar = ini1,ini2 => t,T
process t_dec_had = t:T => b:B,q:aq,aq:q
process t_dec_lep = t:T => b:B,al:l,n:an
process t_dec_tau = t:T => b:B,E3:e3,n3:N3
process t_dec_gamma_had = t:T => b:B,q:aq,aq:q,A
process t_dec_gamma_lep = t:T => b:B,al:l,n:an,A
process t_dec_gamma_tau = t:T => b:B,E3:e3,n3:N3,A

compile

mtop=172.5 GeV
wtop=1.523 GeV

mH=115 GeV
wH=0.3228E-02 GeV

sqrts = 7 TeV
scale = 2*mtop

beams = p,p => lhapdf {$lhapdf_file = "cteq6ll.LHpdf"}

?rebuild_phase_space = true
?rebuild_grids = true
?rebuild_events = true

?read_color_factors = true

?alpha_s_is_fixed = false
#alphas = 0.10689
?alpha_s_from_lhapdf = true
#alpha_s_order =            0
#alpha_s_nf =            5
#?alpha_s_from_mz = false
# lambda_qcd = 0.168 GeV

cuts = all Pt > 8 GeV [A]
and all Pt > 10 GeV [q]
and all Pt > 10 GeV [aq]
and all Pt > 10 GeV [l]
and all Pt > 10 GeV [al]
and all Pt > 10 GeV [taus]
and all M > 5 GeV [q,aq]
and all M > 5 GeV [A,aq]
and all M > 5 GeV [A,q]
and all M > 5 GeV [A,l]
and all M > 5 GeV [A,al]
and all M > 5 GeV [A,taus]
and all M > 5 GeV [ini1,A]
and all M > 5 GeV [ini2,A]
and all M > 5 GeV [gl,q]
and all M > 5 GeV [gl,aq]

integrate (ttbar)

sample_format = lhef

n_events = 50000
unstable t(t_dec_had)
unstable T(t_dec_gamma_lep,t_dec_gamma_tau)

simulate(ttbar) {$sample="tty_2step_2_1"}