Commit 8d0a1771 authored by Qi Chu's avatar Qi Chu

change hard-coded chisq ratio thresh to be given by option --finalsink-chisq-ratio-thresh

parent 005295a4
......@@ -210,6 +210,7 @@ def parse_command_line():
parser.add_option("--job-tag", help = "Set the string to identify this job and register the resources it provides on a node. Should be 4 digits of the form 0001, 0002, etc. required")
parser.add_option("--finalsink-far-factor", type = "float", default = 1.0, help = "normalization factor for false alarm rate threshold. = number of nodes. ")
parser.add_option("--finalsink-singlefar-veto-thresh", type = "float", default = 0.01, help = "at least two detectors have fars less than this single far veto threshold. ")
parser.add_option("--finalsink-chisq-ratio-thresh", type = "float", default = 8, help = "any two detectors have chisq ratio less than this threshold. ")
parser.add_option("--finalsink-superevent-thresh", type = "float", default = 3.8e-7, help = "the superevent thresh which we will apply all kinds of vetoes. ")
parser.add_option("--finalsink-gracedb-far-threshold", type = "float", default = None, help = "false alarm rate threshold for gracedb (Hz), if not given gracedb events are not sent")
parser.add_option("--finalsink-gracedb-group", default = "Test", help = "gracedb group, default is Test")
......@@ -483,6 +484,7 @@ finalsink = postcoh_finalsink.FinalSink(
fapupdater_output_fname = options.finalsink_fapupdater_output_fname,
fapupdater_collect_walltime_string = options.finalsink_fapupdater_collect_walltime,
singlefar_veto_thresh = options.finalsink_singlefar_veto_thresh,
chisq_ratio_veto_thresh = options.finalsink_chisq_ratio_thresh,
superevent_thresh = options.finalsink_superevent_thresh,
gracedb_far_threshold = options.finalsink_gracedb_far_threshold,
gracedb_group = options.finalsink_gracedb_group,
......
......@@ -326,7 +326,7 @@ class FAPUpdater(object):
collected_fnames.append("%s/%s" % (self.path, one_bank_fname))
class FinalSink(object):
def __init__(self, channel_dict, process_params, pipeline, need_online_perform, path, output_prefix, output_name, far_factor, cluster_window = 0.5, snapshot_interval = None, fapupdater_interval = None, cohfar_accumbackground_output_prefix = None, cohfar_accumbackground_output_name = None, fapupdater_output_fname = None, fapupdater_collect_walltime_string = None, singlefar_veto_thresh = 0.01, gracedb_far_threshold = None, gracedb_group = "Test", gracedb_search = "LowMass", gracedb_pipeline = "spiir", gracedb_service_url = "https://gracedb.ligo.org/api/", gracedb_offline_annote = None, output_skymap = 0, superevent_thresh = 3.8e-7, verbose = False):
def __init__(self, channel_dict, process_params, pipeline, need_online_perform, path, output_prefix, output_name, far_factor, cluster_window = 0.5, snapshot_interval = None, fapupdater_interval = None, cohfar_accumbackground_output_prefix = None, cohfar_accumbackground_output_name = None, fapupdater_output_fname = None, fapupdater_collect_walltime_string = None, singlefar_veto_thresh = 0.01, chisq_ratio_veto_thresh = 8.0, gracedb_far_threshold = None, gracedb_group = "Test", gracedb_search = "LowMass", gracedb_pipeline = "spiir", gracedb_service_url = "https://gracedb.ligo.org/api/", gracedb_offline_annote = None, output_skymap = 0, superevent_thresh = 3.8e-7, verbose = False):
#
# initialize
#
......@@ -345,7 +345,7 @@ class FinalSink(object):
self.need_candidate_check = False
self.cur_event_table = lsctables.New(postcoh_table_def.PostcohInspiralTable)
# FIXME: hard-coded chisq_ratio_thresh to veto
self.chisq_ratio_thresh = 8
self.chisq_ratio_thresh = chisq_ratio_veto_thresh
self.superevent_thresh = superevent_thresh
self.nevent_clustered = 0
self.singlefar_veto_thresh = singlefar_veto_thresh
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment