theoSysW = Systematic("theoSysW", configMgr.weights, 1.0 + theoSysWNumber, 1.0 - theoSysWNumber, "user", "userOverallSys") #diboson theoSysDiboson = Systematic("theoSysDiboson", configMgr.weights, 1.5, 0.5, "user", "userOverallSys") #photon systematics in SR for Z gammaToZSyst = Systematic("gammaToZSyst", configMgr.weights, 1.25, 0.75, "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("Top", kGreen - 9) topSample.setTreeName("Top_SRAll") topSample.setNormFactor("mu_Top", 1., 0., 500.) topSample.setFileList(topFiles) topSample.setStatConfig(useStat) qcdSample = Sample("MCMultijet", kOrange + 2) qcdSample.setTreeName("QCD_SRAll") qcdSample.setNormFactor("mu_MCMultijet", 1., 0., 500.) qcdSample.setFileList(qcdFiles) qcdSample.setStatConfig(useStat)
theoSysTop = Systematic("theoSysTop", configMgr.weights, 1.+theoSysTopNumber,1.-theoSysTopNumber, "user","userOverallSys") theoSysW = Systematic("theoSysW", configMgr.weights, 1.0+theoSysWNumber,1.0-theoSysWNumber, "user","userOverallSys") #diboson theoSysDiboson = Systematic("theoSysDiboson", configMgr.weights, 1.5,0.5, "user","userOverallSys") #photon systematics in SR for Z gammaToZSyst = Systematic("gammaToZSyst", configMgr.weights, 1.25,0.75, "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("Top",kGreen-9) topSample.setTreeName("Top_SRAll") topSample.setNormFactor("mu_Top",1.,0.,50000.) topSample.setFileList(topFiles) topSample.setStatConfig(useStat) qcdSample = Sample("MCMultijet",kOrange+2) qcdSample.setTreeName("QCD_SRAll") qcdSample.setNormFactor("mu_MCMultijet",1.,0.,500.) qcdSample.setFileList(qcdFiles) qcdSample.setStatConfig(useStat)
log.info("Wait 3 seconds for you to panic if these settings are wrong") wait(3) log.info("No panicking detected, continuing...") ####################################################################### # List of samples and their plotting colours ####################################################################### #-------------------------- # Diboson #-------------------------- # NB: note that theoSys on diboson are applied on the level of the region definitions, # since we have one for the SR and one for the CR dibosonSample = Sample(zlFitterConfig.dibosonSampleName, kRed+3) dibosonSample.setTreeName("Diboson_SRAll") dibosonSample.setFileList(dibosonFiles) dibosonSample.setStatConfig(zlFitterConfig.useStat) #-------------------------- # QCD #-------------------------- qcdSample = Sample(zlFitterConfig.qcdSampleName, kOrange+2) if zlFitterConfig.useDDQCDsample:#normWeight is 0 => remove it qcdSample.setTreeName("Data_SRAll") else : qcdSample.setTreeName("QCD_SRAll") qcdSample.setNormFactor("mu_"+zlFitterConfig.qcdSampleName, 1., 0., 50000000.) qcdSample.setFileList(qcdFiles) qcdSample.setStatConfig(zlFitterConfig.useStat)
"weight", "overallNormHistoSys") # QCD theoSysQCD = Systematic("theoSysQCD", configMgr.weights, 1.0 + theoSysQCDNumber,1.0-theoSysQCDNumber, "user", "userOverallSys") QCDGausSys = Systematic("QCDGausSys", "", "_ghi", "_glo", "tree", "overallNormHistoSys") QCDTailSys = Systematic("QCDTailSys", "", "_thi", "_tlo", "tree", "overallNormHistoSys") # Diboson theoSysDiboson = Systematic("theoSysDiboson", configMgr.weights, 1.5, 0.5, "user", "userOverallSys") #------------------------------------------- # List of samples and their plotting colours #------------------------------------------- # Diboson dibosonSample = Sample("Diboson", kRed+3) dibosonSample.setTreeName("Diboson_SRAll") dibosonSample.setFileList(dibosonFiles) dibosonSample.setStatConfig(useStat) # Top topSample = Sample("ttbar", kGreen-9) topSample.setTreeName("Top_SRAll") topSample.setNormFactor("mu_Top", 1., 0., 500.) topSample.setFileList(topFiles) topSample.setStatConfig(useStat) if useTheoSys: topSample.addSystematic(theoSysTop) if useSyst : topSample.addSystematic(pileup)