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 4.991475e+02 # m0
2 -5.568761e+02 # m12
3 4.566321e+01 # tb
5 1.713773e+03 # a
###################################
## INFORMATION FOR EXTPAR
###################################
BLOCK EXTPAR # input parameters - non-minimal models
0 4.188316e+02 # q
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS # mass spectrum
24 8.025192e+01 # w+
25 1.250265e+02 # h
35 7.334289e+03 # h
36 7.333786e+03 # a
37 7.333477e+03 # h+
5 4.823272e+00 # b-quark pole mass calculated from mb(mb)_msbar
1000001 1.755530e+03 # ~d_l
2000001 2.096707e+03 # ~d_r
1000002 1.753846e+03 # ~u_l
2000002 1.218387e+03 # ~u_r
1000003 1.755530e+03 # ~s_l
2000003 2.096707e+03 # ~s_r
1000004 1.753846e+03 # ~c_l
2000004 1.218387e+03 # ~c_r
1000005 5.510043e+02 # ~b_1
2000005 2.048243e+03 # ~b_2
1000006 5.481446e+02 # ~t_1
2000006 1.109599e+03 # ~t_2
1000011 1.275389e+03 # ~e_l
2000011 2.958254e+03 # ~e_r
1000012 1.272538e+03 # ~nu_el
1000013 1.275389e+03 # ~mu_l
2000013 2.958254e+03 # ~mu_r
1000014 1.272538e+03 # ~nu_mul
1000015 7.724033e+02 # ~tau_1
2000015 2.327526e+03 # ~tau_2
1000016 7.681417e+02 # ~nu_taul
1000021 1.697657e+03 # ~g
1000022 -3.730022e+02 # ~chi_10
1000023 3.818391e+02 # ~chi_20
1000025 5.167620e+02 # ~chi_30
1000035 -5.924809e+02 # ~chi_40
1000024 3.749539e+02 # ~chi_1+
1000037 5.928396e+02 # ~chi_2+
###################################
## INFORMATION FOR NMIX
###################################
BLOCK NMIX # neutralino mixing matrix
1 1 -3.321959e-02 # n_11
1 2 -2.723261e-01 # n_12
1 3 6.933703e-01 # n_13
1 4 6.663126e-01 # n_14
2 1 -2.544851e-01 # n_21
2 2 -5.824509e-02 # n_22
2 3 -6.868716e-01 # n_23
2 4 6.782715e-01 # n_24
3 1 9.664595e-01 # n_31
3 2 -1.527817e-02 # n_32
3 3 -1.555364e-01 # n_33
3 4 2.037918e-01 # n_34
4 1 -9.479493e-03 # n_41
4 2 9.603190e-01 # n_42
4 3 1.524906e-01 # n_43
4 4 2.333328e-01 # n_44
###################################
## INFORMATION FOR UMIX
###################################
BLOCK UMIX # chargino mixing matrix u
1 1 -2.179198e-01 # u_11
1 2 9.759667e-01 # u_12
2 1 -9.759667e-01 # u_21
2 2 -2.179198e-01 # u_22
###################################
## INFORMATION FOR VMIX
###################################
BLOCK VMIX # chargino mixing matrix v
1 1 3.333570e-01 # v_11
1 2 9.428007e-01 # v_12
2 1 9.428007e-01 # v_21
2 2 -3.333570e-01 # v_22
###################################
## INFORMATION FOR HMIX
###################################
BLOCK HMIX Q= 4.188315e+02 # drbar higgs parameters
1 3.834256e+02 # mue
2 4.518262e+01 # tb
4 5.365529e+07 # ma02
###################################
## INFORMATION FOR GAUGE
###################################
BLOCK GAUGE Q= 4.188315e+02 # the gauge couplings
1 3.606660e-01 # gprime(q) drbar
2 6.394403e-01 # g(q) drbar
3 1.074037e+00 # g3(q) drbar
###################################
## INFORMATION FOR YU
###################################
BLOCK YU Q= 4.188315e+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.349079e-01 # y_t(q) drbar
###################################
## INFORMATION FOR YD
###################################
BLOCK YD Q= 4.188315e+02 # the yukawa couplings
1 1 0.000000e+00 # y_d(q) drbar
2 2 0.000000e+00 # y_s(q) drbar
3 3 5.437502e-01 # y_b(q) drbar
###################################
## INFORMATION FOR YE
###################################
BLOCK YE Q= 4.188315e+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.997011e-01 # y_tau(q) drbar
###################################
## INFORMATION FOR MSOFT
###################################
BLOCK MSOFT Q= 4.188315e+02 # the soft susy breaking masses at the scale q
1 4.991475e+02 # m1
2 -5.568761e+02 # m2
3 1.713773e+03 # m3
21 5.198419e+07 # m^2_hd
22 -7.649628e+04 # m^2_hu
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.443298e+00 # top decays
1.000000e+00 2 5 24 # br(t -> b w+)
DECAY 1000021 7.638391e+01 # gluino decays
5.765927e-02 2 -2 2000002 # br(~g -> ~u_r ub)
5.765927e-02 2 -2000002 2 # br(~g -> ~u_r* u )
5.765927e-02 2 -4 2000004 # br(~g -> ~c_r cb)
5.765927e-02 2 -2000004 4 # br(~g -> ~c_r* c )
1.608688e-01 2 -5 1000005 # br(~g -> ~b_1 bb)
1.608688e-01 2 -1000005 5 # br(~g -> ~b_1* b )
1.557256e-01 2 -6 1000006 # br(~g -> ~t_1 tb)
1.557256e-01 2 -1000006 6 # br(~g -> ~t_1* t )
6.808707e-02 2 -6 2000006 # br(~g -> ~t_2 tb)
6.808707e-02 2 -2000006 6 # br(~g -> ~t_2* t )
DECAY 1000006 9.666431e-01 # stop1 decays
9.023853e-02 2 6 1000022 # br(~t_1 -> ~chi_10 t )
9.097615e-01 2 5 1000024 # br(~t_1 -> ~chi_1+ b )
DECAY 2000006 3.131326e+01 # stop2 decays
1.529089e-01 2 6 1000022 # br(~t_2 -> ~chi_10 t )
1.519473e-01 2 6 1000023 # br(~t_2 -> ~chi_20 t )
6.462276e-02 2 6 1000025 # br(~t_2 -> ~chi_30 t )
6.980458e-03 2 6 1000035 # br(~t_2 -> ~chi_40 t )
3.289251e-01 2 5 1000024 # br(~t_2 -> ~chi_1+ b )
1.631880e-02 2 5 1000037 # br(~t_2 -> ~chi_2+ b )
1.011417e-01 2 25 1000006 # br(~t_2 -> ~t_1 h )
6.226248e-02 2 23 1000006 # br(~t_2 -> ~t_1 z )
1.148925e-01 2 24 1000005 # br(~t_2 -> ~b_1 w+)
DECAY 1000005 1.062410e+00 # sbottom1 decays
4.544492e-01 2 5 1000022 # br(~b_1 -> ~chi_10 b )
3.643430e-01 2 5 1000023 # br(~b_1 -> ~chi_20 b )
2.550111e-03 2 5 1000025 # br(~b_1 -> ~chi_30 b )
1.786577e-01 2 -1000024 6 # br(~b_1 -> ~chi_1- t )
DECAY 2000005 3.633948e+01 # sbottom2 decays
1.327474e-01 2 5 1000022 # br(~b_2 -> ~chi_10 b )
1.322094e-01 2 5 1000023 # br(~b_2 -> ~chi_20 b )
3.191844e-02 2 5 1000025 # br(~b_2 -> ~chi_30 b )
6.435517e-03 2 5 1000035 # br(~b_2 -> ~chi_40 b )
2.602179e-01 2 -1000024 6 # br(~b_2 -> ~chi_1- t )
1.299905e-02 2 -1000037 6 # br(~b_2 -> ~chi_2- t )
3.975247e-01 2 5 1000021 # br(~b_2 -> ~g b )
4.582559e-03 2 25 1000005 # br(~b_2 -> ~b_1 h )
9.366085e-03 2 23 1000005 # br(~b_2 -> ~b_1 z )
1.192299e-02 2 -24 1000006 # br(~b_2 -> ~t_1 w-)
7.594098e-05 2 -24 2000006 # br(~b_2 -> ~t_2 w-)
DECAY 1000002 1.710586e+01 # sup_l decays
2.809633e-02 2 2 1000022 # br(~u_l -> ~chi_10 u)
4.056948e-03 2 2 1000023 # br(~u_l -> ~chi_20 u)
9.187633e-03 2 2 1000025 # br(~u_l -> ~chi_30 u)
2.870556e-01 2 2 1000035 # br(~u_l -> ~chi_40 u)
8.038992e-02 2 1 1000024 # br(~u_l -> ~chi_1+ d)
5.552448e-01 2 1 1000037 # br(~u_l -> ~chi_2+ d)
3.596872e-02 2 2 1000021 # br(~u_l -> ~g u)
DECAY 2000002 1.905185e+00 # sup_r decays
1.327213e-03 2 2 1000022 # br(~u_r -> ~chi_10 u)
7.712578e-02 2 2 1000023 # br(~u_r -> ~chi_20 u)
9.214701e-01 2 2 1000025 # br(~u_r -> ~chi_30 u)
7.693489e-05 2 2 1000035 # br(~u_r -> ~chi_40 u)
DECAY 1000001 1.704493e+01 # sdown_l decays
2.575117e-02 2 1 1000022 # br(~d_l -> ~chi_10 d)
3.915843e-05 2 1 1000023 # br(~d_l -> ~chi_20 d)
1.293275e-02 2 1 1000025 # br(~d_l -> ~chi_30 d)
2.906178e-01 2 1 1000035 # br(~d_l -> ~chi_40 d)
3.451217e-02 2 -1000024 2 # br(~d_l -> ~chi_1- u)
5.979256e-01 2 -1000037 2 # br(~d_l -> ~chi_2- u)
3.822131e-02 2 1 1000021 # br(~d_l -> ~g d)
DECAY 2000001 2.296599e+01 # sdown_r decays
5.156117e-05 2 1 1000022 # br(~d_r -> ~chi_10 d)
3.016737e-03 2 1 1000023 # br(~d_r -> ~chi_20 d)
4.111949e-02 2 1 1000025 # br(~d_r -> ~chi_30 d)
3.799592e-06 2 1 1000035 # br(~d_r -> ~chi_40 d)
9.558084e-01 2 1 1000021 # br(~d_r -> ~g d)
DECAY 1000004 1.710586e+01 # scharm_l decays
2.809633e-02 2 4 1000022 # br(~c_l -> ~chi_10 c)
4.056948e-03 2 4 1000023 # br(~c_l -> ~chi_20 c)
9.187633e-03 2 4 1000025 # br(~c_l -> ~chi_30 c)
2.870556e-01 2 4 1000035 # br(~c_l -> ~chi_40 c)
8.038992e-02 2 3 1000024 # br(~c_l -> ~chi_1+ s)
5.552448e-01 2 3 1000037 # br(~c_l -> ~chi_2+ s)
3.596872e-02 2 4 1000021 # br(~c_l -> ~g c)
DECAY 2000004 1.905185e+00 # scharm_r decays
1.327213e-03 2 4 1000022 # br(~c_r -> ~chi_10 c)
7.712578e-02 2 4 1000023 # br(~c_r -> ~chi_20 c)
9.214701e-01 2 4 1000025 # br(~c_r -> ~chi_30 c)
7.693489e-05 2 4 1000035 # br(~c_r -> ~chi_40 c)
DECAY 1000003 1.704493e+01 # sstrange_l decays
2.575117e-02 2 3 1000022 # br(~s_l -> ~chi_10 s)
3.915843e-05 2 3 1000023 # br(~s_l -> ~chi_20 s)
1.293275e-02 2 3 1000025 # br(~s_l -> ~chi_30 s)
2.906178e-01 2 3 1000035 # br(~s_l -> ~chi_40 s)
3.451217e-02 2 -1000024 4 # br(~s_l -> ~chi_1- c)
5.979256e-01 2 -1000037 4 # br(~s_l -> ~chi_2- c)
3.822131e-02 2 3 1000021 # br(~s_l -> ~g s)
DECAY 2000003 2.296599e+01 # sstrange_r decays
5.156117e-05 2 3 1000022 # br(~s_r -> ~chi_10 s)
3.016737e-03 2 3 1000023 # br(~s_r -> ~chi_20 s)
4.111949e-02 2 3 1000025 # br(~s_r -> ~chi_30 s)
3.799592e-06 2 3 1000035 # br(~s_r -> ~chi_40 s)
9.558084e-01 2 3 1000021 # br(~s_r -> ~g s)
DECAY 1000011 1.095321e+01 # selectron_l decays
3.355143e-02 2 11 1000022 # br(~e_l -> ~chi_10 e-)
1.598098e-02 2 11 1000023 # br(~e_l -> ~chi_20 e-)
9.287997e-02 2 11 1000025 # br(~e_l -> ~chi_30 e-)
2.656011e-01 2 11 1000035 # br(~e_l -> ~chi_40 e-)
3.754111e-02 2 -1000024 12 # br(~e_l -> ~chi_1- nu_e)
5.544454e-01 2 -1000037 12 # br(~e_l -> ~chi_2- nu_e)
DECAY 2000011 1.441820e+01 # selectron_r decays
1.134915e-03 2 11 1000022 # br(~e_r -> ~chi_10 e-)
6.650078e-02 2 11 1000023 # br(~e_r -> ~chi_20 e-)
9.322764e-01 2 11 1000025 # br(~e_r -> ~chi_30 e-)
8.792372e-05 2 11 1000035 # br(~e_r -> ~chi_40 e-)
DECAY 1000013 1.095321e+01 # smuon_l decays
3.355143e-02 2 13 1000022 # br(~mu_l -> ~chi_10 mu-)
1.598098e-02 2 13 1000023 # br(~mu_l -> ~chi_20 mu-)
9.287997e-02 2 13 1000025 # br(~mu_l -> ~chi_30 mu-)
2.656011e-01 2 13 1000035 # br(~mu_l -> ~chi_40 mu-)
3.754111e-02 2 -1000024 14 # br(~mu_l -> ~chi_1- nu_mu)
5.544454e-01 2 -1000037 14 # br(~mu_l -> ~chi_2- nu_mu)
DECAY 2000013 1.441820e+01 # smuon_r decays
1.134915e-03 2 13 1000022 # br(~mu_r -> ~chi_10 mu-)
6.650078e-02 2 13 1000023 # br(~mu_r -> ~chi_20 mu-)
9.322764e-01 2 13 1000025 # br(~mu_r -> ~chi_30 mu-)
8.792372e-05 2 13 1000035 # br(~mu_r -> ~chi_40 mu-)
DECAY 1000015 3.560180e+00 # stau_1 decays
2.403845e-01 2 15 1000022 # br(~tau_1 -> ~chi_10 tau-)
2.069391e-01 2 15 1000023 # br(~tau_1 -> ~chi_20 tau-)
8.022783e-02 2 15 1000025 # br(~tau_1 -> ~chi_30 tau-)
1.384847e-01 2 15 1000035 # br(~tau_1 -> ~chi_40 tau-)
5.062357e-02 2 -1000024 16 # br(~tau_1 -> ~chi_1- nu_tau)
2.833403e-01 2 -1000037 16 # br(~tau_1 -> ~chi_2- nu_tau)
DECAY 2000015 2.531218e+01 # stau_2 decays
1.332140e-01 2 15 1000022 # br(~tau_2 -> ~chi_10 tau-)
1.594956e-01 2 15 1000023 # br(~tau_2 -> ~chi_20 tau-)
4.085137e-01 2 15 1000025 # br(~tau_2 -> ~chi_30 tau-)
6.503158e-03 2 15 1000035 # br(~tau_2 -> ~chi_40 tau-)
2.629753e-01 2 -1000024 16 # br(~tau_2 -> ~chi_1- nu_tau)
1.320356e-02 2 -1000037 16 # br(~tau_2 -> ~chi_2- nu_tau)
7.995331e-03 2 -24 1000016 # br(~tau_2 -> ~nu_taul w-)
4.040751e-03 2 25 1000015 # br(~tau_2 -> ~tau_1 h)
4.058585e-03 2 23 1000015 # br(~tau_2 -> ~tau_1 z)
DECAY 1000012 1.100335e+01 # snu_el decays
2.527423e-02 2 12 1000022 # br(~nu_el -> ~chi_10 nu_e)
2.833490e-03 2 12 1000023 # br(~nu_el -> ~chi_20 nu_e)
1.030154e-01 2 12 1000025 # br(~nu_el -> ~chi_30 nu_e)
2.690751e-01 2 12 1000035 # br(~nu_el -> ~chi_40 nu_e)
8.717837e-02 2 11 1000024 # br(~nu_el -> ~chi_1+ e-)
5.126234e-01 2 11 1000037 # br(~nu_el -> ~chi_2+ e-)
DECAY 1000014 1.100335e+01 # snu_mul decays
2.527423e-02 2 14 1000022 # br(~nu_mul -> ~chi_10 nu_mu)
2.833490e-03 2 14 1000023 # br(~nu_mul -> ~chi_20 nu_mu)
1.030154e-01 2 14 1000025 # br(~nu_mul -> ~chi_30 nu_mu)
2.690751e-01 2 14 1000035 # br(~nu_mul -> ~chi_40 nu_mu)
8.717837e-02 2 13 1000024 # br(~nu_mul -> ~chi_1+ mu-)
5.126234e-01 2 13 1000037 # br(~nu_mul -> ~chi_2+ mu-)
DECAY 1000016 3.583462e+00 # snu_taul decays
3.274286e-02 2 16 1000022 # br(~nu_taul -> ~chi_10 nu_tau)
3.595324e-03 2 16 1000023 # br(~nu_taul -> ~chi_20 nu_tau)
8.204686e-02 2 16 1000025 # br(~nu_taul -> ~chi_30 nu_tau)
1.333989e-01 2 16 1000035 # br(~nu_taul -> ~chi_40 nu_tau)
4.901808e-01 2 15 1000024 # br(~nu_taul -> ~chi_1+ tau-)
2.580353e-01 2 15 1000037 # br(~nu_taul -> ~chi_2+ tau-)
DECAY 1000024 2.239658e-11 # chargino1+ decays
5.870364e-01 3 -1 2 1000022 # br(~chi_1+ -> ~chi_10 u db)
2.407108e-02 3 -3 4 1000022 # br(~chi_1+ -> ~chi_10 c sb)
1.957357e-01 3 -11 12 1000022 # br(~chi_1+ -> ~chi_10 e+ nu_e)
1.929101e-01 3 -13 14 1000022 # br(~chi_1+ -> ~chi_10 mu+ nu_mu)
2.467857e-04 3 -15 16 1000022 # br(~chi_1+ -> ~chi_10 tau+ nu_tau)
DECAY 1000037 3.239862e+00 # chargino2+ decays
6.704224e-02 2 -5 1000006 # br(~chi_2+ -> ~t_1 bb)
2.499438e-01 2 23 1000024 # br(~chi_2+ -> ~chi_1+ z )
2.624993e-01 2 24 1000022 # br(~chi_2+ -> ~chi_10 w+)
2.489089e-01 2 24 1000023 # br(~chi_2+ -> ~chi_20 w+)
1.716057e-01 2 25 1000024 # br(~chi_2+ -> ~chi_1+ h )
DECAY 1000022 0.000000e+00 # neutralino1 decays
DECAY 1000023 6.819462e-08 # neutralino2 decays
1.858603e-02 2 22 1000022 # br(~chi_20 -> ~chi_10 gam)
7.887979e-02 3 -2 2 1000022 # br(~chi_20 -> ~chi_10 ub u)
1.021606e-01 3 -1 1 1000022 # br(~chi_20 -> ~chi_10 db d)
6.642187e-02 3 -4 4 1000022 # br(~chi_20 -> ~chi_10 cb c)
1.017784e-01 3 -3 3 1000022 # br(~chi_20 -> ~chi_10 sb s)
2.334007e-02 3 -11 11 1000022 # br(~chi_20 -> ~chi_10 e+ e-)
2.332353e-02 3 -13 13 1000022 # br(~chi_20 -> ~chi_10 mu+ mu-)
1.825135e-02 3 -15 15 1000022 # br(~chi_20 -> ~chi_10 tau+ tau-)
4.664169e-02 3 -12 12 1000022 # br(~chi_20 -> ~chi_10 nu_eb nu_e)
4.664169e-02 3 -14 14 1000022 # br(~chi_20 -> ~chi_10 nu_mub nu_mu)
4.662040e-02 3 -16 16 1000022 # br(~chi_20 -> ~chi_10 nu_taub nu_tau)
7.843947e-02 3 -2 1 1000024 # br(~chi_20 -> ~chi_1+ ub d)
7.843947e-02 3 -1000024 -1 2 # br(~chi_20 -> ~chi_1- db u)
6.399146e-02 3 -4 3 1000024 # br(~chi_20 -> ~chi_1+ cb s)
6.399146e-02 3 -1000024 -3 4 # br(~chi_20 -> ~chi_1- sb c)
2.613564e-02 3 -12 11 1000024 # br(~chi_20 -> ~chi_1+ nu_eb e-)
2.613564e-02 3 -1000024 -11 12 # br(~chi_20 -> ~chi_1- nu_e e+)
2.610499e-02 3 -14 13 1000024 # br(~chi_20 -> ~chi_1+ nu_mub mu-)
2.610499e-02 3 -1000024 -13 14 # br(~chi_20 -> ~chi_1- nu_mu mu+)
1.900573e-02 3 -16 15 1000024 # br(~chi_20 -> ~chi_1+ nu_taub tau-)
1.900573e-02 3 -1000024 -15 16 # br(~chi_20 -> ~chi_1- nu_tau tau+)
DECAY 1000025 5.468999e-01 # neutralino3 decays
2.490025e-01 2 23 1000022 # br(~chi_30 -> ~chi_10 z )
1.041787e-02 2 23 1000023 # br(~chi_30 -> ~chi_20 z )
3.033543e-01 2 -24 1000024 # br(~chi_30 -> ~chi_1+ w-)
3.033543e-01 2 -1000024 24 # br(~chi_30 -> ~chi_1- w+)
1.477112e-03 2 25 1000022 # br(~chi_30 -> ~chi_10 h )
1.323939e-01 2 25 1000023 # br(~chi_30 -> ~chi_20 h )
DECAY 1000035 3.216483e+00 # neutralino4 decays
1.369598e-02 2 23 1000022 # br(~chi_40 -> ~chi_10 z )
2.138265e-01 2 23 1000023 # br(~chi_40 -> ~chi_20 z )
2.763433e-01 2 -24 1000024 # br(~chi_40 -> ~chi_1+ w-)
2.763433e-01 2 -1000024 24 # br(~chi_40 -> ~chi_1- w+)
1.564895e-01 2 25 1000022 # br(~chi_40 -> ~chi_10 h )
6.793065e-03 2 25 1000023 # br(~chi_40 -> ~chi_20 h )
2.825421e-02 2 -5 1000005 # br(~chi_40 -> ~b_1 bb)
2.825421e-02 2 -1000005 5 # br(~chi_40 -> ~b_1* b )
DECAY 25 4.213767e-03 # h decays
6.185985e-01 2 -5 5 # br(h -> b bb )
6.027044e-02 2 -15 15 # br(h -> tau+ tau- )
2.133686e-04 2 -13 13 # br(h -> mu+ mu- )
4.138539e-04 2 -3 3 # br(h -> s sb )
1.839168e-02 2 -4 4 # br(h -> c cb )
6.994984e-02 2 21 21 # br(h -> g g )
2.217115e-03 2 22 22 # br(h -> gam gam )
1.478479e-03 2 22 23 # br(h -> z gam )
2.035484e-01 2 -24 24 # br(h -> w+ w- )
2.491828e-02 2 23 23 # br(h -> z z )
DECAY 35 2.881034e+02 # h decays
5.248890e-01 2 -5 5 # br(h -> b bb )
1.076897e-01 2 -15 15 # br(h -> tau+ tau- )
3.807816e-04 2 -13 13 # br(h -> mu+ mu- )
3.784980e-04 2 -3 3 # br(h -> s sb )
4.000197e-09 2 -4 4 # br(h -> c cb )
4.439045e-04 2 -6 6 # br(h -> t tb )
2.025319e-06 2 21 21 # br(h -> g g )
8.299497e-09 2 22 22 # br(h -> gam gam )
2.042808e-09 2 22 23 # br(h -> z gam )
1.903742e-05 2 -24 24 # br(h -> w+ w- )
9.516721e-06 2 23 23 # br(h -> z z )
3.029446e-08 2 25 25 # br(h -> h h )
-5.914165e-25 2 36 36 # br(h -> a a )
1.210942e-16 2 23 36 # br(h -> z a )
3.231447e-15 2 -37 24 # br(h -> w+ h- )
3.231447e-15 2 -24 37 # br(h -> w- h+ )
2.303778e-02 2 -1000024 1000024 # br(h -> ~chi_1+ ~chi_1-)
9.342894e-03 2 -1000037 1000037 # br(h -> ~chi_2+ ~chi_2-)
8.946608e-02 2 -1000037 1000024 # br(h -> ~chi_1+ ~chi_2-)
8.946608e-02 2 -1000024 1000037 # br(h -> ~chi_2+ ~chi_1-)
6.870126e-03 2 1000022 1000022 # br(h -> ~chi_10 ~chi_10)
5.973086e-04 2 1000023 1000023 # br(h -> ~chi_20 ~chi_20)
1.372235e-03 2 1000025 1000025 # br(h -> ~chi_30 ~chi_30)
4.786492e-03 2 1000035 1000035 # br(h -> ~chi_40 ~chi_40)
5.509327e-03 2 1000022 1000023 # br(h -> ~chi_10 ~chi_20)
1.244733e-02 2 1000022 1000025 # br(h -> ~chi_10 ~chi_30)
4.335305e-02 2 1000022 1000035 # br(h -> ~chi_10 ~chi_40)
1.278446e-02 2 1000023 1000025 # br(h -> ~chi_20 ~chi_30)
4.310799e-02 2 1000023 1000035 # br(h -> ~chi_20 ~chi_40)
5.644602e-03 2 1000025 1000035 # br(h -> ~chi_30 ~chi_40)
2.683537e-08 2 -1000002 1000002 # br(h -> ~u_l ~u_l* )
5.325150e-09 2 -2000002 2000002 # br(h -> ~u_r ~u_r* )
2.683537e-08 2 -1000004 1000004 # br(h -> ~c_l ~c_l* )
5.325150e-09 2 -2000004 2000004 # br(h -> ~c_r ~c_r* )
1.037839e-04 2 -1000006 1000006 # br(h -> ~t_1 ~t_1* )
7.216360e-05 2 -2000006 2000006 # br(h -> ~t_2 ~t_2* )
1.706638e-03 2 -2000006 1000006 # br(h -> ~t_1 ~t_2* )
1.706638e-03 2 -1000006 2000006 # br(h -> ~t_2 ~t_1* )
3.959874e-08 2 -1000001 1000001 # br(h -> ~d_l ~d_l* )
1.158003e-09 2 -2000001 2000001 # br(h -> ~d_r ~d_r* )
3.959874e-08 2 -1000003 1000003 # br(h -> ~s_l ~s_l* )
1.158003e-09 2 -2000003 2000003 # br(h -> ~s_r ~s_r* )
8.049235e-06 2 -1000005 1000005 # br(h -> ~b_1 ~b_1* )
1.102549e-07 2 -2000005 2000005 # br(h -> ~b_2 ~b_2* )
5.995349e-03 2 -2000005 1000005 # br(h -> ~b_1 ~b_2* )
5.995349e-03 2 -1000005 2000005 # br(h -> ~b_2 ~b_1* )
5.885905e-09 2 -1000011 1000011 # br(h -> ~e_l- ~e_l+ )
2.502419e-09 2 -2000011 2000011 # br(h -> ~e_r- ~e_r+ )
5.885905e-09 2 -1000013 1000013 # br(h -> ~mu_l- ~mu_l+ )
2.502419e-09 2 -2000013 2000013 # br(h -> ~mu_r- ~mu_r+ )
1.941791e-07 2 -1000015 1000015 # br(h -> ~tau_1- ~tau_1+)
2.140061e-07 2 -2000015 2000015 # br(h -> ~tau_2- ~tau_2+)
1.406505e-03 2 -2000015 1000015 # br(h -> ~tau_1- ~tau_2+)
1.406505e-03 2 -1000015 2000015 # br(h -> ~tau_2- ~tau_1+)
1.952994e-08 2 -1000012 1000012 # br(h -> ~nu_el ~nu_el* )
1.952994e-08 2 -1000014 1000014 # br(h -> ~nu_mul ~nu_mul* )
1.952994e-08 2 -1000016 1000016 # br(h -> ~nu_taul ~nu_taul*)
DECAY 36 2.880849e+02 # a decays
5.248946e-01 2 -5 5 # br(a -> b bb )
1.076883e-01 2 -15 15 # br(a -> tau+ tau- )
3.807764e-04 2 -13 13 # br(a -> mu+ mu- )
3.785082e-04 2 -3 3 # br(a -> s sb )
4.114601e-09 2 -4 4 # br(a -> c cb )
4.528646e-04 2 -6 6 # br(a -> t tb )
4.107392e-06 2 21 21 # br(a -> g g )
1.565458e-08 2 22 22 # br(a -> gam gam )
4.235270e-09 2 22 23 # br(a -> z gam )
1.902296e-05 2 23 25 # br(a -> z h )
2.202057e-02 2 -1000024 1000024 # br(a -> ~chi_1+ ~chi_1-)
8.343360e-03 2 -1000037 1000037 # br(a -> ~chi_2+ ~chi_2-)
9.047338e-02 2 -1000037 1000024 # br(a -> ~chi_1+ ~chi_2-)
9.047338e-02 2 -1000024 1000037 # br(a -> ~chi_2+ ~chi_1-)
6.378084e-03 2 1000022 1000022 # br(a -> ~chi_10 ~chi_10)
6.587471e-04 2 1000023 1000023 # br(a -> ~chi_20 ~chi_20)
1.571383e-03 2 1000025 1000025 # br(a -> ~chi_30 ~chi_30)
4.294465e-03 2 1000035 1000035 # br(a -> ~chi_40 ~chi_40)
5.707571e-03 2 1000022 1000023 # br(a -> ~chi_10 ~chi_20)
1.100683e-02 2 1000022 1000025 # br(a -> ~chi_10 ~chi_30)
4.062299e-02 2 1000022 1000035 # br(a -> ~chi_10 ~chi_40)
1.414045e-02 2 1000023 1000025 # br(a -> ~chi_20 ~chi_30)
4.642505e-02 2 1000023 1000035 # br(a -> ~chi_20 ~chi_40)
5.667846e-03 2 1000025 1000035 # br(a -> ~chi_30 ~chi_40)
1.792101e-03 2 -2000006 1000006 # br(a -> ~t_1 ~t_2* )
1.792101e-03 2 -1000006 2000006 # br(a -> ~t_1* ~t_2 )
6.000081e-03 2 -2000005 1000005 # br(a -> ~b_1 ~b_2* )
6.000081e-03 2 -1000005 2000005 # br(a -> ~b_1* ~b_2 )
1.406686e-03 2 -2000015 1000015 # br(a -> ~tau_1- ~tau_2+)
1.406686e-03 2 -1000015 2000015 # br(a -> ~tau_1+ ~tau_2-)
DECAY 37 2.689170e+02 # h+ decays
7.827961e-04 2 -5 4 # br(h+ -> c bb )
1.153592e-01 2 -15 16 # br(h+ -> tau+ nu_tau )
4.079003e-04 2 -13 14 # br(h+ -> mu+ nu_mu )
5.009894e-06 2 -5 2 # br(h+ -> u bb )
1.944474e-05 2 -3 2 # br(h+ -> u sb )
3.999351e-04 2 -3 4 # br(h+ -> c sb )
4.896592e-01 2 -5 6 # br(h+ -> t bb )
2.037844e-05 2 24 25 # br(h+ -> w+ h )
5.522770e-04 2 1000022 1000024 # br(h+ -> ~chi_1+ ~chi_10)
1.846933e-02 2 1000023 1000024 # br(h+ -> ~chi_1+ ~chi_20)
2.260524e-02 2 1000024 1000025 # br(h+ -> ~chi_1+ ~chi_30)
1.084805e-01 2 1000024 1000035 # br(h+ -> ~chi_1+ ~chi_40)
1.177582e-01 2 1000022 1000037 # br(h+ -> ~chi_2+ ~chi_10)
9.272348e-02 2 1000023 1000037 # br(h+ -> ~chi_2+ ~chi_20)
1.188737e-02 2 1000025 1000037 # br(h+ -> ~chi_2+ ~chi_30)
6.046486e-07 2 1000035 1000037 # br(h+ -> ~chi_2+ ~chi_40)
2.533741e-08 2 -1000011 1000012 # br(h+ -> ~e_l+ ~nu_el )
2.533741e-08 2 -1000013 1000014 # br(h+ -> ~mu_l+ ~nu_mul)
1.792864e-07 2 -1000015 1000016 # br(h+ -> ~tau_1+ ~nu_taul)
2.933851e-03 2 -2000015 1000016 # br(h+ -> ~tau_2+ ~nu_taul)
7.117781e-08 2 -1000001 1000002 # br(h+ -> ~u_l ~d_l* )
7.117781e-08 2 -1000003 1000004 # br(h+ -> ~c_l ~s_l* )
9.221608e-05 2 -1000005 1000006 # br(h+ -> ~t_1 ~b_1* )
4.630592e-04 2 -2000005 2000006 # br(h+ -> ~t_2 ~b_2* )
1.357703e-02 2 -2000005 1000006 # br(h+ -> ~t_1 ~b_2* )
3.802551e-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.937670e-01 # cos(theta_t)
3 6 1.114774e-01 # sin(theta_t)
6 3 -1.114774e-01 # -sin(theta_t)
6 6 9.937670e-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.999451e-01 # cos(theta_b)
3 6 1.047907e-02 # sin(theta_b)
6 3 -1.047907e-02 # -sin(theta_b)
6 6 9.999451e-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.999820e-01 # cos(theta_tau)
3 6 6.008300e-03 # sin(theta_tau)
6 3 -6.008300e-03 # -sin(theta_tau)
6 6 9.999820e-01 # cos(theta_tau)
###################################
## INFORMATION FOR FRALPHA
###################################
BLOCK FRALPHA #
1 -2.192289e-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 -5.068526e+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 -7.857336e+02 # 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 2.983050e+02 # T_b(Q) DRbar
###################################
## INFORMATION FOR MSL2
###################################
BLOCK MSL2 #
1 1 1.554608e+06 # m_el
2 2 1.554608e+06 # m_mul
3 3 2.480042e+05 # m_taul
###################################
## INFORMATION FOR MSE2
###################################
BLOCK MSE2 #
1 1 8.872243e+06 # m_er
2 2 8.872243e+06 # m_mur
3 3 5.003845e+06 # m_taur
###################################
## INFORMATION FOR MSQ2
###################################
BLOCK MSQ2 #
1 1 3.190670e+06 # m_q1l
2 2 3.190670e+06 # m_q2l
3 3 2.144333e+04 # m_q3l
###################################
## INFORMATION FOR MSU2
###################################
BLOCK MSU2 #
1 1 1.516040e+06 # m_ur
2 2 1.516040e+06 # m_cr
3 3 1.232435e+06 # m_tr
###################################
## INFORMATION FOR MSD2
###################################
BLOCK MSD2 #
1 1 4.547427e+06 # m_dr
2 2 4.547427e+06 # m_sr
3 3 3.370756e+06 # m_br
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.357140504752
# Matched Integrated weight (pb) : 0.212577168425
2212 2212 6.500000e+03 6.500000e+03 0 0 247000 247000 -4 3
1.125629e-01 2.362710e-04 3.571405e-01 1
2.108060e-01 1.762719e-04 3.571405e-01 0
3.377163e-02 1.070235e-04 3.571405e-01 2
please cite 1405.0301