#!/usr/bin/python
import ROOT as r
import sys
from array import array
from math import sqrt
from SCALE import *
from binning import bins2D
njets = sys.argv[1]
rebinX = bins2D().rebinX #5 #2GeV/bin
numybins = len(bins2D().absybins) - 1
nominal = "noEW"
f = r.TFile("/afs/cern.ch/user/y/yduh/work/lpcresults/%sunc/%s/pdf/%s.root" %
(njets, nominal, sys.argv[2]))
outfile = "./%s/skim_pdf_%s.root" % (njets, sys.argv[2])
outfile = r.TFile(outfile, "RECREATE")
kint = "3j_RECO/3j_Mtt_delY" if njets == '3j' else "YUKAWA_RECO/yukawa_Mtt_delY"
hnominal = [[0 for x in range(2)] for x in range(numybins)]
#print hnominal
for i in range(len(bins2D().absybins) - 1):
binpL = f.Get(kint).GetYaxis().FindFixBin(bins2D().absybins[i])
binpH = f.Get(kint).GetYaxis().FindFixBin(bins2D().absybins[i + 1])
binmH = f.Get(kint).GetYaxis().FindFixBin(-(bins2D().absybins[i]))
binmL = f.Get(kint).GetYaxis().FindFixBin(-(bins2D().absybins[i + 1]))
#print binpL, binpH, binmL, binmH
hnominal[i][0] = f.Get(kint).ProjectionX("mtt%s0" % i, binpL,
binpH).Rebin(rebinX)
rebinX.append([0, 400, 440, 480, 520, 580, 660, 780, 2000])
elif njets == '4j':
rebinX.append([0, 360, 400, 440, 480, 540, 2000])
rebinX.append([0, 400, 440, 480, 520, 580, 2000])
rebinX.append([0, 500, 560, 620, 700, 2000])
elif njets == '5j':
rebinX.append([0, 380, 420, 480, 2000])
rebinX.append([0, 420, 480, 560, 2000])
rebinX.append([0, 540, 640, 2000])
elif njets == '6j':
rebinX.append([0, 400, 2000])
rebinX.append([0, 460, 2000])
rebinX.append([0, 600, 2000])
###########################################
numxbins = len(rebinX[0]) - 1 + len(rebinX[1]) - 1 + len(rebinX[2]) - 1
numybins = len(bins2D().absybins) - 1
nominal = "noEW"
print numxbins, numybins
def prepared(fileName, scaleFactor):
f = r.TFile("/afs/cern.ch/user/y/yduh/work/lpcresults/%sunc/%s/pdf/%s" %
(njets, nominal, fileName))
outfile = "./%s/%s/skim_pdf_%s_%s.root" % ('CT14' if sys.argv[2] == 'CT14'
else 'CategorizedNNPDF', njets,
fileName.split('.')[0], nominal)
outfile = r.TFile(outfile, "RECREATE")
kint = "3j_RECO/3j_Mtt_delY" if njets == '3j' else "YUKAWA_RECO/yukawa_Mtt_delY"
hsigList = []
]
#uncUptxt = open('err/uncerUp_%s.txt' %njets)
#uncDowntxt = open('err/uncerUp_%s.txt' %njets)
#uncUp = uncUptxt.readlines()
#uncDown = uncDowntxt.readlines()
#uncer = [x.strip() for x in uncUp]
#print uncer
inpath = "/uscms/home/yiting11/nobackup/CMSSW_7_4_12_patch4/src/ttbar_preparation/WrapUp/%s/skimroot/" % njets
inpathsb = "/uscms/home/yiting11/nobackup/CMSSW_7_4_12_patch4/src/ttbar_preparation/WrapUp/%s/skimrootSB/" % njets
inpathsys = "/uscms/home/yiting11/nobackup/CMSSW_7_4_12_patch4/src/ttbar_preparation/WrapUp/%s/skimrootSYS/" % njets
outpath = "/uscms/home/yiting11/nobackup/CMSSW_7_4_12_patch4/src/ttbar_preparation/WrapUp/%s/packroot/" % njets
kint = "mtt_dely_RECO" #vararg2D().xy_savename+"_RECO"
ybins = bins2D().ybins
absybins = bins2D().absybins
xbins = bins2D.xbins_4j if njets != '3j' else bins2D().xbins_3j
if (FixBins):
xbins = [20 * i for i in range(101)]
xfixbins = bins2D().xfixbins_3j
if njets == '4j':
xfixbins = bins2D().xfixbins_4j
elif njets == '5j':
xfixbins = bins2D().xfixbins_5j
elif njets == '6j':
xfixbins = bins2D().xfixbins_6j
#print "xbins = ", bins2D().xbins_4j if njets != '3j' else bins2D().xbins_3j, " ", len(bins2D().xbins_4j)-1 if njets != '3j' else len(bins2D().xbins_3j)-1
#print "ybins = ", bins2D().ybins, " ", len(bins2D().ybins)-1
#!/usr/bin/python
import ROOT as r
import sys
from array import array
import Rebin
#import binning
from binning import bins2D
njets = sys.argv[1]
numin = 10000 #int(sys.argv[2])
xbins = bins2D().xbins_4j if njets != '3j' else bins2D().xbins_3j
kint = "mtt_dely_RECO"
GT = ['1695', '1715', '1735', '1755', 'noEW']
inpath = '/afs/cern.ch/user/y/yduh/CMSSW_7_1_5/src/ttbar_preparation/WrapUp/%s/' % njets
f = r.TFile(inpath + 'packroot/ch%s.root' % njets)
numbins = f.Get("ttsig").GetNbinsX()
fnominal = r.TFile(inpath + 'skimroot/skim_tt_PowhegP8_noEW.root', 'READ')
histnominal = Rebin.Absy(
Rebin.newRebin2D(fnominal.Get(kint), 'ttsig_temp', xbins,
bins2D().ybins), kint, xbins,
bins2D().absybins)
outfile = r.TFile('./%s/ch%s_mt.root' % (njets, njets), 'RECREATE')
for gt in GT:
f = r.TFile(inpath + 'skimroot/skim_mtop%s_tt_PowhegP8_noEW.root' % gt,