Merge branch 'axu-tune-ccbaryon-hadronic' into '2024-patches'

[Charm] Tune doubly charmed baryon hadronic lines

See merge request !3300

Hlt/Hlt2Conf/python/Hlt2Conf/lines/charm/ccbaryon_to_cbaryon_sl.py

@@ -42,6 +42,9 @@ from Hlt2Conf.standard_particles import (
     make_has_rich_long_protons, make_ismuon_long_muon)
 from Hlt2Conf.algorithms_thor import (ParticleFilter, ParticleCombiner)
 from .prefilters import charm_prefilters
+from .ccbaryon_hadronic import create_hc_mon
+
+from .particle_properties import _LC_M, _XICP_M, _XICZ_M, _OMEGAC_M
 
 all_lines = {}
 
@@ -230,6 +233,7 @@ def xiccpp_to_lcpmu_line(name="Hlt2Charm_XiccppToLcpMupNu", prescale=1):
     pions = filter_long_pions()
     protons = filter_long_protons()
     lambdac = make_lc_to_pkpi(protons=protons, kaons=kaons, pions=pions)
+    lambdac_mon = create_hc_mon(lambdac, name, _LC_M)
     line_alg = make_cc_baryon(
         baryon=lambdac,
         particle=muons,
@@ -243,8 +247,9 @@ def xiccpp_to_lcpmu_line(name="Hlt2Charm_XiccppToLcpMupNu", prescale=1):
         bpvdira_min=0.99)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [muons, protons, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [muons, protons, line_alg, lambdac_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -254,6 +259,7 @@ def xiccpp_to_xicpmu_line(name="Hlt2Charm_XiccppToXicpMupNu", prescale=1):
     pions = filter_long_pions()
     protons = filter_long_protons()
     xicp = make_xicp_to_pkpi(protons=protons, kaons=kaons, pions=pions)
+    xicp_mon = create_hc_mon(xicp, name, _XICP_M)
     line_alg = make_cc_baryon(
         baryon=xicp,
         particle=muons,
@@ -267,8 +273,9 @@ def xiccpp_to_xicpmu_line(name="Hlt2Charm_XiccppToXicpMupNu", prescale=1):
         bpvdira_min=0.99)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [muons, protons, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [muons, protons, line_alg, xicp_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -278,6 +285,7 @@ def xiccp_to_xic0mu_line(name="Hlt2Charm_XiccpToXic0MupNu", prescale=1):
     pions = filter_long_pions()
     protons = filter_long_protons()
     xic0 = make_xic0_to_pkkpi(protons=protons, kaons=kaons, pions=pions)
+    xic0_mon = create_hc_mon(xic0, name, _XICZ_M)
     line_alg = make_cc_baryon(
         baryon=xic0,
         particle=muons,
@@ -291,8 +299,9 @@ def xiccp_to_xic0mu_line(name="Hlt2Charm_XiccpToXic0MupNu", prescale=1):
         bpvdira_min=0.99)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [muons, protons, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [muons, protons, line_alg, xic0_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -302,6 +311,7 @@ def xiccpp_to_lcpmu_wsline(name="Hlt2Charm_XiccppToLcpMumNu_WS", prescale=0.5):
     pions = filter_long_pions()
     protons = filter_long_protons()
     lambdac = make_lc_to_pkpi(protons=protons, kaons=kaons, pions=pions)
+    lambdac_mon = create_hc_mon(lambdac, name, _LC_M)
     line_alg = make_cc_baryon(
         baryon=lambdac,
         particle=muons,
@@ -315,8 +325,9 @@ def xiccpp_to_lcpmu_wsline(name="Hlt2Charm_XiccppToLcpMumNu_WS", prescale=0.5):
         bpvdira_min=0.99)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [muons, protons, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [muons, protons, line_alg, lambdac_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -326,6 +337,7 @@ def xiccpp_to_xicpmu_wsline(name="Hlt2Charm_XiccppToXicpMumNu_WS", prescale=1):
     pions = filter_long_pions()
     protons = filter_long_protons()
     xicp = make_xicp_to_pkpi(protons=protons, kaons=kaons, pions=pions)
+    xicp_mon = create_hc_mon(xicp, name, _XICP_M)
     line_alg = make_cc_baryon(
         baryon=xicp,
         particle=muons,
@@ -339,8 +351,9 @@ def xiccpp_to_xicpmu_wsline(name="Hlt2Charm_XiccppToXicpMumNu_WS", prescale=1):
         bpvdira_min=0.99)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [muons, protons, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [muons, protons, line_alg, xicp_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -350,6 +363,7 @@ def xiccp_to_xic0mu_wsline(name="Hlt2Charm_XiccpToXic0MumNu_WS", prescale=1):
     pions = filter_long_pions()
     protons = filter_long_protons()
     xic0 = make_xic0_to_pkkpi(protons=protons, kaons=kaons, pions=pions)
+    xic0_mon = create_hc_mon(xic0, name, _XICZ_M)
     line_alg = make_cc_baryon(
         baryon=xic0,
         particle=muons,
@@ -363,8 +377,9 @@ def xiccp_to_xic0mu_wsline(name="Hlt2Charm_XiccpToXic0MumNu_WS", prescale=1):
         bpvdira_min=0.99)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [muons, protons, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [muons, protons, line_alg, xic0_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -374,6 +389,7 @@ def omegaccp_to_omegac0mu_line(name="Hlt2Charm_OccpToOc0MupNu", prescale=1):
     pions = filter_long_pions()
     protons = filter_long_protons()
     omegac0 = make_omegac0_to_pkkpi(protons=protons, kaons=kaons, pions=pions)
+    omegac0_mon = create_hc_mon(omegac0, name, _OMEGAC_M)
     line_alg = make_cc_baryon(
         baryon=omegac0,
         particle=muons,
@@ -387,8 +403,9 @@ def omegaccp_to_omegac0mu_line(name="Hlt2Charm_OccpToOc0MupNu", prescale=1):
         bpvdira_min=0.99)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [muons, protons, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [muons, protons, line_alg, omegac0_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -398,6 +415,7 @@ def omegaccp_to_xic0mu_line(name="Hlt2Charm_OccpToXic0MupNu", prescale=1):
     pions = filter_long_pions()
     protons = filter_long_protons()
     xic0 = make_xic0_to_pkkpi(protons=protons, kaons=kaons, pions=pions)
+    xic0_mon = create_hc_mon(xic0, name, _XICZ_M)
     line_alg = make_cc_baryon(
         baryon=xic0,
         particle=muons,
@@ -411,8 +429,9 @@ def omegaccp_to_xic0mu_line(name="Hlt2Charm_OccpToXic0MupNu", prescale=1):
         bpvdira_min=0.99)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [muons, protons, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [muons, protons, line_alg, xic0_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -423,6 +442,7 @@ def omegaccp_to_omegac0mu_wsline(name="Hlt2Charm_OccpToOc0MumNu_WS",
     pions = filter_long_pions()
     protons = filter_long_protons()
     omegac0 = make_omegac0_to_pkkpi(protons=protons, kaons=kaons, pions=pions)
+    omegac0_mon = create_hc_mon(omegac0, name, _OMEGAC_M)
     line_alg = make_cc_baryon(
         baryon=omegac0,
         particle=muons,
@@ -436,8 +456,9 @@ def omegaccp_to_omegac0mu_wsline(name="Hlt2Charm_OccpToOc0MumNu_WS",
         bpvdira_min=0.99)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [muons, protons, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [muons, protons, line_alg, omegac0_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -447,6 +468,7 @@ def omegaccp_to_xic0mu_wsline(name="Hlt2Charm_OccpToXic0MumNu_WS", prescale=1):
     pions = filter_long_pions()
     protons = filter_long_protons()
     xic0 = make_xic0_to_pkkpi(protons=protons, kaons=kaons, pions=pions)
+    xic0_mon = create_hc_mon(xic0, name, _XICZ_M)
     line_alg = make_cc_baryon(
         baryon=xic0,
         particle=muons,
@@ -460,5 +482,6 @@ def omegaccp_to_xic0mu_wsline(name="Hlt2Charm_OccpToXic0MumNu_WS", prescale=1):
         bpvdira_min=0.99)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [muons, protons, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [muons, protons, line_alg, xic0_mon],
+        prescale=prescale,
+        monitoring_variables=())

Hlt/Hlt2Conf/python/Hlt2Conf/lines/charm/ccbaryon_to_hyperon_dh.py

@@ -35,6 +35,9 @@ from ...standard_particles import (
     make_has_rich_long_kaons, make_has_rich_long_pions)
 from ...algorithms_thor import ParticleFilter, ParticleCombiner
 from .prefilters import charm_prefilters
+from .ccbaryon_hadronic import create_hc_mon
+
+from .particle_properties import _XIM_M, _DP_M, _DS_M
 
 all_lines = {}
 
@@ -255,9 +258,11 @@ def omegaccp_to_ximdppip_lll_line(name="Hlt2Charm_OccpToDpXimPip_LLL",
         vchi2dof_max=10.,
         m_min=1790 * MeV,
         m_max=1940 * MeV)
+    dp_mon = create_hc_mon(dp, name, _DP_M)
     xim_lambdas = make_lambdall()
     xim_pions = filter_long_pions(mipchi2_min=16.)
     xim = make_xim_to_lambdapi(lambdas=xim_lambdas, pions=xim_pions)
+    xim_mon = create_hc_mon(xim, name, _XIM_M, name="hyp")
     pions = filter_long_pions(mipchi2_min=1.)
     line_alg = make_cc_baryon(
         particle1=xim,
@@ -273,8 +278,9 @@ def omegaccp_to_ximdppip_lll_line(name="Hlt2Charm_OccpToDpXimPip_LLL",
         bpvdira_min=0.)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [xim, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [xim, line_alg, dp_mon, xim_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -293,9 +299,11 @@ def omegaccp_to_ximdppip_ddl_line(name="Hlt2Charm_OccpToDpXimPip_DDL",
         vchi2dof_max=10.,
         m_min=1790 * MeV,
         m_max=1940 * MeV)
+    dp_mon = create_hc_mon(dp, name, _DP_M)
     xim_lambdas = make_lambdadd()
     xim_pions = filter_long_pions(mipchi2_min=16.)
     xim = make_xim_to_lambdapi(lambdas=xim_lambdas, pions=xim_pions)
+    xim_mon = create_hc_mon(xim, name, _XIM_M, name="hyp")
     pions = filter_long_pions(mipchi2_min=1.)
     line_alg = make_cc_baryon(
         particle1=xim,
@@ -311,8 +319,9 @@ def omegaccp_to_ximdppip_ddl_line(name="Hlt2Charm_OccpToDpXimPip_DDL",
         bpvdira_min=0.)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [xim, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [xim, line_alg, dp_mon, xim_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -331,9 +340,11 @@ def xiccp_to_ximdsppip_lll_line(name="Hlt2Charm_XiccpToDspXimPip_LLL",
         vchi2dof_max=10.,
         m_min=1890 * MeV,
         m_max=2040 * MeV)
+    dsp_mon = create_hc_mon(dsp, name, _DS_M)
     xim_lambdas = make_lambdall()
     xim_pions = filter_long_pions(mipchi2_min=16.)
     xim = make_xim_to_lambdapi(lambdas=xim_lambdas, pions=xim_pions)
+    xim_mon = create_hc_mon(xim, name, _XIM_M, name="hyp")
     pions = filter_long_pions(mipchi2_min=1.)
     line_alg = make_cc_baryon(
         particle1=xim,
@@ -349,8 +360,9 @@ def xiccp_to_ximdsppip_lll_line(name="Hlt2Charm_XiccpToDspXimPip_LLL",
         bpvdira_min=0.)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [xim, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [xim, line_alg, dsp_mon, xim_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -369,9 +381,11 @@ def xiccp_to_ximdsppip_ddl_line(name="Hlt2Charm_XiccpToDspXimPip_DDL",
         vchi2dof_max=10.,
         m_min=1890 * MeV,
         m_max=2040 * MeV)
+    dsp_mon = create_hc_mon(dsp, name, _DS_M)
     xim_lambdas = make_lambdadd()
     xim_pions = filter_long_pions(mipchi2_min=16.)
     xim = make_xim_to_lambdapi(lambdas=xim_lambdas, pions=xim_pions)
+    xim_mon = create_hc_mon(xim, name, _XIM_M, name="hyp")
     pions = filter_long_pions(mipchi2_min=1.)
     line_alg = make_cc_baryon(
         particle1=xim,
@@ -387,8 +401,9 @@ def xiccp_to_ximdsppip_ddl_line(name="Hlt2Charm_XiccpToDspXimPip_DDL",
         bpvdira_min=0.)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [xim, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [xim, line_alg, dsp_mon, xim_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -406,9 +421,11 @@ def xiccp_to_ximdpkp_lll_line(name="Hlt2Charm_XiccpToDpXimKp_LLL", prescale=1):
         vchi2dof_max=10.,
         m_min=1790 * MeV,
         m_max=1940 * MeV)
+    dp_mon = create_hc_mon(dp, name, _DP_M)
     xim_lambdas = make_lambdall()
     xim_pions = filter_long_pions(mipchi2_min=16.)
     xim = make_xim_to_lambdapi(lambdas=xim_lambdas, pions=xim_pions)
+    xim_mon = create_hc_mon(xim, name, _XIM_M, name="hyp")
     kaons = filter_long_kaons(mipchi2_min=1.)
     line_alg = make_cc_baryon(
         particle1=xim,
@@ -424,8 +441,9 @@ def xiccp_to_ximdpkp_lll_line(name="Hlt2Charm_XiccpToDpXimKp_LLL", prescale=1):
         bpvdira_min=0.)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [xim, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [xim, line_alg, dp_mon, xim_mon],
+        prescale=prescale,
+        monitoring_variables=())
 
 
 @register_line_builder(all_lines)
@@ -443,9 +461,11 @@ def xiccp_to_ximdpkp_ddl_line(name="Hlt2Charm_XiccpToDpXimKp_DDL", prescale=1):
         vchi2dof_max=10.,
         m_min=1790 * MeV,
         m_max=1940 * MeV)
+    dp_mon = create_hc_mon(dp, name, _DP_M)
     xim_lambdas = make_lambdadd()
     xim_pions = filter_long_pions(mipchi2_min=16.)
     xim = make_xim_to_lambdapi(lambdas=xim_lambdas, pions=xim_pions)
+    xim_mon = create_hc_mon(xim, name, _XIM_M, name="hyp")
     kaons = filter_long_kaons(mipchi2_min=1.)
     line_alg = make_cc_baryon(
         particle1=xim,
@@ -461,5 +481,6 @@ def xiccp_to_ximdpkp_ddl_line(name="Hlt2Charm_XiccpToDpXimKp_DDL", prescale=1):
         bpvdira_min=0.)
     return Hlt2Line(
         name=name,
-        algs=charm_prefilters() + [xim, line_alg],
-        prescale=prescale)
+        algs=charm_prefilters() + [xim, line_alg, dp_mon, xim_mon],
+        prescale=prescale,
+        monitoring_variables=())

Hlt/Hlt2Conf/python/Hlt2Conf/lines/charm/particle_properties.py

@@ -21,6 +21,7 @@ _PION_M = 139.57061 * MeV  # +/- 0.00024
 _KAON_M = 493.677 * MeV  # +/- 0.016
 _KS_M = 497.611 * MeV  # +/- 0.013
 _ETA_M = 547.862 * MeV  # +/- 0.017
+_XIM_M = 1321.71 * MeV  # +/- 0.07
 _D0_M = 1864.84 * MeV  # +/- 0.05
 _DP_M = 1869.66 * MeV  # +/- 0.05
 _DS_M = 1968.35 * MeV  # +/- 0.07