Heavy Quarkonium Production Phenomenology and Automation of One-Loop Scattering Amplitude Computations