{
   TH1D* h = new TH1D("h1", "h1", 40, 0.5, 40.5);
   TGraph* g = new TGraph(40);
   
   for(Int_t i = 1; i <= 40; i++) {
      h->Fill(i, 0.001*i);
      g->SetPoint(i-1, i, 0.001*i);
   }

   TCanvas *c = new TCanvas("c1", "c1", 1);
   h->Draw("LP");
   cout << endl << "************ Histogram Fit... ***********" << endl;
   h->Fit("pol0", "", "");
   
   
   TCanvas *c = new TCanvas("c2", "c2", 1);
   g->Draw("ALP");
   cout << endl << "************* Graph Fit.... ************" << endl;
   g->Fit("pol0", "", "");

}