////////////////////////////////////////////////////////////////////////////////
//
// graph2d
// -------
//
// 12/13/2005 Philipp Schieferdecker <schieferd@lmu.de>
////////////////////////////////////////////////////////////////////////////////

int graph2d()
{
    ifstream fin("graph2d.txt");
    if (!fin) { cout<<"ERROR. Abort."<<endl; return 0;}

    int njes;
    int nmtop;
    fin>>njes>>nmtop;
    
    double*  jes =new double[njes];
    double*  mtop=new double[nmtop];
    double** lnl =new double*[njes];
    for (int i=0;i<njes;i++) lnl[i]=new double[nmtop];
    
    for (int i=0;i<njes;i++)  fin>>jes[i];
    for (int i=0;i<nmtop;i++) fin>>mtop[i];
    
    for (int i=0;i<njes;i++)
	for (int j=0;j<nmtop;j++) fin>>lnl[i][j];

    fin.close();

    
    TGraph2D* g=new TGraph2D(njes*nmtop);
    for (int i=0;i<njes;i++)
	for (int j=0;j<nmtop;j++)
	    g->SetPoint(i*nmtop+j,jes[i],mtop[j],lnl[i][j]);
    

    TCanvas*c00=new TCanvas("c00","c00",0,0,750,500);
    g->Draw("AP");
    g->SetMarkerStyle(kFullCircle);
    
    TF2* fit1=new TF2("fit1","[0]+[1]*x+[2]*y+[3]*x*y+[4]*x*x+[5]*y*y",
		      jes[2],jes[njes-3],mtop[2],mtop[nmtop-3]);
    fit1->SetParameters(1,1,1,1,1,1);
    g->Fit(fit1,"R");
    fit1->Draw("SURFSAME");
    
    TH1F* h1=new TH1F("h1","h1",100,-10,10);
    for (int i=2;i<njes-2;i++)
	for (int j=2;j<nmtop-2;j++)
	    h1->Fill(lnl[i][j]-fit1->Eval(jes[i],mtop[j]));
    
    TCanvas* c01=new TCanvas("c01","c01",800,0,750,500);
    h1->Draw();

    
    return 0;
}