code: def _get_window_thresh(n_times, sfreq, mt_bandwidth, p_value):
      from .time_frequency.multitaper import _compute_mt_params
  
      # figure out what tapers to use
      window_fun, _, _ = _compute_mt_params(
          n_times, sfreq, mt_bandwidth, False, False, verbose=False
      )
  
      # F-stat of 1-p point
      threshold = fstat.ppf(
          1 - p_value / n_times, 2, 2 * len(window_fun) - 2
      )
      return window_fun, threshold
  
imports: from collections import Counter
  from copy import deepcopy
  from functools import partial
  from math import gcd
  import numpy
  from scipy import fft, signal
  from scipy.stats import f
  from ._fiff.pick import _picks_to_idx
  from ._ola import _COLA
  from .cuda import _fft_multiply_repeated, _fft_resample, _setup_cuda_fft_multiply_repeated, _setup_cuda_fft_resample, _smart_pad
  from .parallel import parallel_func
  from .utils import _check_option, _check_preload, _ensure_int, _pl, _validate_type, logger, sum_squared, verbose, warn
label_desc: Use image node_id 18 for calling a global function with example usage: _get_window_thresh(n_times, sfreq, mt_bandwidth, p_value) and returns: window_fun,  threshold