/*****************************************************************/
/*    NAME: Michael Benjamin                                     */
/*    ORGN: Dept of Mechanical Eng / CSAIL, MIT Cambridge MA     */
/*    FILE: ZAIC_PEAK.h                                          */
/*    DATE: Aug 17th 2006                                        */
/*                                                               */
/* This file is part of IvP Helm Core Libs                       */
/*                                                               */
/* IvP Helm Core Libs is free software: you can redistribute it  */
/* and/or modify it under the terms of the Lesser GNU General    */
/* Public License as published by the Free Software Foundation,  */
/* either version 3 of the License, or (at your option) any      */
/* later version.                                                */
/*                                                               */
/* IvP Helm Core Libs is distributed in the hope that it will    */
/* be useful but WITHOUT ANY WARRANTY; without even the implied  */
/* warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR       */
/* PURPOSE. See the Lesser GNU General Public License for more   */
/* details.                                                      */
/*                                                               */
/* You should have received a copy of the Lesser GNU General     */
/* Public License along with MOOS-IvP.  If not, see              */
/* <http://www.gnu.org/licenses/>.                               */
/*****************************************************************/

#ifdef _WIN32
#pragma warning(disable : 4786)
#pragma warning(disable : 4503)
#endif
#ifndef OF_ZAIC_PEAK_HEADER
#define OF_ZAIC_PEAK_HEADER

#include <string>
#include "IvPDomain.h"
#include "IvPFunction.h"

class PDMap;

class ZAIC_PEAK {
public:
  ZAIC_PEAK(IvPDomain g_domain, const std::string& g_varname);
  virtual ~ZAIC_PEAK() {}

  bool setParams(double summit, double pwidth,  double bwidth, 
		 double delta,  double minutil, double maxutil, 
		 unsigned int index=0);

  bool   setSummit(double, unsigned int index=0);
  bool   setBaseWidth(double, unsigned int index=0);
  bool   setPeakWidth(double, unsigned int index=0);
  bool   setSummitDelta(double, unsigned int index=0);
  bool   setMinMaxUtil(double, double, unsigned int index=0);
  void   setSummitInsist(bool v)      {m_summit_insist=v;}
  void   setValueWrap(bool v)         {m_value_wrap=v;}

  int    addComponent();

  double getParam(std::string, unsigned int index=0);
  bool   getValueWrap()     {return(m_value_wrap);}
  bool   getSummitInsist()  {return(m_summit_insist);}
  int    getSummitCount()   {return(v_summit.size());}

  bool         stateOK()     {return(m_state_ok);}
  std::string  getWarnings() {return(m_warning);}
  IvPFunction* extractOF(bool maxval=true);
  IvPFunction* extractIvPFunction(bool maxval=true) 
  {return(extractOF(maxval));}
  
  IvPDomain getIvPDomain() const {return(m_ivp_domain);}

protected:
  double evalPoint(unsigned int pt_ix, bool maxval=true);
  double evalPoint(unsigned int summit_ix, unsigned int pt_ix);

  void   insistSummit(unsigned int summit_ix);
  PDMap* setPDMap(double tolerance = 0.001);
  
protected:
  std::vector<double> v_summit;
  std::vector<double> v_basewidth;
  std::vector<double> v_peakwidth;
  std::vector<double> v_summitdelta;
  std::vector<double> v_minutil;
  std::vector<double> v_maxutil;

  bool   m_summit_insist;
  bool   m_value_wrap;

private:
  unsigned int m_domain_pts;

  int    m_domain_ix;
  double m_domain_high;
  double m_domain_low;
  double m_domain_delta;

  bool   m_state_ok;

  std::string         m_warning;
  std::vector<double> m_ptvals;
  IvPDomain           m_ivp_domain;
};
#endif