#!/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,