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)
