[fittedPars, chisq] = pwFit Fits the currently combined couples.
If no custom optimization function is supplied, the Likelihood is maximized which is equivalent for Gaussian errors to minimize the Chi-Square value defined as Chi-Square = sum_i ((yModel(i) - yData(i))/yStd(i))^2 for all available datapoints i.