Re-implementation of CMS single top

[jacoterh]

This PR reimplements CMS single top. All datasets agree with the legacy implementation.

• CMS_SINGLETOP_7TEV: no hepdata entry available. I manually checked legacy against the experimental paper and found perfect agreement.
• CMS_SINGLETOP_8TEV: hepdata entry is available. Filter implemented and found perfect agreement. Validphys report this PR vs legacy.

obs = "TCHANNEL-XSEC"

new_implementation = f"CMS_SINGLETOP_8TEV_{obs}"
old_implementation = f"CMS_SINGLETOP_8TEV_{obs}"

inp1 = {
    "dataset_input": {"dataset": f"{new_implementation}"},
    "theoryid": 40_000_000,
    "use_cuts": "internal",
    "t0pdfset": "NNPDF40_nnlo_as_01180",
    "use_t0": True,
}
inp2 = {
    "dataset_input": {"dataset": f"{old_implementation}", "variant": "legacy"},
    "theoryid": 40_000_000,
    "use_cuts": "internal",
    "t0pdfset": "NNPDF40_nnlo_as_01180",
    "use_t0": True,
}

covmat1 = API.covmat_from_systematics(**inp1)
covmat2 = API.covmat_from_systematics(**inp2)

t0_covmat1 = API.t0_covmat_from_systematics(**inp1)
t0_covmat2 = API.t0_covmat_from_systematics(**inp2)

print(f"Comparison for {new_implementation}")
print(np.all(np.isclose(covmat1, covmat2)))
print(np.all(np.isclose(t0_covmat1, t0_covmat2)))

>>> Comparison for CMS_SINGLETOP_8TEV_TCHANNEL-XSEC
>>> True
>>> True

• CMS_SINGLETOP_13TEV: no hepdata entry available. I manually checked legacy against the experimental paper and found perfect agreement.

All metadata has been updated accordingly.

[jacoterh]

Following up on @scarlehoff suggestion in PR#2185 to compare the overall experimental chi2 with all datasets combined, I'm attaching here the reports:

• CMS_SINGLETOP_7TEV
• CMS_SINGLETOP_8TEV
• CMS_SINGLETOP_13TEV

Regarding the ordering, for each dataset, legacy always appears first and is followed by the new commondata implementation.

[scarlehoff]

Thanks. Please, bump the version number and change the version comment just like in the ATLAS one.

Although, I think you forgot to add some of the filter.py files?

(btw, like in the ATLAS pr, rebase on top of master please)

[scarlehoff]

Please, remove this runcard from here (maybe put in the PR for reference)

(it will be eventually removed alongside all the rawdata anyway)

[jacoterh]

Hi @scarlehoff , also this PR should be ready for merging once the tests pass. The reason some of the filter scripts aren't there is because there's no HEPdata entry for those datasets.

Final validphys report can be found at https://vp.nnpdf.science/2txronVJRuG9xBFIOjwPWQ==/