#include "RooFit.h"
#include "RooCategory.h"
#include "RooRealVar.h"
#include "RooGaussModel.h"
#include "RooArgSet.h"
#include "RooDataSet.h"
#include "RooBinning.h"
#include "RooPlot.h"
#include "RooBDecay.h"
#include "RooMsgService.h"
#include "RooRealConstant.h"
using namespace RooFit;
int bu(  ) {
  RooRealVar t("t","t",-0.1,15,"ps");
  RooCategory tag("tag","tag");
  tag.defineType("unmixed",+1);
  tag.defineType("mixed",-1);

  RooRealVar tau("tau","tau",1.536,1.3,1.7,"ps");
  RooRealVar omega("omega","omega",0.38);
  RooRealVar dmd("dmd","dmd",0.502,0,1);
  RooRealVar DG("DG","DG",0);
  RooRealVar mu("mu","mu",0,-0.1,0.1);
  RooRealVar si("si","si",0.04,0.,0.1);
  
  RooGaussModel gres("gres","gres",t,mu,si);
  RooFormulaVar fcos("fcos","@0*(1-2*@1)",RooArgList(tag,omega));
  RooBDecay sdec("bdec","bdec",t,tau,DG,
                RooRealConstant::value(1),
                RooRealConstant::value(0),
                fcos,
                RooRealConstant::value(0),
                dmd,
                gres,
                RooBDecay::SingleSided);
  
  RooArgSet obs(t,tag);
  
  RooDataSet *data = sdec.generate(obs,10000);
  
  sdec.fitTo(*data);

  RooMsgService::instance().Print() ;
  cout << endl ;

  RooBinning asymBins(0,10);
  int nB=20;
  asymBins.addUniform(nB,0,0);
  for (int i=0;i<nB;++i) 
    {
      asymBins.addBoundary(exp(i*log(1+10)/nB)-1.);
    }
  cout<<"Plotting asymetry"<<endl;
  
  RooPlot *tp = t.frame();
  cout<<"Plotting asymetry of data"<<endl;
  data->plotOn(tp,Asymmetry(tag),Binning(asymBins));
  cout<<"Plotting asymetry of pdf"<<endl;
  RooMsgService::instance().addStream(RooMsgService::DEBUG,Topic(RooMsgService::Plotting)) ;
  sdec.plotOn(tp,Asymmetry(tag),LineColor(kRed)
              ,LineWidth(5));
  
  return 0;
}
//=============================================================================