// $Id$
//-----------------------------------------------------------------------
//       The GSI Online Offline Object Oriented (Go4) Project
//         Experiment Data Processing at EE department, GSI
//-----------------------------------------------------------------------
// Copyright (C) 2000- GSI Helmholtzzentrum fuer Schwerionenforschung GmbH
//                     Planckstr. 1, 64291 Darmstadt, Germany
// Contact:            http://go4.gsi.de
//-----------------------------------------------------------------------
// This software can be used under the license agreements as stated
// in Go4License.txt file which is part of the distribution.
//-----------------------------------------------------------------------

/**
 * Go4Fit Example12
 * @brief Fitting of TGraph
 * @file Example12.cxx
 * @ingroup go4_fit
 */

#ifndef __CINT__

#include "TGraph.h"
#include "TApplication.h"

#include "TGo4Fitter.h"
#include "TGo4FitDataGraph.h"

void Example12();

int main(int argc, char **argv)
{
   TApplication theApp("Application", nullptr, nullptr);

   Example12();

   theApp.Run();

   return 0;
}

#endif

// routine to generate example graph function
TGraph *MakeGraph()
{
   TGraph *gr = new TGraph(40);
   for (Int_t i = 0; i < 40; i++) {
      Double_t x = (i + 1.) / 40.;
      Double_t y = 5 - 0.5 * x + x * x;
      if ((x > 0.3) && (x < 0.7))
         y += 10. * (1. - 25. * (x - .5) * (x - .5));
      gr->SetPoint(i, x, y);
   }
   return gr;
}

void Example12()
{
   // create fitter, select fit function and add standard actions list
   TGo4Fitter fitter("Fitter", TGo4Fitter::ff_ML_Poisson, kTRUE);

   // add histogram to fitter, which should be fitted
   fitter.AddGraph("data1", MakeGraph(), kTRUE);

   // add polynom of 2 order for background approx
   fitter.AddPolynomX("data1", "BackgrPol", 2);

   // add polynom of 2 order for background approx
   fitter.AddPolynomX("data1", "LinePol", 2, 123, 0.3, 0.7);

   // perform all actions
   fitter.DoActions();

   // printing of fitter parameters
   fitter.Print("Pars");

   // Draw data ("-" means without full model), Background and line (has group 123)
   fitter.Draw("#data1-, Background, Group123");
}