  The MadGraph5 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 FeynRules 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 Prospino2 and Resummino, and present illustrative applications
relevant to the 13 TeV LHC, both at the total-rate and differential levels. The
computer programs that we have used to obtain the predictions presented here
are all publicly available.
