run_collector_path.py

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
 
 
 
import basf2
from basf2 import conditions as b2conditions
from basf2 import create_path, process
from basf2 import B2FATAL, B2INFO, B2WARNING
 
# The backend provides the input data files to the job and you can get the list of
# files from this function
from caf.backends import get_input_data
 
 
def run_collectors():
    """
    Runs the CAF collector in a general way by passing serialised basf2 paths, input file
    locations, and local databases from previous calibrations.
    """
    from basf2.pickle_path import get_path_from_file
    import glob
    import os
    import json
 
    # We pass in basf2 and CAF options via the config.json file
    with open("collector_config.json", "r") as config_file:
        config = json.load(config_file)
 
    # Create the database chain to use the necessary central DB global tags and local DBs if they are requested
    # We deliberately override the normal database ordering because we don't want input files GTs to affect
    # the processing. Only explicit GTs and intermediate local DBs made by the CAF should be added here.
    b2conditions.reset()
    b2conditions.override_globaltags()
 
    for db_type, database in config["database_chain"]:
        if db_type == "local":
            B2INFO(f"Adding Local Database {database[0]} to head of chain of local databases.")
            b2conditions.prepend_testing_payloads(database[0])
        else:
            B2INFO(f"Using Global Tag {database}.")
            b2conditions.prepend_globaltag(database)
 
    # create a path with all modules needed before calibration path is run.
    collector_path = create_path()
 
    # Grab all the pickled path files. Should at least be one for the collector. And optionally
    # one for the modules to run before the collector.
    pre_collector_path_name = "pre_collector.path"
    pickle_paths = glob.glob("./*.path")
 
    # Remove the pre-collector path from the overall list
    pickle_paths = filter(lambda x: pre_collector_path_name not in x, pickle_paths)
 
    # Check for and add a path to run before the collector
    if os.path.exists(pre_collector_path_name):
        collector_path.add_path(get_path_from_file(pre_collector_path_name))
    else:
        B2INFO("Couldn't find any pickle files for pre-collector setup")
 
    # Check for any other paths to run (should only be collectors) and append them
    if pickle_paths:
        for pickle_path in pickle_paths:
            collector_path.add_path(get_path_from_file(pickle_path))
    else:
        B2FATAL("Couldn't find any pickle files for collector path!")
 
    # Grab the input data. Can't continue without some
    input_data = get_input_data()
 
    # Now we need to create a path that definitely has RootInput as a module.
    main = create_path()
    # Use this utility wrapper to check for RootInput or SeqRootInput and change params if necessary
    from caf.utils import PathExtras
    pe = PathExtras(collector_path)
    if "RootInput" in pe:
        root_input_mod = collector_path.modules()[pe.index("RootInput")]
        root_input_mod.param("inputFileNames", input_data)
        input_entry_sequences = get_entry_sequences(root_input_mod)
        # If entrySquences is set, then we should make sure that it has the same length as the number of input files (if possible)
        # Other parameters can be taken from the input directly, but this one is sensitive to the input file number.
        if len(input_entry_sequences) == 1 and len(input_data) > 1:
            root_input_mod.param("entrySequences", len(input_data)*input_entry_sequences)
    elif "SeqRootInput" in pe:
        root_input_mod = collector_path.modules()[pe.index("SeqRootInput")]
        root_input_mod.param("inputFileNames", input_data)
    else:
        main.add_module("RootInput", inputFileNames=input_data)
 
    if "HistoManager" not in pe:
        main.add_module("HistoManager", histoFileName="CollectorOutput.root")
 
    main.add_path(collector_path)
    basf2.print_path(main)
    process(main)
    print(basf2.statistics)
 
 
def get_entry_sequences(root_input_module):
    """
    Finds the value of entrySequences (list) from the RootInput module.
    """
    for module_param_info in root_input_module.available_params():
        if module_param_info.name == "entrySequences":
            entry_sequences = module_param_info.values
            break
    else:
        B2WARNING(("entrySequences ModuleParamInfo couldn't be found! "
                   "Setting to empty list but this indicates that this script is out of date!"))
        entry_sequences = []
    return entry_sequences
 
 
if __name__ == "__main__":
    run_collectors()