v14

Author: rkansal47

src/boostedhh/processors/objects.py

@@ -28,3 +28,34 @@ def jetid_v12(jets: ak.Array) -> ak.Array:
     ) | ((np.abs(jets.eta) > 2.7) & jetidtight)
 
     return jetidtight, jetidtightlepveto
+
+
+def jetid_v14(jets: ak.Array) -> tuple[ak.Array, ak.Array]:
+    """
+    Jet ID fix for NanoAOD v14 copying
+    # https://gitlab.cern.ch/cms-jetmet/coordination/coordination/-/issues/117#note_8880788
+    """
+
+    jetidtight = (
+        (
+            (np.abs(jets.eta) <= 2.6)
+            & (jets.neHEF < 0.99)
+            & (jets.neEmEF < 0.9)
+            & ((jets.chMultiplicity + jets.neMultiplicity) > 1)
+            & (jets.chHEF > 0.01)
+            & (jets.chMultiplicity > 0)
+        )
+        | (
+            ((np.abs(jets.eta) > 2.6) & (np.abs(jets.eta) <= 2.7))
+            & (jets.neHEF < 0.90)
+            & (jets.neEmEF < 0.99)
+        )
+        | (((np.abs(jets.eta) > 2.7) & (np.abs(jets.eta) <= 3.0)) & (jets.neHEF < 0.99))
+        | ((np.abs(jets.eta) > 3.0) & (jets.neMultiplicity >= 2) & (jets.neEmEF < 0.4))
+    )
+
+    jetidtightlepveto = (
+        (np.abs(jets.eta) <= 2.7) & jetidtight & (jets.muEF < 0.8) & (jets.chEmEF < 0.8)
+    ) | ((np.abs(jets.eta) > 2.7) & jetidtight)
+
+    return jetidtight, jetidtightlepveto

src/boostedhh/xsecs.py

@@ -292,6 +292,14 @@
 xsecs["DYto2L-2Jets_MLL-50_1J"] = 1017.0
 xsecs["DYto2L-2Jets_MLL-50_2J"] = 385.5
 
+xsecs["DYto2L-4Jets_MLL-50to120_HT-40to70"] = 316.7
+xsecs["DYto2L-4Jets_MLL-50to120_HT-70to100"] = 140.1
+xsecs["DYto2L-4Jets_MLL-50to120_HT-100to400"] = 179.6
+xsecs["DYto2L-4Jets_MLL-50to120_HT-400to800"] = 6.742
+xsecs["DYto2L-4Jets_MLL-50to120_HT-800to1500"] = 0.693
+xsecs["DYto2L-4Jets_MLL-50to120_HT-1500to2500"] = 0.05047
+xsecs["DYto2L-4Jets_MLL-50to120_HT-2500"] = 0.00346
+
 ########################################################
 # Cross Sections for 13 TeV.
 ########################################################