Skip to content

Add closure-truth vertical lines to posterior plots#211

Open
WopkeTelman wants to merge 4 commits into
mainfrom
feature/posterior-vertical-lines
Open

Add closure-truth vertical lines to posterior plots#211
WopkeTelman wants to merge 4 commits into
mainfrom
feature/posterior-vertical-lines

Conversation

@WopkeTelman
Copy link
Copy Markdown
Collaborator

  • Draw vertical lines inside plot_posteriors
  • Pass closure values from report.py to the plotting layer
  • Backward compatible with existing calls

Related to #209

@ElieHammou ElieHammou self-requested a review February 19, 2026 16:05
@ElieHammou
Copy link
Copy Markdown
Collaborator

Thanks for opening this PR @WopkeTelman ! I have a few questions.
To start off, does this new feature require to add extra arguments in the analysis runcards? If so could you please add a working example in smefit_release/runcard/analyze?

@ElieHammou ElieHammou linked an issue Feb 20, 2026 that may be closed by this pull request
@ElieHammou ElieHammou added the enhancement New feature or request label Feb 20, 2026
"""Compute closure-truth Wilson coefficients from a coefficient config."""
fit_truth = {coeff_name: 0.0 for coeff_name in coeff_names}
if not coeff_cfg:
return fit_truth
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I would add a clear warning to the user here that the default values 0 is being used

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Also what does coeff_cfg stand for? It is not immediately obvious to me

coeff_set = set(coefficients.name)
for coeff_name in coeff_names:
if coeff_name not in coeff_set:
continue
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Again, maybe need a relevant warning here

return fit_truth

def _load_closure_truth_points_from_runcards(
self, closure_truth_runcards, coeff_names
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Are you loading the truth from the report runcard or from the projection runcard? Or both?
In any case, the possible options need to be better documented, both in the source code and in the example runcards

@ElieHammou
Copy link
Copy Markdown
Collaborator

Also did you implement this feature for the 2D corner plots already? I did not see an example runcard for it on our private repo

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

enhancement New feature or request

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Improve fit plotting for BSM closure tests

2 participants