def test_spectrum_ParentPeakMower(self): thisfilter = pyopenms.ParentPeakMower() new_firstspec = self.exp[1] thisfilter.filterSpectrum(new_firstspec) self.assertNotEqual(new_firstspec.size(), 0) self.assertNotEqual(new_firstspec, self.exp[1]) # TODO self.assertEqual(new_firstspec[10].getMZ(), self.exp[1][10].getMZ()) self.assertEqual(new_firstspec[10].getIntensity(), self.exp[1][10].getIntensity())
def test_map_ParentPeakMower(self): thisfilter = pyopenms.ParentPeakMower() self.exp[0].setMSLevel(2) old_firstspec = self.exp[0] old_secondspec = self.exp[1] thisfilter.filterPeakMap(self.exp) # the MS1 should not be processed self.assertNotEqual(self.exp.size(), 0) self.assertEqual(old_firstspec, self.exp[0]) self.assertNotEqual(old_secondspec, self.exp[1]) # TODO self.assertEqual(old_secondspec[10].getMZ(), self.exp[1][10].getMZ()) self.assertEqual(old_secondspec[10].getIntensity(), self.exp[1][10].getIntensity())
def test_map_ParentPeakMower(self): thisfilter = pyopenms.ParentPeakMower() param = thisfilter.getDefaults() param.setValue(b"window_size", 20.0, b"") param.setValue(b"set_to_zero", True, b"") thisfilter.setParameters(param) self.assertEqual(self.exp.size(), 5) self.assertEqual(self.exp[0].getMSLevel(), 1) map = [] for spec in self.exp: if spec.getMSLevel() != 1: map.append(spec) self.exp.setSpectra(map) self.assertEqual(self.exp.size(), 3) self.assertEqual(self.exp[0].getMSLevel(), 2) self.assertEqual(self.exp[1].getMSLevel(), 2) self.assertEqual(self.exp[2].getMSLevel(), 2) old_firstspec = self.exp[0] old_secondspec = self.exp[1] thisfilter.filterPeakMap(self.exp) a = old_firstspec b = self.exp[0] self.assertEqual(a[174].getIntensity(), b[174].getIntensity()) self.assertNotEqual(a[175].getIntensity(), b[175].getIntensity()) self.assertNotEqual(a[201].getIntensity(), b[201].getIntensity()) self.assertEqual(a[202].getIntensity(), b[202].getIntensity()) # All peaks between 174 and 202 should be zero for i in range(175, 202, 1): self.assertEqual(b[i].getIntensity(), 0) self.assertEqual(a[i].getMZ(), b[i].getMZ())