#include<stdio.h>
#include<stdlib.h>
#include<math.h>

#include "TH1F.h"
#include "TF1.h"

#define MAX 0.1
#define MIN -0.1

void fit(){

  ifstream input("dati1.dat",ios::in);
  TH1F *histo= new TH1F("histo","First_histo",1000,0.,MAX);

  while(!input.eof()){
    double variable;
    input>>variable;
    histo->Fill(variable);
  }


  TF1 *Fit_Gauss= new TF1("Fit_Gauss","[0]*exp(-0.5*((x-[1])/[2])^2)/([2]*sqrt(2*M_PI))",0.,MAX);


  Fit_Gauss->SetParameters(1.,0.,1.);
  histo->Fit("Gauss");
  histo->Draw("AP*");

}


//TF1 *Fit_Landau=new TF1("Fit_Landau","[0]*TMath::Landau(x,[1],[2])",MIN,MAX);

//TF1 *Fit_Poisson=new TF1("Fit_Poisson","[0]*x",MIN,MAX);