%%cpp -d
#include "RooRealVar.h"
#include "RooDataSet.h"
#include "RooGaussian.h"
#include "RooPolynomial.h"
#include "RooAddPdf.h"
#include "TCanvas.h"
#include "TAxis.h"
#include "RooPlot.h"
using namespace RooFit;

RooRealVar x("x", "x", -10, 10);
RooRealVar m("m", "m", 0, -10, 10);
RooRealVar s("s", "s", 1, -10, 10);
RooGaussian gauss("g", "g", x, m, s);

RooRealVar p0("p0", "p0", 0.01, 0., 1.);
RooPolynomial poly("p", "p", x, p0);

RooRealVar f("f", "f", 0.5, 0., 1.);
RooAddPdf model("model", "model", RooArgSet(gauss, poly), f);

std::unique_ptr<RooDataSet> d{model.generate(x, 1000)};
model.fitTo(*d, PrintLevel(-1));

std::unique_ptr<RooArgSet> params{model.getParameters(x)};

params->writeToFile("rf505_asciicfg_example.txt");

TString dir1 = gROOT->GetTutorialDir() ;
dir1.Append("/roofit/rf505_asciicfg.txt") ;
TString dir2 = "rf505_asciicfg_example.txt";

params->readFromFile(dir2);
params->Print("v");

params->readFromFile(dir1, 0, "Section2");
params->Print("v");

params->readFromFile(dir1, "READ", "Section3");

cout << "The following parameters of the were _not_ read from Section3: "
     << (*params->selectByAttrib("READ", false)) << endl;

params->readFromFile(dir1, 0, "Section4");
params->Print("v");