* 1187008882.4 (GW170817)
* 933200215 (hardware injection)
* 1132401286.33 (Koi Fish Glitch)
""")
else:
chosen_event = st.sidebar.selectbox('Select Event', eventlist)
t0 = datasets.event_gps(chosen_event)
detectorlist = list(datasets.event_detectors(chosen_event))
detectorlist.sort()
st.subheader(chosen_event)
st.write('GPS:', t0)
# -- Experiment to display masses
try:
jsoninfo = fetch_event_json(chosen_event)
for name, nameinfo in jsoninfo['events'].items():
st.write('Mass 1:', nameinfo['mass_1_source'], 'M$_{\odot}$')
st.write('Mass 2:', nameinfo['mass_2_source'], 'M$_{\odot}$')
#st.write('Distance:', int(nameinfo['luminosity_distance']), 'Mpc')
st.write('Network SNR:', int(nameinfo['network_matched_filter_snr']))
eventurl = 'https://gw-osc.org/eventapi/html/event/{}'.format(chosen_event)
st.markdown('Event page: {}'.format(eventurl))
st.write('\n')
except:
pass
#-- Choose detector as H1, L1, or V1
detector = st.sidebar.selectbox('Detector', detectorlist)
qcenter = st.sidebar.slider('Q-value', 5, 120, 5) # min, max, default
qrange = (int(qcenter * 0.8), int(qcenter * 1.2))
hq = strain.q_transform(outseg=(t0 - dt, t0 + dt), qrange=qrange)
with lock:
fig4 = hq.plot()
ax = fig4.gca()
fig4.colorbar(label="Normalised energy", vmax=vmax, vmin=0)
ax.grid(False)
ax.set_yscale('log')
ax.set_ylim(bottom=15)
st.pyplot(fig4, clear_figure=True)
#-- Try getting pe url
jsoninfo = fetch_event_json(chosen_event, catalog='GWTC-1-confident')
#st.write(jsoninfo)
for name, nameinfo in jsoninfo['events'].items():
for peset, peinfo in nameinfo['parameters'].items():
if 'pe' in peset:
sourceurl = peinfo['data_url']
st.write('PE samples URL: ', sourceurl)
pedata, waveform = load_pe(sourceurl)
st.write('Showing samples for {}'.format(waveform))
#st.write('Got some samples')
#st.write(pedata.dtype.names)
