import h5py
import sys
import numpy as np
import lal
filename = sys.argv[1]
f = h5py.File(filename,'r')
pos = f['samples']
ii =  np.argmax(pos['loglikelihood'])
print('maxL values for tc, ra and dec are:')
print(pos['delta_tc'][ii]+f.attrs['trigger_time'])
print(pos['ra'][ii])
print(pos['dec'][ii])