Matching with MG5 & Pythia8
Hi,
I used MG5_aMC_v2.6.2 to generate the following process
import model MSSM_SLHA2
generate p p > x1+ n2 @1
add process p p > x1- n2 @1
add process p p > x1+ n2 j @2
add process p p > x1- n2 j @2
add process p p > x1+ n2 j j @3
add process p p > x1- n2 j j @3
output x1n2_test
with default "param_card.dat", "run_card.dat" and "pythia8_card.dat".
It can run without error. But I looked at the pythia8 log file, it seems that all the decay branch ratio are set to be 0. And there is no lepton in "delphes_
Can you please tell me if there is some problem, or, I am doing it incorrectly ?
Thanks,
Yang
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- Valentin Hirschi Edit question
- Last query:
- Last reply:
Revision history for this message
|
#1 |
<LesHouchesEvents version="3.0">
<header>
<!--
#******
# *
# MadGraph5_aMC@NLO *
# *
# Going Beyond *
# *
# http://
# http://
# http://
# *
# The MadGraph5_aMC@NLO team *
# *
#......
# *
# This file contains all the information necessary to reproduce *
# the events generated: *
# *
# 1. software version *
# 2. proc_card : code generation info including model *
# 3. param_card : model primary parameters in the LH format *
# 4. run_card : running parameters (collider and cuts) *
# 5. pythia_card : present only if pythia has been run *
# 6. pgs_card : present only if pgs has been run *
# 7. delphes_cards : present only if delphes has been run *
# *
# *
#******
-->
<MGVersion>
2.6.2
</MGVersion>
<MG5ProcCard>
<![CDATA[
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.2 2018-04-29 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set default_
set group_subprocesses Auto
set ignore_
set loop_optimized_
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model MSSM_SLHA2
generate p p > x1+ n2 @1
add process p p > x1- n2 @1
add process p p > x1+ n2 j @2
add process p p > x1- n2 j @2
add process p p > x1+ n2 j j @3
add process p p > x1- n2 j j @3
output x1n2_test
]]>
</MG5ProcCard>
<MGProcCard>
<![CDATA[
#******
# MadGraph/MadEvent *
# http://
# *
# proc_card.dat *
#******
# *
# This Files is generated by MADGRAPH 5 *
# *
# WARNING: This Files is generated for MADEVENT (compatibility issue)*
# This files is NOT a valid MG4 proc_card.dat *
# Running this in MG4 will NEVER reproduce the result of MG5*
# *
#******
#******
# Process(es) requested : mg2 input *
#******
# Begin PROCESS # This is TAG. Do not modify this line
p p > x1+ n2 @1 #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
]]>
</MGProcCard>
<MGRunCard>
<![CDATA[
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
* update full_run_card *
#******
#
#******
# Running parameters
#******
#
#******
# Tag name for the run (one word) *
#******
tag_1 = run_tag ! name of the run
#******
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
# If you want to run Pythia, avoid more than 50k events in a run. *
#******
10000 = nevents ! Number of unweighted events requested
21 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#******
1 = lpp1 ! beam 1 type
1 = lpp2 ! beam 2 type
6500.0 = ebeam1 ! beam 1 total energy in GeV
6500.0 = ebeam2 ! beam 2 total energy in GeV
# To see polarised beam options: type "update beam_pol"
#******
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#******
nn23lo1 = pdlabel ! PDF set
230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# To see heavy ion options: type "update ion_pdf"
#******
# Renormalization and factorization scales *
#******
False = fixed_ren_scale ! if .true. use fixed ren scale
False = fixed_fac_scale ! if .true. use fixed fac scale
91.188 = scale ! fixed ren scale
91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
-1 = dynamical_
1.0 = scalefact ! scale factor for event-by-event scales
#******
# Type and output format
#******
False = gridpack !True = setting up the grid pack
-1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
3.0 = lhe_version ! Change the way clustering information pass to shower.
True = clusinfo ! include clustering tag in output
average = event_norm ! average/sum. Normalization of the weight in the LHEF
#******
# Matching parameter (MLM only)
#******
1 = ickkw ! 0 no matching, 1 MLM
1.0 = alpsfact ! scale factor for QCD emission vx
False = chcluster ! cluster only according to channel diag
4 = asrwgtflavor ! highest quark flavor for a_s reweight
False = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0
50.0 = 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.
#******
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*
#******
15.0 = bwcutoff ! (M+/-bwcutoff*
#******
# Apply pt/E/eta/
# (note that etmiss/
#******
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
0.0 = ptb ! minimum pt for the b
10.0 = pta ! minimum pt for the photons
10.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) *
#******
5.0 = etaj ! max rap for the jets
-1.0 = etab ! max rap for the b
2.5 = etaa ! max rap for the photons
2.5 = 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.0 = drjj ! min distance between jets
0.0 = drbb ! min distance between b's
0.4 = drll ! min distance between leptons
0.4 = draa ! min distance between gammas
0.0 = drbj ! min distance between b and jet
0.4 = draj ! min distance between gamma and jet
0.0 = drjl ! min distance between jet and lepton
0.0 = drab ! min distance between gamma and b
0.0 = 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_
#******
# 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, 1000001, 1000002, 1000003, 1000004, 1000005, 1000006, 1000021, 2000001, 2000002, 2000003, 2000004, 2000005, 2000006 = pdgs_for_
#******
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#******
4 = maxjetflavor ! Maximum jet pdg code
#******
#
#******
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#******
True = use_syst ! Enable systematics studies
#
systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_
# Syscalc is deprecated but to see the associate options type'update syscalc'
]]>
</MGRunCard>
<slha>
#######
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
#######
#######
## INFORMATION FOR DSQMIX
#######
BLOCK DSQMIX #
1 1 1.000000e+00 # rrd1x1
2 2 1.000000e+00 # rrd2x2
3 3 9.387379e-01 # rrd3x3
3 6 3.446319e-01 # rrd3x6
4 4 1.000000e+00 # rrd4x4
5 5 1.000000e+00 # rrd5x5
6 3 -3.446319e-01 # rrd6x3
6 6 9.387379e-01 # rrd6x6
#######
## INFORMATION FOR FRALPHA
#######
BLOCK FRALPHA #
1 -1.138252e-01 # alp
#######
## INFORMATION FOR HMIX
#######
BLOCK HMIX #
1 3.576810e+02 # rmuh
2 9.748624e+00 # tb
4 1.664391e+05 # ma2
#######
## INFORMATION FOR MASS
#######
BLOCK MASS #
5 4.889917e+00 # mb
6 1.750000e+02 # mt
15 1.777000e+00 # mta
23 9.118760e+01 # mz
24 7.982901e+01 # mw
25 1.108991e+02 # mh01
35 3.999601e+02 # mh02
36 3.995839e+02 # ma0
37 4.078790e+02 # mh
1000001 5.684411e+02 # set of param :1*msd1, 1*msd2
1000002 5.611190e+02 # set of param :1*msu1, 1*msu2
1000005 5.130652e+02 # msd3
1000006 3.996685e+02 # msu3
1000011 2.029157e+02 # set of param :1*msl1, 1*msl2
1000012 1.852583e+02 # set of param :1*msn1, 1*msn2
1000015 1.344909e+02 # msl3
1000016 1.847085e+02 # msn3
1000021 6.077137e+02 # mgo
1000022 9.668807e+01 # mneu1
1000023 1.810882e+02 # mneu2
1000024 1.816965e+02 # mch1
1000025 -3.637560e+02 # mneu3
1000035 3.817294e+02 # mneu4
1000037 3.799393e+02 # mch2
2000001 5.452285e+02 # set of param :1*msd4, 1*msd5
2000002 5.492593e+02 # set of param :1*msu4, 1*msu5
2000005 5.437267e+02 # msd6
2000006 5.857858e+02 # msu6
2000011 1.441028e+02 # set of param :1*msl4, 1*msl5
2000015 2.068678e+02 # msl6
1 0.000000e+00 # d : 0.0
2 0.000000e+00 # u : 0.0
3 0.000000e+00 # s : 0.0
4 0.000000e+00 # c : 0.0
11 0.000000e+00 # e- : 0.0
12 0.000000e+00 # ve : 0.0
13 0.000000e+00 # mu- : 0.0
14 0.000000e+00 # vm : 0.0
16 0.000000e+00 # vt : 0.0
21 0.000000e+00 # g : 0.0
22 0.000000e+00 # a : 0.0
1000014 1.852583e+02 # svm : msn1
1000013 2.029157e+02 # mul- : msl1
2000013 1.441028e+02 # mur- : msl4
1000004 5.611190e+02 # cl : msu1
2000004 5.492593e+02 # cr : msu4
1000003 5.684411e+02 # sl : msd1
2000003 5.452285e+02 # sr : msd4
#######
## INFORMATION FOR MSD2
#######
BLOCK MSD2 #
1 1 2.736847e+05 # set of param :1*rmd21x1, 1*rmd22x2
2 2 2.736847e+05 # mg5 will not use this value use instead 1*mdl_rmd21x1
3 3 2.702620e+05 # rmd23x3
#######
## INFORMATION FOR MSE2
#######
BLOCK MSE2 #
1 1 1.863063e+04 # set of param :1*rme21x1, 1*rme22x2
2 2 1.863063e+04 # mg5 will not use this value use instead 1*mdl_rme21x1
3 3 1.796764e+04 # rme23x3
#######
## INFORMATION FOR MSL2
#######
BLOCK MSL2 #
1 1 3.815567e+04 # set of param :1*rml21x1, 1*rml22x2
2 2 3.815567e+04 # mg5 will not use this value use instead 1*mdl_rml21x1
3 3 3.782868e+04 # rml23x3
#######
## INFORMATION FOR MSOFT
#######
BLOCK MSOFT #
1 1.013965e+02 # rmx1
2 1.915042e+02 # rmx2
3 5.882630e+02 # rmx3
21 3.233749e+04 # mhd2
22 -1.288001e+05 # mhu2
#######
## INFORMATION FOR MSQ2
#######
BLOCK MSQ2 #
1 1 2.998367e+05 # set of param :1*rmq21x1, 1*rmq22x2
2 2 2.998367e+05 # mg5 will not use this value use instead 1*mdl_rmq21x1
3 3 2.487654e+05 # rmq23x3
#######
## INFORMATION FOR MSU2
#######
BLOCK MSU2 #
1 1 2.803821e+05 # set of param :1*rmu21x1, 1*rmu22x2
2 2 2.803821e+05 # mg5 will not use this value use instead 1*mdl_rmu21x1
3 3 1.791371e+05 # rmu23x3
#######
## INFORMATION FOR NMIX
#######
BLOCK NMIX #
1 1 9.863644e-01 # rnn1x1
1 2 -5.311036e-02 # rnn1x2
1 3 1.464340e-01 # rnn1x3
1 4 -5.311861e-02 # rnn1x4
2 1 9.935054e-02 # rnn2x1
2 2 9.449493e-01 # rnn2x2
2 3 -2.698467e-01 # rnn2x3
2 4 1.561507e-01 # rnn2x4
3 1 -6.033880e-02 # rnn3x1
3 2 8.770049e-02 # rnn3x2
3 3 6.958775e-01 # rnn3x3
3 4 7.102270e-01 # rnn3x4
4 1 -1.165071e-01 # rnn4x1
4 2 3.107390e-01 # rnn4x2
4 3 6.492260e-01 # rnn4x3
4 4 -6.843778e-01 # rnn4x4
#######
## INFORMATION FOR SELMIX
#######
BLOCK SELMIX #
1 1 1.000000e+00 # rrl1x1
2 2 1.000000e+00 # rrl2x2
3 3 2.824872e-01 # rrl3x3
3 6 9.592711e-01 # rrl3x6
4 4 1.000000e+00 # rrl4x4
5 5 1.000000e+00 # rrl5x5
6 3 9.592711e-01 # rrl6x3
6 6 -2.824872e-01 # rrl6x6
#######
## INFORMATION FOR SMINPUTS
#######
BLOCK SMINPUTS #
1 1.279340e+02 # aewm1
3 1.180000e-01 # as
#######
## INFORMATION FOR SNUMIX
#######
BLOCK SNUMIX #
1 1 1.000000e+00 # rrn1x1
2 2 1.000000e+00 # rrn2x2
3 3 1.000000e+00 # rrn3x3
#######
## INFORMATION FOR TD
#######
BLOCK TD #
3 3 -1.106937e+02 # rtd3x3
#######
## INFORMATION FOR TE
#######
BLOCK TE #
3 3 -2.540197e+01 # rte3x3
#######
## INFORMATION FOR TU
#######
BLOCK TU #
3 3 -4.447525e+02 # rtu3x3
#######
## INFORMATION FOR UMIX
#######
BLOCK UMIX #
1 1 9.168349e-01 # ruu1x1
1 2 -3.992666e-01 # ruu1x2
2 1 3.992666e-01 # ruu2x1
2 2 9.168349e-01 # ruu2x2
#######
## INFORMATION FOR UPMNS
#######
BLOCK UPMNS #
1 1 1.000000e+00 # rmns1x1
2 2 1.000000e+00 # rmns2x2
3 3 1.000000e+00 # rmns3x3
#######
## INFORMATION FOR USQMIX
#######
BLOCK USQMIX #
1 1 1.000000e+00 # rru1x1
2 2 1.000000e+00 # rru2x2
3 3 5.536450e-01 # rru3x3
3 6 8.327528e-01 # rru3x6
4 4 1.000000e+00 # rru4x4
5 5 1.000000e+00 # rru5x5
6 3 8.327528e-01 # rru6x3
6 6 -5.536450e-01 # rru6x6
#######
## INFORMATION FOR VCKM
#######
BLOCK VCKM #
1 1 1.000000e+00 # rckm1x1
2 2 1.000000e+00 # rckm2x2
3 3 1.000000e+00 # rckm3x3
#######
## INFORMATION FOR VMIX
#######
BLOCK VMIX #
1 1 9.725578e-01 # rvv1x1
1 2 -2.326612e-01 # rvv1x2
2 1 2.326612e-01 # rvv2x1
2 2 9.725578e-01 # rvv2x2
#######
## INFORMATION FOR YD
#######
BLOCK YD #
3 3 1.388402e-01 # ryd3x3
#######
## INFORMATION FOR YE
#######
BLOCK YE #
3 3 1.008908e-01 # rye3x3
#######
## INFORMATION FOR YU
#######
BLOCK YU #
3 3 8.928445e-01 # ryu3x3
#######
## INFORMATION FOR DECAY
#######
DECAY 6 1.561950e+00 # wt
DECAY 23 2.411433e+00 # wz
DECAY 24 2.002822e+00 # ww
DECAY 25 1.986108e-03 # wh01
DECAY 35 5.748014e-01 # wh02
DECAY 36 6.321785e-01 # wa0
DECAY 37 5.469628e-01 # wh
DECAY 1000001 5.312788e+00 # wsd1
DECAY 1000002 5.477195e+00 # wsu1
DECAY 1000003 5.312788e+00 # wsd2
DECAY 1000004 5.477195e+00 # wsu2
DECAY 1000005 3.736276e+00 # wsd3
DECAY 1000006 2.021596e+00 # wsu3
DECAY 1000011 2.136822e-01 # wsl1
DECAY 1000012 1.498816e-01 # wsn1
DECAY 1000013 2.136822e-01 # wsl2
DECAY 1000014 1.498816e-01 # wsn2
DECAY 1000015 1.483273e-01 # wsl3
DECAY 1000016 1.475190e-01 # wsn3
DECAY 1000021 5.506754e+00 # wgo
DECAY 1000023 2.077700e-02 # wneu2
DECAY 1000024 1.704145e-02 # wch1
DECAY 1000025 -1.915985e+00 # wneu3
DECAY 1000035 2.585851e+00 # wneu4
DECAY 1000037 2.486895e+00 # wch2
DECAY 2000001 2.858123e-01 # wsd4
DECAY 2000002 1.152973e+00 # wsu4
DECAY 2000003 2.858123e-01 # wsd5
DECAY 2000004 1.152973e+00 # wsu5
DECAY 2000005 8.015663e-01 # wsd6
DECAY 2000006 7.373133e+00 # wsu6
DECAY 2000011 2.161216e-01 # wsl4
DECAY 2000013 2.161216e-01 # wsl5
DECAY 2000015 2.699061e-01 # wsl6
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
DECAY 1000022 0.000000e+00 # n1 : 0.0
#######
## INFORMATION FOR QNUMBERS 1000022
#######
BLOCK QNUMBERS 1000022 # n1
1 0 # 3 times electric charge
2 2 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 1000023
#######
BLOCK QNUMBERS 1000023 # n2
1 0 # 3 times electric charge
2 2 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 1000025
#######
BLOCK QNUMBERS 1000025 # n3
1 0 # 3 times electric charge
2 2 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 1000035
#######
BLOCK QNUMBERS 1000035 # n4
1 0 # 3 times electric charge
2 2 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 1000024
#######
BLOCK QNUMBERS 1000024 # x1+
1 3 # 3 times electric charge
2 2 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000037
#######
BLOCK QNUMBERS 1000037 # x2+
1 3 # 3 times electric charge
2 2 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000021
#######
BLOCK QNUMBERS 1000021 # go
1 0 # 3 times electric charge
2 2 # number of spin states (2s+1)
3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 35
#######
BLOCK QNUMBERS 35 # h2
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 36
#######
BLOCK QNUMBERS 36 # h3
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 37
#######
BLOCK QNUMBERS 37 # h+
1 3 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000012
#######
BLOCK QNUMBERS 1000012 # sve
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000014
#######
BLOCK QNUMBERS 1000014 # svm
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000016
#######
BLOCK QNUMBERS 1000016 # svt
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000011
#######
BLOCK QNUMBERS 1000011 # el-
1 -3 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000013
#######
BLOCK QNUMBERS 1000013 # mul-
1 -3 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000015
#######
BLOCK QNUMBERS 1000015 # ta1-
1 -3 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 2000011
#######
BLOCK QNUMBERS 2000011 # er-
1 -3 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 2000013
#######
BLOCK QNUMBERS 2000013 # mur-
1 -3 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 2000015
#######
BLOCK QNUMBERS 2000015 # ta2-
1 -3 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000002
#######
BLOCK QNUMBERS 1000002 # ul
1 2 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000004
#######
BLOCK QNUMBERS 1000004 # cl
1 2 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000006
#######
BLOCK QNUMBERS 1000006 # t1
1 2 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 2000002
#######
BLOCK QNUMBERS 2000002 # ur
1 2 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 2000004
#######
BLOCK QNUMBERS 2000004 # cr
1 2 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 2000006
#######
BLOCK QNUMBERS 2000006 # t2
1 2 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000001
#######
BLOCK QNUMBERS 1000001 # dl
1 -1 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000003
#######
BLOCK QNUMBERS 1000003 # sl
1 -1 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 1000005
#######
BLOCK QNUMBERS 1000005 # b1
1 -1 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 2000001
#######
BLOCK QNUMBERS 2000001 # dr
1 -1 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 2000003
#######
BLOCK QNUMBERS 2000003 # sr
1 -1 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 2000005
#######
BLOCK QNUMBERS 2000005 # b2
1 -1 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
</slha>
<MGPythiaCard>
<![CDATA[
!
! It is possible to run this card manually with:
! LD_LIBRARY_
!
!
! Pythia8 cmd card automatically generated by MadGraph5_aMC@NLO
! For more information on the use of the MG5aMC / Pythia8 interface, visit
! https:/
!
! ==================
! General parameters
! ==================
!
Main:numberOfEvents = -1
!
! -------
! Specify the HEPMC output of the Pythia8 shower. You can set it to:
! auto : MG5aMC will automatically place it the run_<i> directory
! /dev/null : to turn off the HEPMC output.
! <path> : to select where the HEPMC file must written. It will
! therefore not be placed in the run_<i> directory. The
! specified path, if not absolute, will be relative to
! the Event/run_<i> directory of the process output.
! fifo : to have MG5aMC setup the piping of the PY8 output to
! analysis tools such as MadAnalysis5.
! fifo@<fifo_path> :
! 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_
!
! -------
! 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://
! -------
! If equal to -1.0, MadGraph5_aMC@NLO will set it automatically based
! on the parameter 'xqcut' of the run_card.dat
JetMatching:qCut = 7.5000000000e+01
! 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:
! 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://
! -------
! 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,
! 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 = <set_by_user>
! 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:
!
! =======
! 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 technical parameters set by MG5_aMC.
!
! 1.0 corresponds to HEPMC weight given in [mb]. We choose here the [pb] normalization.
HEPMCoutput:
SysCalc:
! Value of the merging scale below which one does not even write the HepMC event.
SysCalc:
! This parameter is automatically set to True by MG5aMC when doing MLM merging with PY8.
Beams:setProduc
! Tell Pythia8 that an LHEF input is used.
Beams:frameType=4
! Specify one must read inputs from the MadGraph banner.
JetMatching:
JetMatching:
JetMatching:
JetMatching:
! Specifiy if we are merging sample of different multiplicity.
JetMatching:
! Do veto externally (e.g. in SysCalc).
JetMatching:
JetMatching:
! Be more forgiving with momentum mismatches.
Check:epTolErr=
!
! =======
! Subrun definitions
! =======
!
LHEFInputs:
Main:subrun=0
!
! Definition of subrun 0
!
Beams:LHEF=
]]>
</MGPythiaCard>
<MGDelphesCard>
<![CDATA[
#######
# Order of execution of various modules
#######
set ExecutionPath {
ParticlePropa
ChargedHadron
ElectronTrack
MuonTrackingE
ChargedHadron
ElectronMomen
MuonMomentumS
TrackMerger
ECal
HCal
Calorimeter
EFlowMerger
EFlowFilter
PhotonEfficiency
PhotonIsolation
ElectronFilter
ElectronEffic
ElectronIsolation
ChargedHadron
MuonEfficiency
MuonIsolation
MissingET
NeutrinoFilter
GenJetFinder
GenMissingET
FastJetFinder
FatJetFinder
JetEnergyScale
JetFlavorAsso
BTagging
TauTagging
UniqueObjectF
ScalarHT
TreeWriter
}
#######
# Propagate particles in cylinder
#######
module ParticlePropagator ParticlePropagator {
set InputArray Delphes/
set OutputArray stableParticles
set ChargedHadronOu
set ElectronOutputArray electrons
set MuonOutputArray muons
# radius of the magnetic field coverage, in m
set Radius 1.29
# half-length of the magnetic field coverage, in m
set HalfLength 3.00
# magnetic field
set Bz 3.8
}
#######
# Charged hadron tracking efficiency
#######
module Efficiency ChargedHadronTr
set InputArray ParticlePropaga
set OutputArray chargedHadrons
# add EfficiencyFormula {efficiency formula as a function of eta and pt}
# tracking efficiency formula for charged hadrons
set EfficiencyFormula { (pt <= 0.1) * (0.00) +
}
#######
# Electron tracking efficiency
#######
module Efficiency ElectronTrackin
set InputArray ParticlePropaga
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) +
}
#######
# Muon tracking efficiency
#######
module Efficiency MuonTrackingEff
set InputArray ParticlePropaga
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) +
}
#######
# Momentum resolution for charged tracks
#######
module MomentumSmearing ChargedHadronMo
set InputArray ChargedHadronTr
set OutputArray chargedHadrons
# set ResolutionFormula {resolution formula as a function of eta and pt}
# resolution formula for charged hadrons
# based on arXiv:1405.6569
set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
}
#######
# Momentum resolution for electrons
#######
module MomentumSmearing ElectronMomentu
set InputArray ElectronTrackin
set OutputArray electrons
# set ResolutionFormula {resolution formula as a function of eta and energy}
# resolution formula for electrons
# based on arXiv:1405.6569
set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
}
#######
# Momentum resolution for muons
#######
module MomentumSmearing MuonMomentumSme
set InputArray MuonTrackingEff
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) +
}
##############
# Track merger
##############
module Merger TrackMerger {
# add InputArray InputArray
add InputArray ChargedHadronMo
add InputArray ElectronMomentu
add InputArray MuonMomentumSme
set OutputArray tracks
}
#############
# ECAL
#############
module SimpleCalorimeter ECal {
set ParticleInputArray ParticlePropaga
set TrackInputArray TrackMerger/tracks
set TowerOutputArray ecalTowers
set EFlowTrackOutpu
set EFlowTowerOutpu
set IsEcal true
set EnergyMin 0.5
set EnergySignifica
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 (HGCAL- ECAL)
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}
# for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
# set ECalResolutionF
# Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^
}
#############
# HCAL
#############
module SimpleCalorimeter HCal {
set ParticleInputArray ParticlePropaga
set TrackInputArray ECal/eflowTracks
set TowerOutputArray hcalTowers
set EFlowTrackOutpu
set EFlowTowerOutpu
set IsEcal false
set EnergyMin 1.0
set EnergySignifica
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
# 5 degrees towers
set PhiBins {}
for {set i -36} {$i <= 36} {incr i} {
add PhiBins [expr {$i * $pi/36.0}]
}
foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
add EtaPhiBins $eta $PhiBins
}
# 10 degrees towers
set PhiBins {}
for {set i -18} {$i <= 18} {incr i} {
add PhiBins [expr {$i * $pi/18.0}]
}
foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
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 {-5 -4.7 -4.525 4.7 5} {
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}
# set HCalResolutionF
set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^
}
#################
# 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/eflowNeutr
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) +
}
##################
# Photon isolation
##################
module Isolation PhotonIsolation {
set CandidateInputArray PhotonEfficienc
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/
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) +
}
#######
# Electron isolation
#######
module Isolation ElectronIsolation {
set CandidateInputArray ElectronEfficie
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 MuonMomentumSme
set OutputArray muons
# set EfficiencyFormula {efficiency as a function of eta and pt}
# efficiency formula for muons
set EfficiencyFormula { (pt <= 10.0) * (0.00) +
}
################
# Muon isolation
################
module Isolation MuonIsolation {
set CandidateInputArray MuonEfficiency/
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 EFlowMerger/eflow
set MomentumOutputArray momentum
}
##################
# Scalar HT merger
##################
module Merger ScalarHT {
# add InputArray InputArray
add InputArray UniqueObjectFin
add InputArray UniqueObjectFin
add InputArray UniqueObjectFin
add InputArray UniqueObjectFin
set EnergyOutputArray energy
}
#######
# Neutrino Filter
#######
module PdgCodeFilter NeutrinoFilter {
set InputArray Delphes/
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/
set OutputArray jets
# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
set JetAlgorithm 6
set ParameterR 0.5
set JetPTMin 20.0
}
#######
# Gen Missing ET merger
#######
module Merger GenMissingET {
# add InputArray InputArray
add InputArray NeutrinoFilter/
set MomentumOutputArray momentum
}
############
# Jet finder
############
module FastJetFinder FastJetFinder {
# set InputArray Calorimeter/towers
set InputArray EFlowMerger/eflow
set OutputArray jets
# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
set JetAlgorithm 6
set ParameterR 0.5
set JetPTMin 20.0
}
##################
# Fat Jet finder
##################
module FastJetFinder FatJetFinder {
set InputArray EFlowMerger/eflow
set OutputArray jets
# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
set JetAlgorithm 6
set ParameterR 0.8
set ComputeNsubjett
set Beta 1.0
set AxisMode 4
set ComputeTrimming 1
set RTrim 0.2
set PtFracTrim 0.05
set ComputePruning 1
set ZcutPrun 0.1
set RcutPrun 0.5
set RPrun 0.8
set ComputeSoftDrop 1
set BetaSoftDrop 0.0
set SymmetryCutSoftDrop 0.1
set R0SoftDrop 0.8
set JetPTMin 200.0
}
##################
# Jet Energy Scale
##################
module EnergyScale JetEnergyScale {
set InputArray FastJetFinder/jets
set OutputArray jets
# scale formula for jets
set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
}
#######
# Jet Flavor Association
#######
module JetFlavorAssoci
set PartonInputArray Delphes/partons
set ParticleInputArray Delphes/
set ParticleLHEFInp
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 arXiv:1211.4462
# default efficiency formula (misidentification rate)
add EfficiencyFormula {0} {0.01+0.000038*pt}
# efficiency formula for c-jets (misidentification rate)
add EfficiencyFormula {4} {0.25*tanh(
# efficiency formula for b-jets
add EfficiencyFormula {5} {0.85*tanh(
}
#############
# tau-tagging
#############
module TauTagging TauTagging {
set ParticleInputArray Delphes/
set PartonInputArray Delphes/partons
set JetInputArray JetEnergyScale/jets
set DeltaR 0.5
set TauPTMin 1.0
set TauEtaMax 2.5
# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
# default efficiency formula (misidentification rate)
add EfficiencyFormula {0} {0.01}
# efficiency formula for tau-jets
add EfficiencyFormula {15} {0.6}
}
#######
# Find uniquely identified photons/
#######
module UniqueObjectFinder UniqueObjectFinder {
# earlier arrays take precedence over later ones
# add InputArray InputArray OutputArray
add InputArray PhotonIsolation
add InputArray ElectronIsolati
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/
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/eflowNeutr
add Branch GenJetFinder/jets GenJet Jet
add Branch GenMissingET/
add Branch UniqueObjectFin
add Branch UniqueObjectFin
add Branch UniqueObjectFin
add Branch UniqueObjectFin
add Branch FatJetFinder/jets FatJet Jet
add Branch MissingET/momentum MissingET MissingET
add Branch ScalarHT/energy ScalarHT ScalarHT
}
]]>
</MGDelphesCard>
<MGGenerationInfo>
# Number of Events : 10000
# Integrated weight (pb) : 3.09652803909
# Matched Integrated weight (pb) : -1.0
</MGGenerationInfo>
</header>
<init>
2212 2212 6.500000e+03 6.500000e+03 0 0 247000 247000 -4 3
9.016164e-01 3.272751e-03 3.096528e+00 1
1.502443e+00 3.190257e-03 3.096528e+00 3
6.924697e-01 3.532158e-03 3.096528e+00 2
<generator name='MadGraph5
</init>
</LesHouchesEvents>
Revision history for this message
|
#2 |
=======
-> Pythia8 log file for run 7 <-
=======
*-----
| |
| *------
| | | |
| | | |
| | PPP Y Y TTTTT H H III A Welcome to the Lund Monte Carlo! | |
| | P P Y Y T H H I A A This is PYTHIA version 8.235 | |
| | PPP Y T HHHHH I AAAAA Last date of change: 27 Mar 2018 | |
| | P Y T H H I A A | |
| | P Y T H H III A A Now is 18 Jun 2018 at 17:52:31 | |
| | | |
| | Christian Bierlich; Department of Astronomy and Theoretical Physics, | |
| | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | |
| | e-mail: <email address hidden> | |
| | Nishita Desai; Laboratoire Charles Coulomb (L2C), | |
| | CNRS-Universite de Montpellier, 34090 Montpellier, France; | |
| | e-mail: <email address hidden> | |
| | Nadine Fischer; School of Physics, | |
| | Monash University, PO Box 27, 3800 Melbourne, Australia; | |
| | e-mail: <email address hidden> | |
| | Ilkka Helenius; Institute for Theoretical Physics, | |
| | Tuebingen University, Auf der Morgenstelle 14, 72076 Tuebingen, Germany; | |
| | e-mail: <email address hidden> | |
| | Philip Ilten; School of Physics and Astronomy, | |
| | University of Birmingham, Birmingham, B152 2TT, UK; | |
| | e-mail: <email address hidden> | |
| | Leif Lonnblad; Department of Astronomy and Theoretical Physics, | |
| | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | |
| | e-mail: <email address hidden> | |
| | Stephen Mrenna; Computing Division, Simulations Group, | |
| | Fermi National Accelerator Laboratory, MS 234, Batavia, IL 60510, USA; | |
| | e-mail: <email address hidden> | |
| | Stefan Prestel; Theoretical Physics Department, | |
| | Fermi National Accelerator Laboratory, MS 106, Batavia, IL 60510, USA; | |
| | e-mail: <email address hidden> | |
| | Christine O. Rasmussen; Department of Astronomy and Theoretical Physics, | |
| | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | |
| | e-mail: <email address hidden> | |
| | Torbjorn Sjostrand; Department of Astronomy and Theoretical Physics, | |
| | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | |
| | e-mail: <email address hidden> | |
| | Peter Skands; School of Physics, | |
| | Monash University, PO Box 27, 3800 Melbourne, Australia; | |
| | e-mail: <email address hidden> | |
| | | |
| | The main program reference is 'An Introduction to PYTHIA 8.2', | |
| | T. Sjostrand et al, Comput. Phys. Commun. 191 (2015) 159 | |
| | [arXiv:1410.3012 [hep-ph]] | |
| | | |
| | The main physics reference is the 'PYTHIA 6.4 Physics and Manual', | |
| | T. Sjostrand, S. Mrenna and P. Skands, JHEP05 (2006) 026 [hep-ph/0603175] | |
| | | |
| | An archive of program versions and documentation is found on the web: | |
| | http://
| | | |
| | This program is released under the GNU General Public Licence version 2. | |
| | Please respect the MCnet Guidelines for Event Generator Authors and Users. | |
| | | |
| | Disclaimer: this program comes without any guarantees. | |
| | Beware of errors and use common sense when interpreting results. | |
| | | |
| | Copyright (C) 2018 Torbjorn Sjostrand | |
| | | |
| | | |
| *------
| |
*-----
*-------- Madgraph parameters --------*
*-----
*----- Madgraph matching parameters -----*
| qCut | 75.000 |
| nQmatch | 4 |
| clFact | 1.000 |
| Jet algorithm | kT |
| eTjetMin | 20.000 |
| etaJetMax | 1000.000 |
| jetAllow | 1 |
| Mode | exclusive |
*-----
*-----
| Last Change 12 Apr 2017 - P. Skands
| Parsing: events.lhe.gz
| (SLHA::readFile) line 430 - storing non-SLHA(2) block: fralpha
| (SLHA::readFile) line 705 - storing QNUMBERS for id = 1000022 n1
| (SLHA::readFile) line 713 - storing QNUMBERS for id = 1000023 n2
| (SLHA::readFile) line 721 - storing QNUMBERS for id = 1000025 n3
| (SLHA::readFile) line 729 - storing QNUMBERS for id = 1000035 n4
| (SLHA::readFile) line 737 - storing QNUMBERS for id = 1000024 x1+
| (SLHA::readFile) line 745 - storing QNUMBERS for id = 1000037 x2+
| (SLHA::readFile) line 753 - storing QNUMBERS for id = 1000021 go
| (SLHA::readFile) line 761 - storing QNUMBERS for id = 35 h2
| (SLHA::readFile) line 769 - storing QNUMBERS for id = 36 h3
| (SLHA::readFile) line 777 - storing QNUMBERS for id = 37 h+
| (SLHA::readFile) line 785 - storing QNUMBERS for id = 1000012 sve
| (SLHA::readFile) line 793 - storing QNUMBERS for id = 1000014 svm
| (SLHA::readFile) line 801 - storing QNUMBERS for id = 1000016 svt
| (SLHA::readFile) line 809 - storing QNUMBERS for id = 1000011 el-
| (SLHA::readFile) line 817 - storing QNUMBERS for id = 1000013 mul-
| (SLHA::readFile) line 825 - storing QNUMBERS for id = 1000015 ta1-
| (SLHA::readFile) line 833 - storing QNUMBERS for id = 2000011 er-
| (SLHA::readFile) line 841 - storing QNUMBERS for id = 2000013 mur-
| (SLHA::readFile) line 849 - storing QNUMBERS for id = 2000015 ta2-
| (SLHA::readFile) line 857 - storing QNUMBERS for id = 1000002 ul
| (SLHA::readFile) line 865 - storing QNUMBERS for id = 1000004 cl
| (SLHA::readFile) line 873 - storing QNUMBERS for id = 1000006 t1
| (SLHA::readFile) line 881 - storing QNUMBERS for id = 2000002 ur
| (SLHA::readFile) line 889 - storing QNUMBERS for id = 2000004 cr
| (SLHA::readFile) line 897 - storing QNUMBERS for id = 2000006 t2
| (SLHA::readFile) line 905 - storing QNUMBERS for id = 1000001 dl
| (SLHA::readFile) line 913 - storing QNUMBERS for id = 1000003 sl
| (SLHA::readFile) line 921 - storing QNUMBERS for id = 1000005 b1
| (SLHA::readFile) line 929 - storing QNUMBERS for id = 2000001 dr
| (SLHA::readFile) line 937 - storing QNUMBERS for id = 2000003 sr
| (SLHA::readFile) line 945 - storing QNUMBERS for id = 2000005 b2
*-----
PYTHIA Info from SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Info from SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
*------- PYTHIA Process Initialization -------
| |
| We collide p+ with p+ at a CM energy of 1.300e+04 GeV |
| |
|-----
| | |
| Subprocess Code | Estimated |
| | max (mb) |
| | |
|-----
| | |
| Les Houches User Process(es) 9999 | 3.000e-09 |
| |
*------- End PYTHIA Process Initialization -------
*------- PYTHIA Multiparton Interactions Initialization ---------*
| |
| sigmaNonDiffractive = 56.42 mb |
| |
| pT0 = 2.60 gives sigmaInteraction = 321.22 mb: accepted |
| |
*------- End PYTHIA Multiparton Interactions Initialization -----*
*------- PYTHIA Flag + Mode + Parm + Word + FVec + MVec + PVec + WVec Settings (changes only) ------------------*
| |
| Name | Now | Default Min Max |
| | | |
| Beams:frameType | 4 | 1 1 5 |
| Beams:LHEF | events.lhe.gz | void |
| Beams:setProduc
| Check:epTolErr | 0.0100000 | 1.0000e-04 |
| HEPMCoutput:file | events.hepmc | void |
| HEPMCoutput:scaling | 1.2500e+08 | 1.00000 0.0 |
| JetMatching:
| JetMatching:doVeto | off | on |
| JetMatching:
| JetMatching:merge | on | off |
| JetMatching:nJetMax | 2 | -1 -1 |
| JetMatching:nQmatch | 4 | 5 3 6 |
| JetMatching:qCut | 75.00000 | 10.00000 0.0 |
| JetMatching:setMad | off | on |
| LHEFInputs:nSubruns | 1 | 0 0 |
| Main:numberOfEvents | 1250 | 1000 0 |
| Main:subrun | 0 | -999 0 |
| Syscalc:qCutList | 112.50000 | 0.0 |
| | 150.00000 | 0.0 |
| | | 0.0 |
| | | 0.0 |
| | | 0.0 |
| | | 0.0 |
| | | 0.0 |
| | | 0.0 |
| | | 0.0 |
| | | 0.0 |
| SysCalc:qWeed | 50.00000 | 0.0 0.0 |
| |
*------- End PYTHIA Flag + Mode + Parm + Word + FVec + MVec + PVec + WVec Settings -------
-------- PYTHIA Particle Data Table (changed only) -------
id name antiName spn chg col m0 mWidth mMin mMax tau0 res dec ext vis wid
no onMode bRatio meMode products
25 h0 1 0 0 110.89910 0.00287 110.88917 110.90903 9.91890e-11 1 1 0 0 0
0 1 0.0000011 0 1 -1
1 1 0.0000003 0 2 -2
2 1 0.0002781 0 3 -3
3 1 0.0370756 0 4 -4
4 1 0.7397740 0 5 -5
5 1 0.0000000 0 6 -6
6 1 0.0000000 0 11 -11
7 1 0.0002756 0 13 -13
8 1 0.0792479 0 15 -15
9 1 0.0850015 0 21 21
10 1 0.0020250 0 22 22
11 1 0.0004654 0 22 23
12 1 0.0044078 0 23 23
13 1 0.0514476 0 24 -24
14 1 0.0000000 103 1000022 1000022
15 1 0.0000000 103 1000023 1000022
16 1 0.0000000 103 1000023 1000023
17 1 0.0000000 103 1000025 1000022
18 1 0.0000000 103 1000025 1000023
19 1 0.0000000 103 1000025 1000025
20 1 0.0000000 103 1000035 1000022
21 1 0.0000000 103 1000035 1000023
22 1 0.0000000 103 1000035 1000025
23 1 0.0000000 103 1000035 1000035
24 1 0.0000000 103 1000024 -1000024
25 1 0.0000000 103 1000024 -1000037
26 1 0.0000000 103 1000037 -1000024
27 1 0.0000000 103 1000037 -1000037
28 1 0.0000000 103 1000001 -1000001
29 1 0.0000000 103 2000001 -2000001
30 1 0.0000000 103 1000001 -2000001
31 1 0.0000000 103 -1000001 2000001
32 1 0.0000000 103 1000002 -1000002
33 1 0.0000000 103 2000002 -2000002
34 1 0.0000000 103 1000002 -2000002
35 1 0.0000000 103 -1000002 2000002
36 1 0.0000000 103 1000003 -1000003
37 1 0.0000000 103 2000003 -2000003
38 1 0.0000000 103 1000003 -2000003
39 1 0.0000000 103 -1000003 2000003
40 1 0.0000000 103 1000004 -1000004
41 1 0.0000000 103 2000004 -2000004
42 1 0.0000000 103 1000004 -2000004
43 1 0.0000000 103 -1000004 2000004
44 1 0.0000000 103 1000005 -1000005
45 1 0.0000000 103 2000005 -2000005
46 1 0.0000000 103 1000005 -2000005
47 1 0.0000000 103 -1000005 2000005
48 1 0.0000000 103 1000006 -1000006
49 1 0.0000000 103 2000006 -2000006
50 1 0.0000000 103 1000006 -2000006
51 1 0.0000000 103 -1000006 2000006
52 1 0.0000000 103 1000011 -1000011
53 1 0.0000000 103 2000011 -2000011
54 1 0.0000000 103 1000011 -2000011
55 1 0.0000000 103 -1000011 2000011
56 1 0.0000000 103 1000012 -1000012
57 1 0.0000000 103 2000012 -2000012
58 1 0.0000000 103 1000012 -2000012
59 1 0.0000000 103 -1000012 2000012
60 1 0.0000000 103 1000013 -1000013
61 1 0.0000000 103 2000013 -2000013
62 1 0.0000000 103 1000013 -2000013
63 1 0.0000000 103 -1000013 2000013
64 1 0.0000000 103 1000014 -1000014
65 1 0.0000000 103 2000014 -2000014
66 1 0.0000000 103 1000014 -2000014
67 1 0.0000000 103 -1000014 2000014
68 1 0.0000000 103 1000015 -1000015
69 1 0.0000000 103 2000015 -2000015
70 1 0.0000000 103 1000015 -2000015
71 1 0.0000000 103 -1000015 2000015
72 1 0.0000000 103 1000016 -1000016
73 1 0.0000000 103 2000016 -2000016
74 1 0.0000000 103 1000016 -2000016
75 1 0.0000000 103 -1000016 2000016
35 H0 1 0 0 399.96010 0.57480 397.08609 402.83411 3.42727e-13 1 1 0 0 0
0 1 0.0000000 0 1 -1
1 1 0.0000000 0 2 -2
2 1 0.0000000 0 3 -3
3 1 0.0000490 0 4 -4
4 1 0.0007740 0 5 -5
5 1 0.0000000 0 6 -6
6 1 0.0000000 0 11 -11
7 1 0.0000000 0 13 -13
8 1 0.0000740 0 15 -15
9 1 0.0004170 0 21 21
10 1 0.0000150 0 22 22
11 1 0.0000610 0 22 23
12 1 0.3067100 0 23 23
13 1 0.6890110 0 24 -24
14 1 0.0000000 0 23 25
15 1 0.0028890 0 25 25
16 1 0.0000000 0 24 -37
17 1 0.0000000 0 37 -24
18 1 0.0000000 0 23 36
19 1 0.0000000 0 25 36
20 1 0.0000000 0 36 36
21 1 0.0000000 103 1000022 1000022
22 1 0.0000000 103 1000023 1000022
23 1 0.0000000 103 1000023 1000023
24 1 0.0000000 103 1000025 1000022
25 1 0.0000000 103 1000025 1000023
26 1 0.0000000 103 1000025 1000025
27 1 0.0000000 103 1000035 1000022
28 1 0.0000000 103 1000035 1000023
29 1 0.0000000 103 1000035 1000025
30 1 0.0000000 103 1000035 1000035
31 1 0.0000000 103 1000024 -1000024
32 1 0.0000000 103 1000024 -1000037
33 1 0.0000000 103 1000037 -1000024
34 1 0.0000000 103 1000037 -1000037
35 1 0.0000000 103 1000001 -1000001
36 1 0.0000000 103 2000001 -2000001
37 1 0.0000000 103 1000001 -2000001
38 1 0.0000000 103 -1000001 2000001
39 1 0.0000000 103 1000002 -1000002
40 1 0.0000000 103 2000002 -2000002
41 1 0.0000000 103 1000002 -2000002
42 1 0.0000000 103 -1000002 2000002
43 1 0.0000000 103 1000003 -1000003
44 1 0.0000000 103 2000003 -2000003
45 1 0.0000000 103 1000003 -2000003
46 1 0.0000000 103 -1000003 2000003
47 1 0.0000000 103 1000004 -1000004
48 1 0.0000000 103 2000004 -2000004
49 1 0.0000000 103 1000004 -2000004
50 1 0.0000000 103 -1000004 2000004
51 1 0.0000000 103 1000005 -1000005
52 1 0.0000000 103 2000005 -2000005
53 1 0.0000000 103 1000005 -2000005
54 1 0.0000000 103 -1000005 2000005
55 1 0.0000000 103 1000006 -1000006
56 1 0.0000000 103 2000006 -2000006
57 1 0.0000000 103 1000006 -2000006
58 1 0.0000000 103 -1000006 2000006
59 1 0.0000000 103 1000011 -1000011
60 1 0.0000000 103 2000011 -2000011
61 1 0.0000000 103 1000011 -2000011
62 1 0.0000000 103 -1000011 2000011
63 1 0.0000000 103 1000012 -1000012
64 1 0.0000000 103 2000012 -2000012
65 1 0.0000000 103 1000012 -2000012
66 1 0.0000000 103 -1000012 2000012
67 1 0.0000000 103 1000013 -1000013
68 1 0.0000000 103 2000013 -2000013
69 1 0.0000000 103 1000013 -2000013
70 1 0.0000000 103 -1000013 2000013
71 1 0.0000000 103 1000014 -1000014
72 1 0.0000000 103 2000014 -2000014
73 1 0.0000000 103 1000014 -2000014
74 1 0.0000000 103 -1000014 2000014
75 1 0.0000000 103 1000015 -1000015
76 1 0.0000000 103 2000015 -2000015
77 1 0.0000000 103 1000015 -2000015
78 1 0.0000000 103 -1000015 2000015
79 1 0.0000000 103 1000016 -1000016
80 1 0.0000000 103 2000016 -2000016
81 1 0.0000000 103 1000016 -2000016
82 1 0.0000000 103 -1000016 2000016
36 A0 1 0 0 399.58390 0.63218 396.42301 402.74479 3.11621e-13 1 1 0 0 0
0 1 0.0000000 0 1 -1
1 1 0.0000000 0 2 -2
2 1 0.0000010 0 3 -3
3 1 0.0001210 0 4 -4
4 1 0.0019240 0 5 -5
5 1 0.0000000 0 6 -6
6 1 0.0000000 0 11 -11
7 1 0.0000010 0 13 -13
8 1 0.0001840 0 15 -15
9 1 0.0031060 0 21 21
10 1 0.0000150 0 22 22
11 1 0.0000030 0 22 23
12 1 0.0000000 0 23 23
13 1 0.0000000 0 24 -24
14 1 0.9946450 0 23 25
15 1 0.0000000 0 25 25
16 1 0.0000000 0 24 -37
17 1 0.0000000 0 37 -24
18 1 0.0000000 103 1000022 1000022
19 1 0.0000000 103 1000023 1000022
20 1 0.0000000 103 1000023 1000023
21 1 0.0000000 103 1000025 1000022
22 1 0.0000000 103 1000025 1000023
23 1 0.0000000 103 1000025 1000025
24 1 0.0000000 103 1000035 1000022
25 1 0.0000000 103 1000035 1000023
26 1 0.0000000 103 1000035 1000025
27 1 0.0000000 103 1000035 1000035
28 1 0.0000000 103 1000024 -1000024
29 1 0.0000000 103 1000024 -1000037
30 1 0.0000000 103 1000037 -1000024
31 1 0.0000000 103 1000037 -1000037
32 1 0.0000000 103 1000001 -1000001
33 1 0.0000000 103 2000001 -2000001
34 1 0.0000000 103 1000001 -2000001
35 1 0.0000000 103 -1000001 2000001
36 1 0.0000000 103 1000002 -1000002
37 1 0.0000000 103 2000002 -2000002
38 1 0.0000000 103 1000002 -2000002
39 1 0.0000000 103 -1000002 2000002
40 1 0.0000000 103 1000003 -1000003
41 1 0.0000000 103 2000003 -2000003
42 1 0.0000000 103 1000003 -2000003
43 1 0.0000000 103 -1000003 2000003
44 1 0.0000000 103 1000004 -1000004
45 1 0.0000000 103 2000004 -2000004
46 1 0.0000000 103 1000004 -2000004
47 1 0.0000000 103 -1000004 2000004
48 1 0.0000000 103 1000005 -1000005
49 1 0.0000000 103 2000005 -2000005
50 1 0.0000000 103 1000005 -2000005
51 1 0.0000000 103 -1000005 2000005
52 1 0.0000000 103 1000006 -1000006
53 1 0.0000000 103 2000006 -2000006
54 1 0.0000000 103 1000006 -2000006
55 1 0.0000000 103 -1000006 2000006
56 1 0.0000000 103 1000011 -1000011
57 1 0.0000000 103 2000011 -2000011
58 1 0.0000000 103 1000011 -2000011
59 1 0.0000000 103 -1000011 2000011
60 1 0.0000000 103 1000012 -1000012
61 1 0.0000000 103 2000012 -2000012
62 1 0.0000000 103 1000012 -2000012
63 1 0.0000000 103 -1000012 2000012
64 1 0.0000000 103 1000013 -1000013
65 1 0.0000000 103 2000013 -2000013
66 1 0.0000000 103 1000013 -2000013
67 1 0.0000000 103 -1000013 2000013
68 1 0.0000000 103 1000014 -1000014
69 1 0.0000000 103 2000014 -2000014
70 1 0.0000000 103 1000014 -2000014
71 1 0.0000000 103 -1000014 2000014
72 1 0.0000000 103 1000015 -1000015
73 1 0.0000000 103 2000015 -2000015
74 1 0.0000000 103 1000015 -2000015
75 1 0.0000000 103 -1000015 2000015
76 1 0.0000000 103 1000016 -1000016
77 1 0.0000000 103 2000016 -2000016
78 1 0.0000000 103 1000016 -2000016
79 1 0.0000000 103 -1000016 2000016
37 H+ H- 1 3 0 407.87900 0.54696 405.14419 410.61381 3.60171e-13 1 1 0 1 0
0 1 0.0000000 0 -1 2
1 1 0.0000210 0 -3 4
2 1 0.0901350 0 -5 6
3 1 0.0000000 0 -11 12
4 1 0.0000130 0 -13 14
5 1 0.0037140 0 -15 16
6 1 0.9061170 0 24 25
7 1 0.0000000 103 1000022 1000024
8 1 0.0000000 103 1000022 1000037
9 1 0.0000000 103 1000023 1000024
10 1 0.0000000 103 1000023 1000037
11 1 0.0000000 103 1000025 1000024
12 1 0.0000000 103 1000025 1000037
13 1 0.0000000 103 1000035 1000024
14 1 0.0000000 103 1000035 1000037
15 1 0.0000000 103 1000006 -1000005
16 1 0.0000000 103 2000006 -1000005
17 1 0.0000000 103 1000006 -2000005
18 1 0.0000000 103 2000006 -2000005
19 1 0.0000000 103 -1000001 1000002
20 1 0.0000000 103 -1000003 1000004
21 1 0.0000000 103 -1000011 1000012
22 1 0.0000000 103 -1000013 1000014
23 1 0.0000000 103 -1000015 1000016
24 1 0.0000000 103 -2000015 1000016
1000001 ~d_L ~d_Lbar 1 -1 1 568.44110 0.00000 541.87716 595.00504 3.70803e-14 1 0 0 1 0
0 1 0.0000000 103 1000039 1
1 1 0.0000000 103 1000039 3
2 1 0.0000000 103 1000039 5
3 1 0.0000000 0 -1000024 2
4 1 0.0000000 0 -1000037 2
5 1 0.0000000 0 -1000024 4
6 1 0.0000000 0 -1000037 4
7 1 0.0000000 0 -1000024 6
8 1 0.0000000 0 -1000037 6
9 1 0.0000000 0 1000022 1
10 1 0.0000000 0 1000022 3
11 1 0.0000000 0 1000022 5
12 1 0.0000000 0 1000023 1
13 1 0.0000000 0 1000023 3
14 1 0.0000000 0 1000023 5
15 1 0.0000000 0 1000025 1
16 1 0.0000000 0 1000025 3
17 1 0.0000000 0 1000025 5
18 1 0.0000000 0 1000035 1
19 1 0.0000000 0 1000035 3
20 1 0.0000000 0 1000035 5
21 1 0.0000000 0 1000002 -24
22 1 0.0000000 0 1000004 -24
23 1 0.0000000 0 1000006 -24
24 1 0.0000000 0 2000002 -24
25 1 0.0000000 0 2000004 -24
26 1 0.0000000 0 2000006 -24
27 1 0.0000000 0 1000002 -37
28 1 0.0000000 0 1000004 -37
29 1 0.0000000 0 1000006 -37
30 1 0.0000000 0 2000002 -37
31 1 0.0000000 0 2000004 -37
32 1 0.0000000 0 2000006 -37
33 1 0.0000000 0 1000021 1
34 1 0.0000000 0 1000021 2
35 1 0.0000000 0 1000021 5
36 1 0.0000000 0 -12 1
37 1 0.0000000 0 -12 3
38 1 0.0000000 0 -12 5
39 1 0.0000000 0 -14 1
40 1 0.0000000 0 -14 3
41 1 0.0000000 0 -14 5
42 1 0.0000000 0 -16 1
43 1 0.0000000 0 -16 3
44 1 0.0000000 0 -16 5
45 1 0.0000000 0 12 1
46 1 0.0000000 0 11 2
47 1 0.0000000 0 12 3
48 1 0.0000000 0 11 4
49 1 0.0000000 0 12 5
50 1 0.0000000 0 11 6
51 1 0.0000000 0 14 1
52 1 0.0000000 0 13 2
53 1 0.0000000 0 14 3
54 1 0.0000000 0 13 4
55 1 0.0000000 0 14 5
56 1 0.0000000 0 13 6
57 1 0.0000000 0 16 1
58 1 0.0000000 0 15 2
59 1 0.0000000 0 16 3
60 1 0.0000000 0 15 4
61 1 0.0000000 0 16 5
62 1 0.0000000 0 15 6
63 1 0.0000000 0 -2 -3
64 1 0.0000000 0 -2 -5
65 1 0.0000000 0 -4 -3
66 1 0.0000000 0 -4 -5
67 1 0.0000000 0 -6 -3
68 1 0.0000000 0 -6 -5
1000002 ~u_L ~u_Lbar 1 2 1 561.11900 0.00000 533.73302 588.50498 3.59673e-14 1 0 0 1 0
0 1 0.0000000 103 1000039 2
1 1 0.0000000 103 1000039 4
2 1 0.0000000 103 1000039 6
3 1 0.0000000 0 1000024 1
4 1 0.0000000 0 1000024 3
5 1 0.0000000 0 1000024 5
6 1 0.0000000 0 1000037 1
7 1 0.0000000 0 1000037 3
8 1 0.0000000 0 1000037 5
9 1 0.0000000 0 1000022 2
10 1 0.0000000 0 1000022 4
11 1 0.0000000 0 1000022 6
12 1 0.0000000 0 1000023 2
13 1 0.0000000 0 1000023 4
14 1 0.0000000 0 1000023 6
15 1 0.0000000 0 1000025 2
16 1 0.0000000 0 1000025 4
17 1 0.0000000 0 1000025 6
18 1 0.0000000 0 1000035 2
19 1 0.0000000 0 1000035 4
20 1 0.0000000 0 1000035 6
21 1 0.0000000 0 1000001 -24
22 1 0.0000000 0 1000003 -24
23 1 0.0000000 0 1000005 -24
24 1 0.0000000 0 2000001 -24
25 1 0.0000000 0 2000003 -24
26 1 0.0000000 0 2000005 -24
27 1 0.0000000 0 1000001 -37
28 1 0.0000000 0 1000003 -37
29 1 0.0000000 0 1000005 -37
30 1 0.0000000 0 2000001 -37
31 1 0.0000000 0 2000003 -37
32 1 0.0000000 0 2000005 -37
33 1 0.0000000 0 1000021 2
34 1 0.0000000 0 1000021 4
35 1 0.0000000 0 1000021 6
36 1 0.0000000 0 -11 1
37 1 0.0000000 0 -11 3
38 1 0.0000000 0 -11 5
39 1 0.0000000 0 -13 1
40 1 0.0000000 0 -13 3
41 1 0.0000000 0 -13 5
42 1 0.0000000 0 -15 1
43 1 0.0000000 0 -15 3
44 1 0.0000000 0 -15 5
45 1 0.0000000 0 -1 -3
46 1 0.0000000 0 -1 -5
47 1 0.0000000 0 -3 -5
Revision history for this message
|
#3 |
-------- End PYTHIA Particle Data Table -------
-------- LHA initialization information ------------
beam kind energy pdfgrp pdfset
A 2212 6500.000 0 247000
B 2212 6500.000 0 247000
Event weighting strategy = -4
Processes, with strategy-dependent cross section info
number xsec (pb) xerr (pb) xmax (pb)
1 9.0162e-01 3.2728e-03 3.0965e+00
3 1.5024e+00 3.1903e-03 3.0965e+00
2 6.9247e-01 3.5322e-03 3.0965e+00
-------- End LHA initialization information --------
-------- LHA event information and listing -------
process = 3 weight = 3.0965e+00 scale = 7.3707e+02 (GeV)
Participating Particles
no id stat mothers colours p_x p_y p_z e m tau spin
1 1 -1 0 0 501 0 -0.000 0.000 944.379 944.379 0.000 0.000 -1.000
2 2 -1 0 0 502 0 0.000 -0.000 -631.452 631.452 0.000 0.000 -1.000
3 1000001 2 1 2 502 0 -143.474 450.543 394.595 836.048 565.418 0.000 0.000
4 1000002 2 1 2 501 0 143.474 -450.543 -81.668 739.783 563.059 0.000 0.000
5 -1000024 1 3 3 0 0 -339.970 285.919 238.697 536.022 181.696 0.000 -1.000
6 1000023 1 4 4 0 0 38.302 -130.509 -271.139 353.283 181.088 0.000 -1.000
7 2 1 4 4 501 0 105.172 -320.034 189.472 386.500 0.000 0.000 -1.000
8 2 1 3 3 502 0 196.496 164.624 155.898 300.026 0.000 0.000 -1.000
-------- End LHA event information and listing -------
-------- PYTHIA Info Listing -------
Beam A: id = 2212, pz = 6.500e+03, e = 6.500e+03, m = 9.383e-01.
Beam B: id = 2212, pz = -6.500e+03, e = 6.500e+03, m = 9.383e-01.
In 1: id = 1, x = 1.453e-01, pdf = 0.000e+00 at Q2 = 5.433e+05.
In 2: id = 2, x = 9.715e-02, pdf = 0.000e+00 at same Q2.
Process Les Houches User Process(es) with code 9999 is 2 -> 2.
Subprocess user process 3 with code 3 is 2 -> 2.
It has sHat = 2.385e+06, tHat = -5.141e+05, uHat = -1.234e+06,
pTHat = 4.728e+02, m3Hat = 5.654e+02, m4Hat = 5.631e+02,
thetaHat = 1.113e+00, phiHat = 1.879e+00.
alphaEM = 7.817e-03, alphaS = 9.657e-02 at Q2 = 5.433e+05.
Impact parameter b = 2.553e-01 gives enhancement factor = 3.601e+00.
Max pT scale for MPI = 1.300e+04, ISR = 1.300e+04, FSR = 1.300e+04.
Number of MPI = 13, ISR = 33, FSRproc = 122, FSRreson = 9.
-------- End PYTHIA Info Listing -------
-------- PYTHIA Event Listing (hard process) -------
no id name status mothers daughters colours p_x p_y p_z e m
0 90 (system) -11 0 0 0 0 0 0 0.000 0.000 0.000 13000.000 13000.000
1 2212 (p+) -12 0 0 3 0 0 0 0.000 0.000 6500.000 6500.000 0.938
2 2212 (p+) -12 0 0 4 0 0 0 0.000 0.000 -6500.000 6500.000 0.938
3 1 (d) -21 1 0 5 6 501 0 -0.000 0.000 944.379 944.379 0.000
4 2 (u) -21 2 0 5 6 502 0 0.000 -0.000 -631.452 631.452 0.000
5 1000001 (~d_L) -22 3 4 7 8 502 0 -143.474 450.543 394.595 836.048 565.418
6 1000002 (~u_L) -22 3 4 9 10 501 0 143.474 -450.543 -81.668 739.783 563.059
7 -1000024 ~chi_1- 23 5 0 0 0 0 0 -339.970 285.919 238.697 536.022 181.696
8 2 u 23 5 0 0 0 502 0 196.496 164.624 155.898 300.026 0.000
9 1000023 ~chi_20 23 6 0 0 0 0 0 38.302 -130.509 -271.139 353.283 181.088
10 2 u 23 6 0 0 0 501 0 105.172 -320.034 189.472 386.500 0.000
Revision history for this message
|
#4 |
According to this line (and to the param_card that you send)
PYTHIA Warning in SLHAinterface:
You did not specify within the param_card what the Branching ratio are for that BSM particle.
So this is not surprising. If you want the BR to be computed by MG5aMC (at LO without any loop-induced contribution), You can set their width to Auto in the param_card.
Cheers,
Olivier
> On 19 Jun 2018, at 02:46, Zhang Yang <email address hidden> wrote:
>
> Question #670297 on MadGraph5_aMC@NLO changed:
> https:/
>
> Zhang Yang gave more information on the question:
> -------- End PYTHIA Particle Data Table
> -------
>
>
> -------- LHA initialization information ------------
>
> beam kind energy pdfgrp pdfset
> A 2212 6500.000 0 247000
> B 2212 6500.000 0 247000
>
> Event weighting strategy = -4
>
> Processes, with strategy-dependent cross section info
> number xsec (pb) xerr (pb) xmax (pb)
> 1 9.0162e-01 3.2728e-03 3.0965e+00
> 3 1.5024e+00 3.1903e-03 3.0965e+00
> 2 6.9247e-01 3.5322e-03 3.0965e+00
>
> -------- End LHA initialization information --------
>
> -------- LHA event information and listing
> -------
>
> process = 3 weight = 3.0965e+00 scale = 7.3707e+02 (GeV)
> alpha_em = 7.8165e-03 alpha_strong = 9.6566e-02
>
> Participating Particles
> no id stat mothers colours p_x p_y p_z e m tau spin
> 1 1 -1 0 0 501 0 -0.000 0.000 944.379 944.379 0.000 0.000 -1.000
> 2 2 -1 0 0 502 0 0.000 -0.000 -631.452 631.452 0.000 0.000 -1.000
> 3 1000001 2 1 2 502 0 -143.474 450.543 394.595 836.048 565.418 0.000 0.000
> 4 1000002 2 1 2 501 0 143.474 -450.543 -81.668 739.783 563.059 0.000 0.000
> 5 -1000024 1 3 3 0 0 -339.970 285.919 238.697 536.022 181.696 0.000 -1.000
> 6 1000023 1 4 4 0 0 38.302 -130.509 -271.139 353.283 181.088 0.000 -1.000
> 7 2 1 4 4 501 0 105.172 -320.034 189.472 386.500 0.000 0.000 -1.000
> 8 2 1 3 3 502 0 196.496 164.624 155.898 300.026 0.000 0.000 -1.000
>
> -------- End LHA event information and listing
> -------
>
> -------- PYTHIA Info Listing -------
>
> Beam A: id = 2212, pz = 6.500e+03, e = 6.500e+03, m = 9.383e-01.
> Beam B: id = 2212, pz = -6.500e+03, e = 6.500e+03, m = 9.383e-01.
>
> In 1: id = 1, x = 1.453e-01, pdf = 0.000e+00 at Q2 = 5.433e+05.
> In 2: id = 2, x = 9.715e-02, pdf = 0.000e+00 at same Q2.
>
> Process Les Houches User Process(es) with code 9999 is 2 -> 2.
> Subprocess user process 3 with code 3 is 2 -> 2.
> It has sHat = 2.385e+06, tHat = -5.141e+05, uHat = -1.234e+06,
> pTHat = 4.728e+02, m3Hat = 5.654e+02, m4Hat = 5.631e+02,
> thetaHat = 1.113e+00, phiHat = 1.879e+00.
> alphaEM = 7.817e-03, alphaS = 9.657e-02 at Q2 = 5.433e+05.
>
> Impact parameter b = 2.553e-01 gives enhancement factor = 3.601e+00.
> Max pT scale for MPI = 1.300e+04, ISR = 1.300e+04, FSR = 1.300e+04.
> Number of MPI = 13, ISR = 33, FSRproc = 122, FSRreson = 9.
>
> -------- End PYTHIA Info Listing -------
>
> -------- PYTHIA Event Listing (hard process) -------
>
> no id name status mothers daughters colours p_x p_y p_z e m
> 0 90 (system) -11 0 0 0 0 0 0 0.000 0.000 0.000 13000.000 13000.000
> 1 2212 (p+) -12 0 0 3 0 0 0 0.000 0.000 6500.000 6500.000 0.938
> 2 2212 (p+) -12 0 0 4 0 0 0 0.000 0.000 -6500.000 6500.000 0.938
> 3 1 (d) -21 1 0 5 6 501 0 -0.000 0.000 944.379 944.379 0.000
> 4 2 (u) -21 2 0 5 6 502 0 0.000 -0.000 -631.452 631.452 0.000
> 5 1000001 (~d_L) -22 3 4 7 8 502 0 -143.474 450.543 394.595 836.048 565.418
> 6 1000002 (~u_L) -22 3 4 9 10 501 0 143.474 -450.543 -81.668 739.783 563.059
> 7 -1000024 ~chi_1- 23 5 0 0 0 0 0 -339.970 285.919 238.697 536.022 181.696
> 8 2 u 23 5 0 0 0 502 0 196.496 164.624 155.898 300.026 0.000
> 9 1000023 ~chi_20 23 6 0 0 0 0 0 38.302 -130.509 -271.139 353.283 181.088
> 10 2 u 23 6 0 0 0 501 0 105.172 -320.034 189.472 386.500 0.000
> Charge sum: 0.333 Momentum sum: 0.000 0.000 312.927 1575.831 1544.448
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Revision history for this message
|
#5 |
Hi Olivier,
Thanks very much. That solved my question.
I have another question: I want to specify the decay mode of Z boson and W boson. So I add
23:onMode = off
23:onIfAny = 11 13 15
24:onMode = off
24:onIfAny = 1 2 3 4 5
in "pythia8_card.dat", which met an error:
Merging results from the split PY8 runs...
Command "generate_events run_03" interrupted with error:
InvalidCmd : Pythia8 shower failed since it did not accept any event from the MG5aMC event file.
quit
INFO: storing files of previous run
What's the right way to specify the decay mode?
Best,
Yang
Revision history for this message
|
#6 |
I do not know.
I would refer you to the pythia8 manual for such detail.
Cheers,
Olivier
> On 19 Jun 2018, at 11:37, Zhang Yang <email address hidden> wrote:
>
> Question #670297 on MadGraph5_aMC@NLO changed:
> https:/
>
> Zhang Yang posted a new comment:
> Hi Olivier,
>
> Thanks very much. That solved my question.
>
> I have another question: I want to specify the decay mode of Z boson and
> W boson. So I add
>
> 23:onMode = off
> 23:onIfAny = 11 13 15
> 24:onMode = off
> 24:onIfAny = 1 2 3 4 5
>
> in "pythia8_card.dat", which met an error:
>
> Merging results from the split PY8 runs...
> Command "generate_events run_03" interrupted with error:
> InvalidCmd : Pythia8 shower failed since it did not accept any event from the MG5aMC event file.
> quit
> INFO: storing files of previous run
>
> What's the right way to specify the decay mode?
>
> Best,
> Yang
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Revision history for this message
|
#7 |
Hi Olivier,
I found an Error information in pythia8.log:
-------- End LHA initialization information --------
PYTHIA Error in ResonanceDecays
PYTHIA Abort from Pythia::next: processLevel failed; giving up
PYTHIA Abort from Pythia::next: reached end of Les Houches Events File
WARNING in MG5aMC_
The decay channel for id = 1000023 is read correctly in pythia8.log:
1000023 ~chi_20 2 0 0 200.00000 0.02078 199.89612 200.10388 9.48164e-12 1 1 0 0 0
0 1 1.0000000 100 23 1000022
1 1 0.0000000 100 -11 11 1000022
2 1 0.0000000 100 -13 13 1000022
3 1 0.0000000 100 -15 15 1000022
Without specify the decay mode of Z boson and W boson, it can run without error.
Many thanks,
Yang
Revision history for this message
|
#8 |
Thanks Olivier Mattelaer, that solved my question.
Revision history for this message
|
#9 |
Well,
This means that your param_card information is just wrong.
you should have three body decay here and not two body decay.
Cheers,
Olivier
> On 19 Jun 2018, at 11:43, Zhang Yang <email address hidden> wrote:
>
> Question #670297 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Zhang Yang is still having a problem:
> Hi Olivier,
>
> I found an Error information in pythia8.log:
>
> -------- End LHA initialization information --------
> PYTHIA Error in ResonanceDecays
> PYTHIA Abort from Pythia::next: processLevel failed; giving up
> PYTHIA Abort from Pythia::next: reached end of Les Houches Events File
> WARNING in MG5aMC_
>
>
> The decay channel for id = 1000023 is read correctly in pythia8.log:
>
> 1000023 ~chi_20 2 0 0 200.00000 0.02078 199.89612 200.10388 9.48164e-12 1 1 0 0 0
> 0 1 1.0000000 100 23 1000022
> 1 1 0.0000000 100 -11 11 1000022
> 2 1 0.0000000 100 -13 13 1000022
> 3 1 0.0000000 100 -15 15 1000022
>
> Without specify the decay mode of Z boson and W boson, it can run
> without error.
>
> Many thanks,
> Yang
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Revision history for this message
|
#10 |
Hi Olivier,
Sorry to bother you again.
I found the reason why the way I specify the decay mode didn't work.
In "pythia8_card.dat", I added
23:onMode = off
23:onIfAny = 11 13 15
24:onMode = off
24:onIfAny = 1 2 3 4 5
but in "Events/
24:onMode=off
24:onIfAny=1 2 3 4 5
23:onIfAny=11 13 15
23:onMode=off
And the order is important here. Could you tell me where is the codes that write "PY8Card.dat"?
Best,
Yang
Revision history for this message
|
#11 |
Hi,
I have fix the problem in a temporary way. I hard coded the order of the 4 lines in "banner.py", i.e. added
if len(visible_param) == 4 and '24:onMode' in visible_param:
after line 1995 of "banner.py" in MG5_aMC_v2_6_2.
Best,
Yang
Revision history for this message
|
#12 |
Hi Olivier and Zhang,
This issue seems to happen because the ConfigFile class (madgraph/
So that when PY8Card(
This possibly could be fixed using the below code that reorders the visible_param according to the self.user_set that could be put after "visible_param = [p for p in self if p.lower() not in self.hidden_param or p.lower() in self.user_set]"
# Reorder according to self.user_set
for self_param in visible_param:
if self_param.lower() not in self.user_set: ordered_
for user_param in self.user_set:
for self_param in self:
if user_param == self_param.lower():
break
Best regards,
Jaebak
Revision history for this message
|
#13 |
I guess this does not concern 3.4.0 since that version only support python3 where dict are ordered.
For the Long Term Stable Version (2.9.x), python2.7 is still working but officially out of support,...
But I'm pretty sure that such issue has been fixed within that version. Could you confirm that 2.9.12 does not have such issue/does not need such patch.
Cheers,
Olivier
Revision history for this message
|
#14 |
Hi Olivier,
I tried using LTS 2.9.12 with python2.7, but there still seems to be an issue.
Below [1] I show the settings that was used.
Also the fix I showed above will probably not work for all systems because self.user_set is also not ordered. One could make a simple ordered set class from OrderedDict and then used the fix I showed before. At the end of the message [2], a possible fix is shown.
Best regards,
Jaebak Kim
--------
[1]
Events are generated with the below commands.
import model heft
generate p p > h
output HToZG_pythia
launch
And then in the pythia8_card.dat, the below options are added.
25:m0 = 125.0
25:onMode = off
25:onIfMatch = 23 22
23:mMin = 50.0
23:onMode = off
23:onIfAny = 11 13 15
When generating events and looking at Events/
The below mixed-up sequence can be found
!
! Additional general parameters.
!
25:m0=125.0
25:onMode=off
25:onIfMatch=23 22
23:onIfAny=11 13 15
23:mMin=50.0
23:onMode=off
--------
[2]
In madgraph/
class OrderedSet(
def add(self, elem): self[elem] = None
def remove(self, elem): self.pop(elem, None)
def __repr__(self): return 'OrderedSet([%s])' % (', '.join(map(repr, self.keys())))
def __str__(self): return '{%s}' % (', '.join(map(repr, self.keys())))
In "class Config(dict)", change
self.user_set = set()
to
self.user_set = OrderedSet()
In class PY8Card(
# Reorder according to self.user_set
for self_param in visible_param:
if self_param.lower() not in self.user_set: ordered_
for user_param in self.user_set:
for self_param in self:
if user_param == self_param.lower():
break
Below is a "diff -u" for the changes.
--- ../Downloads/
+++ madgraph/
@@ -62,6 +62,12 @@
class UnknownType(str):
pass
+class OrderedSet(
+ def add(self, elem): self[elem] = None
+ def remove(self, elem): self.pop(elem, None)
+ def __repr__(self): return 'OrderedSet([%s])' % (', '.join(map(repr, self.keys())))
+ def __str__(self): return '{%s}' % (', '.join(map(repr, self.keys())))
+
#dict
class Banner(dict):
""" """
@@ -987,7 +993,9 @@
# Initialize it with all the default value
- self.user_set = set()
+ #self.user_set = set()
+ self.user_set = OrderedSet()
+
@@ -1916,6 +1924,18 @@
# First list the visible parameters
+
+ # Reorder according to self.user_set
+ ordered_
+ for self_param in visible_param:
+ if self_param.lower() not in self.user_set: ordered_
+ for user_param in self.user_set:
+ for self_param in self:
+ if user_param == self_param.lower():
+ ordered_
+ break
+ visible_param = ordered_
+
# Filter against list of parameters vetoed for write-out
Can you help with this problem?
Provide an answer of your own, or ask Zhang Yang for more information if necessary.