Skip to content

Implement total phase space for unfolding #661

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
davidwalter2 wants to merge 34 commits into
base: main
Choose a base branch
from
Open

Implement total phase space for unfolding #661

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
34 commits
Select commit Hold shift + click to select a range
36ed5f1
Support total phase space unfolding
davidwalter2 Feb 16, 2026
d8a657c
Always use preFSR for gen mass cut
davidwalter2 Feb 19, 2026
4638eab
Fix previously broken filterProcs argument; use 'is' instead of '==' …
davidwalter2 Feb 19, 2026
d5d9733
Basic support for gen mass in unfolding
davidwalter2 Feb 19, 2026
14c4dd3
Fix full phase space unfolding for Z
davidwalter2 Feb 20, 2026
7b95f3c
Fix filtering procs for datasets
davidwalter2 Feb 20, 2026
7774298
Adapt inclusive xsec plotting script
davidwalter2 Feb 20, 2026
6e9c712
Merge branch '260220_fixProcFilters' of github.com:davidwalter2/WRemn…
davidwalter2 Feb 21, 2026
0d6a0ee
Separate out ROOT dependent input tools to avoid loading root library…
davidwalter2 Feb 21, 2026
baf4c66
Small updates for total phase space unfolding
davidwalter2 Feb 21, 2026
ed0ff51
Merge branch '260223_refactoring' of github.com:davidwalter2/WRemnant…
davidwalter2 Feb 26, 2026
93eb557
Update README and add some docstrings at the top of files to give a s…
davidwalter2 Feb 27, 2026
d10ad8b
Merge branch 'main' of github.com:WMass/WRemnants into 260216_totalPh…
davidwalter2 Feb 27, 2026
e50ccfb
Merge branch 'main' of github.com:WMass/WRemnants into 260216_totalPh…
davidwalter2 Mar 16, 2026
126f111
Add option to add masked channel to extrapolate unfolding to full pha…
davidwalter2 Mar 16, 2026
92eef2c
Merge branch 'main' of github.com:WMass/WRemnants into 260216_totalPh…
davidwalter2 Mar 17, 2026
6e55f1a
Merge branch 'main' of github.com:WMass/WRemnants into 260216_totalPh…
davidwalter2 Mar 19, 2026
e198759
Disable MC inflation for unfolding
davidwalter2 Mar 28, 2026
a96fad4
Merge branch 'main' of github.com:WMass/WRemnants into 260216_totalPh…
davidwalter2 Mar 28, 2026
fb79ae8
Update submodules
davidwalter2 Mar 30, 2026
0145bda
Fix CI
davidwalter2 Mar 31, 2026
b6b494d
Rename 'theory_helpers' to 'helicity_smoothing_helpers' and add optio…
davidwalter2 Apr 6, 2026
d0eb124
Update plotting scripts for unfolding and move them in separate area
davidwalter2 Apr 6, 2026
d2e515a
Edits to 'Rename 'theory_helpers' to 'helicity_smoothing_helpers' and…
davidwalter2 Apr 6, 2026
96b8eee
Merge branch 'main' of github.com:WMass/WRemnants into 260216_totalPh…
davidwalter2 Apr 7, 2026
07293b1
Support MiNNLO QCD scale from event weight hist (default is from heli…
davidwalter2 Apr 7, 2026
c86de65
Fixes to previous commits
davidwalter2 Apr 8, 2026
93c8e05
Replace overflow with np.inf upper edge on ABCD mt/relIso axes
davidwalter2 Apr 11, 2026
f68088f
Fix five bugs identified in code review
davidwalter2 Apr 12, 2026
a4cb03c
Merge branch 'main' of github.com:WMass/WRemnants into 260216_totalPh…
davidwalter2 Apr 12, 2026
75a3d33
Add support for simultaneous (extended)ABCD fit
davidwalter2 Apr 16, 2026
8ef2d11
Add support for simultaneous (extended)ABCD fit
davidwalter2 Apr 16, 2026
587c2bd
Merge branch 'main' of github.com:WMass/WRemnants into 260216_totalPh…
davidwalter2 Apr 16, 2026
79d90c6
Merge branch '260216_totalPhaseSpace' of github.com:davidwalter2/WRem…
davidwalter2 Apr 16, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
27 changes: 13 additions & 14 deletions .github/workflows/main.yml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -341,7 +341,7 @@ jobs:

# - name: bsm plot parameter correlations
# run: >-
# scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_hists_cov.py '$WREMNANTS_OUTDIR/WMass_eta_pt_charge_bsm/fitresults.hdf5'
# scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_cov.py '$WREMNANTS_OUTDIR/WMass_eta_pt_charge_bsm/fitresults.hdf5'
# -o $WEB_DIR/$PLOT_DIR/BSM --params 'WtoNMuMass50' massShiftW100MeV
# --title CMS --subtitle Preliminary --titlePos 0 --config 'wremnants/utilities/styles/styles.py' --correlation --showNumbers

Expand Down Expand Up @@ -475,14 +475,14 @@ jobs:

- name: plot xsec covariance
run: >-
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_hists_cov.py
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_cov.py
$WREMNANTS_OUTDIR/WMass_eta_pt_charge_poiAsNoi/fitresults.hdf5
-o $WEB_DIR/$PLOT_DIR/unfolding_mw --title CMS --subtitle Preliminary --postfix poiAsNoi
-m Project ch0_masked ptGen -m Project ch0_masked absEtaGen --config wremnants/utilities/styles/styles.py

- name: plot xsec corrrelation
run: >-
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_hists_cov.py
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_cov.py
$WREMNANTS_OUTDIR/WMass_eta_pt_charge_poiAsNoi/fitresults.hdf5
-o $WEB_DIR/$PLOT_DIR/unfolding_mw --title CMS --subtitle Preliminary --postfix poiAsNoi
-m Project ch0_masked ptGen -m Project ch0_masked absEtaGen --config wremnants/utilities/styles/styles.py --correlation
Expand Down Expand Up @@ -997,8 +997,7 @@ jobs:
mkdir -p $ALPHAS_OUTDIR
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run_python.sh scripts/rabbit/setupRabbit.py \
-i $ALPHAS_HIST_FILE --fitvar "$ALPHAS_FITVAR" --lumiScale $LUMI_SCALE --realData \
-o $ALPHAS_OUTDIR --noi alphaS --npUnc LatticeEigvars --pdfUncFromCorr \
--axlim ptll 0j 44j
-o $ALPHAS_OUTDIR --noi alphaS --npUnc LatticeEigvars --axlim ptll 0j 44j

- name: alphaS Z reco fit (${{ matrix.mode }})
run: >-
Expand Down Expand Up @@ -1099,14 +1098,14 @@ jobs:

- name: plot xsec covariance
run: >-
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_hists_cov.py
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_cov.py
$WREMNANTS_OUTDIR/ZMassDilepton_ptll_yll_poiAsNoi/fitresults.hdf5
-o $WEB_DIR/$PLOT_DIR/unfolding_dilepton --title CMS --subtitle Preliminary --postfix poiAsNoi
-m Project ch0_masked ptVGen -m Project ch0_masked absYVGen --config wremnants/utilities/styles/styles.py

- name: plot xsec corrrelation
run: >-
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_hists_cov.py
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_cov.py
$WREMNANTS_OUTDIR/ZMassDilepton_ptll_yll_poiAsNoi/fitresults.hdf5
-o $WEB_DIR/$PLOT_DIR/unfolding_dilepton --title CMS --subtitle Preliminary --postfix poiAsNoi
-m Project ch0_masked ptVGen -m Project ch0_masked absYVGen --config wremnants/utilities/styles/styles.py --correlation
Expand Down Expand Up @@ -1187,7 +1186,7 @@ jobs:
- name: unfolding plot pulls
run: >-
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_pulls_and_impacts.py
$WREMNANTS_OUTDIR/ZMassDilepton_ptll_yll_unfolding/fitresults.hdf5 -m ungrouped --sortDescending -s constraint --debug --noImpacts
$WREMNANTS_OUTDIR/ZMassDilepton_ptll_yll_unfolding/fitresults.hdf5 -m ungrouped --sortDescending -s constraint --debug --impactType none
-o $WEB_DIR/$PLOT_DIR/unfolding_dilepton --postfix unfolding -n 50 --otherExtensions png pdf

- name: plot xsec
Expand All @@ -1206,14 +1205,14 @@ jobs:

- name: plot xsec covariance
run: >-
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_hists_cov.py
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_cov.py
$WREMNANTS_OUTDIR/ZMassDilepton_ptll_yll_unfolding/fitresults.hdf5
-o $WEB_DIR/$PLOT_DIR/unfolding_dilepton --postfix unfolding --title CMS --subtitle Preliminary
-m Project ch0_masked ptVGen -m Project ch0_masked absYVGen --config wremnants/utilities/styles/styles.py

- name: plot xsec corrrelation
run: >-
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_hists_cov.py
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_cov.py
$WREMNANTS_OUTDIR/ZMassDilepton_ptll_yll_unfolding/fitresults.hdf5
-o $WEB_DIR/$PLOT_DIR/unfolding_dilepton --postfix unfolding --title CMS --subtitle Preliminary
-m Project ch0_masked ptVGen -m Project ch0_masked absYVGen --config wremnants/utilities/styles/styles.py --correlation
Expand Down Expand Up @@ -1420,7 +1419,7 @@ jobs:
# --fitvar "ptll-yll-cosThetaStarll_quantile-phiStarll_quantile" "eta-pt-charge" \
# --lumiScale $LUMI_SCALE $LUMI_SCALE \
# -o $ALPHAS_OUTDIR --noi alphaS wmass --postfix reco \
# --npUnc LatticeEigvars --pdfUncFromCorr
# --npUnc LatticeEigvars

# - name: alphaS WZ reco fit
# run: >-
Expand Down Expand Up @@ -1469,7 +1468,7 @@ jobs:
--genAxes "ptVGen-absYVGen-helicitySig" \
--scaleNormXsecHistYields "0.05" --allowNegativeExpectation \
--realData --systematicType normal --postfix unfolding \
--npUnc LatticeEigvars --pdfUncFromCorr \
--npUnc LatticeEigvars \
--rebin 2 2 # this is for CI only! needed for max_files=100

- name: alphaS Z unfolding fit
Expand Down Expand Up @@ -1551,14 +1550,14 @@ jobs:

- name: alphaS z gen-fit xsec covariance
run: >-
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_hists_cov.py
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_cov.py
$ALPHAS_Z_GENFIT
-o $ALPHAS_PLOT_DIR --postfix z_genfit --title CMS --subtitle Preliminary
-m Project chSigmaUL ptVGen -m Project chSigmaUL absYVGen --config wremnants/utilities/styles/styles.py

- name: alphaS z gen-fit xsec corrrelation
run: >-
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_hists_cov.py
scripts/ci/run_with_singularity.sh scripts/ci/setup_and_run.sh rabbit_plot_cov.py
$ALPHAS_Z_GENFIT
-o $ALPHAS_PLOT_DIR --postfix z_genfit --title CMS --subtitle Preliminary
-m Project chSigmaUL ptVGen -m Project chSigmaUL absYVGen --config wremnants/utilities/styles/styles.py --correlation
Expand Down
2 changes: 1 addition & 1 deletion rabbit
Open in desktop
Submodule rabbit updated 38 files
+149 −35 .github/workflows/main.yml
+4 −3 README.md
+240 −0 bin/merge_inputs.py
+16 −25 bin/rabbit_debug_inputdata.py
+150 −20 bin/rabbit_fit.py
+3 −3 bin/rabbit_limit.py
+77 −110 bin/rabbit_plot_cov.py
+265 −301 bin/rabbit_plot_hists.py
+20 −133 bin/rabbit_plot_hists_uncertainties.py
+7 −147 bin/rabbit_plot_inputdata.py
+6 −64 bin/rabbit_plot_likelihood_scan.py
+32 −73 bin/rabbit_plot_likelihood_scan2D.py
+202 −278 bin/rabbit_plot_pulls_and_impacts.py
+11 −53 bin/rabbit_print_impacts.py
+6 −19 bin/rabbit_print_pulls_and_constraints.py
+9 −30 bin/rabbit_text2hdf5.py
+1,416 −173 notebooks/tutorial_1_getting_started.ipynb
+273 −217 notebooks/tutorial_2_advanced.ipynb
+904 −0 notebooks/tutorial_3_top_mass.ipynb
+5 −5 pyproject.toml
+492 −572 rabbit/fitter.py
+0 −0 rabbit/impacts/__init__.py
+500 −0 rabbit/impacts/global_impacts.py
+116 −0 rabbit/impacts/nonprofiled_impacts.py
+71 −0 rabbit/impacts/traditional_impacts.py
+3 −0 rabbit/io_tools.py
+8 −5 rabbit/mappings/mapping.py
+243 −11 rabbit/parsing.py
+101 −0 rabbit/poi_models/poi_model.py
+0 −0 rabbit/regularization/__init__.py
+13 −0 rabbit/regularization/helpers.py
+234 −0 rabbit/regularization/lcurve.py
+17 −0 rabbit/regularization/regularizer.py
+92 −0 rabbit/regularization/svd.py
+6 −2 rabbit/tensorwriter.py
+107 −44 rabbit/workspace.py
+103 −31 tests/plot_epoch_loss_time.py
+1 −1 tests/plot_limits.py
Loading
Loading