pyslha comparison

import hepi
print(hepi.__version__)
import smpl
import numpy as np
import hepi.resummino as rs
import hepi.util as util
import matplotlib.pyplot as plt
rs.set_path("~/git/resummino_release")
print (rs.get_path())

0.1.8.9+dirty
~/git/resummino_release

params = [
    "mastercode_with_gm2.in",
]
pss = [
      (1000022,1000022),
     ]

for pa,pb in pss:
    for param in params:
        i = hepi.Input(hepi.Order.NLO,13000,pa,pb,param,"cteq6l1","cteq66",1., 1.,id="5")
        #i = hepi.Input(hepi.Order.NLO,13000,pa,pb,param,"CT14lo","CT14lo",1., 1.,model_path=model_path,id="5")
        li = [i]
        li = hepi.mass_scan([i],pa, np.linspace(2.48758955E+02,2.48758955E+02,1))
        rs_dl = rs.run(li,noskip=False)
rs_dl["NLO"]

Running: 1 jobs

0    (9.721+/-0.008)e-06
Name: NLO, dtype: object

for pa,pb in pss:
    for param in params:
        i = hepi.Input(hepi.Order.NLO,13000,pa,pb,param,"cteq6l1","cteq66",1., 1.,id="5")
        #i = hepi.Input(hepi.Order.NLO,13000,pa,pb,param,"CT14lo","CT14lo",1., 1.,model_path=model_path,id="5")
        li = [i]
        #li = hepi.mass_scan([i],pa, np.linspace(2.4876e2,2.4876e2,1))
        rs_dl = rs.run(li,noskip=False)
rs_dl["NLO"]

Running: 1 jobs

0    (9.721+/-0.008)e-06
Name: NLO, dtype: object