fitting.h

Go to the documentation of this file.

#ifndef ATOMPROBE_FITTING_H
#define ATOMPROBE_FITTING_H


namespace AtomProbe
{


bool fitVoigt(const std::vector<double> &x, const std::vector<double> &y, 
    double &sigma, double &gamma, double &mu, double &amp, bool autoInit=true);


void voigtProfile(const std::vector<double> &x, double sigma, 
    double gamma, double mu, double amp, std::vector<double> &y);

bool fitDoniachSunjic(const std::vector<double> &x, const std::vector<double> &y, 
    double &a, double &f, double &mu, double &amp, bool autoInit=true);

void doniachSunjic(const std::vector<double> &x, double &a, double &F, double &mu,double &amp,
                std::vector<double> &y);

bool fitLikeLogGaussian(const std::vector<double> &x, const std::vector<double> &y, 
    double &lambda, double &sigma, double &xp, double &amp, double &h,bool autoInit=true);


void likeLogGaussian(const std::vector<double> &x, double &xp, double &sigma,
                double &lambda,double &amp, double &h,
                std::vector<double> &y);

bool fitExpNorm(const std::vector<double> &x, const std::vector<double> &y, 
    double &K,  double &mu, double &sigma,double &amp, bool autoInit=true);

/* mu - shift.
   sigma - squeeze
  K - asymmetry
  amp - amplitude
*/
void expNorm(const std::vector<double> &x, double &K, double &mu, double&sigma, double &amp, std::vector<double> &y);

double lsq(const std::vector<double>  &y, const std::vector<double> &yFit);
}
#endif