#include <TROOT.h>
#include <TFile.h>
#include <TCanvas.h>

#include <RooFitResult.h>
#include <RooDataSet.h>
#include <RooPlot.h>
#include <RooArgSet.h>
#include <RooAbsPdf.h>
#include <RooRealVar.h>
#include <RooWorkspace.h>
#include <RooNLLVar.h>

#include "../RooBetaFun.h"
#include "../RooGaussian3.h"
#include "../RooRadiativeTail.h"

int main(int argc, char *argv[]) {

  gROOT->SetStyle("Plain"); 
  bool debug = false;

  TFile f("Simplot_FitRes.root");
  f.ls();

  TCanvas *c1 = new TCanvas("c1","c1",600,600);
  c1->cd();
  // Retrieve workspace from file
  RooWorkspace* w; 
  w =  ((RooWorkspace*)gDirectory->Get("w"));
  //  f.GetObject("w",w) ;
  cout<<" 1 "<<endl;
  // Retrieve x,model and data from workspace
  RooRealVar* mpipi = w->var("mpipi") ;
  RooRealVar* pA = w->var("pA") ;
  RooRealVar* pidp = w->var("pidp") ;
  RooRealVar* pidm = w->var("pidm") ;
  RooRealVar* pidKp = w->var("pidKp") ;
  RooRealVar* pidKm = w->var("pidKm") ;
  RooAbsPdf* pdf = w->pdf("pdf") ;
  RooAbsData* pl_data = w->data("pdfData") ;
  RooPlot *myf, *myf_pA;
  if(mpipi) myf = mpipi->frame();
  if(pA) myf_pA = pA->frame();
  cout<<"umm"<<endl;

  pdf->printCompactTree();

  if(myf) {
    pl_data->plotOn(myf);
    if(pdf) {
      pdf->printCompactTree();
      //      pdf->plotOn(myf);
    }
  }
  if(myf_pA) {
    pl_data->plotOn(myf_pA);
    if(pdf) {
      pdf->printCompactTree();
      pdf->plotOn(myf_pA);
    }
  }

  cout<<" 2 "<<endl;
  //  RooAbsReal * lik = w->function("nll_fit");   
  //  RooFitResult * result = ((RooFitResult*)gDirectory->Get("L"));   
  myf->Draw();
  //  if(debug) pdf->printCompactTree();   
  c1->Print("myplot.eps");
  return 0;
}