generateFlavCfromFlavB.py

Go to the documentation of this file.

#!/usr/bin/env python

from __future__ import print_function
from __future__ import absolute_import
import os
import sys
import itertools
from . import dataLoader
from . import checkBTagCalibrationConsistency as checker

import six

def generate_flav_c(loaded_data):
    flav_b_data = [e for e in loaded_data.entries if e.params.jetFlavor == 0]
    flav_b_data = sorted(flav_b_data, key=lambda e: e.params.operatingPoint)
    flav_b_data = sorted(flav_b_data, key=lambda e: e.params.measurementType)
    flav_b_data = sorted(flav_b_data, key=lambda e: e.params.etaMin)
    flav_b_data = sorted(flav_b_data, key=lambda e: e.params.ptMin)
    flav_b_data = sorted(flav_b_data, key=lambda e: e.params.discrMin)
    sys_groups = itertools.groupby(
        flav_b_data,
        key=lambda e: '%d, %s, %.02f, %.02f, %.02f' % (
            e.params.operatingPoint,
            e.params.measurementType,
            e.params.etaMin,
            e.params.ptMin,
            e.params.discrMin
        )
    )

    def gen_entry_dict(groups):
        for _, grp in groups:
            grp = list(grp)
            entries_by_sys = dict((e.params.sysType, e) for e in grp)
            assert len(grp) == len(entries_by_sys)  # every sysType is unique
            yield entries_by_sys
    sys_dicts = gen_entry_dict(sys_groups)

    def gen_flavb_csv_line(dicts):
        for d in dicts:
            central = d.pop('central')
            central.params.jetFlavor = 1
            yield central.makeCSVLine()
            for e in six.itervalues(d):
                e.params.jetFlavor = 1
                e.formula = '2*(%s)-(%s)' % (e.formula, central.formula)
                yield e.makeCSVLine()
    csv_lines = gen_flavb_csv_line(sys_dicts)

    return list(l for l in csv_lines)


def main():
    if len(sys.argv) < 3:
        print('Need input/output filenames as first/second arguments. Exit.')
        exit(-1)
    if os.path.exists(sys.argv[2]):
        print('Output file exists. Exit.')
        exit(-1)

    print('\nChecking input file consistency...')
    loaders = dataLoader.get_data(sys.argv[1])
    checks = checker.run_check_data(loaders, True, True, False)
    for data in loaders:
        typ = data.meas_type
        if 1 in data.flavs:
            print('FLAV_C already present in input file for %s. Exit.' % typ)
            exit(-1)
    if not any(0 in data.flavs for data in loaders):
        print('FLAV_B not found in input file. Exit.')
        exit(-1)


    print('\nGenerating new csv content...')
    new_csv_data = list(itertools.chain.from_iterable(
        l
        for d in loaders
        for l in generate_flav_c(d)
    ))

    with open(sys.argv[1]) as f:
        old_csv_data = f.readlines()

    with open(sys.argv[2], 'w') as f:
        f.writelines(old_csv_data)
        f.write('\n')
        f.writelines(new_csv_data)

    print('Done.')


if __name__ == '__main__':
    main()