compute_widths command not executed in MadGraph 2.7.2

Asked by Jochen Jens Heinrich on 2020-07-23

Hello,

I am trying to run MadGraph 2.7.2 with a process card that looks like this

import model MODEL
[...some definitions...]
generate PROCESS
launch
set WDP0 Auto
set WDP Auto
set Wro Auto
set Wrop Auto
set lhc 13
set pionmass 300.0
set fermioneta 0.45
set xi 0.173611111111
update dependent
compute_widths dp+ dp0 --body_decay=3 --min_br=0.001
compute_widths rho0 rho+ --body_decay=2 --min_br=0.001

Running this yields the following log output:

[...]
set xi 0.173611111111
INFO: modify param_card information BLOCK frblock with id (5,) set to 0.173611111111
update dependent
INFO: Update the dependent parameter of the param_card.dat
INFO: For consistency, the mass of particle 4900113 (rho0) is changed to 666.666666667.
Command "generate_events run_01" interrupted with error:
ValueError : invalid literal for int() with base 10: 'AutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAutoAuto'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'ME5_debug'.
Please attach this file to your report.
INFO:
INFO:
command not executed: compute_widths dp+ dp0 --body_decay=3 --min_br=0.001
command not executed: compute_widths rho0 rho+ --body_decay=2 --min_br=0.001
quit
INFO:
quit

I am a bit surprised at this, since I have used the exact same process card successfully in MadGraph 2.5.5. What changed in between those two versions that I am missing? How can I fix the error and compute the widths?

Thanks for all the help,
Jochen

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2020-07-23
Last reply:
2020-07-29

looks like the "update dependent" does not auto-trigger first the compuation of the width that you set as "Auto"
The simplest solution is to force the direct computation of those width via the compute_widths command

launch
set lhc 13
set pionmass 300.0
set fermioneta 0.45
set xi 0.173611111111
compute_widths dp+ dp0
compute_widths rho0 rho+
update dependent
compute_widths dp+ dp0 --body_decay=3 --min_br=0.001
compute_widths rho0 rho+ --body_decay=2 --min_br=0.001

So this should work and should be equuivalent to your previous command (but not sure what is the point to compute them twice)
Olivier

Can you help with this problem?

Provide an answer of your own, or ask Jochen Jens Heinrich for more information if necessary.

To post a message you must log in.