qpython 31.08.2016
"""
import facoLib as fl
msg = "testing new function ...\n"
print msg
print fl.combinedRelUncertainty(100, [40., 80], True)
print msg
print fl.combinedUncertainty([40., 80], True)
print msg
x = 4
y = 6
z = 2.1
myDoubleArray = [[0, 1, 2], [x, y, z]]
myHeader = ["this", "stuff", "is great"]
print fl.makeTable("testOutput", myDoubleArray, myHeader, True, "any caption", True)
print msg
print fl.getDictFromJson("WW", "", True)
print msg
print fl.floatToPercent(0.12)
# define the bound of the Signal region
#LowBounds=[0.001,0.02,0.05]
bound1 = 0.01
# template histogram that contains one single bin
hist=rt.TH1D("template","template",1,0,1,)
# hard coded lumi taken from doubleEG. Note that the exact value does not influence the results as we are calculating ratio.
lumivalue = 2610.486
# dictionary for the cross section which depends on the samples and on the systShift
XSWeight_dict = fl.getDictFromJson("CrossSection", "", True)
samples_dict={'NP_overlay_stopTobl_m500_Ctau10': 'central', 'TTJets_SingleLeptFromT': 'central', 'ZZTo2L2Nu': 'central', 'TTJets_Dilept': 'central', 'SingleTop_tW': 'central', 'TTJets_SingleLeptFromTbar': 'central', 'ZG': 'central', 'WJetsToLNu': 'central', 'SingleTop_tbarW': 'central', 'WGToLNuG': 'central', 'ZZTo4Lpowheg': 'central', 'WZTo3LNu': 'central', 'WWTo2l2Nu': 'central', 'WWToLNuQQ': 'central'}
# loop over the channel (lepton in final statue)
for chan in channels:
isElEl=False
isMuMu=False
doubleArray=[]
# get the xmlfile that corresponds to the channel
if "MuMu" in chan:
isMuMu=True
tree = ET.ElementTree(file='../config/MuMuV4.xml')
FinalState="At least two muons"
# dict with the cst systemtic uncertainties
dict_cstSystType={"lumi":0.05, "trigger_sf_electron":0.20, "trigger_sf_muon":0.10, "trk_sf_electron":0.24, "trk_sf_muon":0.12 }
# lepton dependent syst type
muonSystTypes=["XSWeight", "evt_puSF","sf_iso_muon","sf_id_muon","lumi","trigger_sf_muon", "trk_sf_muon"]
electronSystTypes=["XSWeight", "evt_puSF","sf_reco_electron","sf_id_electron","trigger_sf_electron", "trk_sf_electron"]
# loop over the channels
for chan in ["ElEl", "MuMu"]:
# dictionary for the cross section which depends on the samples and on the systShift
Yield_dict = fl.getDictFromJson(chan, "", True)
print "the merged dict is \n"
print Yield_dict
if "ElEl" in chan:
isElEl=True
tree = ET.ElementTree(file='../config/ElElV4.xml')
FinalState="At least two electrons"
print FinalState
systTypes=list(electronSystTypes)
elif "MuMu" in chan:
isMuMu=True
tree = ET.ElementTree(file='../config/MuMuV4.xml')
FinalState="At least two muons"
