It looks like your
makeSwaspRawVisitInfo.py is not significantly modified from the obs_subaru one. In particular, your observatory value is almost certainly wrong (unless your telescope is at the exact same location as Subaru!). I’m assuming you do have all of the other listed fields correctly filled out in your raw files, otherwise singleFrameDriver should fail, but jointcal does need many of those fields correctly filled out for your observatory/observation.
Also, depending on the reference catalog you’re using, jointcal may or may not give a major improvement in astrometric fit at present.
How exactly are you trying to use the output? Are you going through the Butler, or are you trying to read the FITS catalogs directly (not recommended)?
Our catalog positions are Coord objects that contain afw.Angles, which can be converted on the fly via e.g.
coord.getLatitude().asDegrees(). We have plans to eventually make these interoperate with astropy.Quantity, but that’s quite a ways off.
Similarly for fluxes/magnitudes, you want to use the Calib object to convert from calibrated instrument fluxes to magnitudes. Given a calib object from the calexp and a flux (e.g.
catalog['slot_PsfFlux_flux']) from the source catalog, you can get an AB magnitude via
calib.getMagnitude(flux) and a flux in Janskys via
Does this help?