Issues with dipole operator from FeynRules
Dear MadGraph team,
I have taken an existing FeynRules model (heavyN) and added dim-5 effective operators. Everything works fine until I add dipole operators:
****
convert model /home/jjones/
import model heavyNeff5_all_UFO --modelname
compute_widths n4 --body_decay=3
Please note that the automatic computation of the width is only valid in narrow-width approximation and at tree-level.
INFO: load particles
INFO: load vertices
INFO: Get two body decay from FeynRules formula
Command "compute_widths n4 --body_decay=3 " interrupted in sub-command:
"compute_widths n4 --body_decay=3" with error:
SyntaxError : invalid syntax (<string>, line 1)
Please report this bug on https:/
More information is found in 'MG5_debug'.
Please attach this file to your report.
****
The MG5_debug is not very helpful. It includes:
****
Traceback (most recent call last):
File "/home/
return self.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
return self.cmd.
File "/home/
value = eval(expr,
File "<string>", line 1
****
Then comes a long line of code involving the coefficient of the dipole operator, I copy a short part of it:
****
((mdl_mN4**2 - mdl_MZ*
****
It then finishes with:
****
SyntaxError: invalid syntax
Related File: <string>
****
And I am stumped. I have adapted the dipole operator from the BSMC_v2.2.fr file of the BSM Characterisation model in the FeynRules page, so I really do not know what I am doing wrong.
This might be a long shot, but could it be a MadGraph issue related to the Majorana nature of the neutrinos?
Many thanks for your attention.
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Olivier Mattelaer
- Solved:
- 2021-03-25
- Last query:
- 2021-03-25
- Last reply:
- 2021-03-23
Hi,
This is a FeynRules issue. Your file decays.py (part of the UFO) is clearly not correctly formatted (the expression you quote is not a valid python syntax as it should be).
So I'm not really able to help you on that. What you can do is to remove the full file decays.py (since it is optional). This will slow down the computation of the width but it should go trough. The second solution is of course to ask FeynRules author for support.
Cheers,
Olivier
Joel Jones (n-jones-j) said : | #2 |
Thanks Olivier Mattelaer, that solved my question.
Joel Jones (n-jones-j) said : | #3 |
Thanks for the info!