TString dir = gROOT->GetTutorialDir();

TFITSHDU hdu(dir + "/fitsio/sample5.fits[1]");

std::unique_ptr<TVectorD> vec1(hdu.GetTabRealVectorColumn("ra"));
std::unique_ptr<TVectorD> vec2(hdu.GetTabRealVectorColumn("dec"));
std::unique_ptr<TVectorD> vec3(hdu.GetTabRealVectorColumn("flux_g"));
std::unique_ptr<TVectorD> vec4(hdu.GetTabRealVectorColumn("brick_primary"));

for (auto i : ROOT::TSeqI(vec1->GetLwb(), vec1->GetUpb())) {
   const auto bp = (*vec4)[i];
   if (bp) {
      const auto gflux = (*vec3)[i];
      const auto ra = (*vec1)[i];
      const auto dec = (*vec2)[i];
      printf("RA %f DEC %f G-FLUX %f\n", ra, dec, gflux);
   }
}

gROOT->GetListOfCanvases()->Draw()