Polish analysis of a multiclosure test

pyproject.toml

@@ -50,6 +50,7 @@ wiki-upload = "validphys.scripts.wiki_upload:main"
 vp-get = "validphys.scripts.vp_get:main"
 vp-comparefits = "validphys.scripts.vp_comparefits:main"
 vp-fitrename = "validphys.scripts.vp_fitrename:main"
+vp-multiclosure = "validphys.scripts.vp_multiclosure:main"
 vp-checktheory = "validphys.scripts.vp_checktheory:main"
 vp-pdfrename = "validphys.scripts.vp_pdfrename:main"
 vp-pdffromreplicas = "validphys.scripts.vp_pdffromreplicas:main"

validphys2/examples/pca_bias_variance_report.yaml

@@ -0,0 +1,65 @@
+meta:
+  title: PCA report for bias and variance of multiclosure fits.
+  author: Mark N. Costantini
+  keywords: [multiclosure, PCA, bias, variance]
+
+
+dataset_inputs:
+  - {dataset: ATLAS_SINGLETOP_7TEV_TCHANNEL-XSEC, variant: legacy}
+  - {dataset: HERA_NC_225GEV_EP-SIGMARED, variant: legacy}
+  - {dataset: HERA_NC_251GEV_EP-SIGMARED, variant: legacy}
+  - {dataset: HERA_NC_300GEV_EP-SIGMARED, variant: legacy}
+  - {dataset: HERA_NC_318GEV_EP-SIGMARED, variant: legacy}
+  - {dataset: HERA_CC_318GEV_EM-SIGMARED, variant: legacy}
+  - {dataset: HERA_CC_318GEV_EP-SIGMARED, variant: legacy}
+
+
+
+theoryid: 200
+use_cuts: internal
+use_t0: True
+t0pdfset: 210223-mw-000_fakepdf
+
+
+explained_variance_ratio: 0.99
+
+##### Bootstrap parameters
+n_fit_max: 8
+n_fit: 8
+
+n_rep_max: 100
+n_rep: 60
+
+n_boot_multiclosure: 100
+rng_seed_mct_boot: 1
+
+
+fits:
+  - 25_5_2023_19_47_5_dis_pt1_mnc_commit_4d5d473c_filterseed_415295
+  - 25_5_2023_19_49_38_dis_pt1_mnc_commit_4d5d473c_filterseed_120750
+  - 25_5_2023_19_50_27_dis_pt1_mnc_commit_4d5d473c_filterseed_754299
+  - 25_5_2023_19_51_19_dis_pt1_mnc_commit_4d5d473c_filterseed_361864
+  - 25_5_2023_19_52_54_dis_pt1_mnc_commit_4d5d473c_filterseed_516909
+  - 25_5_2023_19_52_7_dis_pt1_mnc_commit_4d5d473c_filterseed_539528
+  - 25_5_2023_19_53_41_dis_pt1_mnc_commit_4d5d473c_filterseed_307907
+  - 25_5_2023_19_54_28_dis_pt1_mnc_commit_4d5d473c_filterseed_279829
+
+
+template_text: |
+
+  PCA report for bias and variance of multiclosure fits.
+  ------------------------------------------------------
+
+  ## Table of bias and variance
+  {@table_bias_variance_datasets@}
+
+
+  ## Bootstrapped table for full dataset  
+  {@bootstrapped_table_bias_variance_data@}
+
+  ## L2 condition number
+  {@plot_l2_condition_number@}
+
+
+actions_:
+  - report(main=true)

validphys2/src/validphys/closuretest/__init__.py

@@ -11,3 +11,5 @@
 from validphys.closuretest.multiclosure_pdf_output import *
 from validphys.closuretest.multiclosure_preprocessing import *
 from validphys.closuretest.multiclosure_pseudodata import *
+from validphys.closuretest.inconsistent_closuretest.multiclosure_inconsistent import *
+from validphys.closuretest.inconsistent_closuretest.multiclosure_inconsistent_output import *