nll=fr.minNll() print '***** nll = ',nll,' ***** \n' print 'total yield: {0:0.0f} +/- {1:0.0f}'.format(totalYield, sqrt(sig2)) # assert(False) cWpJ = TCanvas('cWpJ', 'W+jets shape') pars4 = config.the4BodyConfig(fitterPars, opts.mH, opts.syst, opts.alpha) pars4.mHiggs = opts.mH pars4.wHiggs = HWWSignalShapes.HiggsWidth[opts.mH] pars4.initParamsFile = sigYieldFilename fitter4 = RooWjjMjjFitter(pars4) fitter4.makeFitter((opts.ParamWpJ>=0)) fitter4.loadData() fitter4.make4BodyPdf(theFitter) fitter4.loadParameters(sigYieldFilename) ## assert(False) mf4 = fitter4.stackedPlot(False, RooWjjMjjFitter.mlnujj) mf4.SetName("Mlvjj_Stacked") ## sf4 = fitter4.residualPlot(mf4, "h_background", "dibosonPdf", False) # pf4 = fitter4.residualPlot(mf4, "h_total", "", True) pf4 = pulls.createPull(mf4.getHist('theData'), mf4.getCurve('h_total')) pf4.SetName("Mlvjj_Pull") lf4 = fitter4.stackedPlot(True, RooWjjMjjFitter.mlnujj) lf4.SetName("Mlvjj_log") fitUtils = RooWjjFitterUtils(pars4) iwt = 0
nll=fr.minNll() print '***** nll = ',nll,' ***** \n' print 'total yield: {0:0.0f} +/- {1:0.0f}'.format(totalYield, sqrt(sig2)) # assert(False) cWpJ = TCanvas('cWpJ', 'W+jets shape') pars4 = config.the4BodyConfig(fitterPars, opts.mH, opts.syst, opts.alpha) pars4.mHiggs = opts.mH pars4.wHiggs = HWWSignalShapes.HiggsWidth[opts.mH] pars4.initParamsFile = sigYieldFilename fitter4 = RooWjjMjjFitter(pars4) fitter4.makeFitter((opts.ParamWpJ>=0)) fitter4.loadData() fitter4.make4BodyPdf(theFitter) fitter4.loadParameters(sigYieldFilename) ## assert(False) mf4 = fitter4.stackedPlot(False, RooWjjMjjFitter.mlnujj) mf4.SetName("Mlvjj_Stacked") ## sf4 = fitter4.residualPlot(mf4, "h_background", "dibosonPdf", False) # pf4 = fitter4.residualPlot(mf4, "h_total", "", True) pf4 = pulls.createPull(mf4.getHist('theData'), mf4.getCurve('h_total')) pf4.SetName("Mlvjj_Pull") rf4 = pulls.createResid(mf4.getHist('theData'), mf4.getCurve('h_total')) rf4.SetName("Mlvjj_Residuals") (chi2_4, ndf_4) = pulls.computeChi2(mf4.getHist('theData'), mf4.getCurve('h_total')) lf4 = fitter4.stackedPlot(True, RooWjjMjjFitter.mlnujj)