Skip to content

Commit

Permalink
squash
Browse files Browse the repository at this point in the history
  • Loading branch information
@lukasheinrich
lukasheinrich committed Mar 19, 2021
1 parent 5c05e6e commit fa514ba
Show file tree
Hide file tree
Showing 29 changed files with 1,349,505 additions and 6,168 deletions.
Binary file added 1Lbb-likelihoods.tar.gz
View file
Binary file not shown.
105,998 changes: 105,998 additions & 0 deletions 1Lbb-likelihoods/BkgOnly.json
View file

Large diffs are not rendered by default.

26 changes: 26 additions & 0 deletions 1Lbb-likelihoods/README.md
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,26 @@
# JSON Likelihoods for 1Lbb Analysis

The JSON likelihoods are serialized in this folder. This is done by providing a background-only workspace containing the signal/control channels at `BkgOnly.json` as well as patch files for each mass point on the signal phase-space explored in the analysis.

All [jsonpatches](http://jsonpatch.com/) are contained in the file `patchset.json`. Each patch is identified in `patchset.json` by the metadata field `"name": "C1N2_Wh_hbb_[m1]_[m2]"` where `m1` is the mass of both the lightest chargino and the next-to-lightest neutralino (which are assumed to be nearly mass degenerate), and `m2` is the mass of the lightest neutralino.

## Producing signal workspaces

As an example, we use [python jsonpatch](https://python-json-patch.readthedocs.io/en/latest/) to make the full json likelihood workspace for the signal point `C1N2_Wh_hbb_700_400`:

```
jsonpatch BkgOnly.json <(pyhf patchset extract patchset.json --name "C1N2_Wh_hbb_700_400") > C1N2_Wh_hbb_700_400.json
```

## Computing signal workspaces

For example, with [pyhf](https://scikit-hep.org/pyhf/), you can do any of the following:

```
pyhf cls BkgOnly.json --patch <(pyhf patchset extract patchset.json --name "C1N2_Wh_hbb_700_400")
jsonpatch BkgOnly.json <(pyhf patchset extract patchset.json --name "C1N2_Wh_hbb_700_400") | pyhf cls
pyhf cls C1N2_Wh_hbb_700_400.json
```

1,242,030 changes: 1,242,030 additions & 0 deletions 1Lbb-likelihoods/patchset.json
View file

Large diffs are not rendered by default.

275 changes: 275 additions & 0 deletions aaa
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,275 @@
{
"channels": [
{
"name": "channel1",
"samples": [
{
"data": [
20.0,
10.0
],
"modifiers": [
{
"data": {
"hi": 1.05,
"lo": 0.95
},
"name": "syst1",
"type": "normsys"
},
{
"data": null,
"name": "SigXsecOverSM",
"type": "normfactor"
}
],
"name": "signal"
},
{
"data": [
100.0,
0.0
],
"modifiers": [
{
"data": null,
"name": "lumi",
"type": "lumi"
},
{
"data": [
5.000000074505806,
0.0
],
"name": "staterror_channel1",
"type": "staterror"
},
{
"data": {
"hi": 1.05,
"lo": 0.95
},
"name": "syst2",
"type": "normsys"
}
],
"name": "background1"
},
{
"data": [
0.0,
100.0
],
"modifiers": [
{
"data": null,
"name": "lumi",
"type": "lumi"
},
{
"data": [
0.0,
10.0
],
"name": "staterror_channel1",
"type": "staterror"
},
{
"data": {
"hi": 1.05,
"lo": 0.95
},
"name": "syst3",
"type": "normsys"
}
],
"name": "background2"
}
]
}
],
"measurements": [
{
"config": {
"parameters": [
{
"auxdata": [
1.0
],
"bounds": [
[
0.5,
1.5
]
],
"fixed": true,
"inits": [
1.0
],
"name": "lumi",
"sigmas": [
0.1
]
},
{
"bounds": [
[
0.0,
3.0
]
],
"inits": [
1.0
],
"name": "SigXsecOverSM"
},
{
"fixed": true,
"name": "syst1"
}
],
"poi": "SigXsecOverSM"
},
"name": "GaussExample"
},
{
"config": {
"parameters": [
{
"auxdata": [
1.0
],
"bounds": [
[
0.5,
1.5
]
],
"fixed": true,
"inits": [
1.0
],
"name": "lumi",
"sigmas": [
0.1
]
},
{
"bounds": [
[
0.0,
3.0
]
],
"inits": [
1.0
],
"name": "SigXsecOverSM"
},
{
"fixed": true,
"name": "syst1"
}
],
"poi": "SigXsecOverSM"
},
"name": "GammaExample"
},
{
"config": {
"parameters": [
{
"auxdata": [
1.0
],
"bounds": [
[
0.5,
1.5
]
],
"fixed": true,
"inits": [
1.0
],
"name": "lumi",
"sigmas": [
0.1
]
},
{
"bounds": [
[
0.0,
3.0
]
],
"inits": [
1.0
],
"name": "SigXsecOverSM"
},
{
"fixed": true,
"name": "syst1"
}
],
"poi": "SigXsecOverSM"
},
"name": "LogNormExample"
},
{
"config": {
"parameters": [
{
"auxdata": [
1.0
],
"bounds": [
[
0.5,
1.5
]
],
"fixed": true,
"inits": [
1.0
],
"name": "lumi",
"sigmas": [
0.1
]
},
{
"bounds": [
[
0.0,
3.0
]
],
"inits": [
1.0
],
"name": "SigXsecOverSM"
},
{
"fixed": true,
"name": "syst1"
}
],
"poi": "SigXsecOverSM"
},
"name": "ConstExample"
}
],
"observations": [
{
"data": [
122.0,
112.0
],
"name": "channel1"
}
],
"version": "1.0.0"
}
27 changes: 0 additions & 27 deletions docs/_templates/modifierclass.rst
View file
Original file line number Diff line number Diff line change
@@ -1,30 +1,3 @@
:github_url: https://github.com/scikit-hep/pyhf/blob/master/{{module | replace(".", "/") }}

{{ name | escape | underline}}

.. currentmodule:: {{ module }}

.. autoclass:: {{ name }}
:show-inheritance:

{% block attributes %}
{% if attributes %}
.. rubric:: Attributes

{% for item in attributes %}
.. autoattribute:: {{ name }}.{{ item }}
{%- endfor %}
{% endif %}
{% endblock %}

{% block methods %}
{% if methods %}
.. rubric:: Methods

{% for item in methods %}
{% if item not in inherited_members %}
.. automethod:: {{ name }}.{{ item }}
{% endif %}
{%- endfor %}
{% endif %}
{% endblock %}
17 changes: 0 additions & 17 deletions docs/api.rst
View file
Original file line number Diff line number Diff line change
Expand Up @@ -81,23 +81,6 @@ Optimizers
opt_scipy.scipy_optimizer
opt_minuit.minuit_optimizer

Modifiers
---------

.. currentmodule:: pyhf.modifiers

.. autosummary::
:toctree: _generated/
:nosignatures:
:template: modifierclass.rst

histosys
normfactor
normsys
shapefactor
shapesys
staterror

Interpolators
-------------

Expand Down
Loading

0 comments on commit fa514ba

Please sign in to comment.