Problem generating from model

Asked by Giacomo

Hello everyone;

I'm quite new to madgraph and event generation. I was asked to generate the following process from the model RS_bulk_ktilda:
proton proton collision creating a gravition (y, pdg_id = 39) which decays in two Higgs each of them decaying w+ w-. Now 3 ws decay hadronically and one creates a lepton neutrino pair. This process as written gives errors but if i impose one H to decay into a b b~ pair or a photon pair everything is fine. I'll post here the proc card that i use to generate the 1st event that doesn't work with its error. Am i missing something here or am i doing something wrong? Thanks everyone.

import model RS_bulk_ktilda
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define q = u c d s b u~ c~ d~ s~ b~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define lep = e+ e- mu+ mu-
define nus = ve~ vm~ ve vm
define ws = w+ w-
define vs = ws z
define vis = lep q
generate p p > y, (y > H H, (H > w+ w-, w+ > q q, w- > lep nus), (H > w+ w-, w+ > q q, w- > q q ) ) @1
add process p p > y, (y > H H, (H > w+ w-, w- > q q, w+ > lep nus), (H > w+ w-, w+ > q q, w- > q q ) ) @2
output Bulk_signal
launch

The error here is:
.
.
.
Survey return zero cross section.
   Typical reasons are the following:
   1) A massive s-channel particle has a width set to zero.
   2) The pdf are zero for at least one of the initial state particles
      or you are using maxjetflavor=4 for initial state b:s.
   3) The cuts are too strong.
   Please check/correct your param_card and/or your run_card.
Zero result detected: See https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/FAQ-General-14

note that the following processes work fine:

generate p p > y, (y > H H, H > q q, (H > w+ w-, w+ > q q, w- > lep nus) ) @1
add process p p > y, (y > H H, H > q q, (H > w+ w-, w- > q q, w+ > lep nus) ) @2

generate p p > y, (y > H H, H > a a, (H > w+ w-, w+ > q q, w- > lep nus) ) @1
add process p p > y, (y > H H, H > a a, (H > w+ w-, w- > q q, w+ > lep nus) ) @2

And this does not work either:

generate p p > y, (y > H H, H > z z, (H > w+ w-, w+ > q q, w- > lep nus) ) @1
add process p p > y, (y > H H, H > z z, (H > w+ w-, w- > q q, w+ > lep nus) ) @2

generate p p > y, (y > H H, H > w+ w-, (H > w+ w-, w+ > q q, w- > lep nus) ) @1
add process p p > y, (y > H H, H > w+ w-, (H > w+ w-, w- > q q, w+ > lep nus) ) @2

This doesn't work too but the error this time is different:

generate p p > y, (y > H H, H > w+ w-, H > w+ w-) @1
.
.
.
Zero result detected: No Phase Space. Please check particle masses.

INFO:
quit
INFO:
Generation failed (no results.dat file found)

Let me know if you need more information about the model or so on, i'll provide them. Thanks again for your kind attention. have a nice day.

Giacomo

Question information

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

Hi,

You are using the "onshell" decay syntax.
So you can not ask for
H > w+ w-
since the mass of the H is smaller than twice the mass of the W.

Cheers,

Olivier

Revision history for this message
Giacomo (giacbold) said :
#2

Hi, thank you for your answer

Yes indeed one of the w is off-shell. I thought it was automatic sorry. How should i tell mg5 that one w is off shell? Tried putting a $ in front of one but it does not recognize this command (or it seems like so). Also searched for a variable bw_cut (as written in pag 30 https://cp3.irmp.ucl.ac.be/projects/madgraph/raw-attachment/wiki/Madrid2015_2%20MC/Tutorial.pdf )but couldn't find it. How would you suggest to move on? Thanks in advance.

Giacomo

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

you can either ask directly for
H > q q lep nus

no?

Olivier

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

If you are sure that the following syntax makes sense, (that it is gauge invariant)
you can also ask for
H > w+ w- > q q lep nus
This syntax is quite dangerous (lot of side effects and can lead to unphysical result when a breaking of gauge invariance occurs)
You will find more details about it in the Tutorial that you mention.

Cheers,

Olivier

Revision history for this message
Giacomo (giacbold) said :
#5

This seems right. Thank you very much for your help.

Can you help with this problem?

Provide an answer of your own, or ask Giacomo for more information if necessary.

To post a message you must log in.