void SetupAxesg( TGraph2D *gr, TString mainTitle, TString xtitle, TString ytitle, TString ztitle)
{
   gPad -> Update();
   gr -> SetTitle(mainTitle.Data());
   gr -> GetXaxis() -> SetTitle(xtitle.Data());
   gr -> GetXaxis() -> SetTitleSize(0.03);
   gr -> GetXaxis() -> SetLabelSize(0.02);
   gr -> GetXaxis() -> SetDecimals(kTRUE);
   gr -> GetXaxis() -> CenterTitle(true);
   gr -> GetXaxis() -> SetTitleOffset(1.30);
   gr -> GetYaxis() -> SetDecimals(kTRUE);
   gr -> GetYaxis() -> SetTitle(ytitle.Data());
   gr -> GetYaxis() -> SetTitleSize(0.03);
   gr -> GetYaxis() -> SetLabelSize(0.02);
   gr -> GetYaxis() -> CenterTitle(true);
   gr -> GetYaxis() -> SetTitleOffset(1.30);
   gr -> GetZaxis() -> SetDecimals(kTRUE);
   gr -> GetZaxis() -> SetTitle(ztitle.Data());
   gr -> GetZaxis() -> SetTitleSize(0.03);
   gr -> GetZaxis() -> SetLabelSize(0.02);
   gr -> GetZaxis() -> CenterTitle(true);
   gr -> GetZaxis() -> SetTitleOffset(1.30);
}

void rob()
{

   int nn = 10;
   Float_t rAcc[nn][nn];
   Float_t Emax = 52.8304;
   TGraph2D *g1 = new TGraph2D();

   Int_t pnt = 0;
   for( int i=0; i<nn; i++ )
   {
         double xx = Emax*rX0[i];
         for( int j=0; j<nn; j++ )
         {
                 double yy = Emax*rY0[j];
                 g1 -> SetPoint(pnt,xx,yy,rAcc[i][j]);
                 pnt++;
         }
   }
   TCanvas  *c1 = new TCanvas("c1","Radiative Acceptance",70,60,1200,900);
   c1 -> Divide(2,2);
   c1 -> cd(1);
   g1 -> Draw("surf1");
   SetupAxesg( g1, "Radiative Acceptance", "Minimum Positron Energy [MeV]","Minimum Photon Energy [MeV]","Acceptance");
   TVirtualPad* c1_1 = c1->cd(1);
   c1_1 -> SetTheta(30);
   c1_1 -> SetPhi(210);
   c1_1 -> Modified();
   c1_1 -> Update();
}