Merge pull request #18 from jmduarte/convert_jmr_values

convert JMR values

Author: farakiko

boostedhiggs/corrections.py

@@ -738,12 +738,20 @@ def pad_val(
 
 # jet mass resolution: https://twiki.cern.ch/twiki/bin/view/CMS/JetWtagging
 # nominal, down, up (these are switched in the github!!!)
-jmrValues["msoftdrop"] = {
+jmrValues["msoftdrop_twiki"] = {
     "2016": [1.0, 0.8, 1.2],
     "2017": [1.09, 1.04, 1.14],
     # Use 2017 values for 2018 until 2018 are released
     "2018": [1.09, 1.04, 1.14],
 }
+# convert from SF on sigma to value to smear mass by
+# note: clamped to 1.0 (so 2016 JMR down needs to be interpolated by combine)
+jmrValues["msoftdrop"] = {}
+for year in jmrValues["msoftdrop_twiki"]:
+    rel_res = 0.12 if year == "2016" else 0.10
+    jmrValues["msoftdrop"][year] = [
+        1 + np.sqrt(max(sf_res**2 - 1, 0)) * rel_res for sf_res in jmrValues["msoftdrop_twiki"][year]
+    ]
 
 # jet mass scale
 # W-tagging PUPPI softdrop JMS values: https://twiki.cern.ch/twiki/bin/view/CMS/JetWtagging