text: user3061: I found this documentation on pyOpenMS mass spec library. I want to try different peak picking algorithms. Unfortunately, the docs are not very detailed and don't show how to use the other peak picking algorithms and the codebase itself is also quite cryptic for me.
  Any idea how to use, for example, the FeatureFindingMetabo algorithm insead of the centroided one?
  from urllib.request import urlretrieve
  # from urllib import urlretrieve  # use this code for Python 2.x
  gh = "https://raw.githubusercontent.com/OpenMS/OpenMS/develop"
  urlretrieve (gh +"/src/tests/topp/FeatureFinderCentroided_1_input.mzML", "feature_test.mzML")
  
  from pyopenms import *
  
  # Prepare data loading (save memory by only
  # loading MS1 spectra into memory)
  options = PeakFileOptions()
  options.setMSLevels([1])
  fh = MzMLFile()
  fh.setOptions(options)
  
  # Load data
  input_map = MSExperiment()
  fh.load("feature_test.mzML", input_map)
  input_map.updateRanges()
  
  ff = FeatureFinder()
  ff.setLogType(LogType.CMD)
  
  # Run the feature finder
  name = "centroided"
  features = FeatureMap()
  seeds = FeatureMap()
  params = FeatureFinder().getParameters(name)
  ff.run(name, input_map, features, params, seeds)
  
  features.setUniqueIds()
  fh = FeatureXMLFile()
  fh.store("output.featureXML", features)
  print("Found", features.size(), "features")
  
  
  user8977: This is a draft answer, if the OP still need needs I will extend it later.
  Me too I found the documentation and examples of pyopenms very scarce. At the same time  OpenMS is huge, sometimes you find key information not in various other sources instead of the official documentation. This is especially true if it comes to metabolomics.
  We built a Python OpenMS lipidomics LC MS/MS preprocessing pipeline, contained in this file:
  https://github.com/saezlab/lipyd/blob/master/src/lipyd/msproc.py
  The top class is the one called MSPreprocess, it implements the workflow, the other classes implement its steps.
  Briefly, you start from profile mode mzML files, then
  
  Peak picking with PeakPickerHiRes results centroided data
  FeatureFindingMetabo assembles the traces in a metabolomics
  compatible way
  A MapAlignment algorithm alignes the traces across
  scans
  A FeatureGrouping algorithm alignes the features across multiple experiments
  
  Our software mentioned above (lipyd) outputs data frames from the m/z, RT and intensity.
  About OpenMS I can recommend to read the C++ code and docs. The Python classes and methods are automatically built minimal wrappers, everything works the same way as in C++.
  
original_text: 3061: <p>I found <a href="https://pyopenms.readthedocs.io/en/latest/feature_detection.html?highlight=featurefinder" rel="nofollow noreferrer">this</a> documentation on pyOpenMS mass spec library. I want to try different peak picking algorithms. Unfortunately, the docs are not very detailed and don't show how to use the other peak picking algorithms and the codebase itself is also quite cryptic for me.</p>
  <p>Any idea how to use, for example, the FeatureFindingMetabo algorithm insead of the centroided one?</p>
  <pre><code>from urllib.request import urlretrieve
  # from urllib import urlretrieve  # use this code for Python 2.x
  gh = &quot;https://raw.githubusercontent.com/OpenMS/OpenMS/develop&quot;
  urlretrieve (gh +&quot;/src/tests/topp/FeatureFinderCentroided_1_input.mzML&quot;, &quot;feature_test.mzML&quot;)
  
  from pyopenms import *
  
  # Prepare data loading (save memory by only
  # loading MS1 spectra into memory)
  options = PeakFileOptions()
  options.setMSLevels([1])
  fh = MzMLFile()
  fh.setOptions(options)
  
  # Load data
  input_map = MSExperiment()
  fh.load(&quot;feature_test.mzML&quot;, input_map)
  input_map.updateRanges()
  
  ff = FeatureFinder()
  ff.setLogType(LogType.CMD)
  
  # Run the feature finder
  name = &quot;centroided&quot;
  features = FeatureMap()
  seeds = FeatureMap()
  params = FeatureFinder().getParameters(name)
  ff.run(name, input_map, features, params, seeds)
  
  features.setUniqueIds()
  fh = FeatureXMLFile()
  fh.store(&quot;output.featureXML&quot;, features)
  print(&quot;Found&quot;, features.size(), &quot;features&quot;)
  </code></pre>
   8977: <p>This is a draft answer, if the OP still need needs I will extend it later.</p>
  <p>Me too I found the documentation and examples of <code>pyopenms</code> very scarce. At the same time  OpenMS is huge, sometimes you find key information not in various other sources instead of the official documentation. This is especially true if it comes to metabolomics.</p>
  <p>We built a Python OpenMS lipidomics LC MS/MS preprocessing pipeline, contained in this file:
  <a href="https://github.com/saezlab/lipyd/blob/master/src/lipyd/msproc.py" rel="nofollow noreferrer">https://github.com/saezlab/lipyd/blob/master/src/lipyd/msproc.py</a></p>
  <p>The top class is the one called <code>MSPreprocess</code>, it implements the workflow, the other classes implement its steps.</p>
  <p>Briefly, you start from profile mode mzML files, then</p>
  <ol>
  <li>Peak picking with PeakPickerHiRes results centroided data</li>
  <li>FeatureFindingMetabo assembles the traces in a metabolomics
  compatible way</li>
  <li>A MapAlignment algorithm alignes the traces across
  scans</li>
  <li>A FeatureGrouping algorithm alignes the features across multiple experiments</li>
  </ol>
  <p>Our software mentioned above (<code>lipyd</code>) outputs data frames from the m/z, RT and intensity.</p>
  <p>About <code>OpenMS</code> I can recommend to read the <code>C++</code> code and docs. The Python classes and methods are automatically built minimal wrappers, everything works the same way as in <code>C++.</code></p>
  