Merge branch 'ift-line-ion' into 'master'

HLT2 baseline set-up and tests for ion data

See merge request !2117

Hlt/RecoConf/options/hlt2_lead_lead_light_reco_pr_kf.py → Hlt/Hlt1Conf/options/allen_hlt1_PbPb.py

 ###############################################################################
-# (c) Copyright 2019 CERN for the benefit of the LHCb Collaboration           #
+# (c) Copyright 2000-2022 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".   #
@@ -8,13 +8,8 @@
 # granted to it by virtue of its status as an Intergovernmental Organization  #
 # or submit itself to any jurisdiction.                                       #
 ###############################################################################
+from Moore import (options, run_allen)
+from RecoConf.hlt1_allen import (allen_gaudi_config)
 
-from Moore import options, run_reconstruction
-from RecoConf.standalone import standalone_hlt2_light_reco
-from RecoConf.event_filters import require_gec
-from RecoConf.hlt2_tracking import make_hlt2_tracks
-
-# A higher gec cut than 40000 leads to unreasonable processing times for the nightlies.
-with require_gec.bind(cut=40000), \
-     make_hlt2_tracks.bind(use_pr_kf=True):
-    run_reconstruction(options, standalone_hlt2_light_reco)
+with allen_gaudi_config.bind(sequence="hlt1_PbPb_default"):
+    run_allen(options)

Hlt/Hlt1Conf/options/allen_hlt1_PbPb_minimal.py

+###############################################################################
+# (c) Copyright 2000-2022 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_allen)
+from RecoConf.hlt1_allen import (allen_gaudi_config)
+
+with allen_gaudi_config.bind(sequence="hlt1_PbPb_minimal"):
+    run_allen(options)

Hlt/Hlt2Conf/options/hlt2_PbPb_default.py

+###############################################################################
+# (c) Copyright 2021 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_moore
+from RecoConf.reconstruction_objects import reconstruction
+from RecoConf.hlt1_muonid import make_muon_hits
+from RecoConf.hlt1_tracking import (make_reco_pvs, make_PatPV3DFuture_pvs)
+from RecoConf.hlt2_tracking import make_TrackBestTrackCreator_tracks
+from RecoConf.hlt2_global_reco import make_fastest_reconstruction
+from RecoConf.hlt2_global_reco import reconstruction as reco
+from RecoConf.event_filters import require_gec
+from Hlt2Conf.lines.iftPbPb import all_lines as ift_lines
+from PyConf.Algorithms import PrHybridSeeding
+
+from Configurables import LHCb__Det__LbDD4hep__DD4hepSvc as DD4hepSvc
+
+dd4hepSvc = DD4hepSvc()
+dd4hepSvc.DetectorList = [
+    '/world', 'VP', 'FT', 'Magnet', 'Rich1', 'Rich2', 'Ecal', 'Hcal', 'Muon'
+]  #no UT
+dd4hepSvc.ConditionsLocation = 'git:/cvmfs/lhcb.cern.ch/lib/lhcb/git-conddb/lhcb-conditions-database.git'
+
+options.output_file = 'hlt2_ion_data.dst'
+options.output_type = 'ROOT'
+options.output_manifest_file = "hlt2_ion_data.tck.json"
+options.histo_file = "hlt2_ion_data.root"
+public_tools = []
+
+
+def make_lines():
+    lines = []
+    for builder in ift_lines.values():
+        has_gec = False
+        node = builder().node
+        for child in node.children:
+            for grand_child in child.children:
+                if "PrGECFilter" in grand_child.name:
+                    has_gec = True
+        if has_gec:
+            lines.append(builder)
+        else:
+            print("WARNING: line %s has not GEC and will not be processed" %
+                  node.name)
+
+    return [builder() for builder in lines]
+
+with reconstruction.bind(from_file=False),\
+    require_gec.bind(cut=30000,skipUT=True), \
+    make_fastest_reconstruction.bind(skipUT=False),\
+    PrHybridSeeding.bind(RemoveBeamHole=True, RemoveClones_forLead=True), \
+    reco.bind(make_reconstruction=make_fastest_reconstruction),\
+    make_reco_pvs.bind(make_pvs_from_velo_tracks=make_PatPV3DFuture_pvs),\
+    make_TrackBestTrackCreator_tracks.bind(max_ghost_prob=0.8),\
+    make_muon_hits.bind(geometry_version=3):
+    config = run_moore(options, make_lines, public_tools)

Hlt/Hlt2Conf/options/hlt2_PbPb_default_retina.py

+###############################################################################
+# (c) Copyright 2021 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_moore
+from PyConf.Algorithms import VeloRetinaClusterTrackingSIMD
+from RecoConf.reconstruction_objects import reconstruction
+from RecoConf.hlt1_muonid import make_muon_hits
+from RecoConf.hlt1_tracking import (
+    make_VeloClusterTrackingSIMD, make_RetinaClusters, make_reco_pvs,
+    make_PatPV3DFuture_pvs, get_global_measurement_provider)
+from RecoConf.hlt2_tracking import (get_UpgradeGhostId_tool_no_UT,
+                                    make_TrackBestTrackCreator_tracks)
+from RecoConf.hlt2_global_reco import make_fastest_reconstruction
+from RecoConf.hlt2_global_reco import reconstruction as reco
+from RecoConf.event_filters import require_gec
+from Hlt2Conf.lines.iftPbPb import all_lines as ift_lines
+from PyConf.Algorithms import PrHybridSeeding
+
+from Configurables import LHCb__Det__LbDD4hep__DD4hepSvc as DD4hepSvc
+
+dd4hepSvc = DD4hepSvc()
+dd4hepSvc.DetectorList = [
+    '/world', 'VP', 'FT', 'Magnet', 'Rich1', 'Rich2', 'Ecal', 'Hcal', 'Muon'
+]  #no UT
+dd4hepSvc.ConditionsLocation = 'git:/cvmfs/lhcb.cern.ch/lib/lhcb/git-conddb/lhcb-conditions-database.git'
+
+options.output_file = 'hlt2_ion_data.dst'
+options.output_type = 'ROOT'
+options.output_manifest_file = "hlt2_ion_data.tck.json"
+options.histo_file = "hlt2_ion_data.root"
+public_tools = []
+
+
+def make_lines():
+    lines = []
+    for builder in ift_lines.values():
+        has_gec = False
+        node = builder().node
+        for child in node.children:
+            for grand_child in child.children:
+                if "PrGECFilter" in grand_child.name:
+                    has_gec = True
+        if has_gec:
+            lines.append(builder)
+        else:
+            print("WARNING: line %s has not GEC and will not be processed" %
+                  node.name)
+
+    return [builder() for builder in lines]
+
+with reconstruction.bind(from_file=False),\
+    require_gec.bind(cut=30000,skipUT=True), \
+    make_RetinaClusters.bind(detector="VPRetinaCluster"),\
+    make_fastest_reconstruction.bind(skipUT=False),\
+    PrHybridSeeding.bind(RemoveBeamHole=True, RemoveClones_forLead=True),\
+    reco.bind(make_reconstruction=make_fastest_reconstruction),\
+    make_VeloClusterTrackingSIMD.bind(algorithm=VeloRetinaClusterTrackingSIMD, SkipForward=4),\
+    make_reco_pvs.bind(make_pvs_from_velo_tracks=make_PatPV3DFuture_pvs),\
+    get_UpgradeGhostId_tool_no_UT.bind(velo_hits=make_RetinaClusters),\
+    get_global_measurement_provider.bind(velo_hits=make_RetinaClusters),\
+    make_TrackBestTrackCreator_tracks.bind(max_ghost_prob=0.8),\
+    make_muon_hits.bind(geometry_version=3):
+    config = run_moore(options, make_lines, public_tools)

Hlt/Hlt2Conf/options/hlt2_PbPb_without_UT.py

+###############################################################################
+# (c) Copyright 2021 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_moore
+from RecoConf.reconstruction_objects import reconstruction
+from RecoConf.hlt1_muonid import make_muon_hits
+from RecoConf.hlt1_tracking import (make_reco_pvs, make_PatPV3DFuture_pvs)
+from RecoConf.hlt2_tracking import make_TrackBestTrackCreator_tracks
+from RecoConf.hlt2_global_reco import make_fastest_reconstruction
+from RecoConf.hlt2_global_reco import reconstruction as reco
+from RecoConf.event_filters import require_gec
+from Hlt2Conf.lines.iftPbPb import all_lines as ift_lines
+from PyConf.Algorithms import PrHybridSeeding
+
+from Configurables import LHCb__Det__LbDD4hep__DD4hepSvc as DD4hepSvc
+
+dd4hepSvc = DD4hepSvc()
+dd4hepSvc.DetectorList = [
+    '/world', 'VP', 'FT', 'Magnet', 'Rich1', 'Rich2', 'Ecal', 'Hcal', 'Muon'
+]  #no UT
+dd4hepSvc.ConditionsLocation = 'git:/cvmfs/lhcb.cern.ch/lib/lhcb/git-conddb/lhcb-conditions-database.git'
+
+options.output_file = 'hlt2_ion_data.dst'
+options.output_type = 'ROOT'
+options.output_manifest_file = "hlt2_ion_data.tck.json"
+options.histo_file = "hlt2_ion_data.root"
+public_tools = []
+
+
+def make_lines():
+    lines = []
+    for builder in ift_lines.values():
+        has_gec = False
+        node = builder().node
+        for child in node.children:
+            for grand_child in child.children:
+                if "PrGECFilter" in grand_child.name:
+                    has_gec = True
+        if has_gec:
+            lines.append(builder)
+        else:
+            print("WARNING: line %s has not GEC and will not be processed" %
+                  node.name)
+
+    return [builder() for builder in lines]
+
+with reconstruction.bind(from_file=False),\
+    require_gec.bind(cut=30000, skipUT=True), \
+    make_fastest_reconstruction.bind(skipUT=True),\
+    PrHybridSeeding.bind(RemoveBeamHole=True, RemoveClones_forLead=True),\
+    reco.bind(make_reconstruction=make_fastest_reconstruction),\
+    make_reco_pvs.bind(make_pvs_from_velo_tracks=make_PatPV3DFuture_pvs),\
+    make_TrackBestTrackCreator_tracks.bind(max_ghost_prob=0.8),\
+    make_muon_hits.bind(geometry_version=3):
+    config = run_moore(options, make_lines, public_tools)

Hlt/Hlt2Conf/options/hlt2_PbPb_without_UT_retina.py

+###############################################################################
+# (c) Copyright 2021 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_moore
+from PyConf.Algorithms import VeloRetinaClusterTrackingSIMD
+from RecoConf.reconstruction_objects import reconstruction
+from RecoConf.hlt1_muonid import make_muon_hits
+from RecoConf.hlt1_tracking import (
+    make_VeloClusterTrackingSIMD, make_RetinaClusters, make_reco_pvs,
+    make_PatPV3DFuture_pvs, get_global_measurement_provider)
+from RecoConf.hlt2_tracking import (get_UpgradeGhostId_tool_no_UT,
+                                    make_TrackBestTrackCreator_tracks)
+from RecoConf.hlt2_global_reco import make_fastest_reconstruction
+from RecoConf.hlt2_global_reco import reconstruction as reco
+from RecoConf.event_filters import require_gec
+from Hlt2Conf.lines.iftPbPb import all_lines as ift_lines
+from PyConf.Algorithms import PrHybridSeeding
+
+from Configurables import LHCb__Det__LbDD4hep__DD4hepSvc as DD4hepSvc
+
+dd4hepSvc = DD4hepSvc()
+dd4hepSvc.DetectorList = [
+    '/world', 'VP', 'FT', 'Magnet', 'Rich1', 'Rich2', 'Ecal', 'Hcal', 'Muon'
+]  #no UT
+dd4hepSvc.ConditionsLocation = 'git:/cvmfs/lhcb.cern.ch/lib/lhcb/git-conddb/lhcb-conditions-database.git'
+
+options.output_type = 'ROOT'
+public_tools = []
+
+
+def make_lines():
+    lines = []
+    for builder in ift_lines.values():
+        has_gec = False
+        node = builder().node
+        for child in node.children:
+            for grand_child in child.children:
+                if "PrGECFilter" in grand_child.name:
+                    has_gec = True
+        if has_gec:
+            lines.append(builder)
+        else:
+            print("WARNING: line %s has not GEC and will not be processed" %
+                  node.name)
+
+    return [builder() for builder in lines]
+
+with reconstruction.bind(from_file=False),\
+    require_gec.bind(cut=30000, skipUT=True), \
+    make_RetinaClusters.bind(detector="VPRetinaCluster"),\
+    make_fastest_reconstruction.bind(skipUT=True),\
+    PrHybridSeeding.bind(RemoveBeamHole=True, RemoveClones_forLead=True), \
+    reco.bind(make_reconstruction=make_fastest_reconstruction),\
+    make_VeloClusterTrackingSIMD.bind(algorithm=VeloRetinaClusterTrackingSIMD, SkipForward=4),\
+    make_reco_pvs.bind(make_pvs_from_velo_tracks=make_PatPV3DFuture_pvs),\
+    get_UpgradeGhostId_tool_no_UT.bind(velo_hits=make_RetinaClusters),\
+    get_global_measurement_provider.bind(velo_hits=make_RetinaClusters),\
+    make_TrackBestTrackCreator_tracks.bind(max_ghost_prob=0.8),\
+    make_muon_hits.bind(geometry_version=3):
+    config = run_moore(options, make_lines, public_tools)

Hlt/Hlt2Conf/python/Hlt2Conf/lines/ift/builders/__init__.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.                                       #
+###############################################################################
+
+##init

Hlt/RecoConf/options/hlt2_lead_lead_light_reco_pr_kf_gec_50000.py → Hlt/Hlt2Conf/python/Hlt2Conf/lines/iftPbPb/__init__.py

-###############################################################################
 # (c) Copyright 2019 CERN for the benefit of the LHCb Collaboration           #
 #                                                                             #
 # This software is distributed under the terms of the GNU General Public      #
@@ -8,13 +7,16 @@
 # granted to it by virtue of its status as an Intergovernmental Organization  #
 # or submit itself to any jurisdiction.                                       #
 ###############################################################################
+"""
+Submodule that defines all the ift HLT2 lines
+"""
+
+from . import hlt2_iftPbPb
 
-from Moore import options, run_reconstruction
-from RecoConf.standalone import standalone_hlt2_light_reco_ion
-from RecoConf.hlt2_tracking import make_hlt2_tracks_ion
-from RecoConf.event_filters import require_gec
+all_lines = {}
+all_lines.update(hlt2_iftPbPb.all_lines)
 
+from . import spruce_iftPbPb
 
-with require_gec.bind(cut=50000), \
-     make_hlt2_tracks_ion.bind( fast_reco=True, use_pr_kf=True, post_fit_selection='(TrGHOSTPROB<0.9999)'):
-    run_reconstruction(options, standalone_hlt2_light_reco_ion)
+sprucing_lines = {}
+sprucing_lines.update(spruce_iftPbPb.sprucing_lines)

Hlt/Hlt2Conf/python/Hlt2Conf/lines/iftPbPb/builders/charm_builders.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.                                       #
+###############################################################################
+import Functors as F
+from Functors.math import in_range
+from GaudiKernel.SystemOfUnits import GeV, MeV, mm
+
+from PyConf import configurable
+from RecoConf.reconstruction_objects import make_pvs
+from Hlt2Conf.algorithms_thor import (
+    ParticleFilter,
+    ParticleCombiner,
+)
+from Hlt2Conf.standard_particles import (make_has_rich_long_pions,
+                                         make_has_rich_long_kaons,
+                                         make_has_rich_long_protons)
+
+all_lines = {}
+
+################################################################################
+# Prompt lines
+################################################################################
+
+
+@configurable
+def make_charm_pions(bpvipchi2=9.0,
+                     maxPIDKaon=None,
+                     min_pt=250 * MeV,
+                     min_p=3000.0 * MeV,
+                     max_trchi2dof=3):
+    code = F.require_all(F.PT > min_pt, F.P > min_p, F.CHI2DOF < max_trchi2dof,
+                         F.BPVIPCHI2(make_pvs()) > bpvipchi2)
+    if maxPIDKaon is not None:
+        code &= (F.PID_K < maxPIDKaon)
+    return ParticleFilter(make_has_rich_long_pions(), F.FILTER(code))
+
+
+@configurable
+def make_charm_kaons(bpvipchi2=9.0,
+                     minPIDKaon=None,
+                     min_pt=250 * MeV,
+                     min_p=3000.0 * MeV,
+                     max_trchi2dof=3):
+    code = F.require_all(F.PT > min_pt, F.P > min_p, F.CHI2DOF < max_trchi2dof,
+                         F.BPVIPCHI2(make_pvs()) > bpvipchi2)
+    if minPIDKaon is not None:
+        code &= (F.PID_K > minPIDKaon)
+    return ParticleFilter(make_has_rich_long_kaons(), F.FILTER(code))
+
+
+@configurable
+def make_charm_protons_loose(bpvipchi2=9.0,
+                             minPIDProton=None,
+                             min_pt=250 * MeV,
+                             min_p=3000.0 * MeV,
+                             max_trchi2dof=3):
+    code = F.require_all(F.PT > min_pt, F.P > min_p, F.CHI2DOF < max_trchi2dof,
+                         F.BPVIPCHI2(make_pvs()) > bpvipchi2)
+    if minPIDProton is not None:
+        code &= (F.PID_P > minPIDProton)
+    return ParticleFilter(make_has_rich_long_protons(), F.FILTER(code))
+
+
+@configurable
+def make_charm_protons_tight(bpvipchi2=9.0,
+                             minPIDProton=None,
+                             minPIDProtonKaon=None,
+                             pt_min=250 * MeV,
+                             p_min=3000.0 * MeV,
+                             max_trchi2dof=3):
+    code = F.require_all(F.PT > pt_min, F.P > p_min, F.CHI2DOF < max_trchi2dof,
+                         F.BPVIPCHI2(make_pvs()) > bpvipchi2)
+    if (minPIDProton and minPIDProtonKaon) is not None:
+        code &= (F.PID_P > minPIDProton)
+        code &= ((F.PID_P - F.PID_K) < minPIDProtonKaon)
+    return ParticleFilter(make_has_rich_long_protons(), F.FILTER(code))
+
+
+@configurable
+def define_child(particle, Criteria):
+
+    if particle == "kaon":
+        dau = make_charm_kaons(
+            bpvipchi2=Criteria["bpvipchi2"],
+            minPIDKaon=Criteria["minPIDK"],
+            min_pt=Criteria["minpT"],
+            min_p=Criteria["minp"],
+            max_trchi2dof=Criteria["maxtrchi2dof"])
+    elif particle == "pion":
+        dau = make_charm_pions(
+            bpvipchi2=Criteria["bpvipchi2"],
+            maxPIDKaon=Criteria["maxPIDK"],
+            min_pt=Criteria["minpT"],
+            min_p=Criteria["minp"],
+            max_trchi2dof=Criteria["maxtrchi2dof"])
+    elif particle == "proton":
+        if (Criteria["minPIDpK"] is None):
+            dau = make_charm_protons_loose(
+                bpvipchi2=Criteria["bpvipchi2"],
+                minPIDProton=Criteria["minPIDp"],
+                min_pt=Criteria["minpT"],
+                min_p=Criteria["minp"],
+                max_trchi2dof=Criteria["maxtrchi2dof"])
+        else:
+            dau = make_charm_protons_tight(
+                bpvipchi2=Criteria["bpvipchi2"],
+                minPIDProton=Criteria["minPIDp"],
+                minPIDProtonKaon=Criteria["minPIDpK"],
+                min_pt=Criteria["minpT"],
+                min_p=Criteria["minp"],
+                max_trchi2dof=Criteria["maxtrchi2dof"])
+    else:
+        print(particle + "is not correctly defined")
+
+    return dau
+
+
+@configurable
+def make_particle_criteria(
+        bpvipchi2=None,
+        minPIDK=None,
+        maxPIDK=None,
+        minPIDp=None,
+        minPIDpK=None,  #PIDp - PIDK > minPIDpK
+        minpT=0,
+        minp=1 * GeV,
+        maxtrchi2dof=5):
+
+    criteria = {
+        "bpvipchi2": bpvipchi2,
+        "minPIDK": minPIDK,
+        "maxPIDK": maxPIDK,
+        "minPIDp": minPIDp,
+        "minPIDpK": minPIDpK,
+        "minpT": minpT,
+        "minp": minp,
+        "maxtrchi2dof": maxtrchi2dof
+    }
+
+    return criteria
+
+
+@configurable
+def make_comb_criteria(mmin=None,
+                       mmax=None,
+                       min_child_pt=0,
+                       max_doca=None,
+                       vchi2pdof_max=25):
+
+    criteria_comb = {
+        "mmin": mmin,
+        "mmax": mmax,
+        "min_child_pt": min_child_pt,
+        "max_doca": max_doca,
+        "vchi2pdof_max": vchi2pdof_max
+    }
+
+    return criteria_comb
+
+
+@configurable
+def make_charm2hadrons(name="",
+                       process="",
+                       decay="",
+                       particle="",
+                       CriteriaParticle="",
+                       CriteriaCombination=""):
+    assert process in ['hlt2', 'spruce'
+                       ], 'Line must be defined as Hlt2 or Sprucing line!'
+
+    dau1 = define_child(particle[0], CriteriaParticle[0])
+    dau2 = define_child(particle[1], CriteriaParticle[1])
+
+    combination_code = F.require_all(
+        in_range(CriteriaCombination["mmin"], F.MASS,
+                 CriteriaCombination["mmin"]),
+        F.SUM(F.PT > CriteriaCombination["min_child_pt"]) > 0)
+
+    if CriteriaCombination["max_doca"] is not None:
+        combination_code &= F.MAXDOCACHI2CUT(CriteriaCombination["max_doca"])
+
+    vertex_code = F.require_all(
+        F.CHI2DOF < CriteriaCombination["vchi2pdof_max"])
+
+    return ParticleCombiner(
+        name=name,
+        Inputs=[dau1, dau2],
+        DecayDescriptor=decay,
+        CombinationCut=combination_code,
+        CompositeCut=vertex_code)
+
+
+@configurable
+def make_charm3hadrons(name="",
+                       process="",
+                       decay="",
+                       particle="",
+                       CriteriaParticle="",
+                       CriteriaCombination=""):
+    """Return a decay maker for either of the decays:
+    1. D-  -> K- pi+ pi+
+    2. Ds- -> K- K+  pi+
+    3. Lc  -> K pi p
+    """
+
+    assert process in ['hlt2', 'spruce'
+                       ], 'Line must be defined as Hlt2 or Sprucing line!'
+    dau1 = define_child(particle[0], CriteriaParticle[0])
+    dau2 = define_child(particle[1], CriteriaParticle[1])
+    dau3 = define_child(particle[2], CriteriaParticle[2])
+
+    combination_code = F.require_all(
+        in_range(CriteriaCombination["mmin"], F.MASS,
+                 CriteriaCombination["mmin"]),
+        F.SUM(F.PT > CriteriaCombination["min_child_pt"]) > 0)
+
+    if CriteriaCombination["max_doca"] is not None:
+        combination_code &= F.MAXDOCACHI2CUT(CriteriaCombination["max_doca"])
+
+    vertex_code = F.require_all(
+        F.CHI2DOF < CriteriaCombination["vchi2pdof_max"])
+    combination12_code = F.require_all(F.MAXDOCACUT(0.5 * mm))
+
+    return ParticleCombiner(
+        name=name,
+        Inputs=[dau1, dau2, dau3],
+        DecayDescriptor=decay,
+        Combination12Cut=combination12_code,
+        CombinationCut=combination_code,
+        CompositeCut=vertex_code)
+
+
+@configurable
+def make_charm4hadrons(name="",
+                       process="",
+                       decay="",
+                       particle="",
+                       CriteriaParticle="",
+                       CriteriaCombination="",
+                       allchild_ptcut=0.25 * GeV,
+                       twochild_ptcut=0.25 * GeV,
+                       onechild_ptcut=0.25 * GeV):
+    assert process in ['hlt2', 'spruce'
+                       ], 'Line must be defined as Hlt2 or Sprucing line!'
+    daus = [
+        define_child(particle[0], CriteriaParticle[0]),
+        define_child(particle[1], CriteriaParticle[1]),
+        define_child(particle[2], CriteriaParticle[2]),
+        define_child(particle[3], CriteriaParticle[3])
+    ]
+
+    allhaschild_pt = F.SUM(F.PT > allchild_ptcut) > 3
+    twohaschild_pt = F.SUM(F.PT > twochild_ptcut) > 1
+    ahaschild_pt = F.SUM(F.PT > onechild_ptcut) > 0
+
+    combination_code = F.require_all(
+        in_range(CriteriaCombination["mmin"], F.MASS,
+                 CriteriaCombination["mmin"]),
+        F.MAXDOCACHI2CUT(CriteriaCombination["max_doca"]), allhaschild_pt,
+        twohaschild_pt, ahaschild_pt)
+
+    vertex_code = F.require_all(
+        F.CHI2DOF < CriteriaCombination["vchi2pdof_max"])
+    combination12_code = F.require_all(F.MAXDOCACUT(0.5 * mm))
+    return ParticleCombiner(
+        name=name,
+        Inputs=daus,
+        DecayDescriptor=decay,
+        CombinationCut=combination_code,
+        Combination12Cut=combination12_code,
+        CompositeCut=vertex_code)

Hlt/Hlt2Conf/python/Hlt2Conf/lines/iftPbPb/builders/dimuon_builders.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 __future__ import absolute_import, division, print_function
+
+from GaudiKernel.SystemOfUnits import GeV, MeV
+
+from RecoConf.reconstruction_objects import make_pvs
+
+from Hlt2Conf.algorithms_thor import (ParticleCombiner, ParticleFilter)
+import Functors as F
+from Functors.math import in_range
+from PyConf import configurable
+from Hlt2Conf.standard_particles import make_ismuon_long_muon
+
+# MASS WINDOWS
+_MASSWINDOW_JPSI = [2900. * MeV, 3300. * MeV]
+_MASSWINDOW_PSI2S = [3500.0 * MeV, 3900.0 * MeV]
+_MASSMINDOW_UPS = [8000.0 * MeV, 12000.0 * MeV]
+_MASSWINDOW_Z = [60. * GeV, 120.0 * GeV]
+_MASSWINDOW_DY1 = [1.2 * GeV, 3.0 * GeV]
+
+# PID CUT FOR SPRUCE AND HLT2 (HLT2 SET TO NONE FOR NOW)
+
+_SPRUCE_PIDmu_JPSI = None  #"PIDmu > 1"
+_SPRUCE_PIDmu_PSI2S = None  #"PIDmu > 1"
+_SPRUCE_PIDmu_UPSILON = None  #"PIDmu > 1"
+_SPRUCE_PIDmu_Z = None  #"PIDmu > 1"
+_SPRUCE_PIDmu_DY = None  #"PIDmu > 1"
+
+_HLT2_PIDmu_JPSI = None
+_HLT2_PIDmu_PSI2S = None
+_HLT2_PIDmu_UPSILON = None
+_HLT2_PIDmu_Z = None
+
+
+def make_muons(pid=None,
+               trchi2_max=3.,
+               trghostprob_max=1,
+               isMuon=False,
+               pt_min=700 * MeV):
+    code = F.require_all(F.GHOSTPROB < trghostprob_max, F.PT > pt_min,
+                         F.CHI2DOF < trchi2_max)
+    if pid is not None:
+        pid_cuts = F.PID_MU > pid
+        code &= (pid_cuts)
+    if isMuon is not None:
+        ismuon_cuts = F.ISMUON
+        code &= ismuon_cuts
+    return ParticleFilter(make_ismuon_long_muon(), F.FILTER(code))
+
+
+## constructor for the different Dimuon lines
+@configurable
+def make_dimuon(particle1,
+                particle2,
+                descriptors,
+                pvs,
+                pt_mom_min=0 * MeV,
+                vchi2pdof_max=25.):
+
+    combination_code = F.require_all(F.PT > pt_mom_min)
+    vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max)
+
+    return ParticleCombiner([particle1, particle2],
+                            DecayDescriptor=descriptors,
+                            CombinationCut=combination_code,
+                            CompositeCut=vertex_code)
+
+
+@configurable
+def make_jpsi(process="hlt2",
+              mmin=_MASSWINDOW_JPSI[0],
+              mmax=_MASSWINDOW_JPSI[1],
+              pt_mom_min=0 * MeV,
+              vchi2pdof_max=25.,
+              pid_cut=-99999,
+              trchi2_max=3.,
+              trghostprob_max=1,
+              isMuon=True,
+              pt_min=700 * MeV):
+    assert process in ['hlt2', 'spruce'
+                       ], 'Line must be defined as Hlt2 or Sprucing line!'
+    pvs = make_pvs()
+    code = F.require_all(in_range(mmin, F.MASS, mmax))
+    if process == 'spruce':
+        muons = make_muons(
+            pid=pid_cut,
+            trchi2_max=trchi2_max,
+            trghostprob_max=trghostprob_max,
+            isMuon=isMuon,
+            pt_min=pt_min)
+        dimuon = make_dimuon(muons, muons, 'J/psi(1S) -> mu- mu+', pvs=pvs)
+    if process == 'hlt2':
+        muons = make_muons(
+            pid=pid_cut,
+            trchi2_max=trchi2_max,
+            trghostprob_max=trghostprob_max,
+            isMuon=isMuon,
+            pt_min=pt_min)
+        dimuon = make_dimuon(muons, muons, 'J/psi(1S) -> mu- mu+', pvs=pvs)
+
+    return ParticleFilter(dimuon, F.FILTER(code))
+
+
+@configurable
+def make_psi2s(process="hlt2",
+               mmin=_MASSWINDOW_PSI2S[0],
+               mmax=_MASSWINDOW_PSI2S[1],
+               pt_mom_min=0 * MeV,
+               vchi2pdof_max=25.,
+               pid_cut=None,
+               trchi2_max=3.,
+               trghostprob_max=1,
+               isMuon=True,
+               pt_min=700 * MeV):
+    assert process in ['hlt2', 'spruce'
+                       ], 'Line must be defined as Hlt2 or Sprucing line!'
+
+    code = F.require_all(in_range(mmin, F.MASS, mmax))
+    pvs = make_pvs()
+    if process == 'spruce':
+        muons = make_muons(
+            pid=pid_cut,
+            trchi2_max=trchi2_max,
+            trghostprob_max=trghostprob_max,
+            isMuon=isMuon,
+            pt_min=pt_min)
+        dimuon = make_dimuon(muons, muons, 'psi(2S) -> mu- mu+', pvs=pvs)
+    if process == 'hlt2':
+        muons = make_muons(
+            pid=pid_cut,
+            trchi2_max=trchi2_max,
+            trghostprob_max=trghostprob_max,
+            isMuon=isMuon,
+            pt_min=pt_min)
+        dimuon = make_dimuon(muons, muons, 'psi(2S) -> mu- mu+', pvs=pvs)
+
+    return ParticleFilter(dimuon, F.FILTER(code))
+
+
+@configurable
+def make_upsilon(process="hlt2",
+                 mmin=_MASSMINDOW_UPS[0],
+                 mmax=_MASSMINDOW_UPS[1],
+                 pt_mom_min=0 * MeV,
+                 vchi2pdof_max=25.,
+                 pid_cut=None,
+                 trchi2_max=3.,
+                 trghostprob_max=1,
+                 isMuon=True,
+                 pt_min=700 * MeV):
+
+    code = F.require_all(in_range(mmin, F.MASS, mmax))
+    pvs = make_pvs()
+    if process == 'spruce':
+        muons = make_muons(
+            pid=pid_cut,
+            trchi2_max=trchi2_max,
+            trghostprob_max=trghostprob_max,
+            isMuon=isMuon,
+            pt_min=pt_min)
+        dimuon = make_dimuon(muons, muons, 'Upsilon(1S) -> mu- mu+', pvs,
+                             pt_mom_min, vchi2pdof_max)
+    if process == 'hlt2':
+        muons = make_muons(
+            pid=pid_cut,
+            trchi2_max=trchi2_max,
+            trghostprob_max=trghostprob_max,
+            isMuon=isMuon,
+            pt_min=pt_min)
+        dimuon = make_dimuon(muons, muons, 'Upsilon(1S) -> mu- mu+', pvs,
+                             pt_mom_min, vchi2pdof_max)
+
+    return ParticleFilter(dimuon, F.FILTER(code))
+
+
+@configurable
+def make_Zboson(process="hlt2",
+                mmin=_MASSWINDOW_Z[0],
+                mmax=_MASSWINDOW_Z[1],
+                pt_mom_min=0 * MeV,
+                vchi2pdof_max=25.,
+                pid_cut=None,
+                trchi2_max=3.,
+                trghostprob_max=1,
+                isMuon=True,
+                pt_min=700 * MeV):
+
+    code = F.require_all(in_range(mmin, F.MASS, mmax))
+    pvs = make_pvs()
+    if process == 'spruce':
+        muons = make_muons(
+            pid=pid_cut,
+            trchi2_max=trchi2_max,
+            trghostprob_max=trghostprob_max,
+            isMuon=isMuon,
+            pt_min=pt_min)
+        dimuon = make_dimuon(muons, muons, 'Z0 -> mu- mu+', pvs, pt_mom_min,
+                             vchi2pdof_max)
+    if process == 'hlt2':
+        muons = make_muons(
+            pid=pid_cut,
+            trchi2_max=trchi2_max,
+            trghostprob_max=trghostprob_max,
+            isMuon=isMuon,
+            pt_min=pt_min)
+        dimuon = make_dimuon(muons, muons, 'Z0 -> mu- mu+', pvs, pt_mom_min,
+                             vchi2pdof_max)
+
+    return ParticleFilter(dimuon, F.FILTER(code))
+
+
+@configurable
+def make_DY(process="hlt2",
+            mmin=_MASSWINDOW_DY1[0],
+            mmax=_MASSWINDOW_DY1[1],
+            pt_mom_min=0 * MeV,
+            vchi2pdof_max=25.,
+            pid_cut=None,
+            trchi2_max=3.,
+            trghostprob_max=1,
+            isMuon=True,
+            pt_min=700 * MeV):
+
+    code = F.require_all(in_range(mmin, F.MASS, mmax))
+    pvs = make_pvs()
+    if process == 'spruce':
+        muons = make_muons(
+            pid=pid_cut,
+            trchi2_max=trchi2_max,
+            trghostprob_max=trghostprob_max,
+            isMuon=isMuon,
+            pt_min=pt_min)
+        dimuon = make_dimuon(muons, muons, 'Z0 -> mu- mu+', pvs, pt_mom_min,
+                             vchi2pdof_max)
+    if process == 'hlt2':
+        muons = make_muons(
+            pid=pid_cut,
+            trchi2_max=trchi2_max,
+            trghostprob_max=trghostprob_max,
+            isMuon=isMuon,
+            pt_min=pt_min)
+        dimuon = make_dimuon(muons, muons, 'Z0 -> mu- mu+', pvs, pt_mom_min,
+                             vchi2pdof_max)
+
+    return ParticleFilter(dimuon, F.FILTER(code))
+
+
+@configurable
+def make_DYLS(process="hlt2",
+              massWind_DY=_MASSWINDOW_DY1,
+              pt_mom_min=0 * MeV,
+              vchi2pdof_max=25.,
+              pid_cut=None,
+              trchi2_max=3.,
+              trghostprob_max=1,
+              isMuon=True,
+              pt_min=700 * MeV):
+
+    mmax = massWind_DY[1]
+    mmin = massWind_DY[0]
+    code = F.require_all(in_range(mmin, F.MASS, mmax))
+    pvs = make_pvs()
+
+    if process == 'spruce':
+        muons = make_muons(
+            pid=pid_cut,
+            isMuon=isMuon,
+            trchi2_max=trchi2_max,
+            trghostprob_max=trghostprob_max,
+            pt_min=pt_min)
+        dimuon = make_dimuon(muons, muons, '[Z0 -> mu+ mu+]cc', pvs,
+                             pt_mom_min, vchi2pdof_max)
+    if process == 'hlt2':
+        muons = make_muons(
+            pid=pid_cut,
+            isMuon=isMuon,
+            trchi2_max=trchi2_max,
+            trghostprob_max=trghostprob_max,
+            pt_min=pt_min)
+        dimuon = make_dimuon(muons, muons, '[Z0 -> mu+ mu+]cc', pvs,
+                             pt_mom_min, vchi2pdof_max)
+    return ParticleFilter(dimuon, F.FILTER(code))

Hlt/Hlt2Conf/python/Hlt2Conf/lines/iftPbPb/hlt2_iftPbPb.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.config import register_line_builder
+from Moore.lines import Hlt2Line
+
+from RecoConf.reconstruction_objects import make_pvs
+from RecoConf.event_filters import require_pvs
+from PyConf import configurable
+
+from RecoConf.event_filters import require_gec
+
+PbSMOG_hlt1 = [
+    "Hlt1BeamGasDecision", "Hlt1PbSMOGMicroBiasDecision",
+    "Hlt1PbSMOGMBOneTrackDecision"
+]
+
+PbPb_hlt1 = [
+    "Hlt1HeavyIonPbPbMicroBias", "Hlt1HeavyIonPbPbPeripheral",
+    "Hlt1HeavyIonPbPbCentral", "Hlt1HeavyIonPbPbUPCMB", "Hlt1PbPbMBOneTrack"
+]
+
+all_lines = {}
+PROCESS = "hlt2"
+
+CUTGEC = 30000
+
+
+## Mass ranges
+## simple prefilters that requires at least one PV ##
+def prefilters():
+    return [require_pvs(make_pvs())]  ## might be adapted in case of PbPb
+
+
+@register_line_builder(all_lines)
+@configurable
+def passthrough(name='Hlt2IFT_PassThrough',
+                prescale=1,
+                persistreco=True,
+                cutGeC=CUTGEC):
+    algs = [require_gec(cut=cutGeC)]
+    return Hlt2Line(
+        name=name, algs=algs, prescale=prescale, persistreco=persistreco)

Hlt/RecoConf/options/hlt2_lead_lead_light_reco_pr_kf_gec_40000_with_mcchecking.py → Hlt/Hlt2Conf/python/Hlt2Conf/lines/iftPbPb/spruce_iftPbPb.py

@@ -8,20 +8,32 @@
 # granted to it by virtue of its status as an Intergovernmental Organization  #
 # or submit itself to any jurisdiction.                                       #
 ###############################################################################
+from Moore.config import SpruceLine, register_line_builder
 
-from Moore import options, run_reconstruction
-from RecoConf.standalone import standalone_hlt2_light_reco_ion
-from RecoConf.hlt2_tracking import make_hlt2_tracks_ion
-
+from RecoConf.reconstruction_objects import upfront_reconstruction
+from PyConf import configurable
 from RecoConf.event_filters import require_gec
-from RecoConf.mc_checking import check_track_resolution
 
-options.histo_file = 'hlt2_lead_lead_track_monitoring_with_mc_histos.root'
-options.ntuple_file = 'hlt2_lead_lead_track_monitoring_with_mc_ntuples.root'
-options.evt_max = 20
+sprucing_lines = {}
+PROCESS = "spruce"
+
+hlt2_ion_triggers = [
+    'Hlt2IFT_PassThroughDecision', 'Hlt2IFT_PbPbMicroBiasDecision'
+]
+
 
-with require_gec.bind(cut=40000), \
-     standalone_hlt2_light_reco_ion.bind(do_mc_checking=True, do_data_monitoring=True, use_pr_kf=True), \
-     check_track_resolution.bind(per_hit_resolutions=True, split_per_type=True), \
-     make_hlt2_tracks_ion.bind( fast_reco=True, use_pr_kf=True, post_fit_selection='(TrGHOSTPROB<0.9999)'):
-    run_reconstruction(options, standalone_hlt2_light_reco_ion)
+@register_line_builder(sprucing_lines)
+@configurable
+def passthrough_PbPb(name='SpruceIFT_PbPbPassThrough',
+                     prescale=1,
+                     skipUT=True,
+                     persistreco=True,
+                     cutGEC=30000,
+                     hlt2_filter=hlt2_ion_triggers):
+    algs = upfront_reconstruction() + [require_gec(cut=cutGEC, skipUT=skipUT)]
+    return SpruceLine(
+        name=name,
+        algs=algs,
+        hlt2_filter_code=hlt2_filter,
+        prescale=prescale,
+        persistreco=persistreco)

Hlt/RecoConf/options/hlt2_lead_lead_light_reco_pr_kf_gec_60000.py → Hlt/Moore/tests/options/lead_lead_input_and_conds_hlt2_2023_velo_open.py

@@ -8,13 +8,16 @@
 # granted to it by virtue of its status as an Intergovernmental Organization  #
 # or submit itself to any jurisdiction.                                       #
 ###############################################################################
+from Moore import options
 
-from Moore import options, run_reconstruction
-from RecoConf.standalone import standalone_hlt2_light_reco_ion
-from RecoConf.hlt2_tracking import make_hlt2_tracks_ion
-from RecoConf.event_filters import require_gec
+# To be removed once data is updated to use a tag including
+# https://gitlab.cern.ch/lhcb-conddb/DDDB/-/merge_requests/107
+from DDDB.CheckDD4Hep import UseDD4Hep
+if UseDD4Hep:
+    from Configurables import DDDBConf
+    DDDBConf().GeometryVersion = 'before-rich1-geom-update-26052022'
 
+options.set_input_and_conds_from_testfiledb(
+    'upgrade-official-PbPb-EPOS-b-8_22-VELO-23.5mm-digi')
 
-with require_gec.bind(cut=60000), \
-     make_hlt2_tracks_ion.bind( fast_reco=True, use_pr_kf=True, post_fit_selection='(TrGHOSTPROB<0.9999)'):
-    run_reconstruction(options, standalone_hlt2_light_reco_ion)
+options.evt_max = 200

Hlt/RecoConf/CMakeLists.txt

@@ -51,8 +51,8 @@ if(BUILD_TESTING AND USE_DD4HEP)
            RecoConf.hlt1_reco_baseline_with_mcchecking
            RecoConf.performance.hlt1_reco_trackingeff_plots
            RecoConf.hlt1_reco_decode_retina_with_mcchecking_FTv6
-           RecoConf.hlt2_lead_lead_light_reco_pr_kf
-           RecoConf.hlt2_lead_lead_light_reco_pr_kf_with_mcchecking
+           RecoConf.hlt2_lead_lead_fast_reco_pr_kf_gec_40000_with_mcchecking_velo_open_2023
+           RecoConf.hlt2_lead_lead_fast_reco_pr_kf_without_UT_gec_25000_with_mcchecking_velo_open_2023
            RecoConf.hlt2_light_reco_pr_kf_velo_open_without_UT_with_mcchecking
            RecoConf.hlt2_pr_kf_longmuon
            RecoConf.hlt2_protoparticles_baseline

Hlt/RecoConf/options/hlt2_lead_lead_light_reco_pr_kf_gec_40000.py → Hlt/RecoConf/options/examples/lead_lead_input_and_conds_hlt2_2023_velo_open_mdf.py

@@ -8,13 +8,16 @@
 # granted to it by virtue of its status as an Intergovernmental Organization  #
 # or submit itself to any jurisdiction.                                       #
 ###############################################################################
+from Moore import options
 
-from Moore import options, run_reconstruction
-from RecoConf.standalone import standalone_hlt2_light_reco_ion
-from RecoConf.hlt2_tracking import make_hlt2_tracks_ion
-from RecoConf.event_filters import require_gec
+# To be removed once data is updated to use a tag including
+# https://gitlab.cern.ch/lhcb-conddb/DDDB/-/merge_requests/107
+from DDDB.CheckDD4Hep import UseDD4Hep
+if UseDD4Hep:
+    from Configurables import DDDBConf
+    DDDBConf().GeometryVersion = 'before-rich1-geom-update-26052022'
 
+options.set_input_and_conds_from_testfiledb(
+    'upgrade-202305-PbPb-EPOS-b-8_22-VELO-10mm-mdf')
 
-with require_gec.bind(cut=40000), \
-     make_hlt2_tracks_ion.bind( fast_reco=True, use_pr_kf=True, post_fit_selection='(TrGHOSTPROB<0.9999)'):
-    run_reconstruction(options, standalone_hlt2_light_reco_ion)
+options.evt_max = 300

Hlt/RecoConf/options/examples/mdf_input_and_conds_data_ion_2022.py

+###############################################################################
+# (c) Copyright 2022 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.                                       #
+###############################################################################
+#####
+## part I
+#####
+from Moore import options
+from DDDB.CheckDD4Hep import UseDD4Hep
+from PyConf.application import metainfo_repos
+
+options.input_files = [
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-203228-575.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-203342-100.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-203446-741.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-203552-259.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-203715-831.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-203824-972.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-203936-532.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-204049-286.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-204159-306.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-204307-403.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-204418-140.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-204527-058.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-204632-862.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-204742-885.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-204854-007.mdf',
+    'mdf:root://eoslhcb.cern.ch//eos/lhcb/hlt2/LHCb/0000254507/Run_0000254507_HLT24640_20230117-205826-007.mdf',
+]
+
+# Set the tags depending on DetDesc or DD4HEP.
+# Note DetDesc requires simulation to be false.
+options.dddb_tag = "upgrade/master"
+options.conddb_tag = 'md_VP_SciFi_macromicrosurvey_from20220923'
+options.geometry_version = "trunk"
+options.conditions_version = 'master'
+options.simulation = not UseDD4Hep
+metainfo_repos.global_bind(extra_central_tags=['commissioning'])
+
+options.evt_max = 5000
+options.input_type = "MDF"
+
+if UseDD4Hep:
+    # The data have no UT so exclude it here.
+    from Configurables import LHCb__Det__LbDD4hep__DD4hepSvc as DD4hepSvc
+    dd4hepSvc = DD4hepSvc()
+    dd4hepSvc.DetectorList = [
+        '/world', 'VP', 'FT', 'Magnet', 'Rich1', 'Rich2', 'Ecal', 'Hcal',
+        'Muon'
+    ]

Hlt/RecoConf/options/hlt2_lead_lead_fast_reco_pr_kf_gec_40000.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.hlt1_tracking import (make_VeloClusterTrackingSIMD,
+                                    make_reco_pvs, make_PatPV3DFuture_pvs)
+from RecoConf.hlt2_tracking import make_TrackBestTrackCreator_tracks
+
+from RecoConf.event_filters import require_gec
+from RecoConf.hlt1_muonid import make_muon_hits
+
+from RecoConf.standalone import standalone_hlt2_light_reco
+
+from PyConf.Algorithms import PrHybridSeeding
+
+with standalone_hlt2_light_reco.bind(do_mc_checking=False, do_data_monitoring=False, use_pr_kf=False, fast_reco=False),\
+    require_gec.bind(cut=40000,skipUT=False), \
+    PrHybridSeeding.bind(RemoveBeamHole=True, RemoveClones_forLead=True), \
+    make_VeloClusterTrackingSIMD.bind(SkipForward=4),\
+    make_reco_pvs.bind(make_pvs_from_velo_tracks=make_PatPV3DFuture_pvs),\
+    make_TrackBestTrackCreator_tracks.bind(max_ghost_prob=0.8),\
+    make_muon_hits.bind(geometry_version=3):
+    config = run_reconstruction(options, standalone_hlt2_light_reco)

Hlt/RecoConf/options/hlt2_lead_lead_light_reco_pr_kf_with_mcchecking.py → Hlt/RecoConf/options/hlt2_lead_lead_fast_reco_pr_kf_gec_40000_with_mcchecking.py

@@ -8,24 +8,29 @@
 # granted to it by virtue of its status as an Intergovernmental Organization  #
 # or submit itself to any jurisdiction.                                       #
 ###############################################################################
-"""
-This set of options is used for reconstruction development purposes,
-and assumes that the input contains MCHits (i.e. is of `Exended`
-DST/digi type).
-"""
 
 from Moore import options, run_reconstruction
-from RecoConf.standalone import standalone_hlt2_light_reco
+from RecoConf.hlt1_tracking import (make_VeloClusterTrackingSIMD,
+                                    make_reco_pvs, make_PatPV3DFuture_pvs)
+from RecoConf.hlt2_tracking import make_TrackBestTrackCreator_tracks
 
 from RecoConf.event_filters import require_gec
+from RecoConf.hlt1_muonid import make_muon_hits
 from RecoConf.mc_checking import check_track_resolution
 
+from RecoConf.standalone import standalone_hlt2_light_reco
+
+from PyConf.Algorithms import PrHybridSeeding
+
 options.histo_file = 'hlt2_lead_lead_track_monitoring_with_mc_histos.root'
 options.ntuple_file = 'hlt2_lead_lead_track_monitoring_with_mc_ntuples.root'
-options.evt_max = 20
 
-# A higher gec cut than 40000 leads to unreasonable processing times for the nightlies.
-with standalone_hlt2_light_reco.bind(do_mc_checking=True, do_data_monitoring=True, use_pr_kf=True), \
-     check_track_resolution.bind(per_hit_resolutions=False, split_per_type=True), \
-     require_gec.bind(cut=40000):
-    run_reconstruction(options, standalone_hlt2_light_reco)
+with standalone_hlt2_light_reco.bind(do_mc_checking=True, do_data_monitoring=False, use_pr_kf=True, fast_reco=True),\
+    require_gec.bind(cut=40000,skipUT=False), \
+    make_VeloClusterTrackingSIMD.bind(SkipForward=4),\
+    PrHybridSeeding.bind(RemoveBeamHole=True, RemoveClones_forLead=True), \
+    check_track_resolution.bind(per_hit_resolutions=True, split_per_type=True), \
+    make_reco_pvs.bind(make_pvs_from_velo_tracks=make_PatPV3DFuture_pvs),\
+    make_TrackBestTrackCreator_tracks.bind(max_ghost_prob=0.8),\
+    make_muon_hits.bind(geometry_version=3):
+    config = run_reconstruction(options, standalone_hlt2_light_reco)

Hlt/RecoConf/options/hlt2_lead_lead_light_reco_pr_kf_without_UT_gec_40000_with_mcchecking.py → Hlt/RecoConf/options/hlt2_lead_lead_fast_reco_pr_kf_without_UT_gec_25000.py

@@ -10,18 +10,22 @@
 ###############################################################################
 
 from Moore import options, run_reconstruction
-from RecoConf.standalone import standalone_hlt2_light_reco_ion_without_UT
-from RecoConf.hlt2_tracking import make_hlt2_tracks_ion_without_UT
+from RecoConf.hlt1_tracking import (make_VeloClusterTrackingSIMD,
+                                    make_reco_pvs, make_PatPV3DFuture_pvs)
+from RecoConf.hlt2_tracking import make_TrackBestTrackCreator_tracks
 
 from RecoConf.event_filters import require_gec
-from RecoConf.mc_checking import check_track_resolution
+from RecoConf.hlt1_muonid import make_muon_hits
 
-options.histo_file = 'hlt2_lead_lead_track_monitoring_with_mc_histos.root'
-options.ntuple_file = 'hlt2_lead_lead_track_monitoring_with_mc_ntuples.root'
-options.evt_max = 20
+from RecoConf.standalone import standalone_hlt2_light_reco_without_UT
 
-with require_gec.bind(cut=40000), \
-     standalone_hlt2_light_reco_ion_without_UT.bind(do_mc_checking=True, do_data_monitoring=True, fast_reco=True), \
-     check_track_resolution.bind(per_hit_resolutions=True, split_per_type=True), \
-     make_hlt2_tracks_ion_without_UT.bind( fast_reco=True, use_pr_kf=True, post_fit_selection='(TrGHOSTPROB<0.8)'):
-    run_reconstruction(options, standalone_hlt2_light_reco_ion_without_UT)
+from PyConf.Algorithms import PrHybridSeeding
+
+with standalone_hlt2_light_reco_without_UT.bind(do_mc_checking=False, do_data_monitoring=False, use_pr_kf=True, fast_reco=True),\
+    require_gec.bind(cut=25000,skipUT=True), \
+    PrHybridSeeding.bind(RemoveBeamHole=True, RemoveClones_forLead=True), \
+    make_VeloClusterTrackingSIMD.bind(SkipForward=4),\
+    make_reco_pvs.bind(make_pvs_from_velo_tracks=make_PatPV3DFuture_pvs),\
+    make_TrackBestTrackCreator_tracks.bind(max_ghost_prob=0.8),\
+    make_muon_hits.bind(geometry_version=3):
+    config = run_reconstruction(options, standalone_hlt2_light_reco_without_UT)

Hlt/RecoConf/options/hlt2_lead_lead_light_reco_pr_kf_gec_60000_with_mcchecking.py → Hlt/RecoConf/options/hlt2_lead_lead_fast_reco_pr_kf_without_UT_gec_25000_with_mcchecking.py

@@ -10,18 +10,28 @@
 ###############################################################################
 
 from Moore import options, run_reconstruction
-from RecoConf.standalone import standalone_hlt2_light_reco_ion
-from RecoConf.hlt2_tracking import make_hlt2_tracks_ion
+from RecoConf.hlt1_tracking import (make_VeloClusterTrackingSIMD,
+                                    make_reco_pvs, make_PatPV3DFuture_pvs)
+from RecoConf.hlt2_tracking import make_TrackBestTrackCreator_tracks
 
 from RecoConf.event_filters import require_gec
 from RecoConf.mc_checking import check_track_resolution
+from RecoConf.hlt1_muonid import make_muon_hits
+
+from PyConf.Algorithms import PrHybridSeeding
 
 options.histo_file = 'hlt2_lead_lead_track_monitoring_with_mc_histos.root'
 options.ntuple_file = 'hlt2_lead_lead_track_monitoring_with_mc_ntuples.root'
-options.evt_max = 20
 
-with require_gec.bind(cut=60000), \
-     standalone_hlt2_light_reco_ion.bind(do_mc_checking=True, do_data_monitoring=True, use_pr_kf=True), \
-     check_track_resolution.bind(per_hit_resolutions=True, split_per_type=True), \
-     make_hlt2_tracks_ion.bind( fast_reco=True, use_pr_kf=True, post_fit_selection='(TrGHOSTPROB<0.9999)'):
-    run_reconstruction(options, standalone_hlt2_light_reco_ion)
+from RecoConf.standalone import standalone_hlt2_light_reco_without_UT
+
+
+with standalone_hlt2_light_reco_without_UT.bind(do_mc_checking=True, do_data_monitoring=False, use_pr_kf=True, fast_reco=True),\
+    require_gec.bind(cut=25000,skipUT=True), \
+    PrHybridSeeding.bind(RemoveBeamHole=True, RemoveClones_forLead=True), \
+    make_VeloClusterTrackingSIMD.bind(SkipForward=4),\
+    make_reco_pvs.bind(make_pvs_from_velo_tracks=make_PatPV3DFuture_pvs),\
+    check_track_resolution.bind(per_hit_resolutions=True, split_per_type=True), \
+    make_TrackBestTrackCreator_tracks.bind(max_ghost_prob=0.8),\
+    make_muon_hits.bind(geometry_version=3):
+    config = run_reconstruction(options, standalone_hlt2_light_reco_without_UT)