peakList.h

Go to the documentation of this file.

// -*- Mode: C++; tab-width: 2; -*-
// vi: set ts=2:
//
// $Id: peakList.h,v 1.17 2003/08/26 08:04:45 oliver Exp $
//

#ifndef BALL_NMR_PEAKLIST_H
#define BALL_NMR_PEAKLIST_H

#ifndef BALL_NMR_PEAK_H
# include<BALL/NMR/peak.h>
#endif

#ifndef BALL_COMMON_LIMITS_H
# include<BALL/COMMON/limits.h>
#endif

#ifndef BALL_DATATYPE_LIST_H
# include<BALL/DATATYPE/list.h>
#endif

namespace BALL 
{

  template <typename PT>
  class PeakList
    : public List<PT>
  {
    public:
    
    BALL_CREATE(PeakList<PT>)

    
      
    
    typedef typename PT::Position Position;

    typedef PT PeakType;
  
    typedef typename List<PT>::iterator Iterator;

    typedef typename List<PT>::ConstIterator ConstIterator;


    PeakList() 
    {
    }
    
    PeakList(const PeakList& peak_list)
      : List<PT>(peak_list)
    {
    }   

    virtual ~PeakList()
    {
    }
    

      
    void scale(float x)
    {
      Iterator it = List<PT>::begin();
      for (; it != List<PT>::end(); ++it)
      {
        it->setIntensity(it->getIntensity() * x);
      }
    }

    float getMaxIntensity() const
    {
      ConstIterator it = List<PT>::begin();
      float max = -Limits<float>::max();
      for (; it != List<PT>::end(); ++it)
      {
        max = std::max(max, it->getIntensity());
      }
    
      return max;
    }

    float getMinIntensity() const
    {
      ConstIterator it = List<PT>::begin();
      float min = Limits<float>::max();
      for (; it != List<PT>::end(); ++it)
      {
        min = std::min(min, it->getIntensity());
      }
    
      return min;
    }

  };

    
  typedef PeakList<Peak<float> > PeakList1D;

  typedef PeakList<Peak<Vector2> > PeakList2D;

  typedef PeakList<Peak<Vector3> > PeakList3D;
  
  
} // namespace BALL

#endif // BALL_NMR_PEAKLIST_H