qFlipSample.setFileList(dataFiles)
fLepSample.setFileList(dataFiles)
dataSample.setFileList(dataFiles)
for aSample in [
wgammaSample, dibosonSample, tribosonSample, topXSample, higgsSample
]:
aSample.setNormByTheory(True)
aSample.setStatConfig(True)
for aSys in reconSys:
aSample.addSystematic(aSys)
qFlipSample.addSystematic(qFlipSys)
qFlipWeights = ("ElSF", "MuSF", "BtagSF", "weight", "pwt", "qFwt",
dataDrivenBkgScale)
qFlipSample.setWeights(qFlipWeights) #overrides configMgr.weights
qFlipSample.setNormByTheory(False)
qFlipSample.setStatConfig(True)
fLepSample.addSystematic(fLepSys)
fLepWeights = ("ElSF", "MuSF", "BtagSF", "weight", "pwt", "fLwt",
dataDrivenBkgScale)
fLepSample.setWeights(fLepWeights) #overrides configMgr.weights
fLepSample.setNormByTheory(False)
fLepSample.setStatConfig(True)
dataWeights = ("((%s)&&(%s))" % (ssUtil.sigLepCut, ssUtil.ssCut),
dataDrivenBkgScale)
dataSample.setWeights(
dataWeights
) #need to edit python/configManager.setWeightsCutsVariable() in HistFitter for this to take effect
hasHippo = doToys and 'hippo' in userOpts.gridForToys
hasSig = doToys and signame in userOpts.gridForToys
if not doToys or hasHippo or hasSig :
exclusion = configMgr.addFitConfig(tlx, 'TopLvlXML_Exclusion_%s'%s)
userPrint(" >>> !!! --- in LOOP %s"%s)
sigSample = Sample(s, kRed)
sigSample.setStatConfig(useStat)
if userOpts.splitMCSys :
sigSample.addSystematic(sysObj.AR_mcstat_SIG)
sigUncert = userOpts.sigUncert
sigSample.setFileList([userOpts.sigFile])
if sigUncert == "up":
sigSample.setWeights(("eventweight","syst_XSUP","1.0"))
elif sigUncert == "down":
sigSample.setWeights(("eventweight","syst_XSDOWN","1.0"))
sigSample.addSystematic(sysObj.AR_all_JES_MC)
sigSample.addSystematic(sysObj.AR_all_JER_MC)
# sigSample.addSystematic(sysObj.AR_all_JVF_MC)
sigSample.addSystematic(sysObj.AR_all_TES_MC)
sigSample.addSystematic(sysObj.AR_all_BJET_MC)
sigSample.addSystematic(sysObj.AR_all_CJET_MC)
sigSample.addSystematic(sysObj.AR_all_BMISTAG_MC)
sigSample.addSystematic(sysObj.AR_all_EER_MC)
sigSample.addSystematic(sysObj.AR_all_EESLOW_MC)
sigSample.addSystematic(sysObj.AR_all_EESMAT_MC)
sigSample.addSystematic(sysObj.AR_all_EESPS_MC)
sigSample.addSystematic(sysObj.AR_all_EESZ_MC)
hasHippo = doToys and "hippo" in userOpts.gridForToys
hasSig = doToys and signame in userOpts.gridForToys
if not doToys or hasHippo or hasSig:
exclusion = configMgr.addFitConfigClone(tlx, "TopLvlXML_Exclusion_%s" % s)
print "!!!--- in LOOP ", s
sigSample = Sample(s, ROOT.kRed)
sigSample.setStatConfig(useStat)
if not useStat:
sigSample.addSystematic(sysObj.AR_mcstat_SIG)
sigUncert = userOpts.sigUncert
sigSample.setFileList([userOpts.sigFile])
if sigUncert == "up":
# sigSample.setWeights(("eventweight","syst_XSUP","1.0"))
sigSample.setWeights(("eventweight", "1.07", "1.0"))
elif sigUncert == "down":
# sigSample.setWeights(("eventweight","syst_XSDOWN","1.0"))
sigSample.setWeights(("eventweight", "0.93", "1.0"))
#!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
# Add systematics here
#!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
# Set normalization factor
sigSample.setNormFactor("mu_SIG", 1.0, 0.0, 10.0)
# Add the sample to the exclusion
exclusion.addSamples(sigSample)
exclusion.setSignalSample(sigSample)
print s, "Leaving LOOP ---!!!"
