diff --git a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/data_legacy_WPWM-RATIO.yaml b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/data_WPWM-RATIO.yaml similarity index 100% rename from nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/data_legacy_WPWM-RATIO.yaml rename to nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/data_WPWM-RATIO.yaml diff --git a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/data_legacy_WPWM-TOT.yaml b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/data_WPWM-TOT.yaml similarity index 100% rename from nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/data_legacy_WPWM-TOT.yaml rename to nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/data_WPWM-TOT.yaml diff --git a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/filter.py b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/filter.py new file mode 100644 index 0000000000..c0604b80d7 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/filter.py @@ -0,0 +1,19 @@ +''' +Filter script for CMS_WCHARM_7TEV +''' + +import logging +import os + +from filter_utils import Extractor + +logging.basicConfig(level=logging.INFO, format='[%(levelname)s] %(message)s') + +CURRENT_DIR = os.path.dirname(os.path.abspath(__file__)) + +if __name__ == "__main__": + CMS_WCHARM_TOT = Extractor(f"{CURRENT_DIR}/metadata.yaml", "WPWM-TOT", mult_factor=1000) + CMS_WCHARM_TOT.generate_data() + + CMS_WCHARM_RATIO = Extractor(f"{CURRENT_DIR}/metadata.yaml", "WPWM-RATIO", mult_factor=1.0) + CMS_WCHARM_RATIO.generate_data() diff --git a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/filter_utils.py b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/filter_utils.py new file mode 100644 index 0000000000..691d251bdd --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/filter_utils.py @@ -0,0 +1,260 @@ +import functools +import logging +import os + +import numpy as np +import yaml + +from nnpdf_data.filter_utils.utils import prettify_float + +yaml.add_representer(float, prettify_float) + +MW2 = 80.385**2 + +CURRENT_DIR = os.path.dirname(os.path.abspath(__file__)) +ART_LABEL = 'art_corr' +STAT_LABEL = 'stat_uncorr' +TABLE_TOKEN = 'Table' + + +class Extractor: + + def __init__(self, metadata_file, observable, mult_factor=1): + """ + Parameters + ---------- + metadata_file: str + Path to the metadata file + observable: str + The name of the observable for which the data is extracted. The name + must be listed in the metadata file. + mult_factor: float + Multiplication factor to apply to the central data points. This is + useful to convert the data in the metadata file to the desired + units. + """ + # Open metadata and select process + with open(metadata_file, 'r') as file: + metadata = yaml.safe_load(file) + self.metadata = next( + ( + md + for md in metadata["implemented_observables"] + if md['observable_name'] == observable + ), + None, + ) + if self.metadata is None: + raise Exception(f"{observable} is not listed in the metadata file.") + + self.observable = observable + self.mult_factor = mult_factor + + @functools.cache + def _retrieve_table(self, table_id): + """ + Implementation of the loading for the table. + + Parameters + ---------- + table_id: int + Index that specifies the table. + + Return + ------ + The table specified by `table_id`. + """ + with open(f'{CURRENT_DIR}/rawdata/{TABLE_TOKEN}{table_id}.yaml') as tab: + tab_dict = yaml.safe_load(tab) + return tab_dict + + def _generate_kinematics(self): + """ + The function generates the kinematics by reading and processing it from + the referenced table. Kinematics is processed in the format of a list of + dictionaries. The keys in each dictionaries specify the label (i.e. name) + for the kinematic variables. For this dataset, they are 'abs_eta' and 'm_W2'. + The labels are taken from the matadata file. The corresponding values are + 'min', 'mid', and 'max'. + + For this dataset, 'm_W2' is used in the computation of the (x,Q2)-map and + does not have any active role in the fit. For that reason, every bin has the + same value. Moreover, only the mid value is used. + """ + logging.info(f"Generating kinematics for CMS_{self.observable}...") + + table_ID = self.metadata["tables"][0] + tab_dict = self._retrieve_table(table_ID) + + data = tab_dict['independent_variables'][0] + label = self.metadata['kinematic_coverage'] + kinematics = [] + for eta_bin in data['values']: + abs_eta_max = eta_bin['high'] + abs_eta_min = eta_bin['low'] + kin_bin = { + label[0]: { + 'min': abs_eta_min, + 'mid': (abs_eta_max + abs_eta_min) / 2, + 'max': abs_eta_max, + }, + label[1]: {'min': None, 'mid': MW2, 'max': None}, + } + kinematics.append(kin_bin) + + # Check number of data agrees with metadata + ndata = len(kinematics) + if not self.metadata['ndata'] == ndata: + raise ValueError( + f"Mismatch in 'ndata': expected {self.metadata['ndata']}, but got {ndata}" + ) + self.ndata = ndata + return kinematics + + def _generate_data_and_unc(self): + """ + Return a list with central data points and two additional lists with the corresponding + statistical and systematic uncertainties. For this dataset, uncertainties are always + symmetric. Uncertainties are given as absolute values. + + Note that, for the total x-sec, the correlation matrix is provided. The corresponding + covariance matrix is constructed in `_generate_covmat`. + """ + logging.info(f"Generating central data for CMS_{self.observable}...") + dat_central = [] + stat_unc = [] + asy_sys_unc = [] + table_ID = self.metadata['tables'][0] + tab_dict = self._retrieve_table(table_ID) + + # Select data with pT > 25 GeV + tab_dict = tab_dict['dependent_variables'][0]['values'] + + # Loop over bins + for rap_bin in tab_dict: + dat_central.append(rap_bin['value'] * self.mult_factor) + stat_unc.append(rap_bin['errors'][0]['symerror'] * self.mult_factor) + asy_sys_unc.append(rap_bin['errors'][1]['symerror'] * self.mult_factor) + return dat_central, stat_unc, asy_sys_unc + + def _build_unc_definitions(self): + """ + Build the dictionary containing the definitions of the uncertainties to be + used in the uncertainty data file. + """ + unc_definitions = {} + + # Statistical uncertainty + unc_definitions[STAT_LABEL] = { + 'description': f'Statistical uncertainty', + 'treatment': 'ADD', + 'type': 'UNCORR', + } + + if self.observable == 'WPWM-RATIO': + unc_definitions['ART_LABEL'] = { + 'description': f'Correlated systematic uncertainty', + 'treatment': 'MULT', + 'type': 'CORR', + } + elif self.observable == 'WPWM-TOT': + for idx in range(self.ndata): + unc_definitions[f'{ART_LABEL}_{idx+1}'] = { + 'description': f'Correlated systematic uncertainty {idx+1}', + 'treatment': 'ADD', + 'type': 'CORR', + } + + return unc_definitions + + def _generate_covmat(self, diag_uncs): + """ + Generate the covariance matrix for the total x-sec. This function requires + the diagonal systematic uncertainties as argument. The diagonal uncertainties + are used to construct the covariance matrix from the correlation matrix stored + in the HepData table. + + Note that such a correlation matrix exists for the total x-sec only, while the + ratio observable does not provide this information. + """ + if not self.observable == 'WPWM-TOT': + raise ValueError( + "The construction of the covariance matrix is defined for the total x-sec only." + ) + table_ID = self.metadata["tables"][1] + tab_dict = self._retrieve_table(table_ID) + matlist = tab_dict['dependent_variables'][0]['values'] + matlist = [d['value'] for d in matlist] + covmat = np.zeros((self.ndata, self.ndata)) + for i in range(self.ndata): + for j in range(self.ndata): + covmat[i, j] = matlist[i + self.ndata * j] * diag_uncs[i] * diag_uncs[j] + return covmat + + def generate_data(self): + """ + The function collects central data, kinematics, and uncertainties ans save them + into yaml files. + + The systematic uncertainties are given as percentages relative the central data point. + The absolute value of the uncertainty is obtained from the central data point before + the shifts are applied. + """ + # Get central data and kinematics + central_data, stat_unc, sys_unc = self._generate_data_and_unc() + kinematics = self._generate_kinematics() + + # Uncertainty definitions + unc_definitions = self._build_unc_definitions() + sys_artificial = [] # Initialize vector of artificial uncertainties + + if self.observable == 'WPWM-TOT': + # Generate covmat and perform eigen decomposition + covmat = self._generate_covmat(sys_unc) + eigvals, eigvecs = np.linalg.eig(covmat) + art_unc = np.sqrt(eigvals) * eigvecs + + # Loop over bins + for data_idx in range(len(central_data)): + # Statistical uncertainty + unc_dict = {STAT_LABEL: stat_unc[data_idx]} + + # Artificial systematic uncertainties + for sys_idx, art_sys in enumerate(art_unc[data_idx, :]): + unc_dict[f'{ART_LABEL}_{sys_idx+1}'] = float(art_sys) + + # Append to list + sys_artificial.append(unc_dict) + + elif self.observable == 'WPWM-RATIO': + for data_idx in range(len(central_data)): + # Statistical uncertainty + unc_dict = {STAT_LABEL: stat_unc[data_idx]} + + # Systematic uncertainty + unc_dict[f'{ART_LABEL}'] = sys_unc[data_idx] + sys_artificial.append(unc_dict) + + # Save kinematics into file + logging.info("Dumping kinematics to file...") + kinematics_yaml = {'bins': kinematics} + kins_file_name = self.metadata['kinematics']['file'] + with open(CURRENT_DIR + '/' + kins_file_name, 'w') as file: + yaml.dump(kinematics_yaml, file, sort_keys=False) + logging.info("Done!") + + # Save central data into file + logging.info("Dumping kinematics to file...") + dat_central_yaml = {'data_central': central_data} + data_file_name = self.metadata['data_central'] + with open(CURRENT_DIR + '/' + data_file_name, 'w') as file: + yaml.dump(dat_central_yaml, file, sort_keys=False) + logging.info("Done!") + + # Save unertainties + logging.info("Dumping kinematics to file...") + uncertainties_yaml = {'definitions': unc_definitions, 'bins': sys_artificial} + unc_file_name = self.metadata['data_uncertainties'][0] + with open(CURRENT_DIR + '/' + unc_file_name, 'w') as file: + yaml.dump(uncertainties_yaml, file, sort_keys=False) + logging.info("Done!") diff --git a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/kinematics_WPWM-RATIO.yaml b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/kinematics_WPWM-RATIO.yaml index d56da464e5..b3549529e5 100644 --- a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/kinematics_WPWM-RATIO.yaml +++ b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/kinematics_WPWM-RATIO.yaml @@ -1,61 +1,41 @@ bins: -- k1: - min: null +- abs_eta: + min: 0.0 mid: 0.175 - max: null - k2: - min: null - mid: 6463.838404 - max: null - k3: - min: null - mid: 7000.0 - max: null -- k1: - min: null - mid: 0.525 - max: null - k2: + max: 0.35 + m_W2: min: null - mid: 6463.838404 + mid: 6.46174823e+03 max: null - k3: +- abs_eta: + min: 0.35 + mid: 5.25000000e-01 + max: 0.7 + m_W2: min: null - mid: 7000.0 + mid: 6.46174823e+03 max: null -- k1: - min: null +- abs_eta: + min: 0.7 mid: 0.9 - max: null - k2: + max: 1.1 + m_W2: min: null - mid: 6463.838404 + mid: 6.46174823e+03 max: null - k3: - min: null - mid: 7000.0 - max: null -- k1: - min: null +- abs_eta: + min: 1.1 mid: 1.35 - max: null - k2: - min: null - mid: 6463.838404 - max: null - k3: + max: 1.6 + m_W2: min: null - mid: 7000.0 + mid: 6.46174823e+03 max: null -- k1: - min: null +- abs_eta: + min: 1.6 mid: 1.85 - max: null - k2: - min: null - mid: 6463.838404 - max: null - k3: + max: 2.1 + m_W2: min: null - mid: 7000.0 + mid: 6.46174823e+03 max: null diff --git a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/kinematics_WPWM-TOT.yaml b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/kinematics_WPWM-TOT.yaml index d56da464e5..b3549529e5 100644 --- a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/kinematics_WPWM-TOT.yaml +++ b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/kinematics_WPWM-TOT.yaml @@ -1,61 +1,41 @@ bins: -- k1: - min: null +- abs_eta: + min: 0.0 mid: 0.175 - max: null - k2: - min: null - mid: 6463.838404 - max: null - k3: - min: null - mid: 7000.0 - max: null -- k1: - min: null - mid: 0.525 - max: null - k2: + max: 0.35 + m_W2: min: null - mid: 6463.838404 + mid: 6.46174823e+03 max: null - k3: +- abs_eta: + min: 0.35 + mid: 5.25000000e-01 + max: 0.7 + m_W2: min: null - mid: 7000.0 + mid: 6.46174823e+03 max: null -- k1: - min: null +- abs_eta: + min: 0.7 mid: 0.9 - max: null - k2: + max: 1.1 + m_W2: min: null - mid: 6463.838404 + mid: 6.46174823e+03 max: null - k3: - min: null - mid: 7000.0 - max: null -- k1: - min: null +- abs_eta: + min: 1.1 mid: 1.35 - max: null - k2: - min: null - mid: 6463.838404 - max: null - k3: + max: 1.6 + m_W2: min: null - mid: 7000.0 + mid: 6.46174823e+03 max: null -- k1: - min: null +- abs_eta: + min: 1.6 mid: 1.85 - max: null - k2: - min: null - mid: 6463.838404 - max: null - k3: + max: 2.1 + m_W2: min: null - mid: 7000.0 + mid: 6.46174823e+03 max: null diff --git a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/metadata.yaml b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/metadata.yaml index a741e0159d..2e60beeb07 100644 --- a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/metadata.yaml +++ b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/metadata.yaml @@ -1,58 +1,55 @@ setname: CMS_WCHARM_7TEV -version: 1 -version_comment: Port of old commondata + nnpdf_metadata: nnpdf31_process: DY CC experiment: CMS + arXiv: url: https://arxiv.org/abs/1310.1138 journal: JHEP 02 (2014) 013 iNSPIRE: - url: '' + url: https://inspirehep.net/literature/1256938 hepdata: - url: '' - version: -1 + url: https://www.hepdata.net/record/ins1256938 + version: 1 + +version: 2 +version_comment: Implementation in the new format + implemented_observables: - observable_name: WPWM-RATIO observable: description: Jet Rapidity Distribution label: CMS $W+c$ ratio - units: '' - process_type: EWJ_RAP - tables: [] - npoints: [] + units: '[fb]' + process_type: DY_W_ETA + tables: [9] ndata: 5 plotting: kinematics_override: ewj_rap_sqrt_scale dataset_label: CMS $W+c$ ratio y_label: $\sigma(W^+ + \bar{c})/\sigma(W^- + c)$ - plot_x: k1 - kinematic_coverage: - - k1 - - k2 - - k3 + plot_x: abs_eta + kinematic_coverage: [abs_eta, m_W2] kinematics: variables: - k1: - description: Variable k1 - label: k1 - units: '' - k2: - description: Variable k2 - label: k2 - units: '' - k3: - description: Variable k3 - label: k3 + abs_eta: + description: Absolute pseudo-rapidity of the Z boson + label: $|\eta|$ units: '' + m_W2: + description: Mass of the W boson squared + label: $m_W^2$ + units: GeV$^{2}$ file: kinematics_WPWM-RATIO.yaml + data_central: data_WPWM-RATIO.yaml + data_uncertainties: [uncertainties_WPWM-RATIO.yaml] theory: conversion_factor: 1.0 operation: ratio FK_tables: - - CMSWCHARMRAT-CMSWCHARM-WpCb-eta4 - - CMSWCHARMRAT-CMSWCHARM-WmC-eta3 - data_uncertainties: [] variants: legacy: data_uncertainties: @@ -60,40 +57,34 @@ implemented_observables: legacy_10: data_uncertainties: - uncertainties_WPWM-RATIO_sys_10.yaml - data_central: data_legacy_WPWM-RATIO.yaml ported_from: CMSWCHARMRAT + - observable_name: WPWM-TOT observable: - description: Jet Rapidity Distribution + description: identity label: CMS $W+c$ total - units: '' - process_type: EWJ_RAP - tables: [] - npoints: [] + units: '[fb]' + process_type: DY_W_ETA + tables: [5, 6] ndata: 5 plotting: - kinematics_override: ewj_rap_sqrt_scale + kinematics_override: identity dataset_label: CMS $W+c$ total y_label: $d\sigma(W+c)/d|\eta_l|$ (fb) plot_x: k1 - kinematic_coverage: - - k1 - - k2 - - k3 + kinematic_coverage: [abs_eta, m_W2] + data_uncertainties: [uncertainties_WPWM-TOT.yaml] + data_central: data_WPWM-TOT.yaml kinematics: variables: - k1: - description: Variable k1 - label: k1 - units: '' - k2: - description: Variable k2 - label: k2 - units: '' - k3: - description: Variable k3 - label: k3 + abs_eta: + description: Absolute pseudo-rapidity of the Z boson + label: $|\eta|$ units: '' + m_W2: + description: Mass of the W boson squared + label: $m_W^2$ + units: GeV$^{2}$ file: kinematics_WPWM-TOT.yaml theory: conversion_factor: 1.0 @@ -101,7 +92,6 @@ implemented_observables: FK_tables: - - CMSWCHARM-CMSWCHARM-WmC-eta3 - - CMSWCHARM-CMSWCHARM-WpCb-eta4 - data_uncertainties: [] variants: legacy: data_uncertainties: @@ -109,5 +99,4 @@ implemented_observables: legacy_10: data_uncertainties: - uncertainties_WPWM-TOT_sys_10.yaml - data_central: data_legacy_WPWM-TOT.yaml ported_from: CMSWCHARMTOT diff --git a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/rawdata/Table5.yaml b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/rawdata/Table5.yaml new file mode 100755 index 0000000000..2d329707b5 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/rawdata/Table5.yaml @@ -0,0 +1,69 @@ +dependent_variables: +- header: {name: D(SIG)/DABS(ETARAP(LEPTON)), units: PB} + qualifiers: + - {name: ABS(ETARAP(JET)), value: < 2.5} + - {name: JETS, value: 'ANTI-KT, R=1'} + - {name: PT(JET), units: GEV, value: '> 25'} + - {name: PT(LEPTON), units: GEV, value: '> 25'} + - {name: RE, value: P P --> W+ < LEPTON+ NU > CHARMBAR X} + - {name: RE, value: P P --> W- < LEPTON- NUBAR > CHARM X} + - {name: SQRT(S), units: GeV, value: '7000.0'} + values: + - errors: + - {label: stat, symerror: 2.7} + - {label: sys, symerror: 4.6} + value: 68.7 + - errors: + - {label: stat, symerror: 2.5} + - {label: sys, symerror: 4.0} + value: 59.9 + - errors: + - {label: stat, symerror: 2.4} + - {label: sys, symerror: 3.8} + value: 56.7 + - errors: + - {label: stat, symerror: 1.9} + - {label: sys, symerror: 3.2} + value: 44.8 + - errors: + - {label: stat, symerror: 1.7} + - {label: sys, symerror: 2.4} + value: 35.1 +- header: {name: D(SIG)/DABS(ETARAP(LEPTON)), units: PB} + qualifiers: + - {name: ABS(ETARAP(JET)), value: < 2.5} + - {name: JETS, value: 'ANTI-KT, R=1'} + - {name: PT(JET), units: GEV, value: '> 25'} + - {name: PT(LEPTON), units: GEV, value: '> 35'} + - {name: RE, value: P P --> W+ < LEPTON+ NU > CHARMBAR X} + - {name: RE, value: P P --> W- < LEPTON- NUBAR > CHARM X} + - {name: SQRT(S), units: GeV, value: '7000.0'} + values: + - errors: + - {label: stat, symerror: 1.7} + - {label: sys, symerror: 3.2} + value: 52.3 + - errors: + - {label: stat, symerror: 1.6} + - {label: sys, symerror: 3.0} + value: 49.2 + - errors: + - {label: stat, symerror: 1.5} + - {label: sys, symerror: 2.7} + value: 45.5 + - errors: + - {label: stat, symerror: 1.2} + - {label: sys, symerror: 2.1} + value: 34.2 + - errors: + - {label: stat, symerror: 1.0} + - {label: sys, symerror: 1.7} + value: 26.6 +independent_variables: +- header: {name: ABS(ETARAP(LEPTON))} + values: + - {high: 0.35, low: 0.0} + - {high: 0.7, low: 0.35} + - {high: 1.1, low: 0.7} + - {high: 1.6, low: 1.1} + - {high: 2.1, low: 1.6} diff --git a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/rawdata/Table6.yaml b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/rawdata/Table6.yaml new file mode 100755 index 0000000000..5ea79d8856 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/rawdata/Table6.yaml @@ -0,0 +1,91 @@ +dependent_variables: +- header: {name: ''} + qualifiers: + - {name: ABS(ETARAP(JET)), value: < 2.5} + - {name: JETS, value: 'ANTI-KT, R=1'} + - {name: PT(JET), units: GEV, value: '> 25'} + - {name: PT(LEPTON), units: GEV, value: '> 25'} + - {name: RE, value: P P --> W+ < LEPTON+ NU > CHARMBAR X} + - {name: RE, value: P P --> W- < LEPTON- NUBAR > CHARM X} + - {name: SQRT(S), units: GeV, value: '7000.0'} + values: + - {value: 1.0} + - {value: 0.77} + - {value: 0.78} + - {value: 0.76} + - {value: 0.72} + - {value: 0.77} + - {value: 1.0} + - {value: 0.76} + - {value: 0.74} + - {value: 0.7} + - {value: 0.78} + - {value: 0.76} + - {value: 1.0} + - {value: 0.74} + - {value: 0.7} + - {value: 0.76} + - {value: 0.74} + - {value: 0.74} + - {value: 1.0} + - {value: 0.69} + - {value: 0.72} + - {value: 0.7} + - {value: 0.7} + - {value: 0.69} + - {value: 1.0} +independent_variables: +- header: {name: ABS(ETARAP(LEPTON))} + values: + - {high: 0.35, low: 0.0} + - {high: 0.7, low: 0.35} + - {high: 1.1, low: 0.7} + - {high: 1.6, low: 1.1} + - {high: 2.1, low: 1.6} + - {high: 0.35, low: 0.0} + - {high: 0.7, low: 0.35} + - {high: 1.1, low: 0.7} + - {high: 1.6, low: 1.1} + - {high: 2.1, low: 1.6} + - {high: 0.35, low: 0.0} + - {high: 0.7, low: 0.35} + - {high: 1.1, low: 0.7} + - {high: 1.6, low: 1.1} + - {high: 2.1, low: 1.6} + - {high: 0.35, low: 0.0} + - {high: 0.7, low: 0.35} + - {high: 1.1, low: 0.7} + - {high: 1.6, low: 1.1} + - {high: 2.1, low: 1.6} + - {high: 0.35, low: 0.0} + - {high: 0.7, low: 0.35} + - {high: 1.1, low: 0.7} + - {high: 1.6, low: 1.1} + - {high: 2.1, low: 1.6} +- header: {name: ABS(ETARAP(LEPTON))} + values: + - {high: 0.35, low: 0.0} + - {high: 0.35, low: 0.0} + - {high: 0.35, low: 0.0} + - {high: 0.35, low: 0.0} + - {high: 0.35, low: 0.0} + - {high: 0.7, low: 0.35} + - {high: 0.7, low: 0.35} + - {high: 0.7, low: 0.35} + - {high: 0.7, low: 0.35} + - {high: 0.7, low: 0.35} + - {high: 1.1, low: 0.7} + - {high: 1.1, low: 0.7} + - {high: 1.1, low: 0.7} + - {high: 1.1, low: 0.7} + - {high: 1.1, low: 0.7} + - {high: 1.6, low: 1.1} + - {high: 1.6, low: 1.1} + - {high: 1.6, low: 1.1} + - {high: 1.6, low: 1.1} + - {high: 1.6, low: 1.1} + - {high: 2.1, low: 1.6} + - {high: 2.1, low: 1.6} + - {high: 2.1, low: 1.6} + - {high: 2.1, low: 1.6} + - {high: 2.1, low: 1.6} diff --git a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/rawdata/Table9.yaml b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/rawdata/Table9.yaml new file mode 100755 index 0000000000..f260ee1f54 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/rawdata/Table9.yaml @@ -0,0 +1,69 @@ +dependent_variables: +- header: {name: SIG(W+ CHARMBAR) / SIG(W- CHARM)} + qualifiers: + - {name: ABS(ETARAP(JET)), value: < 2.5} + - {name: JETS, value: 'ANTI-KT, R=1'} + - {name: PT(JET), units: GEV, value: '> 25'} + - {name: PT(LEPTON), units: GEV, value: '> 25'} + - {name: RE(W+ CHARMBAR), value: P P --> W+ < LEPTON+ NU > CHARMBAR X} + - {name: RE(W- CHARM), value: P P --> W- < LEPTON- NUBAR > CHARM X} + - {name: SQRT(S), units: GeV, value: '7000.0'} + values: + - errors: + - {label: stat, symerror: 0.052} + - {label: sys, symerror: 0.005} + value: 1.013 + - errors: + - {label: stat, symerror: 0.053} + - {label: sys, symerror: 0.005} + value: 0.96 + - errors: + - {label: stat, symerror: 0.051} + - {label: sys, symerror: 0.008} + value: 0.897 + - errors: + - {label: stat, symerror: 0.061} + - {label: sys, symerror: 0.014} + value: 1.062 + - errors: + - {label: stat, symerror: 0.058} + - {label: sys, symerror: 0.016} + value: 0.776 +- header: {name: SIG(W+ CHARMBAR) / SIG(W- CHARM)} + qualifiers: + - {name: ABS(ETARAP(JET)), value: < 2.5} + - {name: JETS, value: 'ANTI-KT, R=1'} + - {name: PT(JET), units: GEV, value: '> 25'} + - {name: PT(LEPTON), units: GEV, value: '> 35'} + - {name: RE(W+ CHARMBAR), value: P P --> W+ < LEPTON+ NU > CHARMBAR X} + - {name: RE(W- CHARM), value: P P --> W- < LEPTON- NUBAR > CHARM X} + - {name: SQRT(S), units: GeV, value: '7000.0'} + values: + - errors: + - {label: stat, symerror: 0.041} + - {label: sys, symerror: 0.007} + value: 0.993 + - errors: + - {label: stat, symerror: 0.039} + - {label: sys, symerror: 0.007} + value: 0.977 + - errors: + - {label: stat, symerror: 0.04} + - {label: sys, symerror: 0.008} + value: 0.927 + - errors: + - {label: stat, symerror: 0.046} + - {label: sys, symerror: 0.01} + value: 0.948 + - errors: + - {label: stat, symerror: 0.05} + - {label: sys, symerror: 0.011} + value: 0.784 +independent_variables: +- header: {name: ABS(ETARAP(LEPTON))} + values: + - {high: 0.35, low: 0.0} + - {high: 0.7, low: 0.35} + - {high: 1.1, low: 0.7} + - {high: 1.6, low: 1.1} + - {high: 2.1, low: 1.6} diff --git a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/uncertainties_WPWM-RATIO.yaml b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/uncertainties_WPWM-RATIO.yaml new file mode 100644 index 0000000000..64740a4566 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/uncertainties_WPWM-RATIO.yaml @@ -0,0 +1,20 @@ +definitions: + stat_uncorr: + description: Statistical uncertainty + treatment: ADD + type: UNCORR + ART_LABEL: + description: Correlated systematic uncertainty + treatment: MULT + type: CORR +bins: +- stat_uncorr: 0.052 + art_corr: 0.005 +- stat_uncorr: 0.053 + art_corr: 0.005 +- stat_uncorr: 0.051 + art_corr: 0.008 +- stat_uncorr: 0.061 + art_corr: 0.014 +- stat_uncorr: 0.058 + art_corr: 0.016 diff --git a/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/uncertainties_WPWM-TOT.yaml b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/uncertainties_WPWM-TOT.yaml new file mode 100644 index 0000000000..19e422c314 --- /dev/null +++ b/nnpdf_data/nnpdf_data/commondata/CMS_WCHARM_7TEV/uncertainties_WPWM-TOT.yaml @@ -0,0 +1,56 @@ +definitions: + stat_uncorr: + description: Statistical uncertainty + treatment: ADD + type: UNCORR + art_corr_1: + description: Correlated systematic uncertainty 1 + treatment: ADD + type: CORR + art_corr_2: + description: Correlated systematic uncertainty 2 + treatment: ADD + type: CORR + art_corr_3: + description: Correlated systematic uncertainty 3 + treatment: ADD + type: CORR + art_corr_4: + description: Correlated systematic uncertainty 4 + treatment: ADD + type: CORR + art_corr_5: + description: Correlated systematic uncertainty 5 + treatment: ADD + type: CORR +bins: +- stat_uncorr: 2700.0 + art_corr_1: -4.30085960e+03 + art_corr_2: 1.48958525e+02 + art_corr_3: -1.49014534e+03 + art_corr_4: -3.33760715e+02 + art_corr_5: 5.55417584e+02 +- stat_uncorr: 2500.0 + art_corr_1: -3.61128670e+03 + art_corr_2: 1.46245103e+02 + art_corr_3: 1.37172210e+03 + art_corr_4: -3.95171209e+02 + art_corr_5: 9.48387525e+02 +- stat_uncorr: 2400.0 + art_corr_1: -3.42032648e+03 + art_corr_2: 1.58201894e+02 + art_corr_3: 2.94011007e+02 + art_corr_4: -5.54441014e+02 + art_corr_5: -1.52397231e+03 +- stat_uncorr: 1900.0 + art_corr_1: -2.77487954e+03 + art_corr_2: 3.53874801e+02 + art_corr_3: 1.26728625e+02 + art_corr_4: 1.53882437e+03 + art_corr_5: -1.75429767e+02 +- stat_uncorr: 1700.0 + art_corr_1: -1.94086828e+03 + art_corr_2: -1.38692811e+03 + art_corr_3: 5.04690564e+01 + art_corr_4: 2.51872601e+02 + art_corr_5: -5.89390783e+01