Hi Jeremy,
The -inf likelihood problem is finally solved. It turns out that calcNeSqdToNormPerKpc3 returns nan if a float32 redshift is loaded into Cosmology, by chance that a FITS catalog stores everything in float32.
A type conversion in Cosmology.__init__() is necessary to avoid this issue.