parser.add_option('--nameSuffix', dest="histNameSuffix", default="")
o, a = parser.parse_args()
###########
## Setup ##
###########
if o.year == "2016": lumi = "24.3 fb^{-1}"
if path.exists(o.limitFile):
limitFile = TFile.Open(o.limitFile, "UPDATE")
else:
limitFile = TFile.Open(o.limitFile, "RECREATE")
inDir = "hists_" + o.year + "_" + o.weights
files = rootFiles.getFiles(o.iteration, o.nTuple, inDir, o.year)
shapeVariation = "_CR"
#shapeVariation="_CR"
if o.doSyst:
filesCRw = rootFiles.getFiles(o.iteration, o.nTuple,
inDir + shapeVariation, o.year)
muFile = "XhhResolved/data/mu_qcd_FourTag_" + o.weights + "_" + o.year + "_" + o.variation + "_" + o.iteration + ".txt"
mu_qcd_dict = read_mu_qcd_file(muFile)
if o.doSyst:
muFileCRw = "XhhResolved/data/mu_qcd_FourTag_" + o.weights + shapeVariation + "_" + o.year + "_" + o.variation + "_" + o.iteration + ".txt"
mu_qcd_dictCRw = read_mu_qcd_file(muFileCRw)
useAllhad2bShape = True
#get files
blind = True
if o.year == "2016": lumi = "24.3 fb^{-1}"
iteration = o.iteration
if o.threeTag:
tag = "Three"
else:
tag = "Four"
if os.path.exists(o.limitFile):
limitFile = ROOT.TFile.Open(o.limitFile, "UPDATE")
else:
limitFile = ROOT.TFile.Open(o.limitFile, "RECREATE")
files = rootFiles.getFiles(iteration,o.nTuple, o.inDir, o.year)
muFile = "XhhResolved/data/mu_qcd_"+tag+"Tag_"+o.weights+"_"+o.year+"_"+o.variation+"_"+iteration+".txt"
mu_qcd_dict = read_mu_qcd_file(muFile)
rebins={}
#rebins["m4j_cor_l"] = [100,120,140,160,180,200,220,240,260,280,300,320,340,360,380,400,420,440,460,480,
# 500,520,540,560,580,600,640,680,720,820,900,1000,1500]
#rebins["m4j_cor_l"] = [100,250,270,290,330,360,400,450,510,580,680,800,1000,1500]
rebins["m4j_cor_l"] = [100,120,140,160,180,200,220,240,260,280,300,320,340,360,380,400,420,440,460,480,
500,520,540,560,580,600,640,680,720,820,900,1000,1200]
# rebins["m4j_cor_l"] = [100,250,270,350,430,
# 510,590,730,1200]
#rebins["m4j_cor_l"] = [100,250,270,290,330,360,450,580,800,1500]
import sys
sys.path.insert(0, 'XhhResolved/plotting/')
import rootFiles
import ROOT
iteration="7"
nTuple="02-03-04"
year="2015"
inDir="hists_"+year+"_"+weightSet
rootFiles2015 = rootFiles.getFiles(iteration,nTuple,inDir,year)
year="2016"
inDir="hists_"+year+"_"+weightSet
rootFiles2016 = rootFiles.getFiles(iteration,nTuple,inDir,year)
class sample:
def __init__(self, path, tag, scale = 1, error=0):
f = ROOT.TFile(path,"read")
self.path = path
self.n_pm_e = {}
self.n = {}
self.e = {}
self.cuts = ["Inclusive","PassHCdEta","PassAllhadVeto"]
self.regions = ["Sideband","Control","Signal"]
for cut in self.cuts:
def histCleanup(h):
h.SetStats(0)
h.GetXaxis().SetLabelFont(43)
h.GetXaxis().SetLabelSize(16)
h.GetYaxis().SetLabelFont(43)
h.GetYaxis().SetLabelSize(16)
h.GetXaxis().SetTitleFont(43)
h.GetXaxis().SetTitleSize(20)
h.GetYaxis().SetTitleFont(43)
h.GetYaxis().SetTitleSize(20)
h.GetXaxis().SetTitleOffset(2.8)
h.GetYaxis().SetTitleOffset(1.6)
files = rootFiles.getFiles("0", o.nTuple, "hists", o.year)
masses = [
"M300", "M400", "M500", "M600", "M700", "M800", "M900", "M1000", "M1100",
"M1200", "SMNR"
]
#masses=["M300"]
f = {}
c = {}
if o.year == "2015":
L1 = [
"passL1", "L1_J100", "L1_4J20", "L1_3J25.0ETA23", "L1_4J15.0ETA25",
"L1_J75_3J20"
]
parser.add_option('--inDir', dest="inDir", default="hists", help="")
parser.add_option('--out', dest="out", default="", help="")
parser.add_option('-i', '--iter', dest="iteration", default="0", help="")
parser.add_option('-v', '--nTuple', dest="nTuple", default="01-01-01", help="")
parser.add_option('-y', '--year', dest="year", default="2015", help="")
o, a = parser.parse_args()
iteration = o.iteration
nTuple = o.nTuple
inDir = o.inDir
year = o.year
out = o.out
files = rootFiles.getFiles(iteration, nTuple, inDir, year)
L1s = [
'L1_4J15.0ETA25',
'L1_J75_3J20',
'L1_J100',
'L1',
]
HLTs = [
'HLT_2j35_bmv2c2060_split_2j35_L14J15.0ETA25',
'HLT_j100_2j55_bmv2c2060_split',
'HLT_j225_bmv2c2060_split',
'HLT',
'HLT_SF_up',
'HLT_SF_down',
#
condor = True
script = "config_hh4bCombUnique.py"
nTuple = "02-03-04"
tagger = "MV2c10"
var = "Nominal"
#ttbarVariation="ttbar_hard2b_xwt"
ttbarVariation=""
weightRegion="Sideband"
weightSet = "redo"+("_"+ttbarVariation if ttbarVariation else "")+("_CR" if weightRegion == "Control" else "")
#weightSet = "noLowMDR"+("_"+ttbarVariation if ttbarVariation else "")+("_CR" if weightRegion == "Control" else "")
year = "2016"
outdir = "hists_"+year+"_"+weightSet
files = rootFiles.getFiles("0",nTuple, outdir, year)
injectFile = None
injectMu = None
# injectFile = files["H280"]
# injectMu = "0.04"
iterations = [int(i) for i in o.iterations.split(",")] if o.iterations else []
fast_iterations = [int(i) for i in o. fast.split(",")] if o.fast else []
iterations = sorted(iterations+fast_iterations)
plot_iterations = iterations
applyWeight = True
#limitVariable = "m4j_cor_f"
limitVariable = "m4j_cor_Z_f"
#limitCut = "PassAllhadVeto"
limitCut = "PassHCdEta"
#limitCut = "Pass_ggVeto"