Distribute resummino run files to clusters

After running this example one only needs to run the generated output/*.sh files or queue them to your favourite cluster:

$ ./output/*.sh

[1]:

import hepi
print(hepi.__version__)
import smpl
import numpy as np
from hepi.run import resummino as rs
import hepi.util as util
import matplotlib.pyplot as plt
rs.set_path("resummino")
print (rs.get_path())

0.2.9.7
~/git/resummino

Single point

[2]:

params = [
    "slha.in",
]
pss = [
    (1000023,1000024),
    (1000024,-1000024),
    (1000023,1000022),
    (1000002,1000022),
     ]
for pa,pb in pss:
    for param in params:
        i = hepi.Input(hepi.Order.NLO_PLUS_NLL,13000,pa,pb,param,"PDF4LHC21_mc","PDF4LHC21_mc",1., 1.,id="0")
        #i = hepi.Input(hepi.Order.NLO,13000,pa,pb,param,"CT14lo","CT14lo",1., 1.,model_path=model_path,id="5")
        li = [i]
        rs_dl = rs.run(li,noskip=False,run=False,parse=False)

Skipped: 0 Not skipped: 1

Skipped: 0 Not skipped: 1

Skipped: 0 Not skipped: 1

Skipped: 0 Not skipped: 1

Mass scan

[3]:

params = [
    "slha.in",
]
pss = [
    (1000011,-1000011),
     ]
for pa,pb in pss:
    for param in params:
        i = hepi.Input(hepi.Order.NLO_PLUS_NLL,13000,pa,pb,param,"PDF4LHC21_mc","PDF4LHC21_mc",1., 1.,id="0")
        #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(1000,2000,9))
        rs_dl = rs.run(li,noskip=False,run=False,parse=False)

Skipped: 0 Not skipped: 9

Scale scan

[4]:

params = [
    "wino.in",
    "hino.in",
]
pss = [
    (1000022,1000022),
     ]
for pa,pb in pss:
    for param in params:
        i = hepi.Input(hepi.Order.NLO_PLUS_NLL,13000,pa,pb,param,"PDF4LHC21_mc","PDF4LHC21_mc",1., 1.,id="0")
        #i = hepi.Input(hepi.Order.NLO,13000,pa,pb,param,"CT14lo","CT14lo",1., 1.,model_path=model_path,id="5")
        li = [i]
        li = hepi.seven_point_scan(li)
        rs_dl = rs.run(li,noskip=False,run=False,parse=False)

Skipped: 0 Not skipped: 7

Skipped: 0 Not skipped: 7

PDF scan

[5]:

params = [
    "slha.in",
]
pss = [
    (1000011,-1000011),
     ]
for pa,pb in pss:
    for param in params:
        i = hepi.Input(hepi.Order.NLO_PLUS_NLL,13600,pa,pb,param,"PDF4LHC21_mc","PDF4LHC21_mc",1., 1.,id="0")
        #i = hepi.Input(hepi.Order.NLO,13000,pa,pb,param,"CT14lo","CT14lo",1., 1.,model_path=model_path,id="5")
        li = [i]
        li = hepi.pdf_scan(li)
        rs_dl = rs.run(li,noskip=False,run=False,parse=False)

Skipped: 0 Not skipped: 101

[ ]: