# Return matrix element squared of all terms

Dear experts,

I am generating the following process: g g > h1 h1 [QCD] using the 2HDM_NLO model from the feynrules database. It consists of three types of diagrams, let's call them M1, M2 and M3. These diagrams include one loop of fermions. We want to have separately the contribution for each diagram when only top or bottom quarks are allowed to run in this loop. Therefore, each diagram can be seen as a sum of two: M1= M1t+ M1b, M2=M2t+ M2b,M3= M3t+ M3b, where t and b indicate the fermion contributing to the loop. The total matrix element is then given by:

|M|² =|M1t+M1b+

My question is the following: is there a way of generating the entire process at once and then access the value of the matrix element squared for each of the terms described above?

Thank you,

Ana Luisa

## Question information

- Language:
- English Edit question

- Status:
- Open

- Assignee:
- Valentin Hirschi Edit question

- Last query:

- Last reply:

## Can you help with this problem?

Provide an answer of your own, or ask Ana Luisa Moreira de Carvalho for more information if necessary.