{
	TFile f("muon3.root");
	TF1 f1("f1", "[0]+exp([1]+[2]*x)");
	Int_t indf, i = 2;
	Float_t tau, chi2, ndf, cn;
	TH1F *h1[ 10 ];
	char hname[ 100 ];

	TCanvas c1("c1","c1",2);
	c1.SetFillColor( 10 );

	//for( indf = 4; indf < 50; indf += 1 )
	indf = 30;
	{		
	  //sprintf( hname, "h%d", i );
	  sprintf( hname, "TOF" );
		h1[i] = new TH1F(hname, "Binned measurements from muon experiment", indf, 0.6, 10);
		mlout->Project(hname, "time");
		h1[i]->Fit("f1","QEM");
		h1[i]->GetXaxis()->SetTitle("Measured lifetime (uS)");
		h1[i]->GetXaxis()->SetTitleOffset(1.2);
		h1[i]->GetXaxis()->SetTitleSize(0.03);
		h1[i]->GetYaxis()->SetTitle("Number of events                                                                                    ");
		h1[i]->GetYaxis()->SetTitleOffset(1.6);
		h1[i]->GetYaxis()->SetTitleSize(0.03);
		h1[i]->Draw("ep");
		i++;

		tau = -1.0 / f1.GetParameter(2);
		chi2 = f1.GetChisquare();
		ndf = f1.GetNDF();
		cn = chi2 / ndf;

		printf("indf=%d, ndf=%f, tau=%f, chi2=%f, chi2/ndf=%f", indf, ndf, tau, chi2, cn);	
		if( cn < 1.33 ) printf (" ***");
		printf("\n");
	}
}