MadGraph5_aMC@NLO

How do I generate pp → H+ → h W+ → γγ l+ νl

Asked by Amine Alida

Hello, I am new to madgraph (v 3.5.4) and I want to generate the process pp → H+ → h W+ → γγ l+ νl using the model 2HDMtII_NLO.
When I try to do it like this :

MG5_aMC>import model 2HDMtII_NLO/
MG5_aMC>generate p p > h+ > w+ h1, (w+ > l+ vl, h1 > a a)

it gives me this error :

INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=4: WEIGTHED IS QCD+2*QED
INFO: Trying coupling order WEIGHTED<=5: WEIGTHED IS QCD+2*QED
INFO: Trying coupling order WEIGHTED<=6: WEIGTHED IS QCD+2*QED
INFO: Trying coupling order WEIGHTED<=7: WEIGTHED IS QCD+2*QED
INFO: Trying process: u d~ > h+ > w+ h1 WEIGHTED<=8 @1
INFO: Trying process: u s~ > h+ > w+ h1 WEIGHTED<=8 @1
INFO: Trying process: c d~ > h+ > w+ h1 WEIGHTED<=8 @1
INFO: Trying process: c s~ > h+ > w+ h1 WEIGHTED<=8 @1
INFO: Trying process: d~ u > h+ > w+ h1 WEIGHTED<=8 @1
INFO: Trying process: s~ c > h+ > w+ h1 WEIGHTED<=8 @1
Command "generate p p > h+ > w+ h1, (w+ > l+ vl, h1 > a a)" interrupted with error:
NoDiagramException : No amplitudes generated from process Process: g/u/c/d/s/u~/c~/d~/s~ g/u/c/d/s/u~/c~/d~/s~ > h+> w+ h1 WEIGHTED=8 @1
   Decay: w+ > e+/mu+ ve/vm/vt
     Decay: h1 > a a. Please enter a valid process

Regards.
Amine.

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Sihyun Jeon
Solved:
Last query:
Last reply:
Revision history for this message
Best Sihyun Jeon (shjeon) said : 		#1

You are trying to produce h+ from initial state quarks which is not possible in the model definition. "pp > h+", such feynman diagram does not exist but you are trying to enforce h+ to be present in the diagram which results in an error.

Revision history for this message
Amine Alida (amineaal) said (last edit ): 		#2

That solved my question. I get it now.
I tried p p > h+ t~ b, (h+ > w+ h1, w+ > l+ vl, h1 > a a) instead of p p > h+ > w+ h1, (w+ > l+ vl, h1 > a a) and now it works.
Thank you Sihyun.