NLO+NLL SUSY-QCD (coord. Anna Kulesza, Christoph Borschensky)
References
- Squark and Gluino Production at Hadron Colliders, W. Beenakker, R. Höpker, M. Spira, P.M. Zerwas, Nucl. Phys. B492 (1997) 51-103
- Stop Production at Hadron Colliders, W. Beenakker, M. Krämer, T. Plehn, M. Spira, P.M. Zerwas, Nucl.Phys. B515 (1998) 3-14
- Threshold resummation for squark-antisquark and gluino-pair production at the LHC, A. Kulesza, L. Motyka, Phys. Rev. Lett. 102 (2009) 111802
- Soft gluon resummation for the production of gluino-gluino and squark-antisquark pairs at the LHC, A. Kulesza, L. Motyka, Phys. Rev. D80 (2009) 095004
- Soft-gluon resummation for squark and gluino hadroproduction, Wim Beenakker, Silja Brensing, Michael Krämer, Anna Kulesza, Eric Laenen, Irene Niessen, JHEP 0912 (2009) 041
- Supersymmetric top and bottom squark production at hadron colliders, Wim Beenakker, Silja Brensing, Michael Krämer, Anna Kulesza, Eric Laenen, Irene Niessen, JHEP 1008(2010)098
- Squark and gluino hadroproduction, W. Beenakker, S. Brensing, M. Krämer, A. Kulesza, E. Laenen, L. Motyka, I. Niessen, Int. J. Mod. Phys. A26 (2011) 2637-2664
Code
Prospino
Prospino2 is a computer program which computes next-to-leading order cross sections for the production of supersymmetric particles at hadron colliders. The processes currently included are squark, gluino, stop, neutralino/chargino, and slepton pair production. The code is available from the
Prospino web site.
NLL-fast
NLL-fast
is a computer program which computes the squark and gluino hadroproduction cross sections including next-to-leading order supersymmetric QCD corrections and the resummation of soft gluon emission at next-to-leading-logarithmic accuracy.
NLL-fast
also provides an estimate of the theoretical uncertainty due to scale variation and the parton distribution functions. The program reads in grid files with NLO and NLO+NLL results, the scale uncertainty and the pdf and alphas 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 for the LHC collision energy of 13 TeV can be obtained with the latest 3.0 version, whereas earlier 2.1 (1.2) versions deliver cross sections for 8 (7) TeV. Results are obtained with the MSTW2008NLO (MSTW2008LO for the LO results) and the CTEQ6.6M (CTEQ6L1) pdfs.
c
c After compiling the code, it should be called in the following way for all four squark and gluino
c production processes:
c
c ./name_of_the_executable process pdf squark_mass gluino_mass
c
c or, for the case of stop-antistop (sbottom-antisbottom) production:
c
c ./name_of_the_executable st pdf stop_mass
c
c or, for the case of gluino production in the limit of large squark masses (decoupling limit):
c
c ./name_of_the_executable gdcpl pdf gluino_mass
c
c or, for the case of squark production in the limit of large gluino masses (decoupling limit):
c
c ./name_of_the_executable sdcpl pdf squark_mass
c
c The command line attributes above can take the following values:
c - process:
c sb squark-antisquark production,
c ss squark-squark production,
c gg gluino-gluino production,
c sg squark-gluino production,
c - pdf:
c cteq CTEQ6.6M NLO (CTEQ6L1 for LO predictions) parton distribution functions,
c mstw MSTW2008NLO (MSTW2008LO for LO predictions) parton distribution functions,
c - masses for which the NLL+NLO cross section can be calculated:
c for gluino-pair production:
c 200 GeV <= gluino mass <= 3000 GeV and 200 GeV <= squark mass <= 5000 GeV,
c for gluino production in the decoupling limit (large squark mass)
c 200 GeV <= gluino_mass <= 3000 GeV
c for squark-gluino production:
c 200 GeV <= gluino mass <= 3000 GeV and 200 GeV <= squark mass <= 5000 GeV,
c for squark-antisquark and squark-pair production:
c 200 GeV <= squark mass <= 3000 GeV and 200 GeV <= gluino mass <= 3000 GeV
c for squark production in the decoupling limit (large gluino mass)
c 200 GeV <= squark mass <= 3000 GeV
c and for stop-antistop production:
c 100 GeV <= stop mass <= 2000 GeV and stop_mass <= gluino mass
c (other SUSY parameters according to CMSSM benchmark
c point 40.2.5 as defined in arXiv:1109.3859.)
c
c NOTE: option tot (active in versions 1.x) is disabled in this version.
c
c After reading in the mass values, the code returns (in order of appearance) the values of:
c - squark_mass, gluino_mass,
c - LO, NLO and NLL+NLO cross sections in pb,
c - upper and lower error on NLL+NLO due to scale mu variation between m/2 < mu < 2m,
c where m is the average mass of the pair-produced particles, also in pb
c - upper and lower 68% C.L. pdf error in percent (of the NLO value),
c - upper and lower error due to alpha_s uncertainty, in percent (of the NLO value),
c - NLO and NLL K-factors, with the NLL K-factor defined as sigma_NLL+NLO/sigma_NLO.
NEW: NLL-fast for LHC @ 14 TeV and LHC @ 33 TeV
- These versions of the package cover only the cases of stop-pair production and squark and gluino production in the decoupling limits. The main programs and grids can be downloaded from here. For grids for stop production SUSY parameters other that stop masses correspond to CMSSM benchmark point 40.2.5 as defined in arXiv:1109.3859.
NLL-fast, version 3.0 (LHC @ 13 TeV)
- main program and grids in one package nllfast-3.0. For grids for stop/sbottom production SUSY parameters other that stop/sbottom masses correspond to CMSSM benchmark point 40.2.5 as defined in arXiv:1109.3859.
NLL-fast, version 2.1 (LHC @ 8 TeV): files to download
- main program: nllfast-2.1.f.
- grids for inclusive squark and gluino production, stop (sbottom) production and production of squark and gluinos in the decoupling limits: nllfast-2.11grids.tar . For grids for stop/sbottom production SUSY parameters other that stop/sbottom masses correspond to CMSSM benchmark point 40.2.5 as defined in arXiv:1109.3859.
Older versions of NLL-fast
Sample output
# LHC @ 8 TeV, MSTW2008NLO (MSTW2008LO for LO)
# process: gg
# ms[GeV] mg[GeV] LO[pb] NLO[pb] NLL+NLO[pb] d_mu+[pb] d_mu-[pb] d_pdf+[%] d_pdf-[%] d_as+[%] d_as-[%] K_NLO K_NLL
------------------------------------------------------------------------------------------------------------------------------------
1000. 1000. 0.647E-02 0.101E-01 0.119E-01 0.646E-03 -0.828E-03 15.5 -15.4 5.3 -5.1 1.56 1.18
EWK corrections (coord. Wolfgang Hollik, Jonas Lindert)
First results for benchmark #1: tanbeta=40; m0=700 GeV; m1/2=600 GeV; A0=-500 GeV; mu>0
#proc BornQCD[fb] BornEW[fb] NloEW[fb] dBornEW[%] dNloEW[%] dAllEW[%]
LL 0.0055 -0.0003713 -0.0004853 -6.75 -8.82 -15.57
RR 0.0080 -0.0002912 -0.0002600 -3.64 -3.25 -6.89
SUM 0.0135 -0.0006625 -0.0007453 -4.90 -5.52 -10.42
#DETAILS
#proc BornQCD[fb] BornEW[fb] NloEW[fb] dBornEW[%] dNloEW[%] dAllEW[%]
u1u1 0.0023 -0.0003358 -0.0003363 -14.44 -14.46 -28.91
u2u2 0.0032 -0.0002828 -0.0001775 -8.59 -5.39 -13.98
d1d1 0.0012 -0.0000449 -0.0002348 -3.75 -19.62 -23.37
d2d2 0.0018 -0.0000115 -0.0000376 -0.62 -2.03 -2.65
c1c1 0.0010 0.0000102 0.0000500 0.98 4.78 5.76
c2c2 0.0015 0.0000033 -0.0000141 0.21 -0.92 -0.71
s1s1 0.0010 -0.0000008 0.0000358 -0.49 3.49 3.40
s2s2 0.0015 -0.0000002 -0.0000308 -0.01 -1.94 -1.96
NNLL /approx. NNLO / Coulomb resummation
bound state effects & gluonia production
R-parity violating sparticle production
--
MichaelKraemer - 27 Oct 2011