Hi,
I have a problem with fitting using TGraphErrors / TF1. I have the
following code:
Double_t efffun(Double_t *X, Double_t *pars)
{
// par[] = { 0=A, 1=B, 2=C, 3=D, 4=E, 5=F, 6=G, 7=E1, 8=E2, 9=scale }
Double_t x1 = /*log*/(X[0] / pars[7]);
Double_t x2 = /*log*/(X[0] / pars[8]);
if( x1 > exp(-100) ) x1 = log( x1 ); else x1 = -100;
if( x2 > exp(-100) ) x2 = log( x2 ); else x2 = -100;
Double_t g = pars[6];
Double_t f1 = pars[0] + pars[1]*x1 + pars[2]*x1*x1;
Double_t f2 = pars[3] + pars[4]*x2 + pars[5]*x2*x2;
return pars[9]*exp( pow( pow( f1, -g ) + pow( f2, -g ), -1.0 / g));
}
TF1* func = new TF1("eff", efffun, 200, 2500, 10);
func->SetParNames("A", "B", "C", "D", "E", "F", "G", "E1", "E2", "scale");
func->SetParameters(2.9, 2.2, 0, 5.179, -0.302, -0.17, 15, 100, 1000, 1);
func->FixParameter( 2, 0 ); // C=0
//func->FixParameter( 6, 15 ); // G=15
func->FixParameter( 7, 100 ); // E1 = 100
func->FixParameter( 8, 1000 ); // E2 = 1000
g->Fit( func, "" );
The graph has these values:
x[0]=244.454, y[0]=121.118, ex[0]=2.44454, ey[0]=1.21118
x[1]=344.061, y[1]=178.907, ex[1]=3.44061, ey[1]=1.78907
x[2]=367.594, y[2]=189.898, ex[2]=3.67594, ey[2]=1.89898
x[3]=410.915, y[3]=192.414, ex[3]=4.10915, ey[3]=1.92414
x[4]=443.769, y[4]=198.059, ex[4]=4.43769, ey[4]=1.98059
x[5]=688.482, y[5]=185.509, ex[5]=6.88482, ey[5]=1.85509
x[6]=778.677, y[6]=189.777, ex[6]=7.78677, ey[6]=1.89777
x[7]=867.135, y[7]=186.262, ex[7]=8.67135, ey[7]=1.86262
x[8]=963.785, y[8]=182.018, ex[8]=9.63785, ey[8]=1.82018
x[9]=1111.76, y[9]=172.064, ex[9]=11.1176, ey[9]=1.72064
x[10]=1212.64, y[10]=166.903, ex[10]=12.1264, ey[10]=1.66903
x[11]=1298.83, y[11]=157.176, ex[11]=12.9883, ey[11]=1.57176
x[12]=1407.69, y[12]=157.646, ex[12]=14.0769, ey[12]=1.57646
Now the problem is that before fitting, the TF1 fits really well (I
adjusted the initial values), but after fitting it's "worse", I mean it
fits anything but not the graph (it looks a bit like the fitter is trying
to put the TF1 to 0 outside the graph).
How can I solve this problem?
Best wishes,
Alexander
This archive was generated by hypermail 2b29 : Sun Jan 02 2005 - 05:50:10 MET