Add calo efficiency tests using light_reco tracking configuration

Hlt/RecoConf/options/hlt2_light_reco_calo_efficiency.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.                                       #
+###############################################################################
+from Moore import options, run_reconstruction
+from RecoConf.standalone import standalone_hlt2_calo_efficiency
+
+options.evt_max = 1000
+options.ntuple_file = 'outputfile_calo_eff_gamma_light_reco.root'
+
+with standalone_hlt2_calo_efficiency.bind(light_reco=True):
+    run_reconstruction(options, standalone_hlt2_calo_efficiency)

Hlt/RecoConf/options/hlt2_reco_calo_selectivematching.py

@@ -16,5 +16,6 @@ options.evt_max = 1000
 options.ntuple_file = 'outputfile_calo_selectivematching.root'
 
 with make_photons_and_electrons.bind(selective_matching=True), \
-     make_electron_and_brem_match.bind(selective_matching=True):
+     make_electron_and_brem_match.bind(selective_matching=True), \
+     standalone_hlt2_calo_efficiency.bind(light_reco=True):
     run_reconstruction(options, standalone_hlt2_calo_efficiency)

Hlt/RecoConf/python/RecoConf/standalone.py

@@ -528,6 +528,7 @@ def standalone_hlt2_calo_resolution_pi0():
     return Reconstruction('hlt2_calo_resolution_pi0', data, reco_prefilters())
 
 
+@configurable
 def standalone_hlt2_calo_efficiency(light_reco=False):
     # get tracks and pvs, needed for calo reco
     hlt2_tracks = make_hlt2_tracks(light_reco=light_reco)
@@ -536,11 +537,31 @@ def standalone_hlt2_calo_efficiency(light_reco=False):
     pvs_v1 = RecVerticesV1FromV2(
         InputVerticesName=pvs,
         InputTracksName=hlt2_tracks["Velo"]["v1"]).OutputVerticesName
-    # get calo reco
 
-    tracks_for_calo = get_default_tracks_for_calo(
-        hlt2_tracks, light_reco=light_reco)
-    calo = make_calo(tracks_for_calo, pvs_v1)
+    # get tracks for calo
+    if light_reco:
+        out_track_types = get_default_out_track_types_for_light_reco()
+        trackSelections = {}
+        for track_type in out_track_types["Best"] + out_track_types[
+                "Unfitted4Calo"]:
+            trackSelections[track_type] = TracksToSelection(
+                InputLocation=hlt2_tracks[track_type]["v1"]).OutputLocation
+        all_best_tracks = {}
+        all_best_tracks["v1"] = TrackSelectionMerger(InputLocations=[
+            trackSelections[track_type]
+            for track_type in out_track_types["Best"]
+        ]).OutputLocation
+        tracks4calo = [all_best_tracks["v1"]] + [
+            trackSelections[key] for key in out_track_types["Unfitted4Calo"]
+        ]
+        alltracks4calo = TrackSelectionMerger(
+            InputLocations=tracks4calo).OutputLocation
+    else:
+        best_tracks = hlt2_tracks["Best"]
+        alltracks4calo = best_tracks["v1"]
+
+    # get calo reco
+    calo = make_calo(alltracks4calo, pvs_v1)
 
     # add also 'ElectronMatch' and 'BremMatch', in case of selective track matching test
     data = [calo["electronMatch"], calo["bremMatch"]]
@@ -551,7 +572,7 @@ def standalone_hlt2_calo_efficiency(light_reco=False):
 
     # full calo reco with tracking input
     data += check_calo_hypo_efficiency(
-        calo["digitsEcal"], calo["ecalClusters"], tracks_for_calo, pvs_v1)
+        calo["digitsEcal"], calo["ecalClusters"], alltracks4calo, pvs_v1)
 
     return Reconstruction('hlt2_calo_efficiency', data)
 

Hlt/RecoConf/tests/qmtest/hlt2_light_reco_calo_efficiency.qmt

+<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE extension  PUBLIC '-//QM/2.3/Extension//EN'  'http://www.codesourcery.com/qm/dtds/2.3/-//qm/2.3/extension//en.dtd'>
+<!--
+    (c) Copyright 2000-2018 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.
+-->
+<!--
+Check MC efficiency of ECAL reconstruction
+-->
+<extension class="GaudiTest.GaudiExeTest" kind="test">
+<argument name="program"><text>gaudirun.py</text></argument>
+<argument name="args"><set>
+  <text>$MOOREROOT/options/force_functor_cache.py</text>
+  <text>$MOOREROOT/tests/options/xdigi_input_and_conds_calo_res_gamma.py</text>
+  <text>$RECOCONFROOT/options/hlt2_light_reco_calo_efficiency.py</text>
+  <text>--output=hlt2_light_reco_calo_efficiency.opts.py</text>
+  <text>--all-opt</text>
+</set></argument>
+<argument name="use_temp_dir"><enumeral>true</enumeral></argument>
+<argument name="reference"><text>../refs/hlt2_light_reco_calo_efficiency.ref</text></argument>
+<argument name="validator"><text>
+
+from Moore.qmtest.exclusions import ref_preprocessor
+validateWithReference(preproc = ref_preprocessor)
+
+</text></argument>
+</extension>
\ No newline at end of file