EditWorkspace

Edit a workspace using the RooWorkspace.factory function or the editws functionality of RooFitUtils (requires the corresponding library to be available).

Usage:

+EditWorkspaces.mergePOIs {
  +HWWRun2GGF {
    # merge some parameters into one
    <merge.muGGF = {"mu_gg2H_*"}>
    <merge.muVBF = {"mu_VBF_*","mu_*_qq2Hqq_*"}>
    <merge.muVH = {"mu_WH_qq2Hlnu_*","mu_ZH_qqHll_*","mu_ZH_ggHll_*"}>
    # add some POIs
    <addPOIs = {"muGGF"}>
    # set some parameters constant
    <constPars = {"muVBF","muVH"}>
  }
}

Supported edits include:

• merge: merge a set of parameters to one
• lines: lines of string commands that will be forwarded to the RooFit::RooWorkspace.factory function verbatim
• files: files of string commands that will be forwarded to the RooFit::RooWorkspace.factory function verbatim The edits will be applied to the workspace, after which a new, clean workspace will be created, replacing the old workspace of the same name. The flag useRooFitUtils=true can be used to change the processing engine away from the RooFit::RooWorkspace.factory to the more advanced editws engine that is provided within the RooFitUtils extension package. If that package is not present, this flag will trigger a runtime error.

Unfold

Implement Unfolding in a workspace.

Usage:

+Unfold {
  + Unfolded {
    <SystematicsThreshold=0.005>
    <MigrationThreshold=1e-6>
    <mu=false> # decide if you want your POI to be the signal strength (true, default) or the cross-section (false)
    + Measured {
      <workspace="Reco">
    }
    + Migrations {
      <workspace="Truth">
      <region="SR_2D">
    }
    + Truth {
      <workspace="Truth">
      <region="FR_1D">
    }
    + Reconstructed {
      <workspace="Truth">
      <region="SR">
      <obs="obs_x_SR_1D">
    }   
    # option: regularize
    +Regularization {
      <method="curvature"> # Most regularization terms use the curvature of some reference quantity
      <reference="correctionfactor"> 
      # <reference="POI">  # just use the curvature of whatever POI you chose
      # <reference="signalstrength">  # This regularization method corresponds to SVD unfolding
      # <reference="excess">  # This regularization method corresponds to the default settings of TUnfold
      <tau=0.3>
      }
  }
}