{
TFile *myFile = new TFile("twobody.root","READ");
myFile->ls();

TCanvas *myCanvas = new TCanvas("c","twobody",0,0,600,900);
myCanvas->Divide(2,3);

TH1F *th1=h1;
TH1F *th2=h2;
TH1F *th3=h3;

TF1 *f1=new TF1("f1", "[0]*(1.0+[1]*x+[2]*x*x)/(2+[2])",-1.0,1.0);

myCanvas->cd(1);
th1->Fit("f1");
myCanvas->cd(2);
th2->Fit("f1");
myCanvas->cd(3);
th3->Fit("f1");


myCanvas.Update();

myCanvas->Print("twobody4.eps");

}