if fitterPars.includeMuons and fitterPars.includeElectrons:
modeString = ''
elif fitterPars.includeMuons:
modeString = 'Muon'
elif fitterPars.includeElectrons:
modeString = 'Electron'
else:
modeString = ''
pars4 = config.the4BodyConfig(fitterPars, mHbasis, opts.syst,
opts.alpha)
# fitter4 = RooWjjMjjFitter(pars4)
(halfGen, otherdata) = HWWSignalShapes.NgenHiggs(opts.mH, 'HWW')
SignalHist = HWWSignalShapes.makeHiggsHist(opts.mH, pars4, 'HWW')
SignalHist.Scale(otherdata[1]*otherdata[2]*pars4.intLumi* \
opts.factor/otherdata[0])
print "signal only"
SignalHist.Print()
if opts.debug:
print 'drawing signal histogram...'
SignalHist.Draw()
gPad.Update()
gPad.WaitPrimitive()
basisFilename = 'H%i_%s_%iJets_Fit_Shapes.root' % (mHbasis, modeString,
opts.Nj)
print 'opening', basisFilename
bfile = TFile(basisFilename)
bgCurve = bfile.Get('h_total')
fitterPars.smoothingOrder = 0
if fitterPars.includeMuons and fitterPars.includeElectrons:
modeString = ''
elif fitterPars.includeMuons:
modeString = 'Muon'
elif fitterPars.includeElectrons:
modeString = 'Electron'
else:
modeString = ''
pars4 = config.the4BodyConfig(fitterPars, mHbasis, opts.syst, opts.alpha)
# fitter4 = RooWjjMjjFitter(pars4)
(halfGen, otherdata) = HWWSignalShapes.NgenHiggs(opts.mH, 'HWW')
SignalHist = HWWSignalShapes.makeHiggsHist(opts.mH, pars4, 'HWW')
SignalHist.Scale(otherdata[1]*otherdata[2]*pars4.intLumi* \
opts.factor/otherdata[0])
print "signal only"
SignalHist.Print()
if opts.debug:
print 'drawing signal histogram...'
SignalHist.Draw()
gPad.Update()
gPad.WaitPrimitive()
basisFilename = 'H%i_%s_%iJets_Fit_Shapes.root' % (mHbasis, modeString,
opts.Nj)
print 'opening', basisFilename
bfile = TFile(basisFilename)
bgCurve = bfile.Get('h_total')
