I have produced more NLO predictions, using HERAPDF2.0 and ATLAS-epWZ16-EIG PDF sets. As we discussed earlier at the meetings, we do not have proper RT OPT implementation for heavy quarks in CC in available public codes, so I believe it was treated in the ZMVFNS in all earlier HERAPDF-like fits, and I did the same now for the new predictions. Actually I checked that the differences when using ZMVFNS or FONLL-B are of the order of a couple of percent, so probably similarly small differences are expected when using RT OPT or ZMVFNS.

Attached you find the predictions calculated using the nominal HERAPDF2.0 set and four sets with fs variations, as well as using the nominal ATLAS-epWZ16-EIG set (in addition the first two predictions are listed, those are in FF and FONLL which Jae had already, they come with different partonic contributions). Please let me know whether you would like to include them in the paper and will need any other details? (one would need to cite QCDNUM - the package used to calculate predictions in the ZMVFNS)

Regarding the scale variation uncertainties: those could be done in the FF or FONLL schemes only, there is no such implementation in the QCDNUM code. I am not sure whether it will be fine to display them for some sets of the predictions only?



scheme: FFABM, PDF: ABMP16
e-p 200 < Q^2 < 1500 GeV^2  : 4.50[100%] = 0.19[  4%] + 2.31[ 51%] + 2.01[ 45%] + 0.00[  0%]
e-p 1500 < Q^2 < 60000 GeV^2: 1.73[100%] = 0.08[  5%] + 0.85[ 49%] + 0.80[ 46%] + 0.00[  0%]
e+p 200 < Q^2 < 1500 GeV^2  : 4.72[100%] = 0.40[  8%] + 2.31[ 49%] + 2.01[ 43%] + 0.00[  0%]
e+p 1500 < Q^2 < 60000 GeV^2: 1.97[100%] = 0.32[ 16%] + 0.85[ 43%] + 0.80[ 41%] + 0.00[  0%]

scheme: FONLL, PDF: NNPDF3.1
e-p 200 < Q^2 < 1500 GeV^2  : 4.98[100%] = 0.19[  4%] + 2.12[ 43%] + -0.50[-10%] + 3.16[ 64%]
e-p 1500 < Q^2 < 60000 GeV^2: 2.16[100%] = 0.09[  4%] + 0.71[ 33%] + -0.26[-12%] + 1.62[ 75%]
e+p 200 < Q^2 < 1500 GeV^2  : 5.37[100%] = 0.42[  8%] + 2.29[ 43%] + -0.50[ -9%] + 3.16[ 59%]
e+p 1500 < Q^2 < 60000 GeV^2: 2.66[100%] = 0.32[ 12%] + 0.97[ 37%] + -0.26[-10%] + 1.62[ 61%]

scheme: ZMVFNS, PDF: HERAPDF2.0 fs=0.4 [nominal]
e-p 200 < Q^2 < 1500 GeV^2  : 5.41
e-p 1500 < Q^2 < 60000 GeV^2: 2.30
e+p 200 < Q^2 < 1500 GeV^2  : 5.67
e+p 1500 < Q^2 < 60000 GeV^2: 2.57

scheme: ZMVFNS, PDF: HERAPDF2.0 fs=0.3
e-p 200 < Q^2 < 1500 GeV^2  : 5.15
e-p 1500 < Q^2 < 60000 GeV^2: 2.21
e+p 200 < Q^2 < 1500 GeV^2  : 5.40
e+p 1500 < Q^2 < 60000 GeV^2: 2.47

scheme: ZMVFNS, PDF: HERAPDF2.0 fs=0.5
e-p 200 < Q^2 < 1500 GeV^2  : 5.70
e-p 1500 < Q^2 < 60000 GeV^2: 2.37
e+p 200 < Q^2 < 1500 GeV^2  : 5.96
e+p 1500 < Q^2 < 60000 GeV^2: 2.65

scheme: ZMVFNS, PDF: HERAPDF2.0 fs=hermesfs-03
e-p 200 < Q^2 < 1500 GeV^2  : 4.79
e-p 1500 < Q^2 < 60000 GeV^2: 1.89
e+p 200 < Q^2 < 1500 GeV^2  : 5.05
e+p 1500 < Q^2 < 60000 GeV^2: 2.16

scheme: ZMVFNS, PDF: HERAPDF2.0 fs=hermesfs-05
e-p 200 < Q^2 < 1500 GeV^2  : 5.12
e-p 1500 < Q^2 < 60000 GeV^2: 1.93
e+p 200 < Q^2 < 1500 GeV^2  : 5.38
e+p 1500 < Q^2 < 60000 GeV^2: 2.20

scheme: ZMVFNS, PDF: ATLAS-epWZ16-EIG
e-p 200 < Q^2 < 1500 GeV^2  : 6.14
e-p 1500 < Q^2 < 60000 GeV^2: 2.78
e+p 200 < Q^2 < 1500 GeV^2  : 6.41
e+p 1500 < Q^2 < 60000 GeV^2: 3.07