Electroweak superpartner production at 13.6 Tev with Resummino

Due to the greater experimental precision expected from the currently ongoing LHC Run 3, equally accurate theoretical predictions are essential. We update the documentation of the Resummino package, a program dedicated to precision cross section calculations for the production of a pair of sleptons, electroweakinos, and leptons in the presence of extra gauge bosons, and for the production of an associated electroweakino-squark or electroweakino-gluino pair. We detail different additions that have been released since the initial version of the program a decade ago, and then use the code to investigate the impact of threshold resummation corrections at the next-to-next-to-leading-logarithmic accuracy. As an illustration of the code we consider the production of pairs of electroweakinos and sleptons at the LHC for centre-of-mass energies ranging up to 13.6 TeV and in simplified model scenarios. We find slightly increased total cross section values, accompanied by a significant decrease of the associated theoretical uncertainties. Furthermore, we explore the dependence of the results on the squark masses.


Introduction
The Run 3 data taking period of the LHC has started in 2022 at an unprecedented centre-of-mass energy of √ S = 13.6 TeV.The two main general-purpose LHC experiments ATLAS and CMS are expected to collect an integrated luminosity of about 300 fb −1 , which will a E-mail: fiaschi@liverpool.ac.uk b E-mail: fuks@lpthe.jussieu.frc E-mail: michael.klasen@uni-muenster.de d E-mail: alexander.neuwirth@uni-muenster.de complement the 140 fb −1 already collected at 13 TeV during Run 2. This increase in luminosity and centreof-mass energy will thus make it possible to further explore extensions of the Standard Model (SM) of particle physics, such as supersymmetry (SUSY) and its minimal incarnation dubbed the Minimal Supersymmetric Standard Model (MSSM) [1,2].The latter is one of the most appealing options for new physics, as it answers various open questions of the SM and provides explanations to several of its conceptual limitations.These include, among others, the presence of dark matter in the universe, the unification of the strong and electroweak forces at large energy scales, and the stabilisation of the mass of the Higgs boson with respect to radiative corrections.
Experimental searches at the LHC have mainly focused on signatures arising from the production and decay of squarks and gluinos, as these particles can be strongly (and thus generally copiously) produced.For a long time, the signatures of these QCD-sensitive superparticles were consequently expected to be the first visible sign of supersymmetry in LHC data.However, with the associated mass limits being now deeply in the TeV regime, searches for electroweakinos and sleptons received more attention and became equally important.Subsequently, accurate theoretical calculations of signal cross sections and key kinematic distributions for all supersymmetric processes became imperative, and in particular for processes in which at least one nonstrongly interacting superpartner is present in the final state.
The Resummino program [3] has been developed in this context.It consists of a public tool computing precision predictions including soft-gluon radiation resummation effects for the production of a pair of sleptons, electroweakinos, and for the associated produc-tion of one electroweakino and either one squark or one gluino.Moreover, Resummino can also be used to calculate cross section predictions for the neutral-current or charged-current production of a pair of leptons in the presence of extra gauge bosons.So far, its predictions have been used by both the ATLAS and CMS collaborations in order to extract bounds on sleptons and electroweakinos.In particular, the most stringent constraints on simplified models inspired by the MSSM enforce viable slepton and electroweakino masses to be larger than about 700 GeV and 800-1200 GeV respectively, for a not too heavy lightest SUSY particle (see e.g.[4][5][6]).The exact values of these mass limits depend on the details of the search channels, and bounds can always be evaded by either compressing the particle spectrum (thus increasing the mass of the lightest SUSY state) or reducing the branching ratio in the final state of interest (by allowing numerous potential decay modes for a given SUSY particle).
Whereas bounds on viable mass regimes for squarks and gluinos imply that their pair production is now phase-space suppressed, their single production with a typically lighter electroweakino is still relevant and could even provide the best insights on both supersymmetric masses and interactions.Associated signals include in particular the associated production of hard jets with missing transverse energy, a signature well studied at the LHC in the context of searches for dark matter [22,23].Associated production rates are known at LO [8] and NLO [24,25], and the matching of these fixed-order results with threshold resummation has been recently achieved [26,27].Finally, fixed-order predictions for charged-current and neutral-current lepton pairs in the presence of additional gauge bosons are known at next-to-next-to-leading-order (NNLO) [28,29], and matched results at NLO+NLL accuracy are available as well [30,31].
In the remainder of this manuscript, we begin by briefly outlining in section 2 the threshold resummation formalism that is implemented in Resummino, and that is relevant for electroweakino-pair, slepton-pair and squark-electroweakino and gluino-electroweakino total production cross section calculations at the LHC.For interested readers, this short description is further complemented by additional details including self-contained analytical formulas in Appendix A. The installation and running of the Resummino code is described in section 3 and in Appendix B-Appendix E. In section 4 and Appendix F, we make use of Resummino to compute and document for the first time total production rates obtained for the LHC Run 3, operating at a centre-of-mass energy of 13.6 TeV.We consider simplified model scenarios in which all superparticles are decoupled, excepted for those produced in the final state.Our results highlight how theoretical uncertainties are reduced relative to the perturbative order of the fixed-order and resummed component of the matched predictions and how they compare with predictions at a centre-of-mass energy of 13 TeV.Moreover, we explore next-to-minimal scenarios, and discuss the impact of internal squark masses on the predictions for configurations in which squarks are not decoupled but only slightly heavier than the lighter electroweakinos.We summarise our work in section 5.

Threshold resummation at aNNLO+NNLL
Within its first public release, the Resummino package was suitable for NLO+NLL resummation calculations in the threshold regime, as well as at small transverse momentum or in both regimes simultaneously.Since then, the implementation of threshold-resummed cross sections has been updated so that corrections up to aNNLO+NNLL could be calculated.For that reason, we provide below a description of the formalism used for threshold resummation up to aNNLO+NNLL, and we refer instead to [3] for details on p T and joint resummation.
We consider the production of a pair of heavy particles i and j in hadronic collisions through the process AB → ij (in which A and B stand for the initial hadrons).Whereas soft and collinear divergences originating from real and virtual corrections cancel in the perturbative expansion of the associated production cross section, logarithmic terms remain due to the different phase spaces inherent in the different ingredients of the calculation [55,56].These logarithmic contributions encode the effects of soft and collinear emission from either initial-state or final-state coloured particles.As the partonic energy approaches the production energy threshold, they become large and can hence spoil the convergence of fixed-order calculations.
This can be cured through QCD resummation techniques that rely on dynamical and kinematic factorisation of the cross section to account for soft and collinear radiation to all orders [57].In practice, these factorisation properties are exploited in Mellin space.The hadronic differential cross section dσ AB /dM 2 is expressed in terms of the Mellin variable N , that is conjugate to the quantity τ = M 2 /S (with M being the invariant mass of the produced particles in the Born process and S the hadronic centre-of-mass energy).This hadronic cross section can be written as a product of the densities f p/H of parton p in hadron H (in Mellin space, where the corresponding Mellin moments are obtained relative to the momentum fraction x), and of the corresponding partonic cross section σ ab→ij (whose Mellin moments are derived relative to the variable z = M 2 /s with s being the partonic centre-of-mass energy), Here the logarithms depend on the Mellin variable N and they become large in the large-N limit, and µ F and µ R stand for the usual factorisation and renormalisation scales.
After accounting for soft-gluon emission to all orders, the partonic cross section can be re-expressed in a closed exponential form scaled by a hard function The N -independent hard function H can be written in terms of the LO Mellin-transformed cross section σ ab→ij and the hard matching coefficient C ab→ij , where the coefficient C ab→ij can be computed perturbatively, The hard matching coefficients C (n) ab are then derived from fixed-order predictions in Mellin space at a given order in the strong coupling α s .They correspond to the ratio of the finite N -independent pieces of the N n LO correction terms over the LO cross section, The soft and collinear gluon radiation contributions appearing in (2) are included in the so-called Sudakov exponent G ab→ij .They depend on the quark/gluon nature of the initial state and can be written at NNLL accuracy as a sum of leading logarithmic (LL) terms G with L = ln(N e γ E ).For all processes implemented in Resummino, we provide in Appendix A explicit expressions for the resummation coefficients G (n) ab→ij and the hard matching coefficients C (n) ab→ij .In principle, all above expressions should be generalised and include an index referring to a given irreducible colour representation, and the total rate given in (2) should embed a sum over all possible structures emerging from the colour decomposition of the LO partonic cross section in Mellin space [58].However, all the processes considered, namely the production of a pair of colourless particles (electroweakino-pair and slepton-pair production) and that of the associated production of one colourless and one coloured particle (electroweakino-squark and electroweakino-gluino production), only feature a single colour representation.All analytical expressions therefore simplify, and the colour representation indices can be omitted.
In order to obtain meaningful predictions over the entire phase space, the resummed cross section has to be matched to its fixed-order counterpart.This is achieved by subtracting from the sum of the resummed total rate σ res.
ab and the fixed-order one σ f.o.ab their overlap, that is given by the expansion σ exp.ab of the resummed result at the same order in α s as that of the fixed-order calculation, In this expression, fixed-order predictions are computed in physical space, whereas the resummed component and its expansion to given order in α s are calculated in Mellin space.This therefore requires to convolve the associated partonic cross sections with Mellintransformed parton distribution functions (PDFs).
In Resummino, this is achieved by employing the parametrisation of the MSTW collaboration [59] to fit the used PDFs in x-space, and thereby obtain an analytic formula for its expression in Mellin space. 1inally, an inverse Mellin transformation has to be performed in order to go back to physical space.This is achieved by choosing a distorted integration contour in the complex plane inspired by the principal value and minimal prescription procedure [64,65].Along this contour, the Mellin variable N is parameterised by so that all singularities that may appear in the inverse transform process are correctly taken care of.The parameter C is indeed chosen so that the poles in the PDF Mellin moments originating from the Regge singularity are on the left of the contour, and that stemming from the Landau pole of the running of α s is on its right.

Installation
The Resummino package is a high-energy physics program whose source code is written in C++, that is publicly available for download at https://resummino.hepforge.org/,and that is licensed under the European Union Public Licence v1.1.It can be compiled with the GNU compiler collection (GCC) and a working installation of CMake (version 3.0 or more recent) [66].It requires as external dependencies the GNU Scientific Library GSL (version 2.0 or more recent) [67] and the Boost package (version 1.70.0 or more recent) [68].The former allows Resummino to make use of the VEGAS Monte Carlo routines for numerical integration [69], and of the Levenberg-Marquardt algorithm [70,71] to fit non-linear functions like parton distribution functions [72].On the other hand, the Boost package is a dependency necessary to read supersymmetric spectra encoded in the SLHA format [73,74] through the SLHAea library [72].In addition, Resummino must be linked to LHAPDF (version 6.0 or more recent) [75] for PDF handling, and LoopTools (version 2.15 or more recent) [76] for the evaluation of one-loop integrals.
The code can easily be installed by typing the following commands in a shell, once the Resummino zipped sources have been downloaded: unzip Resummino-X.Y.Z.zip cd Resummino* mkdir build cd build cmake ..

-B . [options] make make install
The last step is optional, and it is only relevant to install the code system-wide.The cmake command can be cast with options dictating the behaviour of Resummino relative to its two external dependencies LHAPDF and LoopTools.The program is equipped with LHAPDF version 6.2.3 and LoopTools version 2.15, and the usage of those built-in libraries can be enforced through the cmake options -DBUILD_LHAPDF=ON (default: OFF) and -DBUILD_LOOPTOOLS=ON (default: ON) respectively.Different versions of these packages can be used by providing information about the path to existing installations through the options -DLHAPDF=PATH and -DLOOPTOOLS=PATH if libraries and headers are installed in the same folder, or through LHAPDF_INCLUDE_DIR/LHAPDF_LIB_DIR and LOOPTOOLS_INCLUDE_DIR/LOOPTOOLS_LIB_DIR if not.
The Resummino binary, called resummino, is located in the folder build/bin.After completion of compilation, users can test their local installation by running (from the installation folder) ./bin/resummino --help ./bin/resumminoinput/resummino.in The source files of the code are collected in the folder src whose architecture and structure is described in Appendix B. The folder input includes exemplary input files, whereas the folder external contains necessary external packages such as LoopTools and LHAPDF.Finally, the folder scripts provides a set of docker and apptainer scripts related to the usage of portable installations of Resummino.

Running Resummino and input parameters
Resummino can be run as exemplified at the end of the previous section, by typing in a shell, from the installation folder, the command ./bin/resummino<some-input-file> The keyword <some-input-file> provides the path to a configuration file indicating what to calculate and how.Details on all the options available to write such an input configuration file are briefly provided in the rest of this subsection, and they are additionally extensively documented in Appendix C.Moreover, an example input file input/resummino.in is shipped with the code, and the output obtained from running the code with this input file is described in Appendix D.

Calculation, process and collider settings
A configuration for Resummino includes a list of equalities fixing the values of certain variables acting on how the code functions.This file begins with a definition of the process to be considered, including details on the collider environment.Users must specify the type of colliding particles (protons or antiprotons), the centre-of-mass energy, and the nature of the final-state particles.This is achieved through the self-explanatory variables collider_type (to be set to proton-proton or proton-antiproton), center_of_mass_energy (to be given in GeV), particle1 and particle2 (provided through their Particle Data Group (PDG) identifiers [77]).For example, the following settings, collider_type = proton-proton center_of_mass_energy = 13600 particle1 = 1000024 particle2 = -1000024 would correspond to using Resummino for calculations relevant to LHC proton-proton collisions at a centre-ofmass energy of 13.6 TeV, and leading to the production of a pair of charginos (pp → χ+ ✓ pp → χq [27] pp → χg [26] Table 1 List of processes included in Resummino, given together with the orders of the different calculations supported by the code (defined as the value of the result variable discussed in the text) and the associated references.
joint resummation formalism respectively), or m (differential cross section at fixed invariant mass M according to the threshold resummation formalism).Differential calculations make use of the variables pt and M to get the numerical values to employ for p T and M , respectively, that are provided in GeV.For example, result = total defines a total cross section calculation, whereas estimates dσ/dp T for a final-state transverse momentum of 50 GeV, using the joint resummation formalism (including therefore integration upon the invariant mass M ).In addition, for the Drell-Yan production of a lepton pair, an invariant-mass cut is required to regularise phase-space integration when total rates are evaluated.Information on such a cut is passed through the parameters Minv_min and Minv_max.The list of calculations supported in Resummino is summarised in table 1, together with the key references documenting their implementation.The numerical precision of the performed calculation can be controlled through the input variables precision and max_iters acting on the VEGAS algorithm.The former is related to the relative precision of the numerical integration process, whereas the latter refers to the maximum number of integration iterations (excluding the warm-up phase) allowed before stopping the calculation.The settings precision = 0.01 max_iters = 5 enforces a numerical precision of 1% and imposes that at most five VEGAS iterations are used (even if the desired precision is not reached).
Next, details about the PDF set to be used in LO and higher-order calculations are provided through a few variables.These variables define the name of the LHAPDF set of parton densities to employ (through the self-explanatory variables pdf_lo and pdf_nlo2 ), and the exact identifier of the set member considered (specified as an integer through the variables pdfset_lo and pdfset_nlo).The cross section values spanned after considering all set members included in a given LHAPDF set (obtained by means of multiple runs of Resummino) then allow for PDF error assessments following standard Hessian or Monte Carlo prescriptions [78].For instance, including in the input file enforces the usage of CT14 parton densities [60] for all calculations.The central LO set CT14lo will be used for LO predictions, whereas both higher-order fixed-order predictions and resummed ones will rely on the central NLO set CT14nlo.
As briefly discussed in section 2, Mellin-transformed PDF are obtained from a fit of the PDFs considered to the MSTW parameterisation of [59].This relies on a logarithmic sampling of Bjorken-x values from a minimum value x min to 1, and on a fitting method using weighted least squares as achieved by the Levenberg-Marquardt algorithm (handled through the GSL library).The minimum value x min can be provided through the input variable xmin, and the weights used in the fitting procedure are given through the variables weight_valence, weight_sea and weight_gluon for valence quark, sea quark and gluon PDFs respectively.The default option corresponds to Finally, the value of the renormalisation and factorisation scales is taken to be µ for the computation of total cross sections and p T distributions, with m 1 and m 2 being the masses of the two final-state particles.For invariant-mass distributions, the choice µ F = µ R = M is adopted instead, unless the FIXED_SCALE flag is turned on in the file src/resummino.cc.As required for the evaluation of scale uncertainties, users have the possibility to multiply these values by specific factors via the input variables mu_f and mu_r.For instance, a calculation corresponding to µ F = (m 1 + m 2 )/4 and µ R = (m 1 + m 2 ) would be achieved through mu_f = 0.5 mu_r = 2.0

Model free parameters
Information on the SUSY spectrum and interactions is passed to Resummino through a standard SLHA [73,74] file whose path is specified in the input file through the variable slha.In the case of the charged-current or neutral-current Drell-Yan process in the presence of additional gauge bosons, information on the SM electroweak parameters and on the W ′ and Z ′ properties is provided through a file encoded following an SLHA-like structure, and whose path is passed through the variable zpwp as described in [31].An example of such a file and its structure is provided in input/ssm.in.
By default, the electroweak input parameters are determined from the information included in the provided SLHA file. 3 For SUSY calculations Resummino makes use of tree-level formulas, following the prescription introduced in Prospino [79].In these expressions, M W and M Z stand for the Wboson and Z-boson masses, G F is the Fermi constant, g 2 is the weak coupling constant, and θ w is the electroweak mixing angle.All Yukawa couplings of light quarks (for a number of active quark flavours n f = 5) are set to zero, and the strong coupling constant α s is computed using the value returned by LHAPDF.Consequently, Resummino does not use the values of α s and α included in the SMINPUTS block of the SLHA input file.For computations in the presence of extra gauge bosons, the code uses instead where g 1 stands for the hypercharge coupling constant.
Users have the possibility to change the electroweak parameter scheme through the different macros implemented in the file src/params.cc.For instance, the flag GAUGE_COUPLING can be switched on to use as inputs the values of the SM gauge couplings encoded in the SLHA block GAUGE, or the flag MADGRAPH_COUPLING can be turned on to select the electroweak parameters as in [11].

The command line interface of Resummino
Several of the configuration settings of the code can be fixed through optional arguments when casting the Resummino executable command in a shell, as in ./bin/resummino<input-file> [options] Available options facilitate, for example, iterations over different PDF sets, calculations at different perturbative orders, or the implementation of similar calculations for various processes.Users do not therefore have to manually modify the input file (generically denoted by <input-file> in the above command), as any parameter value passed through the command line supersedes that included in the input file.The full list of options being detailed in Appendix E, we only briefly introduce them below.
The perturbative order of the calculation can be fixed through the options --lo (or -l, for LOaccurate calculations), --nlo (or -n, for NLO-accurate calculations) and --nll (or -s, for NLO+NLLaccurate calculations), as well as through --nnll (or -z) for aNNLO+NNLL-accurate calculations when available (namely for implemented Drell-Yan-like processes).Casting the Resummino command with the --center_of_mass_energy flag (or -e) followed by a double-precision number allows for a modification of the collider hadronic centre-of-mass energy value (in GeV).Similarly, the p T and M values (in GeV) relevant for calculations at fixed transverse momentum or fixed invariant mass can be fixed through the two options --transverse-momentum (or -t) and --invariant-mass (or -m), both followed by a doubleprecision number.Additionally, the PDG identifiers of the final state particles can be modified through the options --particle1 (or -c) and --particle2 (or -d), both followed by an integer.
The evaluation of the theory uncertainties associated with any given calculation can be easily performed through the modification of the PDF set considered and the chosen unphysical scales.This is achieved through the options --pdfset_lo (or -a) and --pdfset_nlo (or -b), both followed by an integer, and --mu_f (or -f) and --mu_r (or -r), both followed by a doubleprecision number.The former two flags allow users to change the identifier of the PDF member set (within a specific LHAPDF collection of parton densities) used for calculations at LO and beyond, whereas the latter two flags provide multiplicative factors to include when fixing the factorisation and renormalisation scales.
In addition, the code can also be started with the options --version (or -v), --help (or -h), --parameter-log (or -p) followed by a string, and --output (or -o) followed by a string.The first of these displays the Resummino release number to the screen, whereas the second of them prints a help message indicating how to run the code.The third possibility makes the code writing all the parameters inherent to the calculation considered in a file (whose path is provided through the included string), and the last one defines the directory (defined through the included string) in which all files created by the code on run time are stored.

Docker containers
We provide convenient 'docker' containers that allow users to employ Resummino without the need to compile it from its source code.This requires to have either Docker or Singularity/Apptainer available on the system, the latter choice being recommended for secured high-performance computing.In this case, Resummino can be run from the scripts localised in the scripts folder, # if Docker is available source scripts/docker_alias.sh # if Singularity/Apptainer is available source scripts/apptainer_alias.shThese scripts automatically detect a local LHAPDF installation, and then enable the usage of any PDF set available within the system.If no installed version of LHAPDF is found, then Resummino is restricted to make use of the default PDF sets included within the docker images.New PDFs can always be added via the command 'lhapdf install', to be cast within a shell.In addition, when it is run from a script, Resummino only accepts input files located in the current directory.Such a behaviour can however be circumvented through modifications of the docker run command's bound directories, through its usual -v option.

Precision predictions for electroweak SUSY processes in simplified models
In this section we provide state-of-the-art cross section predictions for slepton and electroweakino pair production at the LHC Run 3 operating at an increased centreof-mass energy √ S = 13.6 TeV, and we compare our findings with predictions relevant for the LHC Run-2 at √ S = 13 TeV.We employ the PDF4LHC21_40 set of parton distribution functions [80], and we provide results together with the associated PDF and scale uncertainties (the latter being obtained with the seven-point method) added in quadrature.Complementary to the figures shown in this section, the complete collection of numerical predictions are shown in the tables of Appendix F.
We consider simplified SUSY models inspired by the MSSM, and we explore several typical scenarios.In the context of slepton pair production, we focus on a new physics configuration in which all SUSY particles are decoupled by setting their masses at 100 TeV, with the exception of a single slepton species that is taken either left-handed (ẽ L ), right-handed (ẽ R ) or maximally mixed (τ 1 = 1/ √ 2 τL + τR ).This last scenario is representative of models featuring light tau sleptons, as originating from many SUSY scenarios [19].
Electroweakino pair production rates are estimated in similar scenarios, in which all SUSY particles are decoupled with the exception of the lightest electroweakinos [13].We begin our study with scenarios in which the three lightest electroweakinos are all higgsinos.In a first setup, all higgsinos are taken mass-degenerate and the lightest states are defined by Whereas this choice of a real neutralino mixing matrix implies a negative m(χ 0 1 ) eigenvalue, it can always be transformed back to positive a mass eigenvalue through a chiral rotation [73].In a second setup, we introduce some mass splitting between the three electroweakinos, such a splitting being typical of next-to-minimal electroweakino simplified models studied at the LHC, and that turn out to be more realistic in the light of concrete MSSM scenarios [81].In this case, we define the three lightest electroweakinos by Finally, we consider a scenarios in which all lightest electroweakinos are mass-degenerate gauginos, In this last SUSY configuration, we additionally investigate squark mass effects on gaugino pair production.Here, we consider an eight-fold degeneracy of all first-generation and second-generation squarks, the spectrum featuring thus a large number of states reachable at the LHC.TeV and 13 TeV (upper insets), shown together with their ratios to the 13 TeV total rates (lower insets) in which combined scale and PDF uncertainties are included.We consider the production of a pair of left-handed sleptons (top), right-handed sleptons (centre) and maximally-mixed sleptons (bottom), and predictions are presented as a function of the slepton mass m l.

Slepton pair production
In figure 1 we display aNNLO+NNLL cross section predictions for slepton pair production at the LHC as a function of the slepton mass m l.We consider two centre-of-mass energies fixed to √ S = 13.6 TeV (orange) and 13 TeV (blue), and we focus in the upper, central and lower panel of the figure on the respective processes In the figures, we restrict the mass range shown to m l ≲ 1 TeV.This corresponds to cross section values larger than 0.01 fb, to which the LHC Run 3 is in principle sensitive as dozens signal events could populate the signal regions of the relevant ATLAS and CMS analyses. 4hereas cross sections for √ S = 13 TeV and 13.6 TeV are both shown in the upper insets of the three subfigures, the gain in rate at Run 3 is more visible from the ratio plots presented in their lower insets.This indeed illustrates better how cross section increases ranging up to 20% can be obtained in the three classes of scenarios considered, especially for large slepton masses.As expected from the structure of the slepton couplings to the Z-boson (see e.g. in [32]), left-handed sleptons are more easily produced in high-energy hadronic collisions than their righthanded counterparts that only couple through their hypercharge.Consequently, cross sections corresponding to mixed scenarios lie between the two extreme non-mixing cases for a given slepton mass m l.
The different ratio plots of the lower insets of the subfigures also show the dependence of the theoretical systematic uncertainty bands on the slepton mass.The collider energy upgrade achieved at Run 3 naturally leads to a reduction of the PDF uncertainties as the gain in centre-of-mass energy yields a smaller relevant Bjorken-x regime in which parton distribution functions are better fitted.As shown in table 3 (see Appendix F), scale uncertainties contribute to at most 2% of the combined theoretical uncertainty in the entire mass range probed, regardless of the centre-of-mass energies considered.In contrast, PDF errors vary from 3% to 18% at 13.6 TeV, which must be compared to a variation ranging from 3% to 20% at 13 TeV.For a phenomenological study on the reduction of PDF uncertainties in slepton pair production see [20].

Higgsino pair production
We now turn to scenarios in which all SUSY particles are decoupled by setting their masses at 100 TeV, with the exception of all higgsino states.We begin with a calculation of aNNLO+NNLL predictions relevant for higgsino pair production in a scenario in which all higgsinos, defined as in eq. ( 11), are mass-degenerate, i.e. in which Our results are shown in figure 2 for the four processes since in the case of a degenerate spectrum σ( χ0 1 ).In the results displayed, we restrict the mass range considered to m χ ≲ 1.5 TeV, which corresponds to production rates at the LHC larger than 0.01 fb and therefore potentially reachable at Run 3. The largest cross sections are obtained for the charged-current process pp → χ0 1 , such an effect originates from a PDF enhancement related to the ratio of valence and sea quarks in the proton and from the structure of the higgsino gauge couplings.This additionally leads to similar neutral-current higgsino production rates (pp → χ+ 1 ), and the cross section of the negative charged-current process pp → χ0 1 χ− 1 is then smaller.As in the slepton case explored in section 4.1, we find an enhancement of total production cross sections at 13.6 TeV relative to those at 13 TeV thanks to the modest gain in phase space.Rates are indeed found to be 10% to 30% larger at √ S = 13.6 TeV than at √ S = 13 TeV, for low and high electroweakino masses respectively.
Still similarly to the slepton case, theoretical uncertainties get reduced with the increase in centre-of-mass energy.As shown in table 4 (see Appendix F), scale uncertainties negligibly contribute to the total theory errors for both centre-of-mass energies, scale variations indeed leading to errors of about 1% -2% for low higgsino masses and lying in the permille range for m χ ≳ 300 GeV.In contrast, total rates at 13 TeV are plagued with PDF uncertainties varying from a few percent at low masses to more than 20% -40% for higgsino masses larger than about 1.2 TeV.The reduction of the average Bjorken-x value inherent to the larger centre-ofmass energy of 13.6 TeV subsequently leads to smaller PDF errors that are found reduced by about a few permille at low masses, to up to 5% -7% at large masses.For a phenomenological study on the reduction of PDF uncertainties in higgsino pair production see [14].We now move on with a second higgsino scenario in which the three lightest higgsino states are defined as in eq. ( 12).Moreover, their spectrum is enforced to feature a significant level of compression, so that all three higgsino states exhibit a mass splitting of a few percent 5 .In the following, we impose that with all masses being taken positive.Given these mass relations, we present in figure 3 aNNLO+NNLL total cross sections for the process pp → χ0 1 χ+ 1 at the LHC, for centre-of-mass energies of √ S = 13 TeV (left) and 13.6 TeV (right).We consider the mass range in which all higgsinos are lighter than 300 GeV, as this consists of the relevant 5 The masses of the decoupled squarks are set to 4.5 TeV in order to avoid numerical complications.mass configurations in terms of LHC sensitivity to compressed SUSY higgsino scenarios [82][83][84].
In the figures, associated rates are shown logarithmically through a colour code.In this scheme, other higgsino production modes yield almost identical figures, that we therefore omit for brevity.Numerical results are nevertheless provided for all processes in table 5 (see Appendix F), together with separate scale and PDF uncertainties.For the considered mass range, theoretical systematics are found in very good control, the combined uncertainties being of about 5% for all mass configurations explored.Whereas scale uncertainties decrease from 2% -3% in the lightest configurations considered to a few permille for higgsinos of about 200 GeV -300 GeV, PDF errors increase from 2% -3% in the lightest scenarios to 4% -5% in the heavier cases.The combined theory errors are thus similar in size for all scenarios studied.As for the previous calculations achieved, a cross section increase at 13.6 TeV results from the phasespace enhancement inherent to the increased centre-ofmass energy relative to the 13 TeV case, these findings being numerically testified by the results displayed in table 5.Moreover, a rate hierarchy similar to that observed in the mass-degenerate case is obtained, phasespace effects being minimal for compressed scenarios with a non-degenerate spectrum compared to massdegenerate scenarios in which all higgsinos have exactly the same mass.The process pp → χ0 1 χ+ 1 hence dominates, followed by the neutral current modes (pp → χ0 1 .We remind that such a hierarchy is dictated by the structure of the higgsino gauge couplings to the W and Z bosons, and by the PDF ratio of valence and sea quarks in the proton.

Gaugino pair production
We now turn to the analysis of the gaugino scenarios introduced in eq. ( 13), in which the lightest electroweakinos consist of bino and wino eigenstates.In our analysis, we consider that the two wino eigenstates are massdegenerate, with respective masses satisfying In the following, we then focus on the processes and we present aNNLO+NNLL predictions for the associated total rates in figure 4. Other processes are irrelevant as the corresponding cross sections vanish due to the structure of the bino and wino gauge couplings.
As in the previous subsections, we consider results at centre-of-mass energies of 13 TeV (blue) and 13.6 TeV (orange) in the upper insets of the figures.This time, however, we display predictions for wino masses ranging up to 2 TeV, the cross sections being much larger than in the higgsino case by virtue of the weak triplet nature of the winos.Consequently, we can expect a better LHC sensitivity to signatures of wino production and decays, due to the machine being capable to naturally probe a larger mass regime.
In accord with parton density effects, the chargedcurrent process pp → χ0 2 χ+ 1 dominates for a given wino mass, its rate being a factor of 1.5 to 3 larger than that of the other charged-current process pp → χ0 2 χ− 1 in the case of lighter and heavier mass setups respectively.Furthermore, total cross sections for the neutral current process pp → χ+ 1 χ− 1 , mediated by virtual photon and Z-boson exchanges, are usually 1.25 -1.5 smaller than rates corresponding to the charged-current mode pp → χ0 2 χ+ 1 , in which the final state is produced from virtual W -boson exchanges.On the other hand, the increase in cross section observed when the hadronic centre-ofmass energy is modified from 13 TeV to 13.6 TeV can be quite substantial in such mass-degenerate wino scenarios.While for light produced particles the increase is only modest and lies in the 5% -10% range for all three processes, it increases with the wino mass m χ and reaches 35% -40% for wino masses of about 2 TeV.
In the lower insets of the three subfigures, we display the ratio of the total production rates at 13 TeV and 13.6 TeV to that at 13 TeV, including the combined theory systematic error.We remind that numerical values for the cross sections are presented together with separate scale and PDF uncertainties in table 6 (see Appendix F).For the entire mass range considered, scale uncertainties are under good control.They are about 1% -2% at low masses, and then decrease to a few permille for wino masses larger than 300 GeV.In contrast, PDF errors are smaller than 10% for wino masses smaller than about 1 TeV, but quickly increase for heavier mass configurations.In this case, typical Bjorken-x values are large and correspond to phase space regimes in which parton densities are poorly constrained, as already pointed out in [11].This issue will nevertheless be automatically cured with time.Time will indeed allow the LHC collaborations to collect better-quality SM data at large scales, which will consequently help in reducing the PDF errors.
In order to explore the phenomenological consequences of next-to-minimality, we now consider scenarios in which first-generation and second-generation squarks are mass-degenerate, but not decoupled.We introduce the squark mass parameter m q defined by that we then vary between 800 GeV and 4 TeV.For illustrative purposes, we focus on the charged-current process that gives rise to the largest wino production cross sections for a specific mass spectrum.The discussion and the results below are however applicable to other wino production modes as well.This is further numerically depicted in table 7 (see Appendix F).The latter indeed includes a complete set of numerical predictions for all wino pair production processes in the presence of not too heavy squarks, once again together with separate information on the total rate values, and the associated scale and PDF uncertainties presented separately.
In figure 5 we present aNNLO+NNLO total cross sections for the process of eq. ( 21) as a function of the wino and squark masses m χ and m q6 .Results are presented for collider energies of √ S = 13 TeV and 13.6 TeV in the top and bottom row of the figure respectively.In the left subfigures, the rates are represented through a logarithmic colour code, which shows that they exhibit a non-trivial dependence on the squark mass.Due to the destructive interference of s-channel gauge boson exchange diagrams with t/u-channel squark exchange diagrams, the cross section starts by decreasing when the SUSY spectrum is varied from a configuration in which m χ = m q to    one with a larger squark mass value still in the vicinity of the same wino mass m χ.The rate then gets larger and larger with increasing squark masses (the wino mass being constant), and it finally saturates when squarks decouple.This feature is further illustrated in the two upper insets of the right subfigures, that display aNNLO+NNLL production cross sections for χ0 2 χ+ 1 production at the LHC as a function of the squark mass for several choices of wino masses.For a given wino mass m χ, we observe that the cross section always begins by decreasing before quickly reaching a minimum, and then increases for larger and larger squark masses.In the limit of very heavy squarks, the latter decouple and the rates become independent of the squark properties.They hence solely depend on the wino mass.
In the lower inset of these two right subfigures, we present K-factors defined as the ratio of the most pre-cise aNNLO+NNLL rates to the LO ones for a given mass configuration, This further illustrates the squark decoupling at large m q values, the K-factor becoming constant.Moreover, the dependence of the K-factors on the squark mass additionally shows how the minimum of the cross section, for a given wino mass, is shifted by tens of GeV by virtue of the higher-order corrections.In all the insets included in the right subfigures, we additionally include combined theory uncertainties.In general, those uncertainties, that are numerically given in table 7, are under good control, with the exception of the heaviest scenarios in which the PDF errors get very large, as already above-mentioned.On the contrary, scale uncertainties always lie in the percent or permille level for all scenarios considered, the error bars being in fact drastically impacted (and reduced) by QCD resummation.

Conclusion
We have updated the documentation of the Resummino package dedicated to precision calculations of total rates and invariant-mass and transverse-momentum spectra in the context of slepton and electroweakino pair production, the production of a pair of leptons in the presence of additional gauge bosons, and that of an associated electroweakino-squark or electroweakinogluino pair.This update of the documentation is motivated by the significant extensions that have been implemented since the initial release of the code a decade ago, and that have never been collected in a single document.Whereas these new features of Resummino have already been used for various phenomenological studies and experimental searches for particles beyond the Standard Model at the LHC, we have taken the opportunity of this paper to provide a new useful illustration of the capabilities of the program.We have computed and tabulated for the first time precision predictions of total rates relevant to searches for sleptons and electroweakinos at the ongoing Run 3 of the LHC, operating at a centre-of-mass energy of √ S = 13.6 TeV.We have considered several simplified models inspired by the MSSM in which all superpartners are decoupled, with the exception of a few states whose production at the LHC has been explored.We have presented predictions matching fixed-order calculations at approximate NNLO in QCD (in which SUSY twoloop contributions have been neglected) with predictions including the resummation of soft-gluon radiation at NNLL.In our results, we have put emphasis on the gain obtained from updating the centre-of-mass energy from 13 TeV to 13.6 TeV, and on the theory uncertain-ties inherent to all calculations achieved.A brief summary of our findings is given in table 2.
As a rule of thumb valid for all investigated processes, QCD threshold resummation reduces the typical scale dependence to less than 1% at aNNLO+NNLL.These contributions to the combined theory error hence become generally negligible, or at least subleading.For the production of low-mass SUSY states, PDF and scale uncertainties are of the same order of magnitude so that theory systematics are under good control.In our predictions, we have employed the PDF4LHC21 set of parton densities that originates from a combination of variants of the CT18 [61], MSHT20 [62] and NNPDF3.1 [85] global PDF fits.Such a set generally leads to larger PDF uncertainties compared to the individual global sets, particularly at large Bjorken x ≳ 0.4 [80].In this region, yet poorly constrained by experimental data, discrepancies between predictions relying on the different global sets can be significant as the parton density behaviour is mainly driven by the theoretical assumptions made.For instance, recent studies have shown how the recent NNPDF4.0 [63] set predicts (anti)quark distributions falling much faster (slower) with respect to the other sets [86].Moreover, the sea over valence quark ratio is a determining factor for Drell-Yan-like processes [87], and in turn it also has a strong impact on the determination of the gluon density at high x values and large scale [88].Predictions for (differential) cross section can therefore vary substantially with the choice of the PDF set, especially in scenarios featuring heavy SUSY particles.This issue will however be naturally fixed with time, as more data gets collected at high scale by the LHC collaborations.
We have observed that the modest increase in centre-of-mass energy at the LHC Run 3 brings a typical increase in the rates associated with colourless SUSY processes of 20% for light spectra to 40% in the heavier cases considered.In addition, we have explored the impact of next-to-minimality following two aspects.First, we have considered the production of a pair of non-degenerate compressed higgsinos, and next we have studied the impact of squark decoupling on wino pair production.Both setups have led to visible effects, so that their impact will have to be considered for non-minimal excursions in the interpretation of future search results.These findings additionally pave the way to future work.In particular, it will be important to assess the phenomenology of non-minimal model spectra exhibiting several light strongly-interacting and non-stronglyinteracting SUSY states.In this case, rarely-studied semi-strong processes (such as squark-electroweakino associated production) could be relevant, and worth to be explored.
The package Resummino and its current release 3.1.2is ready to achieve related calculations.The source code is publicly available from HEPForge and can be used for aNNLO+NNLL calculations for slepton and electroweakino pair production (as studied in this article), as well as for NLO+NLL calculations for the associated production of a squark-electroweakino and gluino-electroweakino pair.Moreover, nonsupersymmetric models featuring additional W ′ and Z ′ boson could also be investigated at aNNLO+NNLL accuracy, and in addition to total rates Resummino can be used to calculate invariant-mass and transversemomentum distributions for all implemented processes (with an NLO+NLL precision).

Appendix A: Details on the threshold resummation formulas implemented in Resummino
In this appendix we collect analytic expressions for the exponentiated logarithmic contributions and the hard matching coefficients relevant for all processes supported in Resummino.
The logarithmic coefficients of ( 6) are given, for a Born process initiated by the initial-state partons a and b by [35-37, 40, 41] where λ = α s b 0 L with L = ln N and N = N e γ E , and where b 0 is defined from the first coefficient of the QCD beta function.We provide its expression together with that of the next two coefficients of the QCD beta function, that are relevant for the formulas given below.They are normalised according to b n = β n /(2π) n+1 [89,90], which gives for quarks (c = q) and gluons (c = g) only depend on the nature of the initial state.They are given as compact functions of the universal processindependent resummation coefficients A (n) c , the first two of these being written as with C c = C F for quarks and C c = C A for gluons.As in Resummino NNLL precision is only achieved for slepton and electroweakino pair production, the implemented formulas only require the knowledge of the third universal quark coefficient A q , (A.4) The first two resummation coefficients g whereas the third quark coefficient (the only coefficient needed for the processes implemented at aNNLL+NNLL in Resummino) is given by The process dependent term h (2) ab→ij appearing in (A.1) can be expressed in terms of the soft anomalous dimension Γ ab→ij associated with the partonic process ab → ij, from which the Drell-Yan soft anomalous dimension Γ DY ab has been subtracted.The coefficient with For Drell-Yan-like processes computed at the aNNLO+NNLL accuracy (namely electroweakino and slepton pair production), the coefficient h (3) ab→ij must be included too.It is in this case universal, and it is given by (A.9) In the next subsections, we collect the processdependent hard matching coefficients appearing in ( 3) and ( 4), expressions for the modified soft anomalous dimension Γab→ij , as well as the ingredients allowing for the computation of the resummed cross section at a given order in α s , the latter being necessary for the matching of the resummed and fixed-order predictions shown in (7).
Appendix A.1: Supersymmetric Drell-Yan-like processes This subsection collects the formulas relevant for the purely electroweak production of pairs of sleptons and electroweakinos.
These processes do not feature any coloured finalstate particle, and they only involve a quark-antiquark initial state at Born level.The associated hard matching coefficients C (n) q q are hence identical for all Dell-Yanlike processes, the process dependence being factorised in the Born cross section in (3).The first two C (n) q q coefficients can be written in a compact form in Mellin space, C ab→ij therefore vanishes for all Drell-Yan-like processes.
whereas the third coefficient is For SUSY Drell-Yan-like processes, fixed-order predictions are available at aNNLO in QCD.Calculations hence include SM NNLO contributions, but neglect NNLO SUSY-QCD corrections.The latter are however expected to be small, given current bounds on SUSY particles.The matching of the resummed cross sections at NNLL with the aNNLO fixed-order ones is performed by expanding the resummed cross section up to O(α 2 s ).Starting from the expressions given in (2), ( 3) and (4), we get ab + . . .
ab K (1) + . . . . (A.12) In the above expressions, we have omitted the arguments of the various functions from the second equality onward, in order to simplify the notation.The coefficients K (n) of the expanded exponential factor can be organised in powers of the logarithm L = ln(N e γ E ), with the various N -independent contributions being given by Appendix A.2: Gluino-electroweakino associated production The soft anomalous dimension Γ q q→g χ associated with the process q q → g χ can be calculated following standard techniques, from eikonal Feynman rules.The corresponding modified soft anomalous dimension Γqq→g χ is given, after that the Drell-Yan contribution has been subtracted, by For such a process and predictions at NLO+NLL, the hard matching coefficients can be derived from the collinear remainder originating from the P and K insertion operators defined in the language of dipole subtraction [91][92][93].This is achieved by restricting these operators to their N -independent pieces, and by adding to the result the finite part of the virtual corrections, C q q = 2 P+K q,N -ind.
(A. 16) In this expression, the arguments of the various functions are omitted for brevity.Moreover, σ V stands for the virtual corrections contributing at NLO, and σ A consists of the auxiliary cross section shifting the divergences so that integration over phase space can be numerically achieved.Explicit expressions for the insertion operators P and K for an initial quark are where the Mandelstam variables are defined according to the ordering of the two-to-two partonic process q q → g χ.Similar equations hold for u-channel contributions, that must be included as well when computing the ⟨P+ K⟩ q term in (A.16), Appendix A.3: Squark-electroweakino associated production Similarly to Appendix A.2, the modified soft anomalous dimension Γqg→q χ associated with the process qg → q χ is obtained from eikonal Feynman rules, and it is given by The hard-matching coefficients are similarly derived from the N -independent parts of the insertion operators P and K, + P+K g,N -ind.
in which we use the explicit expressions The two quantities Q and G appearing in these formulas are given by where the Mandelstam variables are defined from the two-to-two partonic process qg → q χ.Appendix B: Architecture and structure of the Resummino sources files The source code of Resummino is available from the folder src of any local installation of the code.It includes the following set of files, that we briefly describe.
constants.h:Constants useds throughout the code, such as colour factors and the regulator width used for on-shell subtraction [94].
hxs.cc: Phase space implementations and integration routines relevant for resummed cross sections matched to NLO fixed-order predictions evaluated within the dipole formalism.It thus includes separate components for the partonic LO cross section, the associated virtual corrections (to which integrated dipoles are added), the collinear remainders, the real emission corrections (from which dipole contributions are subtracted), and threshold resummation components.The integration relevant for total rate calculations is performed through the function hadronic_xs.Integrations associated with transverse momentum and invariant mass differential cross sections are similarly implemented in the files hxs_dpt2.ccand hxs_dlnm2.ccrespectively.
kinematics.cc: Calculation of the kinematics inherent to the different partonic configurations relevant for a given process, and momentum reshuffling as required by the regularisation of infrared divergences in the dipole subtraction formalism.
-integration_method.cc: General-purpose wrapper for Monte Carlo integration through the GSL routine VEGAS.
main.cc: Initialisation of the code within its main function.
options.h:Different options for the implementation of on-shell resonance subtraction from the real emission corrections to squark-electroweakino associated production.
params.cc:Values of all SUSY and SM couplings and masses (see also details in section 3.2.2).
pdf.cc: Initialisation of the PDF functions, interface with LHAPDF, and fit of the PDFs in Mellin space.
resummation.cc:Implementation of the resummed exponent defined in (6), and the hard-matching coefficients defined in (5).
resummino.cc:Parsing of the input file and calls to the functions appropriate for the process and calculation chosen (total rate, p T distribution, etc.).
utils.h:Additional functions allowing for internal checks.
Specific functions relevant for each of the supported processes, and including Born, NLO-real and NLOvirtual contributions, are collected in the sub-folders gaugino_gluino ( χg production), gaugino_squark ( χq production), gauginos ( χ χ production), sleptons ( ll production) and leptons (dilepton production including Z ′ and W ′ exchange contributions).

Appendix C: Configuration of the Resummino run
In this appendix, we describe how to write a configuration file for Resummino.Such a file includes a set of equalities allowing to set keywords to user-defined values.A first ensemble of keywords is dedicated to the definition of the hadronic process and the type of calculation to be performed.They are given by: -collider_type = <string>: information on the initial-state hadrons.Allowed values are proton-proton or proton-antiproton.
result = <string>: defines the cross section to compute.Allowed values are total (total cross section using the threshold resummation formalism), pt (differential cross section dσ/dp T using the p T resummation formalism), ptj (differential cross section dσ/dp T using the joint resummation formalism) and m (differential cross section dσ/dM using the threshold resummation formalism).
pt = <double>: transverse momentum in GeV, used for differential cross section calculations at fixed p T .This variable is only relevant when the result parameter is set to pt or ptj.-M = <double>: invariant mass in GeV, used for differential cross section calculations at fixed M .The invariant mass of the final-state system is used by enforcing the specific value auto.
-Minv_min = <double> and Minv_max = <double>: invariant-mass window, given in GeV, that is used for total cross section computations relevant for the production of a Drell-Yan pair of leptons.Such parameters can also be set to auto, which corresponds to M min = (3/4)M Z ′ and M max = (5/4)M Z ′ .
Information on the SUSY properties or on the extra gauge boson properties is provided through files encoded in an SLHA(-like) format.Two variables allow users to provide information about the path to these files: slha = <string>: path to the SLHA input file containing information on the SUSY scenario considered (masses, widths, parameters, etc.).
zpwp = <string>: path to the SLHA-like input file containing information on the properties of extra Z ′ and W ′ bosons.
A third ensemble of keywords are related to the choice of the parton density sets to be used in given calculations, and how their Mellin transform is handled.They are: -pdf_format = <string>: format according to which the LHAPDF grids are encoded.Allowed values include LHpdf (default) and LHgrid (deprecated option relevant for LHAPDF versions anterior to 6).
-pdf_lo = <string> and pdf_nlo = <string>: name of the PDF set to be employed for LO and higher-order calculations respectively, according to the naming scheme of the webpage https://lhapdf.hepforge.org/pdfsets.-pdfset_lo = <int> and pdfset_nlo = <int>: specific PDF set member to be used in LO and higher-order calculations respectively.These integers allow for the selection of the best-fit set or of one of the next-to-best-fit sets provided within a given LHAPDF set, in the context of PDF error estimates.
xmin = <double>: minimum Bjorken-x value to consider in the PDF fitting procedure.When set to auto (default), the code uses x min = M 2 /S where M is the invariant mass of the pair of produced particles.
Renormalisation and factorisation scales are automatically set to specific values for any calculation achieved by Resummino (see section 3.2.1).Those values can however be multiplied by user-defined factors given through two self-explanatory variables: -mu_f = <double> and mu_r = <double>: factors multiplying the factorisation and renormalisation scale values specific to the calculation considered.Default values are 1.0.
Finally, the numerical precision of the Monte Carlo integration achieved with the VEGAS algorithm is controlled through two input variables: precision = <double>: relative precision of the Monte Carlo integration (recommended value: 0.001 ≡ 0.1%).The integration routine stops when the relative numerical error of the result matches the requested precision.For contributions beyond LO, this error is estimated with respect to the tree-level result, which ensures that no computing power is wasted on small contributions.
-max_iters = <int>: maximum number of integration iterations, excluding the warm-up phase (recommended value: 50).After the maximum number of iterations, Monte Carlo integration is stopped, even if the desired precision is not reached.

Appendix D: Resummino output
In this appendix, we detail the (screen) output generated by Resummino in the context of the calculations of the total cross section associated with slepton pairproduction at the LHC, operating at a centre-of-mass energy of 13 TeV.We choose to make use of the LO and NLO sets of CT14 parton densities [60], and to rely on one of the SUSY scenarios whose SLHA input file is provided in the input folder of the Resummino source files.
The used input file including all computational settings is available from the folder input/resummino.in.We remind that several entries (like zpwp or Min_min and Minv_max) are irrelevant for the SUSY process considered.The code is started by typing in a shell ./bin/resumminoinput/resummino.inWhen executed, Resummino begins by printing to the screen its banner (that includes information on the articles to cite when using the code), as well as LHAPDF settings (path to the PDF grids, information on the chosen PDF set, value of the strong coupling, etc.) used for the first calculation achieved by the code, namely the LO rate associated with the process considered.This gives in our case LHAPDF 6.3.0 loading <lhapth>/CT14lo_0000.dat CT14lo PDF set, member #0, version 1; \ LHAPDF ID = 13200 alpha_s(M_Z^2) = 0.1180 for LO PDF set alpha_s(564.015^2)= 0.0927 for NLO PDF set where <lhapth> corresponds to the path to the folder containing the PDF grid CT14lo_0000.dat,that is used for the evaluation of the LO cross section.
Subsequently, Resummino displays to the screen the results of each VEGAS iteration, together with its numerical error, as well as the accumulated result obtained by combining all VEGAS iterations.The input file provided above leads to the screen output shown in figure 6, in which the first five iterations are related to the warmup stage of the numerical integration process, and the last five to the actual calculation.The latter stops when the required precision is reached.
Next, LHAPDF information on the NLO PDF set used is printed to the screen, followed by a similar output as in figure 6 for each component of the full fixedorder calculation.These consist of a recalculation of the LO rate with NLO parton densities, the sum of the virtual corrections and the integrated dipoles (in the language of the dipole subtraction formalism [91][92][93]), the   collinear remainder stemming from the P and K insertion operators (still in the language of the dipole subtraction formalism), and the difference between the real gluon emission corrections and the associated dipoles, and that between the real quark emission corrections and the associated dipoles.In the case in which on-shell resonant contributions appear in the real emission contributions, their integration is handled separately for improved convergence.
The code then moves with the calculation of the resummed components of the cross section, which starts with a fit of the PDFs so that they could be transformed to Mellin space (see section 2).Information on this fit is printed to the screen, as shown in figure 7 for the case of the gluon density.The cross section is then evaluated at NLO+NLL using the collinear-improved resummation formalism of [41][42][43][44], or at aNNLO+NNLL using the standard threshold resummation formalism of [35][36][37][38][39][40] if the code is run with its --nnll option (see Appendix E).
The matched results are finally displayed to the screen, 8Results: LO = (2.56583e-04+-2.73615e-07) pb final-state invariant mass for calculations at fixed invariant mass.
---mu_f (or -f) and --mu_r (or -r), both followed by a double-precision number: this determines the factors multiplying the factorisation and renormalisation scales, relatively to central scale choices.
---output (or -o), followed by a string referring to the path to a specific folder: this defines the folder in which all output files created by Resummino will be stored.
---parameter-log (or -p), followed by a string representing the path to a file: this makes the code writing all the parameters inherent to the calculation considered in the provided file.
---particle1 (or -c) and --particle2 (or -d), followed by integer numbers: modifications of the nature of the first and second final-state particles through a change in their PDG identifiers.
---pdfset_lo (or -a) and --pdfset_nlo (or -b), both followed by an integer: this determines the exact PDF set, within a given collection of LHAPDF parton densities, to employ for calculations at LO and beyond respectively.
---transverse-momentum (or -t), followed by a double-precision number: this sets the value, in GeV, of the final-state transverse momentum, for calculations at fixed p T .
---version (or -v): this displays the Resummino version number to the screen.

Appendix F: Lists of total cross sections
In this appendix we collect total cross section predictions for the different SUSY processes studied in this article, and present them in a tabulated form in which scale and PDF uncertainties are separate.Their digitised version can be found on the webpage https://github.com/APN-Pucky/HEPi/tree/master/hepi/data/json.
Table aNNLO+NNLL total cross sections for the production of a pair of compressed higgsino states at the LHC, for √ S = 13.0TeV (upper) and 13.6 TeV (lower).Higgsinos are defined as in eq. ( 12), and all other SUSY particles are decoupled.

13 Fig. 1
Fig.1Total cross sections for slepton pair production at the LHC, operating at centre-of-mass energies √ S = 13.6 TeV and 13 TeV (upper insets), shown together with their ratios to the 13 TeV total rates (lower insets) in which combined scale and

Fig. 3 2 χ+ 1
Fig. 3 Total cross sections for χ0 2 χ+ 1 production at the LHC, operating at centre-of-mass energies √ S = 13 TeV (left) and 13.6 TeV(right), in a scenario where all SUSY particles are decoupled with the exception of the non-degenerate produced states.Similar figures can be expected for other higgsino production modes (see table5).

13 Fig. 4 2 χ+ 1 (top) and χ0 1 χ− 1 (
Fig. 4 Total cross sections for wino pair production at the LHC, operating at centre-of-mass energies √ S = 13.6 TeV and 13 TeV (upper insets), shown together with their ratios to the 13 TeV total rates (lower insets) in which combined scale and PDF uncertainties are included.We consider the charged-current production of a χ0 2 χ+ 1 (top) and χ0 1 χ− 1 (centre) wino pair, as well as the neutral-current production of a χ+ 1 χ− 1 pair (bottom).Predictions are presented as a function of the wino mass m χ.

KFig. 5 2 χ+ 1
Fig. 5 Total cross sections for wino pair production at the LHC, operating at centre-of-mass energies √ S = 13 TeV (upper) and 13.6 TeV (lower).We consider the charged-current production of a χ0 2 χ+ 1 wino pair, and present predictions as a function of the wino mass m χ and the common first-generation and second-generation squark mass m q .

3 and the number of active quark flavours n f = 5 .
The individual process-independent resummation coefficients g (n) c

Fig. 6
Fig.6Resummino screen output relevant for the calculation of the LO total rate associated with slepton pair production at the LHC, using the SUSY spectrum encoded in the file input/slha.in and the collider settings provided in the file input/resummino.in,both files being shipped with the program.

Fig. 7
Fig. 7 Screen output relevant for the fit of the gluon density in the case of slepton pair production at the LHC operating at √ S = 13 TeV, the details on the calculation being encoded in the file input/resummino.inand the SUSY spectrum in the SLHA file input/slha.in.