Produce f and g plot

michel-decay/f-and-g/f-and-g.py

@@ -416,3 +416,98 @@ GFULL = addplots(GYFS, g3LLsub)
 
 ###########################################################}}}
 ##########################################################################}}}
+## Make plot{{{
+### $f$ plot{{{
+ffig, (ax1, ax2, ax3) = twopanel(
+    labelx="$E_e$",
+    upleft=[flo, FFULL],
+    colupleft=['C2', 'C1'],
+    labupleft="$F$",
+    downleft=[divideplots(fnlo, flo, offset=1)],
+    labdownleft={
+        'ylabel': "$K^{(1)}$",
+        'color': 'C0'
+    },
+    coldownleft=['C0'],
+    downright=[
+        divideplots(fnnlo, flo, offset=1),
+        divideplots(
+            addplots(fLLYFS, f3LLsub),
+            flo,
+            offset=1
+        )
+    ],
+    labdownright={
+        'ylabel': "$K^{(2,\\text{resum})}\\equiv f_i/f_0$",
+        'color': 'C3'
+    },
+    coldownright=['C3', 'C1']
+)
+ax3.tick_params(labelcolor='C1', axis='y')
+
+ax2.set_ylim(0.9, 1.52)
+ax3.set_ylim(0.999, 1.0052)
+ffig.legend(
+    [
+        matplotlib.lines.Line2D([0], [0], color=i) for i in [
+            'C2', 'C0', 'C3', 'C1'
+        ]
+    ],
+    [
+        "$f_0$", "$f_1$", "$f_2$", "$f_\\text{resum}$"
+    ],
+    loc="upper center", ncol=4
+)
+mulify(ffig)
+
+
+ffig.savefig("fplot.pdf")
+###########################################################}}}
+### $g$ plot{{{
+gfig, (ax1, ax2, ax3) = twopanel(
+    labelx="$E_e$",
+    upleft=[glo, GFULL],
+    colupleft=['C2', 'C1'],
+    labupleft="$G$",
+    downleft=[divideplots(gnlo, glo, offset=1)],
+    labdownleft={
+        'ylabel': "$K^{(1)}$",
+        'color': 'C0'
+    },
+    coldownleft=['C0'],
+    downright=[
+        divideplots(gnnlo, glo, offset=1),
+        divideplots(
+            addplots(gLLYFS, g3LLsub),
+            glo,
+            offset=1
+        )
+    ],
+    labdownright={
+        'ylabel': "$K^{(2,\\text{resum})}\\equiv g_i/g_0$",
+        'color': 'C3'
+    },
+    coldownright=['C3', 'C1']
+)
+ax3.tick_params(labelcolor='C1', axis='y')
+
+ax2.set_ylim(0.58, 1.52)
+ax3.set_ylim(0.996, 1.001)
+gfig.legend(
+    [
+        matplotlib.lines.Line2D([0], [0], color=i) for i in [
+            'C2', 'C0', 'C3', 'C1'
+        ]
+    ],
+    [
+        "$g_0$", "$g_1$", "$g_2$", "$g_\\text{resum}$"
+    ],
+    loc="upper center", ncol=4
+)
+mulify(gfig)
+
+
+gfig.savefig("gplot.pdf")
+###########################################################}}}
+##########################################################################}}}
+##########################################################################}}}

michel-decay/f-and-g/makefile

+name=f-and-g
+figures=fplot.pdf gplot.pdf
+include ../../tools/makefile.conf

tools/makefile.conf

@@ -3,11 +3,11 @@ SELF_DIR := $(dir $(lastword $(MAKEFILE_LIST)))
 UNAME_S := $(shell uname -s)
 ifeq ($(UNAME_S),Linux)
     ifneq ("$(wildcard /etc/redhat-release)","")
-        PYTHON=PYTHONPATH=$$PYTHONPATH:/scratch/ulrich/pymule/tools/ /afs/psi.ch/sys/psi.x86_64_slp6/Programming/psi-python27/4.4.0/bin/python
+        PYTHON=PYTHONPATH=$$PYTHONPATH:/scratch/ulrich/pymule/tools/ /afs/psi.ch/sys/psi.x86_64_slp6/Programming/psi-python27/4.4.0/bin/python -B
         PEP8=$(PYTHON) /afs/psi.ch/sys/psi.x86_64_slp6/Programming/psi-python27/4.4.0/bin/pep8
         JUPYTERCONV=$(PYTHON) /afs/psi.ch/sys/psi.x86_64_slp6/Programming/psi-python27/4.4.0/bin/jupyter-nbconvert
     else
-        PYTHON=python
+        PYTHON=python -B
         PEP8=pep8
         JUPYTERCONV=jupyter nbconvert
     endif

tools/py2nb.py

@@ -92,7 +92,9 @@ def executenb(nb, dir):
         pymule + "\nget_ipython().magic('matplotlib inline')"
     )
 
-    ep = nbconvert.preprocessors.ExecutePreprocessor()
+    ep = nbconvert.preprocessors.ExecutePreprocessor(
+        timeout=600
+    )
     ep.preprocess(nb, {'metadata': {'path': dir}})
 
     nb.cells[i]['source'] = old