Commit b7ff1119 authored by Meg Millhouse's avatar Meg Millhouse

Update megasky to read in chains for fullOnly runs

parent ac233364
......@@ -208,8 +208,12 @@ def make_skyview(directory='.', mdc=None, NSIDE=128, inj=None, npost=5000, geo=F
# -- Input skymap data
print("Extracting RA/DEC samples")
filename = './chains/' + 'signal_params_h0.dat.0'
data = np.loadtxt(filename, unpack=True,usecols=(0,1,2))
try:
filename = './chains/' + 'signal_params_h0.dat.0'
data = np.loadtxt(filename, unpack=True,usecols=(0,1,2))
except:
filename = './chains/' + 'full_params_h0.dat.0'
data = np.loadtxt(filename, unpack=True,usecols=(0,1,2))
ralist = data[1]
sin_dec = data[2]
print("Total samples are {0}".format(ralist.size))
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment