Calculation of multiprocesses matrix element

Asked by Nikita

Good day,

Right now I'm trying to calculate the matrix element of VBF Higgs production and I don't quite understand how to combine all the channels together.

When I use the following commands
import model HC_UFO
generate p p > j j x0 $$ w+ w- z a QCD=0
output standalone_cpp VBFME
I got 77 different matrix elements in SubProcesses folder.
Is there some reasonable way to calculate the "combined" matrix element?
I wouldn't actually like to calculate them separately and then weight them by PDF's values myself...

Many thanks in advance,
Nikita

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi Nikita,

You need to define what you call "combined" matrix element before I can comment.
But i guess that in any case you will need to define it yourself.

Olivier

> On 5 Jun 2018, at 12:43, Nikita <email address hidden> wrote:
>
> New question #669983 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/669983
>
> Good day,
>
> Right now I'm trying to calculate the matrix element of VBF Higgs production and I don't quite understand how to combine all the channels together.
>
> When I use the following commands
> import model HC_UFO
> generate p p > j j x0 $$ w+ w- z a QCD=0
> output standalone_cpp VBFME
> I got 77 different matrix elements in SubProcesses folder.
> Is there some reasonable way to calculate the "combined" matrix element?
> I wouldn't actually like to calculate them separately and then weight them by PDF's values myself...
>
> Many thanks in advance,
> Nikita
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Nikita (zerousik) said :
#2

Hi Oliver,

Thanks for your response.
By "combined" matrix element I meant the matrix element that takes into account all the VBF Higgs production diagrams, i.e. the matrix element of the p p > j j x0 process, where p and j can be any allowed quarks or gluons.

Best Regards,
Nikita

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#3

Hi,

You mean just do the sum of them?

If you want to change the output to have a fully optimal way to compute such sum.
You will need to create a quite complex plugin to do that.

Cheers,

Olivier

> On 5 Jun 2018, at 17:42, Nikita <email address hidden> wrote:
>
> Question #669983 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/669983
>
> Status: Answered => Open
>
> Nikita is still having a problem:
> Hi Oliver,
>
> Thanks for your response.
> By "combined" matrix element I meant the matrix element that takes into account all the VBF Higgs production diagrams, i.e. the matrix element of the p p > j j x0 process, where p and j can be any allowed quarks or gluons.
>
> Best Regards,
> Nikita
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Nikita (zerousik) said :
#4

Hi Oliver,

> You mean just do the sum of them?
I'm afraid it's not simply the sum, because if we take into account different subprocesses we will face with the interference effects between different diagrams while calculating the amplitude, so the resulting matrix element will have additional terms.
So, if I understood correctly, there is no "out-of-the-box" way to do that?

Best Regards,
Nikita

Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said :
#5

Hi,

Which kind of interference are you talking about?
The different processes have different final state and therefore are distinguishable (at least in the S-matrix formalism)
So the interference term is zero. I obviously know some situation where this is not the case but this is typically the case when perturbation theory are not valid anymore.

> So, if I understood correctly, there is no "out-of-the-box" way to do that?

Yeah that's for sure (especially since I still do not know how you want to compute such interference)

Cheers,

Olivier

Revision history for this message
Nikita (zerousik) said :
#6

Thanks Olivier Mattelaer, that solved my question.