Source code for has_peak
#!/usr/bin/env python
import pymzml
import os
[docs]
def main():
"""
Testscript to demonstrate functionality of function :py:func:`pymzml.spec.Spectrum.has_peak`
usage:
./has_peak.py
"""
example_file = os.path.join(
os.path.dirname(__file__), os.pardir, "tests", "data", "example.mzML"
)
mz_to_find = 820.7711792
run = pymzml.run.Reader(example_file)
for spectrum in run:
found_peaks = spectrum.has_peak(mz_to_find)
if found_peaks != []:
print("Found peaks: {0} in spectrum {1}".format(found_peaks, spectrum.ID))
if __name__ == "__main__":
main()