code for qcd shape evaluation