infinities inherent in fixed-order calculations. Q is the characteristic scale of the interaction and is related to some physical scale in the interaction, such as the p$_{T}$ of the leading jet. This factorization of hadron-hadron interactions into a short distance interaction (hard interaction) and a large distance interaction (PDFs) is done at an arbitrary energy scale µ$_{F}$ . A parton emitted below the scale µ$_{F}$ is considered to be part of the hadron structure and thus is described by the PDFs. The hadron structure is measured by many different experiments, especially lepton-hadron scattering experiments. The PDFs are determined by fitting these data and parametrized at a scale µ$_{0}$ . The QCD evolution equations [7], currently available at NNLO in perturbation theory, are used to evolve these PDFs to any different scale µ . The cross section is a convolution of the PDFs and the parton level cross section and, to be consistent, both must be calculated to the same order in perturbation theory. A complete calculation, including all orders in the perturbation series, must be independent of µ$_{F}$ , but a fixed order calculation depends on this choice. To evaluate the sensitivity to µ$_{F}$ , it is usually varied up and down by a factor of two. Scale dependencies are expected to decline with the addition of higher order terms in the calculation. In Figure 2, for jet production, examples of leading order (LO) and next-to-leader order (NLO) Feynman diagrams are shown. At NLO, the jet cross section receives contributions from virtual corrections to the two parton final state, and from real corrections from the three parton final state. Both contributions are divergent, but the sum is finite. One sees that a NLO pQCD calculation can describe up to three jets in the final state. For inclusive jet production, many different implementations of NLO pQCD calculations are available [14–17]. All these programs use Monte Carlo integration techniques to calculate the real and virtual contribution to the cross sections. The EKS program [14] calculates the cross sections for user-predefined cuts on transverse momentum and rapidity ranges. JETRAD [15] and NLOJet++[16] generate parton events with weights (both positive and negative) and thus full event kinematics are available to the user for jet clustering, detector acceptance calculations and study of any other distribution. These programs require a huge amount of CPU to reach the desired precision. A large fraction of this CPU is consumed in evaluating the PDFs. In the FastNLO program [17], the convolution in Equation 1 is modified to a product of a perturbatively calculable piece ˜ σ$_{n,i,k,l,m}$ ( µ ), calculated using NLOJet++, the strong coupling constant α n $_{s}$( µ ( m $^{)}$), and a function F$_{i}$ ( x ( k ) a , x ( l ) b , µ ( m $^{)}$) which depends on only the PDFs and the factorization scale µ . The time consuming piece ˜ σ is 5