Re: Fit statistics

From: Woon-Seng Choong (wschoong@lbl.gov)
Date: Mon Feb 14 2000 - 18:22:22 MET


Hi Rene,
  I just checked our system. I am using 2.23/09. I apologized for the
confusion. Is there a quick fix in 2.23/09 without downloading 2.23/10?

	Thanks
	Seng



On Mon, 14 Feb 2000, Rene Brun wrote:

> Hi Seng,
> In your original mail, you claimed that the eps file had been produced
> by 2.23/10. In the header of the file, I see that it has been produced
> by 2.23/09.
> As I already reported to roottalk, the field for the stats box has been
> extended
> from 32 to 64 characters in 2.23/10.
> 
> Rene Brun
> 
> 
> Woon-Seng Choong wrote:
> > 
> > Hi Rene,
> >   Attached is the eps file. You will notice the fitted parameter
> > p1=3.567e-5 +/- 7.088e-05 but the stat box shows p1=3.567e-5 +/- 7.088e-0.
> > 
> >         Seng
> > 
> > On Sun, 13 Feb 2000, Rene Brun wrote:
> > 
> > > Hi Seng,
> > > I cannot reproduce this problem with the truncated parameters.
> > > Could you send me a gif or ps file illustrating the problem ?
> > >
> > > Rene Brun
> > >
> > > On Fri, 11 Feb 2000, Woon-Seng Choong wrote:
> > >
> > > >
> > > >   I am using v2.23/10 on HPUX10.2 running aCC. I have a macro that plot
> > > > several points and do a linear fit. The error on the fitted parameters is
> > > > truncated in the stat box.  Is there a way to fix this?  Below is
> > > > the macro I used.
> > > >
> > > >     Seng
> > > >
> > > >
> > > > void plot_mombias() {
> > > >   gROOT->Reset();
> > > >   c1 = new TCanvas("c1","HMC Output vs Beam Momentum",200,10,700,500);
> > > >   c1->SetGrid();
> > > >
> > > >   // draw a frame to define the range
> > > >   TH1F *hr1 = c1->DrawFrame(140,-0.005.,190,0.03);
> > > >   hr1->SetXTitle("Average Beam Momentum");
> > > >   hr1->SetYTitle("HMC Estimate - Input Generation Value");
> > > >   hr1->SetTitleOffset(1.3, "Y");
> > > >   hr1->SetTitle("HMC Output (Input g=-0.2154) vs Average Beam Momentum");
> > > >
> > > >   const Int_t n = 8;
> > > >   Float_t x[n]  = {147.577700,152.554231,157.480476,162.457551,167.479134,172.467954,
> > > >                    177.447382,182.435490};
> > > >   Float_t ex[n] = {0, 0, 0, 0, 0, 0, 0, 0};
> > > >   Float_t y[n]  = {.12963e-1,.11759e-1,.06063e-1,.08780e-1,.10120e-1,.10231e-1,
> > > >                    .12242e-1,.10814e-1};
> > > >   Float_t ey[n] = {2.82768e-03, 2.43126e-03,2.17597e-03,2.02725e-03,1.97245e-03,1.94164e-03,
> > > >                    1.97880e-03, 2.05305e-03};
> > > >
> > > >   Float_t chi2[n]  = {344.125, 368.906, 390.534, 304.297, 383.312, 359.423,
> > > >                       417.878, 393.309};
> > > >   Float_t dof[n] = {285, 285, 284, 288, 287, 291, 289, 292};
> > > >
> > > >
> > > >   // create graph
> > > >   gr = new TGraphErrors(n,x,y,ex,ey);
> > > >   gr->SetMarkerColor(kBlack);
> > > >   gr->SetLineColor(kBlack);
> > > >   gr->SetMarkerStyle(20);
> > > >   gr->Draw("P");
> > > >
> > > >
> > > >   gStyle->SetOptFit(111);
> > > >   gStyle->SetStatW(0.2);
> > > >   gStyle->SetStatH(0.1);
> > > >   TF1 *f1 = new TF1("f1", "[0]+[1]*(x-155)", 145, 185);
> > > >   gr->Fit("f1");
> > > >
> > > >   c1->Print("hybrid_mombias.eps");
> > > >
> > > >
> > > > }
> > > >
> > >
> > >
> > 
> >   ------------------------------------------------------------------------
> > 
> >                             Name: hybrid_mombias.eps
> >    hybrid_mombias.eps       Type: Postscript Document (APPLICATION/postscript)
> >                         Encoding: BASE64
> 



This archive was generated by hypermail 2b29 : Tue Jan 02 2001 - 11:50:18 MET