Reduce the rates of double-charm lines in B&Q wg

As mentioned in #410, the rates of double-charm lines in B&Q are found to be very high. The rate is reduced to about 10% followed the selection in Charm wg. 50k MB MC samples are used to test the rate of the previous and updated double-charm line.

updated double-charm line:

Hlt2BandQ_DoubleCharmOppositeSign_LineDecisionWithOutput      #=50000   Sum=48          Eff=|(0.09600000 +- 0.0138498)%|
Hlt2BandQ_DoubleCharmSameSign_LineDecisionWithOutput          #=50000   Sum=56          Eff=|(0.1120000 +- 0.0149582)%|

previous double-charm line:

Hlt2BandQ_DoubleCharmOppositeSign_LineDecisionWithOutput      #=50000   Sum=533         Eff=|( 1.066000 +- 0.0459268)%|
Hlt2BandQ_DoubleCharmSameSign_LineDecisionWithOutput          #=50000   Sum=569         Eff=|( 1.138000 +- 0.0474352)%|

Hlt/Hlt2Conf/python/Hlt2Conf/lines/bandq/builders/c_hadrons.py

@@ -11,6 +11,7 @@
 """
 Definition of B&Q c-hadrons
 """
+import math
 import Functors as F
 from Functors.math import in_range
 
@@ -195,3 +196,48 @@ def make_omegac(particles,
         am_max=am_max,
         m_min=m_min,
         m_max=m_max)
+
+
+@configurable
+def make_tightomegac(particles,
+                     descriptor,
+                     name="bandq_tightOmegac",
+                     am_min=2630 * MeV,
+                     am_max=2750 * MeV,
+                     comb_sum_pt_min=3000.0 * MeV,
+                     comb_pt_min_at_least_one=1000.0 * MeV,
+                     comb_pt_min_at_least_two=400.0 * MeV,
+                     comb_mipchi2_min_at_least_one=8.0,
+                     comb_mipchi2_min_at_least_two=6.0,
+                     pt_min=1000.0 * MeV,
+                     vtx_chi2pdof_max=10,
+                     bpvltime_min=0.2 * picosecond,
+                     apt_min=950 * MeV,
+                     bpvdira_min=math.cos(0.01),
+                     achi2_doca_max=25,
+                     m_min=2640 * MeV,
+                     m_max=2740 * MeV):
+
+    pvs = make_pvs()
+
+    combination_code = require_all(
+        in_range(am_min, F.MASS, am_max), F.PT > apt_min,
+        F.DOCACHI2(1, 2) < achi2_doca_max,
+        F.DOCACHI2(1, 3) < achi2_doca_max,
+        F.DOCACHI2(2, 3) < achi2_doca_max,
+        F.SUM(F.PT) > comb_sum_pt_min,
+        F.SUM(F.PT > comb_pt_min_at_least_one) >= 1,
+        F.SUM(F.PT > comb_pt_min_at_least_two) >= 2,
+        F.SUM(F.MINIPCHI2(pvs) > comb_mipchi2_min_at_least_one) >= 1,
+        F.SUM(F.MINIPCHI2(pvs) > comb_mipchi2_min_at_least_two) >= 2)
+    vertex_code = require_all(
+        in_range(m_min, F.MASS, m_max), F.PT > pt_min,
+        F.CHI2DOF < vtx_chi2pdof_max,
+        F.BPVDIRA(pvs) > bpvdira_min,
+        F.BPVLTIME(pvs) > bpvltime_min)
+    return ParticleCombiner(
+        name=name,
+        Inputs=particles,
+        DecayDescriptor=descriptor,
+        CombinationCut=combination_code,
+        CompositeCut=vertex_code)