Update all files

Datacards/python/writeDatacard_SUSYNano19.py

@@ -132,7 +132,10 @@ def sumBkgYields(process, signal, cr_description, yields_dict):
 		crunit = yields_dict[crproc+'_gjets'][cr][0]
 		crunit+=yields[crproc+'_back'][cr][0]
 	#print("crdata: %f, srunit: %f, crunit: %f" %(crdata, srunit, crunit))
-	total += crdata*srunit/crunit
+	if 'znunu' in process:
+		total += srunit*crunit
+	else:
+		total += crdata*srunit/crunit
 	#print("total: %f" %(total))
     return total
 
@@ -200,15 +203,11 @@ def readUncs():
                     if "up" in uncname: 
 			if "nan" in uncval:
 				uncval = 2
-			elif 'znunu' in proc_str and float(uncval) == 0:
-				uncval = 1
 			elif float(uncval) <= 0:
 				uncval = 0.001
 			unc_up = Uncertainty(uncname.strip("up"), unctype, uncval)
                     elif "down" in uncname: 
-			if 'znunu' in proc_str and float(uncval) == 0:
-				uncval = 1
-			elif uncval == "2" or "nan" in uncval or float(uncval) <= 0:
+			if uncval == "2" or "nan" in uncval or float(uncval) <= 0:
 				uncval = 0.001
 			if (unc_up.value > 1 and float(uncval) > 1) or (unc_up.value < 1 and float(uncval) < 1):
 				uncavg = averageUnc(unc_up.value, float(uncval))			

PlotsSMS/config/T1tttt_SUS16005.cfg

@@ -0,0 +1,18 @@
+# AVAILABLE COLORS:
+# kMagenta
+# kBlue
+# kOrange
+# kRed
+#####################################################
+#FORMAT: input root histo-name line-color area-color  
+#####################################################
+HISTOGRAM limit_scan_T1tttt.root hXsec_obs_corr
+EXPECTED limit_scan_T1tttt.root graph_smoothed_Exp graph_smoothed_ExpP graph_smoothed_ExpM kRed kOrange
+OBSERVED limit_scan_T1tttt.root graph_smoothed_Obs graph_smoothed_ObsP graph_smoothed_ObsM kBlack kGray
+# Preliminary Simulation or leave empty
+PRELIMINARY Preliminary
+# Lumi in fb 
+LUMI 2.3
+# Beam energy in TeV
+ENERGY 13
+#

PlotsSMS/python/sms.py

@@ -21,10 +21,10 @@ def T2tt(self):
         self.label = "pp #rightarrow #tilde{t}_{1} #bar{#tilde{t}}_{1}, #tilde{t}_{1} #rightarrow t #tilde{#chi}^{0}_{1}"
         # scan range to plot
         self.Xmin = 150.
-        self.Xmax = 2000.
+        self.Xmax = 1500.
         self.Ymin = 0.
-        self.Ymax = 1100.
-        self.Zmin = 0.000001
+        self.Ymax = 1000.
+        self.Zmin = 0.0001
         self.Zmax = 100.
         # produce sparticle
         self.sParticle = "m_{#tilde{t}_{1}} [GeV]"