  We present a simple method to automatically evaluate arbitrary tree-levelamplitudes involving the production or decay of a heavy quark pair QQbar in ageneric {2S+1}L_J^[1,8] state, i.e., the short distance coefficients appearingin the NRQCD factorization formalism. Our approach is based on extracting therelevant contributions from the open heavy quark-antiquark amplitudes throughan expansion with respect to the quark-antiquark relative momentum and theapplication of suitable color and spin projectors. To illustrate thecapabilities of the method and its implementation in MadGraph a fewapplications to quarkonium collider phenomenology are presented.
Our approach is based on extracting the relevant contributions from the open heavy quark-antiquark amplitudes through an expansion with respect to the quark-antiquark relative momentum and the application of suitable color and spin projectors.