Add a helper function for the one-loop QCD corrections and make
all QCD functions return delta m / m.