Python RooWjjMjjFitter.resetYields示例

示例#1

文件： runHWWFitter.py 项目： cms-analysis/ElectroWeakAnalysis-VPlusJets

#     theFitter.getWorkSpace().var('fSU').setConstant(True)

# theFitter.getWorkSpace().var('nDiboson').setConstant()
#theFitter.getWorkSpace().var('nTTbar').setConstant()
# theFitter.getWorkSpace().var('nSingleTop').setConstant()
# theFitter.getWorkSpace().var('nZjets').setConstant()
#theFitter.getWorkSpace().var('nWjets').setConstant()
theFitter.getWorkSpace().var('nQCD').setConstant()

if (opts.Nj == 3):
    theFitter.getWorkSpace().var('nTTbar').setConstant()
if (opts.mH <= 300) and (opts.Nj == 2):
    theFitter.getWorkSpace().var('nTTbar').setConstant()

theFitter.loadData()
theFitter.resetYields()
if fitterPars.includeMuons:
    theFitter.getWorkSpace().var('nQCD').setConstant()

fr = theFitter.fit()

tries = 1
ndf = Long(fr.floatParsFinal().getSize()-5)
while (fr.covQual() < 2) and (tries < 3):
    print "Fit didn't converge well.  Will try again."
    # theFitter.getWorkSpace().var('nDiboson').setConstant()
    theFitter.getWorkSpace().var('nTTbar').setConstant()
    # theFitter.getWorkSpace().var('nSingleTop').setConstant()
    theFitter.getWorkSpace().var('nQCD').setConstant()
    # theFitter.getWorkSpace().var('nZjets').setConstant()
    #theFitter.getWorkSpace().var('nWjets').setConstant()

示例#2

文件： runHWWFitter.py 项目： VPlusJetsAnalyzers/VPlusJets

#     theFitter.getWorkSpace().var('fSU').setConstant(True)

# theFitter.getWorkSpace().var('nDiboson').setConstant()
#theFitter.getWorkSpace().var('nTTbar').setConstant()
# theFitter.getWorkSpace().var('nSingleTop').setConstant()
# theFitter.getWorkSpace().var('nZjets').setConstant()
#theFitter.getWorkSpace().var('nWjets').setConstant()
theFitter.getWorkSpace().var('nQCD').setConstant()

if (opts.Nj == 3):
    theFitter.getWorkSpace().var('nTTbar').setConstant()
if (opts.mH <= 300) and (opts.Nj == 2):
    theFitter.getWorkSpace().var('nTTbar').setConstant()

theFitter.loadData()
theFitter.resetYields()
if fitterPars.includeMuons:
    theFitter.getWorkSpace().var('nQCD').setConstant()

fr = theFitter.fit()

tries = 1
ndf = Long(fr.floatParsFinal().getSize()-5)
while (fr.covQual() < 2) and (tries < 3):
    print "Fit didn't converge well.  Will try again."
    # theFitter.getWorkSpace().var('nDiboson').setConstant()
    theFitter.getWorkSpace().var('nTTbar').setConstant()
    # theFitter.getWorkSpace().var('nSingleTop').setConstant()
    theFitter.getWorkSpace().var('nQCD').setConstant()
    # theFitter.getWorkSpace().var('nZjets').setConstant()
    #theFitter.getWorkSpace().var('nWjets').setConstant()