QEE: Migrate the single muon HLT2 lines to ThOr

FYI @wbarter @hyin @lugrazet

Shifters please ignore

Now that ThOr has a SUMCONE functor, the HLT2 single, high-pt muon lines can be migrated to ThOr. At the moment, this functor does not remove the signal track from the cone, so the signal is manually removed by removing the track's pT from the sum. When the signal removal is added to SUMCONE, this can be removed. I will open an issue to make sure this is tracked.

Very small change in efficiency firing on 1k Z \to \mu\mu events on the CanRecoChildren denominator from LoKi (top) to ThOr (bottom):

Line:	 Hlt2SingleHighPtMuonIsoDecision               	 Efficiency:	 0.970 +/- 0.015

Line:	 Hlt2SingleHighPtMuonIsoDecision               	 Efficiency:	 0.978 +/- 0.013

Rates now (measured on 50k HLT1-filtered min bias):

Line:	 Hlt2SingleHighPtMuonIsoDecision     	 Incl: 0.13999 +/- 0.0529 kHz,	Excl: 0.12000 +/- 0.0489 kHz
Line:	 Hlt2SingleHighPtMuonDecision        	 Incl: 0.04 +/- 0.0282 kHz,	Excl: 0.02 +/- 0.0199 kHz
Line:	 Hlt2SingleHighPtMuonNoMuIDDecision  	 Incl: 0.24000 +/- 0.0692 kHz,	Excl: 0.24000 +/- 0.0692 kHz
Line:	 Hlt2SingleHighPtMuonPrescaleDecision	 Incl: 0.04 +/- 0.0282 kHz,	Excl: 0.04 +/- 0.0282 kHz

Previous LoKi-based implementation:

Line:	 Hlt2SingleHighPtMuonIsoDecision         Incl: 0.12 +/- 0.04 kHz, Excl: 0.08 +/- 0.03 kHz
Line:	 Hlt2SingleHighPtMuonDecision            Incl: 0.04 +/- 0.02 kHz, Excl: 0.0  +/- 0.0  kHz
Line:	 Hlt2SingleHighPtMuonNoMuIDDecision      Incl: 0.36 +/- 0.08 kHz, Excl: 0.22 +/- 0.06 kHz
Line:	 Hlt2SingleHighPtMuonPrescaleDecision    Incl: 0.02 +/- 0.01 kHz, Excl: 0.0  +/- 0.0  kHz

So compatible within uncertainties.

Hlt/Hlt2Conf/python/Hlt2Conf/lines/qee/single_high_pt_muon.py

@@ -14,12 +14,15 @@ Definition of the HLT2 single high-PT muon line.
 
 from GaudiKernel.SystemOfUnits import GeV
 from PyConf import configurable
+from PyConf.Algorithms import WeightedRelTableAlg
 from Moore.config import register_line_builder
 from Moore.lines import Hlt2Line
-from Hlt2Conf.algorithms import require_all, ParticleFilter
+from Hlt2Conf.algorithms_thor import require_all, ParticleFilter
 from Hlt2Conf.standard_particles import make_ismuon_long_muon, make_long_muons
 from RecoConf.reconstruction_objects import upfront_reconstruction
 from Hlt2Conf.standard_jets import make_onlytrack_particleflow
+import Functors as F
+from Functors.math import in_range
 
 all_lines = {}
 
@@ -29,9 +32,8 @@ def make_highpt_muons(name='HighPtMuonMaker',
                       muons=make_ismuon_long_muon,
                       min_muon_pt=15. * GeV):
 
-    code = require_all('PT > {min_pt}').format(min_pt=min_muon_pt)
-
-    return ParticleFilter(particles=muons(), name=name, Code=code)
+    code = require_all(F.PT > min_muon_pt)
+    return ParticleFilter(muons(), F.FILTER(code), name=name)
 
 
 @configurable
@@ -39,17 +41,22 @@ def make_highpt_isolated_muons(name='HighPtIsolatedMuonMaker',
                                muons=make_ismuon_long_muon,
                                min_muon_pt=12.5 * GeV,
                                max_cone_pt=10.0 * GeV,
-                               pflow_output="HltParticleFlow/Particles"):
+                               pflow_output=make_onlytrack_particleflow()):
+
+    ftAlg = WeightedRelTableAlg(
+        ReferenceParticles=muons(),
+        InputCandidates=pflow_output,
+        Cut=(F.DR2 < (0.5**2)))
+
+    ftAlg_Rels = ftAlg.OutputRelations
 
-    code = require_all(
-        'PT > {min_muon_pt}',
-        'SUMCONE(0.5, PT, TESDATA("{pflow_location}")) < {max_cone_pt}',
-    ).format(
-        min_muon_pt=min_muon_pt,
-        max_cone_pt=max_cone_pt,
-        pflow_location=pflow_output)
+    # FIXME Signal is (for now) manually removed from cone with - F.PT.
+    code = require_all(F.PT > min_muon_pt, (in_range(
+        0,
+        F.VALUE_OR(0) @ F.SUMCONE(Functor=F.PT, Relations=ftAlg_Rels) - F.PT,
+        max_cone_pt)))
 
-    return ParticleFilter(particles=muons(), name=name, Code=code)
+    return ParticleFilter(muons(), F.FILTER(code), name=name)
 
 
 @register_line_builder(all_lines)
@@ -75,7 +82,6 @@ def single_muon_highpt_prescale_line(name='Hlt2SingleHighPtMuonPrescale',
                                      persistreco=True):
     """High PT single muon line with lower pT threshold, prescaled to reduce the rate."""
 
-    # pflow algo populates the HltParticleFlow/Particles path for SUMCONE
     high_pt_muons = make_highpt_muons(min_muon_pt=10. * GeV)
 
     return Hlt2Line(
@@ -92,10 +98,9 @@ def single_muon_highpt_iso_line(name='Hlt2SingleHighPtMuonIso',
                                 persistreco=True):
     """High PT single muon line with lower pT threshold, and isolation to keep the rate down."""
 
-    # pflow algo populates the HltParticleFlow/Particles path for SUMCONE
     pflow = make_onlytrack_particleflow()
     high_pt_isolated_muons = make_highpt_isolated_muons(
-        min_muon_pt=12.5 * GeV, pflow_output=pflow.location)
+        min_muon_pt=12.5 * GeV, pflow_output=pflow)
 
     return Hlt2Line(
         name=name,