Skip to content
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.

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

Adding RG to report #110

Draft
wants to merge 19 commits into
base: main
Choose a base branch
from
Draft

Adding RG to report #110

wants to merge 19 commits into from

Conversation

jacoterh
Copy link
Collaborator

@jacoterh jacoterh commented Dec 16, 2024

This PR includes:

  • The RGE matrix is saved as .pkl file to the result directory
  • Possible to load a precomputed RGE matrix when redoing a fit
  • moves the RGE modules to /rge
  • Report loads RG rotated predictions
  • Rewrote the Fisher heatmap to add support for multiple Fishers in a single heat map
  • Fisher heatmap supports custom ordering of WC and dataset grouping

todo:

  • load the saved RGE matrix when redoing a fit
  • Remove repetitive code
  • Update website with new plot options
  • Fix codefactor issues

Comparison between computing the RG from scratch and reloading it:
coefficient_histo

To enable reloading precomputed RG matrices, one should add its path to the runcard after the key rg_matrix

rge:
  init_scale: 5000.0
  obs_scale: dynamic # float or "dynamic"
  smeft_accuracy: integrate # options: integrate, leadinglog
  yukawa: top # options: top, full or none
  adm_QCD: False # if true, the EW couplings are set to zero
  rg_matrix: <path/to/rge_matrix.pkl>

The fisher heatmap can be added to the report as follows

fisher:

  # heatmap plot
  plot:
    summary_only: True # if True display only the fisher information per dataset group. If False will show the fine grained dataset per dataset
    figsize: [21, 35] # figure size
    title: true # if True display the fit label as title
    column_names:
      - 4H: "$t\\bar{t}t\\bar{t} + t\\bar{t}b\\bar{b}$"
      - Higgs: "$\\rm Higgs$"
      - LEP: "$\\rm LEP$"
      - tt: "$t\\bar{t}$"
      - ttV: "$t\\bar{t}V$"
      - t: "$t$"
      - tV: "$tV$"
      - VV: "$VV$"
      - FCCee-91: "$\\rm FCC-ee\\:91\\:GeV$"
      - FCCee-161: "$\\rm FCC-ee\\:161\\:GeV$"
      - FCCee-240: "$\\rm FCC-ee\\:240\\:GeV$"
      - FCCee-365: "$\\rm FCC-ee\\:365\\:GeV$"
    together: ["fit_without_rg", "with_with_rg"] # list of result IDs to be plotted together

from mpl_toolkits.axes_grid1 import make_axes_locatable
from numpy.f2py.crackfortran import privatepattern
Copy link
Collaborator Author

Choose a reason for hiding this comment

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

I think GitHub copilot added this by mistake haha

@@ -113,10 +126,16 @@ def load_configuration(self):

def load_datasets(self):
"""Load all datasets."""

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

to be improved

Copy link
Collaborator

@LucaMantani LucaMantani left a comment

Choose a reason for hiding this comment

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

Just few comments

src/smefit/analyze/fisher.py Outdated Show resolved Hide resolved
src/smefit/rge/rge.py Show resolved Hide resolved
f"fisher_{fit.name}",
)
links_list.append((f"fisher_{fit.name}", f"Table {fit.label}"))

if plot is not None:
fit_plot = copy.deepcopy(plot)
fit_plot.pop("together")
Copy link
Collaborator

Choose a reason for hiding this comment

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

This line makes it crash in case together is not specified in the runcard. I think it would be better if it is made sure that if it's not there, it gets set to a default value?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants