// CPP includes
#include <cstdlib>
#include <cstdio>
#include <iostream>
#include <cmath>
#include <vector>
#include <string>
#include <iostream>
#include <sstream>
#include <fstream>

// C includes
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>

// ROOT includes
#include <TSystem.h>
#include <TROOT.h>
#include <TApplication.h>
#include <TRint.h>
#include <TDirectory.h>
#include <TFile.h>
#include <TChain.h>
#include <TTree.h>
#include <TNtuple.h>
#include <TBranch.h>
#include <TH1.h>
#include <TAxis.h>
#include <TCanvas.h>
#include <TColor.h>
#include <TPad.h>
#include <TStyle.h>
#include <TKey.h>
#include <TLegend.h>
#include <TText.h>
#include <TPaveText.h>
#include <TLine.h>
#include <TMarker.h>
#include <TLatex.h>
#include <TMath.h>
#include <TH1.h>
#include <TH2.h>
#include <TF1.h>
#include <TF2.h>
#include <TString.h>
#include "TObjString.h"
#include "TRandom3.h"
#include "TObjArray.h"
#include "TSpline.h"
#include "TPad.h"
#include "TGraphSmooth.h"
using namespace std;

void run(){
  TH1::AddDirectory(kFALSE);
  TH1::AddDirectory(kFALSE);
 
  
  FILE *datFile; 
  datFile = fopen("test.dat","r");
  double c,d;
  int cont =0;
  while(1)
         {
          fscanf(datFile, "%lf %lf ",&c,&d);
          cont = cont +1;
          if (feof(datFile)!=NULL) break;
          }
          
          cout << cont <<endl;
  double *t= new double[cont];
  double *y= new double[cont];

         
  rewind (datFile);
  for (int i=0;i<cont;i++){
     fscanf(datFile, "%lf %lf ",&t[i],&y[i]); 
    // cout << t[i] << " " << y[i] << endl;
  
                          }

                          
  TGraph*g = new TGraph(cont,y,t);
  g->Draw("AC*");

  TGraph *grout=new TGraph(cont);
  TGraphSmooth *gs1 = new TGraphSmooth("g");
  grout = gs1->Approx(g,"linear",cont);
  TSpline5 *s = new TSpline5("s",grout);//g);
  s->SetLineColor(kRed);
  cout << " ======== " << s->Eval(5e-02)<<endl;
  s->Draw("same");
  fclose(datFile);
  delete t;
  delete y;
   
}