Bug fix: mergefkswithplot counts up in particle numbers not down

michel-decay/demo-observable/demo-observable.py

@@ -11,7 +11,7 @@ setup(folder='with-cut/out.tar.bz2', flavour='mu-e')
 
 
 loC = scaleset(mergefks(sigma('m2enn0')), 1/gamma0)
-fig, ans = mergefkswithplot([[sigma('m2ennR')], [sigma('m2ennF')]])
+fig, ans = mergefkswithplot([[sigma('m2ennF')], [sigma('m2ennR')]])
 nloC = scaleset(ans, alpha/gamma0)
 nnloC = scaleset(mergefks(
     sigma('m2ennRF'),

mu-e-scattering/muone/muone.py

@@ -64,13 +64,13 @@ aconloEM = scaleset(mergefks(
 
 # Here, we create $\xi_c$ plot
 fignlo, aconlo = mergefkswithplot([
+    [
+        sigma('em2emFEE'), sigma('em2emFMM'), sigma('em2emFEM')
+    ],
     [
         sigma('em2emREE15'), sigma('em2emREE35'),
         sigma('em2emRMM'),
         sigma('em2emREM')
-    ],
-    [
-        sigma('em2emFEE'), sigma('em2emFMM'), sigma('em2emFEM')
     ]
 ], scale=conv*alpha**3)
 
@@ -84,9 +84,9 @@ savefig("plots/xi-nlo.pdf")
 ### Load NNLO{{{
 # Make $\xi_c$ plot
 fignnlo, aconnlo = mergefkswithplot([
-    [sigma('em2emRREE1'), sigma('em2emRREE3')],
+    [sigma('em2emFF', folder='ff.tar.bz2')],
     [sigma('em2emRFEE15'), sigma('em2emRFEE35')],
-    [sigma('em2emFF', folder='ff.tar.bz2')]
+    [sigma('em2emRREE1'), sigma('em2emRREE3')]
 ], conv*alpha**4)
 
 mulify(fignnlo)