import h5py
import sys
import numpy as np
import lal
filename = sys.argv[1]
f = h5py.File(filename,'r')
pos = f['C01:IMRPhenomXPHM']['posterior_samples'][:]
ii =  np.argmax(pos['log_likelihood'])
print('maxL values for tc, ra and dec are:')
print(pos['geocent_time'][ii])
print(pos['ra'][ii])
print(pos['dec'][ii])
print('Final mass is:')
print(pos['final_mass'][ii])
print('10M time is:', 10* pos['final_mass'][ii]*lal.MTSUN_SI)