Abstract

We introduce a new method for the application of one-loop integrand reduction via the Laurent expansion algorithm, as implemented in the public C++ library Ninja. We show how the coefficients of the Laurent expansion can be computed by suitable contractions of the loop numerator tensor with cut-dependent projectors, making it possible to interface Ninja to any one-loop matrix element generator that can provide the components of this tensor. We implemented this technique in the Ninja library and interfaced it to MadLoop, which is part of the public MadGraph5_aMC@NLO framework. We performed a detailed performance study, comparing against other public reduction tools, namely CutTools, Samurai, IREGI, PJFry++ and Golem. We find that Ninja outperforms traditional integrand reduction in both speed and numerical stability, the latter being on par with that of the tensor integral reduction tool Golem which is however more limited and slower than Ninja. We considered many benchmark multi-scale processes of increasing complexity, involving QCD and electro-weak corrections as well as effective non-renormalizable couplings, showing that Ninja's performance scales well with both the rank and multiplicity of the considered process.

Highlights

  • Perturbation theory as a Taylor expansion in the coupling constants

  • We introduce a new method for the application of one-loop integrand reduction via the Laurent expansion algorithm, as implemented in the public C++ library Ninja

  • We show how the coefficients of the Laurent expansion can be computed by suitable contractions of the loop numerator tensor with cut-dependent projectors, making it possible to interface Ninja to any one-loop matrix element generator that can provide the components of this tensor

Read more

Summary

Tensor integrands

A generic one-loop amplitude can be written as a sum of n-point integrals of the form. The numerator N of the integrand is a polynomial in the components of the d-dimensional loop momentum q, with d = 4 − 2 , while the denominators Di correspond to Feynman loop propagators and they have the general quadratic form. [19], addressing the case where a numerical evaluation of the integrand is available but its full polynomial structure is not known, reconstructs the entries of the tensor by sampling the numerator on several values of the loop momentum. We will consider a one-loop tensor integrand to be defined by the entries of the symmetric tensor numerator in eq (2.6), as well as the momenta pi and masses mi appearing in the loop denominators as in eq (2.2). Thanks to the new projection techniques introduced in this paper, these are the only input required by Ninja for performing the corresponding loop reduction

Semi-numerical integrand reduction via Laurent expansion
Notation
Recursive formulas
Tensor projectors for Laurent-expansion terms
Numerator evaluation
The μ2-expansion
The t3-expansion
The t2-expansion
Implementation
Applications
Timing profile
Stability study
Conclusions
A Explicit formulas for expansion terms
B Details of timing performances
C Details of stability performances
D Numerical results for benchmark processes
Findings
SF colour
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call