First iteration of rate and efficiency testing scripts for QEE upgrade

HltEfficiencyChecker/options/hlt2_QEE_eff_example.py

+###############################################################################
+# (c) Copyright 2019 CERN for the benefit of the LHCb Collaboration           #
+#                                                                             #
+# This software is distributed under the terms of the GNU General Public      #
+# Licence version 3 (GPL Version 3), copied verbatim in the file "COPYING".   #
+#                                                                             #
+# In applying this licence, CERN does not waive the privileges and immunities #
+# granted to it by virtue of its status as an Intergovernmental Organization  #
+# or submit itself to any jurisdiction.                                       #
+###############################################################################
+""" Example options file for testing QEE Hlt2 lines on signal MC with HltEfficiencyChecker
+    Run like any other options file:
+
+    ./MooreAnalysis/run gaudirun.py MooreAnalysis/HltEfficiencyChecker/options/hlt2_QEE_eff_example.py
+"""
+
+from Moore import options
+from HltEfficiencyChecker.config import run_moore_with_tuples
+from RecoConf.reco_objects_from_file import stateProvider_with_simplified_geom
+from Hlt2Conf.lines.QEE.SingleHighPtMuon import all_lines
+
+def make_lines():
+    return [builder() for builder in all_lines.values()]
+
+decay = (
+    "${W}[W+ ==> ${mu}mu+ ${nu}nu_mu ]CC"
+)
+
+input_files = [
+    # W+jet 42311011 upgrade sample. See https://its.cern.ch/jira/browse/LHCBGAUSS-1837
+    # Not HLT1 filtered.
+    "root://x509up_u1000@eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/Upgrade/XDST/00100893/0000/00100893_00000010_2.xdst",
+    "root://x509up_u1000@eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/Upgrade/XDST/00100893/0000/00100893_00000011_2.xdst",
+    "root://x509up_u1000@eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/Upgrade/XDST/00100893/0000/00100893_00000012_2.xdst",
+    "root://x509up_u1000@eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/Upgrade/XDST/00100893/0000/00100893_00000013_2.xdst",
+    "root://x509up_u1000@eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/Upgrade/XDST/00100893/0000/00100893_00000014_2.xdst",
+    "root://x509up_u1000@eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/Upgrade/XDST/00100893/0000/00100893_00000015_2.xdst",
+    "root://x509up_u1000@eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/Upgrade/XDST/00100893/0000/00100893_00000016_2.xdst",
+    "root://x509up_u1000@eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/Upgrade/XDST/00100893/0000/00100893_00000017_2.xdst",
+    "root://x509up_u1000@eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/Upgrade/XDST/00100893/0000/00100893_00000018_2.xdst",
+    "root://x509up_u1000@eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/Upgrade/XDST/00100893/0000/00100893_00000020_2.xdst"
+]
+
+options.input_type = 'ROOT'
+options.input_raw_format = 4.3
+# When running from Upgrade MC, must use the post-juggling locations of the raw
+# event
+options.input_files = input_files
+options.data_type = 'Upgrade'
+options.dddb_tag = 'dddb-20190223'
+options.conddb_tag = 'sim-20180530-vc-mu100'
+options.simulation = True
+options.scheduler_legacy_mode = True
+options.evt_max = 1000
+options.lines_maker = make_lines
+
+# Write the output file
+options.ntuple_file = "eff_ntuple.root"
+
+# Uncomment the following to increase the output verbosity
+from Gaudi.Configuration import DEBUG
+#options.output_level = DEBUG
+
+run_moore_with_tuples(
+    options, decay, public_tools=[stateProvider_with_simplified_geom()])
\ No newline at end of file

HltEfficiencyChecker/options/hlt2_QEE_rate_example.yaml

+# HltEfficiencyChecker "wizard" example for Hlt2 rates of QEE rates
+ntuple_path: &NTUPLE rate_ntuple.root
+job:
+    trigger_level: 2
+    evt_max: 10000
+    testfiledb_key: upgrade_minbias_hlt1_filtered
+    input_raw_format: 4.3
+    lines_from: Hlt2Conf.lines.QEE.SingleHighPtMuon
+    options:
+        - $HLTEFFICIENCYCHECKERROOT/options/options_template.py.jinja  # first rendered with jinja2
+analysis:
+    script: $HLTEFFICIENCYCHECKERROOT/scripts/hlt_calculate_rates.py
+    args:
+        null: *NTUPLE
+        level: Hlt2
+        json: Hlt2_rates.json

HltEfficiencyChecker/options/hlt2_QEE_rates_nothlt1filtered.yaml

+# HltEfficiencyChecker "wizard" example for Hlt2 rates of QEE rates, not hlt1 filtered
+ntuple_path: &NTUPLE rate_ntuple.root
+job:
+    trigger_level: 2
+    evt_max: 10000
+    testfiledb_key: upgrade-magdown-sim09c-up02-reco-up01-minbias-ldst
+    input_raw_format: 4.3
+    lines_from: Hlt2Conf.lines.QEE.SingleHighPtMuon
+    #ft_decoding_version: 4
+    options:
+        - $HLTEFFICIENCYCHECKERROOT/options/options_template.py.jinja  # first rendered with jinja2
+analysis:
+    script: $HLTEFFICIENCYCHECKERROOT/scripts/hlt_calculate_rates.py
+    args:
+        null: *NTUPLE
+        level: Hlt2
+        json: Hlt2_rates.json
+        hlt1-output-rate: 30000