# how many diagrams ?

Dear mg5 team

the process that I generate is "generate g g > e+ ve mu- vm~ b b~ b b~ [virt=QCD]"

then I get "INFO: Contributing diagrams generated: 3904 Born, 172416(+53912) loops, 90804 R2, 142528 UV"

so, how many 1-loop diagrams do I have??

Best regards

Huanyu

## Question information

- Language:
- English Edit question

- Status:
- Solved

- Assignee:
- Valentin Hirschi Edit question

- Solved by:
- Olivier Mattelaer

- Solved:

- Last query:

- Last reply:

HuanyuBi (bihuanyu) said : | #1 |

hi Olivier

The question is asked in English. I didnt notice that I choose Haida. but I modified it.

Best regards

Huanyu

在2021-04-18 06:14:02,Olivier Mattelaer<email address hidden>写道：

> Your question #696621 on MadGraph5_aMC@NLO changed:

> https:/

>

> Assignee: None => Valentin Hirschi

>

> WARNING: This question is asked in a language (Haida)

> spoken by none of the registered MadGraph5_aMC@NLO answer contacts.

>

> --

> You received this question notification because you asked the question.

The number of 1-loop diagram is given by 172416(+53912)

the 53k information is probably loop which are identical to another loop and therefore does not need to be computed (typically different mass-less flavor in a loop). So the number of loop that we evaluated is 172k. Which are then interfere with the 4k born amplitude.

Then in order to compute the finite piece of the loop you have to add 90k R2 term.

and then add the 142k UV counter-term.

PS: Valentin correct me if I'm wrong

HuanyuBi (bihuanyu) said : | #3 |

hi Olivier

thanks for your answer, i have further questions.

is it possoble for mg5 to generate (or export/print out) diagrams of specified loop topology.

for example. i need only self-energy diagrams or Octagon-type diagrams.

Best regards

Huanyu

You need to hardcode some filtering:

See:

https:/

HuanyuBi (bihuanyu) said : | #5 |

thanks for your answer.

I find this example in the file, madgraph/

# Ex. 1: Chose the topology, i.e. number of loop line.

# Notice that here particles and antiparticles are not

# differentiated and always the particle PDG is returned.

# In this example, only boxes are selected.

# if len(diag.

# any([i in diag.get_

# valid_diag=False

so if I only need boxes then I uncomment the last 3 lines.

but I am not sure what should I do if I only need Self-energy, Vertex, Pentagon-type, Hexagon-type, Heptagon-type and Octagon-type diagrams.

best regards

Huanyu

I guess that the second line is not needed for what you want to do and this does not "only" select boxes for sure.

If I were you, I would test something like

if len(diag.

if len(diag.

if len(diag.

if len(diag.

if len(diag.

if len(diag.

if len(diag.

and obviously check the drawing of the diagram to check if this correspond to what you want or not

HuanyuBi (bihuanyu) said : | #7 |

Hi Olivier,

Thanks for your answer, and it solves my question.

just one comment:

actually if I only need Octagon-type diagrams, I need

if len(diag.

valid_diag=False

Best regards

Huanyu

HuanyuBi (bihuanyu) said : | #8 |

Thanks Olivier Mattelaer, that solved my question.