cmshistsum with channel masks

Author: Aliya Nigamova

.github/workflows/cvmfs-ci.yml

@@ -243,6 +243,13 @@ jobs:
             text2workspace.py data/ci/template-analysis_shapeNInterp.txt -o ws_template-analysis.root --mass 200 --for-fits --no-wrappers --use-histsum
             combine -M MultiDimFit ws_template-analysis.root --algo singles  --setParameterRanges r=-5,5 --X-rtd FAST_VERTICAL_MORPH
 
+      - uses: ./.github/actions/run-in-cvmfs
+        name: CMSHistSum with channel masks
+        with:
+          script: |
+            text2workspace.py data/ci/htt_multiple_regions.txt -o ws_template-analysis_masks.root --mass 125 --channel-masks --for-fits --no-wrappers --use-histsum
+            combine -M MultiDimFit ws_template-analysis_masks.root --algo singles  --setParameterRanges r=-5,5 --setParameters mask_htt_tt_2_8TeV=1 --X-rtd FAST_VERTICAL_MORPH 
+
       - uses: ./.github/actions/run-in-cvmfs
         name: RooMultiPdf
         if: ${{ !matrix.CMSSW_VERSION || startsWith(matrix.CMSSW_VERSION, 'CMSSW_14') }}

test/CMakeLists.txt

@@ -228,6 +228,16 @@ ADD_COMBINE_TEST(cmshistsum_shapeN
     "combine -M MultiDimFit template-analysis_shapeNInterp.root --algo singles  --setParameterRanges r=-5,5 --X-rtd FAST_VERTICAL_MORPH"
 )
 
+# Template analysis CMSHistFunc with channel masks
+ADD_COMBINE_TEST(cmshistsum_channel_masks
+  COPY_TO_BUILDDIR ${REPO}/data/ci/htt_multiple_regions.txt
+                   ${REPO}/data/ci/htt_multiple_regions.input.root
+  T2W_COMMAND
+    text2workspace.py htt_multiple_regions.txt  -o ws_template-analysis-channel_masks.root --mass 125 --channel-masks -for-fits --no-wrappers --use-histsum
+  COMBINE_COMMANDS
+    "combine -M MultiDimFit ws_template-analysis-channel_masks.root --algo singles  --setParameterRanges r=-5,5 --setParameters mask_htt_tt_2_8TeV=1 --X-rtd FAST_VERTICAL_MORPH"
+)
+
 # Template-analysis datacard -> text2workspace
 COMBINE_ADD_TEST(template_analysis-text2workspace
     COMMAND text2workspace.py template-analysis_shape_autoMCStats.txt -o template-analysis_shape_autoMCStats.root

test/references/cmshistsum_channel_masks.out

@@ -0,0 +1,13 @@
+ <<< Combine >>> 
+ <<< v10.3.3 >>>
+>>> Random number generator seed is 123456
+>>> Method used is MultiDimFit
+Set Range of Parameter r To : (-5,5)
+Set Default Value of Parameter mask_htt_tt_2_8TeV To : 1
+>>> 1 out of 3 channels masked
+
+Doing initial fit: 
+
+ --- MultiDimFit ---
+best fit parameter values and profile-likelihood uncertainties: 
+   r :    +0.943   -0.790/+0.917 (68%)