using namespace RooFit ;
//void fit_qcdANDsignal_toyMC(int niter, int nqcd, int ntop, int nzqq, int nsig)


void toyMC_study_v1()
{
  
 RooRealVar mass("mass","mass",80,200) ;


 RooRealVar meanV_top("meanV_top","meanV",115);
 RooRealVar widthV_top("widthV_sig","widthV",39);
 RooRealVar sigmaV_top("sigmaV_sig","sigmaV",15);
 RooVoigtian pdf2("pdf2","pdf2",mass,meanV_top,widthV_top,sigmaV_top);

 RooDataSet *tmp = pdf2.generate(mass,50000);

 
 RooRealVar meanV_top1("meanV_top1","meanV",115,100,130);
 RooRealVar widthV_top1("widthV_top1","widthV",39,20,50);
 RooRealVar sigmaV_top1("sigmaV_top1","sigmaV",15,5,30);
 RooVoigtian pdf3("pdf3","pdf3",mass,meanV_top1,widthV_top1,sigmaV_top1);

 
 
 mass.setRange("SB1",80,110);
 mass.setRange("SB2",150,200);
 RooFitResult *r = pdf3.fitTo(*tmp,Range("SB1,SB2"),Save());

 
 RooPlot* xframe = mass.frame(Title("signal 125 GeV")) ;
 tmp->plotOn(xframe,Name("sigbkg"));
 pdf3.plotOn(xframe,Name("model"));
 pdf3.plotOn(xframe,Name("model_fullRange"),Range(80,200),LineStyle(kDotted));
 
 TCanvas *c1 = new TCanvas();
 c1->cd();xframe->Draw();


}