add code to find yields of QCD