## 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 if (opts.mH >= 500): iwt = 1 sigHists = HWWSignalShapes.GenHiggsHists(pars4, opts.mH, fitUtils, iwt = iwt) extraFactor = 2. otherdata = HWWSignalShapes.NgenHiggs(opts.mH, 'HWW') SigVisual = TH1D(sigHists['HWW']) SigVisual.Print() SigVisual.SetName('SigVisual') SigVisual.SetLineColor(kBlue) SigVisual.SetLineWidth(3) SigVisual.Scale(pars4.intLumi*extraFactor*otherdata[1]*otherdata[2]/2.) SigVisual.Print() SigVisualLog = TH1D(SigVisual) SigVisualLog.SetName("SigVisualLog") c4body = TCanvas('c4body', '4 body stacked') mf4.addTH1(SigVisual, "hist")