#include <iostream>
#include "TCanvas.h"
#include "TH1.h"
#include "TF1.h"
#include "TStyle.h"
#include "TFile.h"

using namespace std;
double DF_dy_3a_fix_new(TFile*, double, double , double );

void vahe()
{
  TFile *f = new TFile("hist3a.root");
  double DF = DF_dy_3a_fix_new(f,3.22,-60.0,60.0);

}

double DF_dy_3a_fix_new(TFile *f, double center, double fitmin, double fitmax)
{
  gStyle->SetOptStat(000000000);
  gStyle->SetOptFit(1);
  double fit_min =  center - fitmin;
  double fit_max =  center + fitmax;

  TCanvas *c1 = new TCanvas("c1","",200,10,700,500);
  c1 -> cd();

  TH1D *pDeltYDF1    = (TH1D*)f -> Get("pDeltYDF1");
  TH1D *pDeltYPi0DF  = (TH1D*)f -> Get("pDeltYPi0DF");

  TF1 *f1 = new TF1("f1","[0]*(1+[1]*(x-[5])+[2]*(x-[5])**2+[3]*(x-[5])**3+[4]*(x-[5])**4)",fit_min,fit_max);

  f1 -> FixParameter(1,0.0);
  f1 -> FixParameter(3,0.0);

  pDeltYPi0DF -> Fit("f1","","",fitmin,fitmax);
//  pDeltYPi0DF -> Fit("f1","","",-60,60);

  c1 -> Print("kin3a_dy.ps");
  return 0;
}