if not (yields[v1].GetName() == 'nDiboson') \
and not (yields[v2].GetName() == 'nDiboson'):
sig2 += covMatrix(v1, v2)
sig2All += covMatrix(v1, v2)
from plotMjjFit import plot2BodyDist
if (opts.Err >= 0) and (opts.Err < 1):
opts.Err = sqrt(sig2)
leftLegend = True #(not opts.btag)
PrefixString = "Diboson"
if opts.btag:
PrefixString += "_btag"
(plots, cans) = plot2BodyDist(theFitter, fitterPars, chi2, ndf,
opts.Err, Prefix = PrefixString,
Left = leftLegend)
h_total = plots[0].getCurve('h_total')
theData = plots[0].getHist('theData')
## mass.setRange('signal', fitterPars.minTrunc, fitterPars.maxTrunc)
## sigInt = theFitter.makeFitter().createIntegral(iset, 'signal')
## sigFullInt = theFitter.makeFitter().createIntegral(iset)
## dibosonInt = theFitter.makeDibosonPdf().createIntegral(iset, 'signal')
## dibosonFullInt = theFitter.makeDibosonPdf().createIntegral(iset)
## WpJInt = theFitter.makeWpJPdf().createIntegral(iset, 'signal')
## WpJFullInt = theFitter.makeWpJPdf().createIntegral(iset)
## ttbarInt = theFitter.makettbarPdf().createIntegral(iset, 'signal')
## ttbarFullInt = theFitter.makettbarPdf().createIntegral(iset)
## SingleTopInt = theFitter.makeSingleTopPdf().createIntegral(iset, 'signal')
covMatrix = TMatrixDSym(fr.covarianceMatrix())
sig2 = 0.
for v1 in range(0, covMatrix.GetNrows()):
for v2 in range(0, covMatrix.GetNcols()):
if ((yields[v1].GetName())[0] == 'n') and \
((yields[v2].GetName())[0] == 'n'):
#print v1, yields[v1].GetName(),',', v2, yields[v2].GetName()
sig2 += covMatrix(v1, v2)
from plotMjjFit import plot2BodyDist
if (opts.Err >= 0) and (opts.Err < 1):
opts.Err = sqrt(sig2)*1.37
(plots, cans) = plot2BodyDist(theFitter, fitterPars, chi2, ndf,
opts.Err, opts.NP)
h_total = plots[0].getCurve('h_total')
theData = plots[0].getHist('theData')
mass.setRange('signal', fitterPars.minTrunc, fitterPars.maxTrunc)
sigInt = theFitter.makeFitter().createIntegral(iset,iset,'signal')
sigFullInt = theFitter.makeFitter().createIntegral(iset,iset)
dibosonInt = theFitter.makeDibosonPdf().createIntegral(iset,iset,'signal')
dibosonFullInt = theFitter.makeDibosonPdf().createIntegral(iset,iset)
WpJInt = theFitter.makeWpJPdf().createIntegral(iset,iset,'signal')
WpJFullInt = theFitter.makeWpJPdf().createIntegral(iset,iset)
ttbarInt = theFitter.makettbarPdf().createIntegral(iset,iset,'signal')
ttbarFullInt = theFitter.makettbarPdf().createIntegral(iset,iset)
SingleTopInt = theFitter.makeSingleTopPdf().createIntegral(iset,iset,'signal')
SingleTopFullInt = theFitter.makeSingleTopPdf().createIntegral(iset,iset)
