Remove `log_mean_energy` function ... use `scipy.stats.gmean` instead.
Use the safe tempfile.NamedTemporaryFile in SExtractor unit test.
Add unit test for np_to_pha function.