2.6.3.2
t1 t1~ @0
add process p p > t1 t1~ j@1
add process p p > t1 t1~ j j@2
output stop
]]>
t1 t1~ @0 #Process
# Be carefull the coupling are here in MG5 convention
end_coup # End the couplings input
done # this tells MG there are no more procs
# End PROCESS # This is TAG. Do not modify this line
#*********************************************************************
# Model information *
#*********************************************************************
# Begin MODEL # This is TAG. Do not modify this line
MSSM_SLHA2
# End MODEL # This is TAG. Do not modify this line
#*********************************************************************
# Start multiparticle definitions *
#*********************************************************************
# Begin MULTIPARTICLES # This is TAG. Do not modify this line
# End MULTIPARTICLES # This is TAG. Do not modify this line
]]>
0
! (turn off for VBF and single top processes)
100 = xqcut ! minimum kt jet measure between partons
#*********************************************************************
#
#*********************************************************************
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#*********************************************************************
0 = nhel ! using helicities importance sampling or not.
#*********************************************************************
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
#*********************************************************************
None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
{} = bias_parameters ! Specifies the parameters of the module.
#
#*******************************
# Parton level cuts definition *
#*******************************
#
#
#*********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
#*********************************************************************
15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
#*********************************************************************
# Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#*********************************************************************
False = cut_decays ! Cut decay products
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
20.0 = ptj ! minimum pt for the jets
20.0 = ptb ! minimum pt for the b
20.0 = pta ! minimum pt for the photons
15.0 = ptl ! minimum pt for the charged leptons
0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
-1.0 = ptjmax ! maximum pt for the jets
-1.0 = ptbmax ! maximum pt for the b
-1.0 = ptamax ! maximum pt for the photons
-1.0 = ptlmax ! maximum pt for the charged leptons
-1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
{} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
{} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#*********************************************************************
# Minimum and maximum E's (in the center of mass frame) *
#*********************************************************************
0.0 = ej ! minimum E for the jets
0.0 = eb ! minimum E for the b
0.0 = ea ! minimum E for the photons
0.0 = el ! minimum E for the charged leptons
-1.0 = ejmax ! maximum E for the jets
-1.0 = ebmax ! maximum E for the b
-1.0 = eamax ! maximum E for the photons
-1.0 = elmax ! maximum E for the charged leptons
{} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle
{} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50})
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
2.8 = etaj ! max rap for the jets
2.5 = etab ! max rap for the b
2.37 = etaa ! max rap for the photons
2.7 = etal ! max rap for the charged leptons
0.0 = etajmin ! min rap for the jets
0.0 = etabmin ! min rap for the b
0.0 = etaamin ! min rap for the photons
0.0 = etalmin ! main rap for the charged leptons
{} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
{} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
0.4 = drjj ! min distance between jets
0.4 = drbb ! min distance between b's
0.4 = drll ! min distance between leptons
0.4 = draa ! min distance between gammas
0.4 = drbj ! min distance between b and jet
0.4 = draj ! min distance between gamma and jet
0.2 = drjl ! min distance between jet and lepton
0.4 = drab ! min distance between gamma and b
0.4 = drbl ! min distance between b and lepton
0.4 = dral ! min distance between gamma and lepton
-1.0 = drjjmax ! max distance between jets
-1.0 = drbbmax ! max distance between b's
-1.0 = drllmax ! max distance between leptons
-1.0 = draamax ! max distance between gammas
-1.0 = drbjmax ! max distance between b and jet
-1.0 = drajmax ! max distance between gamma and jet
-1.0 = drjlmax ! max distance between jet and lepton
-1.0 = drabmax ! max distance between gamma and b
-1.0 = drblmax ! max distance between b and lepton
-1.0 = dralmax ! maxdistance between gamma and lepton
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
# WARNING: for four lepton final state mmll cut require to have *
# different lepton masses for each flavor! *
#*********************************************************************
0.0 = mmjj ! min invariant mass of a jet pair
0.0 = mmbb ! min invariant mass of a b pair
0.0 = mmaa ! min invariant mass of gamma gamma pair
0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
-1.0 = mmjjmax ! max invariant mass of a jet pair
-1.0 = mmbbmax ! max invariant mass of a b pair
-1.0 = mmaamax ! max invariant mass of gamma gamma pair
-1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
{} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
{'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
#*********************************************************************
# Minimum and maximum invariant mass for all letpons *
#*********************************************************************
0.0 = mmnl ! min invariant mass for all letpons (l+- and vl)
-1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#*********************************************************************
# Minimum and maximum pt for 4-momenta sum of leptons *
#*********************************************************************
0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
-1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
0.0 = ptheavy ! minimum pt for at least one heavy final state
0.0 = xptj ! minimum pt for at least one jet
0.0 = xptb ! minimum pt for at least one b
0.0 = xpta ! minimum pt for at least one photon
0.0 = xptl ! minimum pt for at least one charged lepton
#*********************************************************************
# Control the pt's of the jets sorted by pt *
#*********************************************************************
0.0 = ptj1min ! minimum pt for the leading jet in pt
0.0 = ptj2min ! minimum pt for the second jet in pt
0.0 = ptj3min ! minimum pt for the third jet in pt
0.0 = ptj4min ! minimum pt for the fourth jet in pt
-1.0 = ptj1max ! maximum pt for the leading jet in pt
-1.0 = ptj2max ! maximum pt for the second jet in pt
-1.0 = ptj3max ! maximum pt for the third jet in pt
-1.0 = ptj4max ! maximum pt for the fourth jet in pt
0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#*********************************************************************
# Control the pt's of leptons sorted by pt *
#*********************************************************************
0.0 = ptl1min ! minimum pt for the leading lepton in pt
0.0 = ptl2min ! minimum pt for the second lepton in pt
0.0 = ptl3min ! minimum pt for the third lepton in pt
0.0 = ptl4min ! minimum pt for the fourth lepton in pt
-1.0 = ptl1max ! maximum pt for the leading lepton in pt
-1.0 = ptl2max ! maximum pt for the second lepton in pt
-1.0 = ptl3max ! maximum pt for the third lepton in pt
-1.0 = ptl4max ! maximum pt for the fourth lepton in pt
#*********************************************************************
# Control the Ht(k)=Sum of k leading jets *
#*********************************************************************
0.0 = htjmin ! minimum jet HT=Sum(jet pt)
-1.0 = htjmax ! maximum jet HT=Sum(jet pt)
0.0 = ihtmin !inclusive Ht for all partons (including b)
-1.0 = ihtmax !inclusive Ht for all partons (including b)
0.0 = ht2min ! minimum Ht for the two leading jets
0.0 = ht3min ! minimum Ht for the three leading jets
0.0 = ht4min ! minimum Ht for the four leading jets
-1.0 = ht2max ! maximum Ht for the two leading jets
-1.0 = ht3max ! maximum Ht for the three leading jets
-1.0 = ht4max ! maximum Ht for the four leading jets
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
# When ptgmin>0, pta and draj are not going to be used *
#***********************************************************************
0.0 = ptgmin ! Min photon transverse momentum
0.4 = R0gamma ! Radius of isolation code
1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
True = isoEM ! isolate photons from EM energy (photons and leptons)
#*********************************************************************
# WBF cuts *
#*********************************************************************
0.0 = xetamin ! minimum rapidity for two jets in the WBF case
0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
#***********************************************************************
# Turn on either the ktdurham or ptlund cut to activate *
# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
#***********************************************************************
-1.0 = ktdurham
0.4 = dparameter
-1.0 = ptlund
1, 2, 3, 4, 5, 6, 21 = pdgs_for_merging_cut ! PDGs for two cuts above
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
5 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
False = use_syst ! Enable systematics studies
#
#**************************************
# Parameter of the systematics study
# will be used by SysCalc (if installed)
#**************************************
#
0.5 1 2 = sys_scalefact # factorization/renormalization scale factor
None = sys_alpsfact # \alpha_s emission scale factors
auto = sys_matchscale # variation of merging scale
# PDF sets and number of members (0 or none for all members).
#NNPDF23_lo_as_0130_qed = sys_pdf # list of pdf sets
False = sys_pdf
#
]]>
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR DCINFO
###################################
BLOCK DCINFO # decay program information
1 sdecay/hdecay # decay calculator
2 1.5a /3.4 # version number
###################################
## INFORMATION FOR SPINFO
###################################
BLOCK SPINFO # spectrum calculator information
1 feynhiggs #
2 built on jun 28, 2018 #
###################################
## INFORMATION FOR MODSEL
###################################
BLOCK MODSEL # model selection
0 0 #
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS # standard model inputs
1 1.271180e+02 # alpha_em^-1(m_z)^msbar
3 1.187000e-01 # alpha_s(m_z)^msbar
5 4.180000e+00 # mb(mb)^msbar
###################################
## INFORMATION FOR MINPAR
###################################
BLOCK MINPAR # input parameters - minimal models
1 5.774405e+02 # m0
2 -4.544367e+02 # m12
3 3.670866e+01 # tb
5 1.906052e+03 # a
###################################
## INFORMATION FOR EXTPAR
###################################
BLOCK EXTPAR # input parameters - non-minimal models
0 3.681096e+02 # q
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS # mass spectrum
24 8.025851e+01 # w+
25 1.231921e+02 # h
35 8.146405e+03 # h
36 8.145861e+03 # a
37 8.146249e+03 # h+
5 4.823272e+00 # b-quark pole mass calculated from mb(mb)_msbar
1000001 2.398410e+03 # ~d_l
2000001 1.417294e+03 # ~d_r
1000002 2.397246e+03 # ~u_l
2000002 2.038783e+03 # ~u_r
1000003 2.398410e+03 # ~s_l
2000003 1.417294e+03 # ~s_r
1000004 2.397246e+03 # ~c_l
2000004 2.038783e+03 # ~c_r
1000005 6.304867e+02 # ~b_1
2000005 3.799771e+03 # ~b_2
1000006 6.098075e+02 # ~t_1
2000006 1.095659e+03 # ~t_2
1000011 2.596343e+03 # ~e_l
2000011 3.202957e+03 # ~e_r
1000012 2.594803e+03 # ~nu_el
1000013 2.596343e+03 # ~mu_l
2000013 3.202957e+03 # ~mu_r
1000014 2.594803e+03 # ~nu_mul
1000015 8.084450e+02 # ~tau_1
2000015 2.534642e+03 # ~tau_2
1000016 8.044583e+02 # ~nu_taul
1000021 1.892563e+03 # ~g
1000022 -4.365357e+02 # ~chi_10
1000023 4.996586e+02 # ~chi_20
1000025 -5.515240e+02 # ~chi_30
1000035 6.033017e+02 # ~chi_40
1000024 4.371662e+02 # ~chi_1+
1000037 5.531292e+02 # ~chi_2+
###################################
## INFORMATION FOR NMIX
###################################
BLOCK NMIX # neutralino mixing matrix
1 1 1.519905e-02 # n_11
1 2 8.311307e-01 # n_12
1 3 -4.272709e-01 # n_13
1 4 -3.555703e-01 # n_14
2 1 -3.536438e-01 # n_21
2 2 -5.567011e-02 # n_22
2 3 -6.657189e-01 # n_23
2 4 6.547177e-01 # n_24
3 1 -2.354677e-02 # n_31
3 2 5.528982e-01 # n_32
3 3 5.666454e-01 # n_33
3 4 6.104606e-01 # n_34
4 1 9.349602e-01 # n_41
4 2 -2.064346e-02 # n_42
4 3 -2.305880e-01 # n_43
4 4 2.687981e-01 # n_44
###################################
## INFORMATION FOR UMIX
###################################
BLOCK UMIX # chargino mixing matrix u
1 1 8.041869e-01 # u_11
1 2 -5.943765e-01 # u_12
2 1 5.943765e-01 # u_21
2 2 8.041869e-01 # u_22
###################################
## INFORMATION FOR VMIX
###################################
BLOCK VMIX # chargino mixing matrix v
1 1 -8.688338e-01 # v_11
1 2 -4.951038e-01 # v_12
2 1 -4.951038e-01 # v_21
2 2 8.688338e-01 # v_22
###################################
## INFORMATION FOR HMIX
###################################
BLOCK HMIX Q= 3.681096e+02 # drbar higgs parameters
1 5.021106e+02 # mue
2 3.629302e+01 # tb
4 6.577630e+07 # ma02
###################################
## INFORMATION FOR GAUGE
###################################
BLOCK GAUGE Q= 3.681096e+02 # the gauge couplings
1 3.599444e-01 # gprime(q) drbar
2 6.386130e-01 # g(q) drbar
3 1.073372e+00 # g3(q) drbar
###################################
## INFORMATION FOR YU
###################################
BLOCK YU Q= 3.681096e+02 # the yukawa couplings
1 1 0.000000e+00 # y_u(q) drbar
2 2 0.000000e+00 # y_c(q) drbar
3 3 8.306962e-01 # y_t(q) drbar
###################################
## INFORMATION FOR YD
###################################
BLOCK YD Q= 3.681096e+02 # the yukawa couplings
1 1 0.000000e+00 # y_d(q) drbar
2 2 0.000000e+00 # y_s(q) drbar
3 3 4.933995e-01 # y_b(q) drbar
###################################
## INFORMATION FOR YE
###################################
BLOCK YE Q= 3.681096e+02 # the yukawa couplings
1 1 0.000000e+00 # y_e(q) drbar
2 2 0.000000e+00 # y_mu(q) drbar
3 3 3.330188e-01 # y_tau(q) drbar
###################################
## INFORMATION FOR MSOFT
###################################
BLOCK MSOFT Q= 3.681096e+02 # the soft susy breaking masses at the scale q
1 5.774405e+02 # m1
2 -4.544367e+02 # m2
3 1.906052e+03 # m3
21 6.384664e+07 # m^2_hd
22 -1.515615e+05 # m^2_hu
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.439271e+00 # top decays
1.000000e+00 2 5 24 # br(t -> b w+)
DECAY 1000021 8.314459e+01 # gluino decays
3.912023e-02 2 -1 2000001 # br(~g -> ~d_r db)
3.912023e-02 2 -2000001 1 # br(~g -> ~d_r* d )
3.912023e-02 2 -3 2000003 # br(~g -> ~s_r sb)
3.912023e-02 2 -2000003 3 # br(~g -> ~s_r* s )
1.647008e-01 2 -5 1000005 # br(~g -> ~b_1 bb)
1.647008e-01 2 -1000005 5 # br(~g -> ~b_1* b )
1.609677e-01 2 -6 1000006 # br(~g -> ~t_1 tb)
1.609677e-01 2 -1000006 6 # br(~g -> ~t_1* t )
9.609103e-02 2 -6 2000006 # br(~g -> ~t_2 tb)
9.609103e-02 2 -2000006 6 # br(~g -> ~t_2* t )
DECAY 1000006 1.096631e+00 # stop1 decays
3.292823e-04 2 6 1000022 # br(~t_1 -> ~chi_10 t )
8.997037e-01 2 5 1000024 # br(~t_1 -> ~chi_1+ b )
9.996706e-02 2 5 1000037 # br(~t_1 -> ~chi_2+ b )
DECAY 2000006 2.960240e+01 # stop2 decays
4.996050e-02 2 6 1000022 # br(~t_2 -> ~chi_10 t )
1.131181e-01 2 6 1000023 # br(~t_2 -> ~chi_20 t )
7.490230e-02 2 6 1000025 # br(~t_2 -> ~chi_30 t )
6.232851e-02 2 6 1000035 # br(~t_2 -> ~chi_40 t )
1.065310e-01 2 5 1000024 # br(~t_2 -> ~chi_1+ b )
1.720437e-01 2 5 1000037 # br(~t_2 -> ~chi_2+ b )
1.809620e-01 2 25 1000006 # br(~t_2 -> ~t_1 h )
8.729333e-02 2 23 1000006 # br(~t_2 -> ~t_1 z )
1.528606e-01 2 24 1000005 # br(~t_2 -> ~b_1 w+)
DECAY 1000005 9.932363e-01 # sbottom1 decays
6.362690e-01 2 5 1000022 # br(~b_1 -> ~chi_10 b )
1.703730e-01 2 5 1000023 # br(~b_1 -> ~chi_20 b )
9.095286e-02 2 5 1000025 # br(~b_1 -> ~chi_30 b )
2.159448e-03 2 5 1000035 # br(~b_1 -> ~chi_40 b )
1.002456e-01 2 -1000024 6 # br(~b_1 -> ~chi_1- t )
DECAY 2000005 1.639773e+02 # sbottom2 decays
1.690292e-02 2 5 1000022 # br(~b_2 -> ~chi_10 b )
4.252691e-02 2 5 1000023 # br(~b_2 -> ~chi_20 b )
2.962510e-02 2 5 1000025 # br(~b_2 -> ~chi_30 b )
1.549703e-02 2 5 1000035 # br(~b_2 -> ~chi_40 b )
3.261208e-02 2 -1000024 6 # br(~b_2 -> ~chi_1- t )
5.948367e-02 2 -1000037 6 # br(~b_2 -> ~chi_2- t )
8.000338e-01 2 5 1000021 # br(~b_2 -> ~g b )
5.323226e-04 2 25 1000005 # br(~b_2 -> ~b_1 h )
1.041491e-03 2 23 1000005 # br(~b_2 -> ~b_1 z )
1.709503e-03 2 -24 1000006 # br(~b_2 -> ~t_1 w-)
3.518365e-05 2 -24 2000006 # br(~b_2 -> ~t_2 w-)
DECAY 1000002 5.207067e+01 # sup_l decays
1.153393e-01 2 2 1000022 # br(~u_l -> ~chi_10 u)
2.421794e-03 2 2 1000023 # br(~u_l -> ~chi_20 u)
4.791548e-02 2 2 1000025 # br(~u_l -> ~chi_30 u)
3.745550e-03 2 2 1000035 # br(~u_l -> ~chi_40 u)
2.503106e-01 2 1 1000024 # br(~u_l -> ~chi_1+ d)
7.803853e-02 2 1 1000037 # br(~u_l -> ~chi_2+ d)
5.022288e-01 2 2 1000021 # br(~u_l -> ~g u)
DECAY 2000002 6.723905e+00 # sup_r decays
1.387400e-04 2 2 1000022 # br(~u_r -> ~chi_10 u)
7.293364e-02 2 2 1000023 # br(~u_r -> ~chi_20 u)
3.145366e-04 2 2 1000025 # br(~u_r -> ~chi_30 u)
4.808707e-01 2 2 1000035 # br(~u_r -> ~chi_40 u)
4.457424e-01 2 2 1000021 # br(~u_r -> ~g u)
DECAY 1000001 5.209967e+01 # sdown_l decays
1.137659e-01 2 1 1000022 # br(~d_l -> ~chi_10 d)
1.884572e-05 2 1 1000023 # br(~d_l -> ~chi_20 d)
4.947629e-02 2 1 1000025 # br(~d_l -> ~chi_30 d)
6.006867e-03 2 1 1000035 # br(~d_l -> ~chi_40 d)
2.144486e-01 2 -1000024 2 # br(~d_l -> ~chi_1- u)
1.124764e-01 2 -1000037 2 # br(~d_l -> ~chi_2- u)
5.038072e-01 2 1 1000021 # br(~d_l -> ~g d)
DECAY 2000001 5.508143e-01 # sdown_r decays
2.768362e-04 2 1 1000022 # br(~d_r -> ~chi_10 d)
1.403730e-01 2 1 1000023 # br(~d_r -> ~chi_20 d)
5.846603e-04 2 1 1000025 # br(~d_r -> ~chi_30 d)
8.587655e-01 2 1 1000035 # br(~d_r -> ~chi_40 d)
DECAY 1000004 5.207067e+01 # scharm_l decays
1.153393e-01 2 4 1000022 # br(~c_l -> ~chi_10 c)
2.421794e-03 2 4 1000023 # br(~c_l -> ~chi_20 c)
4.791548e-02 2 4 1000025 # br(~c_l -> ~chi_30 c)
3.745550e-03 2 4 1000035 # br(~c_l -> ~chi_40 c)
2.503106e-01 2 3 1000024 # br(~c_l -> ~chi_1+ s)
7.803853e-02 2 3 1000037 # br(~c_l -> ~chi_2+ s)
5.022288e-01 2 4 1000021 # br(~c_l -> ~g c)
DECAY 2000004 6.723905e+00 # scharm_r decays
1.387400e-04 2 4 1000022 # br(~c_r -> ~chi_10 c)
7.293364e-02 2 4 1000023 # br(~c_r -> ~chi_20 c)
3.145366e-04 2 4 1000025 # br(~c_r -> ~chi_30 c)
4.808707e-01 2 4 1000035 # br(~c_r -> ~chi_40 c)
4.457424e-01 2 4 1000021 # br(~c_r -> ~g c)
DECAY 1000003 5.209967e+01 # sstrange_l decays
1.137659e-01 2 3 1000022 # br(~s_l -> ~chi_10 s)
1.884572e-05 2 3 1000023 # br(~s_l -> ~chi_20 s)
4.947629e-02 2 3 1000025 # br(~s_l -> ~chi_30 s)
6.006867e-03 2 3 1000035 # br(~s_l -> ~chi_40 s)
2.144486e-01 2 -1000024 4 # br(~s_l -> ~chi_1- c)
1.124764e-01 2 -1000037 4 # br(~s_l -> ~chi_2- c)
5.038072e-01 2 3 1000021 # br(~s_l -> ~g s)
DECAY 2000003 5.508143e-01 # sstrange_r decays
2.768362e-04 2 3 1000022 # br(~s_r -> ~chi_10 s)
1.403730e-01 2 3 1000023 # br(~s_r -> ~chi_20 s)
5.846603e-04 2 3 1000025 # br(~s_r -> ~chi_30 s)
8.587655e-01 2 3 1000035 # br(~s_r -> ~chi_40 s)
DECAY 1000011 3.250427e+01 # selectron_l decays
2.157420e-01 2 11 1000022 # br(~e_l -> ~chi_10 e-)
1.953818e-02 2 11 1000023 # br(~e_l -> ~chi_20 e-)
8.603555e-02 2 11 1000025 # br(~e_l -> ~chi_30 e-)
7.434585e-02 2 11 1000035 # br(~e_l -> ~chi_40 e-)
3.956946e-01 2 -1000024 12 # br(~e_l -> ~chi_1- nu_e)
2.086438e-01 2 -1000037 12 # br(~e_l -> ~chi_2- nu_e)
DECAY 2000011 1.540539e+01 # selectron_r decays
2.384824e-04 2 11 1000022 # br(~e_r -> ~chi_10 e-)
1.275975e-01 2 11 1000023 # br(~e_r -> ~chi_20 e-)
5.595371e-04 2 11 1000025 # br(~e_r -> ~chi_30 e-)
8.716045e-01 2 11 1000035 # br(~e_r -> ~chi_40 e-)
DECAY 1000013 3.250427e+01 # smuon_l decays
2.157420e-01 2 13 1000022 # br(~mu_l -> ~chi_10 mu-)
1.953818e-02 2 13 1000023 # br(~mu_l -> ~chi_20 mu-)
8.603555e-02 2 13 1000025 # br(~mu_l -> ~chi_30 mu-)
7.434585e-02 2 13 1000035 # br(~mu_l -> ~chi_40 mu-)
3.956946e-01 2 -1000024 14 # br(~mu_l -> ~chi_1- nu_mu)
2.086438e-01 2 -1000037 14 # br(~mu_l -> ~chi_2- nu_mu)
DECAY 2000013 1.540539e+01 # smuon_r decays
2.384824e-04 2 13 1000022 # br(~mu_r -> ~chi_10 mu-)
1.275975e-01 2 13 1000023 # br(~mu_r -> ~chi_20 mu-)
5.595371e-04 2 13 1000025 # br(~mu_r -> ~chi_30 mu-)
8.716045e-01 2 13 1000035 # br(~mu_r -> ~chi_40 mu-)
DECAY 1000015 5.114682e+00 # stau_1 decays
2.603446e-01 2 15 1000022 # br(~tau_1 -> ~chi_10 tau-)
7.541379e-02 2 15 1000023 # br(~tau_1 -> ~chi_20 tau-)
8.441980e-02 2 15 1000025 # br(~tau_1 -> ~chi_30 tau-)
3.563511e-02 2 15 1000035 # br(~tau_1 -> ~chi_40 tau-)
4.170176e-01 2 -1000024 16 # br(~tau_1 -> ~chi_1- nu_tau)
1.271692e-01 2 -1000037 16 # br(~tau_1 -> ~chi_2- nu_tau)
DECAY 2000015 2.242979e+01 # stau_2 decays
4.175121e-02 2 15 1000022 # br(~tau_2 -> ~chi_10 tau-)
1.689956e-01 2 15 1000023 # br(~tau_2 -> ~chi_20 tau-)
7.403667e-02 2 15 1000025 # br(~tau_2 -> ~chi_30 tau-)
4.654161e-01 2 15 1000035 # br(~tau_2 -> ~chi_40 tau-)
8.054943e-02 2 -1000024 16 # br(~tau_2 -> ~chi_1- nu_tau)
1.485493e-01 2 -1000037 16 # br(~tau_2 -> ~chi_2- nu_tau)
1.016242e-02 2 -24 1000016 # br(~tau_2 -> ~nu_taul w-)
5.376317e-03 2 25 1000015 # br(~tau_2 -> ~tau_1 h)
5.162967e-03 2 23 1000015 # br(~tau_2 -> ~tau_1 z)
DECAY 1000012 3.253223e+01 # snu_el decays
2.067129e-01 2 12 1000022 # br(~nu_el -> ~chi_10 nu_e)
6.191457e-03 2 12 1000023 # br(~nu_el -> ~chi_20 nu_e)
9.455957e-02 2 12 1000025 # br(~nu_el -> ~chi_30 nu_e)
8.682663e-02 2 12 1000035 # br(~nu_el -> ~chi_40 nu_e)
4.611672e-01 2 11 1000024 # br(~nu_el -> ~chi_1+ e-)
1.445422e-01 2 11 1000037 # br(~nu_el -> ~chi_2+ e-)
DECAY 1000014 3.253223e+01 # snu_mul decays
2.067129e-01 2 14 1000022 # br(~nu_mul -> ~chi_10 nu_mu)
6.191457e-03 2 14 1000023 # br(~nu_mul -> ~chi_20 nu_mu)
9.455957e-02 2 14 1000025 # br(~nu_mul -> ~chi_30 nu_mu)
8.682663e-02 2 14 1000035 # br(~nu_mul -> ~chi_40 nu_mu)
4.611672e-01 2 13 1000024 # br(~nu_mul -> ~chi_1+ mu-)
1.445422e-01 2 13 1000037 # br(~nu_mul -> ~chi_2+ mu-)
DECAY 1000016 5.129139e+00 # snu_taul decays
2.142951e-01 2 16 1000022 # br(~nu_taul -> ~chi_10 nu_tau)
4.953687e-03 2 16 1000023 # br(~nu_taul -> ~chi_20 nu_tau)
5.728456e-02 2 16 1000025 # br(~nu_taul -> ~chi_30 nu_tau)
3.653458e-02 2 16 1000035 # br(~nu_taul -> ~chi_40 nu_tau)
5.388663e-01 2 15 1000024 # br(~nu_taul -> ~chi_1+ tau-)
1.480658e-01 2 15 1000037 # br(~nu_taul -> ~chi_2+ tau-)
DECAY 1000024 1.874508e-13 # chargino1+ decays
6.156997e-01 3 -1 2 1000022 # br(~chi_1+ -> ~chi_10 u db)
2.052037e-01 3 -11 12 1000022 # br(~chi_1+ -> ~chi_10 e+ nu_e)
1.790966e-01 3 -13 14 1000022 # br(~chi_1+ -> ~chi_10 mu+ nu_mu)
DECAY 1000037 3.745775e-01 # chargino2+ decays
3.840254e-01 2 23 1000024 # br(~chi_2+ -> ~chi_1+ z )
6.159746e-01 2 24 1000022 # br(~chi_2+ -> ~chi_10 w+)
DECAY 1000022 0.000000e+00 # neutralino1 decays
DECAY 1000023 1.142369e-03 # neutralino2 decays
2.395996e-04 2 22 1000022 # br(~chi_20 -> ~chi_10 gam)
3.238105e-02 3 -2 2 1000022 # br(~chi_20 -> ~chi_10 ub u)
4.192412e-02 3 -1 1 1000022 # br(~chi_20 -> ~chi_10 db d)
3.223674e-02 3 -4 4 1000022 # br(~chi_20 -> ~chi_10 cb c)
4.192192e-02 3 -3 3 1000022 # br(~chi_20 -> ~chi_10 sb s)
4.094287e-02 3 -5 5 1000022 # br(~chi_20 -> ~chi_10 bb b)
9.574011e-03 3 -11 11 1000022 # br(~chi_20 -> ~chi_10 e+ e-)
9.573812e-03 3 -13 13 1000022 # br(~chi_20 -> ~chi_10 mu+ mu-)
9.510830e-03 3 -15 15 1000022 # br(~chi_20 -> ~chi_10 tau+ tau-)
1.913016e-02 3 -12 12 1000022 # br(~chi_20 -> ~chi_10 nu_eb nu_e)
1.913016e-02 3 -14 14 1000022 # br(~chi_20 -> ~chi_10 nu_mub nu_mu)
1.902949e-02 3 -16 16 1000022 # br(~chi_20 -> ~chi_10 nu_taub nu_tau)
1.209853e-01 3 -2 1 1000024 # br(~chi_20 -> ~chi_1+ ub d)
1.209853e-01 3 -1000024 -1 2 # br(~chi_20 -> ~chi_1- db u)
1.206960e-01 3 -4 3 1000024 # br(~chi_20 -> ~chi_1+ cb s)
1.206960e-01 3 -1000024 -3 4 # br(~chi_20 -> ~chi_1- sb c)
4.030061e-02 3 -12 11 1000024 # br(~chi_20 -> ~chi_1+ nu_eb e-)
4.030061e-02 3 -1000024 -11 12 # br(~chi_20 -> ~chi_1- nu_e e+)
4.030009e-02 3 -14 13 1000024 # br(~chi_20 -> ~chi_1+ nu_mub mu-)
4.030009e-02 3 -1000024 -13 14 # br(~chi_20 -> ~chi_1- nu_mu mu+)
3.992060e-02 3 -16 15 1000024 # br(~chi_20 -> ~chi_1+ nu_taub tau-)
3.992060e-02 3 -1000024 -15 16 # br(~chi_20 -> ~chi_1- nu_tau tau+)
DECAY 1000025 4.733107e-01 # neutralino3 decays
5.058443e-03 2 23 1000022 # br(~chi_30 -> ~chi_10 z )
4.974708e-01 2 -24 1000024 # br(~chi_30 -> ~chi_1+ w-)
4.974708e-01 2 -1000024 24 # br(~chi_30 -> ~chi_1- w+)
DECAY 1000035 5.645215e-01 # neutralino4 decays
3.024198e-01 2 23 1000022 # br(~chi_40 -> ~chi_10 z )
2.346536e-03 2 23 1000023 # br(~chi_40 -> ~chi_20 z )
3.455560e-01 2 -24 1000024 # br(~chi_40 -> ~chi_1+ w-)
3.455560e-01 2 -1000024 24 # br(~chi_40 -> ~chi_1- w+)
4.121628e-03 2 25 1000022 # br(~chi_40 -> ~chi_10 h )
DECAY 25 3.802922e-03 # h decays
6.312413e-01 2 -5 5 # br(h -> b bb )
6.533574e-02 2 -15 15 # br(h -> tau+ tau- )
2.313091e-04 2 -13 13 # br(h -> mu+ mu- )
4.500116e-04 2 -3 3 # br(h -> s sb )
2.013893e-02 2 -4 4 # br(h -> c cb )
7.296075e-02 2 21 21 # br(h -> g g )
2.318767e-03 2 22 22 # br(h -> gam gam )
1.385037e-03 2 22 23 # br(h -> z gam )
1.840756e-01 2 -24 24 # br(h -> w+ w- )
2.186249e-02 2 23 23 # br(h -> z z )
DECAY 35 3.183919e+02 # h decays
3.756318e-01 2 -5 5 # br(h -> b bb )
6.983553e-02 2 -15 15 # br(h -> tau+ tau- )
2.469324e-04 2 -13 13 # br(h -> mu+ mu- )
2.423670e-04 2 -3 3 # br(h -> s sb )
6.109934e-09 2 -4 4 # br(h -> c cb )
6.781835e-04 2 -6 6 # br(h -> t tb )
8.946111e-07 2 21 21 # br(h -> g g )
6.652305e-09 2 22 22 # br(h -> gam gam )
1.225714e-09 2 22 23 # br(h -> z gam )
6.927654e-05 2 -24 24 # br(h -> w+ w- )
3.463240e-05 2 23 23 # br(h -> z z )
1.435246e-07 2 25 25 # br(h -> h h )
-3.938638e-25 2 36 36 # br(h -> a a )
1.607980e-16 2 23 36 # br(h -> z a )
6.136426e-19 2 -37 24 # br(h -> w+ h- )
6.136426e-19 2 -24 37 # br(h -> w- h+ )
5.905855e-02 2 -1000024 1000024 # br(h -> ~chi_1+ ~chi_1-)
3.571930e-02 2 -1000037 1000037 # br(h -> ~chi_2+ ~chi_2-)
5.915934e-02 2 -1000037 1000024 # br(h -> ~chi_1+ ~chi_2-)
5.915934e-02 2 -1000024 1000037 # br(h -> ~chi_2+ ~chi_1-)
2.747000e-02 2 1000022 1000022 # br(h -> ~chi_10 ~chi_10)
1.618247e-03 2 1000023 1000023 # br(h -> ~chi_20 ~chi_20)
2.288114e-02 2 1000025 1000025 # br(h -> ~chi_30 ~chi_30)
2.872152e-03 2 1000035 1000035 # br(h -> ~chi_40 ~chi_40)
3.765225e-02 2 1000022 1000023 # br(h -> ~chi_10 ~chi_20)
5.730624e-03 2 1000022 1000025 # br(h -> ~chi_10 ~chi_30)
2.228997e-04 2 1000022 1000035 # br(h -> ~chi_10 ~chi_40)
8.869994e-03 2 1000023 1000025 # br(h -> ~chi_20 ~chi_30)
1.009795e-02 2 1000023 1000035 # br(h -> ~chi_20 ~chi_40)
2.000727e-02 2 1000025 1000035 # br(h -> ~chi_30 ~chi_40)
3.108091e-08 2 -1000002 1000002 # br(h -> ~u_l ~u_l* )
6.139496e-09 2 -2000002 2000002 # br(h -> ~u_r ~u_r* )
3.108091e-08 2 -1000004 1000004 # br(h -> ~c_l ~c_l* )
6.139496e-09 2 -2000004 2000004 # br(h -> ~c_r ~c_r* )
2.427374e-04 2 -1000006 1000006 # br(h -> ~t_1 ~t_1* )
1.871709e-04 2 -2000006 2000006 # br(h -> ~t_2 ~t_2* )
2.217602e-03 2 -2000006 1000006 # br(h -> ~t_1 ~t_2* )
2.217602e-03 2 -1000006 2000006 # br(h -> ~t_2 ~t_1* )
4.585173e-08 2 -1000001 1000001 # br(h -> ~d_l ~d_l* )
1.662163e-09 2 -2000001 2000001 # br(h -> ~d_r ~d_r* )
4.585173e-08 2 -1000003 1000003 # br(h -> ~s_l ~s_l* )
1.662163e-09 2 -2000003 2000003 # br(h -> ~s_r ~s_r* )
6.779099e-06 2 -1000005 1000005 # br(h -> ~b_1 ~b_1* )
3.010460e-07 2 -2000005 2000005 # br(h -> ~b_2 ~b_2* )
9.708559e-02 2 -2000005 1000005 # br(h -> ~b_1 ~b_2* )
9.708559e-02 2 -1000005 2000005 # br(h -> ~b_2 ~b_1* )
6.088463e-09 2 -1000011 1000011 # br(h -> ~e_l- ~e_l+ )
3.285902e-09 2 -2000011 2000011 # br(h -> ~e_r- ~e_r+ )
6.088463e-09 2 -1000013 1000013 # br(h -> ~mu_l- ~mu_l+ )
3.285902e-09 2 -2000013 2000013 # br(h -> ~mu_r- ~mu_r+ )
2.569366e-07 2 -1000015 1000015 # br(h -> ~tau_1- ~tau_1+)
2.061827e-07 2 -2000015 2000015 # br(h -> ~tau_2- ~tau_2+)
1.848563e-03 2 -2000015 1000015 # br(h -> ~tau_1- ~tau_2+)
1.848563e-03 2 -1000015 2000015 # br(h -> ~tau_2- ~tau_1+)
2.018534e-08 2 -1000012 1000012 # br(h -> ~nu_el ~nu_el* )
2.018534e-08 2 -1000014 1000014 # br(h -> ~nu_mul ~nu_mul* )
2.018534e-08 2 -1000016 1000016 # br(h -> ~nu_taul ~nu_taul*)
DECAY 36 3.183755e+02 # a decays
3.756300e-01 2 -5 5 # br(a -> b bb )
6.983313e-02 2 -15 15 # br(a -> tau+ tau- )
2.469239e-04 2 -13 13 # br(a -> mu+ mu- )
2.423727e-04 2 -3 3 # br(a -> s sb )
6.329328e-09 2 -4 4 # br(a -> c cb )
6.965976e-04 2 -6 6 # br(a -> t tb )
2.659044e-06 2 21 21 # br(a -> g g )
1.444836e-08 2 22 22 # br(a -> gam gam )
2.960359e-09 2 22 23 # br(a -> z gam )
6.923299e-05 2 23 25 # br(a -> z h )
5.490614e-02 2 -1000024 1000024 # br(a -> ~chi_1+ ~chi_1-)
3.156795e-02 2 -1000037 1000037 # br(a -> ~chi_2+ ~chi_2-)
6.330911e-02 2 -1000037 1000024 # br(a -> ~chi_1+ ~chi_2-)
6.330911e-02 2 -1000024 1000037 # br(a -> ~chi_2+ ~chi_1-)
2.536705e-02 2 1000022 1000022 # br(a -> ~chi_10 ~chi_10)
1.829744e-03 2 1000023 1000023 # br(a -> ~chi_20 ~chi_20)
2.071274e-02 2 1000025 1000025 # br(a -> ~chi_30 ~chi_30)
3.336356e-03 2 1000035 1000035 # br(a -> ~chi_40 ~chi_40)
4.060423e-02 2 1000022 1000023 # br(a -> ~chi_10 ~chi_20)
5.020129e-03 2 1000022 1000025 # br(a -> ~chi_10 ~chi_30)
5.714729e-05 2 1000022 1000035 # br(a -> ~chi_10 ~chi_40)
1.029036e-02 2 1000023 1000025 # br(a -> ~chi_20 ~chi_30)
1.143167e-02 2 1000023 1000035 # br(a -> ~chi_20 ~chi_40)
1.877808e-02 2 1000025 1000035 # br(a -> ~chi_30 ~chi_40)
2.427005e-03 2 -2000006 1000006 # br(a -> ~t_1 ~t_2* )
2.427005e-03 2 -1000006 2000006 # br(a -> ~t_1* ~t_2 )
9.710398e-02 2 -2000005 1000005 # br(a -> ~b_1 ~b_2* )
9.710398e-02 2 -1000005 2000005 # br(a -> ~b_1* ~b_2 )
1.848618e-03 2 -2000015 1000015 # br(a -> ~tau_1- ~tau_2+)
1.848618e-03 2 -1000015 2000015 # br(a -> ~tau_1+ ~tau_2-)
DECAY 37 3.095257e+02 # h+ decays
5.714056e-04 2 -5 4 # br(h+ -> c bb )
7.183319e-02 2 -15 16 # br(h+ -> tau+ nu_tau )
2.539959e-04 2 -13 14 # br(h+ -> mu+ nu_mu )
3.656995e-06 2 -5 2 # br(h+ -> u bb )
1.195877e-05 2 -3 2 # br(h+ -> u sb )
2.459691e-04 2 -3 4 # br(h+ -> c sb )
3.577984e-01 2 -5 6 # br(h+ -> t bb )
7.122871e-05 2 24 25 # br(h+ -> w+ h )
3.775466e-17 2 24 36 # br(h+ -> w+ a )
1.846467e-05 2 1000022 1000024 # br(h+ -> ~chi_1+ ~chi_10)
8.958250e-02 2 1000023 1000024 # br(h+ -> ~chi_1+ ~chi_20)
1.021440e-01 2 1000024 1000025 # br(h+ -> ~chi_1+ ~chi_30)
9.709563e-05 2 1000024 1000035 # br(h+ -> ~chi_1+ ~chi_40)
1.171574e-01 2 1000022 1000037 # br(h+ -> ~chi_2+ ~chi_10)
1.430620e-02 2 1000023 1000037 # br(h+ -> ~chi_2+ ~chi_20)
1.917374e-04 2 1000025 1000037 # br(h+ -> ~chi_2+ ~chi_30)
3.700837e-02 2 1000035 1000037 # br(h+ -> ~chi_2+ ~chi_40)
2.523714e-08 2 -1000011 1000012 # br(h+ -> ~e_l+ ~nu_el )
2.523714e-08 2 -1000013 1000014 # br(h+ -> ~mu_l+ ~nu_mul)
2.079772e-07 2 -1000015 1000016 # br(h+ -> ~tau_1+ ~nu_taul)
3.314649e-03 2 -2000015 1000016 # br(h+ -> ~tau_2+ ~nu_taul)
7.941443e-08 2 -1000001 1000002 # br(h+ -> ~u_l ~d_l* )
7.941443e-08 2 -1000003 1000004 # br(h+ -> ~c_l ~s_l* )
1.726168e-04 2 -1000005 1000006 # br(h+ -> ~t_1 ~b_1* )
1.948817e-03 2 -2000005 2000006 # br(h+ -> ~t_2 ~b_2* )
1.984324e-01 2 -2000005 1000006 # br(h+ -> ~t_1 ~b_2* )
4.835622e-03 2 -1000005 2000006 # br(h+ -> ~t_2 ~b_1* )
###################################
## INFORMATION FOR USQMIX
###################################
BLOCK USQMIX #
1 1 1.000000e+00 #
2 2 1.000000e+00 #
4 4 1.000000e+00 #
5 5 1.000000e+00 #
3 3 9.886007e-01 # cos(theta_t)
3 6 1.505610e-01 # sin(theta_t)
6 3 -1.505610e-01 # -sin(theta_t)
6 6 9.886007e-01 # cos(theta_t)
###################################
## INFORMATION FOR DSQMIX
###################################
BLOCK DSQMIX #
1 1 1.000000e+00 #
2 2 1.000000e+00 #
4 4 1.000000e+00 #
5 5 1.000000e+00 #
3 3 9.999965e-01 # cos(theta_b)
3 6 2.631347e-03 # sin(theta_b)
6 3 -2.631347e-03 # -sin(theta_b)
6 6 9.999965e-01 # cos(theta_b)
###################################
## INFORMATION FOR SELMIX
###################################
BLOCK SELMIX #
1 1 1.000000e+00 #
2 2 1.000000e+00 #
4 4 1.000000e+00 #
5 5 1.000000e+00 #
3 3 9.999847e-01 # cos(theta_tau)
3 6 5.539834e-03 # sin(theta_tau)
6 3 -5.539834e-03 # -sin(theta_tau)
6 6 9.999847e-01 # cos(theta_tau)
###################################
## INFORMATION FOR FRALPHA
###################################
BLOCK FRALPHA #
1 -2.719383e-02 # mixing angle in the neutral higgs boson sector
###################################
## INFORMATION FOR VCKM
###################################
BLOCK VCKM #
1 1 1.000000e+00 #
2 2 1.000000e+00 #
3 3 1.000000e+00 #
###################################
## INFORMATION FOR SNUMIX
###################################
BLOCK SNUMIX #
1 1 1.000000e+00 #
2 2 1.000000e+00 #
3 3 1.000000e+00 #
###################################
## INFORMATION FOR UPMNS
###################################
BLOCK UPMNS #
1 1 1.000000e+00 #
2 2 1.000000e+00 #
3 3 1.000000e+00 #
###################################
## INFORMATION FOR TE
###################################
BLOCK TE #
1 1 0.000000e+00 # T_e(Q) DRbar
2 2 0.000000e+00 # T_mu(Q) DRbar
3 3 -6.658103e+02 # T_tau(Q) DRbar
###################################
## INFORMATION FOR TU
###################################
BLOCK TU #
1 1 0.000000e+00 # T_u(Q) DRbar
2 2 0.000000e+00 # T_c(Q) DRbar
3 3 -1.078646e+03 # T_t(Q) DRbar
###################################
## INFORMATION FOR TD
###################################
BLOCK TD #
1 1 0.000000e+00 # T_d(Q) DRbar
2 2 0.000000e+00 # T_s(Q) DRbar
3 3 1.652909e+03 # T_b(Q) DRbar
###################################
## INFORMATION FOR MSL2
###################################
BLOCK MSL2 #
1 1 6.618656e+06 # m_el
2 2 6.618656e+06 # m_mul
3 3 2.816987e+05 # m_taul
###################################
## INFORMATION FOR MSE2
###################################
BLOCK MSE2 #
1 1 1.044618e+07 # m_er
2 2 1.044618e+07 # m_mur
3 3 6.096401e+06 # m_taur
###################################
## INFORMATION FOR MSQ2
###################################
BLOCK MSQ2 #
1 1 5.970652e+06 # m_q1l
2 2 5.970652e+06 # m_q2l
3 3 2.505174e+04 # m_q3l
###################################
## INFORMATION FOR MSU2
###################################
BLOCK MSU2 #
1 1 4.207313e+06 # m_ur
2 2 4.207313e+06 # m_cr
3 3 1.237423e+06 # m_tr
###################################
## INFORMATION FOR MSD2
###################################
BLOCK MSD2 #
1 1 2.278543e+06 # m_dr
2 2 2.278543e+06 # m_sr
3 3 1.336830e+07 # m_br
xsection 8.00e+03 2212 2212 2 -1000006 1000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.17823393e-02 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000005 1000005 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.11301517e-03 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000022 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 8.83077481e-03 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000024 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.89356152e-03 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000024 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.86067353e-03 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000025 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.74606328e-03 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000037 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.95016880e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000023 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 8.99342180e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000023 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 6.84074105e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000023 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.97966875e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000037 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.50129525e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000022 1000023 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.49671090e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000037 1000023 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 3.01375305e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000005 1000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.58321690e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000024 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.43970485e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000024 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.43970485e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000024 1000023 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.34403015e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000022 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.29619280e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000037 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.48295785e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000006 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.48295785e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000024 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.14809640e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000024 1000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.10025905e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000024 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.05242170e-04 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000037 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.56747000e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000037 1000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.56747000e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000006 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.56747000e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -2000006 2000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.08909650e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 2000001 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 8.13234950e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000021 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 7.17560250e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000005 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 6.69722900e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000025 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 6.21885550e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000005 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.74048200e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000006 1000005 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.74048200e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000035 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000005 1000023 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 3.34861450e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000005 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 3.34861450e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000037 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.87024100e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 2000001 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.87024100e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000022 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.87024100e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000024 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.39186750e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000021 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.39186750e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000005 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.91349400e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000005 1000023 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.91349400e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000024 1000021 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.91349400e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000021 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.91349400e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000022 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.91349400e-05 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000002 1000005 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.56747000e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000005 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.56747000e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000002 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.56747000e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000024 2000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.56747000e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000021 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.56747000e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000005 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.56747000e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000005 1000021 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.56747000e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000002 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.56747000e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000015 1000016 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.56747000e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 2000001 2000003 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000037 2000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000002 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -2000003 2000003 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000035 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000005 1000005 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000024 1000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000015 1000015 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -2000003 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000005 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000021 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000035 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 2000002 2000003 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000005 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000037 1000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -1000016 1000016 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 1000021 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -2000003 -2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 8.00e+03 2212 2212 2 -2000001 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.78373500e-06 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000006 1000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.38286126e-02 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000005 1000005 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 7.26970429e-02 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000022 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.62131430e-02 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000024 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.82992174e-02 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000024 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.16626272e-02 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000025 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.69247280e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000037 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 3.52655632e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000023 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 3.47102000e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000037 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.58243888e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000023 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.44359808e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000023 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.38806176e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 2000001 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.38806176e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000006 2000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.16591648e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000024 1000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.66608960e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000022 1000023 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.58278512e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000021 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.55501696e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 2000001 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.27733536e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000006 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.24956720e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000037 1000023 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.19403088e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000021 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.16626272e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000005 1000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.11072640e-03 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 2000002 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.99653760e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000037 1000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 9.99653760e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000024 1000023 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 8.33044800e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000005 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 8.05276640e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000024 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 8.05276640e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000006 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 7.49740320e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000022 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 7.21972160e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000021 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 6.94204000e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000037 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 6.10899520e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000002 1000021 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 6.10899520e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000021 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.83131360e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000005 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000006 1000005 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.27595040e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000021 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.27595040e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000024 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.99826880e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000005 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.72058720e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000005 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.44290560e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000024 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 4.44290560e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000005 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 3.88754240e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000024 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 3.60986080e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000005 1000023 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 3.33217920e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000035 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 3.33217920e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000022 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 3.05449760e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000021 1000021 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 3.05449760e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000005 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000024 1000021 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000002 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.49913440e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000001 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.22145280e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000002 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.22145280e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000025 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.94377120e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000001 1000021 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.94377120e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 2000002 2000003 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.66608960e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000005 1000021 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.66608960e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000003 2000003 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.38840800e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000002 1000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.38840800e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000035 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.38840800e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000021 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.38840800e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000022 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.38840800e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000037 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.38840800e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000005 1000023 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.11072640e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000024 1000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.11072640e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000003 1000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.11072640e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 2000001 2000003 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.11072640e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000001 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.11072640e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000005 1000021 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.11072640e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000037 1000021 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 1.11072640e-04 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000024 2000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 8.33044800e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000037 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 8.33044800e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000003 1000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 8.33044800e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000001 1000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 8.33044800e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000024 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 8.33044800e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000001 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 8.33044800e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000002 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000006 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000001 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000015 1000016 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000021 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000001 1000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000001 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000016 1000016 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000002 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000037 2000006 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000002 1000005 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000021 2000003 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000022 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 5.55363200e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000016 1000015 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000023 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000035 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000001 -2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000025 1000025 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000002 1000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000001 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000001 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000004 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000005 1000035 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000003 1000021 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000023 1000023 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000005 1000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000005 1000005 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000006 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000005 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 2000001 2000004 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000003 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000037 1000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000003 -2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000001 1000021 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000005 1000037 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000004 1000021 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000003 2000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000002 1000022 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -1000002 1000024 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000001 1000005 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 1000005 2000002 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
xsection 1.30e+04 2212 2212 2 -2000001 1000001 # 10000 events, [pb], pythia8 for lo
0 0 0 0 0 0 2.77681600e-05 smodelsv1.2.2
directory
! /dev/null : to turn off the HEPMC output.
! : to select where the HEPMC file must written. It will
! therefore not be placed in the run_ directory. The
! specified path, if not absolute, will be relative to
! the Event/run_ directory of the process output.
! fifo : to have MG5aMC setup the piping of the PY8 output to
! analysis tools such as MadAnalysis5.
! fifo@ :
! Same as 'fifo', but selecting a custom path to create the
! fifo pipe. (useful to select a mounted drive that supports
! fifo). Note that the fifo file extension *must* be '.hepmc.fifo'.
! -------------------------------------------------------------------
!
HEPMCoutput:file = tag_1_pythia8_events.hepmc
!
! --------------------------------------------------------------------
! Parameters relevant only when performing MLM merging, which can be
! turned on by setting ickkw to '1' in the run_card and chosing a
! positive value for the parameter xqcut.
! For details, see section 'Jet Matching' on the left-hand menu of
! http://home.thep.lu.se/~torbjorn/pythia81html/Welcome.html
! --------------------------------------------------------------------
! If equal to -1.0, MadGraph5_aMC@NLO will set it automatically based
! on the parameter 'xqcut' of the run_card.dat
JetMatching:qCut = 1.5000000000e+02
! Use default kt-MLM to match parton level jets to those produced by the
! shower. But the other Shower-kt scheme is available too with this option.
JetMatching:doShowerKt = off
! A value of -1 means that it is automatically guessed by MadGraph.
! It is however always safer to explicitly set it.
JetMatching:nJetMax = 2
!
! --------------------------------------------------------------------
! Parameters relevant only when performing CKKW-L merging, which can
! be turned on by setting the parameter 'ptlund' *or* 'ktdurham' to
! a positive value.
! For details, see section 'CKKW-L Merging' on the left-hand menu of
! http://home.thep.lu.se/~torbjorn/pythia81html/Welcome.html
! --------------------------------------------------------------------
! Central merging scale values you want to be used.
! If equal to -1.0, then MadGraph5_aMC@NLO will set this automatically
! based on the parameter 'ktdurham' of the run_card.dat
! The following parameter was forced to be commented out by MG5aMC.
! Merging:TMS = -1.0
! This must be set manually, according to Pythia8 directives.
! An example of possible value is 'pp>LEPTONS,NEUTRINOS'
! Alternatively, from Pythia v8.223 onwards, the value 'guess' can be
! used to instruct Pythia to guess the hard process. The guess would mean
! that all particles apart from light partons will be considered as a part
! of the hard process. This guess is prone to errors if the desired hard
! process is complicated (i.e. contains light partons). The user should
! then be wary of suspicious error messages in the Pythia log file.
! The following parameter was forced to be commented out by MG5aMC.
! Merging:Process =
! A value of -1 means that it is automatically guessed by MadGraph.
! It is however always safer to explicitly set it.
! The following parameter was forced to be commented out by MG5aMC.
! Merging:nJetMax = -1
!
! For all merging schemes, decide whehter you want the merging scale
! variation computed for only the central weights or all other
! PDF and scale variation weights as well
SysCalc:fullCutVariation = off
!
! ==========================
! User customized parameters
! ==========================
!
! By default, Pythia8 generates multi-parton interaction events. This is
! often irrelevant for phenomenology and very slow. You can turn this
! feature off by uncommenting the line below if so desired.
!partonlevel:mpi = off
!
! Additional general parameters.
!
partonlevel:mpi=off
!
! Additional technical parameters set by MG5_aMC.
!
! 1.0 corresponds to HEPMC weight given in [mb]. We choose here the [pb] normalization.
HEPMCoutput:scaling=1.0000000000e+09
! Value of the merging scale below which one does not even write the HepMC event.
SysCalc:qWeed=1.0000000000e+02
! This parameter is automatically set to True by MG5aMC when doing MLM merging with PY8.
Beams:setProductionScalesFromLHEF=on
! Tell Pythia8 that an LHEF input is used.
Beams:frameType=4
! Specify one must read inputs from the MadGraph banner.
JetMatching:setMad=off
JetMatching:nQmatch=5
JetMatching:etaJetMax=1.0000000000e+03
JetMatching:coneRadius=1.0000000000e+00
! Specifiy if we are merging sample of different multiplicity.
JetMatching:merge=on
JetMatching:scheme=1
! Be more forgiving with momentum mismatches.
Check:epTolErr=1.0000000000e-02
!
! ====================
! Subrun definitions
! ====================
!
LHEFInputs:nSubruns=1
Main:subrun=0
!
! Definition of subrun 0
!
Beams:LHEF=unweighted_events.lhe.gz
]]>
0.1 && pt <= 1.0) * (0.70) +
(abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
(abs(eta) > 2.5) * (0.00)}
}
##############################
# Electron tracking efficiency
##############################
module Efficiency ElectronTrackingEfficiency {
set InputArray ParticlePropagator/electrons
set OutputArray electrons
# set EfficiencyFormula {efficiency formula as a function of eta and pt}
# tracking efficiency formula for electrons
set EfficiencyFormula { (pt <= 0.1) * (0.00) +
(abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
(abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
(abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
(abs(eta) > 2.5) * (0.00)}
}
##########################
# Muon tracking efficiency
##########################
module Efficiency MuonTrackingEfficiency {
set InputArray ParticlePropagator/muons
set OutputArray muons
# set EfficiencyFormula {efficiency formula as a function of eta and pt}
# tracking efficiency formula for muons
set EfficiencyFormula { (pt <= 0.1) * (0.00) +
(abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
(abs(eta) <= 1.5) * (pt > 1.0) * (0.99) +
(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) +
(abs(eta) > 2.5) * (0.00)}
}
########################################
# Momentum resolution for charged tracks
########################################
module MomentumSmearing ChargedHadronMomentumSmearing {
set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
set OutputArray chargedHadrons
# set ResolutionFormula {resolution formula as a function of eta and pt}
# resolution formula for charged hadrons
set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
(abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
}
###################################
# Momentum resolution for electrons
###################################
module MomentumSmearing ElectronMomentumSmearing {
set InputArray ElectronTrackingEfficiency/electrons
set OutputArray electrons
# set ResolutionFormula {resolution formula as a function of eta and energy}
# resolution formula for electrons
set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
(abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
}
###############################
# Momentum resolution for muons
###############################
module MomentumSmearing MuonMomentumSmearing {
set InputArray MuonTrackingEfficiency/muons
set OutputArray muons
# set ResolutionFormula {resolution formula as a function of eta and pt}
# resolution formula for muons
set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
(abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
}
##############
# Track merger
##############
module Merger TrackMerger {
# add InputArray InputArray
add InputArray ChargedHadronMomentumSmearing/chargedHadrons
add InputArray ElectronMomentumSmearing/electrons
add InputArray MuonMomentumSmearing/muons
set OutputArray tracks
}
#############
# ECAL
#############
module SimpleCalorimeter ECal {
set ParticleInputArray ParticlePropagator/stableParticles
set TrackInputArray TrackMerger/tracks
set TowerOutputArray ecalTowers
set EFlowTrackOutputArray eflowTracks
set EFlowTowerOutputArray eflowPhotons
set IsEcal true
set EnergyMin 0.5
set EnergySignificanceMin 2.0
set SmearTowerCenter true
set pi [expr {acos(-1)}]
# lists of the edges of each tower in eta and phi
# each list starts with the lower edge of the first tower
# the list ends with the higher edged of the last tower
# assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
set PhiBins {}
for {set i -180} {$i <= 180} {incr i} {
add PhiBins [expr {$i * $pi/180.0}]
}
# 0.02 unit in eta up to eta = 1.5 (barrel)
for {set i -85} {$i <= 86} {incr i} {
set eta [expr {$i * 0.0174}]
add EtaPhiBins $eta $PhiBins
}
# assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0
set PhiBins {}
for {set i -180} {$i <= 180} {incr i} {
add PhiBins [expr {$i * $pi/180.0}]
}
# 0.02 unit in eta up to eta = 3
for {set i 1} {$i <= 84} {incr i} {
set eta [expr { -2.958 + $i * 0.0174}]
add EtaPhiBins $eta $PhiBins
}
for {set i 1} {$i <= 84} {incr i} {
set eta [expr { 1.4964 + $i * 0.0174}]
add EtaPhiBins $eta $PhiBins
}
# take present CMS granularity for HF
# 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
set PhiBins {}
for {set i -18} {$i <= 18} {incr i} {
add PhiBins [expr {$i * $pi/18.0}]
}
foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
add EtaPhiBins $eta $PhiBins
}
add EnergyFraction {0} {0.0}
# energy fractions for e, gamma and pi0
add EnergyFraction {11} {1.0}
add EnergyFraction {22} {1.0}
add EnergyFraction {111} {1.0}
# energy fractions for muon, neutrinos and neutralinos
add EnergyFraction {12} {0.0}
add EnergyFraction {13} {0.0}
add EnergyFraction {14} {0.0}
add EnergyFraction {16} {0.0}
add EnergyFraction {1000022} {0.0}
add EnergyFraction {1000023} {0.0}
add EnergyFraction {1000025} {0.0}
add EnergyFraction {1000035} {0.0}
add EnergyFraction {1000045} {0.0}
# energy fractions for K0short and Lambda
add EnergyFraction {310} {0.3}
add EnergyFraction {3122} {0.3}
# set ResolutionFormula {resolution formula as a function of eta and energy}
# set ECalResolutionFormula {resolution formula as a function of eta and energy}
# http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
# http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
# http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
set ResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
(abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
}
#############
# HCAL
#############
module SimpleCalorimeter HCal {
set ParticleInputArray ParticlePropagator/stableParticles
set TrackInputArray ECal/eflowTracks
set TowerOutputArray hcalTowers
set EFlowTrackOutputArray eflowTracks
set EFlowTowerOutputArray eflowNeutralHadrons
set IsEcal false
set EnergyMin 1.0
set EnergySignificanceMin 2.0
set SmearTowerCenter true
set pi [expr {acos(-1)}]
# lists of the edges of each tower in eta and phi
# each list starts with the lower edge of the first tower
# the list ends with the higher edged of the last tower
# 10 degrees towers
set PhiBins {}
for {set i -18} {$i <= 18} {incr i} {
add PhiBins [expr {$i * $pi/18.0}]
}
foreach eta {-3.2 -2.5 -2.4 -2.3 -2.2 -2.1 -2 -1.9 -1.8 -1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 3.3} {
add EtaPhiBins $eta $PhiBins
}
# 20 degrees towers
set PhiBins {}
for {set i -9} {$i <= 9} {incr i} {
add PhiBins [expr {$i * $pi/9.0}]
}
foreach eta {-4.9 -4.7 -4.5 -4.3 -4.1 -3.9 -3.7 -3.5 -3.3 -3 -2.8 -2.6 2.8 3 3.2 3.5 3.7 3.9 4.1 4.3 4.5 4.7 4.9} {
add EtaPhiBins $eta $PhiBins
}
# default energy fractions {abs(PDG code)} {Fecal Fhcal}
add EnergyFraction {0} {1.0}
# energy fractions for e, gamma and pi0
add EnergyFraction {11} {0.0}
add EnergyFraction {22} {0.0}
add EnergyFraction {111} {0.0}
# energy fractions for muon, neutrinos and neutralinos
add EnergyFraction {12} {0.0}
add EnergyFraction {13} {0.0}
add EnergyFraction {14} {0.0}
add EnergyFraction {16} {0.0}
add EnergyFraction {1000022} {0.0}
add EnergyFraction {1000023} {0.0}
add EnergyFraction {1000025} {0.0}
add EnergyFraction {1000035} {0.0}
add EnergyFraction {1000045} {0.0}
# energy fractions for K0short and Lambda
add EnergyFraction {310} {0.7}
add EnergyFraction {3122} {0.7}
# http://arxiv.org/pdf/hep-ex/0004009v1
# http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
# set HCalResolutionFormula {resolution formula as a function of eta and energy}
set ResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
(abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
(abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
}
#################
# Electron filter
#################
module PdgCodeFilter ElectronFilter {
set InputArray HCal/eflowTracks
set OutputArray electrons
set Invert true
add PdgCode {11}
add PdgCode {-11}
}
######################
# ChargedHadronFilter
######################
module PdgCodeFilter ChargedHadronFilter {
set InputArray HCal/eflowTracks
set OutputArray chargedHadrons
add PdgCode {11}
add PdgCode {-11}
add PdgCode {13}
add PdgCode {-13}
}
###################################################
# Tower Merger (in case not using e-flow algorithm)
###################################################
module Merger Calorimeter {
# add InputArray InputArray
add InputArray ECal/ecalTowers
add InputArray HCal/hcalTowers
set OutputArray towers
}
####################
# Energy flow merger
####################
module Merger EFlowMerger {
# add InputArray InputArray
add InputArray HCal/eflowTracks
add InputArray ECal/eflowPhotons
add InputArray HCal/eflowNeutralHadrons
set OutputArray eflow
}
######################
# EFlowFilter
######################
module PdgCodeFilter EFlowFilter {
set InputArray EFlowMerger/eflow
set OutputArray eflow
add PdgCode {11}
add PdgCode {-11}
add PdgCode {13}
add PdgCode {-13}
}
###################
# Photon efficiency
###################
module Efficiency PhotonEfficiency {
set InputArray ECal/eflowPhotons
set OutputArray photons
# set EfficiencyFormula {efficiency formula as a function of eta and pt}
# efficiency formula for photons
set EfficiencyFormula { (pt <= 10.0) * (0.00) +
(abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
(abs(eta) > 2.5) * (0.00)}
}
##################
# Photon isolation
##################
module Isolation PhotonIsolation {
set CandidateInputArray PhotonEfficiency/photons
set IsolationInputArray EFlowFilter/eflow
set OutputArray photons
set DeltaRMax 0.5
set PTMin 0.5
set PTRatioMax 0.12
}
#####################
# Electron efficiency
#####################
module Efficiency ElectronEfficiency {
set InputArray ElectronFilter/electrons
set OutputArray electrons
# set EfficiencyFormula {efficiency formula as a function of eta and pt}
# efficiency formula for electrons
set EfficiencyFormula { (pt <= 10.0) * (0.00) +
(abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
(abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
(abs(eta) > 2.5) * (0.00)}
}
####################
# Electron isolation
####################
module Isolation ElectronIsolation {
set CandidateInputArray ElectronEfficiency/electrons
set IsolationInputArray EFlowFilter/eflow
set OutputArray electrons
set DeltaRMax 0.5
set PTMin 0.5
set PTRatioMax 0.12
}
#################
# Muon efficiency
#################
module Efficiency MuonEfficiency {
set InputArray MuonMomentumSmearing/muons
set OutputArray muons
# set EfficiencyFormula {efficiency as a function of eta and pt}
# efficiency formula for muons
set EfficiencyFormula { (pt <= 10.0) * (0.00) +
(abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
(abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) +
(abs(eta) > 2.7) * (0.00)}
}
################
# Muon isolation
################
module Isolation MuonIsolation {
set CandidateInputArray MuonEfficiency/muons
set IsolationInputArray EFlowFilter/eflow
set OutputArray muons
set DeltaRMax 0.5
set PTMin 0.5
set PTRatioMax 0.25
}
###################
# Missing ET merger
###################
module Merger MissingET {
# add InputArray InputArray
add InputArray Calorimeter/towers
set MomentumOutputArray momentum
}
##################
# Scalar HT merger
##################
module Merger ScalarHT {
# add InputArray InputArray
add InputArray UniqueObjectFinder/jets
add InputArray UniqueObjectFinder/electrons
add InputArray UniqueObjectFinder/photons
add InputArray UniqueObjectFinder/muons
set EnergyOutputArray energy
}
#####################
# Neutrino Filter
#####################
module PdgCodeFilter NeutrinoFilter {
set InputArray Delphes/stableParticles
set OutputArray filteredParticles
set PTMin 0.0
add PdgCode {12}
add PdgCode {14}
add PdgCode {16}
add PdgCode {-12}
add PdgCode {-14}
add PdgCode {-16}
}
#####################
# MC truth jet finder
#####################
module FastJetFinder GenJetFinder {
set InputArray NeutrinoFilter/filteredParticles
set OutputArray jets
# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
set JetAlgorithm 6
set ParameterR 0.6
set JetPTMin 20.0
}
#########################
# Gen Missing ET merger
########################
module Merger GenMissingET {
# add InputArray InputArray
add InputArray NeutrinoFilter/filteredParticles
set MomentumOutputArray momentum
}
############
# Jet finder
############
module FastJetFinder FastJetFinder {
set InputArray Calorimeter/towers
set OutputArray jets
# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
set JetAlgorithm 6
set ParameterR 0.6
set JetPTMin 20.0
}
##################
# Jet Energy Scale
##################
module EnergyScale JetEnergyScale {
set InputArray FastJetFinder/jets
set OutputArray jets
# scale formula for jets
set ScaleFormula { sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 ) }
}
########################
# Jet Flavor Association
########################
module JetFlavorAssociation JetFlavorAssociation {
set PartonInputArray Delphes/partons
set ParticleInputArray Delphes/allParticles
set ParticleLHEFInputArray Delphes/allParticlesLHEF
set JetInputArray JetEnergyScale/jets
set DeltaR 0.5
set PartonPTMin 1.0
set PartonEtaMax 2.5
}
###########
# b-tagging
###########
module BTagging BTagging {
set JetInputArray JetEnergyScale/jets
set BitNumber 0
# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
# PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
# gluon's PDG code has the lowest priority
# based on ATL-PHYS-PUB-2015-022
# default efficiency formula (misidentification rate)
add EfficiencyFormula {0} {0.002+7.3e-06*pt}
# efficiency formula for c-jets (misidentification rate)
add EfficiencyFormula {4} {0.20*tanh(0.02*pt)*(1/(1+0.0034*pt))}
# efficiency formula for b-jets
add EfficiencyFormula {5} {0.80*tanh(0.003*pt)*(30/(1+0.086*pt))}
}
#############
# tau-tagging
#############
module TrackCountingTauTagging TauTagging {
set ParticleInputArray Delphes/allParticles
set PartonInputArray Delphes/partons
set TrackInputArray TrackMerger/tracks
set JetInputArray JetEnergyScale/jets
set DeltaR 0.2
set DeltaRTrack 0.2
set TrackPTMin 1.0
set TauPTMin 1.0
set TauEtaMax 2.5
# instructions: {n-prongs} {eff}
# 1 - one prong efficiency
# 2 - two or more efficiency
# -1 - one prong mistag rate
# -2 - two or more mistag rate
set BitNumber 0
# taken from ATL-PHYS-PUB-2015-045 (medium working point)
add EfficiencyFormula {1} {0.70}
add EfficiencyFormula {2} {0.60}
add EfficiencyFormula {-1} {0.02}
add EfficiencyFormula {-2} {0.01}
}
#####################################################
# Find uniquely identified photons/electrons/tau/jets
#####################################################
module UniqueObjectFinder UniqueObjectFinder {
# earlier arrays take precedence over later ones
# add InputArray InputArray OutputArray
add InputArray PhotonIsolation/photons photons
add InputArray ElectronIsolation/electrons electrons
add InputArray MuonIsolation/muons muons
add InputArray JetEnergyScale/jets jets
}
##################
# ROOT tree writer
##################
# tracks, towers and eflow objects are not stored by default in the output.
# if needed (for jet constituent or other studies), uncomment the relevant
# "add Branch ..." lines.
module TreeWriter TreeWriter {
# add Branch InputArray BranchName BranchClass
add Branch Delphes/allParticles Particle GenParticle
add Branch TrackMerger/tracks Track Track
add Branch Calorimeter/towers Tower Tower
add Branch HCal/eflowTracks EFlowTrack Track
add Branch ECal/eflowPhotons EFlowPhoton Tower
add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
add Branch GenJetFinder/jets GenJet Jet
add Branch GenMissingET/momentum GenMissingET MissingET
add Branch UniqueObjectFinder/jets Jet Jet
add Branch UniqueObjectFinder/electrons Electron Electron
add Branch UniqueObjectFinder/photons Photon Photon
add Branch UniqueObjectFinder/muons Muon Muon
add Branch MissingET/momentum MissingET MissingET
add Branch ScalarHT/energy ScalarHT ScalarHT
}
]]>
# Number of Events : 100000
# Integrated weight (pb) : 0.194358438071
# Matched Integrated weight (pb) : 0.113206016975
2212 2212 6.500000e+03 6.500000e+03 0 0 247000 247000 -4 3
6.304726e-02 1.049080e-04 1.943584e-01 1
1.117550e-01 1.035728e-04 1.943584e-01 0
1.955617e-02 6.007270e-05 1.943584e-01 2
please cite 1405.0301