diff --git a/buildmaster/LHCB_DY_8TEV_MUON/data.yaml b/buildmaster/LHCB_DY_8TEV_MUON/data.yaml new file mode 100644 index 0000000000..15de29fca1 --- /dev/null +++ b/buildmaster/LHCB_DY_8TEV_MUON/data.yaml @@ -0,0 +1,35 @@ +data_central: +- 1223.0 +- 3263.0 +- 4983.0 +- 6719.0 +- 8051.0 +- 8967.0 +- 9561.0 +- 9822.0 +- 9721.0 +- 9030.0 +- 7748.0 +- 6059.0 +- 4385.0 +- 2724.0 +- 1584.0 +- 749.0 +- 383.0 +- 11.0 +- 236500.0 +- 208400.0 +- 182000.0 +- 153300.0 +- 119500.0 +- 84400.0 +- 86400.0 +- 23000.0 +- 134000.0 +- 119800.0 +- 110600.0 +- 102400.0 +- 92500.0 +- 79900.0 +- 119300.0 +- 60000.0 diff --git a/buildmaster/LHCB_DY_8TEV_MUON/filter.py b/buildmaster/LHCB_DY_8TEV_MUON/filter.py new file mode 100644 index 0000000000..20de5b4125 --- /dev/null +++ b/buildmaster/LHCB_DY_8TEV_MUON/filter.py @@ -0,0 +1,387 @@ +import numpy as np +import pandas as pd +import pathlib +import yaml + +from validphys.commondata_utils import covmat_to_artunc + + +MZ_VALUE = 91.1876 # GeV +MW_VALUE = 80.398 # GeV +SQRT_S = 7_000.0 # GeV +NORM_FACTOR = 1_000.0 # from pb -> fb + +# Correct tables to read values [[Z], [W+, W-]] +TABLES = {2: [0], 1: [0, 2]} # {table_id: [indexes]} +NB_POINTS = [18, 8, 8] # [N_Z, N_W+, N_W-] + +# MAP Boson string to M values +MAP_BOSON = {"Z": MZ_VALUE, "W": MW_VALUE} +MAP_TABLE = {2: "Z", 1: "W"} + + +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-ins1406555-v{version}-Table_{table_id}" + table = pathlib.Path(f"./rawdata/{filename}.yaml") + + return yaml.safe_load(table.read_text()) + + +def read_metadata() -> tuple[int, int, list]: + """Read the version and list of tables from metadata. + + Returns + ------- + tuple(int, list): + data version and list of hepdata tables + + """ + metadata = pathlib.Path("./metadata.yaml") + content = yaml.safe_load(metadata.read_text()) + + version = content["hepdata"]["version"] + nb_datapoints = sum(content["implemented_observables"][0]["npoints"]) + tables = content["implemented_observables"][0]["tables"] + + return version, nb_datapoints, tables + + +def get_kinematics(hepdata: dict, bin_index: list, boson: str = "Z") -> 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 + + Returns + ------- + tuple(int, list): + data version and list of hepdata tables + + """ + rapbins = hepdata["independent_variables"][0]["values"] + + kinematics = [] + for bins in bin_index: + ymin = float(rapbins[bins]["low"]) + ymax = float(rapbins[bins]["high"]) + kin_value = { + "y": { + "min": ymin, + "mid": 0.5 * (ymin + ymax), + "max": ymax, + }, + "M2": {"min": None, "mid": MAP_BOSON[boson] ** 2, "max": None}, + "sqrt_s": {"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 + list of Non-empty bin index + idx: int + 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 [NORM_FACTOR * central[i]["value"] for i in bin_index] + + +def get_errors(hepdata: dict, bin_index: list, indx: int = 0) -> dict: + """Extract the error values from the HepData yaml file. + + Parameters + ---------- + hepdata: dict + dictionary containing all data info + bin_index: list + list of Non-empty bin index + indx: int + index from which the errors will be read + + Returns + ------- + list: + list of dictionaries whose contents are the various + source of uncertainties + + """ + errors = hepdata["dependent_variables"][indx]["values"] + + stat, sys_corr, sys_beam, sys_lumi = [], [], [], [] + for idx in bin_index: + stat.append(NORM_FACTOR * errors[idx]["errors"][0]["symerror"]) + sys_corr.append(NORM_FACTOR * errors[idx]["errors"][1]["symerror"]) + sys_beam.append(NORM_FACTOR * errors[idx]["errors"][2]["symerror"]) + sys_lumi.append(NORM_FACTOR * errors[idx]["errors"][3]["symerror"]) + + return { + "stat": stat, + "sys_corr": sys_corr, + "sys_beam": sys_beam, + "sys_lumi": sys_lumi, + } + + +def read_corrmatrix(nb_datapoints: int) -> np.ndarray: + """Read the matrix and returns a symmetrized verions. + + Parameters + ---------- + nb_datapoints: int + total number of datapoints + + Returns + ------- + np.ndarray: (nb_datapoints, nb_datapoints) + a symmetrized square matrix + + """ + corrmat = pd.read_csv( + "./rawdata/corrmat.corr", + names=[f'{i}' for i in range(nb_datapoints)], + delim_whitespace=True, + ) + corrmat = corrmat.iloc[:, :].values + + # Symmetrize the Correlation Matrix + symmetrised = np.triu(corrmat.T) + np.triu(corrmat.T, 1).T + assert np.allclose(symmetrised, symmetrised.T) + + return symmetrised + + +def multiply_syst(covmat: np.ndarray, sys: list) -> np.ndarray: + """Multiply the covmat with the total systematics. + + Parameters + ---------- + covmat: np.ndarray + a squared matrix representing the covmat + sys: list + a vector containing the correlatdd systematics + + Returns + ------- + np.ndarray: + covariance matrix multiplied by the total systematics + + """ + covmat_sys = np.zeros(covmat.shape) + + for i in range(covmat.shape[0]): + for j in range(covmat.shape[-1]): + covmat_sys[i][j] = covmat[i][j] * sys[i] * sys[j] + + return covmat_sys.flatten() + + +def generate_artificial_unc(**kwargs) -> np.ndarray: + """A wrapper around `covmat_to_artunc` to return an array. + + Returns + ------- + np.ndarray: + a squared matrix with artificial uncertainties + """ + artunc = covmat_to_artunc(**kwargs) + return np.array(artunc) + + +def concatenate_dicts(multidict: list[dict]) -> dict: + """Given a list of dictionaries combined the values. + + Parameters + ---------- + multidict: list[dict] + list of dictionaries with the same keys + + Returns + ------- + dict: + dictionary whose keys are combined + """ + new_dict = {} + for key in multidict[0].keys(): + new_dict[key] = [] + for element in multidict: + new_dict[key] += element[key] + + return new_dict + + +def format_uncertainties(uncs: dict, artunc: np.ndarray) -> list: + """Format the uncertainties to be dumped into the yaml file. + + Parameters + ---------- + uncs: dict + dictionary containing the various source of uncertainties + artunc: no.ndarray + array of the artificial systematic uncertainties + + Returns + ------- + list: + list of ditionaries whose elements are the various errors + + """ + combined_errors = [] + + artunc = artunc.tolist() + for idat in range(len(artunc)): + error_value = {} + for jdat in range(len(artunc[idat])): + error_value[f"sys_corr_{jdat + 1}"] = artunc[idat][jdat] + + error_value["stat"] = uncs["stat"][idat] + error_value["sys_beam"] = uncs["sys_beam"][idat] + error_value["sys_luminosity"] = uncs["sys_lumi"][idat] + combined_errors.append(error_value) + + return combined_errors + + +def dump_commondata(kinematics: list, data: list, errors: list) -> 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: list + list containing the different errors + + """ + error_definition = { + f"sys_corr_{i + 1}": { + "description": "Correlated systematic uncertainties", + "treatment": "ADD", + "type": "CORR", + } + for i in range(len(data)) + } + + error_definition["stat"] = { + "description": "Uncorrelated statistical uncertainties", + "treatment": "ADD", + "type": "UNCORR", + } + + error_definition["sys_beam"] = { + "description": "Systematic Beam uncertainties", + "treatment": "ADD", + "type": "LHCBBEAM7TEV", + } + + error_definition["sys_luminosity"] = { + "description": "Systematic Luminosity uncertainties", + "treatment": "ADD", + "type": "LHCBLUMI7TEV", + } + + with open("data.yaml", "w") as file: + yaml.dump({"data_central": data}, file, sort_keys=False) + + with open("kinematics.yaml", "w") as file: + yaml.dump({"bins": kinematics}, file, sort_keys=False) + + with open("uncertainties.yaml", "w") as file: + yaml.dump( + {"definitions": error_definition, "bins": errors}, + file, + sort_keys=False, + ) + + +def main_filter() -> None: + """Main driver of the LHCB_ZEE_8TEV filter that produces commmondata. + + There are three main different sources of uncertainties. + + 1. Statistical uncertainties: ADD, UNCORR + + 2. Correlated Systematic uncertainties: ADD, CORR: + constructed by symmetrizing the correlation matrix and extracting + the artificial systematic uncertainties from the latter + + 3. Beam Systematic uncertainties: ADD, LHCBBEAM7TEV + + 4. Luminosity Systematic uncertainties: ADD, LHCBLUMI7TEV + + """ + version, _, _ = read_metadata() + nbpoints = sum(NB_POINTS) + + nbp_idx = 0 + comb_kins, comb_data = [], [] + combined_errors = [] + for tabid in TABLES.keys(): # Loop over tables + for idx in TABLES[tabid]: # Loop over Bosons [Z, W+, W-] + bin_index = [i for i in range(NB_POINTS[nbp_idx])] + yaml_content = load_yaml(table_id=tabid, version=version) + + # Extract the kinematic, data, and uncertainties + kinematics = get_kinematics(yaml_content, bin_index, MAP_TABLE[tabid]) + data_central = get_data_values(yaml_content, bin_index, indx=idx) + uncertainties = get_errors(yaml_content, bin_index, indx=idx) + + # Collect all the results from different tables + comb_kins += kinematics + comb_data += data_central + combined_errors.append(uncertainties) + + nbp_idx += 1 + + errors_combined = concatenate_dicts(combined_errors) + # Compute the Artifical Systematics from CovMat + corrmat = read_corrmatrix(nb_datapoints=nbpoints) + covmat = multiply_syst(corrmat, errors_combined["sys_corr"]) + artunc = generate_artificial_unc( + ndata=nbpoints, + covmat_list=covmat.tolist(), + no_of_norm_mat=0, + ) + errors = format_uncertainties(errors_combined, artunc) + + # Generate all the necessary files + dump_commondata(comb_kins, comb_data, errors) + + return + + +if __name__ == "__main__": + main_filter() diff --git a/buildmaster/LHCB_DY_8TEV_MUON/kinematics.yaml b/buildmaster/LHCB_DY_8TEV_MUON/kinematics.yaml new file mode 100644 index 0000000000..462ad559b3 --- /dev/null +++ b/buildmaster/LHCB_DY_8TEV_MUON/kinematics.yaml @@ -0,0 +1,409 @@ +bins: +- y: + min: 2.0 + mid: 2.0625 + max: 2.125 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.125 + mid: 2.1875 + max: 2.25 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.25 + mid: 2.3125 + max: 2.375 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.375 + mid: 2.4375 + max: 2.5 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.5 + mid: 2.5625 + max: 2.625 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.625 + mid: 2.6875 + max: 2.75 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.75 + mid: 2.8125 + max: 2.875 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.875 + mid: 2.9375 + max: 3.0 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.0 + mid: 3.0625 + max: 3.125 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.125 + mid: 3.1875 + max: 3.25 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.25 + mid: 3.3125 + max: 3.375 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.375 + mid: 3.4375 + max: 3.5 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.5 + mid: 3.5625 + max: 3.625 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.625 + mid: 3.6875 + max: 3.75 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.75 + mid: 3.8125 + max: 3.875 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.875 + mid: 3.9375 + max: 4.0 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 4.0 + mid: 4.125 + max: 4.25 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 4.25 + mid: 4.375 + max: 4.5 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.0 + mid: 2.125 + max: 2.25 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.25 + mid: 2.375 + max: 2.5 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.5 + mid: 2.625 + max: 2.75 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.75 + mid: 2.875 + max: 3.0 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.0 + mid: 3.125 + max: 3.25 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.25 + mid: 3.375 + max: 3.5 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.5 + mid: 3.75 + max: 4.0 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 4.0 + mid: 4.25 + max: 4.5 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.0 + mid: 2.125 + max: 2.25 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.25 + mid: 2.375 + max: 2.5 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.5 + mid: 2.625 + max: 2.75 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.75 + mid: 2.875 + max: 3.0 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.0 + mid: 3.125 + max: 3.25 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.25 + mid: 3.375 + max: 3.5 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.5 + mid: 3.75 + max: 4.0 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 4.0 + mid: 4.25 + max: 4.5 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null diff --git a/buildmaster/LHCB_DY_8TEV_MUON/metadata.yaml b/buildmaster/LHCB_DY_8TEV_MUON/metadata.yaml new file mode 100644 index 0000000000..4d4c073767 --- /dev/null +++ b/buildmaster/LHCB_DY_8TEV_MUON/metadata.yaml @@ -0,0 +1,49 @@ +setname: "LHCB_DY_8TEV_MUON" + +version: 1 +version_comment: "First implementation" + +# References +arXiv: + url: "https://arxiv.org/abs/1511.08039" +iNSPIRE: + url: "https://inspirehep.net/literature/1406555" +hepdata: + url: "https://www.hepdata.net/record/ins1406555" + version: 1 + +nnpdf_metadata: + nnpdf31_process: "DY" + experiment: "LHCb" + +implemented_observables: + +- observable_name: "y" + observable: + description: "Inclusive cross-section for DY production (Z, W+/-) in bins of muon pseudorapidity at 8 TeV" + label: r"$d\sigma / d|y|$" + units: "[fb]" + process_type: "EWK_RAP" + tables: [1, 2] + npoints: [16, 18] + plotting: + kinematics_override: ewk_rap_sqrt_scale + dataset_label: "LHCb $W,Z \\to \\mu$ 8 TeV" + plot_x: y + y_label: '$d\sigma/dy$ (fb)' + kinematic_coverage: [y, M2, sqrt_s] + kinematics: + variables: + y: {description: "Muon rapidity", label: "$y$", units: ""} + M2: {description: "W/Z boson Mass", label: "$M^2$", units: "$GeV^2$"} + sqrt_s: {description: "Center of Mass Energy", label: '$\sqrt{s}$', units: "$GeV$"} + file: kinematics.yaml + data_central: data.yaml + data_uncertainties: + - uncertainties.yaml + theory: + FK_tables: + - - LHCB_DY_8TEV + - LHCB_WP_8TEV + - LHCB_WM_8TEV + operation: 'null' diff --git a/buildmaster/LHCB_DY_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_1.yaml b/buildmaster/LHCB_DY_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_1.yaml new file mode 100644 index 0000000000..0332da55b0 --- /dev/null +++ b/buildmaster/LHCB_DY_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_1.yaml @@ -0,0 +1,184 @@ +dependent_variables: +- header: {name: '$\sigma_{W^{+} \to \mu^{+}\nu}$', units: PB} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: P P --> W+ < MU+ NU > X} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {label: stat, symerror: 1.2} + - {label: sys, symerror: 3.2} + - {label: 'sys,beam', symerror: 2.4} + - {label: 'sys,lumi', symerror: 2.7} + value: 236.5 + - errors: + - {label: stat, symerror: 0.9} + - {label: sys, symerror: 2.2} + - {label: 'sys,beam', symerror: 2.1} + - {label: 'sys,lumi', symerror: 2.4} + value: 208.4 + - errors: + - {label: stat, symerror: 0.8} + - {label: sys, symerror: 1.8} + - {label: 'sys,beam', symerror: 1.8} + - {label: 'sys,lumi', symerror: 2.1} + value: 182.0 + - errors: + - {label: stat, symerror: 0.7} + - {label: sys, symerror: 1.6} + - {label: 'sys,beam', symerror: 1.5} + - {label: 'sys,lumi', symerror: 1.8} + value: 153.3 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.3} + - {label: 'sys,beam', symerror: 1.2} + - {label: 'sys,lumi', symerror: 1.4} + value: 119.5 + - errors: + - {label: stat, symerror: 0.5} + - {label: sys, symerror: 1.0} + - {label: 'sys,beam', symerror: 0.8} + - {label: 'sys,lumi', symerror: 1.0} + value: 84.4 + - errors: + - {label: stat, symerror: 0.5} + - {label: sys, symerror: 1.2} + - {label: 'sys,beam', symerror: 0.9} + - {label: 'sys,lumi', symerror: 1.0} + value: 86.4 + - errors: + - {label: stat, symerror: 0.4} + - {label: sys, symerror: 0.7} + - {label: 'sys,beam', symerror: 0.2} + - {label: 'sys,lumi', symerror: 0.3} + value: 23.0 +- header: {name: '$f^{W^+}_{FSR}$'} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: ''} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {symerror: 0.0047} + value: 1.0188 + - errors: + - {symerror: 0.0028} + value: 1.0163 + - errors: + - {symerror: 0.0025} + value: 1.0158 + - errors: + - {symerror: 0.0028} + value: 1.0148 + - errors: + - {symerror: 0.0032} + value: 1.0152 + - errors: + - {symerror: 0.0046} + value: 1.015 + - errors: + - {symerror: 0.0045} + value: 1.0175 + - errors: + - {symerror: 0.0087} + value: 1.0211 +- header: {name: '$\sigma_{W^{-} \to \mu^{-}\overline{\nu}}$', units: PB} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: P P --> W- < MU- NUBAR > X} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {label: stat, symerror: 0.9} + - {label: sys, symerror: 1.8} + - {label: 'sys,beam', symerror: 1.2} + - {label: 'sys,lumi', symerror: 1.6} + value: 134.0 + - errors: + - {label: stat, symerror: 0.7} + - {label: sys, symerror: 1.4} + - {label: 'sys,beam', symerror: 1.0} + - {label: 'sys,lumi', symerror: 1.4} + value: 119.8 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.2} + - {label: 'sys,beam', symerror: 1.0} + - {label: 'sys,lumi', symerror: 1.3} + value: 110.6 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.2} + - {label: 'sys,beam', symerror: 0.9} + - {label: 'sys,lumi', symerror: 1.2} + value: 102.4 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.1} + - {label: 'sys,beam', symerror: 0.8} + - {label: 'sys,lumi', symerror: 1.1} + value: 92.5 + - errors: + - {label: stat, symerror: 0.5} + - {label: sys, symerror: 0.9} + - {label: 'sys,beam', symerror: 0.7} + - {label: 'sys,lumi', symerror: 0.9} + value: 79.9 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.5} + - {label: 'sys,beam', symerror: 1.0} + - {label: 'sys,lumi', symerror: 1.4} + value: 119.3 + - errors: + - {label: stat, symerror: 0.7} + - {label: sys, symerror: 1.6} + - {label: 'sys,beam', symerror: 0.5} + - {label: 'sys,lumi', symerror: 0.7} + value: 60.0 +- header: {name: '$f^{W^-}_{FSR}$'} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: ''} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {symerror: 0.0026} + value: 1.0172 + - errors: + - {symerror: 0.0027} + value: 1.0155 + - errors: + - {symerror: 0.0028} + value: 1.0153 + - errors: + - {symerror: 0.003} + value: 1.0162 + - errors: + - {symerror: 0.0031} + value: 1.016 + - errors: + - {symerror: 0.0033} + value: 1.0176 + - errors: + - {symerror: 0.0033} + value: 1.02 + - errors: + - {symerror: 0.0053} + value: 1.0243 +independent_variables: +- header: {name: '$\eta^{\mu}$'} + values: + - {high: 2.25, low: 2.0} + - {high: 2.5, low: 2.25} + - {high: 2.75, low: 2.5} + - {high: 3.0, low: 2.75} + - {high: 3.25, low: 3.0} + - {high: 3.5, low: 3.25} + - {high: 4.0, low: 3.5} + - {high: 4.5, low: 4.0} diff --git a/buildmaster/LHCB_DY_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_2.yaml b/buildmaster/LHCB_DY_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_2.yaml new file mode 100644 index 0000000000..9a39670c3b --- /dev/null +++ b/buildmaster/LHCB_DY_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_2.yaml @@ -0,0 +1,200 @@ +dependent_variables: +- header: {name: '$\sigma_{Z \to \mu^{+}\mu^{-}}$', units: PB} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: M(MU+ MU-), units: GeV, value: 60.0-120.0} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: P P --> Z0 < MU+ MU- > X} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {label: stat, symerror: 0.033} + - {label: sys, symerror: 0.055} + - {label: 'sys,beam', symerror: 0.014} + - {label: 'sys,lumi', symerror: 0.014} + value: 1.223 + - errors: + - {label: stat, symerror: 0.051} + - {label: sys, symerror: 0.06} + - {label: 'sys,beam', symerror: 0.038} + - {label: 'sys,lumi', symerror: 0.038} + value: 3.263 + - errors: + - {label: stat, symerror: 0.062} + - {label: sys, symerror: 0.064} + - {label: 'sys,beam', symerror: 0.057} + - {label: 'sys,lumi', symerror: 0.058} + value: 4.983 + - errors: + - {label: stat, symerror: 0.07} + - {label: sys, symerror: 0.072} + - {label: 'sys,beam', symerror: 0.077} + - {label: 'sys,lumi', symerror: 0.078} + value: 6.719 + - errors: + - {label: stat, symerror: 0.076} + - {label: sys, symerror: 0.074} + - {label: 'sys,beam', symerror: 0.093} + - {label: 'sys,lumi', symerror: 0.094} + value: 8.051 + - errors: + - {label: stat, symerror: 0.079} + - {label: sys, symerror: 0.074} + - {label: 'sys,beam', symerror: 0.103} + - {label: 'sys,lumi', symerror: 0.105} + value: 8.967 + - errors: + - {label: stat, symerror: 0.081} + - {label: sys, symerror: 0.076} + - {label: 'sys,beam', symerror: 0.11} + - {label: 'sys,lumi', symerror: 0.112} + value: 9.561 + - errors: + - {label: stat, symerror: 0.082} + - {label: sys, symerror: 0.071} + - {label: 'sys,beam', symerror: 0.113} + - {label: 'sys,lumi', symerror: 0.115} + value: 9.822 + - errors: + - {label: stat, symerror: 0.081} + - {label: sys, symerror: 0.074} + - {label: 'sys,beam', symerror: 0.112} + - {label: 'sys,lumi', symerror: 0.114} + value: 9.721 + - errors: + - {label: stat, symerror: 0.078} + - {label: sys, symerror: 0.071} + - {label: 'sys,beam', symerror: 0.104} + - {label: 'sys,lumi', symerror: 0.105} + value: 9.03 + - errors: + - {label: stat, symerror: 0.072} + - {label: sys, symerror: 0.074} + - {label: 'sys,beam', symerror: 0.089} + - {label: 'sys,lumi', symerror: 0.09} + value: 7.748 + - errors: + - {label: stat, symerror: 0.063} + - {label: sys, symerror: 0.051} + - {label: 'sys,beam', symerror: 0.07} + - {label: 'sys,lumi', symerror: 0.071} + value: 6.059 + - errors: + - {label: stat, symerror: 0.054} + - {label: sys, symerror: 0.041} + - {label: 'sys,beam', symerror: 0.05} + - {label: 'sys,lumi', symerror: 0.051} + value: 4.385 + - errors: + - {label: stat, symerror: 0.042} + - {label: sys, symerror: 0.027} + - {label: 'sys,beam', symerror: 0.031} + - {label: 'sys,lumi', symerror: 0.032} + value: 2.724 + - errors: + - {label: stat, symerror: 0.032} + - {label: sys, symerror: 0.02} + - {label: 'sys,beam', symerror: 0.018} + - {label: 'sys,lumi', symerror: 0.019} + value: 1.584 + - errors: + - {label: stat, symerror: 0.022} + - {label: sys, symerror: 0.012} + - {label: 'sys,beam', symerror: 0.009} + - {label: 'sys,lumi', symerror: 0.009} + value: 0.749 + - errors: + - {label: stat, symerror: 0.016} + - {label: sys, symerror: 0.008} + - {label: 'sys,beam', symerror: 0.004} + - {label: 'sys,lumi', symerror: 0.004} + value: 0.383 + - errors: + - {label: stat, symerror: 0.003} + - {label: sys, symerror: 0.001} + - {label: 'sys,beam', symerror: 0.0} + - {label: 'sys,lumi', symerror: 0.0} + value: 0.011 +- header: {name: '$f^{Z}_{FSR}$'} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: M(MU+ MU-), units: GeV, value: 60.0-120.0} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: P P --> Z0 < MU+ MU- > X} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {symerror: 0.0395} + value: 1.0466 + - errors: + - {symerror: 0.0119} + value: 1.0305 + - errors: + - {symerror: 0.0069} + value: 1.0277 + - errors: + - {symerror: 0.0061} + value: 1.0252 + - errors: + - {symerror: 0.0048} + value: 1.0264 + - errors: + - {symerror: 0.0032} + value: 1.0257 + - errors: + - {symerror: 0.0038} + value: 1.0258 + - errors: + - {symerror: 0.0027} + value: 1.0252 + - errors: + - {symerror: 0.0035} + value: 1.0282 + - errors: + - {symerror: 0.003} + value: 1.0264 + - errors: + - {symerror: 0.0066} + value: 1.0261 + - errors: + - {symerror: 0.004} + value: 1.0248 + - errors: + - {symerror: 0.006} + value: 1.0258 + - errors: + - {symerror: 0.0053} + value: 1.0228 + - errors: + - {symerror: 0.0079} + value: 1.018 + - errors: + - {symerror: 0.01} + value: 1.0207 + - errors: + - {symerror: 0.014} + value: 1.0183 + - errors: + - {symerror: 0.0761} + value: 1.0177 +independent_variables: +- header: {name: '$y_{Z}$'} + values: + - {high: 2.125, low: 2.0} + - {high: 2.25, low: 2.125} + - {high: 2.375, low: 2.25} + - {high: 2.5, low: 2.375} + - {high: 2.625, low: 2.5} + - {high: 2.75, low: 2.625} + - {high: 2.875, low: 2.75} + - {high: 3.0, low: 2.875} + - {high: 3.125, low: 3.0} + - {high: 3.25, low: 3.125} + - {high: 3.375, low: 3.25} + - {high: 3.5, low: 3.375} + - {high: 3.625, low: 3.5} + - {high: 3.75, low: 3.625} + - {high: 3.875, low: 3.75} + - {high: 4.0, low: 3.875} + - {high: 4.25, low: 4.0} + - {high: 4.5, low: 4.25} diff --git a/buildmaster/LHCB_DY_8TEV_MUON/rawdata/corrmat.corr b/buildmaster/LHCB_DY_8TEV_MUON/rawdata/corrmat.corr new file mode 100644 index 0000000000..30782754d8 --- /dev/null +++ b/buildmaster/LHCB_DY_8TEV_MUON/rawdata/corrmat.corr @@ -0,0 +1,34 @@ +1.0 +0.19 1.0 +0.17 0.27 1.0 +0.16 0.26 0.28 1.0 +0.16 0.25 0.28 0.29 1.0 +0.15 0.24 0.27 0.29 0.3 1.0 +0.14 0.23 0.26 0.28 0.29 0.3 1.0 +0.14 0.21 0.25 0.27 0.29 0.3 0.3 1.0 +0.13 0.2 0.23 0.25 0.27 0.28 0.29 0.29 1.0 +0.11 0.17 0.2 0.23 0.25 0.26 0.27 0.28 0.27 1.0 +0.09 0.14 0.16 0.18 0.2 0.22 0.22 0.23 0.23 0.23 1.0 +0.08 0.12 0.15 0.17 0.19 0.2 0.21 0.22 0.22 0.22 0.2 1.0 +0.07 0.1 0.12 0.14 0.16 0.17 0.18 0.19 0.19 0.2 0.18 0.19 1.0 +0.06 0.08 0.1 0.11 0.13 0.14 0.15 0.16 0.16 0.17 0.16 0.16 0.15 1.0 +0.05 0.07 0.08 0.09 0.1 0.11 0.11 0.12 0.13 0.13 0.12 0.13 0.12 0.11 1.0 +0.03 0.05 0.06 0.06 0.07 0.08 0.08 0.09 0.09 0.1 0.09 0.1 0.09 0.08 0.07 1.0 +0.03 0.04 0.04 0.05 0.05 0.06 0.06 0.06 0.07 0.07 0.07 0.08 0.07 0.07 0.06 0.05 1.0 +0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.02 0.02 0.01 0.01 0.01 0.01 1.0 +0.23 0.3 0.28 0.27 0.26 0.25 0.24 0.23 0.21 0.18 0.15 0.13 0.11 0.1 0.07 0.05 0.04 0.01 1.0 +0.05 0.15 0.21 0.2 0.2 0.2 0.2 0.19 0.18 0.16 0.14 0.12 0.1 0.08 0.06 0.04 0.03 0.01 0.2 1.0 +0.04 0.07 0.12 0.15 0.16 0.17 0.17 0.17 0.16 0.15 0.13 0.13 0.11 0.08 0.06 0.05 0.03 0.01 0.13 0.12 1.0 +0.05 0.08 0.1 0.13 0.16 0.17 0.17 0.17 0.16 0.16 0.14 0.13 0.12 0.09 0.07 0.05 0.03 0.01 0.06 0.03 0.26 1.0 +0.06 0.08 0.1 0.11 0.14 0.16 0.17 0.17 0.16 0.16 0.14 0.14 0.12 0.1 0.07 0.05 0.03 0.01 0.07 0.03 0.25 0.33 1.0 +0.04 0.06 0.07 0.09 0.1 0.11 0.12 0.13 0.13 0.12 0.11 0.11 0.09 0.08 0.06 0.04 0.03 0.0 0.03 -0.01 0.28 0.35 0.34 1.0 +0.04 0.06 0.07 0.08 0.09 0.1 0.11 0.12 0.12 0.12 0.11 0.11 0.1 0.09 0.07 0.05 0.04 0.0 0.0 -0.06 0.31 0.41 0.4 0.45 1.0 +0.02 0.03 0.04 0.04 0.04 0.04 0.05 0.05 0.06 0.06 0.06 0.07 0.06 0.06 0.05 0.04 0.04 0.01 -0.07 -0.14 0.14 0.32 0.32 0.35 0.40 1.0 +0.21 0.28 0.26 0.25 0.24 0.24 0.22 0.21 0.2 0.17 0.14 0.12 0.11 0.09 0.07 0.05 0.04 0.01 0.67 0.10 0.24 0.22 0.22 0.23 0.26 0.14 1.0 +0.04 0.14 0.19 0.18 0.18 0.18 0.18 0.17 0.16 0.15 0.12 0.11 0.09 0.07 0.05 0.04 0.03 0.0 0.07 0.54 0.23 0.28 0.28 0.28 0.33 0.24 0.21 1.0 +0.04 0.07 0.11 0.14 0.15 0.15 0.15 0.15 0.15 0.14 0.12 0.12 0.1 0.08 0.06 0.04 0.03 0.01 0.05 0.03 0.64 0.27 0.26 0.27 0.32 0.23 0.18 0.22 1.0 +0.05 0.07 0.09 0.12 0.14 0.15 0.15 0.16 0.15 0.14 0.12 0.12 0.11 0.08 0.06 0.04 0.03 0.01 0.04 0.00 0.25 0.70 0.30 0.33 0.39 0.29 0.21 0.25 0.31 1.0 +0.05 0.08 0.09 0.1 0.13 0.15 0.15 0.16 0.15 0.15 0.13 0.13 0.11 0.09 0.07 0.05 0.03 0.01 0.06 0.03 0.26 0.32 0.68 0.32 0.38 0.26 0.22 0.28 0.27 0.32 1.0 +0.04 0.06 0.08 0.09 0.1 0.12 0.13 0.13 0.13 0.13 0.11 0.11 0.1 0.08 0.06 0.04 0.03 0.0 0.07 0.04 0.30 0.28 0.27 0.63 0.36 0.22 0.23 0.23 0.29 0.32 0.28 1.0 +0.04 0.06 0.08 0.09 0.1 0.11 0.12 0.13 0.14 0.14 0.12 0.12 0.11 0.1 0.08 0.06 0.04 0.01 0.14 0.20 -0.13 -0.08 -0.07 -0.17 -0.04 -0.15 -0.06 -0.04 -0.04 -0.11 -0.07 -0.19 1.0 +0.03 0.04 0.04 0.05 0.05 0.06 0.06 0.06 0.07 0.07 0.07 0.08 0.08 0.07 0.06 0.05 0.04 0.01 0.12 0.17 -0.18 -0.17 -0.15 -0.23 -0.31 0.05 -0.09 -0.11 -0.14 -0.19 -0.18 -0.24 0.48 1.0 diff --git a/buildmaster/LHCB_DY_8TEV_MUON/uncertainties.yaml b/buildmaster/LHCB_DY_8TEV_MUON/uncertainties.yaml new file mode 100644 index 0000000000..afce947435 --- /dev/null +++ b/buildmaster/LHCB_DY_8TEV_MUON/uncertainties.yaml @@ -0,0 +1,1408 @@ +definitions: + sys_corr_1: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_2: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_3: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_4: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_5: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_6: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_7: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_8: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_9: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_10: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_11: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_12: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_13: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_14: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_15: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_16: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_17: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_18: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_19: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_20: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_21: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_22: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_23: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_24: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_25: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_26: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_27: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_28: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_29: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_30: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_31: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_32: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_33: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_34: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + stat: + description: Uncorrelated statistical uncertainties + treatment: ADD + type: UNCORR + sys_beam: + description: Systematic Beam uncertainties + treatment: ADD + type: LHCBBEAM7TEV + sys_luminosity: + description: Systematic Luminosity uncertainties + treatment: ADD + type: LHCBLUMI7TEV +bins: +- sys_corr_1: -13.129371357446566 + sys_corr_2: 0.6853465481125789 + sys_corr_3: -1.8136153124139376 + sys_corr_4: -1.6378919565265173 + sys_corr_5: 1.0056447881503754 + sys_corr_6: 0.9250488794408771 + sys_corr_7: -0.6130610634103019 + sys_corr_8: -1.033777871097276 + sys_corr_9: 0.9022978319043655 + sys_corr_10: 0.5049756150293241 + sys_corr_11: -0.8765726686084679 + sys_corr_12: 0.5378206330642871 + sys_corr_13: -0.10895523093285057 + sys_corr_14: -0.4881555518108563 + sys_corr_15: -0.265568053839214 + sys_corr_16: 0.03401879609850996 + sys_corr_17: -11.515944076984733 + sys_corr_18: -9.94493656819381e-05 + sys_corr_19: -0.007574441556960995 + sys_corr_20: 0.0036078709274008384 + sys_corr_21: 0.054497988233331206 + sys_corr_22: -0.08061168467584123 + sys_corr_23: 0.2232461165909764 + sys_corr_24: 5.565016708385234 + sys_corr_25: 0.2609714050586262 + sys_corr_26: -39.974071914597864 + sys_corr_27: -27.943155261002897 + sys_corr_28: -16.32894214901613 + sys_corr_29: -4.436995222440807 + sys_corr_30: -0.8115865284675327 + sys_corr_31: -0.3755659425689436 + sys_corr_32: 1.9696013567201123 + sys_corr_33: -2.6324056596660883 + sys_corr_34: -0.1664720974747847 + stat: 33.0 + sys_beam: 14.0 + sys_luminosity: 14.0 +- sys_corr_1: -19.84410989915685 + sys_corr_2: -0.33428052767614547 + sys_corr_3: 2.6709227371602084 + sys_corr_4: -1.122291532865575 + sys_corr_5: 2.2445396406954448 + sys_corr_6: 1.2318412040028195 + sys_corr_7: -0.9643207591596343 + sys_corr_8: -2.2908607279681825 + sys_corr_9: 1.134280089906541 + sys_corr_10: -0.8548057471102127 + sys_corr_11: -1.796249557177666 + sys_corr_12: -0.07911083910570488 + sys_corr_13: -0.472951502025977 + sys_corr_14: 0.2558186095238292 + sys_corr_15: -0.3882769649810862 + sys_corr_16: -0.06521893670249582 + sys_corr_17: -20.03065356651147 + sys_corr_18: -5.938906840067905e-05 + sys_corr_19: -0.005317597654044067 + sys_corr_20: -0.012280961378499562 + sys_corr_21: 0.04723338963268214 + sys_corr_22: 0.08725074045417443 + sys_corr_23: -0.2200826518235283 + sys_corr_24: 10.334303066975695 + sys_corr_25: -1.8910305617977035 + sys_corr_26: 31.412187227642292 + sys_corr_27: -31.612351485478605 + sys_corr_28: -24.0763415362685 + sys_corr_29: -7.556695045723271 + sys_corr_30: 1.0636142561345379 + sys_corr_31: 1.678562392078886 + sys_corr_32: 3.831782907879968 + sys_corr_33: -4.904694970016218 + sys_corr_34: -1.4462770772783355 + stat: 51.0 + sys_beam: 38.0 + sys_luminosity: 38.0 +- sys_corr_1: -21.06803631144305 + sys_corr_2: -2.8187263609581716 + sys_corr_3: 7.429751053129648 + sys_corr_4: -1.2363854092311348 + sys_corr_5: 0.440020299276015 + sys_corr_6: 0.8285805636548081 + sys_corr_7: -0.1833307684413444 + sys_corr_8: -1.8866832461202852 + sys_corr_9: 1.0748601041255013 + sys_corr_10: -0.9841457536126791 + sys_corr_11: -2.2265738633784835 + sys_corr_12: -0.5655169182098072 + sys_corr_13: -0.4881568704060051 + sys_corr_14: -0.782035779352363 + sys_corr_15: -0.05919509595153284 + sys_corr_16: -1.1097250976141622 + sys_corr_17: -24.699858197476445 + sys_corr_18: -3.105265932964177e-05 + sys_corr_19: 0.0027633726426610577 + sys_corr_20: -0.01714976905989588 + sys_corr_21: 0.03260600131367636 + sys_corr_22: -0.02846572509875826 + sys_corr_23: -0.1467691264227021 + sys_corr_24: 12.565805182397174 + sys_corr_25: 0.625615759576152 + sys_corr_26: -3.5879461032450957 + sys_corr_27: 31.925361874882043 + sys_corr_28: -40.75747797840909 + sys_corr_29: -8.403699350561684 + sys_corr_30: -2.1525840708834774 + sys_corr_31: 1.1118986638848498 + sys_corr_32: 3.4014199295399585 + sys_corr_33: -4.750919333226475 + sys_corr_34: -1.0948035576743735 + stat: 62.0 + sys_beam: 57.0 + sys_luminosity: 58.0 +- sys_corr_1: -23.338793935211896 + sys_corr_2: -5.015020292076225 + sys_corr_3: 8.195207361706466 + sys_corr_4: -3.6167946329492535 + sys_corr_5: -1.5375896951455132 + sys_corr_6: -0.5055645332316372 + sys_corr_7: -0.1656605444022098 + sys_corr_8: -1.3674236995197913 + sys_corr_9: 0.9863174769320691 + sys_corr_10: -0.11689147629752011 + sys_corr_11: -2.76279796378281 + sys_corr_12: -0.7334770006813159 + sys_corr_13: -0.39710202721088006 + sys_corr_14: -0.9513343420516247 + sys_corr_15: -0.5236614050232143 + sys_corr_16: 0.14328103033362025 + sys_corr_17: -33.108074461181786 + sys_corr_18: -1.1676986084828315e-05 + sys_corr_19: -0.005199428393585041 + sys_corr_20: 0.0030204570342489963 + sys_corr_21: 0.04012076569205404 + sys_corr_22: 0.003803307862097004 + sys_corr_23: 0.13823101139063843 + sys_corr_24: 20.18490057862894 + sys_corr_25: 0.783534475773199 + sys_corr_26: -1.5428612589110309 + sys_corr_27: 2.9113392297877394 + sys_corr_28: 19.877343744892453 + sys_corr_29: -18.817522495155952 + sys_corr_30: 5.171011463669989 + sys_corr_31: 3.266548144634222 + sys_corr_32: 45.33543603498649 + sys_corr_33: -6.762197324713448 + sys_corr_34: -10.777895659545207 + stat: 70.0 + sys_beam: 77.0 + sys_luminosity: 78.0 +- sys_corr_1: -23.62879238801978 + sys_corr_2: -6.871934369335336 + sys_corr_3: 8.838236802397901 + sys_corr_4: -5.75372720918565 + sys_corr_5: -1.5402984336354577 + sys_corr_6: -0.3451192706498088 + sys_corr_7: 0.413105816684792 + sys_corr_8: 0.5073510714926635 + sys_corr_9: 1.5788941118756248 + sys_corr_10: 0.7607246969863836 + sys_corr_11: -2.673007365634146 + sys_corr_12: -1.2417474817490068 + sys_corr_13: -0.7432869137172015 + sys_corr_14: -0.598997451623972 + sys_corr_15: -0.06637383592033852 + sys_corr_16: 0.741162658809314 + sys_corr_17: -36.499510974556024 + sys_corr_18: 5.846611289286119e-06 + sys_corr_19: -0.0001752413952545087 + sys_corr_20: -0.008085268741514851 + sys_corr_21: 0.05044391052283411 + sys_corr_22: -0.15328563391070796 + sys_corr_23: 0.49110134937679384 + sys_corr_24: 17.63914440713869 + sys_corr_25: 1.5555770256950772 + sys_corr_26: -0.29419626869461024 + sys_corr_27: 1.4513100444746851 + sys_corr_28: 12.997331532175536 + sys_corr_29: -22.22192546156861 + sys_corr_30: 7.42411526077336 + sys_corr_31: -10.124571665104387 + sys_corr_32: -22.849337902782768 + sys_corr_33: -10.788104366418192 + sys_corr_34: 40.447688777367155 + stat: 76.0 + sys_beam: 93.0 + sys_luminosity: 94.0 +- sys_corr_1: -23.305638822180896 + sys_corr_2: -7.8937366993396365 + sys_corr_3: 9.253795984805592 + sys_corr_4: -7.241905559104537 + sys_corr_5: -1.9123722247234651 + sys_corr_6: 0.5912138549587218 + sys_corr_7: 0.0465224938402271 + sys_corr_8: 0.7199048005800637 + sys_corr_9: 1.828692488405049 + sys_corr_10: 2.1832257020458705 + sys_corr_11: -2.015036523294676 + sys_corr_12: -1.8353749829104327 + sys_corr_13: 0.0015187347231060773 + sys_corr_14: -0.48209896219515624 + sys_corr_15: 1.213798569671312 + sys_corr_16: 1.1211154805244627 + sys_corr_17: -36.92425033443672 + sys_corr_18: 1.4188601727116922e-05 + sys_corr_19: -0.008850403418020465 + sys_corr_20: -0.02487825307330309 + sys_corr_21: 0.09023808536309276 + sys_corr_22: -0.21273521070752027 + sys_corr_23: 0.568346566723835 + sys_corr_24: 9.560426327356453 + sys_corr_25: 1.6835388239353692 + sys_corr_26: -0.6834505876496619 + sys_corr_27: 0.6825047572018638 + sys_corr_28: 7.890847045554724 + sys_corr_29: -14.45530599394527 + sys_corr_30: 10.42790997030499 + sys_corr_31: -2.6562075756183785 + sys_corr_32: -29.80150469584712 + sys_corr_33: 18.833535008874055 + sys_corr_34: -40.257082690578 + stat: 79.0 + sys_beam: 103.0 + sys_luminosity: 105.0 +- sys_corr_1: -23.122870553655478 + sys_corr_2: -8.472992295351233 + sys_corr_3: 9.978884071326496 + sys_corr_4: -8.191892596977262 + sys_corr_5: -2.170985058688932 + sys_corr_6: 1.1927202283582423 + sys_corr_7: 1.027768527951832 + sys_corr_8: 1.3222563308413176 + sys_corr_9: 0.3390668514570481 + sys_corr_10: 2.5599804325753284 + sys_corr_11: -1.6074832724593935 + sys_corr_12: -1.6149874429471374 + sys_corr_13: -0.18163026461144396 + sys_corr_14: -1.317510032846632 + sys_corr_15: 1.2937203505127906 + sys_corr_16: 0.4372359769181125 + sys_corr_17: -39.002034940744274 + sys_corr_18: 1.5486962831135256e-05 + sys_corr_19: -0.006718922393629441 + sys_corr_20: -0.018187853006976003 + sys_corr_21: 0.059586532298960715 + sys_corr_22: -0.29423833975817293 + sys_corr_23: 0.8058671737702604 + sys_corr_24: 7.740143059284684 + sys_corr_25: 2.274379055527933 + sys_corr_26: -0.94388224336938 + sys_corr_27: 0.23546915242985275 + sys_corr_28: 4.67316228570499 + sys_corr_29: 31.63416661772606 + sys_corr_30: 4.348764566808575 + sys_corr_31: 47.20725640418786 + sys_corr_32: -6.675227803548107 + sys_corr_33: -7.705400923277242 + sys_corr_34: 4.070133302521287 + stat: 81.0 + sys_beam: 110.0 + sys_luminosity: 112.0 +- sys_corr_1: -20.85999640061773 + sys_corr_2: -8.242540680871082 + sys_corr_3: 8.961204056550143 + sys_corr_4: -8.614385898951557 + sys_corr_5: -2.306048360098447 + sys_corr_6: 1.3799189252318764 + sys_corr_7: 1.763104057565564 + sys_corr_8: 1.2044387625997617 + sys_corr_9: -0.029223529687805137 + sys_corr_10: 3.0215802371301006 + sys_corr_11: -1.4142992542156707 + sys_corr_12: -1.1889723630536118 + sys_corr_13: 0.9257611979897812 + sys_corr_14: -0.6372331880602672 + sys_corr_15: 0.3263646033902244 + sys_corr_16: 0.8296699988914757 + sys_corr_17: -34.31148751935087 + sys_corr_18: 2.223371222085629e-05 + sys_corr_19: -0.004079579657142941 + sys_corr_20: -0.04109755671287726 + sys_corr_21: 0.12179926772253029 + sys_corr_22: -0.4318081034916117 + sys_corr_23: 1.2383669027745627 + sys_corr_24: -0.555349911332747 + sys_corr_25: 3.847663593630981 + sys_corr_26: 1.214500716621149 + sys_corr_27: -1.0896588789845654 + sys_corr_28: 5.220822383009355 + sys_corr_29: 4.056607732576169 + sys_corr_30: -53.73393248419216 + sys_corr_31: -5.961938629073214 + sys_corr_32: 0.1723904768898964 + sys_corr_33: 13.699964798406059 + sys_corr_34: 2.433988171680721 + stat: 82.0 + sys_beam: 113.0 + sys_luminosity: 115.0 +- sys_corr_1: -20.223983250604604 + sys_corr_2: -8.154175836436359 + sys_corr_3: 9.26327356637858 + sys_corr_4: -9.718752686919935 + sys_corr_5: -2.815203846765515 + sys_corr_6: 1.8702305366091985 + sys_corr_7: 1.6422084351635182 + sys_corr_8: -0.16224786719614528 + sys_corr_9: -0.22452711746931522 + sys_corr_10: 3.572768532873175 + sys_corr_11: -1.807628621364004 + sys_corr_12: -0.9629409377418384 + sys_corr_13: 0.5171356641657681 + sys_corr_14: -0.6176549011257892 + sys_corr_15: 0.29396865202369565 + sys_corr_16: 0.19228110231051582 + sys_corr_17: -35.97332885136173 + sys_corr_18: 2.1030318953479127e-05 + sys_corr_19: -0.01573421063351599 + sys_corr_20: -0.03959349841824993 + sys_corr_21: 0.20298279711011666 + sys_corr_22: -0.43673348436197756 + sys_corr_23: 1.222881890842337 + sys_corr_24: -8.050128059309468 + sys_corr_25: 3.706780158719833 + sys_corr_26: 0.2856448791722523 + sys_corr_27: 0.3134908927527956 + sys_corr_28: 1.8316820476472904 + sys_corr_29: 32.11772042509791 + sys_corr_30: 6.0053293548508035 + sys_corr_31: -35.91989545508118 + sys_corr_32: -0.02011206044083529 + sys_corr_33: -30.62631510390598 + sys_corr_34: -10.705014255279856 + stat: 81.0 + sys_beam: 112.0 + sys_luminosity: 114.0 +- sys_corr_1: -17.04123512385607 + sys_corr_2: -7.918393116785677 + sys_corr_3: 8.440143768234751 + sys_corr_4: -10.024020107965285 + sys_corr_5: -2.6847034507320844 + sys_corr_6: 2.0671669249445555 + sys_corr_7: 1.837279687306533 + sys_corr_8: 0.3650568659958519 + sys_corr_9: -0.4904354192346455 + sys_corr_10: 3.0568684568954443 + sys_corr_11: -1.0638801573487726 + sys_corr_12: -1.3582456192662984 + sys_corr_13: 0.07169286786156893 + sys_corr_14: -0.2199043557753311 + sys_corr_15: 1.5744304454334228 + sys_corr_16: -0.0013292406672338916 + sys_corr_17: -32.11721511436332 + sys_corr_18: 1.563925464777089e-05 + sys_corr_19: -0.01602557329814359 + sys_corr_20: -0.06980990577079062 + sys_corr_21: 0.2132859532473274 + sys_corr_22: -0.6435224332446 + sys_corr_23: 1.9689443448390334 + sys_corr_24: -12.899652788631762 + sys_corr_25: 5.004417879944764 + sys_corr_26: 0.4905581183348719 + sys_corr_27: -0.6767901157597395 + sys_corr_28: -5.192725175071233 + sys_corr_29: 15.90697875227382 + sys_corr_30: 20.16293066600603 + sys_corr_31: -11.987331589786724 + sys_corr_32: 14.255145302365655 + sys_corr_33: 43.78533937999355 + sys_corr_34: 17.535461370942752 + stat: 78.0 + sys_beam: 104.0 + sys_luminosity: 105.0 +- sys_corr_1: -14.949658691661552 + sys_corr_2: -7.165496747473827 + sys_corr_3: 7.808477726717307 + sys_corr_4: -9.56305409326678 + sys_corr_5: -2.670729336856583 + sys_corr_6: 1.870910013082464 + sys_corr_7: 1.1343784857727517 + sys_corr_8: 0.5731050576349122 + sys_corr_9: -0.8237338552005756 + sys_corr_10: 3.4751925833597936 + sys_corr_11: -0.13851422837746447 + sys_corr_12: -0.10023616022750506 + sys_corr_13: -0.297566757573512 + sys_corr_14: -0.01650093308831284 + sys_corr_15: 0.6964491524517722 + sys_corr_16: 0.020654173191060936 + sys_corr_17: -29.914143247343116 + sys_corr_18: -2.595663656276923e-06 + sys_corr_19: -0.020768693638988766 + sys_corr_20: -0.05834119476442283 + sys_corr_21: 0.19884628619594974 + sys_corr_22: -0.6436033819979039 + sys_corr_23: 1.6165938253859384 + sys_corr_24: -56.831154776699634 + sys_corr_25: 4.207755805618636 + sys_corr_26: 0.13546165938685006 + sys_corr_27: 0.16794567880689754 + sys_corr_28: -1.254823525202448 + sys_corr_29: -23.915653026107577 + sys_corr_30: 1.316650496021619 + sys_corr_31: 12.59376213706475 + sys_corr_32: 1.475379610059091 + sys_corr_33: -11.882612801600084 + sys_corr_34: -0.6558219368176063 + stat: 72.0 + sys_beam: 89.0 + sys_luminosity: 90.0 +- sys_corr_1: -9.14904976731458 + sys_corr_2: -4.972284733299394 + sys_corr_3: 4.934501253783746 + sys_corr_4: -7.146715955712768 + sys_corr_5: -2.4806832426488516 + sys_corr_6: 1.6716036517481274 + sys_corr_7: 0.3731925008926304 + sys_corr_8: 0.46906925363558766 + sys_corr_9: -1.0870814502316777 + sys_corr_10: 2.1864260496234604 + sys_corr_11: -0.08527632330590701 + sys_corr_12: -0.012526363709937075 + sys_corr_13: 0.424383407931963 + sys_corr_14: -0.13456579621486817 + sys_corr_15: 0.479004540879783 + sys_corr_16: -0.15340419169042 + sys_corr_17: -14.879810632978513 + sys_corr_18: -0.00024564680501881395 + sys_corr_19: -0.0454435934292143 + sys_corr_20: -0.12517695573643947 + sys_corr_21: 0.4495879379661885 + sys_corr_22: -1.2047889254655135 + sys_corr_23: 6.991666808010768 + sys_corr_24: -5.151973939331355 + sys_corr_25: -45.641593280488614 + sys_corr_26: -1.3784434779206556 + sys_corr_27: 1.565294265899055 + sys_corr_28: 1.4093209031809397 + sys_corr_29: 2.3280091228437603 + sys_corr_30: -0.7987670255568619 + sys_corr_31: -1.492549105845241 + sys_corr_32: 0.1634737023690217 + sys_corr_33: 2.02289736511137 + sys_corr_34: 1.008066369978752 + stat: 63.0 + sys_beam: 70.0 + sys_luminosity: 71.0 +- sys_corr_1: -6.326085258885516 + sys_corr_2: -3.439415506883048 + sys_corr_3: 3.3413277376739843 + sys_corr_4: -5.541898966162248 + sys_corr_5: -1.666176480784771 + sys_corr_6: 0.979375217874185 + sys_corr_7: 0.0002485455875154854 + sys_corr_8: -0.08187806243705073 + sys_corr_9: -0.43331484389360025 + sys_corr_10: 1.7688356519672581 + sys_corr_11: 0.16738920479722036 + sys_corr_12: 0.11055325680962605 + sys_corr_13: 0.4269058132515623 + sys_corr_14: -0.1138800538385906 + sys_corr_15: 0.9613977606946876 + sys_corr_16: -0.013820528956704024 + sys_corr_17: -9.506358516911547 + sys_corr_18: -0.000322491958865001 + sys_corr_19: -0.04499426072423848 + sys_corr_20: -0.13755544896289332 + sys_corr_21: 0.5827172665108897 + sys_corr_22: -1.9470722746536389 + sys_corr_23: -37.513389438295306 + sys_corr_24: -3.666147114225558 + sys_corr_25: -7.529167532066956 + sys_corr_26: -0.6017477237226491 + sys_corr_27: 0.16993176961810133 + sys_corr_28: 0.8185166967637814 + sys_corr_29: 1.5139539639474104 + sys_corr_30: -0.23812496924184945 + sys_corr_31: -0.8615737285431316 + sys_corr_32: 0.09807293782358924 + sys_corr_33: 1.5488833831456272 + sys_corr_34: 0.7674261133314435 + stat: 54.0 + sys_beam: 50.0 + sys_luminosity: 51.0 +- sys_corr_1: -3.581622334817089 + sys_corr_2: -1.5952014747683092 + sys_corr_3: 1.696419758223885 + sys_corr_4: -3.146991276459658 + sys_corr_5: -0.8213111707130737 + sys_corr_6: 0.8751706472028377 + sys_corr_7: 0.2206946879339562 + sys_corr_8: -0.10396281407173584 + sys_corr_9: -0.6365967085276881 + sys_corr_10: 0.9100948097624815 + sys_corr_11: 0.13025444043125295 + sys_corr_12: 0.34852058051345974 + sys_corr_13: 0.000963419400960771 + sys_corr_14: -0.024415407833905774 + sys_corr_15: 0.23683126623679324 + sys_corr_16: -0.213123518257198 + sys_corr_17: -4.887880950067507 + sys_corr_18: -9.666121427649194e-05 + sys_corr_19: -0.08679639349439726 + sys_corr_20: -0.21795914219962417 + sys_corr_21: 1.5906032250675446 + sys_corr_22: 25.605815007759354 + sys_corr_23: -2.2802107816509967 + sys_corr_24: -2.064997634263775 + sys_corr_25: -2.1904506948921627 + sys_corr_26: -0.31575888365863575 + sys_corr_27: 0.13593361557791478 + sys_corr_28: 0.2645398552247651 + sys_corr_29: 0.7254250538696598 + sys_corr_30: -0.14873611967742692 + sys_corr_31: -0.34721727099001154 + sys_corr_32: -0.06366415037975065 + sys_corr_33: 0.7056752901669989 + sys_corr_34: 0.34189479965746694 + stat: 42.0 + sys_beam: 31.0 + sys_luminosity: 32.0 +- sys_corr_1: -1.92704439512093 + sys_corr_2: -0.8906975447977756 + sys_corr_3: 0.9859044565822679 + sys_corr_4: -1.900006590186975 + sys_corr_5: -0.500061597538644 + sys_corr_6: 0.46968476717803787 + sys_corr_7: 0.026567509007107172 + sys_corr_8: -0.24644151318871022 + sys_corr_9: -0.25414422874330955 + sys_corr_10: 0.5977450748892065 + sys_corr_11: 0.12461147105302377 + sys_corr_12: 0.22106489547048022 + sys_corr_13: 0.002600996773599335 + sys_corr_14: 0.18567315840954113 + sys_corr_15: 0.11880756529326222 + sys_corr_16: -0.18567999103818122 + sys_corr_17: -2.8272588693821663 + sys_corr_18: -0.00019905744405662599 + sys_corr_19: -0.11219010309646388 + sys_corr_20: -0.36203400976763017 + sys_corr_21: -19.33064816976645 + sys_corr_22: 1.9792982895260696 + sys_corr_23: -1.076070580397553 + sys_corr_24: -0.9794399369760031 + sys_corr_25: -1.2690397198581644 + sys_corr_26: -0.11077568336747423 + sys_corr_27: -0.0455017730958743 + sys_corr_28: 0.014709472995496412 + sys_corr_29: 0.3343664195538066 + sys_corr_30: 0.011541248946823261 + sys_corr_31: -0.3839086612900922 + sys_corr_32: 0.06978441016448338 + sys_corr_33: 0.26523110211595075 + sys_corr_34: 0.06515123861157356 + stat: 32.0 + sys_beam: 18.0 + sys_luminosity: 19.0 +- sys_corr_1: -0.8296601388639502 + sys_corr_2: -0.3805251722782747 + sys_corr_3: 0.4256454984305268 + sys_corr_4: -0.8640659174103269 + sys_corr_5: -0.3029421108714545 + sys_corr_6: 0.2176950458024982 + sys_corr_7: -0.0513015110495244 + sys_corr_8: -0.12919103018667238 + sys_corr_9: -0.10721047160309029 + sys_corr_10: 0.1529847335090339 + sys_corr_11: 0.17487604901206083 + sys_corr_12: -0.01248019049416586 + sys_corr_13: 0.04226511417819689 + sys_corr_14: 0.07817088485442444 + sys_corr_15: 0.0736083183745245 + sys_corr_16: -0.10963219724155511 + sys_corr_17: -1.2120042253538708 + sys_corr_18: -0.00046489040093820644 + sys_corr_19: -0.25803633958738176 + sys_corr_20: 11.795834334234783 + sys_corr_21: -0.5401002869203217 + sys_corr_22: 0.47787211458447243 + sys_corr_23: -0.4012872851787313 + sys_corr_24: -0.4773980799989664 + sys_corr_25: -0.5578977015713101 + sys_corr_26: 0.015910316662032797 + sys_corr_27: 0.039315340097579664 + sys_corr_28: -0.034916993021117164 + sys_corr_29: 0.15079334377683432 + sys_corr_30: -0.025448127463030965 + sys_corr_31: -0.12947869114049182 + sys_corr_32: 0.003045495350807366 + sys_corr_33: 0.21539792088791174 + sys_corr_34: 0.050432772372592215 + stat: 22.0 + sys_beam: 9.0 + sys_luminosity: 9.0 +- sys_corr_1: -0.42393728550091253 + sys_corr_2: -0.14894396859548853 + sys_corr_3: 0.1981748571052791 + sys_corr_4: -0.39930231793851473 + sys_corr_5: -0.11765515578109007 + sys_corr_6: 0.10694488123136513 + sys_corr_7: -0.0737090810991067 + sys_corr_8: -0.13441755182880608 + sys_corr_9: -0.11536049109165246 + sys_corr_10: 0.04875973725724547 + sys_corr_11: 0.04521689153192965 + sys_corr_12: 0.10130389568753956 + sys_corr_13: 0.04322115716786899 + sys_corr_14: 0.028136185353375684 + sys_corr_15: 0.004731362627832487 + sys_corr_16: -0.11516055102024063 + sys_corr_17: -0.6185243817055507 + sys_corr_18: -0.0008188847906671834 + sys_corr_19: 7.911038971966661 + sys_corr_20: 0.37528005595773883 + sys_corr_21: -0.26672538386835976 + sys_corr_22: 0.3020038060076616 + sys_corr_23: -0.21369897814298744 + sys_corr_24: -0.2553303274366322 + sys_corr_25: -0.33152188617632095 + sys_corr_26: -0.03180232046108036 + sys_corr_27: -0.045026544656933135 + sys_corr_28: 0.015861015519060512 + sys_corr_29: 0.06645220360768585 + sys_corr_30: 0.039901597979006935 + sys_corr_31: -0.053250101359537255 + sys_corr_32: 0.028197119479147756 + sys_corr_33: 0.04831256818640388 + sys_corr_34: -0.01828330836598994 + stat: 16.0 + sys_beam: 4.0 + sys_luminosity: 4.0 +- sys_corr_1: -0.013414184810947182 + sys_corr_2: -0.0037960802643431835 + sys_corr_3: 0.007050325110555244 + sys_corr_4: -0.012793678842011953 + sys_corr_5: -0.006988237873416711 + sys_corr_6: -0.0011527303922321177 + sys_corr_7: -0.0038282744129991384 + sys_corr_8: 0.00580351666510989 + sys_corr_9: 0.010110892923237483 + sys_corr_10: 0.006099701507322768 + sys_corr_11: -0.002713842573751903 + sys_corr_12: 0.0038844748413824986 + sys_corr_13: -0.00021929479123481416 + sys_corr_14: -5.673249487341189e-05 + sys_corr_15: -0.005283439658346962 + sys_corr_16: -0.010413435144737736 + sys_corr_17: -0.009022711855041334 + sys_corr_18: 0.9993144938282417 + sys_corr_19: 0.006305972548180469 + sys_corr_20: 0.005629159396434368 + sys_corr_21: -0.003870372244719926 + sys_corr_22: 0.0024480503883864435 + sys_corr_23: -0.011281196460354895 + sys_corr_24: -0.0016038955252630446 + sys_corr_25: -0.015001223022070295 + sys_corr_26: -0.0028588158298082783 + sys_corr_27: -0.0032024298307163834 + sys_corr_28: -0.00378958229068439 + sys_corr_29: -0.0012810170230193785 + sys_corr_30: 0.00020952138697761106 + sys_corr_31: -0.0002487670121694488 + sys_corr_32: 0.0015993114173559364 + sys_corr_33: 1.83609916547212e-05 + sys_corr_34: 0.00045379693968826006 + stat: 3.0 + sys_beam: 0.0 + sys_luminosity: 0.0 +- sys_corr_1: -3038.5793019795487 + sys_corr_2: 768.151619729029 + sys_corr_3: -485.10950751290466 + sys_corr_4: 53.08972426745508 + sys_corr_5: 7.850087868619188 + sys_corr_6: -88.07520308797234 + sys_corr_7: 137.25556071652076 + sys_corr_8: 283.50818110253323 + sys_corr_9: -239.42618180677908 + sys_corr_10: -114.9282210891482 + sys_corr_11: 8.42245032181839 + sys_corr_12: 31.98627693875795 + sys_corr_13: 1.50717484391528 + sys_corr_14: 6.896281410681029 + sys_corr_15: -7.595533001147195 + sys_corr_16: -7.633120097776864 + sys_corr_17: 1.0195549272803348 + sys_corr_18: 7.284609808103727e-07 + sys_corr_19: -1.6778355582638205e-05 + sys_corr_20: 0.00025537746071577504 + sys_corr_21: -0.0011339103675613135 + sys_corr_22: -0.003415602574189624 + sys_corr_23: -0.005932843196956676 + sys_corr_24: -0.1280796093801708 + sys_corr_25: -7.564443428452589e-06 + sys_corr_26: 0.016222470669920764 + sys_corr_27: 0.06514473289220417 + sys_corr_28: 0.10371048272556845 + sys_corr_29: 0.06929157468645339 + sys_corr_30: 0.018776298638962904 + sys_corr_31: -0.021881712529226115 + sys_corr_32: -0.036065189677757854 + sys_corr_33: 0.03731793154064317 + sys_corr_34: -0.006212574243613175 + stat: 1200.0 + sys_beam: 2400.0 + sys_luminosity: 2700.0 +- sys_corr_1: -757.4168560908336 + sys_corr_2: 75.90309841401225 + sys_corr_3: 1987.0588271704808 + sys_corr_4: 380.50217714972797 + sys_corr_5: 186.6540628119187 + sys_corr_6: -98.92177817595513 + sys_corr_7: 19.386151565620303 + sys_corr_8: 90.24316654080108 + sys_corr_9: 97.86762228902114 + sys_corr_10: 257.88630047369355 + sys_corr_11: 27.80542496923231 + sys_corr_12: 182.1786676656212 + sys_corr_13: -40.1069768929091 + sys_corr_14: 11.34420533632285 + sys_corr_15: -37.39982404709859 + sys_corr_16: -32.34179686209255 + sys_corr_17: 1.0943227332297263 + sys_corr_18: -6.383721004848815e-06 + sys_corr_19: -3.428582414792927e-05 + sys_corr_20: 0.00046593140222169854 + sys_corr_21: 0.0016542604394633518 + sys_corr_22: -0.0033912447442104652 + sys_corr_23: 0.004148891256763143 + sys_corr_24: -0.009165474667733913 + sys_corr_25: 0.008144508683747427 + sys_corr_26: -0.026150151611025615 + sys_corr_27: -0.061315511041762666 + sys_corr_28: 0.012687205157747163 + sys_corr_29: 0.018693498085528087 + sys_corr_30: 0.01425778424606732 + sys_corr_31: -0.006241271517890433 + sys_corr_32: -0.006706266150301394 + sys_corr_33: 0.025389499587724858 + sys_corr_34: 0.00543093342266935 + stat: 900.0 + sys_beam: 2100.0 + sys_luminosity: 2400.0 +- sys_corr_1: -539.1116598566842 + sys_corr_2: -1082.196953898048 + sys_corr_3: 139.26538263419215 + sys_corr_4: 259.0007882801631 + sys_corr_5: -1235.3638860478686 + sys_corr_6: -147.8992484251472 + sys_corr_7: -30.971965139079433 + sys_corr_8: 66.12783421119019 + sys_corr_9: 79.27645225016443 + sys_corr_10: 6.1083395365405675 + sys_corr_11: 298.8256052265815 + sys_corr_12: -160.33052814602755 + sys_corr_13: 123.74448329511169 + sys_corr_14: -32.43971184140428 + sys_corr_15: -11.811037446007145 + sys_corr_16: -23.921899411063976 + sys_corr_17: 0.5767811158596198 + sys_corr_18: -2.8777510060249715e-06 + sys_corr_19: -1.7746259062853234e-05 + sys_corr_20: -0.0015668674597882877 + sys_corr_21: 0.00047893537625031795 + sys_corr_22: 0.00013796288062932144 + sys_corr_23: 0.015193558356304633 + sys_corr_24: 0.08368951181574025 + sys_corr_25: 0.03252795209784002 + sys_corr_26: 0.012074739801121198 + sys_corr_27: -0.07568283304928537 + sys_corr_28: -0.10839682755241084 + sys_corr_29: -0.06354115477152299 + sys_corr_30: 0.022137948768943465 + sys_corr_31: -0.00040771647319118586 + sys_corr_32: 0.029190918104867233 + sys_corr_33: -0.045338286183720544 + sys_corr_34: 0.009281120599104669 + stat: 800.0 + sys_beam: 1800.0 + sys_luminosity: 2100.0 +- sys_corr_1: -346.6272090627615 + sys_corr_2: -1030.181994954958 + sys_corr_3: 2.3905943298608827 + sys_corr_4: -672.7845398398682 + sys_corr_5: 431.24137151214705 + sys_corr_6: -736.2297158517016 + sys_corr_7: -0.9211964079113855 + sys_corr_8: 167.18955777019715 + sys_corr_9: 125.27532107231038 + sys_corr_10: 0.6952330910211038 + sys_corr_11: 153.6669942345208 + sys_corr_12: -132.51248962046776 + sys_corr_13: -319.0057886061013 + sys_corr_14: 92.20109904875524 + sys_corr_15: 52.134955258314406 + sys_corr_16: 11.019143373701814 + sys_corr_17: 0.8035534599422568 + sys_corr_18: -1.9494090664902555e-06 + sys_corr_19: 0.0004555292374753845 + sys_corr_20: -0.0005602017511280195 + sys_corr_21: 0.0017116083298735739 + sys_corr_22: -0.0012111825214759626 + sys_corr_23: 0.013093740239201866 + sys_corr_24: 0.0663328708802538 + sys_corr_25: 0.005741342860060582 + sys_corr_26: -0.005243893280854078 + sys_corr_27: -0.014898575132749711 + sys_corr_28: -0.07832385774732964 + sys_corr_29: -0.008103685999407343 + sys_corr_30: -0.028478571390293136 + sys_corr_31: -0.02209874007861173 + sys_corr_32: 0.06915915493569391 + sys_corr_33: -0.05225712733244364 + sys_corr_34: -0.021446400930781855 + stat: 700.0 + sys_beam: 1500.0 + sys_luminosity: 1800.0 +- sys_corr_1: -274.8755548245817 + sys_corr_2: -712.8587225878063 + sys_corr_3: 9.424364874672388 + sys_corr_4: -310.81026743552223 + sys_corr_5: 162.63611732393898 + sys_corr_6: 728.7752580649696 + sys_corr_7: -76.45444548329307 + sys_corr_8: 508.0244085137416 + sys_corr_9: 158.27207582811076 + sys_corr_10: 47.22082110864786 + sys_corr_11: 34.177741491179496 + sys_corr_12: 13.152729580561742 + sys_corr_13: -117.7021162614959 + sys_corr_14: -377.1740405636019 + sys_corr_15: 48.51723655535921 + sys_corr_16: 31.56153128981498 + sys_corr_17: 0.9822152491275007 + sys_corr_18: -2.4619998369348916e-06 + sys_corr_19: 0.0005612970734051026 + sys_corr_20: 0.0005386362464870477 + sys_corr_21: -0.001815465250724656 + sys_corr_22: -0.006347847483481793 + sys_corr_23: 0.017227109060982387 + sys_corr_24: 0.09711895842625598 + sys_corr_25: 0.03514572275551089 + sys_corr_26: 0.04299205294753696 + sys_corr_27: -0.011167609711057146 + sys_corr_28: -0.04891378723110805 + sys_corr_29: -0.09846940785881475 + sys_corr_30: 0.02063997187136929 + sys_corr_31: -0.03877240980632347 + sys_corr_32: 0.06848351774573039 + sys_corr_33: -0.04952056140314971 + sys_corr_34: -0.003110697817064795 + stat: 600.0 + sys_beam: 1200.0 + sys_luminosity: 1400.0 +- sys_corr_1: -167.11074525470343 + sys_corr_2: -569.573314600226 + sys_corr_3: -47.732461471837524 + sys_corr_4: -61.80663852923446 + sys_corr_5: 100.02781413206776 + sys_corr_6: 138.327542900984 + sys_corr_7: -68.1364178650373 + sys_corr_8: -175.18878405870703 + sys_corr_9: -517.694932637838 + sys_corr_10: 375.78574105417636 + sys_corr_11: -22.358913549372588 + sys_corr_12: -171.54256789293123 + sys_corr_13: -77.7987263976552 + sys_corr_14: -20.62483812309415 + sys_corr_15: -333.83570107460514 + sys_corr_16: 141.46991248825637 + sys_corr_17: 0.48158322751296384 + sys_corr_18: 4.342515998397655e-06 + sys_corr_19: 0.00012465687701839452 + sys_corr_20: 0.0002705298702991312 + sys_corr_21: 0.0015800747762540462 + sys_corr_22: -0.006142676190768974 + sys_corr_23: -0.010540544502378445 + sys_corr_24: 0.09906325673279748 + sys_corr_25: 0.02814902488221871 + sys_corr_26: 0.004982507679248855 + sys_corr_27: 0.029057225727090522 + sys_corr_28: -0.06471540581684256 + sys_corr_29: -0.04520273973785426 + sys_corr_30: 0.05994446708871603 + sys_corr_31: 0.05655868692291631 + sys_corr_32: -0.01759802246321701 + sys_corr_33: 0.034334691922246696 + sys_corr_34: -0.009203113866061742 + stat: 500.0 + sys_beam: 800.0 + sys_luminosity: 1000.0 +- sys_corr_1: -188.20228910192097 + sys_corr_2: -777.7968807961512 + sys_corr_3: -83.523149571827 + sys_corr_4: -219.69021728485177 + sys_corr_5: 118.43402492877796 + sys_corr_6: 252.7137844343597 + sys_corr_7: 294.355545808125 + sys_corr_8: -370.94543280886097 + sys_corr_9: -209.10455100615636 + sys_corr_10: -130.42393782916918 + sys_corr_11: 365.6282260301084 + sys_corr_12: 480.6203850587949 + sys_corr_13: 20.079947758739483 + sys_corr_14: 67.06685990860872 + sys_corr_15: 65.32932344673472 + sys_corr_16: 87.47120757016614 + sys_corr_17: -0.09826105456667454 + sys_corr_18: 8.386171914222564e-06 + sys_corr_19: -0.0010224700247338724 + sys_corr_20: -0.0019943997996906787 + sys_corr_21: 0.0084136973091293 + sys_corr_22: -0.02061594320152601 + sys_corr_23: 0.03529448123666306 + sys_corr_24: 0.18401404863525153 + sys_corr_25: 0.059450908402464545 + sys_corr_26: 0.005920112777482507 + sys_corr_27: 0.038432895713877376 + sys_corr_28: 0.01778891141365446 + sys_corr_29: -0.04235602286263034 + sys_corr_30: 0.019926253823197022 + sys_corr_31: 0.017823211055328757 + sys_corr_32: -0.028407360893841153 + sys_corr_33: -0.015168614203903776 + sys_corr_34: -0.00999576404121969 + stat: 500.0 + sys_beam: 900.0 + sys_luminosity: 1000.0 +- sys_corr_1: -20.418282214869066 + sys_corr_2: -294.128388854442 + sys_corr_3: -58.60919723330269 + sys_corr_4: -198.41896719160766 + sys_corr_5: 53.836537785331174 + sys_corr_6: 94.56972392250604 + sys_corr_7: -235.80701606356234 + sys_corr_8: -98.75510182567812 + sys_corr_9: -143.68649816865596 + sys_corr_10: -129.97796340553987 + sys_corr_11: 34.287470713688585 + sys_corr_12: 106.63233272961433 + sys_corr_13: -66.03659873895782 + sys_corr_14: -27.168441642658422 + sys_corr_15: -173.44340111205508 + sys_corr_16: -444.7463117440774 + sys_corr_17: -0.24157992028818057 + sys_corr_18: -1.9035352315276237e-05 + sys_corr_19: -0.0020680211971272366 + sys_corr_20: -0.002768681927241517 + sys_corr_21: 0.008251447625955019 + sys_corr_22: -0.017362115654713663 + sys_corr_23: 0.0005594231669970876 + sys_corr_24: 0.00048017663069296006 + sys_corr_25: 0.041775962924260265 + sys_corr_26: -0.0019150653015100835 + sys_corr_27: -0.03489947235145362 + sys_corr_28: 0.16155295845034465 + sys_corr_29: 0.045016400588073914 + sys_corr_30: -0.031051031050705183 + sys_corr_31: 0.007999837636393914 + sys_corr_32: -0.1345409740234758 + sys_corr_33: 0.08854643213947719 + sys_corr_34: -0.03573913245167926 + stat: 400.0 + sys_beam: 200.0 + sys_luminosity: 300.0 +- sys_corr_1: -1418.5660286961804 + sys_corr_2: -300.83618976028646 + sys_corr_3: -372.17857782986243 + sys_corr_4: 110.85229077349138 + sys_corr_5: 213.07294567190755 + sys_corr_6: 192.03560835810046 + sys_corr_7: -315.0880721533769 + sys_corr_8: -688.389231807664 + sys_corr_9: 537.9263076607414 + sys_corr_10: 188.2713964977463 + sys_corr_11: -46.3728491490399 + sys_corr_12: -58.69340052788875 + sys_corr_13: -14.504876419343034 + sys_corr_14: -14.744707792059215 + sys_corr_15: 3.8807055752414743 + sys_corr_16: 8.66452881632407 + sys_corr_17: 0.8337970308503326 + sys_corr_18: -3.098814478249894e-06 + sys_corr_19: -8.957681649067556e-05 + sys_corr_20: -0.0002853717212049436 + sys_corr_21: 0.0010659257978352338 + sys_corr_22: 0.0037092495827079847 + sys_corr_23: 0.006878015330544949 + sys_corr_24: -0.19054727226453752 + sys_corr_25: -0.024486657565609662 + sys_corr_26: 0.0053534274545727146 + sys_corr_27: 0.08616988285704642 + sys_corr_28: 0.14996497858208171 + sys_corr_29: 0.11390900690940575 + sys_corr_30: -0.03560935114680982 + sys_corr_31: 0.012681825895916057 + sys_corr_32: -0.03089693251724615 + sys_corr_33: 0.05081019719485488 + sys_corr_34: 0.052787600581333 + stat: 900.0 + sys_beam: 1200.0 + sys_luminosity: 1600.0 +- sys_corr_1: -379.7758627502127 + sys_corr_2: -614.702742631885 + sys_corr_3: 806.9723200388605 + sys_corr_4: 129.34897826869812 + sys_corr_5: 329.5697346931343 + sys_corr_6: 180.1001765648686 + sys_corr_7: -83.25772608032091 + sys_corr_8: -249.97342476735562 + sys_corr_9: -217.19952916877614 + sys_corr_10: -625.126074096715 + sys_corr_11: -89.05430339410013 + sys_corr_12: -320.1639682551919 + sys_corr_13: 70.13671643344883 + sys_corr_14: -28.02170633911008 + sys_corr_15: 56.046822735985636 + sys_corr_16: 47.589118060404786 + sys_corr_17: 0.4556251546296614 + sys_corr_18: 1.1495839610964414e-05 + sys_corr_19: 0.00013258230364422347 + sys_corr_20: 0.0001381739901380239 + sys_corr_21: -0.005874689172628592 + sys_corr_22: 0.010501500264532582 + sys_corr_23: -0.01794623114827708 + sys_corr_24: -0.21644468180010468 + sys_corr_25: -0.03488192226887612 + sys_corr_26: -0.06857859750272072 + sys_corr_27: -0.049221237602540484 + sys_corr_28: 0.07731091052149563 + sys_corr_29: 0.07604399926457621 + sys_corr_30: -0.025643183602270476 + sys_corr_31: -0.027707930572605736 + sys_corr_32: -0.039392965654304765 + sys_corr_33: 0.008744923900361117 + sys_corr_34: 0.005484652230440729 + stat: 700.0 + sys_beam: 1000.0 + sys_luminosity: 1400.0 +- sys_corr_1: -239.61349262110903 + sys_corr_2: -678.8206626587207 + sys_corr_3: 43.04893366922258 + sys_corr_4: -55.83800780748894 + sys_corr_5: -634.499154269975 + sys_corr_6: -25.947554268123866 + sys_corr_7: 41.427904422469275 + sys_corr_8: -54.053308811602285 + sys_corr_9: -16.362379243606593 + sys_corr_10: -129.38273499463483 + sys_corr_11: -592.7873743872755 + sys_corr_12: 256.3760680903285 + sys_corr_13: -265.5698980460149 + sys_corr_14: 43.29483884631226 + sys_corr_15: -1.6211538751018917 + sys_corr_16: 48.76481505905378 + sys_corr_17: 0.809829264135287 + sys_corr_18: -2.8081644487259427e-06 + sys_corr_19: 0.00010808868172066936 + sys_corr_20: 0.001635589344699784 + sys_corr_21: -0.00017695688461365108 + sys_corr_22: -0.0018627281863186546 + sys_corr_23: -0.007389095578738441 + sys_corr_24: -0.09727117032053585 + sys_corr_25: -0.004215152100620265 + sys_corr_26: -0.00531890952950511 + sys_corr_27: -0.013228979704683317 + sys_corr_28: 0.002784867494607145 + sys_corr_29: 0.07434987123378306 + sys_corr_30: -0.033407111180509784 + sys_corr_31: 0.020745837975215568 + sys_corr_32: -0.03494504746873362 + sys_corr_33: 0.04376355950975352 + sys_corr_34: -0.02172061471681362 + stat: 600.0 + sys_beam: 1000.0 + sys_luminosity: 1300.0 +- sys_corr_1: -222.80016800556484 + sys_corr_2: -742.5368422506577 + sys_corr_3: -36.8272748602642 + sys_corr_4: -408.192048198211 + sys_corr_5: 263.337600274069 + sys_corr_6: -456.734990608118 + sys_corr_7: -7.855868883315059 + sys_corr_8: 78.7247284309053 + sys_corr_9: 44.68970834520459 + sys_corr_10: 42.11339132489611 + sys_corr_11: -238.17127568567935 + sys_corr_12: 147.6273626390351 + sys_corr_13: 522.6255945867229 + sys_corr_14: -164.03665899925022 + sys_corr_15: -66.06062046337215 + sys_corr_16: 10.446280909012556 + sys_corr_17: 0.7081767021411348 + sys_corr_18: -4.695676677119722e-06 + sys_corr_19: -0.0001868136902926309 + sys_corr_20: 0.000835749485750196 + sys_corr_21: -0.0014647779094920056 + sys_corr_22: 0.0005052468688517873 + sys_corr_23: 0.01511537237976092 + sys_corr_24: -0.004813289205308594 + sys_corr_25: 0.024600067592098613 + sys_corr_26: 0.02837806872979211 + sys_corr_27: -0.002254224212950835 + sys_corr_28: -0.053259457136311486 + sys_corr_29: -0.04791596442701525 + sys_corr_30: 0.0681923382744364 + sys_corr_31: 0.04485150499439751 + sys_corr_32: 2.8937518308521544e-05 + sys_corr_33: 0.00247769781776522 + sys_corr_34: 0.01624354351478134 + stat: 600.0 + sys_beam: 900.0 + sys_luminosity: 1200.0 +- sys_corr_1: -221.03406480331822 + sys_corr_2: -598.7026972313603 + sys_corr_3: 7.879832997776385 + sys_corr_4: -218.47633151977615 + sys_corr_5: 119.79805532288195 + sys_corr_6: 540.0535716310895 + sys_corr_7: -14.335697968361968 + sys_corr_8: 370.33754643475316 + sys_corr_9: 156.56323982872624 + sys_corr_10: 48.85089017217374 + sys_corr_11: -72.25805431862382 + sys_corr_12: -60.236085152793216 + sys_corr_13: 166.24096100027964 + sys_corr_14: 496.04191921051756 + sys_corr_15: -42.90701380970265 + sys_corr_16: -16.03726295711185 + sys_corr_17: 0.4785205369484593 + sys_corr_18: -4.710614020048994e-06 + sys_corr_19: -9.20596815872463e-05 + sys_corr_20: -0.0012894448180824792 + sys_corr_21: 0.005305379878287117 + sys_corr_22: -0.004502053808523869 + sys_corr_23: 0.009383386501251282 + sys_corr_24: 0.1241067869378707 + sys_corr_25: 0.03652483327628769 + sys_corr_26: -0.026124954404075508 + sys_corr_27: 0.027004757545574683 + sys_corr_28: -0.022177739726991944 + sys_corr_29: -0.04456292446841822 + sys_corr_30: 0.02875750635448535 + sys_corr_31: 0.05203761053449261 + sys_corr_32: 0.0586078660639764 + sys_corr_33: -0.05914759281192592 + sys_corr_34: 0.0033480429937501 + stat: 600.0 + sys_beam: 800.0 + sys_luminosity: 1100.0 +- sys_corr_1: -175.95234830648135 + sys_corr_2: -458.29719201587636 + sys_corr_3: -25.221580090315847 + sys_corr_4: 34.36989661893009 + sys_corr_5: 36.823678063795654 + sys_corr_6: 68.23139353147354 + sys_corr_7: -45.66847954320662 + sys_corr_8: -111.32792227020646 + sys_corr_9: -424.8650642158514 + sys_corr_10: 408.5006539031898 + sys_corr_11: -128.01638358132988 + sys_corr_12: -105.17863116654347 + sys_corr_13: 49.58952197883063 + sys_corr_14: 17.210940422189474 + sys_corr_15: 392.3340808181992 + sys_corr_16: -123.20759730685347 + sys_corr_17: 0.7580998432848071 + sys_corr_18: 5.007363639324624e-06 + sys_corr_19: -0.00036884888231146186 + sys_corr_20: -0.0012978386362287946 + sys_corr_21: 0.005410694702775865 + sys_corr_22: -0.013313942732979372 + sys_corr_23: 0.05115695988736416 + sys_corr_24: 0.23458858072676236 + sys_corr_25: 0.06360147645469767 + sys_corr_26: 0.013309433725921567 + sys_corr_27: -0.04071310626280849 + sys_corr_28: -0.028177852122879966 + sys_corr_29: -0.1799505617141715 + sys_corr_30: -0.02169242324782223 + sys_corr_31: -0.007009667327721166 + sys_corr_32: 0.03717354123857821 + sys_corr_33: -0.12043966335792472 + sys_corr_34: 0.01543371410455415 + stat: 500.0 + sys_beam: 700.0 + sys_luminosity: 900.0 +- sys_corr_1: -183.14713236431768 + sys_corr_2: 567.3613770804149 + sys_corr_3: 468.56312146099094 + sys_corr_4: -942.5696139985873 + sys_corr_5: -309.77329328945194 + sys_corr_6: 193.68459039286637 + sys_corr_7: 711.2139889472766 + sys_corr_8: -281.5128856875479 + sys_corr_9: 108.21391439053171 + sys_corr_10: 77.90499127010314 + sys_corr_11: -58.041274291234444 + sys_corr_12: -198.5639555860381 + sys_corr_13: 23.54786254952854 + sys_corr_14: -43.82890699205783 + sys_corr_15: -20.50657414285486 + sys_corr_16: -67.43375805284023 + sys_corr_17: 1.4190064875049369 + sys_corr_18: -3.892991343602838e-06 + sys_corr_19: -0.00017797671767751046 + sys_corr_20: -0.0017744229331906653 + sys_corr_21: 0.005908187636415936 + sys_corr_22: -0.014350998864706447 + sys_corr_23: 0.02762243269222977 + sys_corr_24: 0.18750360514629275 + sys_corr_25: 0.05462605794528981 + sys_corr_26: 0.010284578419506411 + sys_corr_27: -0.040815268044172244 + sys_corr_28: -0.0392171074296804 + sys_corr_29: -0.1371256172369289 + sys_corr_30: 0.017759500556972105 + sys_corr_31: 0.06695235138005146 + sys_corr_32: -0.00374669242826696 + sys_corr_33: -0.061967289524856445 + sys_corr_34: -0.019287334663180073 + stat: 600.0 + sys_beam: 1000.0 + sys_luminosity: 1400.0 +- sys_corr_1: -108.57118702106705 + sys_corr_2: 839.8291501072293 + sys_corr_3: 459.8390870183657 + sys_corr_4: -925.6968612905044 + sys_corr_5: -400.4040360456309 + sys_corr_6: 42.68917936257094 + sys_corr_7: -752.2174879814365 + sys_corr_8: -53.07108189474481 + sys_corr_9: -128.09208022343634 + sys_corr_10: -32.98357281980944 + sys_corr_11: 74.02432402617801 + sys_corr_12: 78.5267575038786 + sys_corr_13: 36.42357651871526 + sys_corr_14: 39.17038903485109 + sys_corr_15: 53.70701716655165 + sys_corr_16: 94.08338628291777 + sys_corr_17: 0.6888256080273181 + sys_corr_18: 4.136991764242596e-07 + sys_corr_19: -0.000857665038759362 + sys_corr_20: -0.0026926244135948014 + sys_corr_21: 0.0075394244035477806 + sys_corr_22: -0.01576189784794336 + sys_corr_23: 0.04972357343844988 + sys_corr_24: 0.14135531791323716 + sys_corr_25: 0.08199691279107951 + sys_corr_26: 0.010188403737656443 + sys_corr_27: 0.02352486740748513 + sys_corr_28: -0.04847116685701623 + sys_corr_29: -0.03684786732394954 + sys_corr_30: -0.0070126705255896095 + sys_corr_31: 0.01275766566589463 + sys_corr_32: 0.022736517156611544 + sys_corr_33: -0.04112496966862067 + sys_corr_34: 0.018462409568420584 + stat: 700.0 + sys_beam: 500.0 + sys_luminosity: 700.0 diff --git a/buildmaster/LHCB_WPWM_8TEV_MUON/data.yaml b/buildmaster/LHCB_WPWM_8TEV_MUON/data.yaml new file mode 100644 index 0000000000..d87f6fe500 --- /dev/null +++ b/buildmaster/LHCB_WPWM_8TEV_MUON/data.yaml @@ -0,0 +1,17 @@ +data_central: +- 236500.0 +- 208400.0 +- 182000.0 +- 153300.0 +- 119500.0 +- 84400.0 +- 86400.0 +- 23000.0 +- 134000.0 +- 119800.0 +- 110600.0 +- 102400.0 +- 92500.0 +- 79900.0 +- 119300.0 +- 60000.0 diff --git a/buildmaster/LHCB_WPWM_8TEV_MUON/filter.py b/buildmaster/LHCB_WPWM_8TEV_MUON/filter.py new file mode 100644 index 0000000000..dfd2869217 --- /dev/null +++ b/buildmaster/LHCB_WPWM_8TEV_MUON/filter.py @@ -0,0 +1,390 @@ +import numpy as np +import pandas as pd +import pathlib +import yaml + +from validphys.commondata_utils import covmat_to_artunc + + +MZ_VALUE = 91.1876 # GeV +MW_VALUE = 80.398 # GeV +SQRT_S = 7_000.0 # GeV +NORM_FACTOR = 1_000.0 # from pb -> fb + +# Correct tables to read values [[Z], [W+, W-]] +TABLES = {2: [0], 1: [0, 2]} # {table_id: [indexes]} +NB_POINTS = [18, 8, 8] # [N_Z, N_W+, N_W-] + +# MAP Boson string to M values +MAP_BOSON = {"Z": MZ_VALUE, "W": MW_VALUE} +MAP_TABLE = {2: "Z", 1: "W"} + + +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-ins1406555-v{version}-Table_{table_id}" + table = pathlib.Path(f"./rawdata/{filename}.yaml") + + return yaml.safe_load(table.read_text()) + + +def read_metadata() -> tuple[int, int, list]: + """Read the version and list of tables from metadata. + + Returns + ------- + tuple(int, list): + data version and list of hepdata tables + + """ + metadata = pathlib.Path("./metadata.yaml") + content = yaml.safe_load(metadata.read_text()) + + version = content["hepdata"]["version"] + nb_datapoints = sum(content["implemented_observables"][0]["npoints"]) + tables = content["implemented_observables"][0]["tables"] + + return version, nb_datapoints, tables + + +def get_kinematics(hepdata: dict, bin_index: list, boson: str = "Z") -> 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 + + Returns + ------- + tuple(int, list): + data version and list of hepdata tables + + """ + rapbins = hepdata["independent_variables"][0]["values"] + + kinematics = [] + for bins in bin_index: + ymin = float(rapbins[bins]["low"]) + ymax = float(rapbins[bins]["high"]) + kin_value = { + "y": { + "min": ymin, + "mid": 0.5 * (ymin + ymax), + "max": ymax, + }, + "M2": {"min": None, "mid": MAP_BOSON[boson] ** 2, "max": None}, + "sqrt_s": {"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 + list of Non-empty bin index + idx: int + 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 [NORM_FACTOR * central[i]["value"] for i in bin_index] + + +def get_errors(hepdata: dict, bin_index: list, indx: int = 0) -> dict: + """Extract the error values from the HepData yaml file. + + Parameters + ---------- + hepdata: dict + dictionary containing all data info + bin_index: list + list of Non-empty bin index + indx: int + index from which the errors will be read + + Returns + ------- + list: + list of dictionaries whose contents are the various + source of uncertainties + + """ + errors = hepdata["dependent_variables"][indx]["values"] + + stat, sys_corr, sys_beam, sys_lumi = [], [], [], [] + for idx in bin_index: + stat.append(NORM_FACTOR * errors[idx]["errors"][0]["symerror"]) + sys_corr.append(NORM_FACTOR * errors[idx]["errors"][1]["symerror"]) + sys_beam.append(NORM_FACTOR * errors[idx]["errors"][2]["symerror"]) + sys_lumi.append(NORM_FACTOR * errors[idx]["errors"][3]["symerror"]) + + return { + "stat": stat, + "sys_corr": sys_corr, + "sys_beam": sys_beam, + "sys_lumi": sys_lumi, + } + + +def read_corrmatrix(nb_datapoints: int) -> np.ndarray: + """Read the matrix and returns a symmetrized verions. + + Parameters + ---------- + nb_datapoints: int + total number of datapoints + + Returns + ------- + np.ndarray: (nb_datapoints, nb_datapoints) + a symmetrized square matrix + + """ + corrmat = pd.read_csv( + "./rawdata/corrmat.corr", + names=[f'{i}' for i in range(nb_datapoints)], + delim_whitespace=True, + ) + corrmat = corrmat.iloc[:, :].values + + # Symmetrize the Correlation Matrix + symmetrised = np.triu(corrmat.T) + np.triu(corrmat.T, 1).T + assert np.allclose(symmetrised, symmetrised.T) + + return symmetrised + + +def multiply_syst(covmat: np.ndarray, sys: list) -> np.ndarray: + """Multiply the covmat with the total systematics. + + Parameters + ---------- + covmat: np.ndarray + a squared matrix representing the covmat + sys: list + a vector containing the correlatdd systematics + + Returns + ------- + np.ndarray: + covariance matrix multiplied by the total systematics + + """ + covmat_sys = np.zeros(covmat.shape) + + for i in range(covmat.shape[0]): + for j in range(covmat.shape[-1]): + covmat_sys[i][j] = covmat[i][j] * sys[i] * sys[j] + + return covmat_sys.flatten() + + +def generate_artificial_unc(**kwargs) -> np.ndarray: + """A wrapper around `covmat_to_artunc` to return an array. + + Returns + ------- + np.ndarray: + a squared matrix with artificial uncertainties + """ + artunc = covmat_to_artunc(**kwargs) + return np.array(artunc) + + +def concatenate_dicts(multidict: list[dict]) -> dict: + """Given a list of dictionaries combined the values. + + Parameters + ---------- + multidict: list[dict] + list of dictionaries with the same keys + + Returns + ------- + dict: + dictionary whose keys are combined + """ + new_dict = {} + for key in multidict[0].keys(): + new_dict[key] = [] + for element in multidict: + new_dict[key] += element[key] + + return new_dict + + +def format_uncertainties(uncs: dict, artunc: np.ndarray) -> list: + """Format the uncertainties to be dumped into the yaml file. + + Parameters + ---------- + uncs: dict + dictionary containing the various source of uncertainties + artunc: no.ndarray + array of the artificial systematic uncertainties + + Returns + ------- + list: + list of ditionaries whose elements are the various errors + + """ + combined_errors = [] + + artunc = artunc.tolist() + for idat in range(len(artunc)): + error_value = {} + for jdat in range(len(artunc[idat])): + error_value[f"sys_corr_{jdat + 1}"] = artunc[idat][jdat] + + error_value["stat"] = uncs["stat"][idat] + error_value["sys_beam"] = uncs["sys_beam"][idat] + error_value["sys_luminosity"] = uncs["sys_lumi"][idat] + combined_errors.append(error_value) + + return combined_errors + + +def dump_commondata(kinematics: list, data: list, errors: list) -> 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: list + list containing the different errors + + """ + error_definition = { + f"sys_corr_{i + 1}": { + "description": "Correlated systematic uncertainties", + "treatment": "ADD", + "type": "CORR", + } + for i in range(len(data)) + } + + error_definition["stat"] = { + "description": "Uncorrelated statistical uncertainties", + "treatment": "ADD", + "type": "UNCORR", + } + + error_definition["sys_beam"] = { + "description": "Systematic Beam uncertainties", + "treatment": "ADD", + "type": "LHCBBEAM7TEV", + } + + error_definition["sys_luminosity"] = { + "description": "Systematic Luminosity uncertainties", + "treatment": "ADD", + "type": "LHCBLUMI7TEV", + } + + with open("data.yaml", "w") as file: + yaml.dump({"data_central": data}, file, sort_keys=False) + + with open("kinematics.yaml", "w") as file: + yaml.dump({"bins": kinematics}, file, sort_keys=False) + + with open("uncertainties.yaml", "w") as file: + yaml.dump( + {"definitions": error_definition, "bins": errors}, + file, + sort_keys=False, + ) + + +def main_filter(boson: str = "Z") -> None: + """Main driver of the LHCB_ZEE_8TEV filter that produces commmondata. + + There are three main different sources of uncertainties. + + 1. Statistical uncertainties: ADD, UNCORR + + 2. Correlated Systematic uncertainties: ADD, CORR: + constructed by symmetrizing the correlation matrix and extracting + the artificial systematic uncertainties from the latter + + 3. Beam Systematic uncertainties: ADD, LHCBBEAM7TEV + + 4. Luminosity Systematic uncertainties: ADD, LHCBLUMI7TEV + + """ + version, _, _ = read_metadata() + nbpoints = sum(NB_POINTS) + + # Select only the bins corresponding to the correct boson + bslice = slice(18) if boson == "Z" else slice(18, nbpoints + 1) + + nbp_idx = 0 + comb_kins, comb_data = [], [] + combined_errors = [] + for tabid in TABLES.keys(): # Loop over tables + for idx in TABLES[tabid]: # Loop over Bosons [Z, W+, W-] + bin_index = [i for i in range(NB_POINTS[nbp_idx])] + yaml_content = load_yaml(table_id=tabid, version=version) + + # Extract the kinematic, data, and uncertainties + kinematics = get_kinematics(yaml_content, bin_index, MAP_TABLE[tabid]) + data_central = get_data_values(yaml_content, bin_index, indx=idx) + uncertainties = get_errors(yaml_content, bin_index, indx=idx) + + # Collect all the results from different tables + comb_kins += kinematics + comb_data += data_central + combined_errors.append(uncertainties) + + nbp_idx += 1 + + errors_combined = concatenate_dicts(combined_errors) + # Compute the Artifical Systematics from CovMat + corrmat = read_corrmatrix(nb_datapoints=nbpoints) + covmat = multiply_syst(corrmat, errors_combined["sys_corr"]) + artunc = generate_artificial_unc( + ndata=nbpoints, + covmat_list=covmat.tolist(), + no_of_norm_mat=0, + ) + errors = format_uncertainties(errors_combined, artunc) + + # Generate all the necessary files + dump_commondata(comb_kins[bslice], comb_data[bslice], errors[bslice]) + + return + + +if __name__ == "__main__": + main_filter(boson="W") diff --git a/buildmaster/LHCB_WPWM_8TEV_MUON/kinematics.yaml b/buildmaster/LHCB_WPWM_8TEV_MUON/kinematics.yaml new file mode 100644 index 0000000000..278ea1b812 --- /dev/null +++ b/buildmaster/LHCB_WPWM_8TEV_MUON/kinematics.yaml @@ -0,0 +1,193 @@ +bins: +- y: + min: 2.0 + mid: 2.125 + max: 2.25 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.25 + mid: 2.375 + max: 2.5 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.5 + mid: 2.625 + max: 2.75 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.75 + mid: 2.875 + max: 3.0 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.0 + mid: 3.125 + max: 3.25 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.25 + mid: 3.375 + max: 3.5 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.5 + mid: 3.75 + max: 4.0 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 4.0 + mid: 4.25 + max: 4.5 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.0 + mid: 2.125 + max: 2.25 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.25 + mid: 2.375 + max: 2.5 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.5 + mid: 2.625 + max: 2.75 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.75 + mid: 2.875 + max: 3.0 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.0 + mid: 3.125 + max: 3.25 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.25 + mid: 3.375 + max: 3.5 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.5 + mid: 3.75 + max: 4.0 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 4.0 + mid: 4.25 + max: 4.5 + M2: + min: null + mid: 6463.838403999999 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null diff --git a/buildmaster/LHCB_WPWM_8TEV_MUON/metadata.yaml b/buildmaster/LHCB_WPWM_8TEV_MUON/metadata.yaml new file mode 100644 index 0000000000..a63d7ed7c4 --- /dev/null +++ b/buildmaster/LHCB_WPWM_8TEV_MUON/metadata.yaml @@ -0,0 +1,48 @@ +setname: "LHCB_WPWM_8TEV_MUON" + +version: 1 +version_comment: "First implementation" + +# References +arXiv: + url: "https://arxiv.org/abs/1511.08039" +iNSPIRE: + url: "https://inspirehep.net/literature/1406555" +hepdata: + url: "https://www.hepdata.net/record/ins1406555" + version: 1 + +nnpdf_metadata: + nnpdf31_process: "DY CC" + experiment: "LHCb" + +implemented_observables: + +- observable_name: "y" + observable: + description: "Inclusive cross-section for W+ and W- boson production in bins of muon pseudorapidity at 8 TeV" + label: r"$d\sigma / d|y|$" + units: "[fb]" + process_type: "EWK_RAP" + tables: [1] + npoints: [16] + plotting: + kinematics_override: ewk_rap_sqrt_scale + dataset_label: "LHCb $W \\to \\mu$ 8 TeV" + plot_x: y + y_label: '$d\sigma/dy$ (fb)' + kinematic_coverage: [y, M2, sqrt_s] + kinematics: + variables: + y: {description: "Muon rapidity", label: "$y$", units: ""} + M2: {description: "W boson Mass", label: "$M^2$", units: "$GeV^2$"} + sqrt_s: {description: "Center of Mass Energy", label: '$\sqrt{s}$', units: "$GeV$"} + file: kinematics.yaml + data_central: data.yaml + data_uncertainties: + - uncertainties.yaml + theory: + FK_tables: + - - LHCB_WP_8TEV + - LHCB_WM_8TEV + operation: 'null' diff --git a/buildmaster/LHCB_WPWM_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_1.yaml b/buildmaster/LHCB_WPWM_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_1.yaml new file mode 100644 index 0000000000..0332da55b0 --- /dev/null +++ b/buildmaster/LHCB_WPWM_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_1.yaml @@ -0,0 +1,184 @@ +dependent_variables: +- header: {name: '$\sigma_{W^{+} \to \mu^{+}\nu}$', units: PB} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: P P --> W+ < MU+ NU > X} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {label: stat, symerror: 1.2} + - {label: sys, symerror: 3.2} + - {label: 'sys,beam', symerror: 2.4} + - {label: 'sys,lumi', symerror: 2.7} + value: 236.5 + - errors: + - {label: stat, symerror: 0.9} + - {label: sys, symerror: 2.2} + - {label: 'sys,beam', symerror: 2.1} + - {label: 'sys,lumi', symerror: 2.4} + value: 208.4 + - errors: + - {label: stat, symerror: 0.8} + - {label: sys, symerror: 1.8} + - {label: 'sys,beam', symerror: 1.8} + - {label: 'sys,lumi', symerror: 2.1} + value: 182.0 + - errors: + - {label: stat, symerror: 0.7} + - {label: sys, symerror: 1.6} + - {label: 'sys,beam', symerror: 1.5} + - {label: 'sys,lumi', symerror: 1.8} + value: 153.3 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.3} + - {label: 'sys,beam', symerror: 1.2} + - {label: 'sys,lumi', symerror: 1.4} + value: 119.5 + - errors: + - {label: stat, symerror: 0.5} + - {label: sys, symerror: 1.0} + - {label: 'sys,beam', symerror: 0.8} + - {label: 'sys,lumi', symerror: 1.0} + value: 84.4 + - errors: + - {label: stat, symerror: 0.5} + - {label: sys, symerror: 1.2} + - {label: 'sys,beam', symerror: 0.9} + - {label: 'sys,lumi', symerror: 1.0} + value: 86.4 + - errors: + - {label: stat, symerror: 0.4} + - {label: sys, symerror: 0.7} + - {label: 'sys,beam', symerror: 0.2} + - {label: 'sys,lumi', symerror: 0.3} + value: 23.0 +- header: {name: '$f^{W^+}_{FSR}$'} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: ''} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {symerror: 0.0047} + value: 1.0188 + - errors: + - {symerror: 0.0028} + value: 1.0163 + - errors: + - {symerror: 0.0025} + value: 1.0158 + - errors: + - {symerror: 0.0028} + value: 1.0148 + - errors: + - {symerror: 0.0032} + value: 1.0152 + - errors: + - {symerror: 0.0046} + value: 1.015 + - errors: + - {symerror: 0.0045} + value: 1.0175 + - errors: + - {symerror: 0.0087} + value: 1.0211 +- header: {name: '$\sigma_{W^{-} \to \mu^{-}\overline{\nu}}$', units: PB} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: P P --> W- < MU- NUBAR > X} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {label: stat, symerror: 0.9} + - {label: sys, symerror: 1.8} + - {label: 'sys,beam', symerror: 1.2} + - {label: 'sys,lumi', symerror: 1.6} + value: 134.0 + - errors: + - {label: stat, symerror: 0.7} + - {label: sys, symerror: 1.4} + - {label: 'sys,beam', symerror: 1.0} + - {label: 'sys,lumi', symerror: 1.4} + value: 119.8 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.2} + - {label: 'sys,beam', symerror: 1.0} + - {label: 'sys,lumi', symerror: 1.3} + value: 110.6 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.2} + - {label: 'sys,beam', symerror: 0.9} + - {label: 'sys,lumi', symerror: 1.2} + value: 102.4 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.1} + - {label: 'sys,beam', symerror: 0.8} + - {label: 'sys,lumi', symerror: 1.1} + value: 92.5 + - errors: + - {label: stat, symerror: 0.5} + - {label: sys, symerror: 0.9} + - {label: 'sys,beam', symerror: 0.7} + - {label: 'sys,lumi', symerror: 0.9} + value: 79.9 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.5} + - {label: 'sys,beam', symerror: 1.0} + - {label: 'sys,lumi', symerror: 1.4} + value: 119.3 + - errors: + - {label: stat, symerror: 0.7} + - {label: sys, symerror: 1.6} + - {label: 'sys,beam', symerror: 0.5} + - {label: 'sys,lumi', symerror: 0.7} + value: 60.0 +- header: {name: '$f^{W^-}_{FSR}$'} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: ''} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {symerror: 0.0026} + value: 1.0172 + - errors: + - {symerror: 0.0027} + value: 1.0155 + - errors: + - {symerror: 0.0028} + value: 1.0153 + - errors: + - {symerror: 0.003} + value: 1.0162 + - errors: + - {symerror: 0.0031} + value: 1.016 + - errors: + - {symerror: 0.0033} + value: 1.0176 + - errors: + - {symerror: 0.0033} + value: 1.02 + - errors: + - {symerror: 0.0053} + value: 1.0243 +independent_variables: +- header: {name: '$\eta^{\mu}$'} + values: + - {high: 2.25, low: 2.0} + - {high: 2.5, low: 2.25} + - {high: 2.75, low: 2.5} + - {high: 3.0, low: 2.75} + - {high: 3.25, low: 3.0} + - {high: 3.5, low: 3.25} + - {high: 4.0, low: 3.5} + - {high: 4.5, low: 4.0} diff --git a/buildmaster/LHCB_WPWM_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_2.yaml b/buildmaster/LHCB_WPWM_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_2.yaml new file mode 100644 index 0000000000..9a39670c3b --- /dev/null +++ b/buildmaster/LHCB_WPWM_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_2.yaml @@ -0,0 +1,200 @@ +dependent_variables: +- header: {name: '$\sigma_{Z \to \mu^{+}\mu^{-}}$', units: PB} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: M(MU+ MU-), units: GeV, value: 60.0-120.0} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: P P --> Z0 < MU+ MU- > X} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {label: stat, symerror: 0.033} + - {label: sys, symerror: 0.055} + - {label: 'sys,beam', symerror: 0.014} + - {label: 'sys,lumi', symerror: 0.014} + value: 1.223 + - errors: + - {label: stat, symerror: 0.051} + - {label: sys, symerror: 0.06} + - {label: 'sys,beam', symerror: 0.038} + - {label: 'sys,lumi', symerror: 0.038} + value: 3.263 + - errors: + - {label: stat, symerror: 0.062} + - {label: sys, symerror: 0.064} + - {label: 'sys,beam', symerror: 0.057} + - {label: 'sys,lumi', symerror: 0.058} + value: 4.983 + - errors: + - {label: stat, symerror: 0.07} + - {label: sys, symerror: 0.072} + - {label: 'sys,beam', symerror: 0.077} + - {label: 'sys,lumi', symerror: 0.078} + value: 6.719 + - errors: + - {label: stat, symerror: 0.076} + - {label: sys, symerror: 0.074} + - {label: 'sys,beam', symerror: 0.093} + - {label: 'sys,lumi', symerror: 0.094} + value: 8.051 + - errors: + - {label: stat, symerror: 0.079} + - {label: sys, symerror: 0.074} + - {label: 'sys,beam', symerror: 0.103} + - {label: 'sys,lumi', symerror: 0.105} + value: 8.967 + - errors: + - {label: stat, symerror: 0.081} + - {label: sys, symerror: 0.076} + - {label: 'sys,beam', symerror: 0.11} + - {label: 'sys,lumi', symerror: 0.112} + value: 9.561 + - errors: + - {label: stat, symerror: 0.082} + - {label: sys, symerror: 0.071} + - {label: 'sys,beam', symerror: 0.113} + - {label: 'sys,lumi', symerror: 0.115} + value: 9.822 + - errors: + - {label: stat, symerror: 0.081} + - {label: sys, symerror: 0.074} + - {label: 'sys,beam', symerror: 0.112} + - {label: 'sys,lumi', symerror: 0.114} + value: 9.721 + - errors: + - {label: stat, symerror: 0.078} + - {label: sys, symerror: 0.071} + - {label: 'sys,beam', symerror: 0.104} + - {label: 'sys,lumi', symerror: 0.105} + value: 9.03 + - errors: + - {label: stat, symerror: 0.072} + - {label: sys, symerror: 0.074} + - {label: 'sys,beam', symerror: 0.089} + - {label: 'sys,lumi', symerror: 0.09} + value: 7.748 + - errors: + - {label: stat, symerror: 0.063} + - {label: sys, symerror: 0.051} + - {label: 'sys,beam', symerror: 0.07} + - {label: 'sys,lumi', symerror: 0.071} + value: 6.059 + - errors: + - {label: stat, symerror: 0.054} + - {label: sys, symerror: 0.041} + - {label: 'sys,beam', symerror: 0.05} + - {label: 'sys,lumi', symerror: 0.051} + value: 4.385 + - errors: + - {label: stat, symerror: 0.042} + - {label: sys, symerror: 0.027} + - {label: 'sys,beam', symerror: 0.031} + - {label: 'sys,lumi', symerror: 0.032} + value: 2.724 + - errors: + - {label: stat, symerror: 0.032} + - {label: sys, symerror: 0.02} + - {label: 'sys,beam', symerror: 0.018} + - {label: 'sys,lumi', symerror: 0.019} + value: 1.584 + - errors: + - {label: stat, symerror: 0.022} + - {label: sys, symerror: 0.012} + - {label: 'sys,beam', symerror: 0.009} + - {label: 'sys,lumi', symerror: 0.009} + value: 0.749 + - errors: + - {label: stat, symerror: 0.016} + - {label: sys, symerror: 0.008} + - {label: 'sys,beam', symerror: 0.004} + - {label: 'sys,lumi', symerror: 0.004} + value: 0.383 + - errors: + - {label: stat, symerror: 0.003} + - {label: sys, symerror: 0.001} + - {label: 'sys,beam', symerror: 0.0} + - {label: 'sys,lumi', symerror: 0.0} + value: 0.011 +- header: {name: '$f^{Z}_{FSR}$'} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: M(MU+ MU-), units: GeV, value: 60.0-120.0} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: P P --> Z0 < MU+ MU- > X} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {symerror: 0.0395} + value: 1.0466 + - errors: + - {symerror: 0.0119} + value: 1.0305 + - errors: + - {symerror: 0.0069} + value: 1.0277 + - errors: + - {symerror: 0.0061} + value: 1.0252 + - errors: + - {symerror: 0.0048} + value: 1.0264 + - errors: + - {symerror: 0.0032} + value: 1.0257 + - errors: + - {symerror: 0.0038} + value: 1.0258 + - errors: + - {symerror: 0.0027} + value: 1.0252 + - errors: + - {symerror: 0.0035} + value: 1.0282 + - errors: + - {symerror: 0.003} + value: 1.0264 + - errors: + - {symerror: 0.0066} + value: 1.0261 + - errors: + - {symerror: 0.004} + value: 1.0248 + - errors: + - {symerror: 0.006} + value: 1.0258 + - errors: + - {symerror: 0.0053} + value: 1.0228 + - errors: + - {symerror: 0.0079} + value: 1.018 + - errors: + - {symerror: 0.01} + value: 1.0207 + - errors: + - {symerror: 0.014} + value: 1.0183 + - errors: + - {symerror: 0.0761} + value: 1.0177 +independent_variables: +- header: {name: '$y_{Z}$'} + values: + - {high: 2.125, low: 2.0} + - {high: 2.25, low: 2.125} + - {high: 2.375, low: 2.25} + - {high: 2.5, low: 2.375} + - {high: 2.625, low: 2.5} + - {high: 2.75, low: 2.625} + - {high: 2.875, low: 2.75} + - {high: 3.0, low: 2.875} + - {high: 3.125, low: 3.0} + - {high: 3.25, low: 3.125} + - {high: 3.375, low: 3.25} + - {high: 3.5, low: 3.375} + - {high: 3.625, low: 3.5} + - {high: 3.75, low: 3.625} + - {high: 3.875, low: 3.75} + - {high: 4.0, low: 3.875} + - {high: 4.25, low: 4.0} + - {high: 4.5, low: 4.25} diff --git a/buildmaster/LHCB_WPWM_8TEV_MUON/rawdata/corrmat.corr b/buildmaster/LHCB_WPWM_8TEV_MUON/rawdata/corrmat.corr new file mode 100644 index 0000000000..30782754d8 --- /dev/null +++ b/buildmaster/LHCB_WPWM_8TEV_MUON/rawdata/corrmat.corr @@ -0,0 +1,34 @@ +1.0 +0.19 1.0 +0.17 0.27 1.0 +0.16 0.26 0.28 1.0 +0.16 0.25 0.28 0.29 1.0 +0.15 0.24 0.27 0.29 0.3 1.0 +0.14 0.23 0.26 0.28 0.29 0.3 1.0 +0.14 0.21 0.25 0.27 0.29 0.3 0.3 1.0 +0.13 0.2 0.23 0.25 0.27 0.28 0.29 0.29 1.0 +0.11 0.17 0.2 0.23 0.25 0.26 0.27 0.28 0.27 1.0 +0.09 0.14 0.16 0.18 0.2 0.22 0.22 0.23 0.23 0.23 1.0 +0.08 0.12 0.15 0.17 0.19 0.2 0.21 0.22 0.22 0.22 0.2 1.0 +0.07 0.1 0.12 0.14 0.16 0.17 0.18 0.19 0.19 0.2 0.18 0.19 1.0 +0.06 0.08 0.1 0.11 0.13 0.14 0.15 0.16 0.16 0.17 0.16 0.16 0.15 1.0 +0.05 0.07 0.08 0.09 0.1 0.11 0.11 0.12 0.13 0.13 0.12 0.13 0.12 0.11 1.0 +0.03 0.05 0.06 0.06 0.07 0.08 0.08 0.09 0.09 0.1 0.09 0.1 0.09 0.08 0.07 1.0 +0.03 0.04 0.04 0.05 0.05 0.06 0.06 0.06 0.07 0.07 0.07 0.08 0.07 0.07 0.06 0.05 1.0 +0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.02 0.02 0.01 0.01 0.01 0.01 1.0 +0.23 0.3 0.28 0.27 0.26 0.25 0.24 0.23 0.21 0.18 0.15 0.13 0.11 0.1 0.07 0.05 0.04 0.01 1.0 +0.05 0.15 0.21 0.2 0.2 0.2 0.2 0.19 0.18 0.16 0.14 0.12 0.1 0.08 0.06 0.04 0.03 0.01 0.2 1.0 +0.04 0.07 0.12 0.15 0.16 0.17 0.17 0.17 0.16 0.15 0.13 0.13 0.11 0.08 0.06 0.05 0.03 0.01 0.13 0.12 1.0 +0.05 0.08 0.1 0.13 0.16 0.17 0.17 0.17 0.16 0.16 0.14 0.13 0.12 0.09 0.07 0.05 0.03 0.01 0.06 0.03 0.26 1.0 +0.06 0.08 0.1 0.11 0.14 0.16 0.17 0.17 0.16 0.16 0.14 0.14 0.12 0.1 0.07 0.05 0.03 0.01 0.07 0.03 0.25 0.33 1.0 +0.04 0.06 0.07 0.09 0.1 0.11 0.12 0.13 0.13 0.12 0.11 0.11 0.09 0.08 0.06 0.04 0.03 0.0 0.03 -0.01 0.28 0.35 0.34 1.0 +0.04 0.06 0.07 0.08 0.09 0.1 0.11 0.12 0.12 0.12 0.11 0.11 0.1 0.09 0.07 0.05 0.04 0.0 0.0 -0.06 0.31 0.41 0.4 0.45 1.0 +0.02 0.03 0.04 0.04 0.04 0.04 0.05 0.05 0.06 0.06 0.06 0.07 0.06 0.06 0.05 0.04 0.04 0.01 -0.07 -0.14 0.14 0.32 0.32 0.35 0.40 1.0 +0.21 0.28 0.26 0.25 0.24 0.24 0.22 0.21 0.2 0.17 0.14 0.12 0.11 0.09 0.07 0.05 0.04 0.01 0.67 0.10 0.24 0.22 0.22 0.23 0.26 0.14 1.0 +0.04 0.14 0.19 0.18 0.18 0.18 0.18 0.17 0.16 0.15 0.12 0.11 0.09 0.07 0.05 0.04 0.03 0.0 0.07 0.54 0.23 0.28 0.28 0.28 0.33 0.24 0.21 1.0 +0.04 0.07 0.11 0.14 0.15 0.15 0.15 0.15 0.15 0.14 0.12 0.12 0.1 0.08 0.06 0.04 0.03 0.01 0.05 0.03 0.64 0.27 0.26 0.27 0.32 0.23 0.18 0.22 1.0 +0.05 0.07 0.09 0.12 0.14 0.15 0.15 0.16 0.15 0.14 0.12 0.12 0.11 0.08 0.06 0.04 0.03 0.01 0.04 0.00 0.25 0.70 0.30 0.33 0.39 0.29 0.21 0.25 0.31 1.0 +0.05 0.08 0.09 0.1 0.13 0.15 0.15 0.16 0.15 0.15 0.13 0.13 0.11 0.09 0.07 0.05 0.03 0.01 0.06 0.03 0.26 0.32 0.68 0.32 0.38 0.26 0.22 0.28 0.27 0.32 1.0 +0.04 0.06 0.08 0.09 0.1 0.12 0.13 0.13 0.13 0.13 0.11 0.11 0.1 0.08 0.06 0.04 0.03 0.0 0.07 0.04 0.30 0.28 0.27 0.63 0.36 0.22 0.23 0.23 0.29 0.32 0.28 1.0 +0.04 0.06 0.08 0.09 0.1 0.11 0.12 0.13 0.14 0.14 0.12 0.12 0.11 0.1 0.08 0.06 0.04 0.01 0.14 0.20 -0.13 -0.08 -0.07 -0.17 -0.04 -0.15 -0.06 -0.04 -0.04 -0.11 -0.07 -0.19 1.0 +0.03 0.04 0.04 0.05 0.05 0.06 0.06 0.06 0.07 0.07 0.07 0.08 0.08 0.07 0.06 0.05 0.04 0.01 0.12 0.17 -0.18 -0.17 -0.15 -0.23 -0.31 0.05 -0.09 -0.11 -0.14 -0.19 -0.18 -0.24 0.48 1.0 diff --git a/buildmaster/LHCB_WPWM_8TEV_MUON/uncertainties.yaml b/buildmaster/LHCB_WPWM_8TEV_MUON/uncertainties.yaml new file mode 100644 index 0000000000..8be0ba9477 --- /dev/null +++ b/buildmaster/LHCB_WPWM_8TEV_MUON/uncertainties.yaml @@ -0,0 +1,670 @@ +definitions: + sys_corr_1: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_2: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_3: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_4: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_5: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_6: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_7: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_8: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_9: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_10: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_11: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_12: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_13: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_14: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_15: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_16: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + stat: + description: Uncorrelated statistical uncertainties + treatment: ADD + type: UNCORR + sys_beam: + description: Systematic Beam uncertainties + treatment: ADD + type: LHCBBEAM7TEV + sys_luminosity: + description: Systematic Luminosity uncertainties + treatment: ADD + type: LHCBLUMI7TEV +bins: +- sys_corr_1: -3038.5793019795487 + sys_corr_2: 768.151619729029 + sys_corr_3: -485.10950751290466 + sys_corr_4: 53.08972426745508 + sys_corr_5: 7.850087868619188 + sys_corr_6: -88.07520308797234 + sys_corr_7: 137.25556071652076 + sys_corr_8: 283.50818110253323 + sys_corr_9: -239.42618180677908 + sys_corr_10: -114.9282210891482 + sys_corr_11: 8.42245032181839 + sys_corr_12: 31.98627693875795 + sys_corr_13: 1.50717484391528 + sys_corr_14: 6.896281410681029 + sys_corr_15: -7.595533001147195 + sys_corr_16: -7.633120097776864 + sys_corr_17: 1.0195549272803348 + sys_corr_18: 7.284609808103727e-07 + sys_corr_19: -1.6778355582638205e-05 + sys_corr_20: 0.00025537746071577504 + sys_corr_21: -0.0011339103675613135 + sys_corr_22: -0.003415602574189624 + sys_corr_23: -0.005932843196956676 + sys_corr_24: -0.1280796093801708 + sys_corr_25: -7.564443428452589e-06 + sys_corr_26: 0.016222470669920764 + sys_corr_27: 0.06514473289220417 + sys_corr_28: 0.10371048272556845 + sys_corr_29: 0.06929157468645339 + sys_corr_30: 0.018776298638962904 + sys_corr_31: -0.021881712529226115 + sys_corr_32: -0.036065189677757854 + sys_corr_33: 0.03731793154064317 + sys_corr_34: -0.006212574243613175 + stat: 1200.0 + sys_beam: 2400.0 + sys_luminosity: 2700.0 +- sys_corr_1: -757.4168560908336 + sys_corr_2: 75.90309841401225 + sys_corr_3: 1987.0588271704808 + sys_corr_4: 380.50217714972797 + sys_corr_5: 186.6540628119187 + sys_corr_6: -98.92177817595513 + sys_corr_7: 19.386151565620303 + sys_corr_8: 90.24316654080108 + sys_corr_9: 97.86762228902114 + sys_corr_10: 257.88630047369355 + sys_corr_11: 27.80542496923231 + sys_corr_12: 182.1786676656212 + sys_corr_13: -40.1069768929091 + sys_corr_14: 11.34420533632285 + sys_corr_15: -37.39982404709859 + sys_corr_16: -32.34179686209255 + sys_corr_17: 1.0943227332297263 + sys_corr_18: -6.383721004848815e-06 + sys_corr_19: -3.428582414792927e-05 + sys_corr_20: 0.00046593140222169854 + sys_corr_21: 0.0016542604394633518 + sys_corr_22: -0.0033912447442104652 + sys_corr_23: 0.004148891256763143 + sys_corr_24: -0.009165474667733913 + sys_corr_25: 0.008144508683747427 + sys_corr_26: -0.026150151611025615 + sys_corr_27: -0.061315511041762666 + sys_corr_28: 0.012687205157747163 + sys_corr_29: 0.018693498085528087 + sys_corr_30: 0.01425778424606732 + sys_corr_31: -0.006241271517890433 + sys_corr_32: -0.006706266150301394 + sys_corr_33: 0.025389499587724858 + sys_corr_34: 0.00543093342266935 + stat: 900.0 + sys_beam: 2100.0 + sys_luminosity: 2400.0 +- sys_corr_1: -539.1116598566842 + sys_corr_2: -1082.196953898048 + sys_corr_3: 139.26538263419215 + sys_corr_4: 259.0007882801631 + sys_corr_5: -1235.3638860478686 + sys_corr_6: -147.8992484251472 + sys_corr_7: -30.971965139079433 + sys_corr_8: 66.12783421119019 + sys_corr_9: 79.27645225016443 + sys_corr_10: 6.1083395365405675 + sys_corr_11: 298.8256052265815 + sys_corr_12: -160.33052814602755 + sys_corr_13: 123.74448329511169 + sys_corr_14: -32.43971184140428 + sys_corr_15: -11.811037446007145 + sys_corr_16: -23.921899411063976 + sys_corr_17: 0.5767811158596198 + sys_corr_18: -2.8777510060249715e-06 + sys_corr_19: -1.7746259062853234e-05 + sys_corr_20: -0.0015668674597882877 + sys_corr_21: 0.00047893537625031795 + sys_corr_22: 0.00013796288062932144 + sys_corr_23: 0.015193558356304633 + sys_corr_24: 0.08368951181574025 + sys_corr_25: 0.03252795209784002 + sys_corr_26: 0.012074739801121198 + sys_corr_27: -0.07568283304928537 + sys_corr_28: -0.10839682755241084 + sys_corr_29: -0.06354115477152299 + sys_corr_30: 0.022137948768943465 + sys_corr_31: -0.00040771647319118586 + sys_corr_32: 0.029190918104867233 + sys_corr_33: -0.045338286183720544 + sys_corr_34: 0.009281120599104669 + stat: 800.0 + sys_beam: 1800.0 + sys_luminosity: 2100.0 +- sys_corr_1: -346.6272090627615 + sys_corr_2: -1030.181994954958 + sys_corr_3: 2.3905943298608827 + sys_corr_4: -672.7845398398682 + sys_corr_5: 431.24137151214705 + sys_corr_6: -736.2297158517016 + sys_corr_7: -0.9211964079113855 + sys_corr_8: 167.18955777019715 + sys_corr_9: 125.27532107231038 + sys_corr_10: 0.6952330910211038 + sys_corr_11: 153.6669942345208 + sys_corr_12: -132.51248962046776 + sys_corr_13: -319.0057886061013 + sys_corr_14: 92.20109904875524 + sys_corr_15: 52.134955258314406 + sys_corr_16: 11.019143373701814 + sys_corr_17: 0.8035534599422568 + sys_corr_18: -1.9494090664902555e-06 + sys_corr_19: 0.0004555292374753845 + sys_corr_20: -0.0005602017511280195 + sys_corr_21: 0.0017116083298735739 + sys_corr_22: -0.0012111825214759626 + sys_corr_23: 0.013093740239201866 + sys_corr_24: 0.0663328708802538 + sys_corr_25: 0.005741342860060582 + sys_corr_26: -0.005243893280854078 + sys_corr_27: -0.014898575132749711 + sys_corr_28: -0.07832385774732964 + sys_corr_29: -0.008103685999407343 + sys_corr_30: -0.028478571390293136 + sys_corr_31: -0.02209874007861173 + sys_corr_32: 0.06915915493569391 + sys_corr_33: -0.05225712733244364 + sys_corr_34: -0.021446400930781855 + stat: 700.0 + sys_beam: 1500.0 + sys_luminosity: 1800.0 +- sys_corr_1: -274.8755548245817 + sys_corr_2: -712.8587225878063 + sys_corr_3: 9.424364874672388 + sys_corr_4: -310.81026743552223 + sys_corr_5: 162.63611732393898 + sys_corr_6: 728.7752580649696 + sys_corr_7: -76.45444548329307 + sys_corr_8: 508.0244085137416 + sys_corr_9: 158.27207582811076 + sys_corr_10: 47.22082110864786 + sys_corr_11: 34.177741491179496 + sys_corr_12: 13.152729580561742 + sys_corr_13: -117.7021162614959 + sys_corr_14: -377.1740405636019 + sys_corr_15: 48.51723655535921 + sys_corr_16: 31.56153128981498 + sys_corr_17: 0.9822152491275007 + sys_corr_18: -2.4619998369348916e-06 + sys_corr_19: 0.0005612970734051026 + sys_corr_20: 0.0005386362464870477 + sys_corr_21: -0.001815465250724656 + sys_corr_22: -0.006347847483481793 + sys_corr_23: 0.017227109060982387 + sys_corr_24: 0.09711895842625598 + sys_corr_25: 0.03514572275551089 + sys_corr_26: 0.04299205294753696 + sys_corr_27: -0.011167609711057146 + sys_corr_28: -0.04891378723110805 + sys_corr_29: -0.09846940785881475 + sys_corr_30: 0.02063997187136929 + sys_corr_31: -0.03877240980632347 + sys_corr_32: 0.06848351774573039 + sys_corr_33: -0.04952056140314971 + sys_corr_34: -0.003110697817064795 + stat: 600.0 + sys_beam: 1200.0 + sys_luminosity: 1400.0 +- sys_corr_1: -167.11074525470343 + sys_corr_2: -569.573314600226 + sys_corr_3: -47.732461471837524 + sys_corr_4: -61.80663852923446 + sys_corr_5: 100.02781413206776 + sys_corr_6: 138.327542900984 + sys_corr_7: -68.1364178650373 + sys_corr_8: -175.18878405870703 + sys_corr_9: -517.694932637838 + sys_corr_10: 375.78574105417636 + sys_corr_11: -22.358913549372588 + sys_corr_12: -171.54256789293123 + sys_corr_13: -77.7987263976552 + sys_corr_14: -20.62483812309415 + sys_corr_15: -333.83570107460514 + sys_corr_16: 141.46991248825637 + sys_corr_17: 0.48158322751296384 + sys_corr_18: 4.342515998397655e-06 + sys_corr_19: 0.00012465687701839452 + sys_corr_20: 0.0002705298702991312 + sys_corr_21: 0.0015800747762540462 + sys_corr_22: -0.006142676190768974 + sys_corr_23: -0.010540544502378445 + sys_corr_24: 0.09906325673279748 + sys_corr_25: 0.02814902488221871 + sys_corr_26: 0.004982507679248855 + sys_corr_27: 0.029057225727090522 + sys_corr_28: -0.06471540581684256 + sys_corr_29: -0.04520273973785426 + sys_corr_30: 0.05994446708871603 + sys_corr_31: 0.05655868692291631 + sys_corr_32: -0.01759802246321701 + sys_corr_33: 0.034334691922246696 + sys_corr_34: -0.009203113866061742 + stat: 500.0 + sys_beam: 800.0 + sys_luminosity: 1000.0 +- sys_corr_1: -188.20228910192097 + sys_corr_2: -777.7968807961512 + sys_corr_3: -83.523149571827 + sys_corr_4: -219.69021728485177 + sys_corr_5: 118.43402492877796 + sys_corr_6: 252.7137844343597 + sys_corr_7: 294.355545808125 + sys_corr_8: -370.94543280886097 + sys_corr_9: -209.10455100615636 + sys_corr_10: -130.42393782916918 + sys_corr_11: 365.6282260301084 + sys_corr_12: 480.6203850587949 + sys_corr_13: 20.079947758739483 + sys_corr_14: 67.06685990860872 + sys_corr_15: 65.32932344673472 + sys_corr_16: 87.47120757016614 + sys_corr_17: -0.09826105456667454 + sys_corr_18: 8.386171914222564e-06 + sys_corr_19: -0.0010224700247338724 + sys_corr_20: -0.0019943997996906787 + sys_corr_21: 0.0084136973091293 + sys_corr_22: -0.02061594320152601 + sys_corr_23: 0.03529448123666306 + sys_corr_24: 0.18401404863525153 + sys_corr_25: 0.059450908402464545 + sys_corr_26: 0.005920112777482507 + sys_corr_27: 0.038432895713877376 + sys_corr_28: 0.01778891141365446 + sys_corr_29: -0.04235602286263034 + sys_corr_30: 0.019926253823197022 + sys_corr_31: 0.017823211055328757 + sys_corr_32: -0.028407360893841153 + sys_corr_33: -0.015168614203903776 + sys_corr_34: -0.00999576404121969 + stat: 500.0 + sys_beam: 900.0 + sys_luminosity: 1000.0 +- sys_corr_1: -20.418282214869066 + sys_corr_2: -294.128388854442 + sys_corr_3: -58.60919723330269 + sys_corr_4: -198.41896719160766 + sys_corr_5: 53.836537785331174 + sys_corr_6: 94.56972392250604 + sys_corr_7: -235.80701606356234 + sys_corr_8: -98.75510182567812 + sys_corr_9: -143.68649816865596 + sys_corr_10: -129.97796340553987 + sys_corr_11: 34.287470713688585 + sys_corr_12: 106.63233272961433 + sys_corr_13: -66.03659873895782 + sys_corr_14: -27.168441642658422 + sys_corr_15: -173.44340111205508 + sys_corr_16: -444.7463117440774 + sys_corr_17: -0.24157992028818057 + sys_corr_18: -1.9035352315276237e-05 + sys_corr_19: -0.0020680211971272366 + sys_corr_20: -0.002768681927241517 + sys_corr_21: 0.008251447625955019 + sys_corr_22: -0.017362115654713663 + sys_corr_23: 0.0005594231669970876 + sys_corr_24: 0.00048017663069296006 + sys_corr_25: 0.041775962924260265 + sys_corr_26: -0.0019150653015100835 + sys_corr_27: -0.03489947235145362 + sys_corr_28: 0.16155295845034465 + sys_corr_29: 0.045016400588073914 + sys_corr_30: -0.031051031050705183 + sys_corr_31: 0.007999837636393914 + sys_corr_32: -0.1345409740234758 + sys_corr_33: 0.08854643213947719 + sys_corr_34: -0.03573913245167926 + stat: 400.0 + sys_beam: 200.0 + sys_luminosity: 300.0 +- sys_corr_1: -1418.5660286961804 + sys_corr_2: -300.83618976028646 + sys_corr_3: -372.17857782986243 + sys_corr_4: 110.85229077349138 + sys_corr_5: 213.07294567190755 + sys_corr_6: 192.03560835810046 + sys_corr_7: -315.0880721533769 + sys_corr_8: -688.389231807664 + sys_corr_9: 537.9263076607414 + sys_corr_10: 188.2713964977463 + sys_corr_11: -46.3728491490399 + sys_corr_12: -58.69340052788875 + sys_corr_13: -14.504876419343034 + sys_corr_14: -14.744707792059215 + sys_corr_15: 3.8807055752414743 + sys_corr_16: 8.66452881632407 + sys_corr_17: 0.8337970308503326 + sys_corr_18: -3.098814478249894e-06 + sys_corr_19: -8.957681649067556e-05 + sys_corr_20: -0.0002853717212049436 + sys_corr_21: 0.0010659257978352338 + sys_corr_22: 0.0037092495827079847 + sys_corr_23: 0.006878015330544949 + sys_corr_24: -0.19054727226453752 + sys_corr_25: -0.024486657565609662 + sys_corr_26: 0.0053534274545727146 + sys_corr_27: 0.08616988285704642 + sys_corr_28: 0.14996497858208171 + sys_corr_29: 0.11390900690940575 + sys_corr_30: -0.03560935114680982 + sys_corr_31: 0.012681825895916057 + sys_corr_32: -0.03089693251724615 + sys_corr_33: 0.05081019719485488 + sys_corr_34: 0.052787600581333 + stat: 900.0 + sys_beam: 1200.0 + sys_luminosity: 1600.0 +- sys_corr_1: -379.7758627502127 + sys_corr_2: -614.702742631885 + sys_corr_3: 806.9723200388605 + sys_corr_4: 129.34897826869812 + sys_corr_5: 329.5697346931343 + sys_corr_6: 180.1001765648686 + sys_corr_7: -83.25772608032091 + sys_corr_8: -249.97342476735562 + sys_corr_9: -217.19952916877614 + sys_corr_10: -625.126074096715 + sys_corr_11: -89.05430339410013 + sys_corr_12: -320.1639682551919 + sys_corr_13: 70.13671643344883 + sys_corr_14: -28.02170633911008 + sys_corr_15: 56.046822735985636 + sys_corr_16: 47.589118060404786 + sys_corr_17: 0.4556251546296614 + sys_corr_18: 1.1495839610964414e-05 + sys_corr_19: 0.00013258230364422347 + sys_corr_20: 0.0001381739901380239 + sys_corr_21: -0.005874689172628592 + sys_corr_22: 0.010501500264532582 + sys_corr_23: -0.01794623114827708 + sys_corr_24: -0.21644468180010468 + sys_corr_25: -0.03488192226887612 + sys_corr_26: -0.06857859750272072 + sys_corr_27: -0.049221237602540484 + sys_corr_28: 0.07731091052149563 + sys_corr_29: 0.07604399926457621 + sys_corr_30: -0.025643183602270476 + sys_corr_31: -0.027707930572605736 + sys_corr_32: -0.039392965654304765 + sys_corr_33: 0.008744923900361117 + sys_corr_34: 0.005484652230440729 + stat: 700.0 + sys_beam: 1000.0 + sys_luminosity: 1400.0 +- sys_corr_1: -239.61349262110903 + sys_corr_2: -678.8206626587207 + sys_corr_3: 43.04893366922258 + sys_corr_4: -55.83800780748894 + sys_corr_5: -634.499154269975 + sys_corr_6: -25.947554268123866 + sys_corr_7: 41.427904422469275 + sys_corr_8: -54.053308811602285 + sys_corr_9: -16.362379243606593 + sys_corr_10: -129.38273499463483 + sys_corr_11: -592.7873743872755 + sys_corr_12: 256.3760680903285 + sys_corr_13: -265.5698980460149 + sys_corr_14: 43.29483884631226 + sys_corr_15: -1.6211538751018917 + sys_corr_16: 48.76481505905378 + sys_corr_17: 0.809829264135287 + sys_corr_18: -2.8081644487259427e-06 + sys_corr_19: 0.00010808868172066936 + sys_corr_20: 0.001635589344699784 + sys_corr_21: -0.00017695688461365108 + sys_corr_22: -0.0018627281863186546 + sys_corr_23: -0.007389095578738441 + sys_corr_24: -0.09727117032053585 + sys_corr_25: -0.004215152100620265 + sys_corr_26: -0.00531890952950511 + sys_corr_27: -0.013228979704683317 + sys_corr_28: 0.002784867494607145 + sys_corr_29: 0.07434987123378306 + sys_corr_30: -0.033407111180509784 + sys_corr_31: 0.020745837975215568 + sys_corr_32: -0.03494504746873362 + sys_corr_33: 0.04376355950975352 + sys_corr_34: -0.02172061471681362 + stat: 600.0 + sys_beam: 1000.0 + sys_luminosity: 1300.0 +- sys_corr_1: -222.80016800556484 + sys_corr_2: -742.5368422506577 + sys_corr_3: -36.8272748602642 + sys_corr_4: -408.192048198211 + sys_corr_5: 263.337600274069 + sys_corr_6: -456.734990608118 + sys_corr_7: -7.855868883315059 + sys_corr_8: 78.7247284309053 + sys_corr_9: 44.68970834520459 + sys_corr_10: 42.11339132489611 + sys_corr_11: -238.17127568567935 + sys_corr_12: 147.6273626390351 + sys_corr_13: 522.6255945867229 + sys_corr_14: -164.03665899925022 + sys_corr_15: -66.06062046337215 + sys_corr_16: 10.446280909012556 + sys_corr_17: 0.7081767021411348 + sys_corr_18: -4.695676677119722e-06 + sys_corr_19: -0.0001868136902926309 + sys_corr_20: 0.000835749485750196 + sys_corr_21: -0.0014647779094920056 + sys_corr_22: 0.0005052468688517873 + sys_corr_23: 0.01511537237976092 + sys_corr_24: -0.004813289205308594 + sys_corr_25: 0.024600067592098613 + sys_corr_26: 0.02837806872979211 + sys_corr_27: -0.002254224212950835 + sys_corr_28: -0.053259457136311486 + sys_corr_29: -0.04791596442701525 + sys_corr_30: 0.0681923382744364 + sys_corr_31: 0.04485150499439751 + sys_corr_32: 2.8937518308521544e-05 + sys_corr_33: 0.00247769781776522 + sys_corr_34: 0.01624354351478134 + stat: 600.0 + sys_beam: 900.0 + sys_luminosity: 1200.0 +- sys_corr_1: -221.03406480331822 + sys_corr_2: -598.7026972313603 + sys_corr_3: 7.879832997776385 + sys_corr_4: -218.47633151977615 + sys_corr_5: 119.79805532288195 + sys_corr_6: 540.0535716310895 + sys_corr_7: -14.335697968361968 + sys_corr_8: 370.33754643475316 + sys_corr_9: 156.56323982872624 + sys_corr_10: 48.85089017217374 + sys_corr_11: -72.25805431862382 + sys_corr_12: -60.236085152793216 + sys_corr_13: 166.24096100027964 + sys_corr_14: 496.04191921051756 + sys_corr_15: -42.90701380970265 + sys_corr_16: -16.03726295711185 + sys_corr_17: 0.4785205369484593 + sys_corr_18: -4.710614020048994e-06 + sys_corr_19: -9.20596815872463e-05 + sys_corr_20: -0.0012894448180824792 + sys_corr_21: 0.005305379878287117 + sys_corr_22: -0.004502053808523869 + sys_corr_23: 0.009383386501251282 + sys_corr_24: 0.1241067869378707 + sys_corr_25: 0.03652483327628769 + sys_corr_26: -0.026124954404075508 + sys_corr_27: 0.027004757545574683 + sys_corr_28: -0.022177739726991944 + sys_corr_29: -0.04456292446841822 + sys_corr_30: 0.02875750635448535 + sys_corr_31: 0.05203761053449261 + sys_corr_32: 0.0586078660639764 + sys_corr_33: -0.05914759281192592 + sys_corr_34: 0.0033480429937501 + stat: 600.0 + sys_beam: 800.0 + sys_luminosity: 1100.0 +- sys_corr_1: -175.95234830648135 + sys_corr_2: -458.29719201587636 + sys_corr_3: -25.221580090315847 + sys_corr_4: 34.36989661893009 + sys_corr_5: 36.823678063795654 + sys_corr_6: 68.23139353147354 + sys_corr_7: -45.66847954320662 + sys_corr_8: -111.32792227020646 + sys_corr_9: -424.8650642158514 + sys_corr_10: 408.5006539031898 + sys_corr_11: -128.01638358132988 + sys_corr_12: -105.17863116654347 + sys_corr_13: 49.58952197883063 + sys_corr_14: 17.210940422189474 + sys_corr_15: 392.3340808181992 + sys_corr_16: -123.20759730685347 + sys_corr_17: 0.7580998432848071 + sys_corr_18: 5.007363639324624e-06 + sys_corr_19: -0.00036884888231146186 + sys_corr_20: -0.0012978386362287946 + sys_corr_21: 0.005410694702775865 + sys_corr_22: -0.013313942732979372 + sys_corr_23: 0.05115695988736416 + sys_corr_24: 0.23458858072676236 + sys_corr_25: 0.06360147645469767 + sys_corr_26: 0.013309433725921567 + sys_corr_27: -0.04071310626280849 + sys_corr_28: -0.028177852122879966 + sys_corr_29: -0.1799505617141715 + sys_corr_30: -0.02169242324782223 + sys_corr_31: -0.007009667327721166 + sys_corr_32: 0.03717354123857821 + sys_corr_33: -0.12043966335792472 + sys_corr_34: 0.01543371410455415 + stat: 500.0 + sys_beam: 700.0 + sys_luminosity: 900.0 +- sys_corr_1: -183.14713236431768 + sys_corr_2: 567.3613770804149 + sys_corr_3: 468.56312146099094 + sys_corr_4: -942.5696139985873 + sys_corr_5: -309.77329328945194 + sys_corr_6: 193.68459039286637 + sys_corr_7: 711.2139889472766 + sys_corr_8: -281.5128856875479 + sys_corr_9: 108.21391439053171 + sys_corr_10: 77.90499127010314 + sys_corr_11: -58.041274291234444 + sys_corr_12: -198.5639555860381 + sys_corr_13: 23.54786254952854 + sys_corr_14: -43.82890699205783 + sys_corr_15: -20.50657414285486 + sys_corr_16: -67.43375805284023 + sys_corr_17: 1.4190064875049369 + sys_corr_18: -3.892991343602838e-06 + sys_corr_19: -0.00017797671767751046 + sys_corr_20: -0.0017744229331906653 + sys_corr_21: 0.005908187636415936 + sys_corr_22: -0.014350998864706447 + sys_corr_23: 0.02762243269222977 + sys_corr_24: 0.18750360514629275 + sys_corr_25: 0.05462605794528981 + sys_corr_26: 0.010284578419506411 + sys_corr_27: -0.040815268044172244 + sys_corr_28: -0.0392171074296804 + sys_corr_29: -0.1371256172369289 + sys_corr_30: 0.017759500556972105 + sys_corr_31: 0.06695235138005146 + sys_corr_32: -0.00374669242826696 + sys_corr_33: -0.061967289524856445 + sys_corr_34: -0.019287334663180073 + stat: 600.0 + sys_beam: 1000.0 + sys_luminosity: 1400.0 +- sys_corr_1: -108.57118702106705 + sys_corr_2: 839.8291501072293 + sys_corr_3: 459.8390870183657 + sys_corr_4: -925.6968612905044 + sys_corr_5: -400.4040360456309 + sys_corr_6: 42.68917936257094 + sys_corr_7: -752.2174879814365 + sys_corr_8: -53.07108189474481 + sys_corr_9: -128.09208022343634 + sys_corr_10: -32.98357281980944 + sys_corr_11: 74.02432402617801 + sys_corr_12: 78.5267575038786 + sys_corr_13: 36.42357651871526 + sys_corr_14: 39.17038903485109 + sys_corr_15: 53.70701716655165 + sys_corr_16: 94.08338628291777 + sys_corr_17: 0.6888256080273181 + sys_corr_18: 4.136991764242596e-07 + sys_corr_19: -0.000857665038759362 + sys_corr_20: -0.0026926244135948014 + sys_corr_21: 0.0075394244035477806 + sys_corr_22: -0.01576189784794336 + sys_corr_23: 0.04972357343844988 + sys_corr_24: 0.14135531791323716 + sys_corr_25: 0.08199691279107951 + sys_corr_26: 0.010188403737656443 + sys_corr_27: 0.02352486740748513 + sys_corr_28: -0.04847116685701623 + sys_corr_29: -0.03684786732394954 + sys_corr_30: -0.0070126705255896095 + sys_corr_31: 0.01275766566589463 + sys_corr_32: 0.022736517156611544 + sys_corr_33: -0.04112496966862067 + sys_corr_34: 0.018462409568420584 + stat: 700.0 + sys_beam: 500.0 + sys_luminosity: 700.0 diff --git a/buildmaster/LHCB_Z0_8TEV_MUON/data.yaml b/buildmaster/LHCB_Z0_8TEV_MUON/data.yaml new file mode 100644 index 0000000000..8b8f5759cc --- /dev/null +++ b/buildmaster/LHCB_Z0_8TEV_MUON/data.yaml @@ -0,0 +1,19 @@ +data_central: +- 1223.0 +- 3263.0 +- 4983.0 +- 6719.0 +- 8051.0 +- 8967.0 +- 9561.0 +- 9822.0 +- 9721.0 +- 9030.0 +- 7748.0 +- 6059.0 +- 4385.0 +- 2724.0 +- 1584.0 +- 749.0 +- 383.0 +- 11.0 diff --git a/buildmaster/LHCB_Z0_8TEV_MUON/filter.py b/buildmaster/LHCB_Z0_8TEV_MUON/filter.py new file mode 100644 index 0000000000..37509a2e5c --- /dev/null +++ b/buildmaster/LHCB_Z0_8TEV_MUON/filter.py @@ -0,0 +1,390 @@ +import numpy as np +import pandas as pd +import pathlib +import yaml + +from validphys.commondata_utils import covmat_to_artunc + + +MZ_VALUE = 91.1876 # GeV +MW_VALUE = 80.398 # GeV +SQRT_S = 7_000.0 # GeV +NORM_FACTOR = 1_000.0 # from pb -> fb + +# Correct tables to read values [[Z], [W+, W-]] +TABLES = {2: [0], 1: [0, 2]} # {table_id: [indexes]} +NB_POINTS = [18, 8, 8] # [N_Z, N_W+, N_W-] + +# MAP Boson string to M values +MAP_BOSON = {"Z": MZ_VALUE, "W": MW_VALUE} +MAP_TABLE = {2: "Z", 1: "W"} + + +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-ins1406555-v{version}-Table_{table_id}" + table = pathlib.Path(f"./rawdata/{filename}.yaml") + + return yaml.safe_load(table.read_text()) + + +def read_metadata() -> tuple[int, int, list]: + """Read the version and list of tables from metadata. + + Returns + ------- + tuple(int, list): + data version and list of hepdata tables + + """ + metadata = pathlib.Path("./metadata.yaml") + content = yaml.safe_load(metadata.read_text()) + + version = content["hepdata"]["version"] + nb_datapoints = sum(content["implemented_observables"][0]["npoints"]) + tables = content["implemented_observables"][0]["tables"] + + return version, nb_datapoints, tables + + +def get_kinematics(hepdata: dict, bin_index: list, boson: str = "Z") -> 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 + + Returns + ------- + tuple(int, list): + data version and list of hepdata tables + + """ + rapbins = hepdata["independent_variables"][0]["values"] + + kinematics = [] + for bins in bin_index: + ymin = float(rapbins[bins]["low"]) + ymax = float(rapbins[bins]["high"]) + kin_value = { + "y": { + "min": ymin, + "mid": 0.5 * (ymin + ymax), + "max": ymax, + }, + "M2": {"min": None, "mid": MAP_BOSON[boson] ** 2, "max": None}, + "sqrt_s": {"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 + list of Non-empty bin index + idx: int + 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 [NORM_FACTOR * central[i]["value"] for i in bin_index] + + +def get_errors(hepdata: dict, bin_index: list, indx: int = 0) -> dict: + """Extract the error values from the HepData yaml file. + + Parameters + ---------- + hepdata: dict + dictionary containing all data info + bin_index: list + list of Non-empty bin index + indx: int + index from which the errors will be read + + Returns + ------- + list: + list of dictionaries whose contents are the various + source of uncertainties + + """ + errors = hepdata["dependent_variables"][indx]["values"] + + stat, sys_corr, sys_beam, sys_lumi = [], [], [], [] + for idx in bin_index: + stat.append(NORM_FACTOR * errors[idx]["errors"][0]["symerror"]) + sys_corr.append(NORM_FACTOR * errors[idx]["errors"][1]["symerror"]) + sys_beam.append(NORM_FACTOR * errors[idx]["errors"][2]["symerror"]) + sys_lumi.append(NORM_FACTOR * errors[idx]["errors"][3]["symerror"]) + + return { + "stat": stat, + "sys_corr": sys_corr, + "sys_beam": sys_beam, + "sys_lumi": sys_lumi, + } + + +def read_corrmatrix(nb_datapoints: int) -> np.ndarray: + """Read the matrix and returns a symmetrized verions. + + Parameters + ---------- + nb_datapoints: int + total number of datapoints + + Returns + ------- + np.ndarray: (nb_datapoints, nb_datapoints) + a symmetrized square matrix + + """ + corrmat = pd.read_csv( + "./rawdata/corrmat.corr", + names=[f'{i}' for i in range(nb_datapoints)], + delim_whitespace=True, + ) + corrmat = corrmat.iloc[:, :].values + + # Symmetrize the Correlation Matrix + symmetrised = np.triu(corrmat.T) + np.triu(corrmat.T, 1).T + assert np.allclose(symmetrised, symmetrised.T) + + return symmetrised + + +def multiply_syst(covmat: np.ndarray, sys: list) -> np.ndarray: + """Multiply the covmat with the total systematics. + + Parameters + ---------- + covmat: np.ndarray + a squared matrix representing the covmat + sys: list + a vector containing the correlatdd systematics + + Returns + ------- + np.ndarray: + covariance matrix multiplied by the total systematics + + """ + covmat_sys = np.zeros(covmat.shape) + + for i in range(covmat.shape[0]): + for j in range(covmat.shape[-1]): + covmat_sys[i][j] = covmat[i][j] * sys[i] * sys[j] + + return covmat_sys.flatten() + + +def generate_artificial_unc(**kwargs) -> np.ndarray: + """A wrapper around `covmat_to_artunc` to return an array. + + Returns + ------- + np.ndarray: + a squared matrix with artificial uncertainties + """ + artunc = covmat_to_artunc(**kwargs) + return np.array(artunc) + + +def concatenate_dicts(multidict: list[dict]) -> dict: + """Given a list of dictionaries combined the values. + + Parameters + ---------- + multidict: list[dict] + list of dictionaries with the same keys + + Returns + ------- + dict: + dictionary whose keys are combined + """ + new_dict = {} + for key in multidict[0].keys(): + new_dict[key] = [] + for element in multidict: + new_dict[key] += element[key] + + return new_dict + + +def format_uncertainties(uncs: dict, artunc: np.ndarray) -> list: + """Format the uncertainties to be dumped into the yaml file. + + Parameters + ---------- + uncs: dict + dictionary containing the various source of uncertainties + artunc: no.ndarray + array of the artificial systematic uncertainties + + Returns + ------- + list: + list of ditionaries whose elements are the various errors + + """ + combined_errors = [] + + artunc = artunc.tolist() + for idat in range(len(artunc)): + error_value = {} + for jdat in range(len(artunc[idat])): + error_value[f"sys_corr_{jdat + 1}"] = artunc[idat][jdat] + + error_value["stat"] = uncs["stat"][idat] + error_value["sys_beam"] = uncs["sys_beam"][idat] + error_value["sys_luminosity"] = uncs["sys_lumi"][idat] + combined_errors.append(error_value) + + return combined_errors + + +def dump_commondata(kinematics: list, data: list, errors: list) -> 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: list + list containing the different errors + + """ + error_definition = { + f"sys_corr_{i + 1}": { + "description": "Correlated systematic uncertainties", + "treatment": "ADD", + "type": "CORR", + } + for i in range(len(data)) + } + + error_definition["stat"] = { + "description": "Uncorrelated statistical uncertainties", + "treatment": "ADD", + "type": "UNCORR", + } + + error_definition["sys_beam"] = { + "description": "Systematic Beam uncertainties", + "treatment": "ADD", + "type": "LHCBBEAM7TEV", + } + + error_definition["sys_luminosity"] = { + "description": "Systematic Luminosity uncertainties", + "treatment": "ADD", + "type": "LHCBLUMI7TEV", + } + + with open("data.yaml", "w") as file: + yaml.dump({"data_central": data}, file, sort_keys=False) + + with open("kinematics.yaml", "w") as file: + yaml.dump({"bins": kinematics}, file, sort_keys=False) + + with open("uncertainties.yaml", "w") as file: + yaml.dump( + {"definitions": error_definition, "bins": errors}, + file, + sort_keys=False, + ) + + +def main_filter(boson: str = "Z") -> None: + """Main driver of the LHCB_ZEE_8TEV filter that produces commmondata. + + There are three main different sources of uncertainties. + + 1. Statistical uncertainties: ADD, UNCORR + + 2. Correlated Systematic uncertainties: ADD, CORR: + constructed by symmetrizing the correlation matrix and extracting + the artificial systematic uncertainties from the latter + + 3. Beam Systematic uncertainties: ADD, LHCBBEAM7TEV + + 4. Luminosity Systematic uncertainties: ADD, LHCBLUMI7TEV + + """ + version, _, _ = read_metadata() + nbpoints = sum(NB_POINTS) + + # Select only the bins corresponding to the correct boson + bslice = slice(18) if boson == "Z" else slice(18, nbpoints + 1) + + nbp_idx = 0 + comb_kins, comb_data = [], [] + combined_errors = [] + for tabid in TABLES.keys(): # Loop over tables + for idx in TABLES[tabid]: # Loop over Bosons [Z, W+, W-] + bin_index = [i for i in range(NB_POINTS[nbp_idx])] + yaml_content = load_yaml(table_id=tabid, version=version) + + # Extract the kinematic, data, and uncertainties + kinematics = get_kinematics(yaml_content, bin_index, MAP_TABLE[tabid]) + data_central = get_data_values(yaml_content, bin_index, indx=idx) + uncertainties = get_errors(yaml_content, bin_index, indx=idx) + + # Collect all the results from different tables + comb_kins += kinematics + comb_data += data_central + combined_errors.append(uncertainties) + + nbp_idx += 1 + + errors_combined = concatenate_dicts(combined_errors) + # Compute the Artifical Systematics from CovMat + corrmat = read_corrmatrix(nb_datapoints=nbpoints) + covmat = multiply_syst(corrmat, errors_combined["sys_corr"]) + artunc = generate_artificial_unc( + ndata=nbpoints, + covmat_list=covmat.tolist(), + no_of_norm_mat=0, + ) + errors = format_uncertainties(errors_combined, artunc) + + # Generate all the necessary files + dump_commondata(comb_kins[bslice], comb_data[bslice], errors[bslice]) + + return + + +if __name__ == "__main__": + main_filter(boson="Z") diff --git a/buildmaster/LHCB_Z0_8TEV_MUON/kinematics.yaml b/buildmaster/LHCB_Z0_8TEV_MUON/kinematics.yaml new file mode 100644 index 0000000000..634a97d47d --- /dev/null +++ b/buildmaster/LHCB_Z0_8TEV_MUON/kinematics.yaml @@ -0,0 +1,217 @@ +bins: +- y: + min: 2.0 + mid: 2.0625 + max: 2.125 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.125 + mid: 2.1875 + max: 2.25 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.25 + mid: 2.3125 + max: 2.375 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.375 + mid: 2.4375 + max: 2.5 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.5 + mid: 2.5625 + max: 2.625 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.625 + mid: 2.6875 + max: 2.75 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.75 + mid: 2.8125 + max: 2.875 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 2.875 + mid: 2.9375 + max: 3.0 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.0 + mid: 3.0625 + max: 3.125 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.125 + mid: 3.1875 + max: 3.25 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.25 + mid: 3.3125 + max: 3.375 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.375 + mid: 3.4375 + max: 3.5 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.5 + mid: 3.5625 + max: 3.625 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.625 + mid: 3.6875 + max: 3.75 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.75 + mid: 3.8125 + max: 3.875 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 3.875 + mid: 3.9375 + max: 4.0 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 4.0 + mid: 4.125 + max: 4.25 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null +- y: + min: 4.25 + mid: 4.375 + max: 4.5 + M2: + min: null + mid: 8315.178393760001 + max: null + sqrt_s: + min: null + mid: 7000.0 + max: null diff --git a/buildmaster/LHCB_Z0_8TEV_MUON/metadata.yaml b/buildmaster/LHCB_Z0_8TEV_MUON/metadata.yaml new file mode 100644 index 0000000000..f9e005d77d --- /dev/null +++ b/buildmaster/LHCB_Z0_8TEV_MUON/metadata.yaml @@ -0,0 +1,47 @@ +setname: "LHCB_Z0_8TEV_MUON" + +version: 1 +version_comment: "First implementation" + +# References +arXiv: + url: "https://arxiv.org/abs/1511.08039" +iNSPIRE: + url: "https://inspirehep.net/literature/1406555" +hepdata: + url: "https://www.hepdata.net/record/ins1406555" + version: 1 + +nnpdf_metadata: + nnpdf31_process: "DY NC" + experiment: "LHCb" + +implemented_observables: + +- observable_name: "y" + observable: + description: "Inclusive cross-section for Z boson production in bins of muon pseudorapidity at 8 TeV" + label: r"$d\sigma / d|y|$" + units: "[fb]" + process_type: "EWK_RAP" + tables: [2] + npoints: [18] + plotting: + kinematics_override: ewk_rap_sqrt_scale + dataset_label: "LHCb $Z \\to \\mu$ 8 TeV" + plot_x: y + y_label: '$d\sigma/dy$ (fb)' + kinematic_coverage: [y, M2, sqrt_s] + kinematics: + variables: + y: {description: "Muon rapidity", label: "$y$", units: ""} + M2: {description: "Z boson Mass", label: "$M^2$", units: "$GeV^2$"} + sqrt_s: {description: "Center of Mass Energy", label: '$\sqrt{s}$', units: "$GeV$"} + file: kinematics.yaml + data_central: data.yaml + data_uncertainties: + - uncertainties.yaml + theory: + FK_tables: + - - LHCBZMU8TEV + operation: 'null' diff --git a/buildmaster/LHCB_Z0_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_1.yaml b/buildmaster/LHCB_Z0_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_1.yaml new file mode 100644 index 0000000000..0332da55b0 --- /dev/null +++ b/buildmaster/LHCB_Z0_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_1.yaml @@ -0,0 +1,184 @@ +dependent_variables: +- header: {name: '$\sigma_{W^{+} \to \mu^{+}\nu}$', units: PB} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: P P --> W+ < MU+ NU > X} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {label: stat, symerror: 1.2} + - {label: sys, symerror: 3.2} + - {label: 'sys,beam', symerror: 2.4} + - {label: 'sys,lumi', symerror: 2.7} + value: 236.5 + - errors: + - {label: stat, symerror: 0.9} + - {label: sys, symerror: 2.2} + - {label: 'sys,beam', symerror: 2.1} + - {label: 'sys,lumi', symerror: 2.4} + value: 208.4 + - errors: + - {label: stat, symerror: 0.8} + - {label: sys, symerror: 1.8} + - {label: 'sys,beam', symerror: 1.8} + - {label: 'sys,lumi', symerror: 2.1} + value: 182.0 + - errors: + - {label: stat, symerror: 0.7} + - {label: sys, symerror: 1.6} + - {label: 'sys,beam', symerror: 1.5} + - {label: 'sys,lumi', symerror: 1.8} + value: 153.3 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.3} + - {label: 'sys,beam', symerror: 1.2} + - {label: 'sys,lumi', symerror: 1.4} + value: 119.5 + - errors: + - {label: stat, symerror: 0.5} + - {label: sys, symerror: 1.0} + - {label: 'sys,beam', symerror: 0.8} + - {label: 'sys,lumi', symerror: 1.0} + value: 84.4 + - errors: + - {label: stat, symerror: 0.5} + - {label: sys, symerror: 1.2} + - {label: 'sys,beam', symerror: 0.9} + - {label: 'sys,lumi', symerror: 1.0} + value: 86.4 + - errors: + - {label: stat, symerror: 0.4} + - {label: sys, symerror: 0.7} + - {label: 'sys,beam', symerror: 0.2} + - {label: 'sys,lumi', symerror: 0.3} + value: 23.0 +- header: {name: '$f^{W^+}_{FSR}$'} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: ''} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {symerror: 0.0047} + value: 1.0188 + - errors: + - {symerror: 0.0028} + value: 1.0163 + - errors: + - {symerror: 0.0025} + value: 1.0158 + - errors: + - {symerror: 0.0028} + value: 1.0148 + - errors: + - {symerror: 0.0032} + value: 1.0152 + - errors: + - {symerror: 0.0046} + value: 1.015 + - errors: + - {symerror: 0.0045} + value: 1.0175 + - errors: + - {symerror: 0.0087} + value: 1.0211 +- header: {name: '$\sigma_{W^{-} \to \mu^{-}\overline{\nu}}$', units: PB} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: P P --> W- < MU- NUBAR > X} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {label: stat, symerror: 0.9} + - {label: sys, symerror: 1.8} + - {label: 'sys,beam', symerror: 1.2} + - {label: 'sys,lumi', symerror: 1.6} + value: 134.0 + - errors: + - {label: stat, symerror: 0.7} + - {label: sys, symerror: 1.4} + - {label: 'sys,beam', symerror: 1.0} + - {label: 'sys,lumi', symerror: 1.4} + value: 119.8 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.2} + - {label: 'sys,beam', symerror: 1.0} + - {label: 'sys,lumi', symerror: 1.3} + value: 110.6 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.2} + - {label: 'sys,beam', symerror: 0.9} + - {label: 'sys,lumi', symerror: 1.2} + value: 102.4 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.1} + - {label: 'sys,beam', symerror: 0.8} + - {label: 'sys,lumi', symerror: 1.1} + value: 92.5 + - errors: + - {label: stat, symerror: 0.5} + - {label: sys, symerror: 0.9} + - {label: 'sys,beam', symerror: 0.7} + - {label: 'sys,lumi', symerror: 0.9} + value: 79.9 + - errors: + - {label: stat, symerror: 0.6} + - {label: sys, symerror: 1.5} + - {label: 'sys,beam', symerror: 1.0} + - {label: 'sys,lumi', symerror: 1.4} + value: 119.3 + - errors: + - {label: stat, symerror: 0.7} + - {label: sys, symerror: 1.6} + - {label: 'sys,beam', symerror: 0.5} + - {label: 'sys,lumi', symerror: 0.7} + value: 60.0 +- header: {name: '$f^{W^-}_{FSR}$'} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: ''} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {symerror: 0.0026} + value: 1.0172 + - errors: + - {symerror: 0.0027} + value: 1.0155 + - errors: + - {symerror: 0.0028} + value: 1.0153 + - errors: + - {symerror: 0.003} + value: 1.0162 + - errors: + - {symerror: 0.0031} + value: 1.016 + - errors: + - {symerror: 0.0033} + value: 1.0176 + - errors: + - {symerror: 0.0033} + value: 1.02 + - errors: + - {symerror: 0.0053} + value: 1.0243 +independent_variables: +- header: {name: '$\eta^{\mu}$'} + values: + - {high: 2.25, low: 2.0} + - {high: 2.5, low: 2.25} + - {high: 2.75, low: 2.5} + - {high: 3.0, low: 2.75} + - {high: 3.25, low: 3.0} + - {high: 3.5, low: 3.25} + - {high: 4.0, low: 3.5} + - {high: 4.5, low: 4.0} diff --git a/buildmaster/LHCB_Z0_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_2.yaml b/buildmaster/LHCB_Z0_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_2.yaml new file mode 100644 index 0000000000..9a39670c3b --- /dev/null +++ b/buildmaster/LHCB_Z0_8TEV_MUON/rawdata/HEPData-ins1406555-v1-Table_2.yaml @@ -0,0 +1,200 @@ +dependent_variables: +- header: {name: '$\sigma_{Z \to \mu^{+}\mu^{-}}$', units: PB} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: M(MU+ MU-), units: GeV, value: 60.0-120.0} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: P P --> Z0 < MU+ MU- > X} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {label: stat, symerror: 0.033} + - {label: sys, symerror: 0.055} + - {label: 'sys,beam', symerror: 0.014} + - {label: 'sys,lumi', symerror: 0.014} + value: 1.223 + - errors: + - {label: stat, symerror: 0.051} + - {label: sys, symerror: 0.06} + - {label: 'sys,beam', symerror: 0.038} + - {label: 'sys,lumi', symerror: 0.038} + value: 3.263 + - errors: + - {label: stat, symerror: 0.062} + - {label: sys, symerror: 0.064} + - {label: 'sys,beam', symerror: 0.057} + - {label: 'sys,lumi', symerror: 0.058} + value: 4.983 + - errors: + - {label: stat, symerror: 0.07} + - {label: sys, symerror: 0.072} + - {label: 'sys,beam', symerror: 0.077} + - {label: 'sys,lumi', symerror: 0.078} + value: 6.719 + - errors: + - {label: stat, symerror: 0.076} + - {label: sys, symerror: 0.074} + - {label: 'sys,beam', symerror: 0.093} + - {label: 'sys,lumi', symerror: 0.094} + value: 8.051 + - errors: + - {label: stat, symerror: 0.079} + - {label: sys, symerror: 0.074} + - {label: 'sys,beam', symerror: 0.103} + - {label: 'sys,lumi', symerror: 0.105} + value: 8.967 + - errors: + - {label: stat, symerror: 0.081} + - {label: sys, symerror: 0.076} + - {label: 'sys,beam', symerror: 0.11} + - {label: 'sys,lumi', symerror: 0.112} + value: 9.561 + - errors: + - {label: stat, symerror: 0.082} + - {label: sys, symerror: 0.071} + - {label: 'sys,beam', symerror: 0.113} + - {label: 'sys,lumi', symerror: 0.115} + value: 9.822 + - errors: + - {label: stat, symerror: 0.081} + - {label: sys, symerror: 0.074} + - {label: 'sys,beam', symerror: 0.112} + - {label: 'sys,lumi', symerror: 0.114} + value: 9.721 + - errors: + - {label: stat, symerror: 0.078} + - {label: sys, symerror: 0.071} + - {label: 'sys,beam', symerror: 0.104} + - {label: 'sys,lumi', symerror: 0.105} + value: 9.03 + - errors: + - {label: stat, symerror: 0.072} + - {label: sys, symerror: 0.074} + - {label: 'sys,beam', symerror: 0.089} + - {label: 'sys,lumi', symerror: 0.09} + value: 7.748 + - errors: + - {label: stat, symerror: 0.063} + - {label: sys, symerror: 0.051} + - {label: 'sys,beam', symerror: 0.07} + - {label: 'sys,lumi', symerror: 0.071} + value: 6.059 + - errors: + - {label: stat, symerror: 0.054} + - {label: sys, symerror: 0.041} + - {label: 'sys,beam', symerror: 0.05} + - {label: 'sys,lumi', symerror: 0.051} + value: 4.385 + - errors: + - {label: stat, symerror: 0.042} + - {label: sys, symerror: 0.027} + - {label: 'sys,beam', symerror: 0.031} + - {label: 'sys,lumi', symerror: 0.032} + value: 2.724 + - errors: + - {label: stat, symerror: 0.032} + - {label: sys, symerror: 0.02} + - {label: 'sys,beam', symerror: 0.018} + - {label: 'sys,lumi', symerror: 0.019} + value: 1.584 + - errors: + - {label: stat, symerror: 0.022} + - {label: sys, symerror: 0.012} + - {label: 'sys,beam', symerror: 0.009} + - {label: 'sys,lumi', symerror: 0.009} + value: 0.749 + - errors: + - {label: stat, symerror: 0.016} + - {label: sys, symerror: 0.008} + - {label: 'sys,beam', symerror: 0.004} + - {label: 'sys,lumi', symerror: 0.004} + value: 0.383 + - errors: + - {label: stat, symerror: 0.003} + - {label: sys, symerror: 0.001} + - {label: 'sys,beam', symerror: 0.0} + - {label: 'sys,lumi', symerror: 0.0} + value: 0.011 +- header: {name: '$f^{Z}_{FSR}$'} + qualifiers: + - {name: ETARAP(MU), value: 2.0-4.5} + - {name: M(MU+ MU-), units: GeV, value: 60.0-120.0} + - {name: PT(MU), units: GEV, value: '> 20'} + - {name: RE, value: P P --> Z0 < MU+ MU- > X} + - {name: SQRT(S), units: GeV, value: '8000.0'} + values: + - errors: + - {symerror: 0.0395} + value: 1.0466 + - errors: + - {symerror: 0.0119} + value: 1.0305 + - errors: + - {symerror: 0.0069} + value: 1.0277 + - errors: + - {symerror: 0.0061} + value: 1.0252 + - errors: + - {symerror: 0.0048} + value: 1.0264 + - errors: + - {symerror: 0.0032} + value: 1.0257 + - errors: + - {symerror: 0.0038} + value: 1.0258 + - errors: + - {symerror: 0.0027} + value: 1.0252 + - errors: + - {symerror: 0.0035} + value: 1.0282 + - errors: + - {symerror: 0.003} + value: 1.0264 + - errors: + - {symerror: 0.0066} + value: 1.0261 + - errors: + - {symerror: 0.004} + value: 1.0248 + - errors: + - {symerror: 0.006} + value: 1.0258 + - errors: + - {symerror: 0.0053} + value: 1.0228 + - errors: + - {symerror: 0.0079} + value: 1.018 + - errors: + - {symerror: 0.01} + value: 1.0207 + - errors: + - {symerror: 0.014} + value: 1.0183 + - errors: + - {symerror: 0.0761} + value: 1.0177 +independent_variables: +- header: {name: '$y_{Z}$'} + values: + - {high: 2.125, low: 2.0} + - {high: 2.25, low: 2.125} + - {high: 2.375, low: 2.25} + - {high: 2.5, low: 2.375} + - {high: 2.625, low: 2.5} + - {high: 2.75, low: 2.625} + - {high: 2.875, low: 2.75} + - {high: 3.0, low: 2.875} + - {high: 3.125, low: 3.0} + - {high: 3.25, low: 3.125} + - {high: 3.375, low: 3.25} + - {high: 3.5, low: 3.375} + - {high: 3.625, low: 3.5} + - {high: 3.75, low: 3.625} + - {high: 3.875, low: 3.75} + - {high: 4.0, low: 3.875} + - {high: 4.25, low: 4.0} + - {high: 4.5, low: 4.25} diff --git a/buildmaster/LHCB_Z0_8TEV_MUON/rawdata/corrmat.corr b/buildmaster/LHCB_Z0_8TEV_MUON/rawdata/corrmat.corr new file mode 100644 index 0000000000..30782754d8 --- /dev/null +++ b/buildmaster/LHCB_Z0_8TEV_MUON/rawdata/corrmat.corr @@ -0,0 +1,34 @@ +1.0 +0.19 1.0 +0.17 0.27 1.0 +0.16 0.26 0.28 1.0 +0.16 0.25 0.28 0.29 1.0 +0.15 0.24 0.27 0.29 0.3 1.0 +0.14 0.23 0.26 0.28 0.29 0.3 1.0 +0.14 0.21 0.25 0.27 0.29 0.3 0.3 1.0 +0.13 0.2 0.23 0.25 0.27 0.28 0.29 0.29 1.0 +0.11 0.17 0.2 0.23 0.25 0.26 0.27 0.28 0.27 1.0 +0.09 0.14 0.16 0.18 0.2 0.22 0.22 0.23 0.23 0.23 1.0 +0.08 0.12 0.15 0.17 0.19 0.2 0.21 0.22 0.22 0.22 0.2 1.0 +0.07 0.1 0.12 0.14 0.16 0.17 0.18 0.19 0.19 0.2 0.18 0.19 1.0 +0.06 0.08 0.1 0.11 0.13 0.14 0.15 0.16 0.16 0.17 0.16 0.16 0.15 1.0 +0.05 0.07 0.08 0.09 0.1 0.11 0.11 0.12 0.13 0.13 0.12 0.13 0.12 0.11 1.0 +0.03 0.05 0.06 0.06 0.07 0.08 0.08 0.09 0.09 0.1 0.09 0.1 0.09 0.08 0.07 1.0 +0.03 0.04 0.04 0.05 0.05 0.06 0.06 0.06 0.07 0.07 0.07 0.08 0.07 0.07 0.06 0.05 1.0 +0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.02 0.02 0.01 0.01 0.01 0.01 1.0 +0.23 0.3 0.28 0.27 0.26 0.25 0.24 0.23 0.21 0.18 0.15 0.13 0.11 0.1 0.07 0.05 0.04 0.01 1.0 +0.05 0.15 0.21 0.2 0.2 0.2 0.2 0.19 0.18 0.16 0.14 0.12 0.1 0.08 0.06 0.04 0.03 0.01 0.2 1.0 +0.04 0.07 0.12 0.15 0.16 0.17 0.17 0.17 0.16 0.15 0.13 0.13 0.11 0.08 0.06 0.05 0.03 0.01 0.13 0.12 1.0 +0.05 0.08 0.1 0.13 0.16 0.17 0.17 0.17 0.16 0.16 0.14 0.13 0.12 0.09 0.07 0.05 0.03 0.01 0.06 0.03 0.26 1.0 +0.06 0.08 0.1 0.11 0.14 0.16 0.17 0.17 0.16 0.16 0.14 0.14 0.12 0.1 0.07 0.05 0.03 0.01 0.07 0.03 0.25 0.33 1.0 +0.04 0.06 0.07 0.09 0.1 0.11 0.12 0.13 0.13 0.12 0.11 0.11 0.09 0.08 0.06 0.04 0.03 0.0 0.03 -0.01 0.28 0.35 0.34 1.0 +0.04 0.06 0.07 0.08 0.09 0.1 0.11 0.12 0.12 0.12 0.11 0.11 0.1 0.09 0.07 0.05 0.04 0.0 0.0 -0.06 0.31 0.41 0.4 0.45 1.0 +0.02 0.03 0.04 0.04 0.04 0.04 0.05 0.05 0.06 0.06 0.06 0.07 0.06 0.06 0.05 0.04 0.04 0.01 -0.07 -0.14 0.14 0.32 0.32 0.35 0.40 1.0 +0.21 0.28 0.26 0.25 0.24 0.24 0.22 0.21 0.2 0.17 0.14 0.12 0.11 0.09 0.07 0.05 0.04 0.01 0.67 0.10 0.24 0.22 0.22 0.23 0.26 0.14 1.0 +0.04 0.14 0.19 0.18 0.18 0.18 0.18 0.17 0.16 0.15 0.12 0.11 0.09 0.07 0.05 0.04 0.03 0.0 0.07 0.54 0.23 0.28 0.28 0.28 0.33 0.24 0.21 1.0 +0.04 0.07 0.11 0.14 0.15 0.15 0.15 0.15 0.15 0.14 0.12 0.12 0.1 0.08 0.06 0.04 0.03 0.01 0.05 0.03 0.64 0.27 0.26 0.27 0.32 0.23 0.18 0.22 1.0 +0.05 0.07 0.09 0.12 0.14 0.15 0.15 0.16 0.15 0.14 0.12 0.12 0.11 0.08 0.06 0.04 0.03 0.01 0.04 0.00 0.25 0.70 0.30 0.33 0.39 0.29 0.21 0.25 0.31 1.0 +0.05 0.08 0.09 0.1 0.13 0.15 0.15 0.16 0.15 0.15 0.13 0.13 0.11 0.09 0.07 0.05 0.03 0.01 0.06 0.03 0.26 0.32 0.68 0.32 0.38 0.26 0.22 0.28 0.27 0.32 1.0 +0.04 0.06 0.08 0.09 0.1 0.12 0.13 0.13 0.13 0.13 0.11 0.11 0.1 0.08 0.06 0.04 0.03 0.0 0.07 0.04 0.30 0.28 0.27 0.63 0.36 0.22 0.23 0.23 0.29 0.32 0.28 1.0 +0.04 0.06 0.08 0.09 0.1 0.11 0.12 0.13 0.14 0.14 0.12 0.12 0.11 0.1 0.08 0.06 0.04 0.01 0.14 0.20 -0.13 -0.08 -0.07 -0.17 -0.04 -0.15 -0.06 -0.04 -0.04 -0.11 -0.07 -0.19 1.0 +0.03 0.04 0.04 0.05 0.05 0.06 0.06 0.06 0.07 0.07 0.07 0.08 0.08 0.07 0.06 0.05 0.04 0.01 0.12 0.17 -0.18 -0.17 -0.15 -0.23 -0.31 0.05 -0.09 -0.11 -0.14 -0.19 -0.18 -0.24 0.48 1.0 diff --git a/buildmaster/LHCB_Z0_8TEV_MUON/uncertainties.yaml b/buildmaster/LHCB_Z0_8TEV_MUON/uncertainties.yaml new file mode 100644 index 0000000000..e0e69fb3b3 --- /dev/null +++ b/buildmaster/LHCB_Z0_8TEV_MUON/uncertainties.yaml @@ -0,0 +1,752 @@ +definitions: + sys_corr_1: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_2: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_3: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_4: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_5: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_6: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_7: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_8: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_9: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_10: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_11: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_12: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_13: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_14: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_15: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_16: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_17: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + sys_corr_18: + description: Correlated systematic uncertainties + treatment: ADD + type: CORR + stat: + description: Uncorrelated statistical uncertainties + treatment: ADD + type: UNCORR + sys_beam: + description: Systematic Beam uncertainties + treatment: ADD + type: LHCBBEAM7TEV + sys_luminosity: + description: Systematic Luminosity uncertainties + treatment: ADD + type: LHCBLUMI7TEV +bins: +- sys_corr_1: -13.129371357446566 + sys_corr_2: 0.6853465481125789 + sys_corr_3: -1.8136153124139376 + sys_corr_4: -1.6378919565265173 + sys_corr_5: 1.0056447881503754 + sys_corr_6: 0.9250488794408771 + sys_corr_7: -0.6130610634103019 + sys_corr_8: -1.033777871097276 + sys_corr_9: 0.9022978319043655 + sys_corr_10: 0.5049756150293241 + sys_corr_11: -0.8765726686084679 + sys_corr_12: 0.5378206330642871 + sys_corr_13: -0.10895523093285057 + sys_corr_14: -0.4881555518108563 + sys_corr_15: -0.265568053839214 + sys_corr_16: 0.03401879609850996 + sys_corr_17: -11.515944076984733 + sys_corr_18: -9.94493656819381e-05 + sys_corr_19: -0.007574441556960995 + sys_corr_20: 0.0036078709274008384 + sys_corr_21: 0.054497988233331206 + sys_corr_22: -0.08061168467584123 + sys_corr_23: 0.2232461165909764 + sys_corr_24: 5.565016708385234 + sys_corr_25: 0.2609714050586262 + sys_corr_26: -39.974071914597864 + sys_corr_27: -27.943155261002897 + sys_corr_28: -16.32894214901613 + sys_corr_29: -4.436995222440807 + sys_corr_30: -0.8115865284675327 + sys_corr_31: -0.3755659425689436 + sys_corr_32: 1.9696013567201123 + sys_corr_33: -2.6324056596660883 + sys_corr_34: -0.1664720974747847 + stat: 33.0 + sys_beam: 14.0 + sys_luminosity: 14.0 +- sys_corr_1: -19.84410989915685 + sys_corr_2: -0.33428052767614547 + sys_corr_3: 2.6709227371602084 + sys_corr_4: -1.122291532865575 + sys_corr_5: 2.2445396406954448 + sys_corr_6: 1.2318412040028195 + sys_corr_7: -0.9643207591596343 + sys_corr_8: -2.2908607279681825 + sys_corr_9: 1.134280089906541 + sys_corr_10: -0.8548057471102127 + sys_corr_11: -1.796249557177666 + sys_corr_12: -0.07911083910570488 + sys_corr_13: -0.472951502025977 + sys_corr_14: 0.2558186095238292 + sys_corr_15: -0.3882769649810862 + sys_corr_16: -0.06521893670249582 + sys_corr_17: -20.03065356651147 + sys_corr_18: -5.938906840067905e-05 + sys_corr_19: -0.005317597654044067 + sys_corr_20: -0.012280961378499562 + sys_corr_21: 0.04723338963268214 + sys_corr_22: 0.08725074045417443 + sys_corr_23: -0.2200826518235283 + sys_corr_24: 10.334303066975695 + sys_corr_25: -1.8910305617977035 + sys_corr_26: 31.412187227642292 + sys_corr_27: -31.612351485478605 + sys_corr_28: -24.0763415362685 + sys_corr_29: -7.556695045723271 + sys_corr_30: 1.0636142561345379 + sys_corr_31: 1.678562392078886 + sys_corr_32: 3.831782907879968 + sys_corr_33: -4.904694970016218 + sys_corr_34: -1.4462770772783355 + stat: 51.0 + sys_beam: 38.0 + sys_luminosity: 38.0 +- sys_corr_1: -21.06803631144305 + sys_corr_2: -2.8187263609581716 + sys_corr_3: 7.429751053129648 + sys_corr_4: -1.2363854092311348 + sys_corr_5: 0.440020299276015 + sys_corr_6: 0.8285805636548081 + sys_corr_7: -0.1833307684413444 + sys_corr_8: -1.8866832461202852 + sys_corr_9: 1.0748601041255013 + sys_corr_10: -0.9841457536126791 + sys_corr_11: -2.2265738633784835 + sys_corr_12: -0.5655169182098072 + sys_corr_13: -0.4881568704060051 + sys_corr_14: -0.782035779352363 + sys_corr_15: -0.05919509595153284 + sys_corr_16: -1.1097250976141622 + sys_corr_17: -24.699858197476445 + sys_corr_18: -3.105265932964177e-05 + sys_corr_19: 0.0027633726426610577 + sys_corr_20: -0.01714976905989588 + sys_corr_21: 0.03260600131367636 + sys_corr_22: -0.02846572509875826 + sys_corr_23: -0.1467691264227021 + sys_corr_24: 12.565805182397174 + sys_corr_25: 0.625615759576152 + sys_corr_26: -3.5879461032450957 + sys_corr_27: 31.925361874882043 + sys_corr_28: -40.75747797840909 + sys_corr_29: -8.403699350561684 + sys_corr_30: -2.1525840708834774 + sys_corr_31: 1.1118986638848498 + sys_corr_32: 3.4014199295399585 + sys_corr_33: -4.750919333226475 + sys_corr_34: -1.0948035576743735 + stat: 62.0 + sys_beam: 57.0 + sys_luminosity: 58.0 +- sys_corr_1: -23.338793935211896 + sys_corr_2: -5.015020292076225 + sys_corr_3: 8.195207361706466 + sys_corr_4: -3.6167946329492535 + sys_corr_5: -1.5375896951455132 + sys_corr_6: -0.5055645332316372 + sys_corr_7: -0.1656605444022098 + sys_corr_8: -1.3674236995197913 + sys_corr_9: 0.9863174769320691 + sys_corr_10: -0.11689147629752011 + sys_corr_11: -2.76279796378281 + sys_corr_12: -0.7334770006813159 + sys_corr_13: -0.39710202721088006 + sys_corr_14: -0.9513343420516247 + sys_corr_15: -0.5236614050232143 + sys_corr_16: 0.14328103033362025 + sys_corr_17: -33.108074461181786 + sys_corr_18: -1.1676986084828315e-05 + sys_corr_19: -0.005199428393585041 + sys_corr_20: 0.0030204570342489963 + sys_corr_21: 0.04012076569205404 + sys_corr_22: 0.003803307862097004 + sys_corr_23: 0.13823101139063843 + sys_corr_24: 20.18490057862894 + sys_corr_25: 0.783534475773199 + sys_corr_26: -1.5428612589110309 + sys_corr_27: 2.9113392297877394 + sys_corr_28: 19.877343744892453 + sys_corr_29: -18.817522495155952 + sys_corr_30: 5.171011463669989 + sys_corr_31: 3.266548144634222 + sys_corr_32: 45.33543603498649 + sys_corr_33: -6.762197324713448 + sys_corr_34: -10.777895659545207 + stat: 70.0 + sys_beam: 77.0 + sys_luminosity: 78.0 +- sys_corr_1: -23.62879238801978 + sys_corr_2: -6.871934369335336 + sys_corr_3: 8.838236802397901 + sys_corr_4: -5.75372720918565 + sys_corr_5: -1.5402984336354577 + sys_corr_6: -0.3451192706498088 + sys_corr_7: 0.413105816684792 + sys_corr_8: 0.5073510714926635 + sys_corr_9: 1.5788941118756248 + sys_corr_10: 0.7607246969863836 + sys_corr_11: -2.673007365634146 + sys_corr_12: -1.2417474817490068 + sys_corr_13: -0.7432869137172015 + sys_corr_14: -0.598997451623972 + sys_corr_15: -0.06637383592033852 + sys_corr_16: 0.741162658809314 + sys_corr_17: -36.499510974556024 + sys_corr_18: 5.846611289286119e-06 + sys_corr_19: -0.0001752413952545087 + sys_corr_20: -0.008085268741514851 + sys_corr_21: 0.05044391052283411 + sys_corr_22: -0.15328563391070796 + sys_corr_23: 0.49110134937679384 + sys_corr_24: 17.63914440713869 + sys_corr_25: 1.5555770256950772 + sys_corr_26: -0.29419626869461024 + sys_corr_27: 1.4513100444746851 + sys_corr_28: 12.997331532175536 + sys_corr_29: -22.22192546156861 + sys_corr_30: 7.42411526077336 + sys_corr_31: -10.124571665104387 + sys_corr_32: -22.849337902782768 + sys_corr_33: -10.788104366418192 + sys_corr_34: 40.447688777367155 + stat: 76.0 + sys_beam: 93.0 + sys_luminosity: 94.0 +- sys_corr_1: -23.305638822180896 + sys_corr_2: -7.8937366993396365 + sys_corr_3: 9.253795984805592 + sys_corr_4: -7.241905559104537 + sys_corr_5: -1.9123722247234651 + sys_corr_6: 0.5912138549587218 + sys_corr_7: 0.0465224938402271 + sys_corr_8: 0.7199048005800637 + sys_corr_9: 1.828692488405049 + sys_corr_10: 2.1832257020458705 + sys_corr_11: -2.015036523294676 + sys_corr_12: -1.8353749829104327 + sys_corr_13: 0.0015187347231060773 + sys_corr_14: -0.48209896219515624 + sys_corr_15: 1.213798569671312 + sys_corr_16: 1.1211154805244627 + sys_corr_17: -36.92425033443672 + sys_corr_18: 1.4188601727116922e-05 + sys_corr_19: -0.008850403418020465 + sys_corr_20: -0.02487825307330309 + sys_corr_21: 0.09023808536309276 + sys_corr_22: -0.21273521070752027 + sys_corr_23: 0.568346566723835 + sys_corr_24: 9.560426327356453 + sys_corr_25: 1.6835388239353692 + sys_corr_26: -0.6834505876496619 + sys_corr_27: 0.6825047572018638 + sys_corr_28: 7.890847045554724 + sys_corr_29: -14.45530599394527 + sys_corr_30: 10.42790997030499 + sys_corr_31: -2.6562075756183785 + sys_corr_32: -29.80150469584712 + sys_corr_33: 18.833535008874055 + sys_corr_34: -40.257082690578 + stat: 79.0 + sys_beam: 103.0 + sys_luminosity: 105.0 +- sys_corr_1: -23.122870553655478 + sys_corr_2: -8.472992295351233 + sys_corr_3: 9.978884071326496 + sys_corr_4: -8.191892596977262 + sys_corr_5: -2.170985058688932 + sys_corr_6: 1.1927202283582423 + sys_corr_7: 1.027768527951832 + sys_corr_8: 1.3222563308413176 + sys_corr_9: 0.3390668514570481 + sys_corr_10: 2.5599804325753284 + sys_corr_11: -1.6074832724593935 + sys_corr_12: -1.6149874429471374 + sys_corr_13: -0.18163026461144396 + sys_corr_14: -1.317510032846632 + sys_corr_15: 1.2937203505127906 + sys_corr_16: 0.4372359769181125 + sys_corr_17: -39.002034940744274 + sys_corr_18: 1.5486962831135256e-05 + sys_corr_19: -0.006718922393629441 + sys_corr_20: -0.018187853006976003 + sys_corr_21: 0.059586532298960715 + sys_corr_22: -0.29423833975817293 + sys_corr_23: 0.8058671737702604 + sys_corr_24: 7.740143059284684 + sys_corr_25: 2.274379055527933 + sys_corr_26: -0.94388224336938 + sys_corr_27: 0.23546915242985275 + sys_corr_28: 4.67316228570499 + sys_corr_29: 31.63416661772606 + sys_corr_30: 4.348764566808575 + sys_corr_31: 47.20725640418786 + sys_corr_32: -6.675227803548107 + sys_corr_33: -7.705400923277242 + sys_corr_34: 4.070133302521287 + stat: 81.0 + sys_beam: 110.0 + sys_luminosity: 112.0 +- sys_corr_1: -20.85999640061773 + sys_corr_2: -8.242540680871082 + sys_corr_3: 8.961204056550143 + sys_corr_4: -8.614385898951557 + sys_corr_5: -2.306048360098447 + sys_corr_6: 1.3799189252318764 + sys_corr_7: 1.763104057565564 + sys_corr_8: 1.2044387625997617 + sys_corr_9: -0.029223529687805137 + sys_corr_10: 3.0215802371301006 + sys_corr_11: -1.4142992542156707 + sys_corr_12: -1.1889723630536118 + sys_corr_13: 0.9257611979897812 + sys_corr_14: -0.6372331880602672 + sys_corr_15: 0.3263646033902244 + sys_corr_16: 0.8296699988914757 + sys_corr_17: -34.31148751935087 + sys_corr_18: 2.223371222085629e-05 + sys_corr_19: -0.004079579657142941 + sys_corr_20: -0.04109755671287726 + sys_corr_21: 0.12179926772253029 + sys_corr_22: -0.4318081034916117 + sys_corr_23: 1.2383669027745627 + sys_corr_24: -0.555349911332747 + sys_corr_25: 3.847663593630981 + sys_corr_26: 1.214500716621149 + sys_corr_27: -1.0896588789845654 + sys_corr_28: 5.220822383009355 + sys_corr_29: 4.056607732576169 + sys_corr_30: -53.73393248419216 + sys_corr_31: -5.961938629073214 + sys_corr_32: 0.1723904768898964 + sys_corr_33: 13.699964798406059 + sys_corr_34: 2.433988171680721 + stat: 82.0 + sys_beam: 113.0 + sys_luminosity: 115.0 +- sys_corr_1: -20.223983250604604 + sys_corr_2: -8.154175836436359 + sys_corr_3: 9.26327356637858 + sys_corr_4: -9.718752686919935 + sys_corr_5: -2.815203846765515 + sys_corr_6: 1.8702305366091985 + sys_corr_7: 1.6422084351635182 + sys_corr_8: -0.16224786719614528 + sys_corr_9: -0.22452711746931522 + sys_corr_10: 3.572768532873175 + sys_corr_11: -1.807628621364004 + sys_corr_12: -0.9629409377418384 + sys_corr_13: 0.5171356641657681 + sys_corr_14: -0.6176549011257892 + sys_corr_15: 0.29396865202369565 + sys_corr_16: 0.19228110231051582 + sys_corr_17: -35.97332885136173 + sys_corr_18: 2.1030318953479127e-05 + sys_corr_19: -0.01573421063351599 + sys_corr_20: -0.03959349841824993 + sys_corr_21: 0.20298279711011666 + sys_corr_22: -0.43673348436197756 + sys_corr_23: 1.222881890842337 + sys_corr_24: -8.050128059309468 + sys_corr_25: 3.706780158719833 + sys_corr_26: 0.2856448791722523 + sys_corr_27: 0.3134908927527956 + sys_corr_28: 1.8316820476472904 + sys_corr_29: 32.11772042509791 + sys_corr_30: 6.0053293548508035 + sys_corr_31: -35.91989545508118 + sys_corr_32: -0.02011206044083529 + sys_corr_33: -30.62631510390598 + sys_corr_34: -10.705014255279856 + stat: 81.0 + sys_beam: 112.0 + sys_luminosity: 114.0 +- sys_corr_1: -17.04123512385607 + sys_corr_2: -7.918393116785677 + sys_corr_3: 8.440143768234751 + sys_corr_4: -10.024020107965285 + sys_corr_5: -2.6847034507320844 + sys_corr_6: 2.0671669249445555 + sys_corr_7: 1.837279687306533 + sys_corr_8: 0.3650568659958519 + sys_corr_9: -0.4904354192346455 + sys_corr_10: 3.0568684568954443 + sys_corr_11: -1.0638801573487726 + sys_corr_12: -1.3582456192662984 + sys_corr_13: 0.07169286786156893 + sys_corr_14: -0.2199043557753311 + sys_corr_15: 1.5744304454334228 + sys_corr_16: -0.0013292406672338916 + sys_corr_17: -32.11721511436332 + sys_corr_18: 1.563925464777089e-05 + sys_corr_19: -0.01602557329814359 + sys_corr_20: -0.06980990577079062 + sys_corr_21: 0.2132859532473274 + sys_corr_22: -0.6435224332446 + sys_corr_23: 1.9689443448390334 + sys_corr_24: -12.899652788631762 + sys_corr_25: 5.004417879944764 + sys_corr_26: 0.4905581183348719 + sys_corr_27: -0.6767901157597395 + sys_corr_28: -5.192725175071233 + sys_corr_29: 15.90697875227382 + sys_corr_30: 20.16293066600603 + sys_corr_31: -11.987331589786724 + sys_corr_32: 14.255145302365655 + sys_corr_33: 43.78533937999355 + sys_corr_34: 17.535461370942752 + stat: 78.0 + sys_beam: 104.0 + sys_luminosity: 105.0 +- sys_corr_1: -14.949658691661552 + sys_corr_2: -7.165496747473827 + sys_corr_3: 7.808477726717307 + sys_corr_4: -9.56305409326678 + sys_corr_5: -2.670729336856583 + sys_corr_6: 1.870910013082464 + sys_corr_7: 1.1343784857727517 + sys_corr_8: 0.5731050576349122 + sys_corr_9: -0.8237338552005756 + sys_corr_10: 3.4751925833597936 + sys_corr_11: -0.13851422837746447 + sys_corr_12: -0.10023616022750506 + sys_corr_13: -0.297566757573512 + sys_corr_14: -0.01650093308831284 + sys_corr_15: 0.6964491524517722 + sys_corr_16: 0.020654173191060936 + sys_corr_17: -29.914143247343116 + sys_corr_18: -2.595663656276923e-06 + sys_corr_19: -0.020768693638988766 + sys_corr_20: -0.05834119476442283 + sys_corr_21: 0.19884628619594974 + sys_corr_22: -0.6436033819979039 + sys_corr_23: 1.6165938253859384 + sys_corr_24: -56.831154776699634 + sys_corr_25: 4.207755805618636 + sys_corr_26: 0.13546165938685006 + sys_corr_27: 0.16794567880689754 + sys_corr_28: -1.254823525202448 + sys_corr_29: -23.915653026107577 + sys_corr_30: 1.316650496021619 + sys_corr_31: 12.59376213706475 + sys_corr_32: 1.475379610059091 + sys_corr_33: -11.882612801600084 + sys_corr_34: -0.6558219368176063 + stat: 72.0 + sys_beam: 89.0 + sys_luminosity: 90.0 +- sys_corr_1: -9.14904976731458 + sys_corr_2: -4.972284733299394 + sys_corr_3: 4.934501253783746 + sys_corr_4: -7.146715955712768 + sys_corr_5: -2.4806832426488516 + sys_corr_6: 1.6716036517481274 + sys_corr_7: 0.3731925008926304 + sys_corr_8: 0.46906925363558766 + sys_corr_9: -1.0870814502316777 + sys_corr_10: 2.1864260496234604 + sys_corr_11: -0.08527632330590701 + sys_corr_12: -0.012526363709937075 + sys_corr_13: 0.424383407931963 + sys_corr_14: -0.13456579621486817 + sys_corr_15: 0.479004540879783 + sys_corr_16: -0.15340419169042 + sys_corr_17: -14.879810632978513 + sys_corr_18: -0.00024564680501881395 + sys_corr_19: -0.0454435934292143 + sys_corr_20: -0.12517695573643947 + sys_corr_21: 0.4495879379661885 + sys_corr_22: -1.2047889254655135 + sys_corr_23: 6.991666808010768 + sys_corr_24: -5.151973939331355 + sys_corr_25: -45.641593280488614 + sys_corr_26: -1.3784434779206556 + sys_corr_27: 1.565294265899055 + sys_corr_28: 1.4093209031809397 + sys_corr_29: 2.3280091228437603 + sys_corr_30: -0.7987670255568619 + sys_corr_31: -1.492549105845241 + sys_corr_32: 0.1634737023690217 + sys_corr_33: 2.02289736511137 + sys_corr_34: 1.008066369978752 + stat: 63.0 + sys_beam: 70.0 + sys_luminosity: 71.0 +- sys_corr_1: -6.326085258885516 + sys_corr_2: -3.439415506883048 + sys_corr_3: 3.3413277376739843 + sys_corr_4: -5.541898966162248 + sys_corr_5: -1.666176480784771 + sys_corr_6: 0.979375217874185 + sys_corr_7: 0.0002485455875154854 + sys_corr_8: -0.08187806243705073 + sys_corr_9: -0.43331484389360025 + sys_corr_10: 1.7688356519672581 + sys_corr_11: 0.16738920479722036 + sys_corr_12: 0.11055325680962605 + sys_corr_13: 0.4269058132515623 + sys_corr_14: -0.1138800538385906 + sys_corr_15: 0.9613977606946876 + sys_corr_16: -0.013820528956704024 + sys_corr_17: -9.506358516911547 + sys_corr_18: -0.000322491958865001 + sys_corr_19: -0.04499426072423848 + sys_corr_20: -0.13755544896289332 + sys_corr_21: 0.5827172665108897 + sys_corr_22: -1.9470722746536389 + sys_corr_23: -37.513389438295306 + sys_corr_24: -3.666147114225558 + sys_corr_25: -7.529167532066956 + sys_corr_26: -0.6017477237226491 + sys_corr_27: 0.16993176961810133 + sys_corr_28: 0.8185166967637814 + sys_corr_29: 1.5139539639474104 + sys_corr_30: -0.23812496924184945 + sys_corr_31: -0.8615737285431316 + sys_corr_32: 0.09807293782358924 + sys_corr_33: 1.5488833831456272 + sys_corr_34: 0.7674261133314435 + stat: 54.0 + sys_beam: 50.0 + sys_luminosity: 51.0 +- sys_corr_1: -3.581622334817089 + sys_corr_2: -1.5952014747683092 + sys_corr_3: 1.696419758223885 + sys_corr_4: -3.146991276459658 + sys_corr_5: -0.8213111707130737 + sys_corr_6: 0.8751706472028377 + sys_corr_7: 0.2206946879339562 + sys_corr_8: -0.10396281407173584 + sys_corr_9: -0.6365967085276881 + sys_corr_10: 0.9100948097624815 + sys_corr_11: 0.13025444043125295 + sys_corr_12: 0.34852058051345974 + sys_corr_13: 0.000963419400960771 + sys_corr_14: -0.024415407833905774 + sys_corr_15: 0.23683126623679324 + sys_corr_16: -0.213123518257198 + sys_corr_17: -4.887880950067507 + sys_corr_18: -9.666121427649194e-05 + sys_corr_19: -0.08679639349439726 + sys_corr_20: -0.21795914219962417 + sys_corr_21: 1.5906032250675446 + sys_corr_22: 25.605815007759354 + sys_corr_23: -2.2802107816509967 + sys_corr_24: -2.064997634263775 + sys_corr_25: -2.1904506948921627 + sys_corr_26: -0.31575888365863575 + sys_corr_27: 0.13593361557791478 + sys_corr_28: 0.2645398552247651 + sys_corr_29: 0.7254250538696598 + sys_corr_30: -0.14873611967742692 + sys_corr_31: -0.34721727099001154 + sys_corr_32: -0.06366415037975065 + sys_corr_33: 0.7056752901669989 + sys_corr_34: 0.34189479965746694 + stat: 42.0 + sys_beam: 31.0 + sys_luminosity: 32.0 +- sys_corr_1: -1.92704439512093 + sys_corr_2: -0.8906975447977756 + sys_corr_3: 0.9859044565822679 + sys_corr_4: -1.900006590186975 + sys_corr_5: -0.500061597538644 + sys_corr_6: 0.46968476717803787 + sys_corr_7: 0.026567509007107172 + sys_corr_8: -0.24644151318871022 + sys_corr_9: -0.25414422874330955 + sys_corr_10: 0.5977450748892065 + sys_corr_11: 0.12461147105302377 + sys_corr_12: 0.22106489547048022 + sys_corr_13: 0.002600996773599335 + sys_corr_14: 0.18567315840954113 + sys_corr_15: 0.11880756529326222 + sys_corr_16: -0.18567999103818122 + sys_corr_17: -2.8272588693821663 + sys_corr_18: -0.00019905744405662599 + sys_corr_19: -0.11219010309646388 + sys_corr_20: -0.36203400976763017 + sys_corr_21: -19.33064816976645 + sys_corr_22: 1.9792982895260696 + sys_corr_23: -1.076070580397553 + sys_corr_24: -0.9794399369760031 + sys_corr_25: -1.2690397198581644 + sys_corr_26: -0.11077568336747423 + sys_corr_27: -0.0455017730958743 + sys_corr_28: 0.014709472995496412 + sys_corr_29: 0.3343664195538066 + sys_corr_30: 0.011541248946823261 + sys_corr_31: -0.3839086612900922 + sys_corr_32: 0.06978441016448338 + sys_corr_33: 0.26523110211595075 + sys_corr_34: 0.06515123861157356 + stat: 32.0 + sys_beam: 18.0 + sys_luminosity: 19.0 +- sys_corr_1: -0.8296601388639502 + sys_corr_2: -0.3805251722782747 + sys_corr_3: 0.4256454984305268 + sys_corr_4: -0.8640659174103269 + sys_corr_5: -0.3029421108714545 + sys_corr_6: 0.2176950458024982 + sys_corr_7: -0.0513015110495244 + sys_corr_8: -0.12919103018667238 + sys_corr_9: -0.10721047160309029 + sys_corr_10: 0.1529847335090339 + sys_corr_11: 0.17487604901206083 + sys_corr_12: -0.01248019049416586 + sys_corr_13: 0.04226511417819689 + sys_corr_14: 0.07817088485442444 + sys_corr_15: 0.0736083183745245 + sys_corr_16: -0.10963219724155511 + sys_corr_17: -1.2120042253538708 + sys_corr_18: -0.00046489040093820644 + sys_corr_19: -0.25803633958738176 + sys_corr_20: 11.795834334234783 + sys_corr_21: -0.5401002869203217 + sys_corr_22: 0.47787211458447243 + sys_corr_23: -0.4012872851787313 + sys_corr_24: -0.4773980799989664 + sys_corr_25: -0.5578977015713101 + sys_corr_26: 0.015910316662032797 + sys_corr_27: 0.039315340097579664 + sys_corr_28: -0.034916993021117164 + sys_corr_29: 0.15079334377683432 + sys_corr_30: -0.025448127463030965 + sys_corr_31: -0.12947869114049182 + sys_corr_32: 0.003045495350807366 + sys_corr_33: 0.21539792088791174 + sys_corr_34: 0.050432772372592215 + stat: 22.0 + sys_beam: 9.0 + sys_luminosity: 9.0 +- sys_corr_1: -0.42393728550091253 + sys_corr_2: -0.14894396859548853 + sys_corr_3: 0.1981748571052791 + sys_corr_4: -0.39930231793851473 + sys_corr_5: -0.11765515578109007 + sys_corr_6: 0.10694488123136513 + sys_corr_7: -0.0737090810991067 + sys_corr_8: -0.13441755182880608 + sys_corr_9: -0.11536049109165246 + sys_corr_10: 0.04875973725724547 + sys_corr_11: 0.04521689153192965 + sys_corr_12: 0.10130389568753956 + sys_corr_13: 0.04322115716786899 + sys_corr_14: 0.028136185353375684 + sys_corr_15: 0.004731362627832487 + sys_corr_16: -0.11516055102024063 + sys_corr_17: -0.6185243817055507 + sys_corr_18: -0.0008188847906671834 + sys_corr_19: 7.911038971966661 + sys_corr_20: 0.37528005595773883 + sys_corr_21: -0.26672538386835976 + sys_corr_22: 0.3020038060076616 + sys_corr_23: -0.21369897814298744 + sys_corr_24: -0.2553303274366322 + sys_corr_25: -0.33152188617632095 + sys_corr_26: -0.03180232046108036 + sys_corr_27: -0.045026544656933135 + sys_corr_28: 0.015861015519060512 + sys_corr_29: 0.06645220360768585 + sys_corr_30: 0.039901597979006935 + sys_corr_31: -0.053250101359537255 + sys_corr_32: 0.028197119479147756 + sys_corr_33: 0.04831256818640388 + sys_corr_34: -0.01828330836598994 + stat: 16.0 + sys_beam: 4.0 + sys_luminosity: 4.0 +- sys_corr_1: -0.013414184810947182 + sys_corr_2: -0.0037960802643431835 + sys_corr_3: 0.007050325110555244 + sys_corr_4: -0.012793678842011953 + sys_corr_5: -0.006988237873416711 + sys_corr_6: -0.0011527303922321177 + sys_corr_7: -0.0038282744129991384 + sys_corr_8: 0.00580351666510989 + sys_corr_9: 0.010110892923237483 + sys_corr_10: 0.006099701507322768 + sys_corr_11: -0.002713842573751903 + sys_corr_12: 0.0038844748413824986 + sys_corr_13: -0.00021929479123481416 + sys_corr_14: -5.673249487341189e-05 + sys_corr_15: -0.005283439658346962 + sys_corr_16: -0.010413435144737736 + sys_corr_17: -0.009022711855041334 + sys_corr_18: 0.9993144938282417 + sys_corr_19: 0.006305972548180469 + sys_corr_20: 0.005629159396434368 + sys_corr_21: -0.003870372244719926 + sys_corr_22: 0.0024480503883864435 + sys_corr_23: -0.011281196460354895 + sys_corr_24: -0.0016038955252630446 + sys_corr_25: -0.015001223022070295 + sys_corr_26: -0.0028588158298082783 + sys_corr_27: -0.0032024298307163834 + sys_corr_28: -0.00378958229068439 + sys_corr_29: -0.0012810170230193785 + sys_corr_30: 0.00020952138697761106 + sys_corr_31: -0.0002487670121694488 + sys_corr_32: 0.0015993114173559364 + sys_corr_33: 1.83609916547212e-05 + sys_corr_34: 0.00045379693968826006 + stat: 3.0 + sys_beam: 0.0 + sys_luminosity: 0.0 diff --git a/buildmaster/dataset_names.yml b/buildmaster/dataset_names.yml index e56b12e016..3c1785a980 100644 --- a/buildmaster/dataset_names.yml +++ b/buildmaster/dataset_names.yml @@ -6,3 +6,6 @@ LHCB_Z0_13TEV_DIMUON: LHCB_Z_13TEV_DIMUON LHCB_Z0_7TEV_MUON: LHCBZMU7TEV LHCB_WPWM_7TEV_MUON: LHCBWMU7TEV LHCB_DY_7TEV_MUON: LHCBWZMU7TEV +LHCB_Z0_8TEV_MUON: LHCBZMU8TEV +LHCB_WPWM_8TEV_MUON: LHCBWMU8TEV +LHCB_DY_8TEV_MUON: LHCBWZMU8TEV