Update HLT1 PbPb sequence for 2024 ion run

The main update consists of adding and tuning the SMOG2 line to reduce the Pb-SMOG2 rate.

configuration/python/AllenConf/HLT1_PbPb.py

@@ -17,7 +17,6 @@ from AllenConf.hlt1_calibration_lines import (
 from AllenConf.hlt1_monitoring_lines import (
     make_velo_micro_bias_line,
     make_beam_gas_line,
-    make_velo_clusters_micro_bias_line,
     make_odin_event_type_with_decoding_line,
 )
 from AllenConf.hlt1_heavy_ions_lines import (
@@ -29,6 +28,7 @@ from AllenConf.hlt1_heavy_ions_lines import (
 from AllenConf.hlt1_inclusive_hadron_lines import make_kstopipi_line, make_lambda2ppi_line
 from AllenConf.hlt1_charm_lines import make_d2kk_line, make_d2pipi_line, make_d2kpi_line
 from AllenConf.hlt1_muon_lines import make_one_muon_track_line, make_di_muon_mass_line
+from AllenConf.hlt1_smog2_lines import *
 from AllenConf.velo_reconstruction import decode_velo
 from AllenConf.calo_reconstruction import decode_calo
 from AllenConf.validators import rate_validation
@@ -40,6 +40,146 @@ from AllenConf.enum_types import TrackingType, includes_matching
 from .HLT1 import default_bgi_activity_lines
 
 
+def default_SMOG2_lines(reconstructed_objects,
+                        with_muon=True,
+                        with_v0s=True,
+                        min_z=-541.,
+                        max_z=-341.,
+                        enable_tupling=False):
+
+    velo_tracks = reconstructed_objects["velo_tracks"]
+    long_tracks = reconstructed_objects["long_tracks"]
+    long_track_particles = reconstructed_objects["long_track_particles"]
+    dihadrons = reconstructed_objects["dihadron_secondary_vertices"]
+    prompt_dihadrons = reconstructed_objects[
+        "prompt_dihadron_secondary_vertices"]
+    v0s = reconstructed_objects["v0_secondary_vertices"]
+    dileptons = reconstructed_objects["dilepton_secondary_vertices"]
+
+    lines = [
+        make_SMOG2_ditrack_line(
+            dihadrons,
+            m1=139.57,
+            m2=493.68,
+            mMother=1864.83,
+            mWindow=100.,
+            min_z=min_z,
+            max_z=max_z,
+            minEitherTrackPt=800.,
+            minTrackPt=500.,
+            minTrackIPCHI2=7.,
+            minFDCHI2=25.,
+            name="Hlt1PbSMOGD2Kpi",
+            pre_scaler=1.,
+            enable_tupling=enable_tupling),
+        make_SMOG2_ditrack_line(
+            prompt_dihadrons,
+            m1=938.27,
+            m2=938.27,
+            mMother=3000.,
+            mWindow=200.,
+            min_z=min_z,
+            max_z=max_z,
+            minTrackIPCHI2=0.,
+            maxTrackIPCHI2=5.,
+            maxFDCHI2=20.,
+            minTrackP=25000.,
+            minTrackPt=1000.,
+            minEitherTrackPt=1200.,
+            name="Hlt1PbSMOGetacTopp",
+            pre_scaler=1.,
+            enable_tupling=enable_tupling),
+        make_SMOG2_kstopipi_line(
+            dihadrons,
+            min_z=min_z,
+            max_z=max_z,
+            name="Hlt1PbSMOGKsTopipi",
+            minTrackPt=250.,
+            minMass=450.,
+            pre_scaler=0.3,
+            enable_tupling=enable_tupling),
+        make_SMOG2_ditrack_line(
+            dihadrons,
+            minTrackPt=500.,
+            minEitherTrackPt=800.,
+            min_z=min_z,
+            max_z=max_z,
+            minMdipion=1300,
+            minFDCHI2=25.,
+            minTrackIPCHI2=7.,
+            name="Hlt1PbSMOG2BodyGeneric",
+            enable_monitoring=False,
+            enable_tupling=False,
+            pre_scaler=1.),
+        make_SMOG2_ditrack_line(
+            prompt_dihadrons,
+            minTrackPt=400.,
+            minEitherTrackPt=400.,
+            min_z=min_z,
+            max_z=max_z,
+            minTrackIPCHI2=0.,
+            enable_monitoring=False,
+            enable_tupling=False,
+            name="Hlt1PbSMOG2BodyGenericPrompt",
+            pre_scaler=0.01),
+        make_SMOG2_singletrack_line(
+            long_tracks,
+            long_track_particles,
+            name="Hlt1PbSMOGSingleTrackVeryHighPt",
+            minPt=5000.,
+            pre_scaler=1,
+            min_z=min_z,
+            max_z=max_z),
+        make_SMOG2_singletrack_line(
+            long_tracks,
+            long_track_particles,
+            name="Hlt1PbSMOGSingleTrackHighPt",
+            minPt=3000.,
+            pre_scaler=1,
+            min_z=min_z,
+            max_z=max_z)
+    ]
+
+    if with_muon:
+        muonid = reconstructed_objects["muonID"]
+        lines += [
+            make_SMOG2_dimuon_highmass_line(
+                dileptons,
+                long_tracks,
+                muonid,
+                maxChi2Corr=9999.,
+                enable_tupling=enable_tupling,
+                name="Hlt1PbSMOGDiMuonHighMass"),
+            make_SMOG2_single_muon_line(
+                long_tracks,
+                long_track_particles,
+                muonid,
+                maxChi2Corr=1.8,
+                MinPt=700,
+                name="Hlt1PbSMOGSingleMuon",
+                pre_scaler=1)
+        ]
+
+    if with_v0s:
+        lines += [
+            make_lambda2ppi_line(
+                v0s,
+                name="Hlt1PbSMOGL0Toppi",
+                minPVZ=min_z,
+                maxPVZ=max_z,
+                minVZ=min_z,
+                maxVtxChi2=10.,
+                minDIRA=0.99985,
+                minpipchi2=16.,
+                minpiipchi2=42.,
+                minpipt=150.,
+                enable_monitoring=True,
+                enable_tupling=enable_tupling)
+        ]
+
+    return [line_maker(line) for line in lines]
+
+
 def default_physics_lines(reconstructed_objects, prescale, reco_particles,
                           with_muon):
 
@@ -47,7 +187,6 @@ def default_physics_lines(reconstructed_objects, prescale, reco_particles,
     long_tracks = reconstructed_objects["long_tracks"]
     long_track_particles = reconstructed_objects["long_track_particles"]
     decoded_calo = reconstructed_objects["decoded_calo"]
-    ecal_clusters = reconstructed_objects["ecal_clusters"]
     pvs = reconstructed_objects["pvs"]
     dihadrons = reconstructed_objects["dihadron_secondary_vertices"]
     dileptons = reconstructed_objects["dilepton_secondary_vertices"]
@@ -321,6 +460,8 @@ def setup_hlt1_node(withMCChecking=False,
                     with_AC_split=False,
                     prescale=False,
                     with_muon=True,
+                    with_v0s=True,
+                    enableTupling=False,
                     reco_particles=False,
                     bx_type=None,
                     tae_passthrough=True,
@@ -416,6 +557,11 @@ def setup_hlt1_node(withMCChecking=False,
         with line_maker.bind(prefilter=prefilters):
             physics_lines = default_physics_lines(
                 reconstructed_objects, prescale, reco_particles, with_muon)
+            physics_lines += default_SMOG2_lines(
+                reconstructed_objects,
+                with_muon,
+                with_v0s,
+                enable_tupling=enableTupling)
         with line_maker.bind(prefilter=prefilter_photon_velo_upc):
             physics_lines += upc_physics_lines(reconstructed_objects)