SqrtMower.h

Go to the documentation of this file.

// Copyright (c) 2002-2023, The OpenMS Team -- EKU Tuebingen, ETH Zurich, and FU Berlin
// SPDX-License-Identifier: BSD-3-Clause
//
// --------------------------------------------------------------------------
// $Maintainer: Mathias Walzer $
// $Authors: $
// --------------------------------------------------------------------------
//
#pragma once
 
#include <OpenMS/KERNEL/StandardTypes.h>
#include <OpenMS/DATASTRUCTURES/DefaultParamHandler.h>
#include <OpenMS/KERNEL/MSSpectrum.h>
#include <OpenMS/KERNEL/MSExperiment.h>
 
#include <cmath>
 
namespace OpenMS
{
  class OPENMS_DLLAPI SqrtMower :
    public DefaultParamHandler
  {
public:
 
    // @name Constructors and Destructors
    // @{
    SqrtMower();
    ~SqrtMower() override;
 
    SqrtMower(const SqrtMower & source);
    SqrtMower & operator=(const SqrtMower & source);
    // @}
 
    template <typename SpectrumType>
    void filterSpectrum(SpectrumType & spectrum)
    {
      bool warning = false;
      for (typename SpectrumType::Iterator it = spectrum.begin(); it != spectrum.end(); ++it)
      {
        double intens = it->getIntensity();
        if (intens < 0)
        {
          intens = 0;
          warning = true;
        }
        it->setIntensity(std::sqrt(intens));
      }
      if (warning)
      {
        std::cerr << "Warning negative intensities were set to zero" << std::endl;
      }
      return;
    }
 
    void filterPeakSpectrum(PeakSpectrum & spectrum);
 
    void filterPeakMap(PeakMap & exp);
 
    //TODO reimplement DefaultParamHandler::updateMembers_()
 
    // @}
 
  };
 
}