From 0de5f637647a15d7b12feb44bd60e01849c30ee2 Mon Sep 17 00:00:00 2001 From: jacoterh <54140851+jacoterh@users.noreply.github.com> Date: Tue, 29 Oct 2024 11:33:52 +0000 Subject: [PATCH 1/2] re-implementing ATLAS single top production in new commondata format --- .../data_TCHANNEL-XSEC.yaml | 2 + .../ATLAS_SINGLETOP_13TEV/filter.py | 273 +++++++++++ .../kinematics_TCHANNEL-XSEC.yaml | 14 +- .../ATLAS_SINGLETOP_13TEV/metadata.yaml | 49 +- .../HEPData-ins1486394-v1-Table_3.yaml | 52 +++ .../HEPData-ins1486394-v1-Table_4.yaml | 249 ++++++++++ .../uncertainties_TCHANNEL-XSEC.yaml | 107 +++++ .../ATLAS_SINGLETOP_7TEV/data_T-Y-NORM.yaml | 4 + .../data_TBAR-Y-NORM.yaml | 4 + .../data_TCHANNEL-XSEC.yaml | 2 + .../commondata/ATLAS_SINGLETOP_7TEV/filter.py | 321 +++++++++++++ .../kinematics_T-Y-NORM.yaml | 38 +- .../kinematics_TBAR-Y-NORM.yaml | 38 +- .../kinematics_TCHANNEL-XSEC.yaml | 14 +- .../ATLAS_SINGLETOP_7TEV/metadata.yaml | 139 +++--- .../HEPData-ins1303905-v2-Table_17.yaml | 58 +++ .../HEPData-ins1303905-v2-Table_18.yaml | 58 +++ .../HEPData-ins1303905-v2-Table_2.yaml | 427 ++++++++++++++++++ .../HEPData-ins1303905-v2-Table_26.yaml | 172 +++++++ .../HEPData-ins1303905-v2-Table_27.yaml | 159 +++++++ .../HEPData-ins1303905-v2-Table_3.yaml | 41 ++ .../HEPData-ins1303905-v2-Table_34.yaml | 56 +++ .../HEPData-ins1303905-v2-Table_35.yaml | 56 +++ .../uncertainties_T-Y-NORM.yaml | 72 +++ .../uncertainties_TBAR-Y-NORM.yaml | 65 +++ .../uncertainties_TCHANNEL-XSEC.yaml | 147 ++++++ .../ATLAS_SINGLETOP_8TEV/data_T-RAP-NORM.yaml | 4 + .../data_TBAR-RAP-NORM.yaml | 4 + .../data_TCHANNEL-XSEC.yaml | 2 + .../commondata/ATLAS_SINGLETOP_8TEV/filter.py | 327 ++++++++++++++ .../kinematics_T-RAP-NORM.yaml | 36 +- .../kinematics_TBAR-RAP-NORM.yaml | 36 +- .../kinematics_TCHANNEL-XSEC.yaml | 14 +- .../ATLAS_SINGLETOP_8TEV/metadata.yaml | 137 +++--- .../HEPData-ins1512776-v2-Table_28.yaml | 52 +++ .../HEPData-ins1512776-v2-Table_29.yaml | 52 +++ .../HEPData-ins1512776-v2-Table_52.yaml | 58 +++ .../HEPData-ins1512776-v2-Table_53.yaml | 58 +++ .../HEPData-ins1512776-v2-Table_77.yaml | 282 ++++++++++++ .../HEPData-ins1512776-v2-Table_79.yaml | 282 ++++++++++++ .../HEPData-ins1512776-v2-Table_8.yaml | 84 ++++ .../uncertainties_T-RAP-NORM.yaml | 135 ++++++ .../uncertainties_TBAR-RAP-NORM.yaml | 135 ++++++ .../uncertainties_TCHANNEL-XSEC.yaml | 12 + validphys2/src/validphys/commondataparser.py | 1 + validphys2/src/validphys/process_options.py | 24 +- 46 files changed, 4097 insertions(+), 255 deletions(-) create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/data_TCHANNEL-XSEC.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/filter.py create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/rawdata/HEPData-ins1486394-v1-Table_3.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/rawdata/HEPData-ins1486394-v1-Table_4.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/uncertainties_TCHANNEL-XSEC.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_T-Y-NORM.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_TBAR-Y-NORM.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_TCHANNEL-XSEC.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/filter.py create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_17.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_18.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_2.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_26.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_27.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_3.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_34.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_35.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/uncertainties_T-Y-NORM.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/uncertainties_TBAR-Y-NORM.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/uncertainties_TCHANNEL-XSEC.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_T-RAP-NORM.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_TBAR-RAP-NORM.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_TCHANNEL-XSEC.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/filter.py create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_28.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_29.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_52.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_53.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_77.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_79.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_8.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/uncertainties_T-RAP-NORM.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/uncertainties_TBAR-RAP-NORM.yaml create mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/uncertainties_TCHANNEL-XSEC.yaml diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/data_TCHANNEL-XSEC.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/data_TCHANNEL-XSEC.yaml new file mode 100644 index 0000000000..df9c7c392f --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/data_TCHANNEL-XSEC.yaml @@ -0,0 +1,2 @@ +data_central: +- 1.72 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/filter.py b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/filter.py new file mode 100644 index 0000000000..6e8f8ae2c0 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/filter.py @@ -0,0 +1,273 @@ +import pathlib + +import numpy as np +import pandas as pd +import yaml + +from nnpdf_data.filter_utils.utils import prettify_float + +yaml.add_representer(float, prettify_float) + +NB_POINTS = 4 +MT_VALUE = 172.5 +SQRT_S = 13_000.0 + +from nnpdf_data.filter_utils.utils import cormat_to_covmat, covmat_to_artunc +from nnpdf_data.filter_utils.utils import symmetrize_errors as se + + +def load_yaml(table_id: int, version: int = 1) -> dict: + """Load the HEP data table in yaml format. + + Parameters + ---------- + table_id: int + table ID number + + Returns + ------- + dict: + ditionary containing the table contents + + """ + filename = f"HEPData-ins1486394-v{version}-Table_{table_id}" + table = pathlib.Path(f"./rawdata/{filename}.yaml") + + return yaml.safe_load(table.read_text()) + + +def get_kinematics(hepdata: dict, bin_index: list = [], indx: int = 0) -> list: + """Read the version and list of tables from metadata. + + Parameters + ---------- + hepdata: dict + dictionary containing all data info + bin_index: list + list of Non-empty bin index + indx: int + Column index from which to read, default=0 + + Returns + ------- + kinematics: list + kinematic info + + """ + if len(hepdata["independent_variables"]) > 0: + bins = hepdata["independent_variables"][indx]["values"] + if len(bin_index) > 0: + bins = [bins[i] for i in bin_index] + else: + bins = [] + + kinematics = [] + if len(bins) > 1: # differential case + for i in bin_index: + ymin, ymax = [float(value) for value in bins[i]["value"].split('-')] + kin_value = { + "y_t": {"min": ymin, "mid": (ymin + ymax) / 2, "max": ymax}, + "m_t2": {"min": None, "mid": MT_VALUE**2, "max": None}, + "sqrts": {"min": None, "mid": SQRT_S, "max": None}, + } + kinematics.append(kin_value) + else: # inclusive case + kin_value = { + "zero": {"min": 0, "mid": 0, "max": 0}, + "m_t2": {"min": None, "mid": MT_VALUE**2, "max": None}, + "sqrts": {"min": None, "mid": SQRT_S, "max": None}, + } + kinematics.append(kin_value) + + return kinematics + + +def get_data_values(hepdata: dict, bin_index: list, indx: int = 0) -> list: + """Extract the central values from the HepData yaml file. + + Parameters + ---------- + hepdata: dict + dictionary containing all data info + bin_index: list + Bin indices that must be parsed + indx: int + Column index from which to read the central value, default=0 + + Returns + ------- + list: + list of dictionaries whose contents are the central values + + """ + central = hepdata["dependent_variables"][indx]["values"] + return np.array([central[i]["value"] for i in bin_index]) + + +def get_errors(hepdata: dict, bin_index: list) -> dict: + """ + Extract the uncertainties from hepdata and computes the shift of the central value in case of + asymmetric uncertainties + + Parameters + ---------- + hepdata: dict + Hepdata yaml file loaded as dictionary + bin_index: list + Bin indices that must be parsed + + Returns + ------- + dict: + Dictionary containing the errors (as pandas DataFrame) and shifts of central values + """ + # parse the systematics + errors = [] + shifts_cv = [] # relevant for asymmetric uncertainties + for i in bin_index: + bin = hepdata["dependent_variables"][i] + + # loop over sources of uncertainty + error_sources = [] + shift_cv = 0 + for source in bin["values"]: + if "errors" not in source: + error_sources.append(None) + elif "symerror" in source["errors"][0]: + error_sources.append(source["errors"][0]["symerror"]) + elif "asymerror" in source["errors"][0]: + + delta_min = source["errors"][0]["asymerror"]["minus"] + delta_plus = source["errors"][0]["asymerror"]["plus"] + se_delta, se_sigma = se(delta_plus, delta_min) + + # shift central value + shift_cv += se_delta + + error_sources.append(se_sigma) + shifts_cv.append(shift_cv) + errors.append(error_sources) + errors = np.array(errors) + + # get the description of the uncertainty from hepdata + error_names = [ + error_name["value"] for error_name in hepdata["independent_variables"][0]["values"] + ] + errors = pd.DataFrame( + errors, columns=error_names, index=[f"bin {i}" for i in range(errors.shape[0])] + ) + + return {"errors": errors, "shifts_cv": np.array(shifts_cv)} + + +def format_uncertainties(uncs: dict) -> list: + """Format the uncertainties to be dumped into the yaml file. + + Parameters + ---------- + uncs: dict + Dictionary containing the various source of uncertainties + + Returns + ------- + list: + list of dictionaries whose elements are the various errors + + """ + + combined_errors = [] + n_bins = uncs["systematics"].index.str.startswith("bin").sum() + for i in range(n_bins): + errors = {} + if "statistics" in uncs: + errors["stat"] = float(uncs["statistics"].iloc[i, 0]) + for j, unc in enumerate(uncs["systematics"].loc[f"bin {i}"].values): + errors[f"sys_corr_{j + 1}"] = float(unc) + + combined_errors.append(errors) + + return combined_errors + + +def dump_commondata(kinematics: list, data: list, errors: dict, obs: str) -> None: + """Function that generates and writes the commondata files. + + Parameters + ---------- + kinematics: list + list containing the kinematic values + data: list + list containing the central values + errors: dict + Dictionary containing the different errors + obs: str + Name to append to the file names + """ + + if "statistics" in errors: + error_definition = { + "stat": { + "description": "Uncorrelated statistical uncertainties", + "treatment": "ADD", + "type": "UNCORR", + } + } + else: + error_definition = {} + + n_sys = errors["systematics"].shape[1] + for i in range(n_sys): + + error_definition[f"sys_corr_{i + 1}"] = { + "description": errors["systematics"].columns[i], + "treatment": errors["systematics"].loc["treatment"].iloc[i], + "type": errors["systematics"].loc["type"].iloc[i], + } + + errors_formatted = format_uncertainties(errors) + with open(f"data_{obs}.yaml", "w") as file: + yaml.dump({"data_central": data.tolist()}, file, sort_keys=False) + + with open(f"kinematics_{obs}.yaml", "w") as file: + yaml.dump({"bins": kinematics}, file, sort_keys=False) + + with open(f"uncertainties_{obs}.yaml", "w") as file: + yaml.dump( + {"definitions": error_definition, "bins": errors_formatted}, file, sort_keys=False + ) + + +def main_filter() -> None: + """ + This filter produces the commondata for the following three observables: + 1) T-Y-NORM + 2) TBAR-Y-NORM + 3) TCHANNEL-XSEC (ratio tq/tqbar) + 1) and 2) follow the same procedure. + """ + + # TCHANNEL-XSEC RATIO + yaml_content_data = load_yaml(table_id=3, version=1) + yaml_uncertainties = load_yaml(table_id=4, version=1) + + data_central = get_data_values(yaml_content_data, bin_index=[0], indx=2) + kinematics = get_kinematics(yaml_content_data) + uncertainties = get_errors(yaml_uncertainties, bin_index=[2]) + + uncertainties_absolute = uncertainties["errors"] * data_central.reshape(1, -1) * 1e-2 + stat_abs = uncertainties_absolute.iloc[[0], [0]] + sys_abs = uncertainties_absolute.iloc[[0], 1:-2] + errors = {"statistics": stat_abs, "systematics": sys_abs} + + n_sys = sys_abs.shape[1] + sys_types = {"treatment": ["MULT"] * n_sys, "type": ["UNCORR"] * (n_sys - 1) + ["ATLASLUMI15"]} + sys_types_df = pd.DataFrame(sys_types, index=sys_abs.columns).T + errors["systematics"] = pd.concat([sys_types_df, sys_abs]) + + dump_commondata(kinematics, data_central, errors, obs="TCHANNEL-XSEC") + + return + + +if __name__ == "__main__": + main_filter() diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/kinematics_TCHANNEL-XSEC.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/kinematics_TCHANNEL-XSEC.yaml index 3980a84f4b..b6a60bc58d 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/kinematics_TCHANNEL-XSEC.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/kinematics_TCHANNEL-XSEC.yaml @@ -1,13 +1,13 @@ bins: -- k1: +- zero: + min: 0 + mid: 0 + max: 0 + m_t2: min: null - mid: 0.0 + mid: 29756.25 max: null - k2: - min: null - mid: 30032.89 - max: null - k3: + sqrts: min: null mid: 13000.0 max: null diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/metadata.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/metadata.yaml index 4a3cbf6e3a..74e215a83b 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/metadata.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/metadata.yaml @@ -1,6 +1,6 @@ setname: ATLAS_SINGLETOP_13TEV -version: 1 -version_comment: Port of old commondata +version: 2 +version_comment: Implementation of new data nnpdf_metadata: nnpdf31_process: SINGLETOP experiment: ATLAS @@ -8,10 +8,10 @@ arXiv: url: https://arxiv.org/abs/1609.03920 journal: JHEP 04 (2017) 086 iNSPIRE: - url: '' + url: https://inspirehep.net/literature/1486394 hepdata: - url: '' - version: -1 + url: https://www.hepdata.net/record/ins1486394 + version: 1 implemented_observables: - observable_name: TCHANNEL-XSEC observable: @@ -19,31 +19,31 @@ implemented_observables: label: ATLAS single top $R_{t}$ 13 TeV units: '' process_type: INC - tables: [] + tables: [3, 4] npoints: [] ndata: 1 plotting: - kinematics_override: inc_sqrt_scale + kinematics_override: identity dataset_label: ATLAS single top $R_{t}$ 13 TeV plot_x: idat kinematic_coverage: - - k1 - - k2 - - k3 + - zero + - m_t2 + - sqrts kinematics: variables: - k1: - description: Variable k1 - label: k1 - units: '' - k2: - description: Variable k2 - label: k2 - units: '' - k3: - description: Variable k3 - label: k3 - units: '' + zero: + description: "" + label: $$ + units: "" + m_t2: + description: "top mass squared" + label: $M^2$ + units: $GeV^2$ + sqrts: + description: "Center of mass energy" + label: $\sqrt{s}$ + units: $GeV$ file: kinematics_TCHANNEL-XSEC.yaml theory: conversion_factor: 1.0 @@ -51,10 +51,11 @@ implemented_observables: FK_tables: - - ATLAS_SINGLETOP_TCH_R_13TEV-ATLAS_SINGLETOP_TCH_R_13TEV_T - - ATLAS_SINGLETOP_TCH_R_13TEV-ATLAS_SINGLETOP_TCH_R_13TEV_TB - data_uncertainties: [] + data_uncertainties: [uncertainties_TCHANNEL-XSEC.yaml] variants: legacy: data_uncertainties: - uncertainties_legacy_TCHANNEL-XSEC.yaml - data_central: data_legacy_TCHANNEL-XSEC.yaml + data_central: data_legacy_TCHANNEL-XSEC.yaml + data_central: data_TCHANNEL-XSEC.yaml ported_from: ATLAS_SINGLETOP_TCH_R_13TEV diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/rawdata/HEPData-ins1486394-v1-Table_3.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/rawdata/HEPData-ins1486394-v1-Table_3.yaml new file mode 100644 index 0000000000..14ac8812ec --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/rawdata/HEPData-ins1486394-v1-Table_3.yaml @@ -0,0 +1,52 @@ +independent_variables: [] +dependent_variables: +- header: {name: '$\sigma(tq)$', units: pb} + qualifiers: + - {name: RE, value: P P --> TOP X} + - {name: SQRT(S), units: GeV, value: '13000.0'} + values: + - errors: + - {label: stat., symerror: 5} + - {label: syst., symerror: 27} + - {label: lumi., symerror: 3} + value: 156 +- header: {name: '$\sigma(\bar{t}q)$', units: pb} + qualifiers: + - {name: RE, value: P P --> TOPBAR X} + - {name: SQRT(S), units: GeV, value: '13000.0'} + values: + - errors: + - {label: stat., symerror: 4} + - {label: syst., symerror: 18} + - {label: lumi., symerror: 2} + value: 91 +- header: {name: '$R_t$'} + qualifiers: + - {name: RE, value: P P --> TOP X} + - {name: SQRT(S), units: GeV, value: '13000.0'} + values: + - errors: + - {label: stat., symerror: 0.09} + - {label: syst., symerror: 0.18} + value: 1.72 +- header: {name: '$\sigma(tq+\bar{t}q)$', units: pb} + qualifiers: + - {name: RE, value: P P --> TOP X} + - {name: SQRT(S), units: GeV, value: '13000.0'} + values: + - errors: + - {label: stat., symerror: 6} + - {label: syst., symerror: 45} + - {label: lumi., symerror: 5} + value: 247 +- header: {name: '$f_\mathrm{LV}\cdot |V_{tb}|$'} + qualifiers: + - {name: RE, value: P P --> TOP X} + - {name: SQRT(S), units: GeV, value: '13000.0'} + values: + - errors: + - {label: stat., symerror: 0.01} + - {label: syst., symerror: 0.09} + - {label: theor., symerror: 0.02} + - {label: lumi., symerror: 0.01} + value: 1.07 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/rawdata/HEPData-ins1486394-v1-Table_4.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/rawdata/HEPData-ins1486394-v1-Table_4.yaml new file mode 100644 index 0000000000..496bc62850 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/rawdata/HEPData-ins1486394-v1-Table_4.yaml @@ -0,0 +1,249 @@ +independent_variables: +- header: {name: Source} + values: + - {value: Data statistics} + - {value: Monte Carlo statistics} + - {value: Muon uncertainties} + - {value: Electron uncertainties} + - {value: JES} + - {value: Jet energy resolution} + - {value: '$E_\mathrm{T}^{\mathrm{miss}}$ modelling'} + - {value: '$b$-tagging efficiency'} + - {value: $c$-tagging efficiency} + - {value: Light-jet tagging efficiency} + - {value: Pile-up reweighting} + - {value: $tq$ parton shower generator} + - {value: $tq$ NLO matching} + - {value: $tq$ radiation} + - {value: '$t\bar{t}$, $Wt$, $t\bar{b}+\bar{t}b$ parton shower generator'} + - {value: '$t\bar{t}$, $Wt$, $t\bar{b}+\bar{t}b$ NLO matching'} + - {value: '$t\bar{t}$, $Wt$, $t\bar{b}+\bar{t}b$ radiation'} + - {value: PDF} + - {value: Multijet normalisation} + - {value: Other background normalisation} + - {value: Luminosity} + - {value: Total systematic uncertainty} + - {value: Total uncertainty} +dependent_variables: +- header: {name: '$\displaystyle \frac{\Delta \sigma(tq)}{\sigma(tq)}$', units: $\%$} + qualifiers: + - {name: RE, value: P P --> TOP X} + - {name: SQRT(S), units: GeV, value: '13000.0'} + values: + - errors: + - {symerror: 2.9, label: stat.} + value: '' + - errors: + - {symerror: 2.8, label: syst.} + value: '' + - errors: + - {symerror: 0.8, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 3.4, label: syst.} + value: '' + - errors: + - {symerror: 3.9, label: syst.} + value: '' + - errors: + - {symerror: 0.9, label: syst.} + value: '' + - errors: + - {symerror: 7.0, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.5, label: syst.} + value: '' + - errors: + - {symerror: 13.0, label: syst.} + value: '' + - errors: + - {symerror: 2.1, label: syst.} + value: '' + - errors: + - {symerror: 3.7, label: syst.} + value: '' + - errors: + - {symerror: 3.2, label: syst.} + value: '' + - errors: + - {symerror: 4.4, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.6, label: syst.} + value: '' + - errors: + - {symerror: 0.3, label: syst.} + value: '' + - errors: + - {symerror: 0.4, label: syst.} + value: '' + - errors: + - {symerror: 2.1, label: syst.} + value: '' + - errors: + - {symerror: 17.5, label: syst.} + value: '' + - errors: + - {symerror: 17.8, label: syst.} + value: '' +- header: {name: '$\displaystyle \frac{\Delta \sigma(\bar{t}q)}{\sigma(\bar{t}q)}$', units: $\%$} + qualifiers: + - {name: RE, value: P P --> TOP X} + - {name: SQRT(S), units: GeV, value: '13000.0'} + values: + - errors: + - {symerror: 4.1, label: stat.} + value: '' + - errors: + - {symerror: 4.2, label: syst.} + value: '' + - errors: + - {symerror: 0.9, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 4.1, label: syst.} + value: '' + - errors: + - {symerror: 3.1, label: syst.} + value: '' + - errors: + - {symerror: 1.2, label: syst.} + value: '' + - errors: + - {symerror: 6.9, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 2.2, label: syst.} + value: '' + - errors: + - {symerror: 14.3, label: syst.} + value: '' + - errors: + - {symerror: 0.7, label: syst.} + value: '' + - errors: + - {symerror: 3.4, label: syst.} + value: '' + - errors: + - {symerror: 4.4, label: syst.} + value: '' + - errors: + - {symerror: 8.6, label: syst.} + value: '' + - errors: + - {symerror: 1.1, label: syst.} + value: '' + - errors: + - {symerror: 0.9, label: syst.} + value: '' + - errors: + - {symerror: 2.0, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 2.1, label: syst.} + value: '' + - errors: + - {symerror: 20.0, label: syst.} + value: '' + - errors: + - {symerror: 20.4, label: syst.} + value: '' +- header: {name: '$\displaystyle \frac{\Delta R_{t}}{R_{t}}$', units: $\%$} + qualifiers: + - {name: RE, value: P P --> TOP X} + - {name: SQRT(S), units: GeV, value: '13000.0'} + values: + - errors: + - {symerror: 5.0, label: stat.} + value: '' + - errors: + - {symerror: 5.1, label: syst.} + value: '' + - errors: + - {symerror: 1.0, label: syst.} + value: '' + - errors: + - {symerror: 0.7, label: syst.} + value: '' + - errors: + - {symerror: 1.2, label: syst.} + value: '' + - errors: + - {symerror: 1.1, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.6, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 3.8, label: syst.} + value: '' + - errors: + - {symerror: 1.9, label: syst.} + value: '' + - errors: + - {symerror: 2.8, label: syst.} + value: '' + - errors: + - {symerror: 3.7, label: syst.} + value: '' + - errors: + - {symerror: 1.2, label: syst.} + value: '' + - errors: + - {symerror: 4.6, label: syst.} + value: '' + - errors: + - {symerror: 0.7, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.8, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 10.2, label: syst.} + value: '' + - errors: + - {symerror: 11.4, label: syst.} + value: '' diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/uncertainties_TCHANNEL-XSEC.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/uncertainties_TCHANNEL-XSEC.yaml new file mode 100644 index 0000000000..4ed50ad05c --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/uncertainties_TCHANNEL-XSEC.yaml @@ -0,0 +1,107 @@ +definitions: + stat: + description: Uncorrelated statistical uncertainties + treatment: ADD + type: UNCORR + sys_corr_1: + description: Monte Carlo statistics + treatment: MULT + type: UNCORR + sys_corr_2: + description: Muon uncertainties + treatment: MULT + type: UNCORR + sys_corr_3: + description: Electron uncertainties + treatment: MULT + type: UNCORR + sys_corr_4: + description: JES + treatment: MULT + type: UNCORR + sys_corr_5: + description: Jet energy resolution + treatment: MULT + type: UNCORR + sys_corr_6: + description: $E_\mathrm{T}^{\mathrm{miss}}$ modelling + treatment: MULT + type: UNCORR + sys_corr_7: + description: $b$-tagging efficiency + treatment: MULT + type: UNCORR + sys_corr_8: + description: $c$-tagging efficiency + treatment: MULT + type: UNCORR + sys_corr_9: + description: Light-jet tagging efficiency + treatment: MULT + type: UNCORR + sys_corr_10: + description: Pile-up reweighting + treatment: MULT + type: UNCORR + sys_corr_11: + description: $tq$ parton shower generator + treatment: MULT + type: UNCORR + sys_corr_12: + description: $tq$ NLO matching + treatment: MULT + type: UNCORR + sys_corr_13: + description: $tq$ radiation + treatment: MULT + type: UNCORR + sys_corr_14: + description: $t\bar{t}$, $Wt$, $t\bar{b}+\bar{t}b$ parton shower generator + treatment: MULT + type: UNCORR + sys_corr_15: + description: $t\bar{t}$, $Wt$, $t\bar{b}+\bar{t}b$ NLO matching + treatment: MULT + type: UNCORR + sys_corr_16: + description: $t\bar{t}$, $Wt$, $t\bar{b}+\bar{t}b$ radiation + treatment: MULT + type: UNCORR + sys_corr_17: + description: PDF + treatment: MULT + type: UNCORR + sys_corr_18: + description: Multijet normalisation + treatment: MULT + type: UNCORR + sys_corr_19: + description: Other background normalisation + treatment: MULT + type: UNCORR + sys_corr_20: + description: Luminosity + treatment: MULT + type: ATLASLUMI15 +bins: +- stat: 0.086 + sys_corr_1: 8.77200000e-02 + sys_corr_2: 0.0172 + sys_corr_3: 0.01204 + sys_corr_4: 2.06400000e-02 + sys_corr_5: 1.89200000e-02 + sys_corr_6: 0.0086 + sys_corr_7: 0.0086 + sys_corr_8: 1.03200000e-02 + sys_corr_9: 0.0086 + sys_corr_10: 0.06536 + sys_corr_11: 0.03268 + sys_corr_12: 0.04816 + sys_corr_13: 0.06364 + sys_corr_14: 2.06400000e-02 + sys_corr_15: 0.07912 + sys_corr_16: 0.01204 + sys_corr_17: 0.0086 + sys_corr_18: 0.03096 + sys_corr_19: 0.0086 + sys_corr_20: 0.0086 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_T-Y-NORM.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_T-Y-NORM.yaml new file mode 100644 index 0000000000..142729748e --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_T-Y-NORM.yaml @@ -0,0 +1,4 @@ +data_central: +- 5.92655000e-01 +- 0.5757 +- 0.45977 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_TBAR-Y-NORM.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_TBAR-Y-NORM.yaml new file mode 100644 index 0000000000..ac04d744e7 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_TBAR-Y-NORM.yaml @@ -0,0 +1,4 @@ +data_central: +- 0.75 +- 6.56700000e-01 +- 5.54167500e-01 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_TCHANNEL-XSEC.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_TCHANNEL-XSEC.yaml new file mode 100644 index 0000000000..da87d41bb0 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_TCHANNEL-XSEC.yaml @@ -0,0 +1,2 @@ +data_central: +- 2.04 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/filter.py b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/filter.py new file mode 100644 index 0000000000..98fdb3e463 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/filter.py @@ -0,0 +1,321 @@ +import pathlib + +import numpy as np +import pandas as pd +import yaml + +from nnpdf_data.filter_utils.utils import prettify_float + +yaml.add_representer(float, prettify_float) + +NB_POINTS = 4 +MT_VALUE = 172.5 +SQRT_S = 7_000.0 + +from nnpdf_data.filter_utils.utils import cormat_to_covmat, covmat_to_artunc +from nnpdf_data.filter_utils.utils import symmetrize_errors as se + + +def load_yaml(table_id: int, version: int = 1) -> dict: + """Load the HEP data table in yaml format. + + Parameters + ---------- + table_id: int + table ID number + + Returns + ------- + dict: + ditionary containing the table contents + + """ + filename = f"HEPData-ins1303905-v{version}-Table_{table_id}" + table = pathlib.Path(f"./rawdata/{filename}.yaml") + + return yaml.safe_load(table.read_text()) + + +def get_kinematics(hepdata: dict, bin_index: list = [], indx: int = 0) -> list: + """Read the version and list of tables from metadata. + + Parameters + ---------- + hepdata: dict + dictionary containing all data info + bin_index: list + list of Non-empty bin index + indx: int + Column index from which to read, default=0 + + Returns + ------- + kinematics: list + kinematic info + + """ + bins = hepdata["independent_variables"][indx]["values"] + + kinematics = [] + if len(bins) > 1: # differential case + for i in bin_index: + ymin, ymax = [float(value) for value in bins[i]["value"].split('-')] + kin_value = { + "y_t": {"min": ymin, "mid": (ymin + ymax) / 2, "max": ymax}, + "m_t2": {"min": None, "mid": MT_VALUE**2, "max": None}, + "sqrts": {"min": None, "mid": SQRT_S, "max": None}, + } + kinematics.append(kin_value) + else: # inclusive case + kin_value = { + "zero": {"min": 0, "mid": 0, "max": 0}, + "m_t2": {"min": None, "mid": MT_VALUE**2, "max": None}, + "sqrts": {"min": None, "mid": SQRT_S, "max": None}, + } + kinematics.append(kin_value) + + return kinematics + + +def get_data_values(hepdata: dict, bin_index: list, indx: int = 0) -> list: + """Extract the central values from the HepData yaml file. + + Parameters + ---------- + hepdata: dict + dictionary containing all data info + bin_index: list + Bin indices that must be parsed + indx: int + Column index from which to read the central value, default=0 + + Returns + ------- + list: + list of dictionaries whose contents are the central values + + """ + central = hepdata["dependent_variables"][indx]["values"] + return np.array([central[i]["value"] for i in bin_index]) + + +def get_errors(hepdata: dict, data_central: list, bin_index: list) -> dict: + """ + Extract the uncertainties from hepdata and computes the shift of the central value in case of + asymmetric uncertainties + + Parameters + ---------- + hepdata: dict + Hepdata yaml file loaded as dictionary + bin_index: list + Bin indices that must be parsed + + Returns + ------- + dict: + Dictionary containing the errors (as pandas DataFrame) and shifts of central values + """ + # parse the systematics + errors = [] + shifts_cv = [] # relevant for asymmetric uncertainties + for i, j in enumerate(bin_index): + bin = hepdata["dependent_variables"][j] + + # loop over sources of uncertainty + error_sources = [] + shift_cv = 0 + for source in bin["values"]: + if "errors" not in source: + error_sources.append(None) + elif "symerror" in source["errors"][0]: + # convert to absolute uncertainties + abs_sym = source["errors"][0]["symerror"] * data_central[i] + error_sources.append(abs_sym) + elif "asymerror" in source["errors"][0]: + delta_min = source["errors"][0]["asymerror"]["minus"] + delta_plus = source["errors"][0]["asymerror"]["plus"] + + # convert to absolute uncertainties + delta_min_abs = delta_min * data_central[i] * 1e-2 + delta_plus_abs = delta_plus * data_central[i] * 1e-2 + + # symmeterise uncertainties and shift central value + se_delta, se_sigma = se(delta_plus_abs, delta_min_abs) + shift_cv += se_delta + + error_sources.append(se_sigma) + shifts_cv.append(shift_cv) + errors.append(error_sources) + errors = np.array(errors) + + # get the description of the uncertainty from hepdata + error_names = [ + error_name["value"] for error_name in hepdata["independent_variables"][0]["values"] + ] + errors = pd.DataFrame( + errors, columns=error_names, index=[f"bin {i}" for i in range(errors.shape[0])] + ) + + return {"errors": errors, "shifts_cv": np.array(shifts_cv)} + + +def format_uncertainties(uncs: dict) -> list: + """Format the uncertainties to be dumped into the yaml file. + + Parameters + ---------- + uncs: dict + Dictionary containing the various source of uncertainties + + Returns + ------- + list: + list of dictionaries whose elements are the various errors + + """ + + combined_errors = [] + n_bins = uncs["systematics"].index.str.startswith("bin").sum() + for i in range(n_bins): + errors = {} + if "statistics" in uncs: + errors["stat"] = float(uncs["statistics"].iloc[i, 0]) + for j, unc in enumerate(uncs["systematics"].loc[f"bin {i}"].values): + errors[f"sys_corr_{j + 1}"] = float(unc) + + combined_errors.append(errors) + + return combined_errors + + +def dump_commondata(kinematics: list, data: list, errors: dict, obs: str) -> None: + """Function that generates and writes the commondata files. + + Parameters + ---------- + kinematics: list + list containing the kinematic values + data: list + list containing the central values + errors: dict + Dictionary containing the different errors + obs: str + Name to append to the file names + """ + + if "statistics" in errors: + error_definition = { + "stat": { + "description": "Uncorrelated statistical uncertainties", + "treatment": "ADD", + "type": "UNCORR", + } + } + else: + error_definition = {} + + n_sys = errors["systematics"].shape[1] + for i in range(n_sys): + + error_definition[f"sys_corr_{i + 1}"] = { + "description": errors["systematics"].columns[i], + "treatment": errors["systematics"].loc["treatment"].iloc[i], + "type": errors["systematics"].loc["type"].iloc[i], + } + + errors_formatted = format_uncertainties(errors) + with open(f"data_{obs}.yaml", "w") as file: + yaml.dump({"data_central": data.tolist()}, file, sort_keys=False) + + with open(f"kinematics_{obs}.yaml", "w") as file: + yaml.dump({"bins": kinematics}, file, sort_keys=False) + + with open(f"uncertainties_{obs}.yaml", "w") as file: + yaml.dump( + {"definitions": error_definition, "bins": errors_formatted}, file, sort_keys=False + ) + + +def main_filter() -> None: + """ + This filter produces the commondata for the following three observables: + 1) T-Y-NORM + 2) TBAR-Y-NORM + 3) TCHANNEL-XSEC (ratio tq/tqbar) + 1) and 2) follow the same procedure. + """ + + observable_tables = {"T-Y-NORM": [17, 34, 26], "TBAR-Y-NORM": [18, 35, 27]} + + for obs, tables in observable_tables.items(): + yaml_content_data = load_yaml(table_id=tables[0], version=2) + yaml_stat_corr = load_yaml(table_id=tables[1], version=2) + yaml_sys_sources = load_yaml(table_id=tables[2], version=2) + + kinematics = get_kinematics(yaml_content_data, bin_index=range(NB_POINTS - 1), indx=0) + data_central = get_data_values(yaml_content_data, bin_index=range(NB_POINTS - 1), indx=0) + uncertainties = get_errors(yaml_sys_sources, data_central, bin_index=range(NB_POINTS - 1)) + + # TODO: do we multiply relative uncertainties by the shifted central value or the unshifted one? + data_central += uncertainties["shifts_cv"] + uncertainties_absolute = uncertainties["errors"].dropna( + axis=1 + ) # some systematics are empty on hepdata + + # statistical correlated uncertainties + stat_cor_dict = yaml_stat_corr["dependent_variables"][0]["values"] + stat_cor = np.array([corr["value"] for corr in stat_cor_dict]).reshape(NB_POINTS, NB_POINTS) + stat_cor = stat_cor[:-1, :-1].flatten() + stat_unc = uncertainties_absolute.iloc[:, 0] + stat_covmat = cormat_to_covmat(err_list=list(stat_unc), cormat_list=stat_cor) + + # convert stat covmat to artificial systematics + stat_art = np.array(covmat_to_artunc(NB_POINTS - 1, stat_covmat)) + stat_art_df = pd.DataFrame( + stat_art, columns=[f"Statistical uncertainty {i}" for i in range(len(stat_art))] + ) + + systematics_df = uncertainties_absolute.iloc[:, 1:-2] + n_sys = systematics_df.shape[1] + n_stat = stat_art_df.shape[1] + n_sys_combined = n_stat + n_sys + # combine stat and sys uncertainties + sys_all_df = pd.concat([stat_art_df, systematics_df.reset_index(drop=True)], axis=1) + sys_all_df.index = uncertainties_absolute.iloc[:, 1:-2].index + + sys_types = { + "treatment": ["ADD"] * len(stat_art) + ["MULT"] * n_sys, + "type": ["CORR"] * n_sys_combined, + } + sys_types_df = pd.DataFrame(sys_types, index=sys_all_df.columns).T + sys_all_df = pd.concat([sys_types_df, sys_all_df]) + errors = {"systematics": sys_all_df} + + # Generate all the necessary files + dump_commondata(kinematics, data_central, errors, obs=obs) + + # TCHANNEL-XSEC RATIO + yaml_content_data = load_yaml(table_id=3, version=2) + yaml_uncertainties = load_yaml(table_id=2, version=2) + + data_central = get_data_values(yaml_content_data, bin_index=[0], indx=2) + kinematics = get_kinematics(yaml_content_data) + uncertainties = get_errors(yaml_uncertainties, data_central, bin_index=[2]) + + stat_abs = uncertainties["errors"].iloc[[0], [0]] + sys_abs = uncertainties["errors"].iloc[[0], 1:-2] + errors = {"statistics": stat_abs, "systematics": sys_abs} + + n_sys = sys_abs.shape[1] + sys_types = {"treatment": ["MULT"] * n_sys, "type": ["UNCORR"] * (n_sys - 1) + ["ATLASLUMI11"]} + sys_types_df = pd.DataFrame(sys_types, index=sys_abs.columns).T + errors["systematics"] = pd.concat([sys_types_df, sys_abs]) + + dump_commondata(kinematics, data_central, errors, obs="TCHANNEL-XSEC") + + return + + +if __name__ == "__main__": + main_filter() diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/kinematics_T-Y-NORM.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/kinematics_T-Y-NORM.yaml index bc6ff90b9b..8c13565883 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/kinematics_T-Y-NORM.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/kinematics_T-Y-NORM.yaml @@ -1,37 +1,37 @@ bins: -- k1: - min: null +- y_t: + min: 0.0 mid: 0.1 - max: null - k2: + max: 0.2 + m_t2: min: null - mid: 30032.89 + mid: 29756.25 max: null - k3: + sqrts: min: null mid: 7000.0 max: null -- k1: - min: null +- y_t: + min: 0.2 mid: 0.4 - max: null - k2: + max: 0.6 + m_t2: min: null - mid: 30032.89 + mid: 29756.25 max: null - k3: + sqrts: min: null mid: 7000.0 max: null -- k1: - min: null - mid: 0.85 - max: null - k2: +- y_t: + min: 0.6 + mid: 8.50000000e-01 + max: 1.1 + m_t2: min: null - mid: 30032.89 + mid: 29756.25 max: null - k3: + sqrts: min: null mid: 7000.0 max: null diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/kinematics_TBAR-Y-NORM.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/kinematics_TBAR-Y-NORM.yaml index bc6ff90b9b..8c13565883 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/kinematics_TBAR-Y-NORM.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/kinematics_TBAR-Y-NORM.yaml @@ -1,37 +1,37 @@ bins: -- k1: - min: null +- y_t: + min: 0.0 mid: 0.1 - max: null - k2: + max: 0.2 + m_t2: min: null - mid: 30032.89 + mid: 29756.25 max: null - k3: + sqrts: min: null mid: 7000.0 max: null -- k1: - min: null +- y_t: + min: 0.2 mid: 0.4 - max: null - k2: + max: 0.6 + m_t2: min: null - mid: 30032.89 + mid: 29756.25 max: null - k3: + sqrts: min: null mid: 7000.0 max: null -- k1: - min: null - mid: 0.85 - max: null - k2: +- y_t: + min: 0.6 + mid: 8.50000000e-01 + max: 1.1 + m_t2: min: null - mid: 30032.89 + mid: 29756.25 max: null - k3: + sqrts: min: null mid: 7000.0 max: null diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/kinematics_TCHANNEL-XSEC.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/kinematics_TCHANNEL-XSEC.yaml index 49b5ce05d2..b8395388a3 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/kinematics_TCHANNEL-XSEC.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/kinematics_TCHANNEL-XSEC.yaml @@ -1,13 +1,13 @@ bins: -- k1: +- zero: + min: 0 + mid: 0 + max: 0 + m_t2: min: null - mid: 0.0 + mid: 29756.25 max: null - k2: - min: null - mid: 30032.89 - max: null - k3: + sqrts: min: null mid: 7000.0 max: null diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/metadata.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/metadata.yaml index 186295afbd..a94a3a7ca5 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/metadata.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/metadata.yaml @@ -1,6 +1,6 @@ setname: ATLAS_SINGLETOP_7TEV -version: 1 -version_comment: Port of old commondata +version: 2 +version_comment: Implementation of new data nnpdf_metadata: nnpdf31_process: SINGLETOP experiment: ATLAS @@ -8,45 +8,45 @@ arXiv: url: https://arxiv.org/abs/1406.7844 journal: Phys. Rev. D 90 (2014) 112006 iNSPIRE: - url: '' + url: https://inspirehep.net/literature/1303905 hepdata: - url: '' - version: -1 + url: https://www.hepdata.net/record/ins1303905 + version: 2 implemented_observables: - observable_name: TBAR-Y-NORM observable: description: Heavy Quarks Production Single Quark Rapidity Distribution label: ATLAS single antitop $y$ (normalised) units: '' - process_type: HQP_YQ - tables: [] + process_type: SINGLETOP + tables: [18, 35, 27] npoints: [] ndata: 3 plotting: - kinematics_override: hqp_yq_sqrt_scale + kinematics_override: identity dataset_label: ATLAS single antitop $y$ (normalised) y_label: $(1/\sigma_{\bar t})d\sigma_{\bar t}/d|y_{\bar t}|$ figure_by: - - k3 - plot_x: k1 + - sqrts + plot_x: y_t kinematic_coverage: - - k1 - - k2 - - k3 + - y_t + - m_t2 + - sqrts kinematics: variables: - k1: - description: Variable k1 - label: k1 - units: '' - k2: - description: Variable k2 - label: k2 - units: '' - k3: - description: Variable k3 - label: k3 - units: '' + y_t: + description: "Antitop rapidity" + label: $y$ + units: "" + m_t2: + description: "top mass squared" + label: $M^2$ + units: $GeV^2$ + sqrts: + description: "Center of mass energy" + label: $\sqrt{s}$ + units: $GeV$ file: kinematics_TBAR-Y-NORM.yaml theory: normalization: @@ -60,47 +60,48 @@ implemented_observables: - ATLAS_SINGLETOP_TCH_DIFF_7TEV_TBAR_RAP_NORM-aMCfast_obs_1 - ATLAS_SINGLETOP_TCH_DIFF_7TEV_TBAR_RAP_NORM-aMCfast_obs_2 - - ATLAS_SINGLETOP_TCH_DIFF_7TEV_TBAR_RAP_NORM-ATLAS_SINGLETOP_TCH_R_7TEV_TB - data_uncertainties: [] + data_uncertainties: [uncertainties_TBAR-Y-NORM.yaml] variants: legacy: data_uncertainties: - uncertainties_legacy_TBAR-Y-NORM.yaml - data_central: data_legacy_TBAR-Y-NORM.yaml + data_central: data_legacy_TBAR-Y-NORM.yaml + data_central: data_TBAR-Y-NORM.yaml ported_from: ATLAS_SINGLETOP_TCH_DIFF_7TEV_TBAR_RAP_NORM - observable_name: T-Y-NORM observable: description: Heavy Quarks Production Single Quark Rapidity Distribution label: ATLAS single top $y_t$ (normalised) units: '' - process_type: HQP_YQ - tables: [] + process_type: SINGLETOP + tables: [17, 26, 34] npoints: [] ndata: 3 plotting: - kinematics_override: hqp_yq_sqrt_scale + kinematics_override: identity dataset_label: ATLAS single top $y_t$ (normalised) y_label: $(1/\sigma_t)d\sigma_t/d|y_t|$ figure_by: - - k3 - plot_x: k1 + - sqrts + plot_x: y_t kinematic_coverage: - - k1 - - k2 - - k3 + - y_t + - m_t2 + - sqrts kinematics: variables: - k1: - description: Variable k1 - label: k1 - units: '' - k2: - description: Variable k2 - label: k2 - units: '' - k3: - description: Variable k3 - label: k3 - units: '' + y_t: + description: "top rapidity" + label: $y$ + units: "" + m_t2: + description: "top mass squared" + label: $M^2$ + units: $GeV^2$ + sqrts: + description: "Center of mass energy" + label: $\sqrt{s}$ + units: $GeV$ file: kinematics_T-Y-NORM.yaml theory: normalization: @@ -114,12 +115,13 @@ implemented_observables: - ATLAS_SINGLETOP_TCH_DIFF_7TEV_T_RAP_NORM-aMCfast_obs_1 - ATLAS_SINGLETOP_TCH_DIFF_7TEV_T_RAP_NORM-aMCfast_obs_2 - - ATLAS_SINGLETOP_TCH_DIFF_7TEV_T_RAP_NORM-ATLAS_SINGLETOP_TCH_R_7TEV_T - data_uncertainties: [] + data_uncertainties: [uncertainties_T-Y-NORM.yaml] variants: legacy: data_uncertainties: - uncertainties_legacy_T-Y-NORM.yaml - data_central: data_legacy_T-Y-NORM.yaml + data_central: data_legacy_T-Y-NORM.yaml + data_central: data_T-Y-NORM.yaml ported_from: ATLAS_SINGLETOP_TCH_DIFF_7TEV_T_RAP_NORM - observable_name: TCHANNEL-XSEC observable: @@ -127,31 +129,31 @@ implemented_observables: label: ATLAS single top $R_{t}$ 7 TeV units: '' process_type: INC - tables: [] + tables: [2, 3] npoints: [] ndata: 1 plotting: - kinematics_override: inc_sqrt_scale + kinematics_override: identity dataset_label: ATLAS single top $R_{t}$ 7 TeV plot_x: idat kinematic_coverage: - - k1 - - k2 - - k3 + - zero + - m_t2 + - sqrts kinematics: variables: - k1: - description: Variable k1 - label: k1 - units: '' - k2: - description: Variable k2 - label: k2 - units: '' - k3: - description: Variable k3 - label: k3 - units: '' + zero: + description: "" + label: $$ + units: "" + m_t2: + description: "top mass squared" + label: $M^2$ + units: $GeV^2$ + sqrts: + description: "Center of mass energy" + label: $\sqrt{s}$ + units: $GeV$ file: kinematics_TCHANNEL-XSEC.yaml theory: conversion_factor: 1.0 @@ -159,10 +161,11 @@ implemented_observables: FK_tables: - - ATLAS_SINGLETOP_TCH_R_7TEV-ATLAS_SINGLETOP_TCH_R_7TEV_T - - ATLAS_SINGLETOP_TCH_R_7TEV-ATLAS_SINGLETOP_TCH_R_7TEV_TB - data_uncertainties: [] + data_uncertainties: [uncertainties_TCHANNEL-XSEC.yaml] variants: legacy: data_uncertainties: - uncertainties_legacy_TCHANNEL-XSEC.yaml - data_central: data_legacy_TCHANNEL-XSEC.yaml + data_central: data_legacy_TCHANNEL-XSEC.yaml + data_central: data_TCHANNEL-XSEC.yaml ported_from: ATLAS_SINGLETOP_TCH_R_7TEV diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_17.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_17.yaml new file mode 100644 index 0000000000..8b0158ab18 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_17.yaml @@ -0,0 +1,58 @@ +independent_variables: +- header: {name: '$|y(t)|$', units: ''} + values: + - {value: '0-0.2'} + - {value: '0.2-0.6'} + - {value: '0.6-1.1'} + - {value: '1.1-3.0'} +dependent_variables: +- header: {name: '$\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(t)|}$', units: ''} + values: + - value: 0.59 + errors: + - {symerror: 0.09} + - value: 0.57 + errors: + - {symerror: 0.05} + - value: 0.46 + errors: + - {symerror: 0.05} + - value: 0.223 + errors: + - {symerror: 0.019} +- header: {name: Total relative uncertainties, units: '%'} + values: + - value: '' + errors: + - {label: total, symerror: 15} + - value: '' + errors: + - {label: total, symerror: 9.0} + - value: '' + errors: + - {label: total, symerror: 9.7} + - value: '' + errors: + - {label: total, symerror: 8.5} +- header: {name: Relative uncertainties, units: '%'} + values: + - value: '' + errors: + - {label: total , symerror: ' '} + - {label: stat, symerror: 9.0} + - {label: syst, symerror: 11} + - value: '' + errors: + - {label: total , symerror: ' '} + - {label: stat, symerror: 6.4} + - {label: syst, symerror: 6.3} + - value: '' + errors: + - {label: total , symerror: ' '} + - {label: stat, symerror: 7.5} + - {label: syst, symerror: 6.2} + - value: '' + errors: + - {label: total , symerror: ' '} + - {label: stat, symerror: 4.9} + - {label: syst, symerror: 6.9} \ No newline at end of file diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_18.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_18.yaml new file mode 100644 index 0000000000..d1e0e057c5 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_18.yaml @@ -0,0 +1,58 @@ +independent_variables: +- header: {name: '$|y(\bar t)|$', units: ''} + values: + - {value: '0-0.2'} + - {value: '0.2-0.6'} + - {value: '0.6-1.1'} + - {value: '1.1-3.0'} +dependent_variables: +- header: {name: '$\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(\bar t)|}$', units: ''} + values: + - value: 0.75 + errors: + - {symerror: 0.14} + - value: 0.66 + errors: + - {symerror: 0.11} + - value: 0.555 + errors: + - {symerror: 0.095} + - value: 0.163 + errors: + - {symerror: 0.030} +- header: {name: Total relative uncertainties, units: '%'} + values: + - value: '' + errors: + - {label: total, symerror: 19} + - value: '' + errors: + - {label: total, symerror: 17} + - value: '' + errors: + - {label: total, symerror: 17} + - value: '' + errors: + - {label: total, symerror: 18} +- header: {name: Relative uncertainties, units: '%'} + values: + - value: '' + errors: + - {label: total , symerror: ' '} + - {label: stat, symerror: 13} + - {label: syst, symerror: 13} + - value: '' + errors: + - {label: total , symerror: ' '} + - {label: stat, symerror: 9.1} + - {label: syst, symerror: 14} + - value: '' + errors: + - {label: total , symerror: ' '} + - {label: stat, symerror: 11} + - {label: syst, symerror: 13} + - value: '' + errors: + - {label: total , symerror: ' '} + - {label: stat, symerror: 11} + - {label: syst, symerror: 15} diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_2.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_2.yaml new file mode 100644 index 0000000000..bb4460a575 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_2.yaml @@ -0,0 +1,427 @@ +independent_variables: +- header: {name: Source} + values: + - {value: Data statistical} + - {value: Monte Carlo statistical} + - {value: Multijet normalization} + - {value: Other background normalization} + - {value: JES detector} + - {value: JES statistical} + - {value: JES physics modeling} + - {value: 'JES $\eta$ intercalibration'} + - {value: JES mixed detector and modeling} + - {value: JES close-by jets} + - {value: JES pile-up} + - {value: JES flavor composition} + - {value: JES flavor response} + - {value: '$b$-JES'} + - {value: Jet energy resolution} + - {value: Jet vertex fraction} + - {value: '$b$-tagging efficiency'} + - {value: '$c$-tagging efficiency'} + - {value: Mistag efficiency} + - {value: '$b/\bar{b}$ acceptance'} + - {value: '$E_\mathrm{T}^{\mathrm{miss}}$ modeling'} + - {value: Lepton uncertainties} + - {value: PDF} + - {value: '$W+$jets shape variation'} + - {value: '$tq$ generator + parton shower'} + - {value: $tq$ scale variations} + - {value: '$t\bar{t}$ generator + parton shower'} + - {value: '$t\bar{t}$ ISR / FSR'} + - {value: Luminosity} + - {value: Total systematic} + - {value: Total} +dependent_variables: +- header: {name: '$\frac{\Delta \sigma(tq)}{\sigma(tq)}$', units: $\%$} + qualifiers: + - {name: RE, value: P P --> TOP X} + - {name: SQRT(S), units: GeV, value: '7000.0'} + values: + - errors: + - {symerror: 3.1, label: stat.} + value: '' + - errors: + - {symerror: 1.9, label: syst.} + value: '' + - errors: + - {symerror: 1.1, label: syst.} + value: '' + - errors: + - {symerror: 1.1, label: syst.} + value: '' + - errors: + - {symerror: 1.6, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 6.9, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.4, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 2.1, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 3.8, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.0, label: syst.} + value: '' + - errors: + - {symerror: 2.3, label: syst.} + value: '' + - errors: + - {symerror: 2.8, label: syst.} + value: '' + - errors: + - {symerror: 3.2, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.9, label: syst.} + value: '' + - errors: + - {symerror: 2.6, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.8, label: syst.} + value: '' + - errors: + - {symerror: 12.0, label: syst.} + value: '' + - errors: + - {symerror: 12.4, label: total} + value: '' +- header: {name: '$\frac{\Delta \sigma(\bar{t} q)}{\sigma(\bar{t}q)}$', units: $\%$} + qualifiers: + - {name: RE, value: P P --> TOPBAR X} + - {name: SQRT(S), units: GeV, value: '7000.0'} + values: + - errors: + - {symerror: 5.4, label: stat.} + value: '' + - errors: + - {symerror: 3.2, label: syst.} + value: '' + - errors: + - {symerror: 2.0, label: syst.} + value: '' + - errors: + - {symerror: 2.8, label: syst.} + value: '' + - errors: + - {symerror: 1.4, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 8.4, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.4, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.6, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 4.1, label: syst.} + value: '' + - errors: + - {symerror: 1.4, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 3.4, label: syst.} + value: '' + - errors: + - {symerror: 3.0, label: syst.} + value: '' + - errors: + - {symerror: 5.8, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.6, label: syst.} + value: '' + - errors: + - {symerror: 3.0, label: syst.} + value: '' + - errors: + - {symerror: 2.1, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.8, label: syst.} + value: '' + - errors: + - {symerror: 14.9, label: syst.} + value: '' + - errors: + - {symerror: 15.9, label: total} + value: '' +- header: {name: '$\frac{\Delta R_{t}}{R_{t}}$', units: $\%$} + qualifiers: + - {name: RE, value: P P --> TOP X} + - {name: SQRT(S), units: GeV, value: '7000.0'} + values: + - errors: + - {symerror: 6.2, label: stat.} + value: '' + - errors: + - {symerror: 3.6, label: syst.} + value: '' + - errors: + - {symerror: 1.6, label: syst.} + value: '' + - errors: + - {symerror: 1.9, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.8, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.2, label: syst.} + value: '' + - errors: + - {symerror: 1.0, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.0, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.6, label: syst.} + value: '' + - errors: + - {symerror: 1.0, label: syst.} + value: '' + - errors: + - {symerror: 2.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.6, label: syst.} + value: '' + - errors: + - {symerror: 1.0, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 6.1, label: syst.} + value: '' + - errors: + - {symerror: 8.7, label: total} + value: '' +- header: {name: '$\frac{\Delta \sigma(tq+\bar{t} q)}{\sigma(tq+\bar{t}q)}$', units: $\%$} + qualifiers: + - {name: RE, value: P P --> TOP X} + - {name: SQRT(S), units: GeV, value: '7000.0'} + values: + - errors: + - {symerror: 2.7, label: stat.} + value: '' + - errors: + - {symerror: 1.9, label: syst.} + value: '' + - errors: + - {symerror: 1.4, label: syst.} + value: '' + - errors: + - {symerror: 1.6, label: syst.} + value: '' + - errors: + - {symerror: 1.4, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 7.3, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.6, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.9, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 3.9, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 2.6, label: syst.} + value: '' + - errors: + - {symerror: 2.8, label: syst.} + value: '' + - errors: + - {symerror: 3.2, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.9, label: syst.} + value: '' + - errors: + - {symerror: 2.6, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 0.5, label: syst.} + value: '' + - errors: + - {symerror: 1.8, label: syst.} + value: '' + - errors: + - {symerror: 12.1, label: syst.} + value: '' + - errors: + - {symerror: 12.4, label: total} + value: '' diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_26.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_26.yaml new file mode 100644 index 0000000000..ef8d6d48c0 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_26.yaml @@ -0,0 +1,172 @@ +dependent_variables: +- header: {name: '[0,0.2]'} + qualifiers: + - {name: '$\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(t)|}$', value: '$|y(t)|$ bins'} + values: + - value: '' + errors: + - {symerror: 9.0, label: 'stat. [%]'} #1 + - value: '' + errors: + - {symerror: 5.9, label: 'stat. [%]'} #2 + - {value: ''} #blank line + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #1 + - value: '' + errors: + - asymerror: {minus: -1.5, plus: +1.6} + label: syst. [%] #2 + - errors: + - {symerror: 1.2, label: 'syst. [%]'} + value: '' #3 + - {value: ''} #blank line + - errors: + - {symerror: 1.7, label: 'syst. [%]'} + value: '' #1 + - value: '' + errors: + - asymerror: {minus: +9.8, plus: -9.0} + label: syst. [%] + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #3 + - {value: ''} #blank line + - value: '' #1 + errors: + - {symerror: 11, label: 'syst. [%]'} + - value: '' #2 + errors: + - {symerror: 15, label: 'total [%]'} +- header: {name: '[0.2,0.6]'} + qualifiers: + - {name: '$\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(t)|}$', value: '$|y(t)|$ bins'} + values: + - value: '' + errors: + - {symerror: 6.4, label: 'stat. [%]'} #1 + - value: '' + errors: + - {symerror: 4.8, label: 'stat. [%]'} #2 + - {value: ''} #blank line + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #1 + - value: '' + errors: + - asymerror: {minus: +2.3, plus: -0.5} + label: syst. [%] #2 + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #3 + - {value: ''} #blank line + - errors: + - {symerror: 1.8, label: 'syst. [%]'} + value: '' #1 + - value: '' + errors: + - asymerror: {minus: +3.0, plus: -2.8} + label: syst. [%] #2 + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #3 + - {value: ''} #blank line + - value: '' #1 + errors: + - {symerror: 6.3, label: 'syst. [%]'} + - value: '' #2 + errors: + - {symerror: 9.0, label: 'total [%]'} +- header: {name: '[0.6,1.1]'} + qualifiers: + - {name: '$\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(t)|}$', value: '$|y(t)|$ bins'} + values: + - value: '' + errors: + - {symerror: 7.5, label: 'stat. [%]'} #1 + - value: '' + errors: + - {symerror: 4.9, label: 'stat. [%]'} #2 + - {value: ''} #blank line + - errors: + - {symerror: 1.1, label: 'syst. [%]'} + value: '' #1 + - value: '' + errors: + - asymerror: {minus: -1.5, plus: +1.4} + label: syst. [%] #2 + - value: '' + errors: + - asymerror: {minus: +1.6, plus: -1.6} + label: syst. [%] #3 + - {value: ''} #blank line + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #1 + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #2 + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #3 + - {value: ''} #blank line + - value: '' #1 + errors: + - {symerror: 6.2, label: 'syst. [%]'} + - value: '' #2 + errors: + - {symerror: 9.7, label: 'total [%]'} +- header: {name: '[1.1,3.0]'} + qualifiers: + - {name: '$\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(t)|}$', value: '$|y(t)|$ bins'} + values: + - value: '' + errors: + - {symerror: 5.0, label: 'stat. [%]'} #1 + - value: '' + errors: + - {symerror: 3.2, label: 'stat. [%]'} #2 + - {value: ''} #blank line + - errors: + - {symerror: 1.0, label: 'syst. [%]'} + value: '' #1 + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #2 + - errors: + - {symerror: 1.0, label: 'syst. [%]'} + value: '' #3 + - {value: ''} #blank line + - errors: + - {symerror: 2.3, label: 'syst. [%]'} + value: '' #1 + - value: '' + errors: + - asymerror: {minus: -5.2, plus: +4.8} + label: syst. [%] #2 + - errors: + - {symerror: 1.5, label: 'syst. [%]'} + value: '' #3 + - {value: ''} #blank line + - value: '' #1 + errors: + - {symerror: 6.9, label: 'syst. [%]'} + - value: '' #2 + errors: + - {symerror: 8.5, label: 'total [%]'} +independent_variables: +- header: {name: Source} + values: + - {value: 'Data statistical'} + - {value: 'Monte Carlo statistical'} + - {value: ''} + - {value: 'Background normalization'} + - {value: 'JES $\eta$ intercalibration'} + - {value: 'Jet energy resolution'} + - {value: ''} + - {value: 'PDF'} + - {value: '$tq$ generator + parton shower'} + - {value: '$tq$ scale variation'} + - {value: ''} + - {value: 'Total systematic'} + - {value: 'Total'} diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_27.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_27.yaml new file mode 100644 index 0000000000..52400f68b4 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_27.yaml @@ -0,0 +1,159 @@ +dependent_variables: +- header: {name: '[0,0.2]'} + qualifiers: + - {name: '$\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(\bar t)|}$', value: '$|y(\bar t)|$ bins'} + values: + - value: '' + errors: + - {symerror: 13, label: 'stat. [%]'} #1 + - value: '' + errors: + - {symerror: 12, label: 'stat. [%]'} #2 + - {value: ''} #blank line + - errors: + - {symerror: 3.4, label: 'syst. [%]'} + value: '' #1 + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #2 + - {value: ''} #blank line + - errors: + - {symerror: 1.6, label: 'syst. [%]'} + value: '' #1 + - value: '' + errors: + - asymerror: {minus: +1.4, plus: -1.4} + label: syst. [%] #2 + - errors: + - {symerror: 1.9, label: 'syst. [%]'} + value: '' #3 + - {value: ''} #blank line + - value: '' #1 + errors: + - {symerror: 13, label: 'syst. [%]'} + - value: '' #2 + errors: + - {symerror: 19, label: 'total [%]'} +- header: {name: '[0.2,0.6]'} + qualifiers: + - {name: '$\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(\bar t)|}$', value: '$|y(\bar t)|$ bins'} + values: + - value: '' + errors: + - {symerror: 9.1, label: 'stat. [%]'} #1 + - value: '' + errors: + - {symerror: 11, label: 'stat. [%]'} #2 + - {value: ''} #blank line + - errors: + - {symerror: 2.4, label: 'syst. [%]'} + value: '' #1 + - value: '' + errors: + - asymerror: {minus: -1.9, plus: +0.5} + label: syst. [%] #2 + - {value: ''} #blank line + - errors: + - {symerror: 1.0, label: 'syst. [%]'} + value: '' #1 + - value: '' + errors: + - asymerror: {minus: +8.2, plus: -7.8} + label: syst. [%] #2 + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #3 + - {value: ''} #blank line + - value: '' #1 + errors: + - {symerror: 14, label: 'syst. [%]'} + - value: '' #2 + errors: + - {symerror: 17, label: 'total [%]'} +- header: {name: '[0.6,1.1]'} + qualifiers: + - {name: '$\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(\bar t)|}$', value: '$|y(\bar t)|$ bins'} + values: + - value: '' + errors: + - {symerror: 11, label: 'stat. [%]'} #1 + - value: '' + errors: + - {symerror: 12, label: 'stat. [%]'} #2 + - {value: ''} #blank line + - errors: + - {symerror: 1.1, label: 'syst. [%]'} + value: '' #1 + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #2 + - {value: ''} #blank line + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #1 + - value: '' + errors: + - asymerror: {minus: -4.3, plus: +4.0} + label: syst. [%] #2 + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #3 + - {value: ''} #blank line + - value: '' #1 + errors: + - {symerror: 13, label: 'syst. [%]'} + - value: '' #2 + errors: + - {symerror: 17, label: 'total [%]'} +- header: {name: '[1.1,3.0]'} + qualifiers: + - {name: '$\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(\bar t)|}$', value: '$|y(\bar t)|$ bins'} + values: + - value: '' + errors: + - {symerror: 11, label: 'stat. [%]'} #1 + - value: '' + errors: + - {symerror: 14, label: 'stat. [%]'} #2 + - {value: ''} #blank line + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #1 + - value: '' + errors: + - asymerror: {minus: -0.8, plus: +1.5} + label: syst. [%] #2 + - {value: ''} #blank line + - errors: + - {symerror: 1.8, label: 'syst. [%]'} + value: '' #1 + - value: '' + errors: + - asymerror: {minus: -3.9, plus: +3.8} + label: syst. [%] #2 + - errors: + - {symerror: 0.5, label: 'syst. [%]'} + value: '' #3 + - {value: ''} #blank line + - value: '' #1 + errors: + - {symerror: 15, label: 'syst. [%]'} + - value: '' #2 + errors: + - {symerror: 18, label: 'total [%]'} + +independent_variables: +- header: {name: Source} + values: + - {value: 'Data statistical'} + - {value: 'Monte Carlo statistical'} + - {value: ''} + - {value: 'Background normalization'} + - {value: 'JES $\eta$ intercalibration'} + - {value: ''} + - {value: 'PDF'} + - {value: '$tq$ generator + parton shower'} + - {value: '$tq$ scale variation'} + - {value: ''} + - {value: 'Total systematic'} + - {value: 'Total'} diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_3.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_3.yaml new file mode 100644 index 0000000000..60387b4461 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_3.yaml @@ -0,0 +1,41 @@ +dependent_variables: +- header: {name: '$\sigma(tq)$', units: '$\mathrm{pb}$'} + qualifiers: + - {name: RE(TOP), value: P P --> TOP X} + - {name: RE(TOPBAR), value: P P --> TOPBAR X} + values: + - errors: + - {label: stat, symerror: 1} + - {label: sys, symerror: 6} + value: 46 +- header: {name: '$\sigma(\bar tq)$', units: '$\mathrm{pb}$'} + qualifiers: + - {name: RE(TOP), value: P P --> TOP X} + - {name: RE(TOPBAR), value: P P --> TOPBAR X} + values: + - errors: + - {label: stat, symerror: 1} + - {label: sys, symerror: 3} + value: 23 +- header: {name: '$R_t$'} + qualifiers: + - {name: RE(TOP), value: P P --> TOP X} + - {name: RE(TOPBAR), value: P P --> TOPBAR X} + values: + - errors: + - {label: stat, symerror: 0.13} + - {label: sys, symerror: 0.12} + value: 2.04 +- header: {name: '$\sigma_{t}(tq+\bar tq)$', units: '$\mathrm{pb}$'} + qualifiers: + - {name: RE(TOP), value: P P --> TOP X} + - {name: RE(TOPBAR), value: P P --> TOPBAR X} + values: + - errors: + - {label: stat, symerror: 2} + - {label: sys, symerror: 8} + value: 68 +independent_variables: +- header: {name: '$\sqrt{s}$', units: '$\mathrm{GeV}$'} + values: + - {value: 7000.0} diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_34.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_34.yaml new file mode 100644 index 0000000000..3b808e431b --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_34.yaml @@ -0,0 +1,56 @@ +independent_variables: +- header: {name: '$|y(t)|$', units: ''} + values: + - {high: 0.2, low: 0.0} + - {high: 0.6, low: 0.2} + - {high: 1.1, low: 0.6} + - {high: 3.0, low: 1.1} + - {high: 0.2, low: 0.0} + - {high: 0.6, low: 0.2} + - {high: 1.1, low: 0.6} + - {high: 3.0, low: 1.1} + - {high: 0.2, low: 0.0} + - {high: 0.6, low: 0.2} + - {high: 1.1, low: 0.6} + - {high: 3.0, low: 1.1} + - {high: 0.2, low: 0.0} + - {high: 0.6, low: 0.2} + - {high: 1.1, low: 0.6} + - {high: 3.0, low: 1.1} +- header: {name: '$|y(t)|$', units: ''} + values: + - {high: 0.2, low: 0.0} + - {high: 0.2, low: 0.0} + - {high: 0.2, low: 0.0} + - {high: 0.2, low: 0.0} + - {high: 0.6, low: 0.2} + - {high: 0.6, low: 0.2} + - {high: 0.6, low: 0.2} + - {high: 0.6, low: 0.2} + - {high: 1.1, low: 0.6} + - {high: 1.1, low: 0.6} + - {high: 1.1, low: 0.6} + - {high: 1.1, low: 0.6} + - {high: 3.0, low: 1.1} + - {high: 3.0, low: 1.1} + - {high: 3.0, low: 1.1} + - {high: 3.0, low: 1.1} +dependent_variables: +- header: {name: Correlation} + values: + - {value: 1.00} + - {value: -0.11} + - {value: -0.14} + - {value: -0.31} + - {value: -0.11} + - {value: 1.00} + - {value: -0.18} + - {value: -0.51} + - {value: -0.14} + - {value: -0.18} + - {value: 1.00} + - {value: -0.60} + - {value: -0.31} + - {value: -0.51} + - {value: -0.60} + - {value: 1.00} diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_35.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_35.yaml new file mode 100644 index 0000000000..60d59cab21 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/rawdata/HEPData-ins1303905-v2-Table_35.yaml @@ -0,0 +1,56 @@ +independent_variables: +- header: {name: '$|y(\bar t)|$', units: ''} + values: + - {high: 0.2, low: 0.0} + - {high: 0.6, low: 0.2} + - {high: 1.1, low: 0.6} + - {high: 3.0, low: 1.1} + - {high: 0.2, low: 0.0} + - {high: 0.6, low: 0.2} + - {high: 1.1, low: 0.6} + - {high: 3.0, low: 1.1} + - {high: 0.2, low: 0.0} + - {high: 0.6, low: 0.2} + - {high: 1.1, low: 0.6} + - {high: 3.0, low: 1.1} + - {high: 0.2, low: 0.0} + - {high: 0.6, low: 0.2} + - {high: 1.1, low: 0.6} + - {high: 3.0, low: 1.1} +- header: {name: '$|y(\bar t)|$', units: ''} + values: + - {high: 0.2, low: 0.0} + - {high: 0.2, low: 0.0} + - {high: 0.2, low: 0.0} + - {high: 0.2, low: 0.0} + - {high: 0.6, low: 0.2} + - {high: 0.6, low: 0.2} + - {high: 0.6, low: 0.2} + - {high: 0.6, low: 0.2} + - {high: 1.1, low: 0.6} + - {high: 1.1, low: 0.6} + - {high: 1.1, low: 0.6} + - {high: 1.1, low: 0.6} + - {high: 3.0, low: 1.1} + - {high: 3.0, low: 1.1} + - {high: 3.0, low: 1.1} + - {high: 3.0, low: 1.1} +dependent_variables: +- header: {name: Correlation} + values: + - {value: 1.00} + - {value: -0.08} + - {value: -0.14} + - {value: -0.28} + - {value: -0.08} + - {value: 1.00} + - {value: -0.03} + - {value: -0.58} + - {value: -0.14} + - {value: -0.03} + - {value: 1.00} + - {value: -0.68} + - {value: -0.28} + - {value: -0.58} + - {value: -0.68} + - {value: 1.00} diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/uncertainties_T-Y-NORM.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/uncertainties_T-Y-NORM.yaml new file mode 100644 index 0000000000..0c4eaf543a --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/uncertainties_T-Y-NORM.yaml @@ -0,0 +1,72 @@ +definitions: + sys_corr_1: + description: Statistical uncertainty 0 + treatment: ADD + type: CORR + sys_corr_2: + description: Statistical uncertainty 1 + treatment: ADD + type: CORR + sys_corr_3: + description: Statistical uncertainty 2 + treatment: ADD + type: CORR + sys_corr_4: + description: Monte Carlo statistical + treatment: MULT + type: CORR + sys_corr_5: + description: Background normalization + treatment: MULT + type: CORR + sys_corr_6: + description: JES $\eta$ intercalibration + treatment: MULT + type: CORR + sys_corr_7: + description: Jet energy resolution + treatment: MULT + type: CORR + sys_corr_8: + description: PDF + treatment: MULT + type: CORR + sys_corr_9: + description: $tq$ generator + parton shower + treatment: MULT + type: CORR + sys_corr_10: + description: $tq$ scale variation + treatment: MULT + type: CORR +bins: +- sys_corr_1: 5.28140740e+00 + sys_corr_2: -5.48184643e-01 + sys_corr_3: -4.82642618e-02 + sys_corr_4: 3.481 + sys_corr_5: 0.295 + sys_corr_6: 9.15451118e-03 + sys_corr_7: 0.708 + sys_corr_8: 1.003 + sys_corr_9: 5.55603348e-02 + sys_corr_10: 0.295 +- sys_corr_1: -6.21517163e-01 + sys_corr_2: -1.82842573e+00 + sys_corr_3: -3.09491192e+00 + sys_corr_4: 2.73600000e+00 + sys_corr_5: 0.285 + sys_corr_6: 1.07849061e-02 + sys_corr_7: 0.285 + sys_corr_8: 1.026 + sys_corr_9: 1.65496435e-02 + sys_corr_10: 0.285 +- sys_corr_1: -7.18475744e-01 + sys_corr_2: -2.44794411e+00 + sys_corr_3: 2.32246900e+00 + sys_corr_4: 2.25400000e+00 + sys_corr_5: 5.06000000e-01 + sys_corr_6: 6.67792632e-03 + sys_corr_7: 7.36000000e-03 + sys_corr_8: 0.23 + sys_corr_9: 0.23 + sys_corr_10: 0.23 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/uncertainties_TBAR-Y-NORM.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/uncertainties_TBAR-Y-NORM.yaml new file mode 100644 index 0000000000..bfd3018af7 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/uncertainties_TBAR-Y-NORM.yaml @@ -0,0 +1,65 @@ +definitions: + sys_corr_1: + description: Statistical uncertainty 0 + treatment: ADD + type: CORR + sys_corr_2: + description: Statistical uncertainty 1 + treatment: ADD + type: CORR + sys_corr_3: + description: Statistical uncertainty 2 + treatment: ADD + type: CORR + sys_corr_4: + description: Monte Carlo statistical + treatment: MULT + type: CORR + sys_corr_5: + description: Background normalization + treatment: MULT + type: CORR + sys_corr_6: + description: JES $\eta$ intercalibration + treatment: MULT + type: CORR + sys_corr_7: + description: PDF + treatment: MULT + type: CORR + sys_corr_8: + description: $tq$ generator + parton shower + treatment: MULT + type: CORR + sys_corr_9: + description: $tq$ scale variation + treatment: MULT + type: CORR +bins: +- sys_corr_1: 9.70671334e+00 + sys_corr_2: -8.57805914e-01 + sys_corr_3: 3.26167250e-01 + sys_corr_4: 9.0 + sys_corr_5: 2.55 + sys_corr_6: 0.375 + sys_corr_7: 1.20000000e+00 + sys_corr_8: 1.05000000e-02 + sys_corr_9: 1.42500000e+00 +- sys_corr_1: -7.27052324e-01 + sys_corr_2: -4.32582320e+00 + sys_corr_3: -4.10252173e+00 + sys_corr_4: 7.26000000e+00 + sys_corr_5: 1.584 + sys_corr_6: 1.02671905e-02 + sys_corr_7: 0.66 + sys_corr_8: 5.28329897e-02 + sys_corr_9: 0.33 +- sys_corr_1: -1.35053932e+00 + sys_corr_2: -3.83652385e+00 + sys_corr_3: 4.55281818e+00 + sys_corr_4: 6.66 + sys_corr_5: 6.10500000e-01 + sys_corr_6: 0.2775 + sys_corr_7: 0.2775 + sys_corr_8: 2.30625707e-02 + sys_corr_9: 0.2775 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/uncertainties_TCHANNEL-XSEC.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/uncertainties_TCHANNEL-XSEC.yaml new file mode 100644 index 0000000000..a042395bee --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/uncertainties_TCHANNEL-XSEC.yaml @@ -0,0 +1,147 @@ +definitions: + stat: + description: Uncorrelated statistical uncertainties + treatment: ADD + type: UNCORR + sys_corr_1: + description: Monte Carlo statistical + treatment: MULT + type: UNCORR + sys_corr_2: + description: Multijet normalization + treatment: MULT + type: UNCORR + sys_corr_3: + description: Other background normalization + treatment: MULT + type: UNCORR + sys_corr_4: + description: JES detector + treatment: MULT + type: UNCORR + sys_corr_5: + description: JES statistical + treatment: MULT + type: UNCORR + sys_corr_6: + description: JES physics modeling + treatment: MULT + type: UNCORR + sys_corr_7: + description: JES $\eta$ intercalibration + treatment: MULT + type: UNCORR + sys_corr_8: + description: JES mixed detector and modeling + treatment: MULT + type: UNCORR + sys_corr_9: + description: JES close-by jets + treatment: MULT + type: UNCORR + sys_corr_10: + description: JES pile-up + treatment: MULT + type: UNCORR + sys_corr_11: + description: JES flavor composition + treatment: MULT + type: UNCORR + sys_corr_12: + description: JES flavor response + treatment: MULT + type: UNCORR + sys_corr_13: + description: $b$-JES + treatment: MULT + type: UNCORR + sys_corr_14: + description: Jet energy resolution + treatment: MULT + type: UNCORR + sys_corr_15: + description: Jet vertex fraction + treatment: MULT + type: UNCORR + sys_corr_16: + description: $b$-tagging efficiency + treatment: MULT + type: UNCORR + sys_corr_17: + description: $c$-tagging efficiency + treatment: MULT + type: UNCORR + sys_corr_18: + description: Mistag efficiency + treatment: MULT + type: UNCORR + sys_corr_19: + description: $b/\bar{b}$ acceptance + treatment: MULT + type: UNCORR + sys_corr_20: + description: $E_\mathrm{T}^{\mathrm{miss}}$ modeling + treatment: MULT + type: UNCORR + sys_corr_21: + description: Lepton uncertainties + treatment: MULT + type: UNCORR + sys_corr_22: + description: PDF + treatment: MULT + type: UNCORR + sys_corr_23: + description: $W+$jets shape variation + treatment: MULT + type: UNCORR + sys_corr_24: + description: $tq$ generator + parton shower + treatment: MULT + type: UNCORR + sys_corr_25: + description: $tq$ scale variations + treatment: MULT + type: UNCORR + sys_corr_26: + description: $t\bar{t}$ generator + parton shower + treatment: MULT + type: UNCORR + sys_corr_27: + description: $t\bar{t}$ ISR / FSR + treatment: MULT + type: UNCORR + sys_corr_28: + description: Luminosity + treatment: MULT + type: ATLASLUMI11 +bins: +- stat: 1.26480000e+01 + sys_corr_1: 7.344 + sys_corr_2: 3.26400000e+00 + sys_corr_3: 3.876 + sys_corr_4: 1.02 + sys_corr_5: 1.02 + sys_corr_6: 1.02 + sys_corr_7: 3.672 + sys_corr_8: 1.02 + sys_corr_9: 1.02 + sys_corr_10: 1.02 + sys_corr_11: 2.448 + sys_corr_12: 2.04 + sys_corr_13: 1.02 + sys_corr_14: 2.04 + sys_corr_15: 1.02 + sys_corr_16: 1.02 + sys_corr_17: 1.02 + sys_corr_18: 1.02 + sys_corr_19: 1.02 + sys_corr_20: 3.26400000e+00 + sys_corr_21: 2.04 + sys_corr_22: 5.1 + sys_corr_23: 1.02 + sys_corr_24: 1.02 + sys_corr_25: 1.02 + sys_corr_26: 3.26400000e+00 + sys_corr_27: 2.04 + sys_corr_28: 1.02 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_T-RAP-NORM.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_T-RAP-NORM.yaml new file mode 100644 index 0000000000..f7fce50017 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_T-RAP-NORM.yaml @@ -0,0 +1,4 @@ +data_central: +- 6.52567800e-01 +- 6.05705300e-01 +- 4.89343200e-01 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_TBAR-RAP-NORM.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_TBAR-RAP-NORM.yaml new file mode 100644 index 0000000000..cadfeb8148 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_TBAR-RAP-NORM.yaml @@ -0,0 +1,4 @@ +data_central: +- 7.00112700e-01 +- 6.26031300e-01 +- 5.41941150e-01 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_TCHANNEL-XSEC.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_TCHANNEL-XSEC.yaml new file mode 100644 index 0000000000..df9c7c392f --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_TCHANNEL-XSEC.yaml @@ -0,0 +1,2 @@ +data_central: +- 1.72 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/filter.py b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/filter.py new file mode 100644 index 0000000000..52260c8120 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/filter.py @@ -0,0 +1,327 @@ +import pathlib + +import numpy as np +import pandas as pd +import yaml + +from nnpdf_data.filter_utils.utils import prettify_float + +yaml.add_representer(float, prettify_float) + +NB_POINTS = 4 +MT_VALUE = 172.5 +SQRT_S = 8_000.0 + +from nnpdf_data.filter_utils.utils import cormat_to_covmat, covmat_to_artunc +from nnpdf_data.filter_utils.utils import symmetrize_errors as se + + +def load_yaml(table_id: int, version: int = 1) -> dict: + """Load the HEP data table in yaml format. + + Parameters + ---------- + table_id: int + table ID number + + Returns + ------- + dict: + ditionary containing the table contents + + """ + filename = f"HEPData-ins1512776-v{version}-Table_{table_id}" + table = pathlib.Path(f"./rawdata/{filename}.yaml") + + return yaml.safe_load(table.read_text()) + + +def get_kinematics(hepdata: dict, bin_index: list = [], indx: int = 0) -> list: + """Read the version and list of tables from metadata. + + Parameters + ---------- + hepdata: dict + dictionary containing all data info + bin_index: list + list of Non-empty bin index + indx: int + Column index from which to read, default=0 + + Returns + ------- + kinematics: list + kinematic info + + """ + if len(hepdata["independent_variables"]) > 0: + bins = hepdata["independent_variables"][indx]["values"] + if len(bin_index) > 0: + bins = [bins[i] for i in bin_index] + else: + bins = [] + + kinematics = [] + if len(bins) > 1: # differential case + for i in bin_index: + ymin, ymax = [float(value) for value in bins[i]["value"].split('-')] + kin_value = { + "y_t": {"min": ymin, "mid": (ymin + ymax) / 2, "max": ymax}, + "m_t2": {"min": None, "mid": MT_VALUE**2, "max": None}, + "sqrts": {"min": None, "mid": SQRT_S, "max": None}, + } + kinematics.append(kin_value) + else: # inclusive case + kin_value = { + "zero": {"min": 0, "mid": 0, "max": 0}, + "m_t2": {"min": None, "mid": MT_VALUE**2, "max": None}, + "sqrts": {"min": None, "mid": SQRT_S, "max": None}, + } + kinematics.append(kin_value) + + return kinematics + + +def get_data_values(hepdata: dict, bin_index: list, indx: int = 0) -> list: + """Extract the central values from the HepData yaml file. + + Parameters + ---------- + hepdata: dict + dictionary containing all data info + bin_index: list + Bin indices that must be parsed + indx: int + Column index from which to read the central value, default=0 + + Returns + ------- + list: + list of dictionaries whose contents are the central values + + """ + central = hepdata["dependent_variables"][indx]["values"] + return np.array([central[i]["value"] for i in bin_index]) + + +def get_errors(hepdata: dict, bin_index: list) -> dict: + """ + Extract the uncertainties from hepdata and computes the shift of the central value in case of + asymmetric uncertainties + + Parameters + ---------- + hepdata: dict + Hepdata yaml file loaded as dictionary + bin_index: list + Bin indices that must be parsed + + Returns + ------- + dict: + Dictionary containing the errors (as pandas DataFrame) and shifts of central values + """ + # parse the systematics + errors = [] + shifts_cv = [] # relevant for asymmetric uncertainties + for i in bin_index: + bin = hepdata["dependent_variables"][i] + + # loop over sources of uncertainty + error_sources = [] + shift_cv = 0 + for source in bin["values"]: + if "errors" not in source: + error_sources.append(None) + elif "symerror" in source["errors"][0]: + error_sources.append(source["errors"][0]["symerror"]) + elif "asymerror" in source["errors"][0]: + + delta_min = source["errors"][0]["asymerror"]["minus"] + delta_plus = source["errors"][0]["asymerror"]["plus"] + se_delta, se_sigma = se(delta_plus, delta_min) + + # shift central value + shift_cv += se_delta + + error_sources.append(se_sigma) + shifts_cv.append(shift_cv) + errors.append(error_sources) + errors = np.array(errors) + + # get the description of the uncertainty from hepdata + error_names = [ + error_name["value"] for error_name in hepdata["independent_variables"][0]["values"] + ] + errors = pd.DataFrame( + errors, columns=error_names, index=[f"bin {i}" for i in range(errors.shape[0])] + ) + + return {"errors": errors, "shifts_cv": np.array(shifts_cv)} + + +def format_uncertainties(uncs: dict) -> list: + """Format the uncertainties to be dumped into the yaml file. + + Parameters + ---------- + uncs: dict + Dictionary containing the various source of uncertainties + + Returns + ------- + list: + list of dictionaries whose elements are the various errors + + """ + + combined_errors = [] + n_bins = uncs["systematics"].index.str.startswith("bin").sum() + for i in range(n_bins): + errors = {} + if "statistics" in uncs: + errors["stat"] = float(uncs["statistics"].iloc[i, 0]) + for j, unc in enumerate(uncs["systematics"].loc[f"bin {i}"].values): + errors[f"sys_corr_{j + 1}"] = float(unc) + + combined_errors.append(errors) + + return combined_errors + + +def dump_commondata(kinematics: list, data: list, errors: dict, obs: str) -> None: + """Function that generates and writes the commondata files. + + Parameters + ---------- + kinematics: list + list containing the kinematic values + data: list + list containing the central values + errors: dict + Dictionary containing the different errors + obs: str + Name to append to the file names + """ + + if "statistics" in errors: + error_definition = { + "stat": { + "description": "Uncorrelated statistical uncertainties", + "treatment": "ADD", + "type": "UNCORR", + } + } + else: + error_definition = {} + + n_sys = errors["systematics"].shape[1] + for i in range(n_sys): + + error_definition[f"sys_corr_{i + 1}"] = { + "description": errors["systematics"].columns[i], + "treatment": errors["systematics"].loc["treatment"].iloc[i], + "type": errors["systematics"].loc["type"].iloc[i], + } + + errors_formatted = format_uncertainties(errors) + with open(f"data_{obs}.yaml", "w") as file: + yaml.dump({"data_central": data.tolist()}, file, sort_keys=False) + + with open(f"kinematics_{obs}.yaml", "w") as file: + yaml.dump({"bins": kinematics}, file, sort_keys=False) + + with open(f"uncertainties_{obs}.yaml", "w") as file: + yaml.dump( + {"definitions": error_definition, "bins": errors_formatted}, file, sort_keys=False + ) + + +def main_filter() -> None: + """ + This filter produces the commondata for the following three observables: + 1) T-Y-NORM + 2) TBAR-Y-NORM + 3) TCHANNEL-XSEC (ratio tq/tqbar) + 1) and 2) follow the same procedure. + """ + + observable_tables = {"T-RAP-NORM": [28, 52, 77], "TBAR-RAP-NORM": [29, 53, 79]} + + for obs, tables in observable_tables.items(): + yaml_content_data = load_yaml(table_id=tables[0], version=2) + yaml_stat_corr = load_yaml(table_id=tables[1], version=2) + yaml_sys_sources = load_yaml(table_id=tables[2], version=2) + + kinematics = get_kinematics(yaml_content_data, bin_index=range(NB_POINTS - 1), indx=0) + + data_central = ( + get_data_values(yaml_content_data, bin_index=range(NB_POINTS - 1), indx=1) * 1e-3 + ) + uncertainties = get_errors(yaml_sys_sources, bin_index=range(NB_POINTS - 1)) + + # TODO: do we multiply relative uncertainties by the shifted central value or the unshifted one? + data_central += uncertainties["shifts_cv"] * data_central * 1e-2 + uncertainties_absolute = ( + uncertainties["errors"] * data_central.reshape(NB_POINTS - 1, -1) * 1e-2 + ) + uncertainties_absolute = uncertainties_absolute.dropna( + axis=1 + ) # some systematics are empty on hepdata + + # statistical correlated uncertainties + stat_cor_dict = yaml_stat_corr["dependent_variables"][0]["values"] + stat_cor = np.array([corr["value"] for corr in stat_cor_dict]).reshape(NB_POINTS, NB_POINTS) + stat_cor = stat_cor[:-1, :-1].flatten() + stat_unc = uncertainties_absolute.iloc[:, 0] + stat_covmat = cormat_to_covmat(err_list=list(stat_unc), cormat_list=stat_cor) + + # convert stat covmat to artificial systematics + stat_art = np.array(covmat_to_artunc(NB_POINTS - 1, stat_covmat)) + stat_art_df = pd.DataFrame( + stat_art, columns=[f"Statistical uncertainty {i}" for i in range(len(stat_art))] + ) + + systematics_df = uncertainties_absolute.iloc[:, 1:-2] + n_sys = systematics_df.shape[1] + n_stat = stat_art_df.shape[1] + n_sys_combined = n_stat + n_sys + # combine stat and sys uncertainties + sys_all_df = pd.concat([stat_art_df, systematics_df.reset_index(drop=True)], axis=1) + sys_all_df.index = uncertainties_absolute.iloc[:, 1:-2].index + + sys_types = { + "treatment": ["ADD"] * len(stat_art) + ["MULT"] * n_sys, + "type": ["CORR"] * n_sys_combined, + } + sys_types_df = pd.DataFrame(sys_types, index=sys_all_df.columns).T + sys_all_df = pd.concat([sys_types_df, sys_all_df]) + errors = {"systematics": sys_all_df} + + # Generate all the necessary files + dump_commondata(kinematics, data_central, errors, obs=obs) + + # TCHANNEL-XSEC RATIO + yaml_content_data = load_yaml(table_id=8, version=2) + + data_central = get_data_values(yaml_content_data, bin_index=[6], indx=0) + kinematics = get_kinematics(yaml_content_data, bin_index=[6]) + + # Eq 13 in 1702.02859 (no hepdata entry) + sys_abs = pd.DataFrame( + np.array([0.07]).reshape(-1, 1), columns=["Sys uncertainty idx: 1"], index=["bin 0"] + ) + stat_abs = pd.DataFrame(np.array([0.05]).reshape(-1, 1)) + errors = {"statistics": stat_abs, "systematics": sys_abs} + + sys_types = {"treatment": ["MULT"], "type": ["UNCORR"]} + sys_types_df = pd.DataFrame(sys_types, index=["Sys uncertainty idx: 1"]).T + errors["systematics"] = pd.concat([sys_types_df, sys_abs]) + + dump_commondata(kinematics, data_central, errors, obs="TCHANNEL-XSEC") + + return + + +if __name__ == "__main__": + main_filter() diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/kinematics_T-RAP-NORM.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/kinematics_T-RAP-NORM.yaml index 5aaae912e9..3f94b802bc 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/kinematics_T-RAP-NORM.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/kinematics_T-RAP-NORM.yaml @@ -1,37 +1,37 @@ bins: -- k1: - min: null +- y_t: + min: 0.0 mid: 0.15 - max: null - k2: + max: 0.3 + m_t2: min: null - mid: 30032.89 + mid: 29756.25 max: null - k3: + sqrts: min: null mid: 8000.0 max: null -- k1: - min: null +- y_t: + min: 0.3 mid: 0.5 - max: null - k2: + max: 0.7 + m_t2: min: null - mid: 30032.89 + mid: 29756.25 max: null - k3: + sqrts: min: null mid: 8000.0 max: null -- k1: - min: null +- y_t: + min: 0.7 mid: 1.0 - max: null - k2: + max: 1.3 + m_t2: min: null - mid: 30032.89 + mid: 29756.25 max: null - k3: + sqrts: min: null mid: 8000.0 max: null diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/kinematics_TBAR-RAP-NORM.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/kinematics_TBAR-RAP-NORM.yaml index 5aaae912e9..3f94b802bc 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/kinematics_TBAR-RAP-NORM.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/kinematics_TBAR-RAP-NORM.yaml @@ -1,37 +1,37 @@ bins: -- k1: - min: null +- y_t: + min: 0.0 mid: 0.15 - max: null - k2: + max: 0.3 + m_t2: min: null - mid: 30032.89 + mid: 29756.25 max: null - k3: + sqrts: min: null mid: 8000.0 max: null -- k1: - min: null +- y_t: + min: 0.3 mid: 0.5 - max: null - k2: + max: 0.7 + m_t2: min: null - mid: 30032.89 + mid: 29756.25 max: null - k3: + sqrts: min: null mid: 8000.0 max: null -- k1: - min: null +- y_t: + min: 0.7 mid: 1.0 - max: null - k2: + max: 1.3 + m_t2: min: null - mid: 30032.89 + mid: 29756.25 max: null - k3: + sqrts: min: null mid: 8000.0 max: null diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/kinematics_TCHANNEL-XSEC.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/kinematics_TCHANNEL-XSEC.yaml index 04c429a38a..ea0fa3227b 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/kinematics_TCHANNEL-XSEC.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/kinematics_TCHANNEL-XSEC.yaml @@ -1,13 +1,13 @@ bins: -- k1: +- zero: + min: 0 + mid: 0 + max: 0 + m_t2: min: null - mid: 0.0 + mid: 29756.25 max: null - k2: - min: null - mid: 30032.89 - max: null - k3: + sqrts: min: null mid: 8000.0 max: null diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/metadata.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/metadata.yaml index 255be33ed6..c0080f554c 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/metadata.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/metadata.yaml @@ -1,6 +1,6 @@ setname: ATLAS_SINGLETOP_8TEV -version: 1 -version_comment: Port of old commondata +version: 2 +version_comment: Implementation of new data nnpdf_metadata: nnpdf31_process: SINGLETOP experiment: ATLAS @@ -8,45 +8,45 @@ arXiv: url: https://arxiv.org/abs/1702.02859 journal: Eur. Phys. J. C 77 (2017) 531 iNSPIRE: - url: '' + url: https://inspirehep.net/literature/1512776 hepdata: - url: '' - version: -1 + url: https://www.hepdata.net/record/ins1512776 + version: 2 implemented_observables: - observable_name: TBAR-RAP-NORM observable: description: Heavy Quarks Production Single Quark Rapidity Distribution label: ATLAS single antitop $y$ (normalised) units: '' - process_type: HQP_YQ - tables: [] + process_type: SINGLETOP + tables: [29, 53, 79 ] npoints: [] ndata: 3 plotting: - kinematics_override: hqp_yq_sqrt_scale + kinematics_override: identity dataset_label: ATLAS single antitop $y$ (normalised) y_label: $(1/\sigma_{\bar t})d\sigma_{\bar t}/d|y_{\bar t}|$ figure_by: - - k3 - plot_x: k1 + - sqrts + plot_x: y_t kinematic_coverage: - - k1 - - k2 - - k3 + - y_t + - m_t2 + - sqrts kinematics: variables: - k1: - description: Variable k1 - label: k1 - units: '' - k2: - description: Variable k2 - label: k2 - units: '' - k3: - description: Variable k3 - label: k3 - units: '' + y_t: + description: "Antitop rapidity" + label: $y$ + units: "" + m_t2: + description: "top mass squared" + label: $M^2$ + units: $GeV^2$ + sqrts: + description: "Center of mass energy" + label: $\sqrt{s}$ + units: $GeV$ file: kinematics_TBAR-RAP-NORM.yaml theory: normalization: @@ -60,47 +60,48 @@ implemented_observables: - ATLAS_SINGLETOP_TCH_DIFF_8TEV_TBAR_RAP_NORM-aMCfast_obs_1 - ATLAS_SINGLETOP_TCH_DIFF_8TEV_TBAR_RAP_NORM-aMCfast_obs_2 - - ATLAS_SINGLETOP_TCH_DIFF_8TEV_TBAR_RAP_NORM-ATLAS_SINGLETOP_TCH_R_8TEV_TB - data_uncertainties: [] + data_uncertainties: [uncertainties_TBAR-RAP-NORM.yaml] variants: legacy: data_uncertainties: - uncertainties_legacy_TBAR-RAP-NORM.yaml - data_central: data_legacy_TBAR-RAP-NORM.yaml + data_central: data_legacy_TBAR-RAP-NORM.yaml + data_central: data_TBAR-RAP-NORM.yaml ported_from: ATLAS_SINGLETOP_TCH_DIFF_8TEV_TBAR_RAP_NORM - observable_name: T-RAP-NORM observable: description: Heavy Quarks Production Single Quark Rapidity Distribution label: ATLAS single top $y_t$ (normalised) units: '' - process_type: HQP_YQ - tables: [] + process_type: SINGLETOP + tables: [28, 52, 77] npoints: [] ndata: 3 plotting: - kinematics_override: hqp_yq_sqrt_scale + kinematics_override: identity dataset_label: ATLAS single top $y_t$ (normalised) y_label: $(1/\sigma_t)d\sigma_t/d|y_t|$ figure_by: - - k3 - plot_x: k1 + - sqrts + plot_x: y_t kinematic_coverage: - - k1 - - k2 - - k3 + - y_t + - m_t2 + - sqrts kinematics: variables: - k1: - description: Variable k1 - label: k1 - units: '' - k2: - description: Variable k2 - label: k2 - units: '' - k3: - description: Variable k3 - label: k3 - units: '' + y_t: + description: "top rapidity" + label: $y$ + units: "" + m_t2: + description: "top mass squared" + label: $M^2$ + units: $GeV^2$ + sqrts: + description: "Center of mass energy" + label: $\sqrt{s}$ + units: $GeV$ file: kinematics_T-RAP-NORM.yaml theory: normalization: @@ -114,12 +115,13 @@ implemented_observables: - ATLAS_SINGLETOP_TCH_DIFF_8TEV_T_RAP_NORM-aMCfast_obs_1 - ATLAS_SINGLETOP_TCH_DIFF_8TEV_T_RAP_NORM-aMCfast_obs_2 - - ATLAS_SINGLETOP_TCH_DIFF_8TEV_T_RAP_NORM-ATLAS_SINGLETOP_TCH_R_8TEV_T - data_uncertainties: [] + data_uncertainties: [uncertainties_T-RAP-NORM.yaml] variants: legacy: data_uncertainties: - uncertainties_legacy_T-RAP-NORM.yaml - data_central: data_legacy_T-RAP-NORM.yaml + data_central: data_legacy_T-RAP-NORM.yaml + data_central: data_T-RAP-NORM.yaml ported_from: ATLAS_SINGLETOP_TCH_DIFF_8TEV_T_RAP_NORM - observable_name: TCHANNEL-XSEC observable: @@ -131,27 +133,27 @@ implemented_observables: npoints: [] ndata: 1 plotting: - kinematics_override: inc_sqrt_scale + kinematics_override: identity dataset_label: ATLAS single top $R_{t}$ 8 TeV plot_x: idat kinematic_coverage: - - k1 - - k2 - - k3 + - zero + - m_t2 + - sqrts kinematics: variables: - k1: - description: Variable k1 - label: k1 - units: '' - k2: - description: Variable k2 - label: k2 - units: '' - k3: - description: Variable k3 - label: k3 - units: '' + zero: + description: "" + label: $$ + units: "" + m_t2: + description: "top mass squared" + label: $M^2$ + units: $GeV^2$ + sqrts: + description: "Center of mass energy" + label: $\sqrt{s}$ + units: $GeV$ file: kinematics_TCHANNEL-XSEC.yaml theory: conversion_factor: 1.0 @@ -159,10 +161,11 @@ implemented_observables: FK_tables: - - ATLAS_SINGLETOP_TCH_R_8TEV-ATLAS_SINGLETOP_TCH_R_8TEV_T - - ATLAS_SINGLETOP_TCH_R_8TEV-ATLAS_SINGLETOP_TCH_R_8TEV_TB - data_uncertainties: [] + data_uncertainties: [uncertainties_TCHANNEL-XSEC.yaml] variants: legacy: data_uncertainties: - uncertainties_legacy_TCHANNEL-XSEC.yaml - data_central: data_legacy_TCHANNEL-XSEC.yaml + data_central: data_legacy_TCHANNEL-XSEC.yaml + data_central: data_TCHANNEL-XSEC.yaml ported_from: ATLAS_SINGLETOP_TCH_R_8TEV diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_28.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_28.yaml new file mode 100644 index 0000000000..d4c7059f9e --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_28.yaml @@ -0,0 +1,52 @@ +independent_variables: +- header: {name: '$|y(t)|$'} + values: + - {value: '0.0-0.3'} + - {value: '0.3-0.7'} + - {value: '0.7-1.3'} + - {value: '1.3-2.2'} +dependent_variables: +- header: {name: '$\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(t)|}$', units: '$\mathrm{pb}$'} + values: + - errors: + - {label: stat, symerror: 1.8} + - asymerror: {minus: -2.1, plus: 2.5} + label: syst + value: 32.7 + - errors: + - {label: stat, symerror: 1.8} + - asymerror: {minus: -2.4, plus: 2.2} + label: syst + value: 31.5 + - errors: + - {label: stat, symerror: 1.3} + - asymerror: {minus: -1.9, plus: 1.9} + label: syst + value: 25.3 + - errors: + - {label: stat, symerror: 0.9} + - asymerror: {minus: -1.2, plus: 1.2} + label: syst + value: 15.4 +- header: {name: '$\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(t)|}$', units: '$\mathrm{10^{-3}}$'} + values: + - errors: + - {label: stat, symerror: 35} + - asymerror: {minus: -39, plus: 47} + label: syst + value: 636 + - errors: + - {label: stat, symerror: 34} + - asymerror: {minus: -33, plus: 31} + label: syst + value: 613 + - errors: + - {label: stat, symerror: 24} + - asymerror: {minus: -27, plus: 26} + label: syst + value: 492 + - errors: + - {label: stat, symerror: 14} + - asymerror: {minus: -15, plus: 14} + label: syst + value: 299 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_29.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_29.yaml new file mode 100644 index 0000000000..b9d9d69a5c --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_29.yaml @@ -0,0 +1,52 @@ +independent_variables: +- header: {name: '$|y(t)|$'} + values: + - {value: '0.0-0.3'} + - {value: '0.3-0.7'} + - {value: '0.7-1.3'} + - {value: '1.3-2.2'} +dependent_variables: +- header: {name: '$\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(t)|}$', units: '$\mathrm{pb}$'} + values: + - errors: + - {label: stat, symerror: 1.7} + - asymerror: {minus: -1.9, plus: 1.8} + label: syst + value: 21.5 + - errors: + - {label: stat, symerror: 1.6} + - asymerror: {minus: -1.7, plus: 1.7} + label: syst + value: 18.8 + - errors: + - {label: stat, symerror: 1.2} + - asymerror: {minus: -1.6, plus: 1.6} + label: syst + value: 16.3 + - errors: + - {label: stat, symerror: 0.8} + - asymerror: {minus: -1.1, plus: 1.2} + label: syst + value: 7.0 +- header: {name: '$\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(t)|}$', units: '$\mathrm{10^{-3}}$'} + values: + - errors: + - {label: stat, symerror: 55} + - asymerror: {minus: -46, plus: 41} + label: syst + value: 714 + - errors: + - {label: stat, symerror: 53} + - asymerror: {minus: -46, plus: 46} + label: syst + value: 626 + - errors: + - {label: stat, symerror: 37} + - asymerror: {minus: -43, plus: 44} + label: syst + value: 543 + - errors: + - {label: stat, symerror: 23} + - asymerror: {minus: -29, plus: 30} + label: syst + value: 233 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_52.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_52.yaml new file mode 100644 index 0000000000..cc331fadce --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_52.yaml @@ -0,0 +1,58 @@ +independent_variables: +- header: {name: '$|y(t)|$'} + values: + - {high: 0.30, low: 0.00} + - {high: 0.70, low: 0.30} + - {high: 1.30, low: 0.70} + - {high: 2.20, low: 1.30} + - {high: 0.30, low: 0.00} + - {high: 0.70, low: 0.30} + - {high: 1.30, low: 0.70} + - {high: 2.20, low: 1.30} + - {high: 0.30, low: 0.00} + - {high: 0.70, low: 0.30} + - {high: 1.30, low: 0.70} + - {high: 2.20, low: 1.30} + - {high: 0.30, low: 0.00} + - {high: 0.70, low: 0.30} + - {high: 1.30, low: 0.70} + - {high: 2.20, low: 1.30} + +- header: {name: '$|y(t)|$'} + values: + - {high: 0.30, low: 0.00} + - {high: 0.30, low: 0.00} + - {high: 0.30, low: 0.00} + - {high: 0.30, low: 0.00} + - {high: 0.70, low: 0.30} + - {high: 0.70, low: 0.30} + - {high: 0.70, low: 0.30} + - {high: 0.70, low: 0.30} + - {high: 1.30, low: 0.70} + - {high: 1.30, low: 0.70} + - {high: 1.30, low: 0.70} + - {high: 1.30, low: 0.70} + - {high: 2.20, low: 1.30} + - {high: 2.20, low: 1.30} + - {high: 2.20, low: 1.30} + - {high: 2.20, low: 1.30} + +dependent_variables: +- header: {name: Correlation} + values: + - {value: 1.00} + - {value: -0.53} + - {value: -0.08} + - {value: -0.16} + - {value: -0.53} + - {value: 1.00} + - {value: -0.38} + - {value: -0.19} + - {value: -0.08} + - {value: -0.38} + - {value: 1.00} + - {value: -0.63} + - {value: -0.16} + - {value: -0.19} + - {value: -0.63} + - {value: 1.00} diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_53.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_53.yaml new file mode 100644 index 0000000000..ebbae65f07 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_53.yaml @@ -0,0 +1,58 @@ +independent_variables: +- header: {name: '$|y(t)|$'} + values: + - {high: 0.30, low: 0.00} + - {high: 0.70, low: 0.30} + - {high: 1.30, low: 0.70} + - {high: 2.20, low: 1.30} + - {high: 0.30, low: 0.00} + - {high: 0.70, low: 0.30} + - {high: 1.30, low: 0.70} + - {high: 2.20, low: 1.30} + - {high: 0.30, low: 0.00} + - {high: 0.70, low: 0.30} + - {high: 1.30, low: 0.70} + - {high: 2.20, low: 1.30} + - {high: 0.30, low: 0.00} + - {high: 0.70, low: 0.30} + - {high: 1.30, low: 0.70} + - {high: 2.20, low: 1.30} + +- header: {name: '$|y(t)|$'} + values: + - {high: 0.30, low: 0.00} + - {high: 0.30, low: 0.00} + - {high: 0.30, low: 0.00} + - {high: 0.30, low: 0.00} + - {high: 0.70, low: 0.30} + - {high: 0.70, low: 0.30} + - {high: 0.70, low: 0.30} + - {high: 0.70, low: 0.30} + - {high: 1.30, low: 0.70} + - {high: 1.30, low: 0.70} + - {high: 1.30, low: 0.70} + - {high: 1.30, low: 0.70} + - {high: 2.20, low: 1.30} + - {high: 2.20, low: 1.30} + - {high: 2.20, low: 1.30} + - {high: 2.20, low: 1.30} + +dependent_variables: +- header: {name: Correlation} + values: + - {value: 1.00} + - {value: -0.54} + - {value: -0.08} + - {value: -0.16} + - {value: -0.54} + - {value: 1.00} + - {value: -0.40} + - {value: -0.17} + - {value: -0.08} + - {value: -0.40} + - {value: 1.00} + - {value: -0.63} + - {value: -0.16} + - {value: -0.17} + - {value: -0.63} + - {value: 1.00} diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_77.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_77.yaml new file mode 100644 index 0000000000..adb9f0e597 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_77.yaml @@ -0,0 +1,282 @@ +dependent_variables: +- header: {name: '$\mathrm{Bin~1} \\ [\%]$'} + qualifiers: + - {name: '', value: '$\mathrm{Uncertainties~on}~\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(t)|}~\mathrm{at~parton~level}$'} + values: + - {value: ''} + - errors: + - {asymerror: {plus: +5.5, minus: -5.5}, label: stat.} + value: '' + - errors: + - {asymerror: {plus: +0.7, minus: -0.7}, label: stat.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +2.7, minus: -2.7}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.9, minus: -0.9}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.01, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.3, minus: -0.2}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +2.5, minus: -2.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.8, minus: -1.6}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.1, minus: -0.1}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.0, minus: -0.6}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +2.2, minus: -2.2}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.4, minus: -0.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +3.7, minus: -3.7}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +4.5, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.3, minus: -0.3}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.4, minus: -0.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.0, minus: -0.3}, label: syst.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +7.4, minus: -6.1}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +9.2, minus: -8.2}, label: total} + value: '' +- header: {name: '$\mathrm{Bin~2} \\ [\%]$'} + qualifiers: + - {name: '', value: '$\mathrm{Uncertainties~on}~\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(t)|}~\mathrm{at~parton~level}$'} + values: + - {value: ''} + - errors: + - {asymerror: {plus: +5.5, minus: -5.5}, label: stat.} + value: '' + - errors: + - {asymerror: {plus: +0.6, minus: -0.6}, label: stat.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +1.2, minus: -1.2}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.5, minus: -0.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.01, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.3, minus: -0.3}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +2.3, minus: -2.3}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.2, minus: -0.8}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.01, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.5, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.6, minus: -0.6}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +3.5, minus: -3.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.3, minus: -0.3}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.2, minus: -2.1}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.9, minus: -0.9}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.5, minus: -1.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.0, minus: -0.7}, label: syst.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +5.0, minus: -5.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +7.4, minus: -7.7}, label: total} + value: '' +- header: {name: '$\mathrm{Bin~3} \\ [\%]$'} + qualifiers: + - {name: '', value: '$\mathrm{Uncertainties~on}~\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(t)|}~\mathrm{at~parton~level}$'} + values: + - {value: ''} + - errors: + - {asymerror: {plus: +4.8, minus: -4.8}, label: stat.} + value: '' + - errors: + - {asymerror: {plus: +0.6, minus: -0.6}, label: stat.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +0.2, minus: -0.2}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.3, minus: -0.3}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.01, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.2, minus: -0.2}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.1, minus: -1.1}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.8, minus: -0.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.01, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.0, minus: -0.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +2.2, minus: -2.2}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +4.0, minus: -4.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.6, minus: -1.6}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.0, minus: -1.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.8, minus: -0.8}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.8, minus: -0.8}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.6, minus: -0.0}, label: syst.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +5.2, minus: -5.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +7.1, minus: -7.2}, label: total} + value: '' +- header: {name: '$\mathrm{Bin~4} \\ [\%]$'} + qualifiers: + - {name: '', value: '$\mathrm{Uncertainties~on}~\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(t)|}~\mathrm{at~parton~level}$'} + values: + - {value: ''} + - errors: + - {asymerror: {plus: +4.8, minus: -4.8}, label: stat.} + value: '' + - errors: + - {asymerror: {plus: +0.7, minus: -0.7}, label: stat.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +3.4, minus: -3.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.8, minus: -0.8}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.01, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.5, minus: -0.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.9, minus: -0.9}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.2, minus: -1.1}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.1, minus: -0.1}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.5, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.5, minus: -1.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.7, minus: -1.7}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.7, minus: -0.7}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.0, minus: -1.3}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.4, minus: -1.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.1, minus: -0.1}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.4, minus: -0.5}, label: syst.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +4.8, minus: -4.9}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +6.8, minus: -6.9}, label: total} + value: '' +independent_variables: +- header: {name: '$\mathrm{Source~of~uncertainty}$'} + values: + - {value: ''} + - {value: '$\mathrm{Data~statistics}$'} + - {value: '$\mathrm{Monte~Carlo~statistics}$'} + - {value: ''} + - {value: '$\mathrm{Unfolding}$'} + - {value: '$\mathrm{Background~normalisation}$'} + - {value: '$\mathrm{Background~modelling}$'} + - {value: '$\mathrm{Lepton~reconstruction}$'} + - {value: '$\mathrm{Jet~reconstruction}$'} + - {value: '$\mathrm{Jet~energy~scale}$'} + - {value: '$\mathrm{Flavour~tagging}$'} + - {value: '$E_{\mathrm{T}}^{\mathrm{miss}}~\mathrm{modelling}$'} + - {value: '$\mathrm{PDF}$'} + - {value: '$tq~(\bar tq)~\mathrm{NLO~matching}$'} + - {value: '$tq~(\bar tq)~\mathrm{parton~shower}$'} + - {value: '$tq~(\bar tq)~\mathrm{scale~variations}$'} + - {value: '$t \bar t~\mathrm{NLO~matching}$'} + - {value: '$t \bar t~\mathrm{parton~shower}$'} + - {value: '$t \bar t~\mathrm{scale~variations}$'} + - {value: ''} + - {value: '$\mathrm{Total~systematic}$'} + - {value: '$\mathrm{Total~(stat.~+~syst.)}$'} diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_79.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_79.yaml new file mode 100644 index 0000000000..2a443d223c --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_79.yaml @@ -0,0 +1,282 @@ +dependent_variables: +- header: {name: '$\mathrm{Bin~1} \\ [\%]$'} + qualifiers: + - {name: '', value: '$\mathrm{Uncertainties~on}~\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(t)|}~\mathrm{at~parton~level}$'} + values: + - {value: ''} + - errors: + - {asymerror: {plus: +7.7, minus: -7.7}, label: stat.} + value: '' + - errors: + - {asymerror: {plus: +0.8, minus: -0.8}, label: stat.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +0.5, minus: -0.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.7, minus: -1.7}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.8, minus: -0.8}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.3, minus: -0.2}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.3, minus: -0.3}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.6, minus: -0.9}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.01, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.3, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.9, minus: -0.9}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +5.0, minus: -5.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.5, minus: -0.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.0, minus: -2.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.7, minus: -0.7}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.4, minus: -1.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.0, minus: -0.9}, label: syst.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +5.8, minus: -6.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +10, minus: -10}, label: total} + value: '' +- header: {name: '$\mathrm{Bin~2} \\ [\%]$'} + qualifiers: + - {name: '', value: '$\mathrm{Uncertainties~on}~\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(t)|}~\mathrm{at~parton~level}$'} + values: + - {value: ''} + - errors: + - {asymerror: {plus: +8.5, minus: -8.5}, label: stat.} + value: '' + - errors: + - {asymerror: {plus: +0.7, minus: -0.7}, label: stat.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +0.1, minus: -0.1}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.5, minus: -0.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.4, minus: -0.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.6, minus: -0.3}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.1, minus: -0.1}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.0, minus: -0.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.01, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.3, minus: -0.6}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +3.5, minus: -3.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +4.5, minus: -4.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +3.0, minus: -3.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.2, minus: -0.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.0, minus: -1.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +3.2, minus: -3.2}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.0, minus: -0.2}, label: syst.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +7.4, minus: -7.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +11, minus: -11}, label: total} + value: '' +- header: {name: '$\mathrm{Bin~3} \\ [\%]$'} + qualifiers: + - {name: '', value: '$\mathrm{Uncertainties~on}~\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(t)|}~\mathrm{at~parton~level}$'} + values: + - {value: ''} + - errors: + - {asymerror: {plus: +6.9, minus: -6.9}, label: stat.} + value: '' + - errors: + - {asymerror: {plus: +0.6, minus: -0.6}, label: stat.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +3.1, minus: -3.1}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.4, minus: -0.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +2.0, minus: -2.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.2, minus: -0.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.7, minus: -0.7}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.1, minus: -1.3}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.01, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.1, minus: -0.6}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.1, minus: -1.2}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +5.7, minus: -5.7}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.6, minus: -1.6}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +3.3, minus: -1.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.8, minus: -1.8}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.0, minus: -1.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.0, minus: -1.9}, label: syst.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +8.1, minus: -7.9}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +11, minus: -10}, label: total} + value: '' +- header: {name: '$\mathrm{Bin~4} \\ [\%]$'} + qualifiers: + - {name: '', value: '$\mathrm{Uncertainties~on}~\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(t)|}~\mathrm{at~parton~level}$'} + values: + - {value: ''} + - errors: + - {asymerror: {plus: +10, minus: -10}, label: stat.} + value: '' + - errors: + - {asymerror: {plus: +0.9, minus: -0.9}, label: stat.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +4.7, minus: -4.7}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +3.1, minus: -3.1}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +3.4, minus: -3.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.4, minus: -0.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.6, minus: -0.6}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +1.7, minus: -0.7}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.1, minus: -0.1}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.9, minus: -0.0}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +5.2, minus: -5.2}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +8.4, minus: -8.4}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.8, minus: -0.8}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +3.2, minus: -2.5}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.9, minus: -0.9}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +0.8, minus: -0.8}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +3.7, minus: -0.0}, label: syst.} + value: '' + - {value: ''} + - errors: + - {asymerror: {plus: +13, minus: -12}, label: syst.} + value: '' + - errors: + - {asymerror: {plus: +16, minus: -16}, label: total} + value: '' +independent_variables: +- header: {name: '$\mathrm{Source~of~uncertainty}$'} + values: + - {value: ''} + - {value: '$\mathrm{Data~statistics}$'} + - {value: '$\mathrm{Monte~Carlo~statistics}$'} + - {value: ''} + - {value: '$\mathrm{Unfolding}$'} + - {value: '$\mathrm{Background~normalisation}$'} + - {value: '$\mathrm{Background~modelling}$'} + - {value: '$\mathrm{Lepton~reconstruction}$'} + - {value: '$\mathrm{Jet~reconstruction}$'} + - {value: '$\mathrm{Jet~energy~scale}$'} + - {value: '$\mathrm{Flavour~tagging}$'} + - {value: '$E_{\mathrm{T}}^{\mathrm{miss}}~\mathrm{modelling}$'} + - {value: '$\mathrm{PDF}$'} + - {value: '$tq~(\bar tq)~\mathrm{NLO~matching}$'} + - {value: '$tq~(\bar tq)~\mathrm{parton~shower}$'} + - {value: '$tq~(\bar tq)~\mathrm{scale~variations}$'} + - {value: '$t \bar t~\mathrm{NLO~matching}$'} + - {value: '$t \bar t~\mathrm{parton~shower}$'} + - {value: '$t \bar t~\mathrm{scale~variations}$'} + - {value: ''} + - {value: '$\mathrm{Total~systematic}$'} + - {value: '$\mathrm{Total~(stat.~+~syst.)}$'} diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_8.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_8.yaml new file mode 100644 index 0000000000..a0936ca8af --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/rawdata/HEPData-ins1512776-v2-Table_8.yaml @@ -0,0 +1,84 @@ +dependent_variables: +- header: {name: '$~~$'} + values: + - errors: + - {label: stat., symerror: 0.9} + - {label: exp., symerror: 2.7} + - asymerror: {minus: -1.7, plus: 2.7} + label: scale + - {label: PDF, symerror: 0.4} + - {label: NLO - matching method, symerror: 1.0} + - {label: parton shower, symerror: 1.1} + - {label: lumi. pb, symerror: 1.1} + value: 56.7 + - errors: + - asymerror: {minus: -3.8, plus: 4.3} + label: pb + value: 56.7 + - errors: + - {label: stat., symerror: 0.8} + - {label: exp., symerror: 2.3} + - asymerror: {minus: -0.8, plus: 1.4} + label: scale + - {label: PDF, symerror: 0.3} + - {label: '', symerror: ' '} + - asymerror: {minus: -0.6, plus: 0.7} + label: NLO-matching method + - {label: parton shower, symerror: 0.6} + - {label: lumi. pb, symerror: 0.6} + value: 32.9 + - errors: + - asymerror: {minus: -2.7, plus: 3.0} + label: pb + value: 32.9 + - errors: + - {label: stat., symerror: 1.2} + - {label: exp., symerror: 5.1} + - asymerror: {minus: -2.5, plus: 4.1} + label: scale + - {label: PDF, symerror: 0.7} + - {label: '', symerror: ''} + - asymerror: {minus: -1.6, plus: 1.7} + label: NLO-matching method + - {label: parton shower, symerror: 1.6} + - {label: lumi. pb, symerror: 1.7} + value: 89.6 + - errors: + - asymerror: {minus: -6.3, plus: 7.1} + label: pb + value: 89.6 + - errors: + - {label: stat., symerror: 0.05} + - {label: exp., symerror: 0.07} + value: 1.72 + - errors: + - {label: '', symerror: 0.09} + value: 1.72 + - errors: + - {label: stat., symerror: 0.007} + - {label: exp., symerror: 0.029} + - asymerror: {minus: -0.014, plus: 0.023} + label: scale + - {label: PDF, symerror: 0.004} + - {label: NLO-matching method, symerror: 0.010} + - {label: parton shower, symerror: 0.009} + - {label: lumi., symerror: 0.010} + - {label: '$m_{t}$', symerror: 0.005} + - {label: theor., symerror: 0.024} + value: 1.029 + - errors: + - {label: '', symerror: 0.048} + value: 1.029 +independent_variables: +- header: {name: '$~~$'} + values: + - {value: '$\sigma_{\mathrm{tot}}(tq)$'} + - {value: '$\sigma_{\mathrm{tot}}(tq)$'} + - {value: '$\sigma_{\mathrm{tot}}(\bar tq)$'} + - {value: '$\sigma_{\mathrm{tot}}(\bar tq)$'} + - {value: '$\sigma_{\mathrm{tot}}(tq ~$+$~\bar tq)$'} + - {value: '$\sigma_{\mathrm{tot}}(tq ~$+$~\bar tq)$'} + - {value: '$\mathrm{R}_{\mathrm{t}}~=~\dfrac{\sigma_{\mathrm{tot}}(tq)}{\sigma_{\mathrm{tot}}(\bar tq)}$'} + - {value: '$\mathrm{R}_{\mathrm{t}}~=~\dfrac{\sigma_{\mathrm{tot}}(tq)}{\sigma_{\mathrm{tot}}(\bar tq)}$'} + - {value: '$f_{\mathrm{LV}}~\dot~|V_{tb}|$'} + - {value: '$f_{\mathrm{LV}}~\dot~|V_{tb}|$'} diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/uncertainties_T-RAP-NORM.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/uncertainties_T-RAP-NORM.yaml new file mode 100644 index 0000000000..2f72698a9c --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/uncertainties_T-RAP-NORM.yaml @@ -0,0 +1,135 @@ +definitions: + sys_corr_1: + description: Statistical uncertainty 0 + treatment: ADD + type: CORR + sys_corr_2: + description: Statistical uncertainty 1 + treatment: ADD + type: CORR + sys_corr_3: + description: Statistical uncertainty 2 + treatment: ADD + type: CORR + sys_corr_4: + description: $\mathrm{Monte~Carlo~statistics}$ + treatment: MULT + type: CORR + sys_corr_5: + description: $\mathrm{Unfolding}$ + treatment: MULT + type: CORR + sys_corr_6: + description: $\mathrm{Background~normalisation}$ + treatment: MULT + type: CORR + sys_corr_7: + description: $\mathrm{Background~modelling}$ + treatment: MULT + type: CORR + sys_corr_8: + description: $\mathrm{Lepton~reconstruction}$ + treatment: MULT + type: CORR + sys_corr_9: + description: $\mathrm{Jet~reconstruction}$ + treatment: MULT + type: CORR + sys_corr_10: + description: $\mathrm{Jet~energy~scale}$ + treatment: MULT + type: CORR + sys_corr_11: + description: $\mathrm{Flavour~tagging}$ + treatment: MULT + type: CORR + sys_corr_12: + description: $E_{\mathrm{T}}^{\mathrm{miss}}~\mathrm{modelling}$ + treatment: MULT + type: CORR + sys_corr_13: + description: $\mathrm{PDF}$ + treatment: MULT + type: CORR + sys_corr_14: + description: $tq~(\bar tq)~\mathrm{NLO~matching}$ + treatment: MULT + type: CORR + sys_corr_15: + description: $tq~(\bar tq)~\mathrm{parton~shower}$ + treatment: MULT + type: CORR + sys_corr_16: + description: $tq~(\bar tq)~\mathrm{scale~variations}$ + treatment: MULT + type: CORR + sys_corr_17: + description: $t \bar t~\mathrm{NLO~matching}$ + treatment: MULT + type: CORR + sys_corr_18: + description: $t \bar t~\mathrm{parton~shower}$ + treatment: MULT + type: CORR + sys_corr_19: + description: $t \bar t~\mathrm{scale~variations}$ + treatment: MULT + type: CORR +bins: +- sys_corr_1: -3.16176678e-02 + sys_corr_2: 1.55853185e-02 + sys_corr_3: 6.75286960e-03 + sys_corr_4: 4.56797460e-03 + sys_corr_5: 1.76193306e-02 + sys_corr_6: 5.87311020e-03 + sys_corr_7: 5.65140292e-05 + sys_corr_8: 1.69542088e-03 + sys_corr_9: 1.63141950e-02 + sys_corr_10: 8.65729017e-03 + sys_corr_11: 6.52567800e-04 + sys_corr_12: 3.39084175e-03 + sys_corr_13: 1.43564916e-02 + sys_corr_14: 2.61027120e-03 + sys_corr_15: 2.41450086e-02 + sys_corr_16: 2.54313132e-02 + sys_corr_17: 1.95770340e-03 + sys_corr_18: 2.61027120e-03 + sys_corr_19: 1.69542088e-03 +- sys_corr_1: 2.88358828e-02 + sys_corr_2: 1.37373476e-02 + sys_corr_3: 9.46498030e-03 + sys_corr_4: 3.63423180e-03 + sys_corr_5: 7.26846360e-03 + sys_corr_6: 3.02852650e-03 + sys_corr_7: 5.24556177e-05 + sys_corr_8: 1.81711590e-03 + sys_corr_9: 1.39312219e-02 + sys_corr_10: 6.29467412e-03 + sys_corr_11: 5.24556177e-05 + sys_corr_12: 2.62278089e-03 + sys_corr_13: 3.63423180e-03 + sys_corr_14: 2.11996855e-02 + sys_corr_15: 1.81711590e-03 + sys_corr_16: 1.07117403e-02 + sys_corr_17: 5.45134770e-03 + sys_corr_18: 9.08557950e-03 + sys_corr_19: 3.67189324e-03 +- sys_corr_1: -4.94051825e-03 + sys_corr_2: -1.95612839e-02 + sys_corr_3: 1.20272958e-02 + sys_corr_4: 2.93605920e-03 + sys_corr_5: 9.78686400e-04 + sys_corr_6: 1.46802960e-03 + sys_corr_7: 4.23783642e-05 + sys_corr_8: 9.78686400e-04 + sys_corr_9: 5.38277520e-03 + sys_corr_10: 3.24593558e-03 + sys_corr_11: 4.23783642e-05 + sys_corr_12: 1.69513457e-03 + sys_corr_13: 1.07655504e-02 + sys_corr_14: 1.95737280e-02 + sys_corr_15: 7.82949120e-03 + sys_corr_16: 5.93297099e-03 + sys_corr_17: 3.91474560e-03 + sys_corr_18: 3.91474560e-03 + sys_corr_19: 2.54270185e-03 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/uncertainties_TBAR-RAP-NORM.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/uncertainties_TBAR-RAP-NORM.yaml new file mode 100644 index 0000000000..ffb7fe8bfd --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/uncertainties_TBAR-RAP-NORM.yaml @@ -0,0 +1,135 @@ +definitions: + sys_corr_1: + description: Statistical uncertainty 0 + treatment: ADD + type: CORR + sys_corr_2: + description: Statistical uncertainty 1 + treatment: ADD + type: CORR + sys_corr_3: + description: Statistical uncertainty 2 + treatment: ADD + type: CORR + sys_corr_4: + description: $\mathrm{Monte~Carlo~statistics}$ + treatment: MULT + type: CORR + sys_corr_5: + description: $\mathrm{Unfolding}$ + treatment: MULT + type: CORR + sys_corr_6: + description: $\mathrm{Background~normalisation}$ + treatment: MULT + type: CORR + sys_corr_7: + description: $\mathrm{Background~modelling}$ + treatment: MULT + type: CORR + sys_corr_8: + description: $\mathrm{Lepton~reconstruction}$ + treatment: MULT + type: CORR + sys_corr_9: + description: $\mathrm{Jet~reconstruction}$ + treatment: MULT + type: CORR + sys_corr_10: + description: $\mathrm{Jet~energy~scale}$ + treatment: MULT + type: CORR + sys_corr_11: + description: $\mathrm{Flavour~tagging}$ + treatment: MULT + type: CORR + sys_corr_12: + description: $E_{\mathrm{T}}^{\mathrm{miss}}~\mathrm{modelling}$ + treatment: MULT + type: CORR + sys_corr_13: + description: $\mathrm{PDF}$ + treatment: MULT + type: CORR + sys_corr_14: + description: $tq~(\bar tq)~\mathrm{NLO~matching}$ + treatment: MULT + type: CORR + sys_corr_15: + description: $tq~(\bar tq)~\mathrm{parton~shower}$ + treatment: MULT + type: CORR + sys_corr_16: + description: $tq~(\bar tq)~\mathrm{scale~variations}$ + treatment: MULT + type: CORR + sys_corr_17: + description: $t \bar t~\mathrm{NLO~matching}$ + treatment: MULT + type: CORR + sys_corr_18: + description: $t \bar t~\mathrm{parton~shower}$ + treatment: MULT + type: CORR + sys_corr_19: + description: $t \bar t~\mathrm{scale~variations}$ + treatment: MULT + type: CORR +bins: +- sys_corr_1: -4.58704797e-02 + sys_corr_2: 2.60138756e-02 + sys_corr_3: 1.11947720e-02 + sys_corr_4: 5.60090160e-03 + sys_corr_5: 3.50056350e-03 + sys_corr_6: 1.19019159e-02 + sys_corr_7: 5.60090160e-03 + sys_corr_8: 1.81894615e-03 + sys_corr_9: 2.10033810e-03 + sys_corr_10: 5.45683845e-03 + sys_corr_11: 6.06315384e-05 + sys_corr_12: 1.81894615e-03 + sys_corr_13: 6.30101430e-03 + sys_corr_14: 3.50056350e-02 + sys_corr_15: 3.50056350e-03 + sys_corr_16: 1.51578846e-02 + sys_corr_17: 4.90078890e-03 + sys_corr_18: 9.80157780e-03 + sys_corr_19: 5.45683845e-03 +- sys_corr_1: 4.77687935e-02 + sys_corr_2: 1.84667968e-02 + sys_corr_3: 1.44466958e-02 + sys_corr_4: 4.38221910e-03 + sys_corr_5: 6.26031300e-04 + sys_corr_6: 3.13015650e-03 + sys_corr_7: 2.50412520e-03 + sys_corr_8: 3.11446639e-03 + sys_corr_9: 6.26031300e-04 + sys_corr_10: 5.19077732e-03 + sys_corr_11: 5.42159009e-05 + sys_corr_12: 3.11446639e-03 + sys_corr_13: 2.19110955e-02 + sys_corr_14: 2.81714085e-02 + sys_corr_15: 1.87809390e-02 + sys_corr_16: 2.56214352e-03 + sys_corr_17: 6.26031300e-03 + sys_corr_18: 2.00330016e-02 + sys_corr_19: 1.08431802e-03 +- sys_corr_1: -9.92291305e-03 + sys_corr_2: -3.13549407e-02 + sys_corr_3: 1.77963533e-02 + sys_corr_4: 3.25164690e-03 + sys_corr_5: 1.68001756e-02 + sys_corr_6: 2.16776460e-03 + sys_corr_7: 1.08388230e-02 + sys_corr_8: 1.79741545e-03 + sys_corr_9: 3.79358805e-03 + sys_corr_10: 5.96135265e-03 + sys_corr_11: 4.69334803e-05 + sys_corr_12: 2.69612318e-03 + sys_corr_13: 6.24409344e-03 + sys_corr_14: 3.08906456e-02 + sys_corr_15: 8.67105840e-03 + sys_corr_16: 1.46098127e-02 + sys_corr_17: 9.75494070e-03 + sys_corr_18: 5.41941150e-03 + sys_corr_19: 8.91736126e-03 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/uncertainties_TCHANNEL-XSEC.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/uncertainties_TCHANNEL-XSEC.yaml new file mode 100644 index 0000000000..02e1abdba8 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/uncertainties_TCHANNEL-XSEC.yaml @@ -0,0 +1,12 @@ +definitions: + stat: + description: Uncorrelated statistical uncertainties + treatment: ADD + type: UNCORR + sys_corr_1: + description: 'Sys uncertainty idx: 1' + treatment: MULT + type: UNCORR +bins: +- stat: 0.05 + sys_corr_1: 0.07 diff --git a/validphys2/src/validphys/commondataparser.py b/validphys2/src/validphys/commondataparser.py index 6ed4779250..eb38a14ede 100644 --- a/validphys2/src/validphys/commondataparser.py +++ b/validphys2/src/validphys/commondataparser.py @@ -112,6 +112,7 @@ def _quick_yaml_load(filepath): "SHP_ASY": ("$\\eta$", "$p_T (GeV)$", "$\\sqrt{s} (GeV)$"), "JET_POL": ("$\\eta$", "$p_T^2 (GeV^2)$", "$\\sqrt{s} (GeV)$"), "DIJET_POL": ("$\\m_{1,2} (GeV)", "$\\eta_1$", "$\\eta_2$"), + "SINGLETOP": ("$y$", "$m_t^2 (GeV^2)$", "$\\sqrt{s} (GeV)$"), } PROCESS_DESCRIPTION_LABEL = { diff --git a/validphys2/src/validphys/process_options.py b/validphys2/src/validphys/process_options.py index 53016409b5..a6ece39ee8 100644 --- a/validphys2/src/validphys/process_options.py +++ b/validphys2/src/validphys/process_options.py @@ -213,7 +213,6 @@ def _inc_xq2map(kin_info): # Compute x, Q2 # k2 necessary to take the mass for DY inclusive cross sections still not migrated mass2 = kin_info.get_one_of(_Vars.m_W2, _Vars.m_Z2, _Vars.m_t2, "k2") - return np.sqrt(mass2) / kin_info[_Vars.sqrts], mass2 @@ -274,6 +273,20 @@ def _dybosonptrap_xq2map(kin_info): x = np.concatenate((x1, x2)) return x, np.concatenate((ET2, ET2)) +def _singletop_xq2map(kin_dict): + + y_t = kin_dict[_Vars.y_t] + sqrts = kin_dict[_Vars.sqrts] + m_t2 = kin_dict[_Vars.m_t2] + + q2 = m_t2 + ratio = np.sqrt(q2) / sqrts + x1 = ratio * np.exp(y_t) + x2 = ratio * np.exp(-y_t) + x = np.concatenate((x1, x2)) + return np.clip(x, a_min=None, a_max=1, out=x), np.concatenate((q2, q2)) + + DIS = _Process( "DIS", @@ -404,6 +417,14 @@ def _dybosonptrap_xq2map(kin_info): "POS_DIS", "Positivity of F2 structure functions", accepted_variables=(_Vars.x, _Vars.Q2) ) +SINGLETOP = _Process( + "SINGLETOP", + "Single top production", + accepted_variables=(_Vars.m_t2, _Vars.sqrts, _Vars.y_t, _Vars.pT_t), + xq2map_function=_singletop_xq2map, +) + + PROCESSES = { "DIS": DIS, "DIS_NC": dataclasses.replace(DIS, name="DIS_NC"), @@ -433,6 +454,7 @@ def _dybosonptrap_xq2map(kin_info): "DY_NC_PTRAP": dataclasses.replace(DY_PT_RAP, name="DY_NC_PTRAP", description="DY Z (ll) + j"), "POS_XPDF": POS_XPDF, "POS_DIS": POS_DIS, + "SINGLETOP": SINGLETOP, } From 51792f18de07abb2fe59d6c8fae65d44b1816ed0 Mon Sep 17 00:00:00 2001 From: jacoterh <54140851+jacoterh@users.noreply.github.com> Date: Mon, 25 Nov 2024 12:01:37 +0000 Subject: [PATCH 2/2] implementing last comments --- .../commondata/ATLAS_SINGLETOP_13TEV/filter.py | 6 +----- .../kinematics_TCHANNEL-XSEC.yaml | 10 +--------- .../commondata/ATLAS_SINGLETOP_13TEV/metadata.yaml | 11 ----------- .../data_legacy_TCHANNEL-XSEC.yaml | 2 -- .../commondata/ATLAS_SINGLETOP_7TEV/metadata.yaml | 1 - .../data_legacy_TCHANNEL-XSEC.yaml | 2 -- .../commondata/ATLAS_SINGLETOP_8TEV/metadata.yaml | 1 - 7 files changed, 2 insertions(+), 31 deletions(-) delete mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_legacy_TCHANNEL-XSEC.yaml delete mode 100644 nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_legacy_TCHANNEL-XSEC.yaml diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/filter.py b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/filter.py index 6e8f8ae2c0..e265ebc4a3 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/filter.py +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/filter.py @@ -72,11 +72,7 @@ def get_kinematics(hepdata: dict, bin_index: list = [], indx: int = 0) -> list: } kinematics.append(kin_value) else: # inclusive case - kin_value = { - "zero": {"min": 0, "mid": 0, "max": 0}, - "m_t2": {"min": None, "mid": MT_VALUE**2, "max": None}, - "sqrts": {"min": None, "mid": SQRT_S, "max": None}, - } + kin_value = {"m_t2": {"min": None, "mid": MT_VALUE**2, "max": None}} kinematics.append(kin_value) return kinematics diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/kinematics_TCHANNEL-XSEC.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/kinematics_TCHANNEL-XSEC.yaml index b6a60bc58d..d73682d60a 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/kinematics_TCHANNEL-XSEC.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/kinematics_TCHANNEL-XSEC.yaml @@ -1,13 +1,5 @@ bins: -- zero: - min: 0 - mid: 0 - max: 0 - m_t2: +- m_t2: min: null mid: 29756.25 max: null - sqrts: - min: null - mid: 13000.0 - max: null diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/metadata.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/metadata.yaml index 74e215a83b..7f73b7bbb6 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/metadata.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_13TEV/metadata.yaml @@ -27,23 +27,13 @@ implemented_observables: dataset_label: ATLAS single top $R_{t}$ 13 TeV plot_x: idat kinematic_coverage: - - zero - m_t2 - - sqrts kinematics: variables: - zero: - description: "" - label: $$ - units: "" m_t2: description: "top mass squared" label: $M^2$ units: $GeV^2$ - sqrts: - description: "Center of mass energy" - label: $\sqrt{s}$ - units: $GeV$ file: kinematics_TCHANNEL-XSEC.yaml theory: conversion_factor: 1.0 @@ -56,6 +46,5 @@ implemented_observables: legacy: data_uncertainties: - uncertainties_legacy_TCHANNEL-XSEC.yaml - data_central: data_legacy_TCHANNEL-XSEC.yaml data_central: data_TCHANNEL-XSEC.yaml ported_from: ATLAS_SINGLETOP_TCH_R_13TEV diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_legacy_TCHANNEL-XSEC.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_legacy_TCHANNEL-XSEC.yaml deleted file mode 100644 index da87d41bb0..0000000000 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/data_legacy_TCHANNEL-XSEC.yaml +++ /dev/null @@ -1,2 +0,0 @@ -data_central: -- 2.04 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/metadata.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/metadata.yaml index a94a3a7ca5..81c9c9fc4e 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/metadata.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_7TEV/metadata.yaml @@ -166,6 +166,5 @@ implemented_observables: legacy: data_uncertainties: - uncertainties_legacy_TCHANNEL-XSEC.yaml - data_central: data_legacy_TCHANNEL-XSEC.yaml data_central: data_TCHANNEL-XSEC.yaml ported_from: ATLAS_SINGLETOP_TCH_R_7TEV diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_legacy_TCHANNEL-XSEC.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_legacy_TCHANNEL-XSEC.yaml deleted file mode 100644 index df9c7c392f..0000000000 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/data_legacy_TCHANNEL-XSEC.yaml +++ /dev/null @@ -1,2 +0,0 @@ -data_central: -- 1.72 diff --git a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/metadata.yaml b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/metadata.yaml index c0080f554c..9b8dc953da 100644 --- a/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/metadata.yaml +++ b/nnpdf_data/nnpdf_data/commondata/ATLAS_SINGLETOP_8TEV/metadata.yaml @@ -166,6 +166,5 @@ implemented_observables: legacy: data_uncertainties: - uncertainties_legacy_TCHANNEL-XSEC.yaml - data_central: data_legacy_TCHANNEL-XSEC.yaml data_central: data_TCHANNEL-XSEC.yaml ported_from: ATLAS_SINGLETOP_TCH_R_8TEV