InvalidCmd : No amplitudes generated from process

Asked by Benedikt Maier

Dear experts,

in the 5FS scheme I am generating:

set group_subprocesses Auto
Set group_subprocesses to Auto
Note that you need to regenerate all processes
set ignore_six_quark_processes False
set loop_optimized_output True
set loop optimized output to True
set gauge unitary
Passing to gauge unitary.
set complex_mass_scheme False
import model loop_sm-no_b_mass
INFO: load particles
INFO: load vertices
INFO: Restrict model loop_sm with file models/loop_sm/restrict_no_b_mass.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Defined multiparticle p = g u c d s u~ c~ d~ s~
Defined multiparticle j = g u c d s u~ c~ d~ s~
Defined multiparticle l+ = e+ mu+
Defined multiparticle l- = e- mu-
Defined multiparticle vl = ve vm vt
Defined multiparticle vl~ = ve~ vm~ vt~
Defined multiparticle all = g gh gh~ d u s c b d~ u~ s~ c~ b~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t t~ z w+ h w- ta- ta+
Switching from interface MG5 to aMC@NLO
define l+ = e+ mu+ ta+
Defined multiparticle l+ = e+ mu+ ta+
define l- = e- mu- ta-
Defined multiparticle l- = e- mu- ta-
define j = j b b~
Defined multiparticle j = g u c d s b u~ c~ d~ s~ b~
define p = j
Defined multiparticle p = g u c d s b u~ c~ d~ s~ b~
define t1 = t t~
Defined multiparticle t1 = t t~
generate p p > t1 h j $$ w+ w- [QCD]

####################################################

All works fine then, I get a nice cross section, in the MadSpin section, however, it fails due to some reason I do not understand. Maybe an incorrect syntax in the decay? here I give the screen output:

####################################################

************************************************************
* *
* W E L C O M E to M A D S P I N *
* *
************************************************************
INFO: Extracting the banner ...
INFO: process: p p > t1 h j
INFO: options: $$ w+ w-
INFO: detected model: loop_sm-no_b_mass. Loading...
WARNING: The UFO model does not include partial widths information.
Impossible to use analytical formula, will use MG5/MadEvent (slower).
#************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph 5 Framework: *
#* The MadGraph Development Team - Please visit us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#Some options (optional)
#
# set seed 1
# set max_weight -1
#
#
# specify the decay for the final state particles
decay h > b b~
decay t~ > w- b~, w- > l- vl~
decay t > w+ b, w+ > l+ vl
# running the actual code
launch
INFO: Will use seed 859922249
INFO: We need to recalculate the branching fractions
INFO: compute_width module not available, use numerical estimates instead
INFO:
INFO: decay channels for t~ :
INFO: BR d1 d2
INFO: 1.000000e+00 b~ w-
INFO:
INFO:
INFO: decay channels for h :
INFO: BR d1 d2
INFO: 1.000000e+00 ta+ ta-
INFO:
INFO:
INFO: decay channels for w- :
INFO: BR d1 d2
INFO: 3.333610e-01 s c~
INFO: 3.333610e-01 d u~
INFO: 1.111195e-01 mu- vm~
INFO: 1.111195e-01 e- ve~
INFO: 1.110390e-01 ta- vt~
INFO:
INFO:
INFO: decay channels for w+ :
INFO: BR d1 d2
INFO: 3.333610e-01 s~ c
INFO: 3.333610e-01 d~ u
INFO: 1.111195e-01 mu+ vm
INFO: 1.111195e-01 e+ ve
INFO: 1.110390e-01 ta+ vt
INFO:
INFO:
INFO: decay channels for t :
INFO: BR d1 d2
INFO: 1.000000e+00 b w+
INFO:
INFO: generating the production square matrix element
INFO: generate p p > t1 h j $$ w+ w- ;add process p p > t1 h j j $$ w+ w- ;
INFO: Done 26.58
INFO: generating the full square matrix element (with decay)
INFO: generate p p > t1 h j $$ w+ w- , (t~ > b~ w- , w- > l- vl~ ), h > b b~ , (t > b w+ , w+ > l+ vl );add process p p > t1 h j j $$ w+ w- , (t~ > b~ w- , w- > l- vl~ ), h > b b~ , (t > b w+ , w+ > l+ vl ) ;
Command "generate p p > t1 h j $$ w+ w- , (t~ > b~ w- , w- > l- vl~ ), h > b b~ , (t > b w+ , w+ > l+ vl )" interrupted with error:
InvalidCmd : No amplitudes generated from process Process: h > b b~ WEIGHTED=4. Please enter a valid process
Command "add process p p > t1 h j j $$ w+ w- , (t~ > b~ w- , w- > l- vl~ ), h > b b~ , (t > b w+ , w+ > l+ vl ) " interrupted in sub-command:
"add process p p > t1 h j j $$ w+ w- , (t~ > b~ w- , w- > l- vl~ ), h > b b~ , (t > b w+ , w+ > l+ vl )" with error:
InvalidCmd : No amplitudes generated from process Process: h > b b~ WEIGHTED=4. Please enter a valid process
Command "launch auto " interrupted with error:
InvalidCmd : No processes generated. Please generate a process first.
quit

quit

#####################################################

Do you have an idea what is going wrong?

Thank you for all your help here, this forum helps us a lot!

ben

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi Ben,

Thanks for reporting.
In fact your are the second to report this Bug (See Bug #1176874).
In that topic, you will find a patch which correct this wrong MS behavior.

Cheers,

Olivier

Revision history for this message
Benedikt Maier (bmaier) said :
#2

Hi Olivier,

since I am also in contact with Alberto Orso who reported the bug you are linking here, I already knew about it. And I had already patched it. Before I had patched it, I also got an error with "Duplicate process..." like him.

But the error I posted above is indeed after the (successful) patch "patch -p1 < bug1176874.patch" ... :/

Do you have another idea?

Thanks in advance

Ben

Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said :
#3

Hi Ben,

Let me first explain you what is the current problem:

With your model (loop_sm-no_b_mass) the b_mass is set to zero.
Therefore it can't couple to Higgs and asking the h > b b~ simply crash
because the code didn't find any valid diagrams.

This being said, what you ask is legitimate, so I would like to have a way to
automatically runs the change of model required in that case. But this needs to be discussed
with the other author of MS. (I'll start that discussion right away).

In the mean time, the only solution that I see is to modify the model/parameter by hand.
1) run in no shower mode
2) edit the lhe events and change the model information from loop_sm-no_b_mass loop_sm)
and put the b_mass to a non zero value. (you need to check that your param_card is indeed a valid one for the new model)
3) run MS+ the shower via the script ./bin/aMCatNLO
(then enter the command
decay_events
shower

Cheers,

Olivier

On May 10, 2013, at 6:36 AM, Benedikt Maier <email address hidden> wrote:

> Question #228419 on MadGraph5 changed:
> https://answers.launchpad.net/madgraph5/+question/228419
>
> Status: Answered => Open
>
> Benedikt Maier is still having a problem:
> Hi Olivier,
>
>
> since I am also in contact with Alberto Orso who reported the bug you are linking here, I already knew about it. And I had already patched it. Before I had patched it, I also got an error with "Duplicate process..." like him.
>
> But the error I posted above is indeed after the (successful) patch
> "patch -p1 < bug1176874.patch" ... :/
>
> Do you have another idea?
>
> Thanks in advance
>
> Ben
>
> --
> You received this question notification because you are a member of
> MadTeam, which is an answer contact for MadGraph5.

Revision history for this message
Benedikt Maier (bmaier) said :
#4

Thanks Olivier Mattelaer, that solved my question.