def _get_intensity_threshold(self, scan): intensity_array = scan.intensity_array intensity_threshold = tic(scan) * self.percent_tic intensity_array = intensity_array[intensity_array >= intensity_threshold] quantile = array_split(sort(intensity_array), self.q)[self.q - self.k] if len(quantile) > 0: return quantile[0], quantile[-1] else: return float("inf"), float("-inf")
def _get_intensity_threshold(self, scan): return (tic(scan) * self.percent_tic, float("inf"))
def calculate(self, spectra, psm): return tic(spectra)