diff --git a/test/.gitignore b/test/.gitignore
new file mode 100644
index 0000000000000000000000000000000000000000..5b8cd369cc6ea006e0abca96812c1da53e007b4c
--- /dev/null
+++ b/test/.gitignore
@@ -0,0 +1,7 @@
+.DS_Store
+.idea
+*.log
+tmp/
+
+output
+InclusiveNJets.str
diff --git a/test/InclusiveNJets_InclusiveJets_fix.wrm b/test/InclusiveNJets_InclusiveJets_fix.wrm
new file mode 100644
index 0000000000000000000000000000000000000000..5d3bc43e838a590b6ef2c2e0383ad6c85ecc4b58
--- /dev/null
+++ b/test/InclusiveNJets_InclusiveJets_fix.wrm
@@ -0,0 +1,386 @@
+# --- Use emacs in sh mode -*-sh-*- #
+# This is a automatically generated file by fastNLO and holds the values of the warmup run.
+# The values are valid for the scenario InclusiveJets_Example_v23_fix_ref
+# and if calculated with the steerfile: InclusiveNJets
+# but only if no serious changes have been performed since its creation.
+#
+# Delete this file, if you want fastNLO to calculate a new one.
+#
+# This file has been calculated using:
+# 1e+09 contributions.
+# 12962320896 entries.
+# ( Mind: contributions != events. And contributions are not necessarily in phase space region.
+# Please check by eye for reasonability of the values.
+# Number of events per bin are listed below.
+
+Warmup.OrderInAlphasOfWarmupRunWas 2
+Warmup.CheckScaleLimitsAgainstBins true
+Warmup.ScaleDescriptionScale1 "pT_jet_[GeV]"
+Warmup.DifferentialDimension 2
+Warmup.DimensionLabels {
+ "|y|" "pT_[GeV]"
+}
+Warmup.DimensionIsDifferential {
+ "2" "2"
+}
+
+Warmup.Values {{
+ ObsBin x_min x_max pT_jet_[GeV]_min pT_jet_[GeV]_max
+ 0 1.0e-09 1.00e+00 18.000000 21.000000
+ 1 1.0e-09 1.00e+00 21.000000 24.000000
+ 2 1.0e-09 1.00e+00 24.000000 28.000000
+ 3 1.0e-09 1.00e+00 28.000000 32.000000
+ 4 1.0e-09 1.00e+00 32.000000 37.000000
+ 5 1.0e-09 1.00e+00 37.000000 43.000000
+ 6 1.0e-09 1.00e+00 43.000000 49.000000
+ 7 1.0e-09 1.00e+00 49.000000 56.000000
+ 8 1.0e-09 1.00e+00 56.000000 64.000000
+ 9 1.0e-09 1.00e+00 64.000000 74.000000
+ 10 1.0e-09 1.00e+00 74.000000 84.000000
+ 11 1.0e-09 1.00e+00 84.000000 97.000000
+ 12 1.0e-09 1.00e+00 97.000000 114.000000
+ 13 1.0e-09 1.00e+00 114.000000 133.000000
+ 14 1.0e-09 1.00e+00 133.000000 153.000000
+ 15 1.0e-09 1.00e+00 153.000000 174.000000
+ 16 1.0e-09 1.00e+00 174.000000 196.000000
+ 17 1.0e-09 1.00e+00 196.000000 220.000000
+ 18 1.0e-09 1.00e+00 220.000000 245.000000
+ 19 1.0e-09 1.00e+00 245.000000 272.000000
+ 20 1.0e-09 1.00e+00 272.000000 300.000000
+ 21 1.0e-09 1.00e+00 300.000000 330.000000
+ 22 1.0e-09 1.00e+00 330.000000 362.000000
+ 23 1.0e-09 1.00e+00 362.000000 395.000000
+ 24 1.0e-09 1.00e+00 395.000000 430.000000
+ 25 1.0e-09 1.00e+00 430.000000 468.000000
+ 26 1.0e-09 1.00e+00 468.000000 507.000000
+ 27 1.0e-09 1.00e+00 507.000000 548.000000
+ 28 1.0e-09 1.00e+00 548.000000 592.000000
+ 29 1.0e-09 1.00e+00 592.000000 638.000000
+ 30 1.0e-09 1.00e+00 638.000000 686.000000
+ 31 1.0e-09 1.00e+00 686.000000 737.000000
+ 32 1.0e-09 1.00e+00 737.000000 846.000000
+ 33 1.0e-09 1.00e+00 846.000000 1684.000000
+ 34 1.0e-09 1.00e+00 18.000000 21.000000
+ 35 1.0e-09 1.00e+00 21.000000 24.000000
+ 36 1.0e-09 1.00e+00 24.000000 28.000000
+ 37 1.0e-09 1.00e+00 28.000000 32.000000
+ 38 1.0e-09 1.00e+00 32.000000 37.000000
+ 39 1.0e-09 1.00e+00 37.000000 43.000000
+ 40 1.0e-09 1.00e+00 43.000000 49.000000
+ 41 1.0e-09 1.00e+00 49.000000 56.000000
+ 42 1.0e-09 1.00e+00 56.000000 64.000000
+ 43 1.0e-09 1.00e+00 64.000000 74.000000
+ 44 1.0e-09 1.00e+00 74.000000 84.000000
+ 45 1.0e-09 1.00e+00 84.000000 97.000000
+ 46 1.0e-09 1.00e+00 97.000000 114.000000
+ 47 1.0e-09 1.00e+00 114.000000 133.000000
+ 48 1.0e-09 1.00e+00 133.000000 153.000000
+ 49 1.0e-09 1.00e+00 153.000000 174.000000
+ 50 1.0e-09 1.00e+00 174.000000 196.000000
+ 51 1.0e-09 1.00e+00 196.000000 220.000000
+ 52 1.0e-09 1.00e+00 220.000000 245.000000
+ 53 1.0e-09 1.00e+00 245.000000 272.000000
+ 54 1.0e-09 1.00e+00 272.000000 300.000000
+ 55 1.0e-09 1.00e+00 300.000000 330.000000
+ 56 1.0e-09 1.00e+00 330.000000 362.000000
+ 57 1.0e-09 1.00e+00 362.000000 395.000000
+ 58 1.0e-09 1.00e+00 395.000000 430.000000
+ 59 1.0e-09 1.00e+00 430.000000 468.000000
+ 60 1.0e-09 1.00e+00 468.000000 507.000000
+ 61 1.0e-09 1.00e+00 507.000000 548.000000
+ 62 1.0e-09 1.00e+00 548.000000 592.000000
+ 63 1.0e-09 1.00e+00 592.000000 638.000000
+ 64 1.0e-09 1.00e+00 638.000000 686.000000
+ 65 1.0e-09 1.00e+00 686.000000 790.000000
+ 66 1.0e-09 1.00e+00 790.000000 1684.000000
+ 67 1.0e-09 1.00e+00 18.000000 21.000000
+ 68 1.0e-09 1.00e+00 21.000000 24.000000
+ 69 1.0e-09 1.00e+00 24.000000 28.000000
+ 70 1.0e-09 1.00e+00 28.000000 32.000000
+ 71 1.0e-09 1.00e+00 32.000000 37.000000
+ 72 1.0e-09 1.00e+00 37.000000 43.000000
+ 73 1.0e-09 1.00e+00 43.000000 49.000000
+ 74 1.0e-09 1.00e+00 49.000000 56.000000
+ 75 1.0e-09 1.00e+00 56.000000 64.000000
+ 76 1.0e-09 1.00e+00 64.000000 74.000000
+ 77 1.0e-09 1.00e+00 74.000000 84.000000
+ 78 1.0e-09 1.00e+00 84.000000 97.000000
+ 79 1.0e-09 1.00e+00 97.000000 114.000000
+ 80 1.0e-09 1.00e+00 114.000000 133.000000
+ 81 1.0e-09 1.00e+00 133.000000 153.000000
+ 82 1.0e-09 1.00e+00 153.000000 174.000000
+ 83 1.0e-09 1.00e+00 174.000000 196.000000
+ 84 1.0e-09 1.00e+00 196.000000 220.000000
+ 85 1.0e-09 1.00e+00 220.000000 245.000000
+ 86 1.0e-09 1.00e+00 245.000000 272.000000
+ 87 1.0e-09 1.00e+00 272.000000 300.000000
+ 88 1.0e-09 1.00e+00 300.000000 330.000000
+ 89 1.0e-09 1.00e+00 330.000000 362.000000
+ 90 1.0e-09 1.00e+00 362.000000 395.000000
+ 91 1.0e-09 1.00e+00 395.000000 430.000000
+ 92 1.0e-09 1.00e+00 430.000000 468.000000
+ 93 1.0e-09 1.00e+00 468.000000 507.000000
+ 94 1.0e-09 1.00e+00 507.000000 548.000000
+ 95 1.0e-09 1.00e+00 548.000000 592.000000
+ 96 1.0e-09 1.00e+00 592.000000 638.000000
+ 97 1.0e-09 1.00e+00 638.000000 686.000000
+ 98 1.0e-09 1.00e+00 686.000000 1410.000000
+ 99 1.0e-09 1.00e+00 18.000000 21.000000
+ 100 1.0e-09 1.00e+00 21.000000 24.000000
+ 101 1.0e-09 1.00e+00 24.000000 28.000000
+ 102 1.0e-09 1.00e+00 28.000000 32.000000
+ 103 1.0e-09 1.00e+00 32.000000 37.000000
+ 104 1.0e-09 1.00e+00 37.000000 43.000000
+ 105 1.0e-09 1.00e+00 43.000000 49.000000
+ 106 1.0e-09 1.00e+00 49.000000 56.000000
+ 107 1.0e-09 1.00e+00 56.000000 64.000000
+ 108 1.0e-09 1.00e+00 64.000000 74.000000
+ 109 1.0e-09 1.00e+00 74.000000 84.000000
+ 110 1.0e-09 1.00e+00 84.000000 97.000000
+ 111 1.0e-09 1.00e+00 97.000000 114.000000
+ 112 1.0e-09 1.00e+00 114.000000 133.000000
+ 113 1.0e-09 1.00e+00 133.000000 153.000000
+ 114 1.0e-09 1.00e+00 153.000000 174.000000
+ 115 1.0e-09 1.00e+00 174.000000 196.000000
+ 116 1.0e-09 1.00e+00 196.000000 220.000000
+ 117 1.0e-09 1.00e+00 220.000000 245.000000
+ 118 1.0e-09 1.00e+00 245.000000 272.000000
+ 119 1.0e-09 1.00e+00 272.000000 300.000000
+ 120 1.0e-09 1.00e+00 300.000000 330.000000
+ 121 1.0e-09 1.00e+00 330.000000 362.000000
+ 122 1.0e-09 1.00e+00 362.000000 395.000000
+ 123 1.0e-09 1.00e+00 395.000000 430.000000
+ 124 1.0e-09 1.00e+00 430.000000 468.000000
+ 125 1.0e-09 1.00e+00 468.000000 507.000000
+ 126 1.0e-09 1.00e+00 507.000000 548.000000
+ 127 1.0e-09 1.00e+00 548.000000 1032.000000
+ 128 1.0e-09 1.00e+00 18.000000 21.000000
+ 129 1.0e-09 1.00e+00 21.000000 24.000000
+ 130 1.0e-09 1.00e+00 24.000000 28.000000
+ 131 1.0e-09 1.00e+00 28.000000 32.000000
+ 132 1.0e-09 1.00e+00 32.000000 37.000000
+ 133 1.0e-09 1.00e+00 37.000000 43.000000
+ 134 1.0e-09 1.00e+00 43.000000 49.000000
+ 135 1.0e-09 1.00e+00 49.000000 56.000000
+ 136 1.0e-09 1.00e+00 56.000000 64.000000
+ 137 1.0e-09 1.00e+00 64.000000 74.000000
+ 138 1.0e-09 1.00e+00 74.000000 84.000000
+ 139 1.0e-09 1.00e+00 84.000000 97.000000
+ 140 1.0e-09 1.00e+00 97.000000 114.000000
+ 141 1.0e-09 1.00e+00 114.000000 133.000000
+ 142 1.0e-09 1.00e+00 133.000000 153.000000
+ 143 1.0e-09 1.00e+00 153.000000 174.000000
+ 144 1.0e-09 1.00e+00 174.000000 196.000000
+ 145 1.0e-09 1.00e+00 196.000000 220.000000
+ 146 1.0e-09 1.00e+00 220.000000 245.000000
+ 147 1.0e-09 1.00e+00 245.000000 272.000000
+ 148 1.0e-09 1.00e+00 272.000000 300.000000
+ 149 1.0e-09 1.00e+00 300.000000 330.000000
+ 150 1.0e-09 1.00e+00 330.000000 362.000000
+ 151 1.0e-09 1.00e+00 362.000000 395.000000
+ 152 1.0e-09 1.00e+00 395.000000 430.000000
+ 153 1.0e-09 1.00e+00 430.000000 737.000000
+ 154 1.0e-09 1.00e+00 18.000000 21.000000
+ 155 1.0e-09 1.00e+00 21.000000 24.000000
+ 156 1.0e-09 1.00e+00 24.000000 28.000000
+ 157 1.0e-09 1.00e+00 28.000000 32.000000
+ 158 1.0e-09 1.00e+00 32.000000 37.000000
+ 159 1.0e-09 1.00e+00 37.000000 43.000000
+ 160 1.0e-09 1.00e+00 43.000000 49.000000
+ 161 1.0e-09 1.00e+00 49.000000 56.000000
+ 162 1.0e-09 1.00e+00 56.000000 64.000000
+ 163 1.0e-09 1.00e+00 64.000000 74.000000
+ 164 1.0e-09 1.00e+00 74.000000 84.000000
+ 165 1.0e-09 1.00e+00 84.000000 97.000000
+ 166 1.0e-09 1.00e+00 97.000000 114.000000
+ 167 1.0e-09 1.00e+00 114.000000 133.000000
+ 168 1.0e-09 1.00e+00 133.000000 153.000000
+ 169 1.0e-09 1.00e+00 153.000000 174.000000
+ 170 1.0e-09 1.00e+00 174.000000 196.000000
+ 171 1.0e-09 1.00e+00 196.000000 220.000000
+ 172 1.0e-09 1.00e+00 220.000000 245.000000
+ 173 1.0e-09 1.00e+00 245.000000 272.000000
+ 174 1.0e-09 1.00e+00 272.000000 300.000000
+ 175 1.0e-09 1.00e+00 300.000000 468.000000
+}}
+
+
+Warmup.Binning {{
+ ObsBin |y|_Lo |y|_Up pT_[GeV]_Lo pT_[GeV]_Up BinSize EventCount
+ 0 0.00000 0.500000 18.0000 21.0000 3.00000 3244885
+ 1 0.00000 0.500000 21.0000 24.0000 3.00000 3509086
+ 2 0.00000 0.500000 24.0000 28.0000 4.00000 5089791
+ 3 0.00000 0.500000 28.0000 32.0000 4.00000 5578139
+ 4 0.00000 0.500000 32.0000 37.0000 5.00000 13849486
+ 5 0.00000 0.500000 37.0000 43.0000 6.00000 42220087
+ 6 0.00000 0.500000 43.0000 49.0000 6.00000 63797111
+ 7 0.00000 0.500000 49.0000 56.0000 7.00000 89607455
+ 8 0.00000 0.500000 56.0000 64.0000 8.00000 111378666
+ 9 0.00000 0.500000 64.0000 74.0000 10.0000 142213393
+ 10 0.00000 0.500000 74.0000 84.0000 10.0000 138600203
+ 11 0.00000 0.500000 84.0000 97.0000 13.0000 168761621
+ 12 0.00000 0.500000 97.0000 114.000 17.0000 196475552
+ 13 0.00000 0.500000 114.000 133.000 19.0000 187801775
+ 14 0.00000 0.500000 133.000 153.000 20.0000 169574713
+ 15 0.00000 0.500000 153.000 174.000 21.0000 154746214
+ 16 0.00000 0.500000 174.000 196.000 22.0000 142317245
+ 17 0.00000 0.500000 196.000 220.000 24.0000 137085599
+ 18 0.00000 0.500000 220.000 245.000 25.0000 126718578
+ 19 0.00000 0.500000 245.000 272.000 27.0000 122031084
+ 20 0.00000 0.500000 272.000 300.000 28.0000 113359687
+ 21 0.00000 0.500000 300.000 330.000 30.0000 109234608
+ 22 0.00000 0.500000 330.000 362.000 32.0000 104985678
+ 23 0.00000 0.500000 362.000 395.000 33.0000 97925128
+ 24 0.00000 0.500000 395.000 430.000 35.0000 94194296
+ 25 0.00000 0.500000 430.000 468.000 38.0000 92800792
+ 26 0.00000 0.500000 468.000 507.000 39.0000 86523591
+ 27 0.00000 0.500000 507.000 548.000 41.0000 82887126
+ 28 0.00000 0.500000 548.000 592.000 44.0000 81112535
+ 29 0.00000 0.500000 592.000 638.000 46.0000 77302834
+ 30 0.00000 0.500000 638.000 686.000 48.0000 73685892
+ 31 0.00000 0.500000 686.000 737.000 51.0000 71466808
+ 32 0.00000 0.500000 737.000 846.000 109.000 133404502
+ 33 0.00000 0.500000 846.000 1684.00 838.000 555151128
+ 34 0.500000 1.00000 18.0000 21.0000 3.00000 3652418
+ 35 0.500000 1.00000 21.0000 24.0000 3.00000 4121614
+ 36 0.500000 1.00000 24.0000 28.0000 4.00000 10869243
+ 37 0.500000 1.00000 28.0000 32.0000 4.00000 21632310
+ 38 0.500000 1.00000 32.0000 37.0000 5.00000 39926467
+ 39 0.500000 1.00000 37.0000 43.0000 6.00000 54958708
+ 40 0.500000 1.00000 43.0000 49.0000 6.00000 55178704
+ 41 0.500000 1.00000 49.0000 56.0000 7.00000 63711732
+ 42 0.500000 1.00000 56.0000 64.0000 8.00000 72031890
+ 43 0.500000 1.00000 64.0000 74.0000 10.0000 89189226
+ 44 0.500000 1.00000 74.0000 84.0000 10.0000 88320155
+ 45 0.500000 1.00000 84.0000 97.0000 13.0000 113800743
+ 46 0.500000 1.00000 97.0000 114.000 17.0000 147458304
+ 47 0.500000 1.00000 114.000 133.000 19.0000 161287756
+ 48 0.500000 1.00000 133.000 153.000 20.0000 159145063
+ 49 0.500000 1.00000 153.000 174.000 21.0000 151362484
+ 50 0.500000 1.00000 174.000 196.000 22.0000 140190638
+ 51 0.500000 1.00000 196.000 220.000 24.0000 134795983
+ 52 0.500000 1.00000 220.000 245.000 25.0000 124324263
+ 53 0.500000 1.00000 245.000 272.000 27.0000 119533099
+ 54 0.500000 1.00000 272.000 300.000 28.0000 110665919
+ 55 0.500000 1.00000 300.000 330.000 30.0000 106354192
+ 56 0.500000 1.00000 330.000 362.000 32.0000 101883684
+ 57 0.500000 1.00000 362.000 395.000 33.0000 94775793
+ 58 0.500000 1.00000 395.000 430.000 35.0000 90787466
+ 59 0.500000 1.00000 430.000 468.000 38.0000 89137881
+ 60 0.500000 1.00000 468.000 507.000 39.0000 82816426
+ 61 0.500000 1.00000 507.000 548.000 41.0000 78978655
+ 62 0.500000 1.00000 548.000 592.000 44.0000 76859111
+ 63 0.500000 1.00000 592.000 638.000 46.0000 72925993
+ 64 0.500000 1.00000 638.000 686.000 48.0000 69029065
+ 65 0.500000 1.00000 686.000 790.000 104.000 129318756
+ 66 0.500000 1.00000 790.000 1684.00 894.000 534602418
+ 67 1.00000 1.50000 18.0000 21.0000 3.00000 11359250
+ 68 1.00000 1.50000 21.0000 24.0000 3.00000 16241652
+ 69 1.00000 1.50000 24.0000 28.0000 4.00000 24808483
+ 70 1.00000 1.50000 28.0000 32.0000 4.00000 24478552
+ 71 1.00000 1.50000 32.0000 37.0000 5.00000 30146578
+ 72 1.00000 1.50000 37.0000 43.0000 6.00000 35559468
+ 73 1.00000 1.50000 43.0000 49.0000 6.00000 35040411
+ 74 1.00000 1.50000 49.0000 56.0000 7.00000 40339229
+ 75 1.00000 1.50000 56.0000 64.0000 8.00000 45514427
+ 76 1.00000 1.50000 64.0000 74.0000 10.0000 56121821
+ 77 1.00000 1.50000 74.0000 84.0000 10.0000 55474552
+ 78 1.00000 1.50000 84.0000 97.0000 13.0000 71379882
+ 79 1.00000 1.50000 97.0000 114.000 17.0000 92203391
+ 80 1.00000 1.50000 114.000 133.000 19.0000 101887345
+ 81 1.00000 1.50000 133.000 153.000 20.0000 106140790
+ 82 1.00000 1.50000 153.000 174.000 21.0000 110465635
+ 83 1.00000 1.50000 174.000 196.000 22.0000 114440872
+ 84 1.00000 1.50000 196.000 220.000 24.0000 120024289
+ 85 1.00000 1.50000 220.000 245.000 25.0000 116193861
+ 86 1.00000 1.50000 245.000 272.000 27.0000 113484315
+ 87 1.00000 1.50000 272.000 300.000 28.0000 104595897
+ 88 1.00000 1.50000 300.000 330.000 30.0000 99818012
+ 89 1.00000 1.50000 330.000 362.000 32.0000 94984358
+ 90 1.00000 1.50000 362.000 395.000 33.0000 87580339
+ 91 1.00000 1.50000 395.000 430.000 35.0000 83217708
+ 92 1.00000 1.50000 430.000 468.000 38.0000 80954958
+ 93 1.00000 1.50000 468.000 507.000 39.0000 74328961
+ 94 1.00000 1.50000 507.000 548.000 41.0000 70068754
+ 95 1.00000 1.50000 548.000 592.000 44.0000 67359915
+ 96 1.00000 1.50000 592.000 638.000 46.0000 62922713
+ 97 1.00000 1.50000 638.000 686.000 48.0000 58677829
+ 98 1.00000 1.50000 686.000 1410.00 724.000 412025222
+ 99 1.50000 2.00000 18.0000 21.0000 3.00000 12728744
+ 100 1.50000 2.00000 21.0000 24.0000 3.00000 12423418
+ 101 1.50000 2.00000 24.0000 28.0000 4.00000 16214380
+ 102 1.50000 2.00000 28.0000 32.0000 4.00000 15861020
+ 103 1.50000 2.00000 32.0000 37.0000 5.00000 19408144
+ 104 1.50000 2.00000 37.0000 43.0000 6.00000 22853397
+ 105 1.50000 2.00000 43.0000 49.0000 6.00000 22425802
+ 106 1.50000 2.00000 49.0000 56.0000 7.00000 25707941
+ 107 1.50000 2.00000 56.0000 64.0000 8.00000 28978887
+ 108 1.50000 2.00000 64.0000 74.0000 10.0000 35647619
+ 109 1.50000 2.00000 74.0000 84.0000 10.0000 35133084
+ 110 1.50000 2.00000 84.0000 97.0000 13.0000 45000130
+ 111 1.50000 2.00000 97.0000 114.000 17.0000 57979131
+ 112 1.50000 2.00000 114.000 133.000 19.0000 63884730
+ 113 1.50000 2.00000 133.000 153.000 20.0000 66388812
+ 114 1.50000 2.00000 153.000 174.000 21.0000 68976530
+ 115 1.50000 2.00000 174.000 196.000 22.0000 71512518
+ 116 1.50000 2.00000 196.000 220.000 24.0000 77389235
+ 117 1.50000 2.00000 220.000 245.000 25.0000 79905868
+ 118 1.50000 2.00000 245.000 272.000 27.0000 85600669
+ 119 1.50000 2.00000 272.000 300.000 28.0000 86526419
+ 120 1.50000 2.00000 300.000 330.000 30.0000 86624748
+ 121 1.50000 2.00000 330.000 362.000 32.0000 82373851
+ 122 1.50000 2.00000 362.000 395.000 33.0000 74632103
+ 123 1.50000 2.00000 395.000 430.000 35.0000 69485661
+ 124 1.50000 2.00000 430.000 468.000 38.0000 66043047
+ 125 1.50000 2.00000 468.000 507.000 39.0000 59062346
+ 126 1.50000 2.00000 507.000 548.000 41.0000 53929043
+ 127 1.50000 2.00000 548.000 1032.00 484.000 277256812
+ 128 2.00000 2.50000 18.0000 21.0000 3.00000 8446858
+ 129 2.00000 2.50000 21.0000 24.0000 3.00000 8203545
+ 130 2.00000 2.50000 24.0000 28.0000 4.00000 10645292
+ 131 2.00000 2.50000 28.0000 32.0000 4.00000 10380489
+ 132 2.00000 2.50000 32.0000 37.0000 5.00000 12652584
+ 133 2.00000 2.50000 37.0000 43.0000 6.00000 14817859
+ 134 2.00000 2.50000 43.0000 49.0000 6.00000 14486416
+ 135 2.00000 2.50000 49.0000 56.0000 7.00000 16570393
+ 136 2.00000 2.50000 56.0000 64.0000 8.00000 18575046
+ 137 2.00000 2.50000 64.0000 74.0000 10.0000 22788563
+ 138 2.00000 2.50000 74.0000 84.0000 10.0000 22357286
+ 139 2.00000 2.50000 84.0000 97.0000 13.0000 28559405
+ 140 2.00000 2.50000 97.0000 114.000 17.0000 36634892
+ 141 2.00000 2.50000 114.000 133.000 19.0000 40209204
+ 142 2.00000 2.50000 133.000 153.000 20.0000 41636119
+ 143 2.00000 2.50000 153.000 174.000 21.0000 43092427
+ 144 2.00000 2.50000 174.000 196.000 22.0000 44605666
+ 145 2.00000 2.50000 196.000 220.000 24.0000 48055329
+ 146 2.00000 2.50000 220.000 245.000 25.0000 49464863
+ 147 2.00000 2.50000 245.000 272.000 27.0000 52872778
+ 148 2.00000 2.50000 272.000 300.000 28.0000 54323871
+ 149 2.00000 2.50000 300.000 330.000 30.0000 57670655
+ 150 2.00000 2.50000 330.000 362.000 32.0000 59613575
+ 151 2.00000 2.50000 362.000 395.000 33.0000 52548181
+ 152 2.00000 2.50000 395.000 430.000 35.0000 46075141
+ 153 2.00000 2.50000 430.000 737.000 307.000 151430405
+ 154 2.50000 3.00000 18.0000 21.0000 3.00000 5656238
+ 155 2.50000 3.00000 21.0000 24.0000 3.00000 5485445
+ 156 2.50000 3.00000 24.0000 28.0000 4.00000 7081781
+ 157 2.50000 3.00000 28.0000 32.0000 4.00000 6863745
+ 158 2.50000 3.00000 32.0000 37.0000 5.00000 8337637
+ 159 2.50000 3.00000 37.0000 43.0000 6.00000 9712710
+ 160 2.50000 3.00000 43.0000 49.0000 6.00000 9449307
+ 161 2.50000 3.00000 49.0000 56.0000 7.00000 10749788
+ 162 2.50000 3.00000 56.0000 64.0000 8.00000 11987227
+ 163 2.50000 3.00000 64.0000 74.0000 10.0000 14638834
+ 164 2.50000 3.00000 74.0000 84.0000 10.0000 14286650
+ 165 2.50000 3.00000 84.0000 97.0000 13.0000 18191222
+ 166 2.50000 3.00000 97.0000 114.000 17.0000 23206827
+ 167 2.50000 3.00000 114.000 133.000 19.0000 25287087
+ 168 2.50000 3.00000 133.000 153.000 20.0000 26054749
+ 169 2.50000 3.00000 153.000 174.000 21.0000 26825953
+ 170 2.50000 3.00000 174.000 196.000 22.0000 27583745
+ 171 2.50000 3.00000 196.000 220.000 24.0000 29590589
+ 172 2.50000 3.00000 220.000 245.000 25.0000 30301929
+ 173 2.50000 3.00000 245.000 272.000 27.0000 32195240
+ 174 2.50000 3.00000 272.000 300.000 28.0000 32609430
+ 175 2.50000 3.00000 300.000 468.000 168.000 89452916
+}}
diff --git a/test/InclusiveNJets_InclusiveJets_flex.wrm b/test/InclusiveNJets_InclusiveJets_flex.wrm
new file mode 100644
index 0000000000000000000000000000000000000000..7c489b9b4009e8fca7d5d9d54c6e97a89f0a26d5
--- /dev/null
+++ b/test/InclusiveNJets_InclusiveJets_flex.wrm
@@ -0,0 +1,387 @@
+# --- Use emacs in sh mode -*-sh-*- #
+# This is a automatically generated file by fastNLO and holds the values of the warmup run.
+# The values are valid for the scenario InclusiveJets_flex
+# and if calculated with the steerfile: InclusiveNJets
+# but only if no serious changes have been performed since its creation.
+#
+# Delete this file, if you want fastNLO to calculate a new one.
+#
+# This file has been calculated using:
+# 100000 contributions.
+# 1296561 entries.
+# ( Mind: contributions != events. And contributions are not necessarily in phase space region.
+# Please check by eye for reasonability of the values.
+# Number of events per bin are listed below.
+
+Warmup.OrderInAlphasOfWarmupRunWas 2
+Warmup.CheckScaleLimitsAgainstBins true
+Warmup.ScaleDescriptionScale1 "pT_jet_[GeV]"
+Warmup.ScaleDescriptionScale2 "pT_max_[GeV]"
+Warmup.DifferentialDimension 2
+Warmup.DimensionLabels {
+ "|y|" "pT_[GeV]"
+}
+Warmup.DimensionIsDifferential {
+ "2" "2"
+}
+
+Warmup.Values {{
+ ObsBin x_min x_max pT_jet_[GeV]_min pT_jet_[GeV]_max pT_max_[GeV]_min pT_max_[GeV]_max
+ 0 1.0e-09 1.00e+00 18.000000 21.000000 18.000000 21.000000
+ 1 1.0e-09 1.00e+00 21.000000 24.000000 21.000000 24.000000
+ 2 1.0e-09 1.00e+00 24.000000 28.000000 24.000000 28.000000
+ 3 1.0e-09 1.00e+00 28.000000 32.000000 28.000000 32.000000
+ 4 1.0e-09 1.00e+00 32.000000 37.000000 32.000000 37.000000
+ 5 1.0e-09 1.00e+00 37.000000 43.000000 37.000000 43.000000
+ 6 1.0e-09 1.00e+00 43.000000 49.000000 43.000000 49.000000
+ 7 1.0e-09 1.00e+00 49.000000 56.000000 49.000000 56.000000
+ 8 1.0e-09 1.00e+00 56.000000 64.000000 56.000000 64.000000
+ 9 1.0e-09 1.00e+00 64.000000 74.000000 64.000000 74.000000
+ 10 1.0e-09 1.00e+00 74.000000 84.000000 74.000000 84.000000
+ 11 1.0e-09 1.00e+00 84.000000 97.000000 84.000000 97.000000
+ 12 1.0e-09 1.00e+00 97.000000 114.000000 97.000000 114.000000
+ 13 1.0e-09 1.00e+00 114.000000 133.000000 114.000000 133.000000
+ 14 1.0e-09 1.00e+00 133.000000 153.000000 133.000000 153.000000
+ 15 1.0e-09 1.00e+00 153.000000 174.000000 153.000000 174.000000
+ 16 1.0e-09 1.00e+00 174.000000 196.000000 174.000000 196.000000
+ 17 1.0e-09 1.00e+00 196.000000 220.000000 196.000000 220.000000
+ 18 1.0e-09 1.00e+00 220.000000 245.000000 220.000000 245.000000
+ 19 1.0e-09 1.00e+00 245.000000 272.000000 245.000000 272.000000
+ 20 1.0e-09 1.00e+00 272.000000 300.000000 272.000000 300.000000
+ 21 1.0e-09 1.00e+00 300.000000 330.000000 300.000000 330.000000
+ 22 1.0e-09 1.00e+00 330.000000 362.000000 330.000000 362.000000
+ 23 1.0e-09 1.00e+00 362.000000 395.000000 362.000000 395.000000
+ 24 1.0e-09 1.00e+00 395.000000 430.000000 395.000000 430.000000
+ 25 1.0e-09 1.00e+00 430.000000 468.000000 430.000000 468.000000
+ 26 1.0e-09 1.00e+00 468.000000 507.000000 468.000000 507.000000
+ 27 1.0e-09 1.00e+00 507.000000 548.000000 507.000000 548.000000
+ 28 1.0e-09 1.00e+00 548.000000 592.000000 548.000000 592.000000
+ 29 1.0e-09 1.00e+00 592.000000 638.000000 592.000000 638.000000
+ 30 1.0e-09 1.00e+00 638.000000 686.000000 638.000000 686.000000
+ 31 1.0e-09 1.00e+00 686.000000 737.000000 686.000000 737.000000
+ 32 1.0e-09 1.00e+00 737.000000 846.000000 737.000000 846.000000
+ 33 1.0e-09 1.00e+00 846.000000 1684.000000 846.000000 1684.000000
+ 34 1.0e-09 1.00e+00 18.000000 21.000000 18.000000 21.000000
+ 35 1.0e-09 1.00e+00 21.000000 24.000000 21.000000 24.000000
+ 36 1.0e-09 1.00e+00 24.000000 28.000000 24.000000 28.000000
+ 37 1.0e-09 1.00e+00 28.000000 32.000000 28.000000 32.000000
+ 38 1.0e-09 1.00e+00 32.000000 37.000000 32.000000 37.000000
+ 39 1.0e-09 1.00e+00 37.000000 43.000000 37.000000 43.000000
+ 40 1.0e-09 1.00e+00 43.000000 49.000000 43.000000 49.000000
+ 41 1.0e-09 1.00e+00 49.000000 56.000000 49.000000 56.000000
+ 42 1.0e-09 1.00e+00 56.000000 64.000000 56.000000 64.000000
+ 43 1.0e-09 1.00e+00 64.000000 74.000000 64.000000 74.000000
+ 44 1.0e-09 1.00e+00 74.000000 84.000000 74.000000 84.000000
+ 45 1.0e-09 1.00e+00 84.000000 97.000000 84.000000 97.000000
+ 46 1.0e-09 1.00e+00 97.000000 114.000000 97.000000 114.000000
+ 47 1.0e-09 1.00e+00 114.000000 133.000000 114.000000 133.000000
+ 48 1.0e-09 1.00e+00 133.000000 153.000000 133.000000 153.000000
+ 49 1.0e-09 1.00e+00 153.000000 174.000000 153.000000 174.000000
+ 50 1.0e-09 1.00e+00 174.000000 196.000000 174.000000 196.000000
+ 51 1.0e-09 1.00e+00 196.000000 220.000000 196.000000 220.000000
+ 52 1.0e-09 1.00e+00 220.000000 245.000000 220.000000 245.000000
+ 53 1.0e-09 1.00e+00 245.000000 272.000000 245.000000 272.000000
+ 54 1.0e-09 1.00e+00 272.000000 300.000000 272.000000 300.000000
+ 55 1.0e-09 1.00e+00 300.000000 330.000000 300.000000 330.000000
+ 56 1.0e-09 1.00e+00 330.000000 362.000000 330.000000 362.000000
+ 57 1.0e-09 1.00e+00 362.000000 395.000000 362.000000 395.000000
+ 58 1.0e-09 1.00e+00 395.000000 430.000000 395.000000 430.000000
+ 59 1.0e-09 1.00e+00 430.000000 468.000000 430.000000 468.000000
+ 60 1.0e-09 1.00e+00 468.000000 507.000000 468.000000 507.000000
+ 61 1.0e-09 1.00e+00 507.000000 548.000000 507.000000 548.000000
+ 62 1.0e-09 1.00e+00 548.000000 592.000000 548.000000 592.000000
+ 63 1.0e-09 1.00e+00 592.000000 638.000000 592.000000 638.000000
+ 64 1.0e-09 1.00e+00 638.000000 686.000000 638.000000 686.000000
+ 65 1.0e-09 1.00e+00 686.000000 790.000000 686.000000 790.000000
+ 66 1.0e-09 1.00e+00 790.000000 1684.000000 790.000000 1684.000000
+ 67 1.0e-09 1.00e+00 18.000000 21.000000 18.000000 21.000000
+ 68 1.0e-09 1.00e+00 21.000000 24.000000 21.000000 24.000000
+ 69 1.0e-09 1.00e+00 24.000000 28.000000 24.000000 28.000000
+ 70 1.0e-09 1.00e+00 28.000000 32.000000 28.000000 32.000000
+ 71 1.0e-09 1.00e+00 32.000000 37.000000 32.000000 37.000000
+ 72 1.0e-09 1.00e+00 37.000000 43.000000 37.000000 43.000000
+ 73 1.0e-09 1.00e+00 43.000000 49.000000 43.000000 49.000000
+ 74 1.0e-09 1.00e+00 49.000000 56.000000 49.000000 56.000000
+ 75 1.0e-09 1.00e+00 56.000000 64.000000 56.000000 64.000000
+ 76 1.0e-09 1.00e+00 64.000000 74.000000 64.000000 74.000000
+ 77 1.0e-09 1.00e+00 74.000000 84.000000 74.000000 84.000000
+ 78 1.0e-09 1.00e+00 84.000000 97.000000 84.000000 97.000000
+ 79 1.0e-09 1.00e+00 97.000000 114.000000 97.000000 114.000000
+ 80 1.0e-09 1.00e+00 114.000000 133.000000 114.000000 133.000000
+ 81 1.0e-09 1.00e+00 133.000000 153.000000 133.000000 153.000000
+ 82 1.0e-09 1.00e+00 153.000000 174.000000 153.000000 174.000000
+ 83 1.0e-09 1.00e+00 174.000000 196.000000 174.000000 196.000000
+ 84 1.0e-09 1.00e+00 196.000000 220.000000 196.000000 220.000000
+ 85 1.0e-09 1.00e+00 220.000000 245.000000 220.000000 245.000000
+ 86 1.0e-09 1.00e+00 245.000000 272.000000 245.000000 272.000000
+ 87 1.0e-09 1.00e+00 272.000000 300.000000 272.000000 300.000000
+ 88 1.0e-09 1.00e+00 300.000000 330.000000 300.000000 330.000000
+ 89 1.0e-09 1.00e+00 330.000000 362.000000 330.000000 362.000000
+ 90 1.0e-09 1.00e+00 362.000000 395.000000 362.000000 395.000000
+ 91 1.0e-09 1.00e+00 395.000000 430.000000 395.000000 430.000000
+ 92 1.0e-09 1.00e+00 430.000000 468.000000 430.000000 468.000000
+ 93 1.0e-09 1.00e+00 468.000000 507.000000 468.000000 507.000000
+ 94 1.0e-09 1.00e+00 507.000000 548.000000 507.000000 548.000000
+ 95 1.0e-09 1.00e+00 548.000000 592.000000 548.000000 592.000000
+ 96 1.0e-09 1.00e+00 592.000000 638.000000 592.000000 638.000000
+ 97 1.0e-09 1.00e+00 638.000000 686.000000 638.000000 686.000000
+ 98 1.0e-09 1.00e+00 686.000000 1410.000000 686.000000 1410.000000
+ 99 1.0e-09 1.00e+00 18.000000 21.000000 18.000000 21.000000
+ 100 1.0e-09 1.00e+00 21.000000 24.000000 21.000000 24.000000
+ 101 1.0e-09 1.00e+00 24.000000 28.000000 24.000000 28.000000
+ 102 1.0e-09 1.00e+00 28.000000 32.000000 28.000000 32.000000
+ 103 1.0e-09 1.00e+00 32.000000 37.000000 32.000000 37.000000
+ 104 1.0e-09 1.00e+00 37.000000 43.000000 37.000000 43.000000
+ 105 1.0e-09 1.00e+00 43.000000 49.000000 43.000000 49.000000
+ 106 1.0e-09 1.00e+00 49.000000 56.000000 49.000000 56.000000
+ 107 1.0e-09 1.00e+00 56.000000 64.000000 56.000000 64.000000
+ 108 1.0e-09 1.00e+00 64.000000 74.000000 64.000000 74.000000
+ 109 1.0e-09 1.00e+00 74.000000 84.000000 74.000000 84.000000
+ 110 1.0e-09 1.00e+00 84.000000 97.000000 84.000000 97.000000
+ 111 1.0e-09 1.00e+00 97.000000 114.000000 97.000000 114.000000
+ 112 1.0e-09 1.00e+00 114.000000 133.000000 114.000000 133.000000
+ 113 1.0e-09 1.00e+00 133.000000 153.000000 133.000000 153.000000
+ 114 1.0e-09 1.00e+00 153.000000 174.000000 153.000000 174.000000
+ 115 1.0e-09 1.00e+00 174.000000 196.000000 174.000000 196.000000
+ 116 1.0e-09 1.00e+00 196.000000 220.000000 196.000000 220.000000
+ 117 1.0e-09 1.00e+00 220.000000 245.000000 220.000000 245.000000
+ 118 1.0e-09 1.00e+00 245.000000 272.000000 245.000000 272.000000
+ 119 1.0e-09 1.00e+00 272.000000 300.000000 272.000000 300.000000
+ 120 1.0e-09 1.00e+00 300.000000 330.000000 300.000000 330.000000
+ 121 1.0e-09 1.00e+00 330.000000 362.000000 330.000000 362.000000
+ 122 1.0e-09 1.00e+00 362.000000 395.000000 362.000000 395.000000
+ 123 1.0e-09 1.00e+00 395.000000 430.000000 395.000000 430.000000
+ 124 1.0e-09 1.00e+00 430.000000 468.000000 430.000000 468.000000
+ 125 1.0e-09 1.00e+00 468.000000 507.000000 468.000000 507.000000
+ 126 1.0e-09 1.00e+00 507.000000 548.000000 507.000000 548.000000
+ 127 1.0e-09 1.00e+00 548.000000 1032.000000 548.000000 1032.000000
+ 128 1.0e-09 1.00e+00 18.000000 21.000000 18.000000 21.000000
+ 129 1.0e-09 1.00e+00 21.000000 24.000000 21.000000 24.000000
+ 130 1.0e-09 1.00e+00 24.000000 28.000000 24.000000 28.000000
+ 131 1.0e-09 1.00e+00 28.000000 32.000000 28.000000 32.000000
+ 132 1.0e-09 1.00e+00 32.000000 37.000000 32.000000 37.000000
+ 133 1.0e-09 1.00e+00 37.000000 43.000000 37.000000 43.000000
+ 134 1.0e-09 1.00e+00 43.000000 49.000000 43.000000 49.000000
+ 135 1.0e-09 1.00e+00 49.000000 56.000000 49.000000 56.000000
+ 136 1.0e-09 1.00e+00 56.000000 64.000000 56.000000 64.000000
+ 137 1.0e-09 1.00e+00 64.000000 74.000000 64.000000 74.000000
+ 138 1.0e-09 1.00e+00 74.000000 84.000000 74.000000 84.000000
+ 139 1.0e-09 1.00e+00 84.000000 97.000000 84.000000 97.000000
+ 140 1.0e-09 1.00e+00 97.000000 114.000000 97.000000 114.000000
+ 141 1.0e-09 1.00e+00 114.000000 133.000000 114.000000 133.000000
+ 142 1.0e-09 1.00e+00 133.000000 153.000000 133.000000 153.000000
+ 143 1.0e-09 1.00e+00 153.000000 174.000000 153.000000 174.000000
+ 144 1.0e-09 1.00e+00 174.000000 196.000000 174.000000 196.000000
+ 145 1.0e-09 1.00e+00 196.000000 220.000000 196.000000 220.000000
+ 146 1.0e-09 1.00e+00 220.000000 245.000000 220.000000 245.000000
+ 147 1.0e-09 1.00e+00 245.000000 272.000000 245.000000 272.000000
+ 148 1.0e-09 1.00e+00 272.000000 300.000000 272.000000 300.000000
+ 149 1.0e-09 1.00e+00 300.000000 330.000000 300.000000 330.000000
+ 150 1.0e-09 1.00e+00 330.000000 362.000000 330.000000 362.000000
+ 151 1.0e-09 1.00e+00 362.000000 395.000000 362.000000 395.000000
+ 152 1.0e-09 1.00e+00 395.000000 430.000000 395.000000 430.000000
+ 153 1.0e-09 1.00e+00 430.000000 737.000000 430.000000 737.000000
+ 154 1.0e-09 1.00e+00 18.000000 21.000000 18.000000 21.000000
+ 155 1.0e-09 1.00e+00 21.000000 24.000000 21.000000 24.000000
+ 156 1.0e-09 1.00e+00 24.000000 28.000000 24.000000 28.000000
+ 157 1.0e-09 1.00e+00 28.000000 32.000000 28.000000 32.000000
+ 158 1.0e-09 1.00e+00 32.000000 37.000000 32.000000 37.000000
+ 159 1.0e-09 1.00e+00 37.000000 43.000000 37.000000 43.000000
+ 160 1.0e-09 1.00e+00 43.000000 49.000000 43.000000 49.000000
+ 161 1.0e-09 1.00e+00 49.000000 56.000000 49.000000 56.000000
+ 162 1.0e-09 1.00e+00 56.000000 64.000000 56.000000 64.000000
+ 163 1.0e-09 1.00e+00 64.000000 74.000000 64.000000 74.000000
+ 164 1.0e-09 1.00e+00 74.000000 84.000000 74.000000 84.000000
+ 165 1.0e-09 1.00e+00 84.000000 97.000000 84.000000 97.000000
+ 166 1.0e-09 1.00e+00 97.000000 114.000000 97.000000 114.000000
+ 167 1.0e-09 1.00e+00 114.000000 133.000000 114.000000 133.000000
+ 168 1.0e-09 1.00e+00 133.000000 153.000000 133.000000 153.000000
+ 169 1.0e-09 1.00e+00 153.000000 174.000000 153.000000 174.000000
+ 170 1.0e-09 1.00e+00 174.000000 196.000000 174.000000 196.000000
+ 171 1.0e-09 1.00e+00 196.000000 220.000000 196.000000 220.000000
+ 172 1.0e-09 1.00e+00 220.000000 245.000000 220.000000 245.000000
+ 173 1.0e-09 1.00e+00 245.000000 272.000000 245.000000 272.000000
+ 174 1.0e-09 1.00e+00 272.000000 300.000000 272.000000 300.000000
+ 175 1.0e-09 1.00e+00 300.000000 468.000000 300.000000 468.000000
+}}
+
+
+Warmup.Binning {{
+ ObsBin |y|_Lo |y|_Up pT_[GeV]_Lo pT_[GeV]_Up BinSize EventCount
+ 0 0.00000 0.500000 18.0000 21.0000 3.00000 308
+ 1 0.00000 0.500000 21.0000 24.0000 3.00000 399
+ 2 0.00000 0.500000 24.0000 28.0000 4.00000 511
+ 3 0.00000 0.500000 28.0000 32.0000 4.00000 609
+ 4 0.00000 0.500000 32.0000 37.0000 5.00000 1463
+ 5 0.00000 0.500000 37.0000 43.0000 6.00000 4291
+ 6 0.00000 0.500000 43.0000 49.0000 6.00000 6615
+ 7 0.00000 0.500000 49.0000 56.0000 7.00000 9079
+ 8 0.00000 0.500000 56.0000 64.0000 8.00000 10969
+ 9 0.00000 0.500000 64.0000 74.0000 10.0000 14035
+ 10 0.00000 0.500000 74.0000 84.0000 10.0000 12929
+ 11 0.00000 0.500000 84.0000 97.0000 13.0000 16940
+ 12 0.00000 0.500000 97.0000 114.000 17.0000 19691
+ 13 0.00000 0.500000 114.000 133.000 19.0000 18921
+ 14 0.00000 0.500000 133.000 153.000 20.0000 17773
+ 15 0.00000 0.500000 153.000 174.000 21.0000 16065
+ 16 0.00000 0.500000 174.000 196.000 22.0000 14035
+ 17 0.00000 0.500000 196.000 220.000 24.0000 14273
+ 18 0.00000 0.500000 220.000 245.000 25.0000 12579
+ 19 0.00000 0.500000 245.000 272.000 27.0000 12159
+ 20 0.00000 0.500000 272.000 300.000 28.0000 11200
+ 21 0.00000 0.500000 300.000 330.000 30.0000 10556
+ 22 0.00000 0.500000 330.000 362.000 32.0000 9506
+ 23 0.00000 0.500000 362.000 395.000 33.0000 10388
+ 24 0.00000 0.500000 395.000 430.000 35.0000 9324
+ 25 0.00000 0.500000 430.000 468.000 38.0000 8883
+ 26 0.00000 0.500000 468.000 507.000 39.0000 8106
+ 27 0.00000 0.500000 507.000 548.000 41.0000 8358
+ 28 0.00000 0.500000 548.000 592.000 44.0000 7973
+ 29 0.00000 0.500000 592.000 638.000 46.0000 7665
+ 30 0.00000 0.500000 638.000 686.000 48.0000 7630
+ 31 0.00000 0.500000 686.000 737.000 51.0000 7413
+ 32 0.00000 0.500000 737.000 846.000 109.000 14644
+ 33 0.00000 0.500000 846.000 1684.00 838.000 56553
+ 34 0.500000 1.00000 18.0000 21.0000 3.00000 364
+ 35 0.500000 1.00000 21.0000 24.0000 3.00000 385
+ 36 0.500000 1.00000 24.0000 28.0000 4.00000 994
+ 37 0.500000 1.00000 28.0000 32.0000 4.00000 2114
+ 38 0.500000 1.00000 32.0000 37.0000 5.00000 3801
+ 39 0.500000 1.00000 37.0000 43.0000 6.00000 5705
+ 40 0.500000 1.00000 43.0000 49.0000 6.00000 5264
+ 41 0.500000 1.00000 49.0000 56.0000 7.00000 6468
+ 42 0.500000 1.00000 56.0000 64.0000 8.00000 7175
+ 43 0.500000 1.00000 64.0000 74.0000 10.0000 8778
+ 44 0.500000 1.00000 74.0000 84.0000 10.0000 8519
+ 45 0.500000 1.00000 84.0000 97.0000 13.0000 11802
+ 46 0.500000 1.00000 97.0000 114.000 17.0000 14868
+ 47 0.500000 1.00000 114.000 133.000 19.0000 16394
+ 48 0.500000 1.00000 133.000 153.000 20.0000 16247
+ 49 0.500000 1.00000 153.000 174.000 21.0000 15547
+ 50 0.500000 1.00000 174.000 196.000 22.0000 14203
+ 51 0.500000 1.00000 196.000 220.000 24.0000 14105
+ 52 0.500000 1.00000 220.000 245.000 25.0000 11788
+ 53 0.500000 1.00000 245.000 272.000 27.0000 12047
+ 54 0.500000 1.00000 272.000 300.000 28.0000 10794
+ 55 0.500000 1.00000 300.000 330.000 30.0000 10227
+ 56 0.500000 1.00000 330.000 362.000 32.0000 9842
+ 57 0.500000 1.00000 362.000 395.000 33.0000 9366
+ 58 0.500000 1.00000 395.000 430.000 35.0000 9303
+ 59 0.500000 1.00000 430.000 468.000 38.0000 9205
+ 60 0.500000 1.00000 468.000 507.000 39.0000 7644
+ 61 0.500000 1.00000 507.000 548.000 41.0000 8162
+ 62 0.500000 1.00000 548.000 592.000 44.0000 7805
+ 63 0.500000 1.00000 592.000 638.000 46.0000 7504
+ 64 0.500000 1.00000 638.000 686.000 48.0000 6853
+ 65 0.500000 1.00000 686.000 790.000 104.000 13265
+ 66 0.500000 1.00000 790.000 1684.00 894.000 53438
+ 67 1.00000 1.50000 18.0000 21.0000 3.00000 1148
+ 68 1.00000 1.50000 21.0000 24.0000 3.00000 1617
+ 69 1.00000 1.50000 24.0000 28.0000 4.00000 2345
+ 70 1.00000 1.50000 28.0000 32.0000 4.00000 2338
+ 71 1.00000 1.50000 32.0000 37.0000 5.00000 2758
+ 72 1.00000 1.50000 37.0000 43.0000 6.00000 3661
+ 73 1.00000 1.50000 43.0000 49.0000 6.00000 3395
+ 74 1.00000 1.50000 49.0000 56.0000 7.00000 4151
+ 75 1.00000 1.50000 56.0000 64.0000 8.00000 4613
+ 76 1.00000 1.50000 64.0000 74.0000 10.0000 5411
+ 77 1.00000 1.50000 74.0000 84.0000 10.0000 5418
+ 78 1.00000 1.50000 84.0000 97.0000 13.0000 7252
+ 79 1.00000 1.50000 97.0000 114.000 17.0000 9149
+ 80 1.00000 1.50000 114.000 133.000 19.0000 10311
+ 81 1.00000 1.50000 133.000 153.000 20.0000 10822
+ 82 1.00000 1.50000 153.000 174.000 21.0000 10920
+ 83 1.00000 1.50000 174.000 196.000 22.0000 11760
+ 84 1.00000 1.50000 196.000 220.000 24.0000 12467
+ 85 1.00000 1.50000 220.000 245.000 25.0000 11557
+ 86 1.00000 1.50000 245.000 272.000 27.0000 11571
+ 87 1.00000 1.50000 272.000 300.000 28.0000 9898
+ 88 1.00000 1.50000 300.000 330.000 30.0000 10136
+ 89 1.00000 1.50000 330.000 362.000 32.0000 9618
+ 90 1.00000 1.50000 362.000 395.000 33.0000 8505
+ 91 1.00000 1.50000 395.000 430.000 35.0000 8043
+ 92 1.00000 1.50000 430.000 468.000 38.0000 8302
+ 93 1.00000 1.50000 468.000 507.000 39.0000 6783
+ 94 1.00000 1.50000 507.000 548.000 41.0000 6825
+ 95 1.00000 1.50000 548.000 592.000 44.0000 6867
+ 96 1.00000 1.50000 592.000 638.000 46.0000 6496
+ 97 1.00000 1.50000 638.000 686.000 48.0000 5929
+ 98 1.00000 1.50000 686.000 1410.00 724.000 41097
+ 99 1.50000 2.00000 18.0000 21.0000 3.00000 1267
+ 100 1.50000 2.00000 21.0000 24.0000 3.00000 1274
+ 101 1.50000 2.00000 24.0000 28.0000 4.00000 1498
+ 102 1.50000 2.00000 28.0000 32.0000 4.00000 1547
+ 103 1.50000 2.00000 32.0000 37.0000 5.00000 1869
+ 104 1.50000 2.00000 37.0000 43.0000 6.00000 2527
+ 105 1.50000 2.00000 43.0000 49.0000 6.00000 2359
+ 106 1.50000 2.00000 49.0000 56.0000 7.00000 2478
+ 107 1.50000 2.00000 56.0000 64.0000 8.00000 2870
+ 108 1.50000 2.00000 64.0000 74.0000 10.0000 3514
+ 109 1.50000 2.00000 74.0000 84.0000 10.0000 3507
+ 110 1.50000 2.00000 84.0000 97.0000 13.0000 4480
+ 111 1.50000 2.00000 97.0000 114.000 17.0000 5866
+ 112 1.50000 2.00000 114.000 133.000 19.0000 6685
+ 113 1.50000 2.00000 133.000 153.000 20.0000 6545
+ 114 1.50000 2.00000 153.000 174.000 21.0000 6867
+ 115 1.50000 2.00000 174.000 196.000 22.0000 6363
+ 116 1.50000 2.00000 196.000 220.000 24.0000 7707
+ 117 1.50000 2.00000 220.000 245.000 25.0000 7287
+ 118 1.50000 2.00000 245.000 272.000 27.0000 8673
+ 119 1.50000 2.00000 272.000 300.000 28.0000 8316
+ 120 1.50000 2.00000 300.000 330.000 30.0000 8743
+ 121 1.50000 2.00000 330.000 362.000 32.0000 7875
+ 122 1.50000 2.00000 362.000 395.000 33.0000 7511
+ 123 1.50000 2.00000 395.000 430.000 35.0000 6720
+ 124 1.50000 2.00000 430.000 468.000 38.0000 6419
+ 125 1.50000 2.00000 468.000 507.000 39.0000 5586
+ 126 1.50000 2.00000 507.000 548.000 41.0000 5355
+ 127 1.50000 2.00000 548.000 1032.00 484.000 28231
+ 128 2.00000 2.50000 18.0000 21.0000 3.00000 854
+ 129 2.00000 2.50000 21.0000 24.0000 3.00000 861
+ 130 2.00000 2.50000 24.0000 28.0000 4.00000 1218
+ 131 2.00000 2.50000 28.0000 32.0000 4.00000 875
+ 132 2.00000 2.50000 32.0000 37.0000 5.00000 1211
+ 133 2.00000 2.50000 37.0000 43.0000 6.00000 1498
+ 134 2.00000 2.50000 43.0000 49.0000 6.00000 1470
+ 135 2.00000 2.50000 49.0000 56.0000 7.00000 1631
+ 136 2.00000 2.50000 56.0000 64.0000 8.00000 1673
+ 137 2.00000 2.50000 64.0000 74.0000 10.0000 2310
+ 138 2.00000 2.50000 74.0000 84.0000 10.0000 2205
+ 139 2.00000 2.50000 84.0000 97.0000 13.0000 2961
+ 140 2.00000 2.50000 97.0000 114.000 17.0000 3619
+ 141 2.00000 2.50000 114.000 133.000 19.0000 4032
+ 142 2.00000 2.50000 133.000 153.000 20.0000 4326
+ 143 2.00000 2.50000 153.000 174.000 21.0000 4585
+ 144 2.00000 2.50000 174.000 196.000 22.0000 4571
+ 145 2.00000 2.50000 196.000 220.000 24.0000 4592
+ 146 2.00000 2.50000 220.000 245.000 25.0000 5327
+ 147 2.00000 2.50000 245.000 272.000 27.0000 5509
+ 148 2.00000 2.50000 272.000 300.000 28.0000 5124
+ 149 2.00000 2.50000 300.000 330.000 30.0000 5747
+ 150 2.00000 2.50000 330.000 362.000 32.0000 5733
+ 151 2.00000 2.50000 362.000 395.000 33.0000 5362
+ 152 2.00000 2.50000 395.000 430.000 35.0000 4802
+ 153 2.00000 2.50000 430.000 737.000 307.000 14854
+ 154 2.50000 3.00000 18.0000 21.0000 3.00000 504
+ 155 2.50000 3.00000 21.0000 24.0000 3.00000 427
+ 156 2.50000 3.00000 24.0000 28.0000 4.00000 602
+ 157 2.50000 3.00000 28.0000 32.0000 4.00000 672
+ 158 2.50000 3.00000 32.0000 37.0000 5.00000 854
+ 159 2.50000 3.00000 37.0000 43.0000 6.00000 917
+ 160 2.50000 3.00000 43.0000 49.0000 6.00000 910
+ 161 2.50000 3.00000 49.0000 56.0000 7.00000 1211
+ 162 2.50000 3.00000 56.0000 64.0000 8.00000 1211
+ 163 2.50000 3.00000 64.0000 74.0000 10.0000 1449
+ 164 2.50000 3.00000 74.0000 84.0000 10.0000 1470
+ 165 2.50000 3.00000 84.0000 97.0000 13.0000 1841
+ 166 2.50000 3.00000 97.0000 114.000 17.0000 2282
+ 167 2.50000 3.00000 114.000 133.000 19.0000 2590
+ 168 2.50000 3.00000 133.000 153.000 20.0000 2674
+ 169 2.50000 3.00000 153.000 174.000 21.0000 2779
+ 170 2.50000 3.00000 174.000 196.000 22.0000 2730
+ 171 2.50000 3.00000 196.000 220.000 24.0000 3283
+ 172 2.50000 3.00000 220.000 245.000 25.0000 2919
+ 173 2.50000 3.00000 245.000 272.000 27.0000 3227
+ 174 2.50000 3.00000 272.000 300.000 28.0000 3255
+ 175 2.50000 3.00000 300.000 468.000 168.000 8883
+}}
diff --git a/test/compare_tables.py b/test/compare_tables.py
new file mode 100644
index 0000000000000000000000000000000000000000..b88e19588cfffe2d2bc8ab0fc28cc01465e94886
--- /dev/null
+++ b/test/compare_tables.py
@@ -0,0 +1,51 @@
+#!/usr/bin/env python3
+
+import os
+
+import fastnlo as fnlo
+
+TYPES = ["fix", "flex"]
+NAMES = ["NodesPerBin", "NodeDensity"]
+SEEDS = [1000, 2000]
+
+# fnlo.SetGlobalVerbosity(fnlo.DEBUG)
+
+os.system("rm -rf output")
+
+
+def generate_table(vartype: str, name: str, seed: int):
+ os.system(f"ln -sf {vartype}_{name}.str InclusiveNJets.str")
+ os.system(
+ "nlojet++ "
+ "--calculate "
+ "-c born "
+ f"-s {seed} "
+ f"--name {vartype}_{name}_{seed:06d} "
+ " -u /usr/lib/fastnlo_interface_nlojet/libInclusiveNJets.la "
+ "--max-event 100000 "
+ "--save-after 100000 "
+ )
+ os.system(f"mv output/run-hhc-born-2jet.tab.gz output/{vartype}_{name}_{seed:06d}.tab.gz")
+
+
+for vartype in TYPES:
+ for steering_file in NAMES:
+ for seed in SEEDS:
+ generate_table(vartype, steering_file, seed)
+ os.system(f"fnlo-tk-merge2 output/{vartype}_{steering_file}_??????.tab.gz "
+ f"output/{vartype}_{steering_file}_merged.tab.gz")
+
+ for seed in SEEDS:
+ filename_0 = f"output/{vartype}_{NAMES[0]}_{seed:06d}.tab.gz"
+ filename_1 = f"output/{vartype}_{NAMES[1]}_{seed:06d}.tab.gz"
+ table_0 = fnlo.fastNLOTable(filename_0)
+ table_1 = fnlo.fastNLOTable(filename_1)
+ print(f"Asserting equal: {filename_0} <-> {filename_1}")
+ assert table_0.IsEquivalent(table_1, 1e-8)
+
+ filename_0 = f"output/{vartype}_{NAMES[0]}_merged.tab.gz"
+ filename_1 = f"output/{vartype}_{NAMES[1]}_merged.tab.gz"
+ table_0_merged = fnlo.fastNLOTable(filename_0)
+ table_1_merged = fnlo.fastNLOTable(filename_1)
+ print(f"Asserting equal: {filename_0} <-> {filename_1}")
+ assert table_0_merged.IsEquivalent(table_1_merged, 1e-8)
diff --git a/test/fix_NodeDensity.str b/test/fix_NodeDensity.str
new file mode 100644
index 0000000000000000000000000000000000000000..8d2bfbd28fe3a70230e602d356942005c7c73529
--- /dev/null
+++ b/test/fix_NodeDensity.str
@@ -0,0 +1,207 @@
+# -*-sh-*-
+# ==================================================================== #
+#
+# A steering file for creating a fastNLO table (with NLOJet++)
+#
+# ==================================================================== #
+
+# ==================================================================== #
+#
+# Scenario specific settings
+#
+# ==================================================================== #
+
+# -------------------------------------------------------------------- #
+# Description of scenario
+# (Write each line of description in quotation marks)
+# -------------------------------------------------------------------- #
+ScenarioName InclusiveJets_fix
+ScenarioDescription {
+ "d2sigma-jet_dpT_dy_[pb_GeV]"
+ "CMS_Collaboration"
+ "Inclusive_Jet_pT"
+ "anti-kT_R=0.5"
+ "CMS-PAPER-QCD-10-011, arXiv:1106.0208, Phys. Rev. Lett. 107 (2011) 132001."
+ "RIVET_ID=CMS_2011_S9086218/D01-x01-y01,SPIRES_ID=9086218,INSPIRE_RECORD=902309"
+ "THEORY_LABEL=NLOJet++"
+ "provided by:"
+ "fastNLO_2.5.0"
+ "If you use this table, please cite:"
+ " D. Britzger, T. Kluge, K. Rabbertz, F. Stober, M. Wobisch, arXiv:1109.1310"
+}
+PublicationUnits 12 # Unit of data cross sections (negative power of 10, e.g. 12->pb, 15->fb)
+
+# -------------------------------------------------------------------- #
+# Global scenario settings
+# -------------------------------------------------------------------- #
+JetAlgo 2 # fastjet jet algorithm
+# # clustering jet algos: 0 = kT, 1 = CA, 2 = anti-kT
+# # cone jet algos: 10 = SISCone, 11 = CDFMidPointCone, 12 = D0RunIICone
+Rjet 0.5 # Jet size parameter: Required for all jets
+#OvThr 0.5 # Overlap threshold: Required for cone jets (Normally 10: 0.75; 11,12: 0.50)
+# # Overall jet phase space
+ptjmin 18. # Minimal jet pT (minimum set in interface to fastjet: 1 GeV)
+yjmin 0.0 # Minimal jet rapidity
+yjmax 3.0 # Maximal jet rapidity
+# Use either rapidity or pseudorapidity but not both! If both are set, only rapidity cuts are applied.
+#etajmin 0.0 # Minimal jet pseudorapidity
+#etajmax 3.0 # Maximal jet pseudorapidity
+#Njetmin 1 # Minimal number of jets required to be in preselected phase space
+# # (Normally one, i.e. single-jet inclusive)
+#obs0min x.y #
+#obs0max x.y # By default these values are derived from the binning in these
+#obs1min x.y # observables and are set automatically. If required the min/max
+#obs1max x.y # values for each of the maximally three observables can be set
+#obs2min x.y # by hand.
+#obs2max x.y #
+
+# -------------------------------------------------------------------- #
+# Observable definition, dimensioning, and normalization
+# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+# Specify, how many differential dimensions are measured.
+# Number must be consistent with DimensionLabels and DimensionIsDifferential, see below
+# 1 : read in 'SingleDifferentialBinning'
+# 2 : read in 'DoubleDifferentialBinning'
+# 3 : read in 'TripleDifferentialBinning'
+# -------------------------------------------------------------------- #
+LeadingOrder 2 # Number of jets for the LO process
+DifferentialDimension 2 # Dimensionality of binning
+DimensionLabels { # Labels (symbol and unit) for the measurement dimensions (from outer to inner "loop")
+ "|y|" # This defines the observables to be calculated!
+ "pT_[GeV]" # Implemented observables are: "|y|", "pT_[GeV]", "phi"
+}
+DimensionIsDifferential { # Specify for each dimension whether
+ 2 2 # 0 : the cross section is NOT differential, i.e. there are two bin borders,
+} # but NO division (normalization) by bin width
+ # 1 : the cross section is point-wise differential, i.e. only one point is given
+ # 2 : the cross section is bin-wise differential, i.e. there are two bin borders
+ # and division by bin width
+CalculateBinSize true # Calculate bin width from lower and upper bin boundaries
+BinSizeFactor 2. # Possibility to provide additional normalization factor, e.g. of 2 for bins in |y|
+#BinSize { x1 x2 x3 ... } # If 'CalculateBinSize' is 'false' provide table with bin widths 'by hand' for normalization
+ #
+ # If the calculation should not be divided by bin width, then use 'DimensionIsDifferential' equal '0'
+ # and set 'CalculateBinSize' 'true' for each dimension.
+ #
+ScaleDescriptionScale1 "pT_jet_[GeV]" # This defines the scale to be used (Note: The 1st scale must always be in units of [GeV]!)
+ScaleDescriptionScale2 "pT_max_[GeV]" # Specify 2nd scale name and unit (ONLY for flexible-scale tables)
+ # Implemented scales are: "pT_jet_[GeV]", "pT_max_[GeV]"
+
+# -------------------------------------------------------------------- #
+# Observable binning
+# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+# Use either 'SingleDifferentialBinning' or
+# 'DoubleDifferentialBinning' or
+# 'TripleDifferentialBinning'
+# in accord with 'DifferentialDimension' above
+# -------------------------------------------------------------------- #
+# Example for SingleDifferentialBinning:
+#SingleDifferentialBinning {
+# 7. 11. 18. 30. 50. 100.
+#}
+
+# Example for DoubleDifferentialBinning:
+DoubleDifferentialBinning {{
+ 1stDimLo 1stDimUp "----- Array of bin-grid for 2nd dimension -----"
+ 0.0 0.5 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 592. 638. 686. 737. 846. 1684.
+ 0.5 1.0 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 592. 638. 686. 790. 1684.
+ 1.0 1.5 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 592. 638. 686. 1410.
+ 1.5 2.0 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 1032.
+ 2.0 2.5 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 737.
+ 2.5 3.0 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 468.
+}}
+
+# Example for TripleDifferentialBinning:
+#TripleDifferentialBinning {{
+# 1stDimLo 1stDimUp 2ndDimLo 2ndDimUp "----- Array of bin-grid for 3rd dimension -----"
+# 100. 150. -1.0 0.0 7. 11. 18. 30. 50.
+# 150. 200. -1.0 0.0 7. 11. 18. 30. 50.
+# 200. 400. 0.0 1.5 7. 11. 18. 30. 50.
+# 400. 700. 0.0 1.5 7. 11. 18. 30. 50.
+# 700. 5000. 1.5 2.5 7. 11. 18. 30. 50.
+# 5000. 15000. 1.5 2.5 7. 11. 18. 30. 50.
+#}}
+
+# ==================================================================== #
+#
+# Process specific and therefore at least partially theory-code
+# dependent settings
+#
+# ==================================================================== #
+CenterOfMassEnergy 7000. # Center-of-mass energy in GeV.
+# Some examples from the past: 200. # RHIC
+# 1800. # TeV Run I
+# 1960. # TeV Run II
+# 900. # LHC Injection Run
+# 2360. # LHC Initial Run
+# 2760. # LHC HIpp base Run
+# 5020. # LHC HIpp base Run
+# 7000. # LHC First Run (2010/2011)
+# 8000. # LHC Second Run (2012)
+# 13000. # LHC Next Run ? (2015)
+# 14000. # LHC Design Run
+PDF1 2212 # PDF of 1st hadron (following PDG convention).
+PDF2 2212 # PDF of 2nd hadron (following PDG convention).
+NPDF 2
+NPDFDim 1
+
+# ==================================================================== #
+#
+# General fastNLO settings
+#
+# ==================================================================== #
+#GlobalVerbosity DEBUG # Global output verbosity of fastNLO toolkit, def.=WARNING,
+GlobalVerbosity INFO # Global output verbosity of fastNLO toolkit, def.=WARNING,
+ # (and of user scenario code if implemented).
+ # Possible values are: DEBUG, MANUAL, INFO, WARNING, ERROR, SILENT
+OutputFilename InclusiveNJets.tab # Filename of fastNLO output table
+OutputPrecision 8 # Number of decimal digits to store in output table, def.=8.
+FlexibleScaleTable false # Create table fully flexible in mu_f (larger size), true, or
+ # table with fixed number of mu_f scale factors, def.=false.
+ScaleVariationFactors { # Factorization scale variations (only needed for fixed-scale tables)
+ 0.5 1.0 2.0 # List of scale factors must include factor '1'
+} # Scale factors will be ordered according to fastNLO convention: (1, min, ... , max)
+ReadBinningFromSteering true # Specify if binning is read from this file. MUST be true for now!
+IgnoreWarmupBinningCheck true # Ignore crosscheck of binning between steering and warmup
+ # Avoids problems with precision of bin borders e.g. in the case of Pi
+#>>fastNLO-warmup.txt # Specify alternative filename to read warm-up results from, def.=SteerFileBasename_${ScenarioName}_warmup.txt
+
+# ==================================================================== #
+#
+# Specific fastNLO settings (normally do not need to be changed)
+# Only change when you know what you are doing ...!
+#
+# ==================================================================== #
+ApplyPDFReweighting true # Apply reweighting of pdfs for an optimized interpolation, def.=true.
+CheckScaleLimitsAgainstBins true # Set limits for scale nodes to bin borders, if possible
+
+# -------------------------------------------------------------------- #
+# Choose fastNLO interpolation kernels and distance measures
+# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+# Currently implemented interpolation kernels
+# Catmull
+# Lagrange
+# OneNode
+# Linear
+#
+# Currently implemented distance measures
+# linear
+# loglog025 eq. to (log(log(4*x)))
+# log10
+# sqrtlog10 eq. to sqrt(log_10(x))
+# -------------------------------------------------------------------- #
+X_Kernel Lagrange
+X_DistanceMeasure "sqrtlog10"
+X_NNodes 4
+X_NoOfNodesPerMagnitude false
+X_NNodeCounting "NodeDensity"
+
+Mu1_Kernel Lagrange
+Mu1_DistanceMeasure "loglog025"
+Mu1_NNodes 6
+Mu1_NNodeCounting "NodesPerDim2"
+
+Mu2_Kernel Lagrange # Scale2 not used for fixed-scale tables
+Mu2_DistanceMeasure "loglog025"
+Mu2_NNodes 6
+Mu2_NNodeCounting "NodesPerDim2"
diff --git a/test/fix_NodesPerBin.str b/test/fix_NodesPerBin.str
new file mode 100644
index 0000000000000000000000000000000000000000..2a45d07baa6a7f6fd66ba8a0678dbf1b001a1799
--- /dev/null
+++ b/test/fix_NodesPerBin.str
@@ -0,0 +1,207 @@
+# -*-sh-*-
+# ==================================================================== #
+#
+# A steering file for creating a fastNLO table (with NLOJet++)
+#
+# ==================================================================== #
+
+# ==================================================================== #
+#
+# Scenario specific settings
+#
+# ==================================================================== #
+
+# -------------------------------------------------------------------- #
+# Description of scenario
+# (Write each line of description in quotation marks)
+# -------------------------------------------------------------------- #
+ScenarioName InclusiveJets_fix
+ScenarioDescription {
+ "d2sigma-jet_dpT_dy_[pb_GeV]"
+ "CMS_Collaboration"
+ "Inclusive_Jet_pT"
+ "anti-kT_R=0.5"
+ "CMS-PAPER-QCD-10-011, arXiv:1106.0208, Phys. Rev. Lett. 107 (2011) 132001."
+ "RIVET_ID=CMS_2011_S9086218/D01-x01-y01,SPIRES_ID=9086218,INSPIRE_RECORD=902309"
+ "THEORY_LABEL=NLOJet++"
+ "provided by:"
+ "fastNLO_2.5.0"
+ "If you use this table, please cite:"
+ " D. Britzger, T. Kluge, K. Rabbertz, F. Stober, M. Wobisch, arXiv:1109.1310"
+}
+PublicationUnits 12 # Unit of data cross sections (negative power of 10, e.g. 12->pb, 15->fb)
+
+# -------------------------------------------------------------------- #
+# Global scenario settings
+# -------------------------------------------------------------------- #
+JetAlgo 2 # fastjet jet algorithm
+# # clustering jet algos: 0 = kT, 1 = CA, 2 = anti-kT
+# # cone jet algos: 10 = SISCone, 11 = CDFMidPointCone, 12 = D0RunIICone
+Rjet 0.5 # Jet size parameter: Required for all jets
+#OvThr 0.5 # Overlap threshold: Required for cone jets (Normally 10: 0.75; 11,12: 0.50)
+# # Overall jet phase space
+ptjmin 18. # Minimal jet pT (minimum set in interface to fastjet: 1 GeV)
+yjmin 0.0 # Minimal jet rapidity
+yjmax 3.0 # Maximal jet rapidity
+# Use either rapidity or pseudorapidity but not both! If both are set, only rapidity cuts are applied.
+#etajmin 0.0 # Minimal jet pseudorapidity
+#etajmax 3.0 # Maximal jet pseudorapidity
+#Njetmin 1 # Minimal number of jets required to be in preselected phase space
+# # (Normally one, i.e. single-jet inclusive)
+#obs0min x.y #
+#obs0max x.y # By default these values are derived from the binning in these
+#obs1min x.y # observables and are set automatically. If required the min/max
+#obs1max x.y # values for each of the maximally three observables can be set
+#obs2min x.y # by hand.
+#obs2max x.y #
+
+# -------------------------------------------------------------------- #
+# Observable definition, dimensioning, and normalization
+# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+# Specify, how many differential dimensions are measured.
+# Number must be consistent with DimensionLabels and DimensionIsDifferential, see below
+# 1 : read in 'SingleDifferentialBinning'
+# 2 : read in 'DoubleDifferentialBinning'
+# 3 : read in 'TripleDifferentialBinning'
+# -------------------------------------------------------------------- #
+LeadingOrder 2 # Number of jets for the LO process
+DifferentialDimension 2 # Dimensionality of binning
+DimensionLabels { # Labels (symbol and unit) for the measurement dimensions (from outer to inner "loop")
+ "|y|" # This defines the observables to be calculated!
+ "pT_[GeV]" # Implemented observables are: "|y|", "pT_[GeV]", "phi"
+}
+DimensionIsDifferential { # Specify for each dimension whether
+ 2 2 # 0 : the cross section is NOT differential, i.e. there are two bin borders,
+} # but NO division (normalization) by bin width
+ # 1 : the cross section is point-wise differential, i.e. only one point is given
+ # 2 : the cross section is bin-wise differential, i.e. there are two bin borders
+ # and division by bin width
+CalculateBinSize true # Calculate bin width from lower and upper bin boundaries
+BinSizeFactor 2. # Possibility to provide additional normalization factor, e.g. of 2 for bins in |y|
+#BinSize { x1 x2 x3 ... } # If 'CalculateBinSize' is 'false' provide table with bin widths 'by hand' for normalization
+ #
+ # If the calculation should not be divided by bin width, then use 'DimensionIsDifferential' equal '0'
+ # and set 'CalculateBinSize' 'true' for each dimension.
+ #
+ScaleDescriptionScale1 "pT_jet_[GeV]" # This defines the scale to be used (Note: The 1st scale must always be in units of [GeV]!)
+ScaleDescriptionScale2 "pT_max_[GeV]" # Specify 2nd scale name and unit (ONLY for flexible-scale tables)
+ # Implemented scales are: "pT_jet_[GeV]", "pT_max_[GeV]"
+
+# -------------------------------------------------------------------- #
+# Observable binning
+# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+# Use either 'SingleDifferentialBinning' or
+# 'DoubleDifferentialBinning' or
+# 'TripleDifferentialBinning'
+# in accord with 'DifferentialDimension' above
+# -------------------------------------------------------------------- #
+# Example for SingleDifferentialBinning:
+#SingleDifferentialBinning {
+# 7. 11. 18. 30. 50. 100.
+#}
+
+# Example for DoubleDifferentialBinning:
+DoubleDifferentialBinning {{
+ 1stDimLo 1stDimUp "----- Array of bin-grid for 2nd dimension -----"
+ 0.0 0.5 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 592. 638. 686. 737. 846. 1684.
+ 0.5 1.0 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 592. 638. 686. 790. 1684.
+ 1.0 1.5 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 592. 638. 686. 1410.
+ 1.5 2.0 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 1032.
+ 2.0 2.5 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 737.
+ 2.5 3.0 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 468.
+}}
+
+# Example for TripleDifferentialBinning:
+#TripleDifferentialBinning {{
+# 1stDimLo 1stDimUp 2ndDimLo 2ndDimUp "----- Array of bin-grid for 3rd dimension -----"
+# 100. 150. -1.0 0.0 7. 11. 18. 30. 50.
+# 150. 200. -1.0 0.0 7. 11. 18. 30. 50.
+# 200. 400. 0.0 1.5 7. 11. 18. 30. 50.
+# 400. 700. 0.0 1.5 7. 11. 18. 30. 50.
+# 700. 5000. 1.5 2.5 7. 11. 18. 30. 50.
+# 5000. 15000. 1.5 2.5 7. 11. 18. 30. 50.
+#}}
+
+# ==================================================================== #
+#
+# Process specific and therefore at least partially theory-code
+# dependent settings
+#
+# ==================================================================== #
+CenterOfMassEnergy 7000. # Center-of-mass energy in GeV.
+# Some examples from the past: 200. # RHIC
+# 1800. # TeV Run I
+# 1960. # TeV Run II
+# 900. # LHC Injection Run
+# 2360. # LHC Initial Run
+# 2760. # LHC HIpp base Run
+# 5020. # LHC HIpp base Run
+# 7000. # LHC First Run (2010/2011)
+# 8000. # LHC Second Run (2012)
+# 13000. # LHC Next Run ? (2015)
+# 14000. # LHC Design Run
+PDF1 2212 # PDF of 1st hadron (following PDG convention).
+PDF2 2212 # PDF of 2nd hadron (following PDG convention).
+NPDF 2
+NPDFDim 1
+
+# ==================================================================== #
+#
+# General fastNLO settings
+#
+# ==================================================================== #
+#GlobalVerbosity DEBUG # Global output verbosity of fastNLO toolkit, def.=WARNING,
+GlobalVerbosity INFO # Global output verbosity of fastNLO toolkit, def.=WARNING,
+ # (and of user scenario code if implemented).
+ # Possible values are: DEBUG, MANUAL, INFO, WARNING, ERROR, SILENT
+OutputFilename InclusiveNJets.tab # Filename of fastNLO output table
+OutputPrecision 8 # Number of decimal digits to store in output table, def.=8.
+FlexibleScaleTable false # Create table fully flexible in mu_f (larger size), true, or
+ # table with fixed number of mu_f scale factors, def.=false.
+ScaleVariationFactors { # Factorization scale variations (only needed for fixed-scale tables)
+ 0.5 1.0 2.0 # List of scale factors must include factor '1'
+} # Scale factors will be ordered according to fastNLO convention: (1, min, ... , max)
+ReadBinningFromSteering true # Specify if binning is read from this file. MUST be true for now!
+IgnoreWarmupBinningCheck true # Ignore crosscheck of binning between steering and warmup
+ # Avoids problems with precision of bin borders e.g. in the case of Pi
+#>>fastNLO-warmup.txt # Specify alternative filename to read warm-up results from, def.=SteerFileBasename_${ScenarioName}_warmup.txt
+
+# ==================================================================== #
+#
+# Specific fastNLO settings (normally do not need to be changed)
+# Only change when you know what you are doing ...!
+#
+# ==================================================================== #
+ApplyPDFReweighting true # Apply reweighting of pdfs for an optimized interpolation, def.=true.
+CheckScaleLimitsAgainstBins true # Set limits for scale nodes to bin borders, if possible
+
+# -------------------------------------------------------------------- #
+# Choose fastNLO interpolation kernels and distance measures
+# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+# Currently implemented interpolation kernels
+# Catmull
+# Lagrange
+# OneNode
+# Linear
+#
+# Currently implemented distance measures
+# linear
+# loglog025 eq. to (log(log(4*x)))
+# log10
+# sqrtlog10 eq. to sqrt(log_10(x))
+# -------------------------------------------------------------------- #
+X_Kernel Lagrange
+X_DistanceMeasure "sqrtlog10"
+X_NNodes 12
+X_NoOfNodesPerMagnitude false
+X_NNodeCounting "NodesPerBin"
+
+Mu1_Kernel Lagrange
+Mu1_DistanceMeasure "loglog025"
+Mu1_NNodes 6
+Mu1_NNodeCounting "NodesPerBin"
+
+Mu2_Kernel Lagrange # Scale2 not used for fixed-scale tables
+Mu2_DistanceMeasure "loglog025"
+Mu2_NNodes 6
+Mu2_NNodeCounting "NodesPerBin"
diff --git a/test/flex_NodeDensity.str b/test/flex_NodeDensity.str
new file mode 100644
index 0000000000000000000000000000000000000000..a3eb18b3c8de0023877e00c609613f836887b5c9
--- /dev/null
+++ b/test/flex_NodeDensity.str
@@ -0,0 +1,207 @@
+# -*-sh-*-
+# ==================================================================== #
+#
+# A steering file for creating a fastNLO table (with NLOJet++)
+#
+# ==================================================================== #
+
+# ==================================================================== #
+#
+# Scenario specific settings
+#
+# ==================================================================== #
+
+# -------------------------------------------------------------------- #
+# Description of scenario
+# (Write each line of description in quotation marks)
+# -------------------------------------------------------------------- #
+ScenarioName InclusiveJets_flex
+ScenarioDescription {
+ "d2sigma-jet_dpT_dy_[pb_GeV]"
+ "CMS_Collaboration"
+ "Inclusive_Jet_pT"
+ "anti-kT_R=0.5"
+ "CMS-PAPER-QCD-10-011, arXiv:1106.0208, Phys. Rev. Lett. 107 (2011) 132001."
+ "RIVET_ID=CMS_2011_S9086218/D01-x01-y01,SPIRES_ID=9086218,INSPIRE_RECORD=902309"
+ "THEORY_LABEL=NLOJet++"
+ "provided by:"
+ "fastNLO_2.5.0"
+ "If you use this table, please cite:"
+ " D. Britzger, T. Kluge, K. Rabbertz, F. Stober, M. Wobisch, arXiv:1109.1310"
+}
+PublicationUnits 12 # Unit of data cross sections (negative power of 10, e.g. 12->pb, 15->fb)
+
+# -------------------------------------------------------------------- #
+# Global scenario settings
+# -------------------------------------------------------------------- #
+JetAlgo 2 # fastjet jet algorithm
+# # clustering jet algos: 0 = kT, 1 = CA, 2 = anti-kT
+# # cone jet algos: 10 = SISCone, 11 = CDFMidPointCone, 12 = D0RunIICone
+Rjet 0.5 # Jet size parameter: Required for all jets
+#OvThr 0.5 # Overlap threshold: Required for cone jets (Normally 10: 0.75; 11,12: 0.50)
+# # Overall jet phase space
+ptjmin 18. # Minimal jet pT (minimum set in interface to fastjet: 1 GeV)
+yjmin 0.0 # Minimal jet rapidity
+yjmax 3.0 # Maximal jet rapidity
+# Use either rapidity or pseudorapidity but not both! If both are set, only rapidity cuts are applied.
+#etajmin 0.0 # Minimal jet pseudorapidity
+#etajmax 3.0 # Maximal jet pseudorapidity
+#Njetmin 1 # Minimal number of jets required to be in preselected phase space
+# # (Normally one, i.e. single-jet inclusive)
+#obs0min x.y #
+#obs0max x.y # By default these values are derived from the binning in these
+#obs1min x.y # observables and are set automatically. If required the min/max
+#obs1max x.y # values for each of the maximally three observables can be set
+#obs2min x.y # by hand.
+#obs2max x.y #
+
+# -------------------------------------------------------------------- #
+# Observable definition, dimensioning, and normalization
+# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+# Specify, how many differential dimensions are measured.
+# Number must be consistent with DimensionLabels and DimensionIsDifferential, see below
+# 1 : read in 'SingleDifferentialBinning'
+# 2 : read in 'DoubleDifferentialBinning'
+# 3 : read in 'TripleDifferentialBinning'
+# -------------------------------------------------------------------- #
+LeadingOrder 2 # Number of jets for the LO process
+DifferentialDimension 2 # Dimensionality of binning
+DimensionLabels { # Labels (symbol and unit) for the measurement dimensions (from outer to inner "loop")
+ "|y|" # This defines the observables to be calculated!
+ "pT_[GeV]" # Implemented observables are: "|y|", "pT_[GeV]", "phi"
+}
+DimensionIsDifferential { # Specify for each dimension whether
+ 2 2 # 0 : the cross section is NOT differential, i.e. there are two bin borders,
+} # but NO division (normalization) by bin width
+ # 1 : the cross section is point-wise differential, i.e. only one point is given
+ # 2 : the cross section is bin-wise differential, i.e. there are two bin borders
+ # and division by bin width
+CalculateBinSize true # Calculate bin width from lower and upper bin boundaries
+BinSizeFactor 2. # Possibility to provide additional normalization factor, e.g. of 2 for bins in |y|
+#BinSize { x1 x2 x3 ... } # If 'CalculateBinSize' is 'false' provide table with bin widths 'by hand' for normalization
+ #
+ # If the calculation should not be divided by bin width, then use 'DimensionIsDifferential' equal '0'
+ # and set 'CalculateBinSize' 'true' for each dimension.
+ #
+ScaleDescriptionScale1 "pT_jet_[GeV]" # This defines the scale to be used (Note: The 1st scale must always be in units of [GeV]!)
+ScaleDescriptionScale2 "pT_max_[GeV]" # Specify 2nd scale name and unit (ONLY for flexible-scale tables)
+ # Implemented scales are: "pT_jet_[GeV]", "pT_max_[GeV]"
+
+# -------------------------------------------------------------------- #
+# Observable binning
+# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+# Use either 'SingleDifferentialBinning' or
+# 'DoubleDifferentialBinning' or
+# 'TripleDifferentialBinning'
+# in accord with 'DifferentialDimension' above
+# -------------------------------------------------------------------- #
+# Example for SingleDifferentialBinning:
+#SingleDifferentialBinning {
+# 7. 11. 18. 30. 50. 100.
+#}
+
+# Example for DoubleDifferentialBinning:
+DoubleDifferentialBinning {{
+ 1stDimLo 1stDimUp "----- Array of bin-grid for 2nd dimension -----"
+ 0.0 0.5 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 592. 638. 686. 737. 846. 1684.
+ 0.5 1.0 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 592. 638. 686. 790. 1684.
+ 1.0 1.5 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 592. 638. 686. 1410.
+ 1.5 2.0 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 1032.
+ 2.0 2.5 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 737.
+ 2.5 3.0 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 468.
+}}
+
+# Example for TripleDifferentialBinning:
+#TripleDifferentialBinning {{
+# 1stDimLo 1stDimUp 2ndDimLo 2ndDimUp "----- Array of bin-grid for 3rd dimension -----"
+# 100. 150. -1.0 0.0 7. 11. 18. 30. 50.
+# 150. 200. -1.0 0.0 7. 11. 18. 30. 50.
+# 200. 400. 0.0 1.5 7. 11. 18. 30. 50.
+# 400. 700. 0.0 1.5 7. 11. 18. 30. 50.
+# 700. 5000. 1.5 2.5 7. 11. 18. 30. 50.
+# 5000. 15000. 1.5 2.5 7. 11. 18. 30. 50.
+#}}
+
+# ==================================================================== #
+#
+# Process specific and therefore at least partially theory-code
+# dependent settings
+#
+# ==================================================================== #
+CenterOfMassEnergy 7000. # Center-of-mass energy in GeV.
+# Some examples from the past: 200. # RHIC
+# 1800. # TeV Run I
+# 1960. # TeV Run II
+# 900. # LHC Injection Run
+# 2360. # LHC Initial Run
+# 2760. # LHC HIpp base Run
+# 5020. # LHC HIpp base Run
+# 7000. # LHC First Run (2010/2011)
+# 8000. # LHC Second Run (2012)
+# 13000. # LHC Next Run ? (2015)
+# 14000. # LHC Design Run
+PDF1 2212 # PDF of 1st hadron (following PDG convention).
+PDF2 2212 # PDF of 2nd hadron (following PDG convention).
+NPDF 2
+NPDFDim 1
+
+# ==================================================================== #
+#
+# General fastNLO settings
+#
+# ==================================================================== #
+#GlobalVerbosity DEBUG # Global output verbosity of fastNLO toolkit, def.=WARNING,
+GlobalVerbosity INFO # Global output verbosity of fastNLO toolkit, def.=WARNING,
+ # (and of user scenario code if implemented).
+ # Possible values are: DEBUG, MANUAL, INFO, WARNING, ERROR, SILENT
+OutputFilename InclusiveNJets.tab # Filename of fastNLO output table
+OutputPrecision 8 # Number of decimal digits to store in output table, def.=8.
+FlexibleScaleTable true # Create table fully flexible in mu_f (larger size), true, or
+ # table with fixed number of mu_f scale factors, def.=false.
+ScaleVariationFactors { # Factorization scale variations (only needed for fixed-scale tables)
+ 0.5 1.0 2.0 # List of scale factors must include factor '1'
+} # Scale factors will be ordered according to fastNLO convention: (1, min, ... , max)
+ReadBinningFromSteering true # Specify if binning is read from this file. MUST be true for now!
+IgnoreWarmupBinningCheck true # Ignore crosscheck of binning between steering and warmup
+ # Avoids problems with precision of bin borders e.g. in the case of Pi
+#>>fastNLO-warmup.txt # Specify alternative filename to read warm-up results from, def.=SteerFileBasename_${ScenarioName}_warmup.txt
+
+# ==================================================================== #
+#
+# Specific fastNLO settings (normally do not need to be changed)
+# Only change when you know what you are doing ...!
+#
+# ==================================================================== #
+ApplyPDFReweighting true # Apply reweighting of pdfs for an optimized interpolation, def.=true.
+CheckScaleLimitsAgainstBins true # Set limits for scale nodes to bin borders, if possible
+
+# -------------------------------------------------------------------- #
+# Choose fastNLO interpolation kernels and distance measures
+# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+# Currently implemented interpolation kernels
+# Catmull
+# Lagrange
+# OneNode
+# Linear
+#
+# Currently implemented distance measures
+# linear
+# loglog025 eq. to (log(log(4*x)))
+# log10
+# sqrtlog10 eq. to sqrt(log_10(x))
+# -------------------------------------------------------------------- #
+X_Kernel Lagrange
+X_DistanceMeasure "sqrtlog10"
+X_NNodes 4
+X_NoOfNodesPerMagnitude false
+X_NNodeCounting "NodeDensity"
+
+Mu1_Kernel Lagrange
+Mu1_DistanceMeasure "loglog025"
+Mu1_NNodes 6
+Mu1_NNodeCounting "NodesPerDim2"
+
+Mu2_Kernel Lagrange # Scale2 not used for fixed-scale tables
+Mu2_DistanceMeasure "loglog025"
+Mu2_NNodes 6
+Mu2_NNodeCounting "NodesPerDim2"
diff --git a/test/flex_NodesPerBin.str b/test/flex_NodesPerBin.str
new file mode 100644
index 0000000000000000000000000000000000000000..9ec5051a957ebbc440456d0b02b1074f81e5d73a
--- /dev/null
+++ b/test/flex_NodesPerBin.str
@@ -0,0 +1,207 @@
+# -*-sh-*-
+# ==================================================================== #
+#
+# A steering file for creating a fastNLO table (with NLOJet++)
+#
+# ==================================================================== #
+
+# ==================================================================== #
+#
+# Scenario specific settings
+#
+# ==================================================================== #
+
+# -------------------------------------------------------------------- #
+# Description of scenario
+# (Write each line of description in quotation marks)
+# -------------------------------------------------------------------- #
+ScenarioName InclusiveJets_flex
+ScenarioDescription {
+ "d2sigma-jet_dpT_dy_[pb_GeV]"
+ "CMS_Collaboration"
+ "Inclusive_Jet_pT"
+ "anti-kT_R=0.5"
+ "CMS-PAPER-QCD-10-011, arXiv:1106.0208, Phys. Rev. Lett. 107 (2011) 132001."
+ "RIVET_ID=CMS_2011_S9086218/D01-x01-y01,SPIRES_ID=9086218,INSPIRE_RECORD=902309"
+ "THEORY_LABEL=NLOJet++"
+ "provided by:"
+ "fastNLO_2.5.0"
+ "If you use this table, please cite:"
+ " D. Britzger, T. Kluge, K. Rabbertz, F. Stober, M. Wobisch, arXiv:1109.1310"
+}
+PublicationUnits 12 # Unit of data cross sections (negative power of 10, e.g. 12->pb, 15->fb)
+
+# -------------------------------------------------------------------- #
+# Global scenario settings
+# -------------------------------------------------------------------- #
+JetAlgo 2 # fastjet jet algorithm
+# # clustering jet algos: 0 = kT, 1 = CA, 2 = anti-kT
+# # cone jet algos: 10 = SISCone, 11 = CDFMidPointCone, 12 = D0RunIICone
+Rjet 0.5 # Jet size parameter: Required for all jets
+#OvThr 0.5 # Overlap threshold: Required for cone jets (Normally 10: 0.75; 11,12: 0.50)
+# # Overall jet phase space
+ptjmin 18. # Minimal jet pT (minimum set in interface to fastjet: 1 GeV)
+yjmin 0.0 # Minimal jet rapidity
+yjmax 3.0 # Maximal jet rapidity
+# Use either rapidity or pseudorapidity but not both! If both are set, only rapidity cuts are applied.
+#etajmin 0.0 # Minimal jet pseudorapidity
+#etajmax 3.0 # Maximal jet pseudorapidity
+#Njetmin 1 # Minimal number of jets required to be in preselected phase space
+# # (Normally one, i.e. single-jet inclusive)
+#obs0min x.y #
+#obs0max x.y # By default these values are derived from the binning in these
+#obs1min x.y # observables and are set automatically. If required the min/max
+#obs1max x.y # values for each of the maximally three observables can be set
+#obs2min x.y # by hand.
+#obs2max x.y #
+
+# -------------------------------------------------------------------- #
+# Observable definition, dimensioning, and normalization
+# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+# Specify, how many differential dimensions are measured.
+# Number must be consistent with DimensionLabels and DimensionIsDifferential, see below
+# 1 : read in 'SingleDifferentialBinning'
+# 2 : read in 'DoubleDifferentialBinning'
+# 3 : read in 'TripleDifferentialBinning'
+# -------------------------------------------------------------------- #
+LeadingOrder 2 # Number of jets for the LO process
+DifferentialDimension 2 # Dimensionality of binning
+DimensionLabels { # Labels (symbol and unit) for the measurement dimensions (from outer to inner "loop")
+ "|y|" # This defines the observables to be calculated!
+ "pT_[GeV]" # Implemented observables are: "|y|", "pT_[GeV]", "phi"
+}
+DimensionIsDifferential { # Specify for each dimension whether
+ 2 2 # 0 : the cross section is NOT differential, i.e. there are two bin borders,
+} # but NO division (normalization) by bin width
+ # 1 : the cross section is point-wise differential, i.e. only one point is given
+ # 2 : the cross section is bin-wise differential, i.e. there are two bin borders
+ # and division by bin width
+CalculateBinSize true # Calculate bin width from lower and upper bin boundaries
+BinSizeFactor 2. # Possibility to provide additional normalization factor, e.g. of 2 for bins in |y|
+#BinSize { x1 x2 x3 ... } # If 'CalculateBinSize' is 'false' provide table with bin widths 'by hand' for normalization
+ #
+ # If the calculation should not be divided by bin width, then use 'DimensionIsDifferential' equal '0'
+ # and set 'CalculateBinSize' 'true' for each dimension.
+ #
+ScaleDescriptionScale1 "pT_jet_[GeV]" # This defines the scale to be used (Note: The 1st scale must always be in units of [GeV]!)
+ScaleDescriptionScale2 "pT_max_[GeV]" # Specify 2nd scale name and unit (ONLY for flexible-scale tables)
+ # Implemented scales are: "pT_jet_[GeV]", "pT_max_[GeV]"
+
+# -------------------------------------------------------------------- #
+# Observable binning
+# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+# Use either 'SingleDifferentialBinning' or
+# 'DoubleDifferentialBinning' or
+# 'TripleDifferentialBinning'
+# in accord with 'DifferentialDimension' above
+# -------------------------------------------------------------------- #
+# Example for SingleDifferentialBinning:
+#SingleDifferentialBinning {
+# 7. 11. 18. 30. 50. 100.
+#}
+
+# Example for DoubleDifferentialBinning:
+DoubleDifferentialBinning {{
+ 1stDimLo 1stDimUp "----- Array of bin-grid for 2nd dimension -----"
+ 0.0 0.5 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 592. 638. 686. 737. 846. 1684.
+ 0.5 1.0 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 592. 638. 686. 790. 1684.
+ 1.0 1.5 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 592. 638. 686. 1410.
+ 1.5 2.0 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 468. 507. 548. 1032.
+ 2.0 2.5 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 330. 362. 395. 430. 737.
+ 2.5 3.0 18. 21. 24. 28. 32. 37. 43. 49. 56. 64. 74. 84. 97. 114. 133. 153. 174. 196. 220. 245. 272. 300. 468.
+}}
+
+# Example for TripleDifferentialBinning:
+#TripleDifferentialBinning {{
+# 1stDimLo 1stDimUp 2ndDimLo 2ndDimUp "----- Array of bin-grid for 3rd dimension -----"
+# 100. 150. -1.0 0.0 7. 11. 18. 30. 50.
+# 150. 200. -1.0 0.0 7. 11. 18. 30. 50.
+# 200. 400. 0.0 1.5 7. 11. 18. 30. 50.
+# 400. 700. 0.0 1.5 7. 11. 18. 30. 50.
+# 700. 5000. 1.5 2.5 7. 11. 18. 30. 50.
+# 5000. 15000. 1.5 2.5 7. 11. 18. 30. 50.
+#}}
+
+# ==================================================================== #
+#
+# Process specific and therefore at least partially theory-code
+# dependent settings
+#
+# ==================================================================== #
+CenterOfMassEnergy 7000. # Center-of-mass energy in GeV.
+# Some examples from the past: 200. # RHIC
+# 1800. # TeV Run I
+# 1960. # TeV Run II
+# 900. # LHC Injection Run
+# 2360. # LHC Initial Run
+# 2760. # LHC HIpp base Run
+# 5020. # LHC HIpp base Run
+# 7000. # LHC First Run (2010/2011)
+# 8000. # LHC Second Run (2012)
+# 13000. # LHC Next Run ? (2015)
+# 14000. # LHC Design Run
+PDF1 2212 # PDF of 1st hadron (following PDG convention).
+PDF2 2212 # PDF of 2nd hadron (following PDG convention).
+NPDF 2
+NPDFDim 1
+
+# ==================================================================== #
+#
+# General fastNLO settings
+#
+# ==================================================================== #
+#GlobalVerbosity DEBUG # Global output verbosity of fastNLO toolkit, def.=WARNING,
+GlobalVerbosity INFO # Global output verbosity of fastNLO toolkit, def.=WARNING,
+ # (and of user scenario code if implemented).
+ # Possible values are: DEBUG, MANUAL, INFO, WARNING, ERROR, SILENT
+OutputFilename InclusiveNJets.tab # Filename of fastNLO output table
+OutputPrecision 8 # Number of decimal digits to store in output table, def.=8.
+FlexibleScaleTable true # Create table fully flexible in mu_f (larger size), true, or
+ # table with fixed number of mu_f scale factors, def.=false.
+ScaleVariationFactors { # Factorization scale variations (only needed for fixed-scale tables)
+ 0.5 1.0 2.0 # List of scale factors must include factor '1'
+} # Scale factors will be ordered according to fastNLO convention: (1, min, ... , max)
+ReadBinningFromSteering true # Specify if binning is read from this file. MUST be true for now!
+IgnoreWarmupBinningCheck true # Ignore crosscheck of binning between steering and warmup
+ # Avoids problems with precision of bin borders e.g. in the case of Pi
+#>>fastNLO-warmup.txt # Specify alternative filename to read warm-up results from, def.=SteerFileBasename_${ScenarioName}_warmup.txt
+
+# ==================================================================== #
+#
+# Specific fastNLO settings (normally do not need to be changed)
+# Only change when you know what you are doing ...!
+#
+# ==================================================================== #
+ApplyPDFReweighting true # Apply reweighting of pdfs for an optimized interpolation, def.=true.
+CheckScaleLimitsAgainstBins true # Set limits for scale nodes to bin borders, if possible
+
+# -------------------------------------------------------------------- #
+# Choose fastNLO interpolation kernels and distance measures
+# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+# Currently implemented interpolation kernels
+# Catmull
+# Lagrange
+# OneNode
+# Linear
+#
+# Currently implemented distance measures
+# linear
+# loglog025 eq. to (log(log(4*x)))
+# log10
+# sqrtlog10 eq. to sqrt(log_10(x))
+# -------------------------------------------------------------------- #
+X_Kernel Lagrange
+X_DistanceMeasure "sqrtlog10"
+X_NNodes 12
+X_NoOfNodesPerMagnitude false
+X_NNodeCounting "NodesPerBin"
+
+Mu1_Kernel Lagrange
+Mu1_DistanceMeasure "loglog025"
+Mu1_NNodes 6
+Mu1_NNodeCounting "NodesPerBin"
+
+Mu2_Kernel Lagrange # Scale2 not used for fixed-scale tables
+Mu2_DistanceMeasure "loglog025"
+Mu2_NNodes 6
+Mu2_NNodeCounting "NodesPerBin"