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


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

//generating pseduo data
 RooRealVar a01("a01","a01",-7.17382e-01) ;
 RooRealVar a11("a11","a11",-5.46139e-01) ;
 RooRealVar a21("a21","a21",4.79449e-01) ;
 RooRealVar a31("a31","a31",-8.53512e-02) ; 
 RooRealVar a41("a41","a41",-1.10019e-01) ;
 RooRealVar a51("a51","a51",7.16816e-02) ;
 RooChebychev pdf1("pdf1","pdf1",mass,RooArgList(a01,a11,a21,a31,a41,a51)); 

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

  
 //tmp = (RooDataSet*)tmp->reduce(Cut("80<mass && mass<100 || 140<mass && mass<200"));
// fitting pseduo data with p5 with all parameters floating
 RooRealVar a02("a02","a02",-7.17382e-01,-8.17382e-01,-6.17382e-01) ;
 RooRealVar a12("a12","a12",-5.46139e-01,-6.46139e-01,-4.46139e-01) ;
 RooRealVar a22("a22","a22",4.79449e-01,3.79449e-01,5.79449e-01) ;
 RooRealVar a32("a32","a32",-8.53512e-02,-9.53512e-02,-7.53512e-02) ; 
 RooRealVar a42("a42","a42",-1.10019e-01,-2.10019e-01,-0.90019e-01) ;
 RooRealVar a52("a52","a52",7.16816e-02,6.16816e-02,8.16816e-02) ;
 RooChebychev qcd_model_polynomial("qcd_model_polynomial","qcd_model_polynomial",mass,RooArgList(a02,a12,a22,a32,a42,a52));
 
 //mass.setRange("SB2",140,200);
 mass.setRange("SB1",80,100);
 mass.setRange("SB2",140,200);
 RooFitResult *r = qcd_model_polynomial.fitTo(*tmp, Range("SB1,SB2"),Save());
 /*mass.setRange("SB2",80,100);
 mass.setRange("SB2",140,200);
 
 RooFitResult *r1 = qcd_model_polynomial1.fitTo(*tmp, Range("SB1,SB2"),Save()); 
*/
 RooPlot* xframe = mass.frame(Title("signal 125 GeV")) ;
 tmp->plotOn(xframe,Name("sigbkg"));
 qcd_model_polynomial.plotOn(xframe,Name("model"));
 //mass.setRange("SB1",80,100);
 //mass.setRange("SB2",140,200);
 //RooFitResult *r1 = qcd_model_polynomial.fitTo(*tmp,Save());
 //qcd_model_polynomial.plotOn(xframe,Name("model1"), LineColor(kGreen), LineStyle(kDotted));

 TCanvas *c1 = new TCanvas();
 c1->cd();xframe->Draw();


}