def find_doppler(filename_fil, kernels):
r""" Run turboseti doppler search on filename with default params """
print("\n===== find_doppler =====")
print("Searching %s" % filename_fil)
filename_dat = filename_fil.replace('.h5', '.dat')
filename_log = filename_fil.replace('.h5', 'log')
if os.path.exists(filename_dat):
os.remove(filename_dat)
if os.path.exists(filename_log):
os.remove(filename_log)
snr = 5.0
coarse_chans = ''
obs_info = None
n_coarse_chan = 1
max_drift = 1.0
find_seti_event = FindDoppler(filename_fil,
max_drift=max_drift,
snr=snr,
out_dir=HERE,
coarse_chans=coarse_chans,
obs_info=obs_info,
n_coarse_chan=n_coarse_chan,
kernels=kernels)
t0 = time.time()
find_seti_event.search()
t_taken = time.time() - t0
print(
"Time taken for find_seti_event.search() [without flagging]: %2.2fs" %
t_taken)
find_seti_event.flagging = True
t0 = time.time()
for dummy, data_dict in enumerate(find_seti_event.data_handle.data_list):
search_coarse_channel(data_dict,
find_seti_event,
filewriter=None,
logwriter=None)
t_taken = time.time() - t0
print("Time taken for search_coarse_channe() [with flagging]: %2.2fs" %
t_taken)
def find_doppler(filename_fil, kernels):
r""" Run turboseti doppler search with and without flagging """
print("\n===== find_doppler =====")
print("Searching %s" % filename_fil)
filename_dat = filename_fil.replace('.h5', '.dat')
filename_log = filename_fil.replace('.h5', 'log')
if os.path.exists(filename_dat):
os.remove(filename_dat)
if os.path.exists(filename_log):
os.remove(filename_log)
find_seti_event = FindDoppler(filename_fil,
max_drift=MAX_DRIFT,
snr=MIN_SNR,
out_dir=TESTDIR,
obs_info=None,
n_coarse_chan=1,
kernels=kernels)
t0 = time.time()
find_seti_event.search()
t_taken = time.time() - t0
print(
"Time taken for find_seti_event.search() [without flagging]: %2.2fs" %
t_taken)
find_seti_event.flagging = True
t0 = time.time()
for dummy, data_dict in enumerate(find_seti_event.data_handle.data_list):
search_coarse_channel(data_dict,
find_seti_event,
filewriter=None,
logwriter=None)
t_taken = time.time() - t0
print("Time taken for search_coarse_channe() [with flagging]: %2.2fs" %
t_taken)