from lalinference.bayestar import command
import logging
logging.basicConfig(level=logging.DEBUG, format='%(asctime)s %(levelname)s %(message)s')
log = logging.getLogger('BAYESTAR')
methods = '''
toa_phoa_snr
toa_phoa_snr_mcmc
'''.split()
default_method = 'toa_phoa_snr'
command.skymap_parser.add_argument(
'--method', choices=methods, default=default_method,
help='Sky localization method [default: %(default)s]')
parser = command.ArgumentParser(
parents=[command.waveform_parser, command.prior_parser, command.skymap_parser])
parser.add_argument('-N', '--dry-run', default=False, action='store_true',
help='Dry run; do not update GraceDB entry [default: %(default)s]')
parser.add_argument('-o', '--output', metavar='FILE.fits[.gz]',
default='bayestar.fits.gz',
help='Name for uploaded file [default: %(default)s]')
parser.add_argument('graceid', metavar='G123456', nargs='*',
help='Run on this GraceDB ID [default: listen to lvalert]')
opts = parser.parse_args()
#
# Late imports
#
import itertools
#
# You should have received a copy of the GNU General Public License
# along with this program. If not, see <http://www.gnu.org/licenses/>.
#
"""
Create a contours for the credible levels of an all-sky probability map.
The input is a HEALPix probability map.
The output is a GeoJSON FeatureCollection (http://geojson.org/).
"""
# Command line interface
import argparse
from lalinference.bayestar import command
parser = command.ArgumentParser()
parser.add_argument('-o',
'--output',
metavar='FILE.geojson',
default='-',
type=argparse.FileType('w'),
help='name of output file [default: stdout]')
parser.add_argument(
'--contour',
metavar='PERCENT',
type=float,
nargs='+',
required=True,
help='plot contour enclosing this percentage of probability mass')
parser.add_argument(
'-i',
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program. If not, see <http://www.gnu.org/licenses/>.
#
"""
Overplot contours for a large number of sky maps, in geographical coordinates.
This can reveal patterns due to the priors (i.e., the network antenna pattern).
"""
__author__ = "Leo Singer <*****@*****.**>"
# Command line interface
from lalinference.bayestar import command
parser = command.ArgumentParser(parents=[command.figure_parser])
parser.add_argument('--contour',
metavar='PERCENT',
type=float,
default=90,
help='plot contour enclosing this percentage of'
' probability mass [default: %(default)s]')
parser.add_argument(
'--alpha',
metavar='ALPHA',
type=float,
default=0.1,
help='alpha blending for each sky map [default: %(default)s]')
parser.add_argument('fitsfilenames',
metavar='GLOB.fits[.gz]',
nargs='+',
