# Define cuts configMgr.cutsDict["CR"] = "1." configMgr.cutsDict["SR"] = "1." # Define weights configMgr.weights = "1." # Define samples bkgSample = Sample("Bkg",kGreen-9) bkgSample.setNormByTheory(True) bkgSample.buildHisto(nBkgCR,"CR","cuts",0.5) bkgSample.buildHisto(nBkgSR,"SR","cuts",0.5) bkgSample.addSystematic(bg1xsec) ddSample = Sample("DataDriven",kGreen+2) ddSample.addShapeFactor("DDShape") sigSample = Sample("Sig",kPink) sigSample.setNormFactor("mu_Sig",1.,0.2,1.5) sigSample.buildHisto(nSigSR,"SR","cuts",0.5) sigSample.setNormByTheory(True) sigSample.addSystematic(sigxsec) dataSample = Sample("Data",kBlack) dataSample.setData() dataSample.buildHisto(nDataCR,"CR","cuts",0.5) dataSample.buildHisto(nDataSR,"SR","cuts",0.5) # Define top-level ana = configMgr.addFitConfig("SPlusB") ana.addSamples([bkgSample,ddSample,dataSample])
#-------------------------- # Data #-------------------------- dataSample = Sample("Data", kBlack) dataSample.setTreeName("Data_SRAll") dataSample.setData() dataSample.setFileList(dataFiles) ####################################################################### # Shape Factor ####################################################################### if zlFitterConfig.useShapeFit and zlFitterConfig.useShapeFactor: topSample.addShapeFactor("topShape") wSample.addShapeFactor("wShape") zSample.addShapeFactor("zShape") gammaSample.addShapeFactor("gammaShape") qcdSample.addShapeFactor("qcdShape") ####################################################################### # Apply systematics ####################################################################### # Here be systematics. Sometime in future. # sysWeight_theoSysSigUp = myreplace(configMgr.weights, ["normWeightUp"], "normWeight") # sysWeight_theoSysSigDown = myreplace(configMgr.weights, ["normWeightDown"], "normWeight") # theoSysSig = Systematic("SigXSec", configMgr.weights, sysWeight_theoSysSigUp, sysWeight_theoSysSigDown, "weight", "overallSys")
sherpaTFSysTop = Systematic("sherpaTFSys", configMgr.weights, 1.0+TopMap['topSherpa'][level][chn],1.0-TopMap['topSherpa'][level][chn], "user","userOverallSys") #------------------------------------------- # List of samples and their plotting colours #------------------------------------------- dibosonSample = Sample("Diboson",kRed+3) dibosonSample.setTreeName("Diboson_SRAll") dibosonSample.setFileList(dibosonFiles) dibosonSample.setStatConfig(useStat) dibosonSample.addSystematic(theoSysDiboson) topSample = Sample("ttbar",kGreen-9) if useShape: topSample.addShapeFactor("topShape") topSample.setTreeName("Top_SRAll") topSample.setNormFactor("mu_Top",1.,0.,500.) topSample.setFileList(topFiles) topSample.setStatConfig(useStat) if useTheoSys: topSample.addSystematic(theoSysTop) ####topSample.addSystematic(mu1ScaleSysTop) ####topSample.addSystematic(mu2ScaleSysTop) ####topSample.addSystematic(matchScaleSysTop) topSample.setNormRegions([("CRT","meffInc"),("CRW","meffInc")]) qcdSample = Sample("Multijets",kOrange+2) qcdSample.setTreeName("QCDdd") qcdSample.setNormFactor("mu_Multijets",1.,0.,500.) qcdSample.setFileList(qcdFiles)