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

Asked by Jay Sandesara on 2020-07-02


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

English Edit question
MadGraph5_aMC@NLO Edit question
No assignee Edit question
Solved by:
Olivier Mattelaer
Last query:
Last reply:
Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said :

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



Revision history for this message
Jay Sandesara (jaysandesara) said :

Thanks Olivier Mattelaer, that solved my question.