InvalidCmd : No particle b~) in model

Asked by muhammad gul

Dear,

I am using 2HDM+a model:
https://github.com/LHC-DMWG/model-repository/tree/DMSimp_t/models/Pseudoscalar_2HDM
with madgraph 5 v3 using the following commands

python3.8 ./bin/mg5_aMC
import model Pseudoscalar_2HDM --modelname

(it's output is:
MG5_aMC>import model Pseudoscalar_2HDM --modelname
INFO: Restrict model Pseudoscalar_2HDM with file models/Pseudoscalar_2HDM/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
Kept definitions of multiparticles p / j / l+ / l- / vl / vl~ unchanged
Defined multiparticle all = g ghg ghg~ u c d s u~ c~ d~ s~ a gha gha~ ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ ghz ghwp ghwm h+ h1 h2 h3 h4 w- ghz~ ghwp~ ghwm~ h- ta- xd ta+ xd~

 )

>> Following command works fine.
g g > xd xd~ h1 [noborn=QCD]

>> But when i am going to decay h1 like:

generate g g > xd xd~ h1 [noborn=QCD], (h1 > b b~)
Command "generate g g > xd xd~ h1 [noborn=QCD], (h1 > b b~)" interrupted with error:
InvalidCmd : No particle b~) in model

What is going on?
Cheers,
M

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
Best Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

The decay chain syntax (the use of the comma) is not supported for loop-induced processes (same for NLO processes).
So the decay of onshell particle needs to be done via MadSpin (or other method).

In the case of loop-induced, the default MadSpin algorithm is not working (it is for NLO processes). So you need to reduce the accuracy by using the mode "spinmode=none". Note that in the case of a scalar with a very small width decaying (like it is for the SM Higgs), this is not a real physics limitation.

Cheers,

Olivier

Revision history for this message
muhammad gul (mgul) said :
#2

Thanks Olivier Mattelaer, that solved my question.