December Updates

The main changes in this update are improvements to the combineTool.py script, new programs to recreate datacards for run 1 analyses utilising the RooMorphingPdf, and the development of the MSSM physics model to support the corresponding signal processes. More details are given below:

New features

• Set of programs to recreate or modify run 1 datacards for BSM Higgs analyses. These have been added to a new Run1BSMComb sub-package.
  • MorphingMSSMUpdate.cpp: MSSM H->tautau update analysis
  • MorphingMSSMLegacy.cpp: MSSM H->tautau legacy publication (mu-tau only, just for testing)
  • MorphingHhh.cpp and MorphingAZh.cpp: H->hh and A->Zh analyses from HIG-14-034
  • AdaptHbb.cpp: modifies existing MSSM H->bb cards to conform to standard process naming and normalisation conventions
  • AdaptChargedHiggs.cpp: modifies existing H+->tau nu cards to conform to standard process naming and normalisation conventions
  • Non-morphing versions of some of the cards can also be created with separate programs
• The MSSM.py physics model is now largely complete for building workspaces for the Run 1 analyses. The model files to use are specified as command line options and multiple data-taking periods can be combined. Theory uncertainties, as a function of mA and tan(beta), are also created within the physics model.
• Example plotting scripts for the combine/combineTool.py output:
  • plotMultiDimFit.py: plotting 2D likelihood scans
  • plotLimitGrid.py: for the asymptotic limits calculated for a 2D grid of values using either the AsymptoticGrid or HybridNewGrid methods implemented in combineTool.py
  • limitCompare.py: generic limit plotting and comparison script
• New program PostFitShapesFromWorkspace.cpp with a similar interface to PostFitShapes.cpp but evaluates the shapes using the pdfs in the combine workspace instead of the datacard inputs.
• Preliminary datacard configuration for the first run 2 MSSM analysis in dedicated sub-package MSSM_13TeV
• New sub-package HIG15002 containing scripts used to produce LHC Higgs combination results
• Updated documentation covering the construction of RooMorphingPdfs in CombineHarvester and the reproduction of legacy datacards using the programs listed above.
• ch::Parameter values can now be frozen such that setting a new value does not actually change the value that is stored. This is useful for doing the uncertainty propagation from a covariance matrix restricted to the effect of a subset of parameters.
• Some parts of the ch::CombineHarvester behaviour is now controlled by flags that can be set with the SetFlag/GetFlag methods. At the moment there is no change to the existing behaviour and so the possible flags are not yet documented.
• Significant development of the HybridNewGrid combineTool.py module. Implemented criteria that terminate the toy generation automatically at each point.

Bug fixes:

• plotImpacts.py: the POI uncertainty drawn at the top of the frame was incorrect, giving roughly half of the true value. Note that this is a plotting issue only, the value calculated by the combineTool.py Impacts method is correct.
• BuildRooMorphing: Fixed a couple of issues in the way the FastVerticalInterpHistPdf2 pdfs were being created that was out of sync with the way it's done in text2workspace. This affects datacards with shape systematics having a scale factor != 1.0.