Skip to content

Update script to produce HEPData for the publication #681

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

Merged
davidwalter2 merged 3 commits into from
Apr 7, 2026
Merged

Update script to produce HEPData for the publication #681

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
60a3f5d
update HEPData entry for final publication
cippy Apr 2, 2026
a0358c2
polish HEPData script
cippy Apr 2, 2026
69cfc73
Merge branch 'main' of github.com:WMass/WRemnants into hepdata_final
cippy Apr 2, 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
45 changes: 21 additions & 24 deletions scripts/hepdata/SMP-23-002/wmass_hepdata.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -141,7 +141,6 @@ def make_pulls_and_constraints(sub, logger=None, verbose=4):
if logger == None:
logger = logging.setup_logger(os.path.basename(__file__), verbose)

# TODO: use a common path
basePath = data_dir + "/hepdata/SMP_23_002/"

# tag name : {file name, matrix name inside file}
Expand Down Expand Up @@ -253,7 +252,6 @@ def make_cross_section(sub, logger=None, verbose=4):
if logger == None:
logger = logging.setup_logger(os.path.basename(__file__), verbose)

# TODO: use a common path
basePath = data_dir + "hepdata/SMP_23_002/"
input_dict = {
"W_pt": {
Expand All @@ -262,7 +260,7 @@ def make_cross_section(sub, logger=None, verbose=4):
"observable": "DSIG/DPT",
"isDividedByBinWidth": True,
"addImage": True,
"reference": "Fig A12",
"reference": "Fig A7",
},
##
"Zwlike_y": {
Expand All @@ -271,7 +269,7 @@ def make_cross_section(sub, logger=None, verbose=4):
"observable": "DSIG/DYRAP",
"isDividedByBinWidth": True,
"addImage": True,
"reference": "Fig xx in supplementary material, and Fig A9 in paper",
"reference": "Fig B3 in supplementary material",
},
##
"Zwlike_pt": {
Expand All @@ -280,7 +278,7 @@ def make_cross_section(sub, logger=None, verbose=4):
"observable": "DSIG/DPT",
"isDividedByBinWidth": True,
"addImage": True,
"reference": "Fig 2 and A9",
"reference": "Fig 2",
},
##
"W_pt_helFit": {
Expand All @@ -289,7 +287,7 @@ def make_cross_section(sub, logger=None, verbose=4):
"observable": "DSIG/DPT",
"isDividedByBinWidth": True,
"addImage": True,
"reference": "Fig A14a",
"reference": "Fig B6a in supplementary material",
},
##
"W_y_helFit": {
Expand All @@ -298,7 +296,7 @@ def make_cross_section(sub, logger=None, verbose=4):
"observable": "DSIG/DYRAP",
"isDividedByBinWidth": True,
"addImage": True,
"reference": "Fig A14b",
"reference": "Fig B6b in supplementary material",
},
}

Expand Down Expand Up @@ -388,37 +386,46 @@ def make_mass_summary(sub, logger=None, verbose=4):
if logger == None:
logger = logging.setup_logger(os.path.basename(__file__), verbose)

# TODO: use a common path
basePath = data_dir + "/hepdata/SMP_23_002/"
input_dict = {
"Zwlike_ptZ_modeling": {
"file": "Zwlike/Wlike_modeling_summary_PerBinStatCorrUnc.root",
"hist": "mass_summary",
"desc": "Comparison of the nominal $\\mathit{m}_{Z}$ measurement and its uncertainty (total or only from $\\mathit{p}_{T}^{Z}$ modeling), with the alternative measurements using different approaches to the $\\mathit{p}_{T}^{Z}$ modeling and its uncertainty.",
"indepVarTitle": "$\\mathit{p}_{T}^{Z}$ model",
"observable": "M",
"addImage": True,
"reference": "Fig B5a in supplementary material",
},
#
"W_ptW_modeling": {
"file": "W/Wmass_modeling_summary.root",
"file": "W/Wmass_modeling_summary_PerBinStatCorrUnc.root",
"hist": "mass_summary",
"desc": "Comparison of the nominal $\\mathit{m}_{W}$ measurement and its uncertainty (total or only from $\\mathit{p}_{T}^{W}$ modeling), with the alternative measurements using different approaches to the $\\mathit{p}_{T}^{W}$ modeling and its uncertainty.",
"indepVarTitle": "$\\mathit{p}_{T}^{W}$ model",
"observable": "M",
"addImage": True,
"reference": "Fig A12",
"reference": "Fig B5b in supplementary material",
},
#
"W_pdfUnscaled": {
"file": "W/Wmass_pdf_summary_unscaled.root",
"hist": "mass_summary",
"desc": "Comparison of the nominal $\\mathit{m}_{W}$ measurement and its uncertainty (total or only from the CT18Z PDF set), with the alternative measurements using different PDFs and their uncertainty before the scaling procedure described in the paper text. The measured mass values are those reported in Table A.7 of the paper.",
"desc": "Comparison of the nominal $\\mathit{m}_{W}$ measurement and its uncertainty (total or only from the CT18Z PDF set), with the alternative measurements using different PDFs and their uncertainty before the scaling procedure described in the paper text. The measured mass values are those reported in Table A.3 of the paper.",
"indepVarTitle": "PDF set",
"observable": "M",
"addImage": True,
"reference": "Fig A18a",
"reference": "Fig B7a in supplementary material",
},
#
"W_pdfScaled": {
"file": "W/Wmass_pdf_summary_scaled.root",
"hist": "mass_summary",
"desc": "Comparison of the nominal $\\mathit{m}_{W}$ measurement and its uncertainty (total or only from the CT18Z PDF set), with the alternative measurements using different PDFs and their uncertainty after the scaling procedure described in the paper text. The measured mass values are those reported in Table A.7 of the paper. The scale factors for the PDF uncertainties from each set are reported in Table A.3 of the paper.",
"desc": "Comparison of the nominal $\\mathit{m}_{W}$ measurement and its uncertainty (total or only from the CT18Z PDF set), with the alternative measurements using different PDFs and their uncertainty after the scaling procedure described in the paper text. The measured mass values and scale factors for the PDF uncertainties from each set are those reported in Table A.3 of the paper.",
"indepVarTitle": "PDF set",
"observable": "M",
"addImage": True,
"reference": "Fig A18b",
"reference": "Fig B7b in supplementary material",
},
}

Expand Down Expand Up @@ -537,16 +544,6 @@ def make_mass_summary(sub, logger=None, verbose=4):
if outdir and not os.path.exists(outdir):
os.makedirs(outdir)

logger.warning("\n\n")
logger.warning(
r"To save PNG images for some plots and convert them within hepdata_lib"
)
logger.warning(r"you need to work within the following singularity image")
logger.warning(
r"/cvmfs/unpacked.cern.ch/gitlab-registry.cern.ch/bendavid/cmswmassdocker/wmassdevrolling\:v38_patch0"
)
logger.warning("\n\n")

if not args.skipCovariance:
logger.warning("Preparing covariance matrices as root files")
make_covariance_matrix(logger=logger, verbose=args.verbose)
Expand Down