Improve plotting script for kfactor with options for scale choice

tools/plotting/fastnnlo_kfactor.py

@@ -44,6 +44,14 @@ def main():
                         help='Output filename (optional).')
     parser.add_argument('-r', '--reftab', default=None, type=str,
                         help='Reference table to compare to (optional).')
+    parser.add_argument('-mur', default=0, type=int,
+                        help='mur index numerator, default is 0.')
+    parser.add_argument('-muf', default=0, type=int,
+                        help='muf index numerator, default is 0.')
+    parser.add_argument('-murref', default=0, type=int,
+                        help='mur index reference, default is 0.')
+    parser.add_argument('-mufref', default=0, type=int,
+                        help='muf index reference, default is 0.')
 
     # Parse arguemnts
     args = vars(parser.parse_args())
@@ -154,9 +162,6 @@ def main():
     xs_ref    = []
     dxsu_ref  = []
     dxsl_ref  = []
-    # This should be an option ...
-    fnlo.SetMuRFunctionalForm(fastnlo.kScale1)
-    fnlo.SetMuFFunctionalForm(fastnlo.kScale1)
 
     for n in orders:
         for j in range(0, 3):
@@ -165,6 +170,8 @@ def main():
                 fref.SetContributionON(fastnlo.kFixedOrder, j, orders[n][j])
         print('\n')
         print('Calc XS for order: %s' % n, '\n')
+        fnlo.SetMuRFunctionalForm(args['mur'])
+        fnlo.SetMuFFunctionalForm(args['muf'])
         fnlo.CalcCrossSection()
         xsunc = np.array(fnlo.GetAddUncertaintyVec(fastnlo.kAddStat))
         xs  = xsunc[0,:]
@@ -179,6 +186,8 @@ def main():
         if reftab:
             print('\n')
             print('Ref XS for order: %s' % n, '\n')
+            fnlo.SetMuRFunctionalForm(args['murref'])
+            fnlo.SetMuFFunctionalForm(args['mufref'])
             fref.CalcCrossSection()
             xsunc = np.array(fref.GetAddUncertaintyVec(fastnlo.kAddStat))
             xs  = xsunc[0,:]