coinc_file, psd_file, opts.waveform, opts.f_low,
opts.min_distance, opts.max_distance, opts.prior_distance_power,
nside=opts.nside, f_high_truncate=opts.f_high_truncate,
method=opts.method, chain_dump=chain_dump,
enable_snr_series=opts.enable_snr_series))
sky_map.meta['objid'] = str(graceid)
sky_map.meta['url'] = 'https://gracedb.ligo.org/events/{0}'.format(graceid)
log.info("sky localization complete")
# upload FITS file
with command.TemporaryDirectory() as fitsdir:
fitspath = os.path.join(fitsdir, opts.output)
fits.write_sky_map(fitspath, sky_map, nest=True)
log.debug('wrote FITS file: %s', opts.output)
if opts.dry_run:
command.rename(fitspath, os.path.join('.', opts.output))
else:
gracedb.writeLog(
graceid, "BAYESTAR rapid sky localization ready",
filename=fitspath, tagname=("sky_loc", "lvem"))
log.debug('uploaded FITS file')
except:
# Produce log message for any otherwise uncaught exception
log.exception("sky localization failed")
if opts.dry_run:
# Then re-raise the exception if we are in dry-run mode
raise
if not opts.dry_run:
# Remove old log handler
logging.root.removeHandler(handler)
phase_convention=opts.phase_convention, nside=opts.nside,
f_high_truncate=opts.f_high_truncate,
method=opts.method, chain_dump=chain_dump)
prob, distmu, distsigma, _ = sky_map
distmean, diststd = distance.parameters_to_marginal_moments(
prob, distmu, distsigma)
log.info("sky localization complete")
# upload FITS file
fitsdir = tempfile.mkdtemp()
try:
fitspath = os.path.join(fitsdir, opts.output)
fits.write_sky_map(fitspath, sky_map, gps_time=float(epoch),
creator=parser.prog, objid=str(graceid),
url='https://gracedb.ligo.org/events/{0}'.format(graceid),
runtime=elapsed_time, instruments=instruments,
distmean=distmean, diststd=diststd,
origin='LIGO/Virgo', nest=True)
if not opts.dry_run:
gracedb.writeLog(graceid, "INFO:BAYESTAR:uploaded sky map",
filename=fitspath, tagname=("sky_loc", "lvem"))
else:
command.rename(fitspath, os.path.join('.', opts.output))
finally:
shutil.rmtree(fitsdir)
except:
# Produce log message for any otherwise uncaught exception
log.exception("sky localization failed")
# Then re-raise the exception
raise