void ToyMCTest() { TFile *fTemplatesBkg = TFile::Open("templates_Bkg_2016.root"); RooWorkspace *wTemplatesBkg = (RooWorkspace *)fTemplatesBkg->Get("w"); RooRealVar mTop("mTop", "mTop", 50., 300.); RooAbsPdf *pdf_qcd_2b = (RooAbsPdf *)wTemplatesBkg->pdf("qcd_pdf"); RooRealVar *kQCD2b = new RooRealVar("kQCD_2b", "kQCD_2b", 1e-3, -1, 1); kQCD2b->setConstant(false); RooRealVar *nFitQCD2b = new RooRealVar("nFitQCD_2b", "nFitQCD_2b", 10000, 0, 10e+4); RooFormulaVar qcdCor_2b("qcdCor_2b", "(1+@0*@1)*@2", RooArgList(mTop, *kQCD2b, *nFitQCD2b)); //RooEffProd pdf_qcdCor_2b("qcdCor_pdf_2b", "qcdCor_pdf_2b", *pdf_qcd_2b, qcdCor_2b); RooAddPdf addPdf("model_2b", "model_2b", RooArgList(*pdf_qcd_2b), RooArgList(qcdCor_2b)); RooDataSet *data_gen = addPdf.generate(mTop, 10000); }