This is an old revision of the document!

Background Polynomial

Chebyshev polynomials of the first kind as described at is used by default in Topas. The recurrences relations of Eq. (26) can be used with the x axis normalized to be between -1 and 1.

A Chebychev Polynomial, say order 5, written as a fit_obj can be coded as:

prm c0  359.79390`
prm c1  140.27403`
prm c2  67.32100`
prm c3  21.05784`
prm c4  11.34432`
prm c5 -4.46399`
prm !xm = X2 - X1;
prm !xp = X2 + X1;
prm x = (2 X - xp) / xm;
prm tn2 = 2 x x - 1;
prm tn3 = 2 x tn2 - x;
prm tn4 = 2 x tn3 - tn2;
prm tn5 = 2 x tn4 - tn3;
fit_obj = 
	c0 +
	c1 x +
	c2 tn2 +
	c3 tn3 +
	c4 tn4 +
	c5 tn5 

Note that X1 and X2 are reserved parameter names that correspond to the start and end of the x-axis.

The attached AAC.INP and data file demonstrates that the fit_obj coefficients after refinement correspond exactly to those of the bkg Chebychev Polynomial.

Personal Tools