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. arXiv:1907.04898v1 [hep-ph] 10 Jul 2019 7 Conclusions 47 A Conventions for one-point and two-point functions 49 B The MoGRe package 50 B.1 The main method MoGRe$Renormalize and its options 50 B.2 Simplifying the procedure 51 B.3 Restrictions 52 B.4 Specifying the renormalisation scheme 52 B.5 Clearing a renormalisation scheme 52 B.6 Technical details on the functioning of the method 53 C Decoupling of heavy SUSY particles in MG5 aMC 55 -1 -which allows for the cancellation of chiral anomalies and to generate masses for all up-type and down-type particles.