increase mumu mass range to include psi two s

Sorry for missing the deadline, personal emergency happened. Making draft merge request to put it on the radar.

I'd like to make possible the study of SL Bc decays to psi(2S) in addition to j/psi.

First draft of this request simply increased the mass range of Hlt2SLB_BcToJpsiMuNu_JpsiToMuMuDecision line (and fake version) to +120 MeV past psi(2S). This increases the rate of the FakeMuon line from 0.0774 kHz to 0.155 kHz with the standard Moore BW script run locally. BW increase from 0.0186GB/s to 0.0386 GB/s. The signal muon line is zero for both selections, so can't compare.

Moving to two mass ranges, the background added will be roughly (240 extra MeV range) / (589 extra MeV range) = 0.4 so depending on the exact psi(2S) signal to background ratio the overall rate would be something like 1.5-1.75 times the old ones probably.

May want to consider getting opinions on expanding this selection to other Bc to Jpsi SL lines.

Hlt/Hlt2Conf/python/Hlt2Conf/lines/semileptonic/builders/charm_hadron_builder.py

@@ -47,6 +47,8 @@ def make_Hc_to_nbody(
     descriptor,
     comb_m_min,
     comb_m_max,
+    comb_m_min2=None,
+    comb_m_max2=None,
     mother_m_min=None,
     mother_m_max=None,
     mother_pt_min=None,
@@ -77,6 +79,10 @@ def make_Hc_to_nbody(
         Lower invariant mass limit for the particle combination.
     comb_m_max :
         Upper invariant mass limit for the particle combination.
+    comb_m_min2 :
+        Lower invariant mass limit for the second particle combination.
+    comb_m_max2 :
+        Upper invariant mass limit for the second particle combination.
     comb_pt_any_min :
         Minimum pt that at least one of the particles in the combination needs to have.
     comb_pt_sum_min :
@@ -97,6 +103,10 @@ def make_Hc_to_nbody(
         Minimum distance of flight in the transverse plane.
     """
     combination_code = in_range(comb_m_min, F.MASS, comb_m_max)
+    if (comb_m_min2 is not None) and (comb_m_max2 is not None):
+        combination_code = F.require_any(
+            combination_code, in_range(comb_m_min2, F.MASS, comb_m_max2)
+        )
     if comb_pt_min is not None:
         combination_code = F.require_all(combination_code, F.PT > comb_pt_min)
     if comb_pt_any_min is not None:
@@ -408,6 +418,8 @@ def make_jpsi_tomumu(
     name="SLB_JpsiToMuMuBuilder_{hash}",
     comb_m_min=3047 * MeV,
     comb_m_max=3147 * MeV,
+    comb_m_min2=None,
+    comb_m_max2=None,
     comb_pt_any_min=1200 * MeV,
     comb_pt_sum_min=3500 * MeV,
     comb_docachi2_max=20.0,
@@ -437,6 +449,8 @@ def make_jpsi_tomumu(
         bpvfdchi2_min=bpvfdchi2_min,
         bpvdira_min=bpvdira_min,
         bpvvdz_min=bpvvdz_min,
+        comb_m_min2=comb_m_min2,
+        comb_m_max2=comb_m_max2,
         comb_pt_min=None,
         comb_pt_any_min=comb_pt_any_min,
         comb_pt_sum_min=comb_pt_sum_min,