TString dir = gROOT->GetTutorialDir();

TFITSHDU hdu(dir + "/fitsio/sample1.fits");

hdu.Print("F+");

printf("Exposure time = %s\n", hdu.GetKeywordValue("EXPTIME").Data());

std::unique_ptr<TMatrixD> mat(hdu.ReadAsMatrix(0));
mat->Print();

TImage * im = (TImage *)hdu.ReadAsImage(0);

TH1 * hist = (TH1 *)hdu.ReadAsHistogram();

auto c = new TCanvas("c1", "FITS tutorial #1", 1400, 800);
c->Divide(2, 1);
c->cd(1);
im->Draw();
c->cd(2);
hist->Draw("COL");

gROOT->GetListOfCanvases()->Draw()