error with topbsm feynrules model

Asked by Minho Kim on 2017-09-13

Dear mg5 team,

Hi, I am minho kim, graduate student in Postech, korea.

recently when i use topBSM model with madgraph5 provided in http://feynrules.irmp.ucl.ac.be/wiki/topBSM,

mg5 shows some error message like

--------------------------------------------------------------------
WARNING: Failed to update dependent parameter. This might create trouble for external program (like MadSpin/shower/...)
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
Traceback (most recent call last):
  File "write_param_card.py", line 205, in <module>
    ParamCardWriter('./param_card.dat', generic=True)
  File "write_param_card.py", line 30, in __init__
    self.write_card(list_of_parameters)
  File "write_param_card.py", line 90, in write_card
    self.write_dep_param_block(lhablock)
  File "write_param_card.py", line 123, in write_dep_param_block
    exec("%s = %s" % (parameter.name, parameter.value))
  File "<string>", line 1, in <module>
ZeroDivisionError: complex division by zero.
---------------------------------------------------------------------------

To trace above error, i add 'print(parameter.name,parameter.value)' in line 123 on 'write_param_card.py',
Then it gives

------------------------------------------------------------------------
.......
('o0fusionAxial', '(G**2*MT**2*o0axial*(-(cmath.pi*complex(0,1)) + cmath.log((1 + cmath.sqrt(1 - (4*MT**2)/MO0**2))/(1 - cmath.sqrt(1 - (4*MT**2)/MO0**2))))**2)/(8.*cmath.pi**2*MO0**2*vev)')
Traceback (most recent call last):
  File "write_param_card.py", line 206, in <module>
    ParamCardWriter('./param_card.dat', generic=True)
  File "write_param_card.py", line 30, in __init__
    self.write_card(list_of_parameters)
  File "write_param_card.py", line 90, in write_card
    self.write_dep_param_block(lhablock)
  File "write_param_card.py", line 124, in write_dep_param_block
    exec("%s = %s" % (parameter.name, parameter.value))
  File "<string>", line 1, in <module>
ZeroDivisionError: complex division by zero
--------------------------------------------------------------------

the method to make UFO file is as follows:

---------------------------------------------------------------------
$FeynRulesPath = SetDirectory["~/Programs/feynrules-current"]; << FeynRules`;
LoadModel["SM.fr", "topBSM.fr"];
FeynmanGauge = False;

WriteUFO[
  LGauge,
  LHiggs,
  LFermions,
  LYukawa,
  LGhost,
  LS0, LO0, LS1, LO1, AddDecays -> False];
---------------------------------------------------------------------

mg5 version is 2.5.5
mathematica version is 11.1 (linux)
feyrules version is 2.3.29

Interesting is even mg5 shows previous error message, it seems to work fine at least process like p p > s0 | s1 | o0 | o1 > t t~.

So why that error happens?

Thanks

Minho Kim

Question information

Language:
English Edit question
Status:
Open
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2017-09-13
Last reply:

Can you help with this problem?

Provide an answer of your own, or ask Minho Kim for more information if necessary.

To post a message you must log in.