Draft: [QEE] Wcharm line

NB: Closure

CLOSED. Superseded by: !3217 (merged)

See: discussion/reasoning !3152 (comment 7774158)

Development of asymmetric dijet lines targeting W->Charm production, but also general use for SVTag + Inclusive line for other future physics.

Specific line rates:

---

MinBias efficiency:

LAZY_AND: Hlt2QEE_DiJetIncSVTag_pT10M40dR20Full  #=25000   Sum=31 Eff=|(0.1240000 +- 0.0222572)%|

Signal MC efficiency:

LAZY_AND: Hlt2QEE_DiJetIncSVTag_pT10M40dR20Full  #=5843    Sum=278 Eff=|( 4.757830 +- 0.278485)%|

pT, di-jet mass, and dR tuning: WcharmInfo.txt

With ALL Changes:

Ran over 25000 events.
Fullstream filesize = 20 MB. Fullstream Events fired = 157

Fullstream Rate = 1.14 MHz * 157 / 25,000 = 7.159 kHz

Fullstream BW = 1.14 MHz * 20 MB / 25,000 = 912 MB/s

Turbo filesize = 8.1 MB. 1758 events passed (overlap 58)

**Note: ** Rates here will have non-negligible overlap with di-jet lines, and likely cannot be reasonably measured independently

With ONLY PV matching on all di-jet lines: Ran over 5000 events.
Fullstream filesize = 3.1 MB. Fullstream Events fired = 27
Fullstream Rate = 1.14 MHz * 27 / 5,000 = 6.156 kHz
Fullstream BW = 1.14 MHz * 3.1 MB / 5,000 = 706.8 MB/s

Turbo filesize = 1.7 MB. 346 events passed (overlap 1)

Without Changes:

Ran over 5000 events.
Fullstream filesize = 3.6 MB. Fullstream Events fired = 28
Fullstream Rate = 1.14 MHz * 28 / 5,000 = 6.384 kHz
Fullstream BW = 1.14 MHz * 3.6 MB / 5,000 = 843.6 MB/s

Turbo filesize = 1.7 MB. 346 events passed (overlap 1)

Hlt/Hlt2Conf/python/Hlt2Conf/lines/qee/jets.py

@@ -12,7 +12,7 @@ import Functors as F
 from GaudiKernel.SystemOfUnits import GeV
 
 from PyConf import configurable
-
+from math import sqrt
 from Moore.config import register_line_builder
 from Moore.lines import Hlt2Line
 
@@ -56,25 +56,39 @@ def make_jets(name='SimpleJets_{hash}',
 
 
 @configurable
-def make_dijets(tagpair=(None, None), prod_pt_min=10 * GeV, min_dphi=0.0):
+def make_dijets(tagpair=(None, None), min_dijet_mass=0 * GeV, prod_pt_min=10 * GeV, min_dphi=0.0, max_dR = 100.0):
     """Make two-jet combinations
     """
-    jets = make_jets(pt_min=prod_pt_min, tags=tagpair[0])
+    jet1 = make_jets(pt_min=prod_pt_min, tags=tagpair[0])
+    jet2 = make_jets(pt_min=prod_pt_min, tags=tagpair[1])
     tagfilt0 = F.ALL
     tagfilt1 = F.ALL
 
-    if (min_dphi > 0.0):
+
+    if (min_dphi > 0.0) and (max_dR < 100.0):
         delta = F.ADJUST_ANGLE @ (F.CHILD(1, F.PHI) - F.CHILD(2, F.PHI))
-        combination_code = F.require_all(tagfilt0, tagfilt1,
+        dR = ( F.SQRT @ ( (F.CHILD(1, F.ETA) - F.CHILD(2, F.ETA)) ** 2 + (F.CHILD(1, F.PHI) - F.CHILD(2, F.PHI)) ** 2) ) 
+        combination_code = F.require_all(tagfilt0, tagfilt1, F.MASS > min_dijet_mass,
+                                         F.ABS @ delta > min_dphi,
+                                         F.ABS @ dR < max_dR)
+    elif (min_dphi > 0.0):
+        delta = F.ADJUST_ANGLE @ (F.CHILD(1, F.PHI) - F.CHILD(2, F.PHI))
+        combination_code = F.require_all(tagfilt0, tagfilt1, F.MASS > min_dijet_mass,
                                          F.ABS @ delta > min_dphi)
+    elif (max_dR < 100.0):
+        dR = ( F.SQRT @ ( (F.CHILD(1, F.ETA) - F.CHILD(2, F.ETA)) ** 2 + (F.CHILD(1, F.PHI) - F.CHILD(2, F.PHI)) ** 2) ) 
+        combination_code = F.require_all(tagfilt0, tagfilt1, F.MASS > min_dijet_mass, F.CHILD(1, F.OWNPV) == F.CHILD(2, F.OWNPV),
+                                         F.ABS @ dR < max_dR)
+
     else:
-        combination_code = F.require_all(tagfilt0, tagfilt1)
+        combination_code = F.require_all(tagfilt0, tagfilt1, F.MASS > min_dijet_mass, F.CHILD(1, F.OWNPV) == F.CHILD(2, F.OWNPV))
 
     return ParticleCombiner(
-        Inputs=[jets, jets],
+        Inputs=[jet1, jet2],
         DecayDescriptor="CLUSjet -> CELLjet CELLjet",
         CombinationCut=combination_code,
         CompositeCut=F.ALL,
+        AllowDiffInputsForSameIDChildren = True,
         ParticleCombiner="ParticleAdder")
 
 
@@ -207,6 +221,22 @@ def diSVTagJet35GeV_line(name='Hlt2QEE_DiSVTagJet35GeVFull',
         calo_clusters=True,
         persistreco=persistreco)
 
+############ Asymmetric dijet lines ################
+@register_line_builder(all_lines)
+@configurable
+def DiJetIncSVTag_line(name='Hlt2QEE_DiJetIncSVTag_pT10M40dR20Full',
+                         prescale=1,
+                         persistreco=True):
+    jets = make_dijets(
+        tagpair=(None,'SV'), min_dijet_mass=40 * GeV, prod_pt_min=20 * GeV, min_dphi=0.0)
+    return Hlt2Line(
+        name=name,
+        algs=upfront_reconstruction() + [require_pvs(make_pvs()), jets],
+        prescale=prescale,
+        calo_clusters=True,
+        persistreco=persistreco)
+
+
 
 ############ Topo Tag ####################