start

NNLL-fast

NNLL-fast is a computer program which computes the squark and gluino hadroproduction cross sections including approximate next-to-next-to-leading order (NNLOApprox) supersymmetric QCD corrections and the resummation of soft gluon emission at next-to-next-to-leading-logarithmic (NNLL) accuracy, including the resummation of Coulomb gluons with a NLO Coulomb potential and bound-state contributions below the production threshold. NNLL-fast also provides an estimate of the theoretical uncertainty due to scale variation, the parton distribution functions, and the value of the strong coupling αs. The program reads in grid files with NLO and NNLOApprox+NNLL results, the scale uncertainty and the PDF and αs error and uses a fast interpolation routine to provide predictions for any value of squark and gluino masses in the mass ranges specified below. The results are provided for the latest LHC collision energy of 13.6 TeV for Run 3, and will be updated for future runs of the LHC. The predictions are created with the recent PDF4LHC21_40_pdfas PDF sets.

NNLL-fast is the successor to NLL-fast, the code which provided predictions at NLO+NLL accuracy for squark and gluino production with previous generations of PDF sets.

Separate versions of NNLL-fast for other models are also available:

  • NNLL-fast-MRSSM for the calculation of squark-antisquark and squark-squark production in the Minimal R-symmetric Supersymmetric Standard Model (MRSSM) at NLO+(N)NLL accuracy,
  • NNLL-fast-LQ for the calculation of scalar leptoquark pair-production at NLO+NNLL accuracy.

NNLL-fast is based on results obtained by many authors over the years. In alphabetical order, the contributions were provided by: Wim Beenakker, Silja Brensing-Thewes, Christoph Borschensky, Michael Krämer, Anna Kulesza, Eric Laenen, Laura Moreno Valero, Leszek Motyka, Irene Niessen, and Daniel Schwartländer. Additionally, NNLL-fast uses NLO results obtained with Prospino. Therefore we ask to cite the following papers when referring to NNLL-fast:

Squark and gluino production:

Stop (sbottom) production:

Squark production in the MRSSM:

NLO numbers for squark production in the MRSSM are obtained on the basis of:

Scalar Leptoquark pair production:

NLO cross sections with t-channel are calculated with the help of MadGraph5_aMC@NLO, please cite

NEW: NNLL-fast-MRSSM (LHC @ 13 TeV and LHC @ 13.6 TeV)

  • Main program and grids in one package: for the LHC @ 13 TeV: nnllfast-mrssm-13; for the LHC @ 13.6 TeV: nnllfast-mrssm-13.6.
  • The grids are generated with the PDF4LHC21 Hessian PDF set (LHAPDF ID 93300).
  • Available for the processes of squark-antisquark production at NLO+NNLL accuracy, and squark-squark (+ antisquark-antisquark) production at NLO+NLL.
  • The theoretical scale uncertainty is determined from the variation of the total cross section when simultaneously varying the renormalisation and factorisation scales around their respective central values up and down by a factor of 2.
  • Please note that the the PDF and αs uncertainties are given as the combined uncertainty from the quadratic sum of the two.
  • NNLL-fast-MRSSM is run in the same way as NNLL-fast, as described below, for the two implemented processes of squark-antisquark and squark-squark production, with the squark and gluino masses as an input.
  • Additionally to the scale and PDF uncertainties, NNLL-fast-MRSSM also outputs an uncertainty with respect to the variation of the sgluon mass from 100 TeV (central value) to 2 TeV.

NEW: NNLL-fast, version 2.0 (LHC @ 13.6 TeV)

  • Main program and grids in one package nnllfast-2.0.
  • The grids are generated for the new LHC Run 3 @ 13.6 TeV with the PDF4LHC21 Hessian PDF set (LHAPDF ID 93300).
  • The theoretical scale uncertainty is determined from the variation of the total cross section when simultaneously varying the renormalisation, factorisation, Coulomb, and Bohr scales around their respective central values up and down by a factor of 2.
  • Please note that the the PDF and αs uncertainties are given as the combined uncertainty from the quadratic sum of the two.
  • For grids for stop production, the light flavour squarks and the heavy stop-2 are considered as very heavy, while the stop mixing angle corresponds to CMSSM benchmark point 40.2.5 as defined in arXiv:1109.3859. Additionally, we provide a parametrisation error for the total cross section when varying the stop mixing angle in the range $-1 \le \sin 2\theta_{\tilde{t}} \le +1$.

NNLL-fast-LQ (LHC @ 13 TeV)

  • Main program and grids in one package for calculation of NLO(QCD)+NNLL nnllfast-LQ.
  • The output can also be used to calculate full NLO (including t-channel contributions) + NNLL.

NNLL-fast, version 1.1 (LHC @ 13 TeV)

  • Main program and grids in one package nnllfast-1.1.
  • The package consists of NNLL-fast version 1.0 package, extended by results for squark production in the large gluino mass limit and gluino production in the large squark mass limit (decoupled cases).

NNLL-fast, version 1.0 (LHC @ 13 TeV)

  • Main program and grids in one package nnllfast-1.0.
  • Upon request, the process of gluino pair-production is also available for gluino masses ranging from 100 to 200 GeV and squark masses ranging from 500 to 3000 GeV.
  • The theoretical scale uncertainty is determined from the variation of the total cross section when simultaneously varying the renormalisation, factorisation, Coulomb, and Bohr scales around their respective central values up and down by a factor of 2.
  • Please note that the the PDF and αs uncertainties are given as the combined uncertainty from the quadratic sum of the two.
  • For grids for stop production, the light flavour squarks and the heavy stop-2 are considered as very heavy, while the stop mixing angle corresponds to CMSSM benchmark point 40.2.5 as defined in arXiv:1109.3859. Additionally, we provide a parametrisation error for the total cross section when varying the stop mixing angle in the range $-1 \le \sin 2\theta_{\tilde{t}} \le +1$.

NNLL-fast-LQ is a python script that runs with python 2.7.X or python 3 with numpy and scipy.

NNLL-fast is written in the Fortran programming language, so most Fortran compilers should be able to compile the code. One that it has been tested with is the free GNU Fortran compiler gfortran. On a Linux terminal, the command to compile the code is:

$ gfortran nnllfast.f -o name_of_the_executable

where nnllfast.f has to be replaced by the correct Fortran source file, depending on the version of NNLL-fast.

NNLL-fast-LQ: can be used in two different ways. NNLLfast_LQ.py can be executed with command line arguments:

$ python NNLLfast_LQ.py PDF m

Here PDF is the parton distribution function and can be either NNPDF31 or CT18. m is the leptoquark mass in GeV between 500 GeV and 2300 GeV. The output is in this case written on screen. NNLLfast_LQ.py can also be executed without command line arguments:

$ python NNLLfast_LQ.py

In this case the script uses PDF and masses, which are defined in the script NNLLfast_LQ.py. The default values can be edited to obtain the output for a list of masses between 500 GeV and 2300 GeV. The output is written in a file with default name output.out.

NNLL-fast: After compiling nnllfast.f, it should be called in the following way for all four squark and gluino production processes:

$ ./name_of_the_executable process squark_mass gluino_mass

or, for the case of stop-antistop (sbottom-antisbottom) production:

$ ./name_of_the_executable st stop_mass gluino_mass

The command line attribute process can take the following values:

sb     squark-antisquark production,
ss     squark-squark production,
gg     gluino-gluino production,
sg     squark-gluino production.

Ranges of gluino and squark masses (squark_mass, gluino_mass, stop_mass) are:

for gluino-pair production:
   500 GeV  <= gluino mass <= 3000 GeV    and   500 GeV <= squark mass <= 3000 GeV,
for squark-gluino  production:
   500 GeV  <= gluino mass <= 3000 GeV    and   500 GeV <= squark mass <= 3000 GeV,
for squark-antisquark and squark-pair production:
   500 GeV  <= squark mass <= 3000 GeV    and   500 GeV <= gluino mass <= 3000 GeV
and for stop-antistop production:
   100 GeV  <=  stop mass  <= 3000 GeV    and   500 GeV <= gluino mass <= 5000 GeV 
   (other SUSY parameters: light flavour squark and stop 2 masses considered as very heavy,
                           and stop mixing angle according to CMSSM benchmark 
                           point 40.2.5 as defined in arXiv:1109.3859.)  

NNLL-fast:

# LHC @ 13.6 TeV, NNLO PDF4LHC21 (LHAPDF ID 93300)
# process: sg   
# ms[GeV]  mg[GeV]   NLO[pb]     NNLL+NNLO_app[pb]   d_mu+[%]    d_mu-[%]   d_pdfas+[%]  d_pdfas-[%]  K_NNLL
--------------------------------------------------------------------------------------------------------------
  1700.    2100.    0.733E-02       0.891E-02          3.69       -5.49       9.86        -9.86        1.21

NNLL-fast-MRSSM:

# LHC @ 13.6 TeV, NNLO PDF4LHC21 (LHAPDF ID 93300)
# process: sb    (in the MRSSM with msgluon = 100 TeV)
# ms[GeV]  mg[GeV]   NLO[pb]     NNLL+NNLO_app[pb]  d_msgluon[%]  d_mu+[%]    d_mu-[%]   d_pdfas+[%]  d_pdfas-[%]  K_NNLL
----------------------------------------------------------------------------------------------------------------------------
  1500.    1800.    0.467E-02       0.508E-02         2.93          3.70       -3.49      16.03       -16.03        1.09

NNLL-fast-LQ:

NNLLfast output for leptoquark pair production
LHC @ 13 TeV, NNPDF31
mLQ = 1000.0 GeV
NLO(QCD) = 5.237E-03 + 6.049E-04 - 7.102E-04 +- 1.990E-04 pb
NLO(QCD)+NNLL = 5.436E-03 + 2.302E-05 - 2.302E-05 +- 2.011E-04 pb
delta NNLL(mu=m) = 5.221E-04 pb
delta NNLL(mu=m/2) = 1.912E-05 pb
delta NNLL(mu=2m) = 1.201E-03 pb

The previous versions of NLL-fast can be downloaded from here.

If you have questions, comments, or want to report a problem regarding the code, please contact Anna Kulesza or Christoph Borschensky.

  • start.txt
  • Last modified: 2024/02/16 08:48
  • by Christoph Borschensky