avoid unwanted feynman diagrams in SMEFTsim

Asked by Neil Warrack on 2018-10-26

Hello experts,
I am new to MadGraph and I'm using the SMEFTsim model (http://feynrules.irmp.ucl.ac.be/wiki/SMEFT using this download: SMEFTsim_A_general_alphaScheme_UFO_v2.tar.gz) I am trying to generate p p > t t~ j j with one or two wilson coeffs set to a non-zero value. This produces lots and lots of unwanted Feynman diagrams that will not contribute to the event generation. Is there a way to avoid this?

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
2018-11-11
Last query:
2018-11-11
Last reply:
2018-11-11

Would that link help you?

Cheers,

Olivier
FAQ #2312: “FR Model much slower than build-in MG model. Why and how to fix?”.

Neil Warrack (nwarrack) said : #2

Hi Oliver,

Thanks for the reply.

I'm not sure that what you suggest answers my question. However, I don't think I fully understood my question when I asked it (still figuring all this MadGraph stuff out - please forgive me).

I'm interested in looking at the process pp->tt~ through the eyes of a SMEFT. If I tell madgraph that I want a dileptonic decay by specifying the output like so:
pp -> b b~ l+ vl u~ d j j
then I suppose it will generate many diagrams that do not contain top quarks. Can I rephrase my question like this: What is the best way to avoid generating the diagrams that I don't want (ie diagrams without tops). Is it still by means of generating my own restrictions card?

Thanks again for your help. Sorry if this is a trivial question.

Neil Warrack (nwarrack) said : #3

Adding a comment:
If this constitutes a new question then let me know and I can post it as one. Thanks.

Hi,

In this case, I guess that what you are looking for is the decay-chain syntax:

generate p p > t t~, (t > w+ b, w+ > l+ vl), ( t~ > w- b~ , w- > l- vl~)

Note that this syntax is associate to a breit-wigner cut to ensure that the top (and W) are "close" to their pole mass.
(You can change the cut in the run_card, but by default "close" means up to 15 times the width of the particle)

Cheers,

Olivier

Neil Warrack (nwarrack) said : #5

Thanks Olivier Mattelaer, that solved my question.

Neil Warrack (nwarrack) said : #6

We were also worried about very off-shell tops too. But you've answered that concern too. Cheers!