Figure 4: Comparison of error bands obtained by the scale-variations for cross sections at NNLO in traditional (upper figure) and optimal (lower figure) TMD implementation. Here, the kinematics bin-integration, etc., is for the Z-boson production measure at ATLAS at 8 TeV [31].