Esempio n. 1
0
        Option("--max-distance", type=float, metavar="Mpc",
            help="Maximum distance of prior in megaparsecs [default: infer from effective distance]"),
        Option("--prior-distance-power", type=int, metavar="-1|2",
            help="Distance prior [-1 for uniform in log, 2 for uniform in volume, default: 2]"),
        Option("--method", choices=methods, metavar='|'.join(methods), default=[], action="append",
            help="Sky localization method [may be specified multiple times, default: "
            + default_method + "]"),
        Option("--chain-dump", default=False, action="store_true",
            help="For MCMC methods, dump the sample chain to disk [default: no]"),
        Option("--keep-going", "-k", default=False, action="store_true",
            help="Keep processing events if a sky map fails to converge [default: no].")
    ]
)
opts, args = parser.parse_args()
infilename = command.get_input_filename(parser, args)
command.check_required_arguments(parser, opts, "f_low", "waveform")

if not opts.method:
    opts.method.append(default_method)


#
# Logging
#

import logging
logging.basicConfig(level=logging.INFO)
log = logging.getLogger('BAYESTAR')

# LIGO-LW XML imports.
import lal.series
        Option("--min-triggers", type=int, default=2,
            help="Emit coincidences only when at least this many triggers "
            "are found [default: %default]"),
        Option("--measurement-error", choices=measurement_error_choices,
            default=measurement_error_choices[0],
            metavar='|'.join(measurement_error_choices),
            help="How to compute the measurement error [default: %default]"),
        Option("--reference-psd", metavar="PSD.xml[.gz]",
            help="Name of PSD file (required)"),
        Option("--f-low", type=float,
            help="Override low frequency cutoff found in sim_inspiral table.")
    ]
)
opts, args = parser.parse_args()
infilename = command.get_input_filename(parser, args)
command.check_required_arguments(parser, opts, 'reference_psd')


# Python standard library imports.
import os

# LIGO-LW XML imports.
from glue.ligolw import ligolw
from glue.ligolw.utils import process as ligolw_process
from glue.ligolw.utils import search_summary as ligolw_search_summary
from glue.ligolw import table as ligolw_table
from pylal import ligolw_thinca
from glue.ligolw import utils as ligolw_utils
from glue.ligolw import lsctables

# glue, LAL and pylal imports.
parser = OptionParser(
    formatter = command.NewlinePreservingHelpFormatter(),
    description = __doc__,
    usage = '%prog [options] --mass1 MSUN --mass2 MSUN [INPUT.xml[.gz]] [-o OUTPUT.xml[.gz]]',
    option_list = [
        Option("-o", "--output", metavar="OUTPUT.xml[.gz]", default="/dev/stdout",
            help="Name of output file [default: %default]"),
        Option("--mass1", type=float, metavar="MSUN"),
        Option("--mass2", type=float, metavar="MSUN"),
        Option("--snr", type=float, metavar="SNR")
    ]
)
opts, args = parser.parse_args()
infilename = command.get_input_filename(parser, args)
command.check_required_arguments(parser, opts, "mass1", "mass2", "snr")


# Python standard library imports.
import os

# LIGO-LW XML imports.
from glue.ligolw.utils import process as ligolw_process
from glue.ligolw import table as ligolw_table
from glue.ligolw import utils as ligolw_utils
from glue.ligolw import lsctables

# lal imports.
import lal
import lalinspiral.sbank.tau0tau3
import lalsimulation
parser = OptionParser(
    formatter = command.NewlinePreservingHelpFormatter(),
    description = __doc__,
    usage = '%prog [options] --mass1 MSUN --mass2 MSUN [INPUT.xml[.gz]] [-o OUTPUT.xml[.gz]]',
    option_list = [
        Option("-o", "--output", metavar="OUTPUT.xml[.gz]", default="/dev/stdout",
            help="Name of output file [default: %default]"),
        Option("--mass1", type=float, metavar="MSUN"),
        Option("--mass2", type=float, metavar="MSUN"),
        Option("--snr", type=float, metavar="SNR")
    ]
)
opts, args = parser.parse_args()
infilename = command.get_input_filename(parser, args)
command.check_required_arguments(parser, opts, "mass1", "mass2", "snr")


# Python standard library imports.
import os

# LIGO-LW XML imports.
from glue.ligolw.utils import process as ligolw_process
from glue.ligolw import table as ligolw_table
from glue.ligolw import utils as ligolw_utils
from glue.ligolw import lsctables

# lal imports.
import lal
import lalinspiral.sbank.tau0tau3
import lalsimulation
            help="Waveform to use for injections"),
        Option("--snr-threshold", type=float, default=4.,
            help="Single-detector SNR threshold [default: %default]"),
        Option("--min-triggers", type=int, default=2,
            help="Emit coincidences only when at least this many triggers are found [default: %default]"),
        Option("-R", "--repeat-first-injection", type=int, default=None,
            help="Instead of performing each injection once, just perform the first injection this many times."),
        Option("--template-bank", metavar="TMPLTBANK.xml[.gz]",
            help="Name of template bank file (required)"),
        Option("--reference-psd", metavar="PSD.xml[.gz]",
            help="Name of PSD file (required)")
    ]
)
opts, args = parser.parse_args()
infilename = command.get_input_filename(parser, args)
command.check_required_arguments(parser, opts, 'waveform', 'template_bank', 'reference_psd')


# Python standard library imports.
import os
import signal
import sys

# LIGO-LW XML imports.
from glue.ligolw import ligolw
from glue.ligolw.utils import process as ligolw_process
from glue.ligolw.utils import search_summary as ligolw_search_summary
from glue.ligolw import table as ligolw_table
from pylal import ligolw_thinca
from glue.ligolw import utils as ligolw_utils
from glue.ligolw import lsctables
        Option("--min-triggers", type=int, default=2,
            help="Emit coincidences only when at least this many triggers "
            "are found [default: %default]"),
        Option("--measurement-error", choices=measurement_error_choices,
            default=measurement_error_choices[0],
            metavar='|'.join(measurement_error_choices),
            help="How to compute the measurement error [default: %default]"),
        Option("--reference-psd", metavar="PSD.xml[.gz]",
            help="Name of PSD file (required)"),
        Option("--f-low", type=float,
            help="Override low frequency cutoff found in sim_inspiral table.")
    ]
)
opts, args = parser.parse_args()
infilename = command.get_input_filename(parser, args)
command.check_required_arguments(parser, opts, 'reference_psd')


# Python standard library imports.
import os

# LIGO-LW XML imports.
from glue.ligolw import ligolw
from glue.ligolw.utils import process as ligolw_process
from glue.ligolw.utils import search_summary as ligolw_search_summary
from glue.ligolw import table as ligolw_table
from pylal import ligolw_thinca
from glue.ligolw import utils as ligolw_utils
from glue.ligolw import lsctables

# glue, LAL and pylal imports.