13: Added some useful constants

pymule/__init__.py

@@ -10,3 +10,5 @@ from loader import importreg, pattern, setup, sigma,                          \
                    addsets, scaleset
 
 from plot import errorband
+
+from constants import pi, alpha, GF, conv, Mmu, Mel, Mtau

pymule/constants.py

+import numpy as np
+import os
+
+
+def read_globaldef(globaldef=""):
+    if len(globaldef) == 0:
+        dir_path = os.path.dirname(os.path.realpath(__file__))
+        globaldef = os.path.abspath(os.path.join(
+            dir_path, os.pardir, os.pardir, 'global_def.f95'
+        ))
+
+    if not os.path.exists(globaldef):
+        raise IOError("No such file or directory %s" % globaldef)
+
+    with open(globaldef) as fp:
+        return [
+            (n,eval(v))
+            for n,v in re.findall(
+                "parameter *:: *([A-Z][A-Z\d]*) = ([\dE\.+/-]+)",
+                fp.read(),
+                re.IGNORECASE
+            )
+        ]
+
+
+pi = np.pi
+
+Mmu = 105.6583715  # MeV
+Mel = 0.510998928  # MeV
+Mtau = 1776.82     # MeV
+
+alpha = 1./137.0359997
+GF = 1.16637E-11   # MeV^-2
+
+conv = 3.893793655569159807686738742e8  # MeV^2 mu barn