Merge branch 'main' of https://github.com/rkansal47/HHbbVV

README.md

@@ -422,6 +422,17 @@ git clone -b v2.0.0 https://github.com/cms-analysis/CombineHarvester.git Combine
 scramv1 b clean; scramv1 b
 ```
 
+### Packages
+
+To create datacards, you need to use the same cmsenv as above + these packages:
+
+```bash
+pip3 install --upgrade pip
+pip3 install rhalphalib
+cd /path/to/your/local/HHbbVV/repo
+pip3 install -e .
+```
+
 I also add this to my .bashrc for convenience:
 
 ```

inference_scans/run_law.sh

@@ -0,0 +1,172 @@
+#!/usr/bin/env bash
+# shellcheck disable=SC2154
+
+####################################################################################################
+# Script for running HH inference 'law' commands
+#
+# Usage:
+# 0) snapshot: --snapshot
+# 1) limits at point: --limpoint (--vbf)  # --vbf runs it for the VBF k2v=0 point
+# 2) 1D kL limit scan: --limkl
+# 3) 1D k2V limit scan: --limc2v
+# 4) impacts: --impacts  (should replot the merged json with plotImpacts.py)
+#
+# General options:
+#   --printdeps:    to print dependencies only
+#   --inject:       inject parameters from preliminary fit
+#   --rmoutput X:   remove outputs up to depth X
+#
+####################################################################################################
+
+####################################################################################################
+# Read options
+####################################################################################################
+
+rmoutput=0
+limits_at_point=0
+pois="r"
+vbfargs=""
+limits_1d_kl=0
+limits_1d_c2v=0
+impacts=0
+snapshot=0
+printdeps=""
+cards="$Cbbww4q"
+unblinded=$UNBLINDED
+
+options=$(getopt -o "ips" --long "limpoint,limkl,limc2v,impacts,printdeps,inject,vbf,unblinded,snapshot,rmoutput:" -- "$@")
+eval set -- "$options"
+
+while true; do
+    case "$1" in
+        -s|--snapshot)
+            snapshot=1
+            ;;
+        --limpoint)
+            limits_at_point=1
+            ;;
+        --limkl)
+            limits_1d_kl=1
+            ;;
+        --limc2v)
+            limits_1d_c2v=1
+            ;;
+        --impacts)
+            impacts=1
+            ;;
+        --vbf)
+            pois="r_qqhh"
+            vbfargs="--x-min 0.1 --x-max 10 --parameter-values C2V=0"
+            ;;
+        --unblinded)
+            unblinded=True
+            ;;
+        -i|--inject)
+            cards="$Cbbww4qInject"
+            ;;
+        -p|--printdeps)
+            printdeps="--print-deps -1"
+            ;;
+        --rmoutput)
+            shift
+            rmoutput=$1
+            ;;
+        --)
+            shift
+            break;;
+        \?)
+            echo "Invalid option: -$OPTARG" >&2
+            exit 1
+            ;;
+        :)
+            echo "Option -$OPTARG requires an argument." >&2
+            exit 1
+            ;;
+    esac
+    shift
+done
+
+# export DHI_CMS_POSTFIX="Supplementary"
+
+common_args="--file-types pdf,png --unblinded $unblinded --version $VERSION $printdeps --remove-output $rmoutput,a,y --campaign run2 --use-snapshot True"
+custom_args="--rMax 200 --setParameterRanges r_qqhh=-40,1000"
+
+
+if [ $snapshot = 1 ]; then
+    law run Snapshot \
+        --datacards $cards \
+        --custom-args=$custom_args
+fi
+
+
+if [ $limits_at_point = 1 ]; then
+    law run PlotUpperLimitsAtPoint \
+        $common_args $vbfargs \
+         --datacard-names bbVV \
+        --multi-datacards $cards \
+        --pois $pois \
+        --show-parameters kl,kt,C2V,CV \
+        --UpperLimits-workflow "htcondor" \
+        --UpperLimits-tasks-per-job 1 \
+        --x-log \
+        --h-lines 1 \
+        --save-hep-data True \
+        --Snapshot-custom-args="$custom_args"
+fi
+
+
+if [ $limits_1d_kl = 1 ]; then
+    law run PlotUpperLimits \
+        $common_args \
+        --version "$VERSION" \
+        --datacards $cards \
+        --xsec fb \
+        --pois r \
+        --scan-parameters kl,-29,-26,2:kl,-22,-18,2:kl,-9,-5,2:kl,-1,0,2:kl,2,5,4:kl,10,11,2:kl,11,30,2 \
+        --UpperLimits-workflow "htcondor" \
+        --UpperLimits-tasks-per-job 1 \
+        --y-log \
+        --show-parameters "kt,C2V,CV" \
+        --br bbww \
+        --save-ranges \
+        --save-hep-data False \
+        --Snapshot-custom-args="$custom_args" \
+        --frozen-groups signal_norm_xsbr
+fi
+
+if [ $limits_1d_c2v = 1 ]; then
+    law run PlotUpperLimits \
+        $common_args \
+        --version "$VERSION" \
+        --datacards $cards \
+        --xsec fb \
+        --pois r \
+        --scan-parameters C2V,-1,-0.2,3:C2V,-0.1,0.3,5:C2V,0.9,1.1,2:C2V,1.5,2.1,7:C2V,2.2,3,3 \
+        --UpperLimits-workflow "htcondor" \
+        --UpperLimits-tasks-per-job 1 \
+        --y-log \
+        --show-parameters "kt,kl,CV" \
+        --br bbww \
+        --save-ranges \
+        --save-hep-data False \
+        --Snapshot-custom-args="$custom_args" \
+        --frozen-groups signal_norm_xsbr
+fi
+
+
+if [ $impacts = 1 ]; then
+    law run PlotPullsAndImpacts \
+        $common_args \
+        --datacards $cards \
+        --pois r \
+        --PullsAndImpacts-workflow "htcondor" \
+        --PullsAndImpacts-tasks-per-job 10 \
+        --PullsAndImpacts-custom-args="--rMin -40 --rMax 200" \
+        --parameters-per-page 40 \
+        --order-by-impact \
+        --labels "nuisance_renames.py" \
+        --skip-parameters "*dataResidual_Bin*" \
+        --page -1 \
+        --pull-range 3 \
+        --Snapshot-custom-args="$custom_args"
+fi