Skip to content
Snippets Groups Projects

update pipe for olib triggers

Merged Meg Millhouse requested to merge meg.millhouse/bayeswave:get-olib-trigs into master
1 file
+ 63
25
Compare changes
  • Side-by-side
  • Inline
@@ -352,36 +352,74 @@ class eventTrigger:
gracedb = GraceDb()
event = gracedb.event(graceid)
event_info = event.json()
# First check for olib info
if event_info['submitter'] == 'oLIB':
try:
self.rho = event_info['extra_attributes']['LalInferenceBurst']['bci']
# put bci in the place of rho for bookkeping
except KeyError:
print >> sys.stderr, \
"graceDB UID %s has no MultiBurst snr attribute"%(graceid)
# Get loudness (for informational, not analysis, purposes)
try:
self.rho = event_info['extra_attributes']['MultiBurst']['snr']
except KeyError:
print("graceDB UID %s has no MultiBurst snr attribute"%(graceid), file=sys.stderr)
# Set time
self.trigger_time = event_info['gpstime']
# Set time
self.trigger_time = event_info['gpstime']
# Set frequency
try:
self.trigger_frequency = \
event_info['extra_attributes']['LalInferenceBurst']['frequency_median']
if self.trigger_frequency < self.frequency_threshold:
self.srate = self.min_srate
self.seglen = self.max_seglen
self.window = self.max_window
self.flow = self.min_flow
else:
self.srate = self.max_srate
self.seglen = self.min_seglen
self.window = self.min_window
self.flow = self.max_flow
except KeyError:
print >> sys.stderr, \
"graceDB UID %s has no MultiBurst central_freq attribute"%(graceid)
print >> sys.stderr, "...using default sample rate"
self.srate = self.default_srate
# Set frequency
try:
self.trigger_frequency = \
event_info['extra_attributes']['MultiBurst']['central_freq']
else:
# Assume it's cwb trigger
# Get loudness (for informational, not analysis, purposes)
try:
self.rho = event_info['extra_attributes']['MultiBurst']['snr']
except KeyError:
print >> sys.stderr, \
"graceDB UID %s has no MultiBurst snr attribute"%(graceid)
if self.trigger_frequency < self.frequency_threshold:
self.srate = self.min_srate
self.seglen = self.max_seglen
self.window = self.max_window
self.flow = self.min_flow
else:
self.srate = self.max_srate
self.seglen = self.min_seglen
self.window = self.min_window
self.flow = self.max_flow
# Set time
self.trigger_time = event_info['gpstime']
except KeyError:
print("graceDB UID %s has no MultiBurst central_freq attribute"%(graceid), file=sys.stderr)
print("...using default sample rate", file=sys.stderr)
self.srate = self.default_srate
# Set frequency
try:
self.trigger_frequency = \
event_info['extra_attributes']['MultiBurst']['central_freq']
if self.trigger_frequency < self.frequency_threshold:
self.srate = self.min_srate
self.seglen = self.max_seglen
self.window = self.max_window
self.flow = self.min_flow
else:
self.srate = self.max_srate
self.seglen = self.min_seglen
self.window = self.min_window
self.flow = self.max_flow
except KeyError:
print >> sys.stderr, \
"graceDB UID %s has no MultiBurst central_freq attribute"%(graceid)
print >> sys.stderr, "...using default sample rate"
self.srate = self.default_srate
Loading