problem with bsm param_card

Asked by nabanita ganguly on 2018-01-09

Hi mg5 team,

I am having a problem when generating events using mg5 in the framework of a bsm model. The problem is as follows :

I have generated the param card for the above mentioned model using mg5 ( v2_4_2) . It looks something like this :

######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR CKMINPUTS
###################################
BLOCK CKMINPUTS #
      1 2.210000e-01 # s12
      2 4.100000e-02 # s23
      3 3.500000e-03 # s13
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      5 4.200000e+00 # mb
      6 1.750000e+02 # mtp
      23 9.118760e+01 # mz
      25 1.250000e+02 # mh
      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
      15 0.000000e+00 # ta- : 0.0
      16 0.000000e+00 # vt : 0.0
      21 0.000000e+00 # g : 0.0
      22 0.000000e+00 # a : 0.0
      24 7.995805e+01 # w+ : cw*zm
      5100001 1.685030e+03 # dd1 : mf1
      5100002 1.685030e+03 # du1 : mf1
      5100003 1.685030e+03 # ds1 : mf1
      5100004 1.685030e+03 # dc1 : mf1
      5100005 1.685030e+03 # db1 : mf1
      5100006 1.694093e+03 # dt1 : cmath.sqrt(mf1__exp__2 + tpm__exp__2)
      5100011 3.235090e+02 # e1l- : ml1
      5100012 3.235090e+02 # ve1 : ml1
      5100013 3.235090e+02 # m1l- : ml1
      5100014 3.235090e+02 # vm1 : ml1
      5100015 3.235090e+02 # t1l- : ml1
      5100016 3.235090e+02 # vt1 : ml1
      5100021 1.781390e+03 # g1 : mgl1
      5100022 2.487074e+02 # a1 : sqrt__massval1
      5100023 4.098216e+02 # z1 : cmath.sqrt(mw1__exp__2 + zm__exp__2)
      5100024 4.114062e+02 # w1+ : cmath.sqrt(massval2 + mw__exp__2)
      5100030 4.186450e+02 # h1 : cmath.sqrt(mh__exp__2 + mw1__exp__2)
      6100001 1.685030e+03 # sd1 : mf1
      6100002 1.685030e+03 # su1 : mf1
      6100003 1.685030e+03 # ss1 : mf1
      6100004 1.685030e+03 # sc1 : mf1
      6100005 1.685030e+03 # sb1 : mf1
      6100006 1.694093e+03 # st1 : cmath.sqrt(mf1__exp__2 + tpm__exp__2)
      6100011 3.235090e+02 # e1r- : ml1
      6100013 3.235090e+02 # m1r- : ml1
      6100015 3.235090e+02 # t1r- : ml1
###################################
## INFORMATION FOR NMUEDINPUTS
###################################
BLOCK NMUEDINPUTS #
      1 1.000000e-03 # r
      2 1.980000e-02 # rb
      3 6.900000e-03 # rw
      4 -3.140000e-03 # rgl
      5 1.110000e-02 # rl
      6 -2.200000e-03 # rf
      7 3.235090e+02 # ml1
      8 1.685030e+03 # mf1
      9 2.470700e+02 # mb1
      10 3.995480e+02 # mw1
      11 1.781390e+03 # mgl1
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 1.280000e+02 # aewm1
      2 1.166390e-05 # gf
      3 1.172000e-01 # as
      4 4.807600e-01 # sw
      5 9.118760e+01 # zm
      6 1.750000e+02 # tpm
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
      1 0.000000e+00 # yme
      2 0.000000e+00 # ymm
      3 1.777000e+00 # ymtau
      4 0.000000e+00 # ymdo
      5 0.000000e+00 # ymup
      6 0.000000e+00 # yms
      7 0.000000e+00 # ymc
      8 4.700000e+00 # ymb
      9 1.720000e+02 # ymt
###################################
## INFORMATION FOR YVALUES
###################################
BLOCK YVALUES #
      1 -1.000000e+00 # y1
      2 -2.000000e+00 # y2
      3 3.333333e-01 # y3
      4 1.333333e+00 # y4
      5 -6.666667e-01 # y5
###################################
## INFORMATION FOR QNUMBERS 5100011
###################################
BLOCK QNUMBERS 5100011 # e1l-
      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/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100013
###################################
BLOCK QNUMBERS 5100013 # m1l-
      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/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100015
###################################
BLOCK QNUMBERS 5100015 # t1l-
      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/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100012
###################################
BLOCK QNUMBERS 5100012 # ve1
      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 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100014
###################################
BLOCK QNUMBERS 5100014 # vm1
      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 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100016
###################################
BLOCK QNUMBERS 5100016 # vt1
      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 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 6100011
###################################
BLOCK QNUMBERS 6100011 # e1r-
      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/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 6100013
###################################
BLOCK QNUMBERS 6100013 # m1r-
      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/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 6100015
###################################
BLOCK QNUMBERS 6100015 # t1r-
      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/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 6100002
###################################
BLOCK QNUMBERS 6100002 # su1
      1 2 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 6100004
###################################
BLOCK QNUMBERS 6100004 # sc1
      1 2 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 6100006
###################################
BLOCK QNUMBERS 6100006 # st1
      1 2 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 6100001
###################################
BLOCK QNUMBERS 6100001 # sd1
      1 -1 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 6100003
###################################
BLOCK QNUMBERS 6100003 # ss1
      1 -1 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 6100005
###################################
BLOCK QNUMBERS 6100005 # sb1
      1 -1 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100002
###################################
BLOCK QNUMBERS 5100002 # du1
      1 2 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100004
###################################
BLOCK QNUMBERS 5100004 # dc1
      1 2 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100006
###################################
BLOCK QNUMBERS 5100006 # dt1
      1 2 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100001
###################################
BLOCK QNUMBERS 5100001 # dd1
      1 -1 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100003
###################################
BLOCK QNUMBERS 5100003 # ds1
      1 -1 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100005
###################################
BLOCK QNUMBERS 5100005 # db1
      1 -1 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100022
###################################
BLOCK QNUMBERS 5100022 # a1
      1 0 # 3 times electric charge
      2 3 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100023
###################################
BLOCK QNUMBERS 5100023 # z1
      1 0 # 3 times electric charge
      2 3 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100024
###################################
BLOCK QNUMBERS 5100024 # w1+
      1 3 # 3 times electric charge
      2 3 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100021
###################################
BLOCK QNUMBERS 5100021 # g1
      1 0 # 3 times electric charge
      2 3 # number of spin states (2s+1)
      3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 5100030
###################################
BLOCK QNUMBERS 5100030 # h1
      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/antiparticle distinction (0=own anti)
#
#*************************
# Decay widths *
#*************************
#
# PDG Width
DECAY 1 0.000000e+00
#
# PDG Width
DECAY 2 0.000000e+00
#
# PDG Width
DECAY 3 0.000000e+00
#
# PDG Width
DECAY 4 0.000000e+00
#
# PDG Width
DECAY 5 0.000000e+00
#
# PDG Width
DECAY 6 1.590000e+00
#
# PDG Width
DECAY 11 0.000000e+00
#
# PDG Width
DECAY 12 0.000000e+00
#
# PDG Width
DECAY 13 0.000000e+00
#
# PDG Width
DECAY 14 0.000000e+00
#
# PDG Width
DECAY 15 0.000000e+00
#
# PDG Width
DECAY 16 0.000000e+00
#
# PDG Width
DECAY 21 0.000000e+00
#
# PDG Width
DECAY 22 0.000000e+00
#
# PDG Width
DECAY 23 2.495200e+00
#
# PDG Width
DECAY 24 2.088950e+00
#
# PDG Width
DECAY 25 6.606090e-03
# BR NDA ID1 ID2 ...
   8.231008e-01 2 -5 5 # 0.00543747796387
   3.948748e-02 2 -15 15 # 0.000260857846753
   1.989074e-02 3 -24 -3 4 # 0.000131400018607
   1.972725e-02 3 -4 3 24 # 0.000130319988953
   1.962583e-02 3 -2 1 24 # 0.000129649999305
   1.961978e-02 3 -24 -1 2 # 0.00012961003246
   7.014739e-03 3 -12 11 24 # 4.63399971605e-05
   7.014739e-03 3 -14 13 24 # 4.63399971605e-05
   7.014739e-03 3 -16 15 24 # 4.63399971605e-05
   6.976896e-03 3 -24 -11 12 # 4.60900028966e-05
   6.976896e-03 3 -24 -13 14 # 4.60900028966e-05
   6.976896e-03 3 -24 -15 16 # 4.60900028966e-05
   1.941920e-03 3 -1 1 23 # 1.28284982928e-05
   1.941920e-03 3 -3 3 23 # 1.28284982928e-05
   1.623350e-03 3 -5 5 23 # 1.07239962015e-05
   1.502855e-03 3 -2 2 23 # 9.92799538695e-06
   1.502855e-03 3 -4 4 23 # 9.92799538695e-06
   1.023026e-03 3 -2 3 24 # 6.75820182834e-06
   1.022390e-03 3 -24 -1 4 # 6.7540003551e-06
   1.021118e-03 3 -24 -3 2 # 6.74559740862e-06
   1.014640e-03 3 -4 1 24 # 6.7028031576e-06
   8.864044e-04 3 -12 12 23 # 5.8556672428e-06
   8.864044e-04 3 -14 14 23 # 5.8556672428e-06
   8.864044e-04 3 -16 16 23 # 5.8556672428e-06
   4.399678e-04 3 -11 11 23 # 2.9064668839e-06
   4.399678e-04 3 -13 13 23 # 2.9064668839e-06
   4.399678e-04 3 -15 15 23 # 2.9064668839e-06
   3.203411e-16 2 22 23 # 2.1162021373e-18
#
# PDG Width
DECAY 5100001 4.474144e+04
# BR NDA ID1 ID2 ...
   9.998562e-01 4 1 5100022 5100030 5100030 # 44735.0061809
   8.814824e-05 2 -5100024 2 # 3.94387919107
   4.463076e-05 2 1 5100023 # 1.99684447069
   5.398576e-06 3 -5100024 1 24 # 0.241540064189
   3.376065e-06 3 -24 2 5100023 # 0.151050009634
   2.261357e-06 2 1 5100022 # 0.101176368534
#
# PDG Width
DECAY 5100002 4.469946e+04
# BR NDA ID1 ID2 ...
   9.998555e-01 4 2 5100022 5100030 5100030 # 44693.000928
   8.392071e-05 2 1 5100024 # 3.75121041982
   4.423294e-05 2 2 5100023 # 1.97718853221
   5.464944e-06 3 -24 2 5100024 # 0.24428004573
   4.309239e-06 2 3 5100024 # 0.192620656311
   3.467946e-06 2 2 5100022 # 0.155015313509
   3.073863e-06 3 1 24 5100023 # 0.137400016214
#
# PDG Width
DECAY 5100003 4.479544e+04
# BR NDA ID1 ID2 ...
   9.998563e-01 4 3 5100022 5100030 5100030 # 44789.0028953
   8.804198e-05 2 -5100024 4 # 3.94387923257
   4.457696e-05 2 3 5100023 # 1.99684453706
   5.413275e-06 3 -5100024 3 24 # 0.242490035466
   3.398784e-06 3 -24 4 5100023 # 0.152250024745
   2.258631e-06 2 3 5100022 # 0.101176369443
#
# PDG Width
DECAY 5100004 4.451845e+04
# BR NDA ID1 ID2 ...
   9.998552e-01 4 4 5100022 5100030 5100030 # 44512.0037284
   8.411583e-05 2 3 5100024 # 3.74470637206
   4.441279e-05 2 4 5100023 # 1.97718857098
   5.426290e-06 3 -24 4 5100024 # 0.241570020051
   4.325016e-06 2 1 5100024 # 0.192543008545
   3.482046e-06 2 4 5100022 # 0.155015290749
   3.071086e-06 3 3 24 5100023 # 0.136719988537
#
# PDG Width
DECAY 5100005 4.466436e+04
# BR NDA ID1 ID2 ...
   9.998353e-01 4 5 5100022 5100030 5100030 # 44657.0037799
   8.543410e-05 2 -5100024 6 # 3.81585939868
   4.470693e-05 2 5 5100023 # 1.99680641601
   1.565252e-05 3 -24 5 5100024 # 0.699109788187
   1.127073e-05 3 -24 6 5100023 # 0.503399942183
   5.400279e-06 3 -5100024 5 24 # 0.241200005356
   2.265217e-06 2 5 5100022 # 0.101174467566
#
# PDG Width
DECAY 5100006 2.503927e+04
# BR NDA ID1 ID2 ...
   9.996700e-01 4 6 5100022 5100030 5100030 # 25031.0070409
   1.589706e-04 2 5 5100024 # 3.98050777546
   7.727449e-05 2 6 5100023 # 1.93489681922
   3.476620e-05 3 5 24 5100023 # 0.870520268674
   2.733387e-05 3 -5100024 6 24 # 0.684420151075
   2.003715e-05 2 6 5100030 # 0.50171560888
   6.441419e-06 2 6 5100022 # 0.161288429524
   5.161494e-06 3 -24 6 5100024 # 0.129240041869
#
# PDG Width
DECAY 5100011 6.339844e-02
# BR NDA ID1 ID2 ...
   1.000000e+00 2 11 5100022 # 0.06339844
#
# PDG Width
DECAY 5100012 5.586162e-02
# BR NDA ID1 ID2 ...
   1.000000e+00 2 12 5100022 # 0.05586162
#
# PDG Width
DECAY 5100013 6.339844e-02
# BR NDA ID1 ID2 ...
   1.000000e+00 2 13 5100022 # 0.06339844
#
# PDG Width
DECAY 5100014 5.586162e-02
# BR NDA ID1 ID2 ...
   1.000000e+00 2 14 5100022 # 0.05586162
#
# PDG Width
DECAY 5100015 6.339844e-02
# BR NDA ID1 ID2 ...
   1.000000e+00 2 15 5100022 # 0.06339844
#
# PDG Width
DECAY 5100016 5.586162e-02
# BR NDA ID1 ID2 ...
   1.000000e+00 2 16 5100022 # 0.05586162
#
# PDG Width
DECAY 5100021 7.240910e-03
# BR NDA ID1 ID2 ...
   5.001029e-02 2 4 -6100004 # 0.00036212
   5.001029e-02 2 1 -6100001 # 0.00036212
   5.001029e-02 2 -1 6100001 # 0.00036212
   5.001029e-02 2 -3 6100003 # 0.00036212
   5.001029e-02 2 -3 5100003 # 0.00036212
   5.001029e-02 2 -4 6100004 # 0.00036212
   5.001029e-02 2 -4 5100004 # 0.00036212
   5.000891e-02 2 4 -5100004 # 0.00036211
   5.000891e-02 2 3 -5100003 # 0.00036211
   5.000891e-02 2 3 -6100003 # 0.00036211
   5.000891e-02 2 2 -5100002 # 0.00036211
   5.000891e-02 2 -1 5100001 # 0.00036211
   5.000891e-02 2 -2 6100002 # 0.00036211
   5.000891e-02 2 -2 5100002 # 0.00036211
   5.000753e-02 2 2 -6100002 # 0.0003621
   5.000615e-02 2 1 -5100001 # 0.00036209
   4.996333e-02 2 5100005 -5 # 0.00036178
   4.996333e-02 2 5 -5100005 # 0.00036178
   4.996333e-02 2 5 -6100005 # 0.00036178
   4.996195e-02 2 6100005 -5 # 0.00036177
#
# PDG Width
DECAY 5100022 0.000000e+00
#
# PDG Width
DECAY 5100023 1.400207e+00
# BR NDA ID1 ID2 ...
   8.389950e-02 2 -12 5100012 # 0.117476667196
   8.389950e-02 2 -14 5100014 # 0.117476667196
   8.389950e-02 2 -16 5100016 # 0.117476667196
   8.389950e-02 2 -5100012 12 # 0.117476667196
   8.389950e-02 2 -5100014 14 # 0.117476667196
   8.389950e-02 2 -5100016 16 # 0.117476667196
   8.115700e-02 2 -11 5100011 # 0.113636599499
   8.115700e-02 2 -13 5100013 # 0.113636599499
   8.115700e-02 2 -15 5100015 # 0.113636599499
   8.115700e-02 2 -5100011 11 # 0.113636599499
   8.115700e-02 2 -5100013 13 # 0.113636599499
   8.115700e-02 2 -5100015 15 # 0.113636599499
   9.524322e-03 2 25 5100022 # 0.0133360223347
   2.278563e-05 2 -11 6100011 # 3.19045986254e-05
   2.278563e-05 2 -13 6100013 # 3.19045986254e-05
   2.278563e-05 2 -15 6100015 # 3.19045986254e-05
   2.278563e-05 2 -6100011 11 # 3.19045986254e-05
   2.278563e-05 2 -6100013 13 # 3.19045986254e-05
   2.278563e-05 2 -6100015 15 # 3.19045986254e-05
#
# PDG Width
DECAY 5100024 1.430156e+00
# BR NDA ID1 ID2 ...
   1.661039e-01 2 -11 5100012 # 0.237554489208
   1.661039e-01 2 -13 5100014 # 0.237554489208
   1.661039e-01 2 -15 5100016 # 0.237554489208
   1.661039e-01 2 -5100011 12 # 0.237554489208
   1.661039e-01 2 -5100013 14 # 0.237554489208
   1.661039e-01 2 -5100015 16 # 0.237554489208
   3.376304e-03 2 24 5100022 # 0.00482863999327
#
# PDG Width
DECAY 5100030 1.647530e-01
# BR NDA ID1 ID2 ...
   9.981434e-01 2 23 5100022 # 0.16444711958
   4.625648e-04 2 -15 6100015 # 7.62089384944e-05
   4.625648e-04 2 -15 5100015 # 7.62089384944e-05
   4.625648e-04 2 -5100015 15 # 7.62089384944e-05
   4.625648e-04 2 -6100015 15 # 7.62089384944e-05
   6.388959e-06 3 -24 24 5100022 # 1.05260016213e-06
#
# PDG Width
DECAY 6100001 5.066624e-01
# BR NDA ID1 ID2 ...
   1.000000e+00 2 1 5100022 # 0.5066624
#
# PDG Width
DECAY 6100002 2.026648e+00
# BR NDA ID1 ID2 ...
   1.000000e+00 2 2 5100022 # 2.026648
#
# PDG Width
DECAY 6100003 5.066624e-01
# BR NDA ID1 ID2 ...
   1.000000e+00 2 3 5100022 # 0.5066624
#
# PDG Width
DECAY 6100004 2.026648e+00
# BR NDA ID1 ID2 ...
   1.000000e+00 2 4 5100022 # 2.026648
#
# PDG Width
DECAY 6100005 5.066530e-01
# BR NDA ID1 ID2 ...
   1.000000e+00 2 5 5100022 # 0.506653
#
# PDG Width
DECAY 6100006 1.222451e+02
# BR NDA ID1 ID2 ...
   9.796716e-01 4 6 5100022 5100030 5100030 # 119.760052709
   1.622425e-02 2 6 5100022 # 1.98333506367
   4.104179e-03 2 6 5100030 # 0.501715772273
#
# PDG Width
DECAY 6100011 2.382817e-01
# BR NDA ID1 ID2 ...
   1.000000e+00 2 11 5100022 # 0.2382817
#
# PDG Width
DECAY 6100013 2.382817e-01
# BR NDA ID1 ID2 ...
   1.000000e+00 2 13 5100022 # 0.2382817
#
# PDG Width
DECAY 6100015 2.382817e-01
# BR NDA ID1 ID2 ...
   1.000000e+00 2 15 5100022 # 0.2382817

My process of interest is (w1+(-) z1) production. After generating the events I pass them through pythia inside mg5 for hadronisation. As you can see from the param_card, w1+(-) and z1 will decay dominantly into leptons. But when I study the .hep files using madanalysis, I find that there is no leptonic event!! Leptons are given 10 GeV PT cut only at the analysis level. Looking at pythia.log file which enlists the history of first 10 events, I find out that z1 is mostly being decayed through the invisible channel although leptonic and invisible modes have got almost equal BRs.

I repeat the same exercise (i.e. same process with same param_card) but this time with mg5 v2_5_5. Analysing .hepmc file using madanalysis , I am now getting nonzero leptonic events.

I have done several analyses in the same model with different param_cards using v2_4_2 before but this kind of problem did not arise there . Apparently, I do not see any inconsistency in this param_card.

Could you please shed some light on it?

Best
Nabanita

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2018-01-09
Last reply:
2018-01-11

Hi,

Which programs is used for the decay?
is it MadSpin, MG5aMC (via the decay chain syntax) or pythia?
Did you use the same version of pythia for 2.4.2 and 2.5.5?
(note that if you use the MG5aMC to pythia interface, the one of 2.4.2 is with pythia6 --not supported anymore neither by us or by pythia author--, while the 2.5.5 is using pythia8)

Cheers,

Olivier

> On 9 Jan 2018, at 12:47, nabanita ganguly <email address hidden> wrote:
>
> New question #662728 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/662728
>
> Hi mg5 team,
>
>
> I am having a problem when generating events using mg5 in the framework of a bsm model. The problem is as follows :
>
> I have generated the param card for the above mentioned model using mg5 ( v2_4_2) . It looks something like this :
>
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
> ######################################################################
> ###################################
> ## INFORMATION FOR CKMINPUTS
> ###################################
> BLOCK CKMINPUTS #
> 1 2.210000e-01 # s12
> 2 4.100000e-02 # s23
> 3 3.500000e-03 # s13
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> BLOCK MASS #
> 5 4.200000e+00 # mb
> 6 1.750000e+02 # mtp
> 23 9.118760e+01 # mz
> 25 1.250000e+02 # mh
> 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
> 15 0.000000e+00 # ta- : 0.0
> 16 0.000000e+00 # vt : 0.0
> 21 0.000000e+00 # g : 0.0
> 22 0.000000e+00 # a : 0.0
> 24 7.995805e+01 # w+ : cw*zm
> 5100001 1.685030e+03 # dd1 : mf1
> 5100002 1.685030e+03 # du1 : mf1
> 5100003 1.685030e+03 # ds1 : mf1
> 5100004 1.685030e+03 # dc1 : mf1
> 5100005 1.685030e+03 # db1 : mf1
> 5100006 1.694093e+03 # dt1 : cmath.sqrt(mf1__exp__2 + tpm__exp__2)
> 5100011 3.235090e+02 # e1l- : ml1
> 5100012 3.235090e+02 # ve1 : ml1
> 5100013 3.235090e+02 # m1l- : ml1
> 5100014 3.235090e+02 # vm1 : ml1
> 5100015 3.235090e+02 # t1l- : ml1
> 5100016 3.235090e+02 # vt1 : ml1
> 5100021 1.781390e+03 # g1 : mgl1
> 5100022 2.487074e+02 # a1 : sqrt__massval1
> 5100023 4.098216e+02 # z1 : cmath.sqrt(mw1__exp__2 + zm__exp__2)
> 5100024 4.114062e+02 # w1+ : cmath.sqrt(massval2 + mw__exp__2)
> 5100030 4.186450e+02 # h1 : cmath.sqrt(mh__exp__2 + mw1__exp__2)
> 6100001 1.685030e+03 # sd1 : mf1
> 6100002 1.685030e+03 # su1 : mf1
> 6100003 1.685030e+03 # ss1 : mf1
> 6100004 1.685030e+03 # sc1 : mf1
> 6100005 1.685030e+03 # sb1 : mf1
> 6100006 1.694093e+03 # st1 : cmath.sqrt(mf1__exp__2 + tpm__exp__2)
> 6100011 3.235090e+02 # e1r- : ml1
> 6100013 3.235090e+02 # m1r- : ml1
> 6100015 3.235090e+02 # t1r- : ml1
> ###################################
> ## INFORMATION FOR NMUEDINPUTS
> ###################################
> BLOCK NMUEDINPUTS #
> 1 1.000000e-03 # r
> 2 1.980000e-02 # rb
> 3 6.900000e-03 # rw
> 4 -3.140000e-03 # rgl
> 5 1.110000e-02 # rl
> 6 -2.200000e-03 # rf
> 7 3.235090e+02 # ml1
> 8 1.685030e+03 # mf1
> 9 2.470700e+02 # mb1
> 10 3.995480e+02 # mw1
> 11 1.781390e+03 # mgl1
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> BLOCK SMINPUTS #
> 1 1.280000e+02 # aewm1
> 2 1.166390e-05 # gf
> 3 1.172000e-01 # as
> 4 4.807600e-01 # sw
> 5 9.118760e+01 # zm
> 6 1.750000e+02 # tpm
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> BLOCK YUKAWA #
> 1 0.000000e+00 # yme
> 2 0.000000e+00 # ymm
> 3 1.777000e+00 # ymtau
> 4 0.000000e+00 # ymdo
> 5 0.000000e+00 # ymup
> 6 0.000000e+00 # yms
> 7 0.000000e+00 # ymc
> 8 4.700000e+00 # ymb
> 9 1.720000e+02 # ymt
> ###################################
> ## INFORMATION FOR YVALUES
> ###################################
> BLOCK YVALUES #
> 1 -1.000000e+00 # y1
> 2 -2.000000e+00 # y2
> 3 3.333333e-01 # y3
> 4 1.333333e+00 # y4
> 5 -6.666667e-01 # y5
> ###################################
> ## INFORMATION FOR QNUMBERS 5100011
> ###################################
> BLOCK QNUMBERS 5100011 # e1l-
> 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/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100013
> ###################################
> BLOCK QNUMBERS 5100013 # m1l-
> 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/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100015
> ###################################
> BLOCK QNUMBERS 5100015 # t1l-
> 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/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100012
> ###################################
> BLOCK QNUMBERS 5100012 # ve1
> 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 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100014
> ###################################
> BLOCK QNUMBERS 5100014 # vm1
> 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 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100016
> ###################################
> BLOCK QNUMBERS 5100016 # vt1
> 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 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 6100011
> ###################################
> BLOCK QNUMBERS 6100011 # e1r-
> 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/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 6100013
> ###################################
> BLOCK QNUMBERS 6100013 # m1r-
> 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/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 6100015
> ###################################
> BLOCK QNUMBERS 6100015 # t1r-
> 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/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 6100002
> ###################################
> BLOCK QNUMBERS 6100002 # su1
> 1 2 # 3 times electric charge
> 2 2 # number of spin states (2s+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 6100004
> ###################################
> BLOCK QNUMBERS 6100004 # sc1
> 1 2 # 3 times electric charge
> 2 2 # number of spin states (2s+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 6100006
> ###################################
> BLOCK QNUMBERS 6100006 # st1
> 1 2 # 3 times electric charge
> 2 2 # number of spin states (2s+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 6100001
> ###################################
> BLOCK QNUMBERS 6100001 # sd1
> 1 -1 # 3 times electric charge
> 2 2 # number of spin states (2s+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 6100003
> ###################################
> BLOCK QNUMBERS 6100003 # ss1
> 1 -1 # 3 times electric charge
> 2 2 # number of spin states (2s+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 6100005
> ###################################
> BLOCK QNUMBERS 6100005 # sb1
> 1 -1 # 3 times electric charge
> 2 2 # number of spin states (2s+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100002
> ###################################
> BLOCK QNUMBERS 5100002 # du1
> 1 2 # 3 times electric charge
> 2 2 # number of spin states (2s+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100004
> ###################################
> BLOCK QNUMBERS 5100004 # dc1
> 1 2 # 3 times electric charge
> 2 2 # number of spin states (2s+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100006
> ###################################
> BLOCK QNUMBERS 5100006 # dt1
> 1 2 # 3 times electric charge
> 2 2 # number of spin states (2s+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100001
> ###################################
> BLOCK QNUMBERS 5100001 # dd1
> 1 -1 # 3 times electric charge
> 2 2 # number of spin states (2s+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100003
> ###################################
> BLOCK QNUMBERS 5100003 # ds1
> 1 -1 # 3 times electric charge
> 2 2 # number of spin states (2s+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100005
> ###################################
> BLOCK QNUMBERS 5100005 # db1
> 1 -1 # 3 times electric charge
> 2 2 # number of spin states (2s+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100022
> ###################################
> BLOCK QNUMBERS 5100022 # a1
> 1 0 # 3 times electric charge
> 2 3 # number of spin states (2s+1)
> 3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 0 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100023
> ###################################
> BLOCK QNUMBERS 5100023 # z1
> 1 0 # 3 times electric charge
> 2 3 # number of spin states (2s+1)
> 3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 0 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100024
> ###################################
> BLOCK QNUMBERS 5100024 # w1+
> 1 3 # 3 times electric charge
> 2 3 # number of spin states (2s+1)
> 3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100021
> ###################################
> BLOCK QNUMBERS 5100021 # g1
> 1 0 # 3 times electric charge
> 2 3 # number of spin states (2s+1)
> 3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 0 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 5100030
> ###################################
> BLOCK QNUMBERS 5100030 # h1
> 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/antiparticle distinction (0=own anti)
> #
> #*************************
> # Decay widths *
> #*************************
> #
> # PDG Width
> DECAY 1 0.000000e+00
> #
> # PDG Width
> DECAY 2 0.000000e+00
> #
> # PDG Width
> DECAY 3 0.000000e+00
> #
> # PDG Width
> DECAY 4 0.000000e+00
> #
> # PDG Width
> DECAY 5 0.000000e+00
> #
> # PDG Width
> DECAY 6 1.590000e+00
> #
> # PDG Width
> DECAY 11 0.000000e+00
> #
> # PDG Width
> DECAY 12 0.000000e+00
> #
> # PDG Width
> DECAY 13 0.000000e+00
> #
> # PDG Width
> DECAY 14 0.000000e+00
> #
> # PDG Width
> DECAY 15 0.000000e+00
> #
> # PDG Width
> DECAY 16 0.000000e+00
> #
> # PDG Width
> DECAY 21 0.000000e+00
> #
> # PDG Width
> DECAY 22 0.000000e+00
> #
> # PDG Width
> DECAY 23 2.495200e+00
> #
> # PDG Width
> DECAY 24 2.088950e+00
> #
> # PDG Width
> DECAY 25 6.606090e-03
> # BR NDA ID1 ID2 ...
> 8.231008e-01 2 -5 5 # 0.00543747796387
> 3.948748e-02 2 -15 15 # 0.000260857846753
> 1.989074e-02 3 -24 -3 4 # 0.000131400018607
> 1.972725e-02 3 -4 3 24 # 0.000130319988953
> 1.962583e-02 3 -2 1 24 # 0.000129649999305
> 1.961978e-02 3 -24 -1 2 # 0.00012961003246
> 7.014739e-03 3 -12 11 24 # 4.63399971605e-05
> 7.014739e-03 3 -14 13 24 # 4.63399971605e-05
> 7.014739e-03 3 -16 15 24 # 4.63399971605e-05
> 6.976896e-03 3 -24 -11 12 # 4.60900028966e-05
> 6.976896e-03 3 -24 -13 14 # 4.60900028966e-05
> 6.976896e-03 3 -24 -15 16 # 4.60900028966e-05
> 1.941920e-03 3 -1 1 23 # 1.28284982928e-05
> 1.941920e-03 3 -3 3 23 # 1.28284982928e-05
> 1.623350e-03 3 -5 5 23 # 1.07239962015e-05
> 1.502855e-03 3 -2 2 23 # 9.92799538695e-06
> 1.502855e-03 3 -4 4 23 # 9.92799538695e-06
> 1.023026e-03 3 -2 3 24 # 6.75820182834e-06
> 1.022390e-03 3 -24 -1 4 # 6.7540003551e-06
> 1.021118e-03 3 -24 -3 2 # 6.74559740862e-06
> 1.014640e-03 3 -4 1 24 # 6.7028031576e-06
> 8.864044e-04 3 -12 12 23 # 5.8556672428e-06
> 8.864044e-04 3 -14 14 23 # 5.8556672428e-06
> 8.864044e-04 3 -16 16 23 # 5.8556672428e-06
> 4.399678e-04 3 -11 11 23 # 2.9064668839e-06
> 4.399678e-04 3 -13 13 23 # 2.9064668839e-06
> 4.399678e-04 3 -15 15 23 # 2.9064668839e-06
> 3.203411e-16 2 22 23 # 2.1162021373e-18
> #
> # PDG Width
> DECAY 5100001 4.474144e+04
> # BR NDA ID1 ID2 ...
> 9.998562e-01 4 1 5100022 5100030 5100030 # 44735.0061809
> 8.814824e-05 2 -5100024 2 # 3.94387919107
> 4.463076e-05 2 1 5100023 # 1.99684447069
> 5.398576e-06 3 -5100024 1 24 # 0.241540064189
> 3.376065e-06 3 -24 2 5100023 # 0.151050009634
> 2.261357e-06 2 1 5100022 # 0.101176368534
> #
> # PDG Width
> DECAY 5100002 4.469946e+04
> # BR NDA ID1 ID2 ...
> 9.998555e-01 4 2 5100022 5100030 5100030 # 44693.000928
> 8.392071e-05 2 1 5100024 # 3.75121041982
> 4.423294e-05 2 2 5100023 # 1.97718853221
> 5.464944e-06 3 -24 2 5100024 # 0.24428004573
> 4.309239e-06 2 3 5100024 # 0.192620656311
> 3.467946e-06 2 2 5100022 # 0.155015313509
> 3.073863e-06 3 1 24 5100023 # 0.137400016214
> #
> # PDG Width
> DECAY 5100003 4.479544e+04
> # BR NDA ID1 ID2 ...
> 9.998563e-01 4 3 5100022 5100030 5100030 # 44789.0028953
> 8.804198e-05 2 -5100024 4 # 3.94387923257
> 4.457696e-05 2 3 5100023 # 1.99684453706
> 5.413275e-06 3 -5100024 3 24 # 0.242490035466
> 3.398784e-06 3 -24 4 5100023 # 0.152250024745
> 2.258631e-06 2 3 5100022 # 0.101176369443
> #
> # PDG Width
> DECAY 5100004 4.451845e+04
> # BR NDA ID1 ID2 ...
> 9.998552e-01 4 4 5100022 5100030 5100030 # 44512.0037284
> 8.411583e-05 2 3 5100024 # 3.74470637206
> 4.441279e-05 2 4 5100023 # 1.97718857098
> 5.426290e-06 3 -24 4 5100024 # 0.241570020051
> 4.325016e-06 2 1 5100024 # 0.192543008545
> 3.482046e-06 2 4 5100022 # 0.155015290749
> 3.071086e-06 3 3 24 5100023 # 0.136719988537
> #
> # PDG Width
> DECAY 5100005 4.466436e+04
> # BR NDA ID1 ID2 ...
> 9.998353e-01 4 5 5100022 5100030 5100030 # 44657.0037799
> 8.543410e-05 2 -5100024 6 # 3.81585939868
> 4.470693e-05 2 5 5100023 # 1.99680641601
> 1.565252e-05 3 -24 5 5100024 # 0.699109788187
> 1.127073e-05 3 -24 6 5100023 # 0.503399942183
> 5.400279e-06 3 -5100024 5 24 # 0.241200005356
> 2.265217e-06 2 5 5100022 # 0.101174467566
> #
> # PDG Width
> DECAY 5100006 2.503927e+04
> # BR NDA ID1 ID2 ...
> 9.996700e-01 4 6 5100022 5100030 5100030 # 25031.0070409
> 1.589706e-04 2 5 5100024 # 3.98050777546
> 7.727449e-05 2 6 5100023 # 1.93489681922
> 3.476620e-05 3 5 24 5100023 # 0.870520268674
> 2.733387e-05 3 -5100024 6 24 # 0.684420151075
> 2.003715e-05 2 6 5100030 # 0.50171560888
> 6.441419e-06 2 6 5100022 # 0.161288429524
> 5.161494e-06 3 -24 6 5100024 # 0.129240041869
> #
> # PDG Width
> DECAY 5100011 6.339844e-02
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 11 5100022 # 0.06339844
> #
> # PDG Width
> DECAY 5100012 5.586162e-02
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 12 5100022 # 0.05586162
> #
> # PDG Width
> DECAY 5100013 6.339844e-02
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 13 5100022 # 0.06339844
> #
> # PDG Width
> DECAY 5100014 5.586162e-02
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 14 5100022 # 0.05586162
> #
> # PDG Width
> DECAY 5100015 6.339844e-02
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 15 5100022 # 0.06339844
> #
> # PDG Width
> DECAY 5100016 5.586162e-02
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 16 5100022 # 0.05586162
> #
> # PDG Width
> DECAY 5100021 7.240910e-03
> # BR NDA ID1 ID2 ...
> 5.001029e-02 2 4 -6100004 # 0.00036212
> 5.001029e-02 2 1 -6100001 # 0.00036212
> 5.001029e-02 2 -1 6100001 # 0.00036212
> 5.001029e-02 2 -3 6100003 # 0.00036212
> 5.001029e-02 2 -3 5100003 # 0.00036212
> 5.001029e-02 2 -4 6100004 # 0.00036212
> 5.001029e-02 2 -4 5100004 # 0.00036212
> 5.000891e-02 2 4 -5100004 # 0.00036211
> 5.000891e-02 2 3 -5100003 # 0.00036211
> 5.000891e-02 2 3 -6100003 # 0.00036211
> 5.000891e-02 2 2 -5100002 # 0.00036211
> 5.000891e-02 2 -1 5100001 # 0.00036211
> 5.000891e-02 2 -2 6100002 # 0.00036211
> 5.000891e-02 2 -2 5100002 # 0.00036211
> 5.000753e-02 2 2 -6100002 # 0.0003621
> 5.000615e-02 2 1 -5100001 # 0.00036209
> 4.996333e-02 2 5100005 -5 # 0.00036178
> 4.996333e-02 2 5 -5100005 # 0.00036178
> 4.996333e-02 2 5 -6100005 # 0.00036178
> 4.996195e-02 2 6100005 -5 # 0.00036177
> #
> # PDG Width
> DECAY 5100022 0.000000e+00
> #
> # PDG Width
> DECAY 5100023 1.400207e+00
> # BR NDA ID1 ID2 ...
> 8.389950e-02 2 -12 5100012 # 0.117476667196
> 8.389950e-02 2 -14 5100014 # 0.117476667196
> 8.389950e-02 2 -16 5100016 # 0.117476667196
> 8.389950e-02 2 -5100012 12 # 0.117476667196
> 8.389950e-02 2 -5100014 14 # 0.117476667196
> 8.389950e-02 2 -5100016 16 # 0.117476667196
> 8.115700e-02 2 -11 5100011 # 0.113636599499
> 8.115700e-02 2 -13 5100013 # 0.113636599499
> 8.115700e-02 2 -15 5100015 # 0.113636599499
> 8.115700e-02 2 -5100011 11 # 0.113636599499
> 8.115700e-02 2 -5100013 13 # 0.113636599499
> 8.115700e-02 2 -5100015 15 # 0.113636599499
> 9.524322e-03 2 25 5100022 # 0.0133360223347
> 2.278563e-05 2 -11 6100011 # 3.19045986254e-05
> 2.278563e-05 2 -13 6100013 # 3.19045986254e-05
> 2.278563e-05 2 -15 6100015 # 3.19045986254e-05
> 2.278563e-05 2 -6100011 11 # 3.19045986254e-05
> 2.278563e-05 2 -6100013 13 # 3.19045986254e-05
> 2.278563e-05 2 -6100015 15 # 3.19045986254e-05
> #
> # PDG Width
> DECAY 5100024 1.430156e+00
> # BR NDA ID1 ID2 ...
> 1.661039e-01 2 -11 5100012 # 0.237554489208
> 1.661039e-01 2 -13 5100014 # 0.237554489208
> 1.661039e-01 2 -15 5100016 # 0.237554489208
> 1.661039e-01 2 -5100011 12 # 0.237554489208
> 1.661039e-01 2 -5100013 14 # 0.237554489208
> 1.661039e-01 2 -5100015 16 # 0.237554489208
> 3.376304e-03 2 24 5100022 # 0.00482863999327
> #
> # PDG Width
> DECAY 5100030 1.647530e-01
> # BR NDA ID1 ID2 ...
> 9.981434e-01 2 23 5100022 # 0.16444711958
> 4.625648e-04 2 -15 6100015 # 7.62089384944e-05
> 4.625648e-04 2 -15 5100015 # 7.62089384944e-05
> 4.625648e-04 2 -5100015 15 # 7.62089384944e-05
> 4.625648e-04 2 -6100015 15 # 7.62089384944e-05
> 6.388959e-06 3 -24 24 5100022 # 1.05260016213e-06
> #
> # PDG Width
> DECAY 6100001 5.066624e-01
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 1 5100022 # 0.5066624
> #
> # PDG Width
> DECAY 6100002 2.026648e+00
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 2 5100022 # 2.026648
> #
> # PDG Width
> DECAY 6100003 5.066624e-01
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 3 5100022 # 0.5066624
> #
> # PDG Width
> DECAY 6100004 2.026648e+00
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 4 5100022 # 2.026648
> #
> # PDG Width
> DECAY 6100005 5.066530e-01
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 5 5100022 # 0.506653
> #
> # PDG Width
> DECAY 6100006 1.222451e+02
> # BR NDA ID1 ID2 ...
> 9.796716e-01 4 6 5100022 5100030 5100030 # 119.760052709
> 1.622425e-02 2 6 5100022 # 1.98333506367
> 4.104179e-03 2 6 5100030 # 0.501715772273
> #
> # PDG Width
> DECAY 6100011 2.382817e-01
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 11 5100022 # 0.2382817
> #
> # PDG Width
> DECAY 6100013 2.382817e-01
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 13 5100022 # 0.2382817
> #
> # PDG Width
> DECAY 6100015 2.382817e-01
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 15 5100022 # 0.2382817
>
>
> My process of interest is (w1+(-) z1) production. After generating the events I pass them through pythia inside mg5 for hadronisation. As you can see from the param_card, w1+(-) and z1 will decay dominantly into leptons. But when I study the .hep files using madanalysis, I find that there is no leptonic event!! Leptons are given 10 GeV PT cut only at the analysis level. Looking at pythia.log file which enlists the history of first 10 events, I find out that z1 is mostly being decayed through the invisible channel although leptonic and invisible modes have got almost equal BRs.
>
> I repeat the same exercise (i.e. same process with same param_card) but this time with mg5 v2_5_5. Analysing .hepmc file using madanalysis , I am now getting nonzero leptonic events.
>
> I have done several analyses in the same model with different param_cards using v2_4_2 before but this kind of problem did not arise there . Apparently, I do not see any inconsistency in this param_card.
>
> Could you please shed some light on it?
>
> Best
> Nabanita
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier,

Thank you for your message. Yes it is pythia6 (the one used in mg5 to pythia interface).

What bothers me is that I have previously used the same version of mg5(2.4.2) for the same model. But this kind of a problem never arose there.

 So my initial guess was there must be something wrong with the param_card. But since v.2.5.5 is working fine this is most probably not the case here.

So my question is : is it really something to do with the card or may be a pythia problem (which is also unlikely as pythia6 inside mg5 worked fine for other points) or I am missing something?

Best
Nabanita

Sorry, I'm confused by your answer.

Can you confirm that you are using pythia6 for both versions of MG5aMC?

Since you are asking PY6 to make the decay (if I understand correctly), you should check the difference between the version of PY6 that you are using. It is likely that you are not using the exact version of py6 in both cases. (i.e. the py6 has been updated correctly).

It sounds unlikely that this is related to our code so can not really help you more than this.

Cheers,

Olivier

PS: Note that PY6 is not supported anymore neither for the interface between MG5aMC and PY6 (since a year), neither by the PY6 author (since three years).

Can you help with this problem?

Provide an answer of your own, or ask nabanita ganguly for more information if necessary.

To post a message you must log in.