Skip to content
Snippets Groups Projects
Commit 299c6a0a authored by Patrick Godwin's avatar Patrick Godwin
Browse files

gstlal_inspiral: avoid issue when calculating expected SNR + frame caching,...

gstlal_inspiral: avoid issue when calculating expected SNR + frame caching, lower sampling rate of expected SNR calculation
parent 3c58956d
No related branches found
No related tags found
1 merge request!41DAG Workflow Overhaul + OSG DAG support
......@@ -772,7 +772,7 @@ for output_file_number, (svd_bank_url_dict, output_url, ranking_stat_output_url,
print("calculating expected SNR for injections...", file=sys.stderr)
# read in injections
xmldoc = ligolw_utils.load_filename(detectors.injection_filename, verbose = options.verbose, contenthandler = LIGOLWContentHandler)
xmldoc = ligolw_utils.load_filename(options.injections, verbose = options.verbose, contenthandler = LIGOLWContentHandler)
sim_inspiral_table = lsctables.SimInspiralTable.get_table(xmldoc)
# calculate expected SNR
......@@ -781,12 +781,12 @@ for output_file_number, (svd_bank_url_dict, output_url, ranking_stat_output_url,
f_low = min(set(bank.f_low for ifo in banks.keys() for bank in banks[ifo]))
max_rate = max(set(rate for ifo in banks.keys() for bank in banks[ifo] for rate in bank.get_rates()))
f_high = 0.8 * (max_rate / 2.)
inspiral.calc_sim_inspiral_table_snrs(sim_inspiral_table, psd, detectors.seg, f_low = f_low, f_high = f_high)
inspiral.calc_sim_inspiral_table_snrs(sim_inspiral_table, psd, detectors.seg, f_low = f_low, f_high = f_high, sample_rate = max_rate)
# write the injection file back out to a temporary directory
f, fname = tempfile.mkstemp(".xml.gz")
ligolw_utils.write_filename(xmldoc, fname, verbose = options.verbose, gz = True)
detectors.injection_filename = fname
options.injections = fname
os.close(f)
......@@ -889,7 +889,7 @@ for output_file_number, (svd_bank_url_dict, output_url, ranking_stat_output_url,
instruments = rankingstat.instruments,
seg = detectors.seg,
offsetvectors = offsetvectors,
injection_filename = detectors.injection_filename,
injection_filename = options.injections,
tmp_path = options.tmp_space,
replace_file = True,
verbose = options.verbose
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment