pyhf xml2json|pyhf cls / pyhf cls wspace.json

[lukasheinrich]

Description

subject to #203

but this would be nice

import pyhf.utils
import json
d = json.load(open('./test.json'))
from pyhf import Model
p = Model({'channels':d['channels']}, poiname='SigXsecOverSM')
result = pyhf.utils.runOnePoint(1.0, d['data']['channel1'] + p.config.auxdata, p)
json.dumps({'CLs_obs': result[-2].tolist()[0], 'CLs_exp': result[-1].ravel().tolist()})

returns

'{"CLs_obs": 0.5250899403989365, "CLs_exp": [0.005398954300520573, 0.02413688765551562, 0.0955643537176678, 0.30007471275202563, 0.6452061679174199]}'

once workspaces are on hepdata we could do

curl http://hepdata.net/ins12345/workspace.json | pyhf cls

and to reinterpret for RECAST

curl http://hepdata.net/ins12345/workspace.json | pyhf cls --patch newsignal.json