Photometry.h

#ifndef Photometry_h
#define Photometry_h
 
/* SPDX-License-Identifier: CC0-1.0 */
 
#include <string>
#include <vector>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_min.h>
#include <gsl/gsl_roots.h>
 
namespace Isis {
  class Pvl;
  class PhotoModel;
  class AtmosModel;
  class NormModel;
  class Photometry {
    public:
      Photometry(Pvl &pvl);
      Photometry() {};
      virtual ~Photometry();
 
      void Compute(double pha, double inc, double ema, double dn,
                   double &albedo, double &mult, double &base);
      void Compute(double pha, double inc, double ema, double deminc,
                   double demema, double dn, double &albedo,
                   double &mult, double &base);
 
      virtual void SetPhotomWl(double wl);
 
      static void minbracket(double &xa, double &xb, double &xc, double &fa,
          double &fb, double &fc, double Func(double par, void *params),
          void *params);
 
      static int brentminimizer(double x_lower, double x_upper, gsl_function *Func, 
          double & x_minimum, double tolerance);
      
      static int brentsolver(double x_lo, double x_hi, gsl_function *Func, double tolerance, double &root);
 
      PhotoModel *GetPhotoModel() const {
        return p_phtPmodel;
      }
 
      AtmosModel *GetAtmosModel() const {
        return p_phtAmodel;
      }
 
      NormModel *GetNormModel() const {
        return p_phtNmodel;
      }
 
    protected:      
      AtmosModel *p_phtAmodel;
      PhotoModel *p_phtPmodel;
      NormModel *p_phtNmodel;
  };
};
 
#endif