#include <iostream>
#include "TH1F.h"

using namespace std;
void test2() {


  TH1F isto("pippo","il mio primo istogramma", 100, 0., 10.);
  int num = 10000;
  for (int i = 0; i < num ; ++i) {
    //isto.Fill(gRandom->Rndm());
    isto.Fill(gRandom->Gaus(5.,2.));
  }

  isto.Fit("gaus");

  TF1* gaus = (TF1*)( isto.GetFunction("gaus"));
  float mean = gaus->GetParameter(1);
  float meanErr = gaus->GetParError(1);
  cout << "mean " << mean << "+/-" << meanErr<< endl;

  // isto.Fill(4.5);
  // isto.Fill(1.1, 0.2);

  gStyle->SetOptFit(111111);
  isto.SetMinimum(0.); // ok
  isto.SetLineColor(4); // ok
  isto.DrawCopy("e"); // ok
  float mean=isto.GetMean();
  // TH1F* pisto = new TH1F("pippo","il mio primo istogramma", 100, 0., 10.);
  // pisto->Draw(); // ok
  //pippo.Draw(); // solo in root

  cout << "Media " << mean  << endl;


  TH2F isto2d("pippo2d","il mio primo istogramma 2d", 100, 0., 10., 100, 0., 10.);
  int numx = 10000;
  for (int i = 0; i < num ; ++i) {
    //isto.Fill(gRandom->Rndm());
    isto2d.Fill(gRandom->Gaus(5.,2.),gRandom->Landau(5.,.2));
  }
  isto2d.DrawCopy("box");

}