if options.small:
prefixes.append('small')
prefix = '_'.join(prefixes)
if prefix!='':
prefix+='_'
print 'maps', options.maps,'prefix',options.prefix, 'input', options.input, 'mode', options.mode
if options.maps=='all':
maps = allMaps
else:
exec("maps = [" +options.maps+ "]")
if options.mode=='all':
modes=["ngoodVertices", "sumPt", "multiplicity"]
else:
modes=[options.mode ]
if os.path.isdir(options.input):
c = getChain(options.input, maxN=1) if options.small else getChain(options.input)
else:
c=None
def getLinSquStr(f):
return "10^{-6}#upoint ("+str(round(10**6*f.GetParameter(0),1))+'#pm '+str(round(10**6*abs(f.GetParError(0)),1))+") #upoint n^{2}+10^{-3}#upoint("+str(round(10**3*f.GetParameter(1),1))+'#pm '+str(round(10**3*abs(f.GetParError(1)),1))+") #upoint n"
def getLinSquStrNoFac(f):
return "("+str(round(f.GetParameter(0),1))+'#pm '+str(round(abs(f.GetParError(0)),1))+") #upoint n^{2}+("+str(round(f.GetParameter(1),1))+'#pm '+str(round(abs(f.GetParError(1)),1))+") #upoint n"
def getSquStr(f):
return "10^{-6}#upoint ("+str(round(10**6*f.GetParameter(0),1))+'#pm '+str(round(10**6*abs(f.GetParError(0)),1))+") #upoint n^{2}"
def getPropStr(f):
return "10^{-3} #upoint ("+str(round(10**3*f.GetParameter(0),1))+'#pm '+str(round(10**3*abs(f.GetParError(0)),1))+") #upoint n"
def getLinStr(f):
return '10^{-2} #upoint ('+str(round(10**2*f.GetParameter(0),1))+'#pm '+str(round(10**2*abs(f.GetParError(0)),1))+") + 10^{-3}#upoint("+str(round(10**3*f.GetParameter(1),1))+'#pm '+str(round(10**3*abs(f.GetParError(1)),1))+") #upoint n"
for mode in modes:
if mode=="ngoodVertices":
import ROOT
from Workspace.HEPHYPythonTools.helpers import getObjFromFile, getChain, getChunks, getObjDict, getVarValue
from Workspace.HEPHYPythonTools.xsec import xsec
ROOT.gROOT.LoadMacro("$CMSSW_BASE/src/Workspace/HEPHYPythonTools/scripts/root/tdrstyle.C")
ROOT.setTDRStyle()
ROOT.gROOT.ProcessLine(".L $CMSSW_BASE/src/Workspace/HEPHYPythonTools/scripts/root/useNiceColorPalette.C+")
ROOT.useNiceColorPalette()
from Workspace.RA4Analysis.cmgTuples_Spring15_150809 import *
lumi=100.
small = False
samples=[
{"name":"DY", "bins":[DYJetsToLL_M50_25ns_CERNRerun], "legendText":"Drell-Yan"}
]
maxN = 10 if small else -1
for s in samples:
totalYield=0
for b in s["bins"]:
chunks, sumWeight = getChunks(b, maxN=maxN)
# print "Chunks:" , chunks
lumiScale = xsec[b['dbsName']]*lumi/float(sumWeight)
b["lumiScale"] = lumiScale
b["chain"] = getChain(chunks, histname="", treeName = b["treeName"])
print b["name"],"xsec",xsec[b['dbsName']],"sumWeight",sumWeight
DYJetsToLL_M50_25ns_CERNRerun['chain'].Draw('Jet_pt/Jet_rawPt:Jet_eta','','COLZ')