format=logging_format, handlers=[RichHandler(rich_tracebacks=True)], ) else: logging.basicConfig( level=logging.INFO, format=logging_format, handlers=[RichHandler(rich_tracebacks=True)], ) logging.info("Input Arguments:-") for arg, value in sorted(vars(values).items()): logging.info("%s: %r", arg, value) your_object = Your(file=values.files) max_delay = your_object.dispersion_delay(dms=np.max(values.dm)) dispersion_delay_samples = np.ceil(max_delay / your_object.your_header.tsamp) logging.info(f"Max Dispersion delay = {max_delay} s") logging.info(f"Max Dispersion delay = {dispersion_delay_samples} samples") if your_object.your_header.nspectra < 2**18: nsamps_gulp = your_object.your_header.nspectra else: nsamps_gulp = int( np.max([(2**np.ceil(np.log2(dispersion_delay_samples))), 2**18])) if values.channel_end > 0: c_max = values.channel_end else: c_max = None
if args.verbose: logging.basicConfig(filename=log_filename, level=logging.DEBUG, format=logging_format) logging.debug("Using debug mode") else: logging.basicConfig(filename=log_filename, level=logging.INFO, format=logging_format) logging.info("Input Arguments:-") for arg, value in sorted(vars(args).items()): logging.info("%s: %r", arg, value) your_object = Your(file=args.files) max_delay = your_object.dispersion_delay(dms=np.max(args.dm)) dispersion_delay_samples = np.ceil(max_delay / your_object.your_header.tsamp) logging.info(f"Max Dispersion delay = {max_delay} s") logging.info(f"Max Dispersion delay = {dispersion_delay_samples} samples") nsamps_gulp = int( np.max([(2**np.ceil(np.log2(dispersion_delay_samples))), 2**18])) your_dada = dada.YourDada(your_object) your_dada.setup() if args.apply_savgol: bandpass = your_object.bandpass(nspectra=8192) mask = savgol_filter(bandpass, your_object.your_header.foff,