Automated simulations beyond the Standard Model: supersymmetry

The Muc(ad)Guc(raph)5_aMC@NLO framework aims to automate all types of leading- and next-to-leading-order-accurate simulations for any user-defined model that stems from a renormalisable Lagrangian. In this paper, we present all of the key ingredients of such models in the context of supersymmetric theories. In order to do so, we extend the Fuc(eyn)Ruc(ules) package by giving it the possibility of dealing with different renormalisation options that are relevant to supersymmetric models. We also show how to deal with the problem posed by the presence of narrow resonances, thus generalising the so-called on-shell subtraction approaches. We extensively compare our total rate results with those of both Puc(rospino)2 and Ruc(esummino), and present illustrative applications relevant to the 13 TeV LHC, both at the total-rate and differential levels. The computer programmes that we have used to obtain the predictions presented here are all publicly available.