add a function to read QCD yields from file
minor updates