test_axis_range() {
  TH1F * h = new TH1F("h", "h", 50, -3, 3);
  h->FillRandom("gaus", 10000);

  cout << "==Before changing the Xaxis range== GetFirst(): "
       << h->GetXaxis()->GetFirst() << endl;
  cout << "==Before changing the Yaxis range== GetFirst(): "
       << h->GetYaxis()->GetFirst() << endl;

  h->GetXaxis()->SetRange(4, 48);
  h->GetYaxis()->SetRangeUser(100, 200);
  
  cout << "==After changing the Xaxis range== GetFirst(): "
       << h->GetXaxis()->GetFirst() << endl;
  cout << "==After changing the Yaxis range== GetFirst(): "
       << h->GetYaxis()->GetFirst() << endl;

  h->Draw();
}