EOY JMR and JMS Update

python/postprocessing/modules/common/pdfWeightProducer.py

@@ -0,0 +1,41 @@
+import ROOT
+import math, os,re
+import numpy as np
+ROOT.PyConfig.IgnoreCommandLineOptions = True
+
+from PhysicsTools.NanoAODTools.postprocessing.framework.datamodel import Collection, Object
+from PhysicsTools.NanoAODTools.postprocessing.framework.eventloop import Module
+
+class pdfWeightProducer(Module):
+    #def __init__(self):
+
+    #def beginJob(self):
+
+    #def endJob(self):
+
+    def beginFile(self, inputFile, outputFile, inputTree, wrappedOutputTree):
+        self.out = wrappedOutputTree
+        self.out.branch("LHEPdfWeightUp", "F")
+        self.out.branch("LHEPdfWeightDown", "F")
+
+    def endFile(self, inputFile, outputFile, inputTree, wrappedOutputTree):
+        pass
+
+    def analyze(self, event):
+        """process event, return True (go to next module) or False (fail, go to next event)"""
+
+        pdfWeightUp = -99.
+        pdfWeightDown = -99.
+
+        if hasattr(event,"LHEPdfWeight"):
+            self.pdfWeightList = [  i for i in event.LHEPdfWeight ]
+
+            if len(self.pdfWeightList) > 1:
+                w = np.std(self.pdfWeightList)
+                pdfWeightUp = 1. + abs(w)
+                pdfWeightDown = 1. - abs(w)
+
+        self.out.fillBranch("LHEPdfWeightUp", pdfWeightUp)
+        self.out.fillBranch("LHEPdfWeightDown", pdfWeightDown)
+
+        return True

python/postprocessing/modules/jme/jetmetHelperRun2.py

@@ -86,9 +86,8 @@
 #nominal, up, down
 jmrValues = {
     '2016': [1.0, 1.2, 0.8],
-    '2017': [1.09, 1.14, 1.04],
-    # Use 2017 values for 2018 until 2018 are released
-    '2018': [1.09, 1.14, 1.04],
+    '2017': [0.99, 1.032, 0.948],
+    '2018': [1.108, 1.142, 1.074],
     'UL2016_preVFP': [1.00, 1.00, 1.00],  # placeholder
     'UL2016': [1.00, 1.00, 1.00],  # placeholder
     'UL2017': [1.00, 1.00, 1.00],  # placeholder
@@ -97,12 +96,10 @@
 
 # jet mass scale
 # W-tagging PUPPI softdrop JMS values: https://twiki.cern.ch/twiki/bin/view/CMS/JetWtagging
-# 2016 values
 jmsValues = {
     '2016': [1.00, 0.9906, 1.0094],  # nominal, down, up
-    '2017': [0.982, 0.978, 0.986],
-    # Use 2017 values for 2018 until 2018 are released
-    '2018': [0.982, 0.978, 0.986],
+    '2017': [1.014, 1.011, 1.017],
+    '2018': [0.999, 0.997, 1.001],
     'UL2016_preVFP': [1.000, 1.000, 1.000],  # placeholder
     'UL2016': [1.000, 1.000, 1.000],  # placeholder
     'UL2017': [1.000, 1.000, 1.000],  # placeholder