Generating loop induced p p > 4lep + 1 jet process

Asked by Jay Sandesara on 2020-07-02

Hello,

Upon your previous recommendation I have started generating p p > l+ l- l+ l- j [noborn=QCD] processes to include all the flavours in initial state. In order to do this I use the following filters:

 if any([abs(pdg) not in range(1,7) for pdg in diag.get_loop_lines_pdgs()]):
   valid_diag = False
 if (23 not in diag.get_pdgs_attached_to_loop(structs)) and (25 not in diag.get_pdgs_attached_to_loop(structs)) and (22 not in diag.get_pdgs_attached_to_loop(structs)):
   valid_diag = False

The idea behind these filters is that they will give me gg-zz, gg-za and gg-aa diagrams as well as gg-h. Doing this gives me about 8k Feynman diagrams, and a few of the following warnings:

WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation. Make sure you did not want to include them.

Would it be safe to ignore these?

Question information

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

If your loop computation is finite (did you check that with the standalone mode) then it should be ok.

Cheers,

Olivier

Jay Sandesara (jaysandesara) said : #2

Thanks Olivier Mattelaer, that solved my question.