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