TString dir = gROOT->GetTutorialDir();

TFITSHDU hdu(dir + "/fitsio/sample2.fits[1][DATAMAX > 2e-15]");

hdu.Print("T");

hdu.Print("T+");

std::unique_ptr<TVectorD> vp(hdu.GetTabRealVectorColumn("DATAMAX"));
const auto &v = *vp;
std::cout << "v[0] = " << v[0] << std::endl;
std::cout << "Does the matched row have DATAMAX > 2e-15? :-)" << std::endl;

gROOT->GetListOfCanvases()->Draw()