Hello guys,
I am trying to fit a function to a histogram. Like this:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Double_t myfunction (Double_t *x, Double_t *par){ int bin = data_8175->GetXaxis()->FindBin(x[0]); Double_t f = 0; f= par[0] * ttbar->GetBinContent(bin) + par[1] * h_w_jets-
>GetBinContent(bin) + par[2] * h_others->GetBinContent(bin);
return f;
}
void myfunc(){
TF1 *f1 = new TF1("myfunc",myfunction,minBin_fit,maxBin_fit,3);
f1->SetParLimits(0,0,2); f1->SetParLimits(1,0,2); f1->FixParameter(2,1); f1->Draw();
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Once my histograms are filled. I do the fit this way:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% myfunc(); TF1 *f1= (TF1*)gROOT->GetFunction("myfunc"); data_8175->Fit("myfunc","ER0"); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
I wonder what are the errors that are entered in the chi square. I believe that the errors are the sqrt( number of entries in a bin )
Can I change these errors to be the theoretical ones, i.e, the sqrt of the value of myfunction per bin .
Many thanks in advance,
Jordi Received on Thu Nov 05 2009 - 17:20:26 CET
This archive was generated by hypermail 2.2.0 : Fri Nov 06 2009 - 17:50:03 CET