1 | | blabla |
| 1 | == HistFactory Input to calculate Combined SR exclusion limits == |
| 2 | |
| 3 | HistFactory JSON files can be attached to the analysis to be able to estimate exclusion limits of combined likelihood profiles and even global profile likelihoods. JSON files need to be under the same folder with info file and required information has to be added into the info file as shown below. |
| 4 | |
| 5 | {{{ |
| 6 | <analysis id="atlas_xyz_00_00"> |
| 7 | <lumi>139.0</lumi> |
| 8 | <pyhf id="<Profile Name>"> |
| 9 | <name>FILE_NAME.json</name> |
| 10 | <regions> |
| 11 | <channel name="<Channel Name1>"> |
| 12 | SRX |
| 13 | SRY |
| 14 | </channel> |
| 15 | <channel name="<Channel Name2>"></channel> |
| 16 | <region> |
| 17 | </pyhf> |
| 18 | </analysis> |
| 19 | }}} |
| 20 | |
| 21 | |
| 22 | where "<Profile Name>" is the identifier of the profile and will be printed in the output file to show the exclusion estimates calculated using this specific likelihood profile. FILE_NAME.json is the name of the JSON file. "<Channel Name>" is the name of the channel as specified in the JSON file, please note that in case of the wrong declaration the profile will be ignored. In the example above, channel 1 has two signal regions declared as SRX and SRY. '''These names are needed to be taken from the PAD analysis and have to be ordered as given in JSON file. In case of wrong ordering, the calculation will be wrong. To make sure, please refer to the analysis description.''' If further help needed, one can use [http://madanalysis.irmp.ucl.ac.be/attachment/wiki/SRCombinations/write_histfactory_info.py write_pyhf_info.py] |
| 23 | |
| 24 | {{{ |
| 25 | python write_histfactory_info.py -i FILE1.json FILE2.json FILE3.json |
| 26 | }}} |
| 27 | where -i refers to the interactive mode, which writes the file by giving you directions. The JSON files has to be named as in info file and placed in the same folder wiith the info file (''~/madanalysis5/tools/<PADofChoice>/Build/SampleAnalyzer/User/Analysis/atlas_xyz_00_00.info'') |
| 28 | |
| 29 | In order to use HistFactory, one needs to install ''pyhf'' package which is automatically installed via following command |
| 30 | |
| 31 | {{{ |
| 32 | install pyhf |
| 33 | }}} |
| 34 | |
| 35 | After these steps all given signal region combinations are automatically calculated and additionally MadAnalysis constructs a global likelihood profile to combine all given HistFactory files with the same parameter of interest. |