def theConfig(Nj, mcdir = '', initFile = '', mH=400, mvaCut = None,
qgCut = None):
optPars = VBFHWWConfig.getOptimalPars(mH, optimalPars2)
mvaVarNames = optPars[1]
mvaCutValues = optPars[2]
VBFHWWConfig.minMlvjj = optPars[3]
VBFHWWConfig.maxMlvjj = optPars[4]
fitterPars = VBFHWWConfig.theConfig(Nj, mcdir, initFile)
fitterPars.includeMuons = True
fitterPars.includeElectrons = False
if mvaCut:
mvaCutValues = mvaCut
fitterPars.cuts += '&& (%s > %f) ' % (mvaVarNames,mvaCutValues)
# if not qgCut:
# qgCut = optPars[11]
# fitterPars.cuts += '&& (%s > %f) ' % ("qgld_Summer11CHS[1]",qgCut)
print '2-body cuts:', fitterPars.cuts
return fitterPars
def theConfig(Nj, mcdir = '', initFile = '', mH=400, mvaCut = None,
qgCut = None):
optPars = VBFHWWConfig.getOptimalPars(mH, optimalPars2)
mvaVarNames = optPars[1]
mvaCutValues = optPars[2]
VBFHWWConfig.minMlvjj = optPars[3]
VBFHWWConfig.maxMlvjj = optPars[4]
fitterPars = VBFHWWConfig.theConfig(Nj, mcdir, initFile)
fitterPars.includeMuons = False
fitterPars.includeElectrons = True
if mvaCut:
mvaCutValues = mvaCut
fitterPars.cuts += '&& (%s > %f) ' % (mvaVarNames,mvaCutValues)
# if not qgCut:
# qgCut = optPars[11]
# fitterPars.cuts += '&& (%s > %f) ' % ("qgld_Summer11CHS[1]",qgCut)
print '2-body cuts:', fitterPars.cuts
return fitterPars
def the4BodyConfig(twoBodyConfig, mH=400, syst=0, alphaOverride = None):
optPars = VBFHWWConfig.getOptimalPars(mH, optimalPars2)
VBFHWWConfig.minMlvjj = optPars[3]
VBFHWWConfig.maxMlvjj = optPars[4]
alpha = optPars[6+syst]
alphaDown = optPars[6+1]
alphaUp = optPars[6+2]
if alphaOverride:
alpha = alphaOverride
alphaDown = alphaOverride
alphaUp = alphaOverride
fitterPars = VBFHWWConfig.the4BodyConfig(twoBodyConfig, alpha, optPars[9],
optPars[10],
alphaRange = (alphaDown,alphaUp))
fitterPars.nbins = optPars[5]
fitterPars.model = 1
if (mH <= 200) and (twoBodyConfig.njets == 2):
fitterPars.model = 4
if (mH == 190) and (twoBodyConfig.njets == 2):
fitterPars.model = 1
return fitterPars
def the4BodyConfig(twoBodyConfig, mH=400, syst=0, alphaOverride = None):
optPars = VBFHWWConfig.getOptimalPars(mH, optimalPars2)
VBFHWWConfig.minMlvjj = optPars[3]
VBFHWWConfig.maxMlvjj = optPars[4]
alpha = optPars[6+syst]
alphaDown = optPars[6+1]
alphaUp = optPars[6+2]
if alphaOverride:
alpha = alphaOverride
alphaDown = alphaOverride
alphaUp = alphaOverride
fitterPars = VBFHWWConfig.the4BodyConfig(twoBodyConfig, alpha, optPars[9],
optPars[10],
alphaRange = (alphaDown,alphaUp))
fitterPars.nbins = optPars[5]
fitterPars.model = 1
if (mH <= 200) and (twoBodyConfig.njets == 2):
fitterPars.model = 4
if (mH == 190) and (twoBodyConfig.njets == 2):
fitterPars.model = 1
return fitterPars
