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)