from gwpy.timeseries import StateVector
# define bitmask
bits = [
'Summary state',
'State 1 damped',
'Stage 1 isolated',
'Stage 2 damped',
'Stage 2 isolated',
'Master switch ON',
'Stage 1 WatchDog OK',
'Stage 2 WatchDog OK',
]
# get data
data = StateVector.fetch('L1:ISI-ETMX_ODC_CHANNEL_OUT_DQ', 'May 22 2014 14:00', 'May 22 15:00', bits=bits)
data = data.resample(16)
# make a plot
plot = data.plot(add_label='inset')
plot.set_title('LLO ETMX internal seismic isolation state')
plot.add_bitmask('0b11101110')
if __name__ == '__main__':
try:
outfile = __file__.replace('.py', '.png')
except NameError:
pass
else:
plot.save(outfile)
print("Example output saved as\n%s" % outfile)
