-
Notifications
You must be signed in to change notification settings - Fork 85
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
1 parent
966742c
commit 139b3fd
Showing
27 changed files
with
1,349,505 additions
and
6,123 deletions.
There are no files selected for viewing
Binary file not shown.
Large diffs are not rendered by default.
Oops, something went wrong.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| 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 | ||
| ``` | ||
|
|
Large diffs are not rendered by default.
Oops, something went wrong.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| 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" | ||
| } |
Oops, something went wrong.