gaussfit.h File Reference

Data Structures

• class gauss1D< T >
• class gauss1D_H< T >
• class gauss1D_2H< T >
• class gauss1D_3H< T >
• class gauss1D_2< T >

Functions

• void gauss1Dfit (ltl::FVector< double, 4 > &parameter, ltl::FVector< double, 4 > &errors, double &chi2, const MArray< double, 1 > &X, const MArray< double, 1 > &Y, const MArray< double, 1 > &Y_error2)
• void gauss1D_H_fit (ltl::FVector< double, 5 > &parameter, ltl::FVector< double, 5 > &errors, double &chi2, const MArray< double, 1 > &X, const MArray< double, 1 > &Y, const MArray< double, 1 > &Y_error2, const ltl::FVector< bool, 5 > &ignore, const int iter=1000, const double tol=1e-7)
• void gauss1D_2H_fit (ltl::FVector< double, 7 > &parameter, ltl::FVector< double, 7 > &errors, double &chi2, const MArray< double, 1 > &X, const MArray< double, 1 > &Y, const MArray< double, 1 > &Y_error2, const ltl::FVector< bool, 7 > &ignore, const int iter=1000, const double tol=1e-7)
• void gauss1D_3H_fit (ltl::FVector< double, 9 > &parameter, ltl::FVector< double, 9 > &errors, double &chi2, const MArray< double, 1 > &X, const MArray< double, 1 > &Y, const MArray< double, 1 > &Y_error2, const ltl::FVector< bool, 9 > &ignore, const int iter=1000, const double tol=1e-7)
• void gauss2Dfit (ltl::FVector< double, 7 > &parameter, ltl::FVector< double, 7 > &errors, double &chi2, MArray< float, 2 > &A, MArray< float, 2 > &A_error, const ltl::FVector< bool, 7 > &ignore, const int iter=200)

---

Function Documentation

void gauss1Dfit	(	ltl::FVector< double, 4 > &	parameter,
		ltl::FVector< double, 4 > &	errors,
		double &	chi2,
		const MArray< double, 1 > &	X,
		const MArray< double, 1 > &	Y,
		const MArray< double, 1 > &	Y_error2
	)

Fit a 1-D gaussian using LTL LMFit.

Referenced by Measure::fit_psf1D().

void gauss1D_H_fit	(	ltl::FVector< double, 5 > &	parameter,
		ltl::FVector< double, 5 > &	errors,
		double &	chi2,
		const MArray< double, 1 > &	X,
		const MArray< double, 1 > &	Y,
		const MArray< double, 1 > &	Y_error2,
		const ltl::FVector< bool, 5 > &	ignore,
		const int	iter = 1000,
		const double	tol = 1e-7
	)

Referenced by Measure::fit_peaks3H_final().

void gauss1D_2H_fit	(	ltl::FVector< double, 7 > &	parameter,
		ltl::FVector< double, 7 > &	errors,
		double &	chi2,
		const MArray< double, 1 > &	X,
		const MArray< double, 1 > &	Y,
		const MArray< double, 1 > &	Y_error2,
		const ltl::FVector< bool, 7 > &	ignore,
		const int	iter = 1000,
		const double	tol = 1e-7
	)

Referenced by Measure::fit_peaks3H(), and Measure::fit_peaks3H_initial().

void gauss1D_3H_fit	(	ltl::FVector< double, 9 > &	parameter,
		ltl::FVector< double, 9 > &	errors,
		double &	chi2,
		const MArray< double, 1 > &	X,
		const MArray< double, 1 > &	Y,
		const MArray< double, 1 > &	Y_error2,
		const ltl::FVector< bool, 9 > &	ignore,
		const int	iter = 1000,
		const double	tol = 1e-7
	)

Referenced by Measure::fit_peaks3H(), and Measure::fit_peaks3H_initial().

void gauss2Dfit	(	ltl::FVector< double, 7 > &	parameter,
		ltl::FVector< double, 7 > &	errors,
		double &	chi2,
		MArray< float, 2 > &	A,
		MArray< float, 2 > &	A_error,
		const ltl::FVector< bool, 7 > &	ignore,
		const int	iter
	)

2-D bivariate gaussian for fitting arc lamp lines in the calibration images

2-D bivariate gaussian for fitting arc lamp lines in the calibration images

Parameters:

	parameter	ltl::FVector of 7 elements: constant, amplitude, sigma_x, sigma_y, relative rotation angle, x_peak, y_peak; needs to be initialised with an initial guess
	errors	ltl::FVector of 7 elements: errors on the seven parameters
	chi2	double chi^2 of the fit
	A	ltl::MArray(): 2D array to fit
	A_error	ltl::MArray(): 2D array with the errors of A
	ignore	ltl::FVector of 7 elements if true don't vary the corresponding parameter
	iter	int maximum number of iterations

Referenced by Measure::fit_gauss2D(), and VirusHealth::validate_peaks().