#include <iostream>
#include <sstream>
#include <fstream>

#include "TROOT.h"
#include "TAxis.h"
#include "TCanvas.h"
#include "TStyle.h"
#include "TApplication.h"
#include "TGraph2D.h"

using namespace std;

int main( int argc, char* argv[] ){
   if(argc < 2){
      cout << "Usage: plot_SVN-MC parameterAndIntegralsFile " << endl;
      return -1;
   }
   int argc2 = 0;
   char* argv2[10];
   TApplication *app = new TApplication("App", &argc2, argv2 );
   string parameterFileName="0";
   string rejectionEfficiencyFileListName="0";
   for(int i=1; i<argc; i++){
   if( i==1 ){parameterFileName = argv[1];
      }
      else{ 
	 cout << "Too many arguments!" << endl; 
	 cout << "Usage: plot_SVN-MC parameterAnIntegralsFile" << endl;
	 return -1;
      }
   }
   
   // Read data from file
   ifstream parameterFile;
   parameterFile.open( parameterFileName.c_str() );

   int sseDiffGraphCount = 0;
   TGraph2D* sseDiffGraph = new TGraph2D();
   sseDiffGraph -> SetTitle( "C, Gamma" );
   
      if( parameterFile.good() ){
	 double c, gamma, diffAll, diffSse, diffMse;

	 while( (parameterFile >> c >> gamma >> diffAll >> diffSse >> diffMse) ){
	    if( diffSse < 10000 && diffSse >-10000 ){	cout << c << "	" << gamma << "	" << diffAll << "	" << diffSse << "	" << diffMse << endl;	}
	    sseDiffGraph -> SetPoint( sseDiffGraphCount, gamma, c, diffSse );
	    sseDiffGraphCount ++;
	 }
      }
      
   // Draw Graph and save as .png and .root
   gStyle -> SetPalette(1);
   gStyle -> SetTitleOffset(1.9, "xy");
   gStyle -> SetTitleOffset(1.35, "z");

   if( 0 != sseDiffGraph -> GetN() ){
      TCanvas* sseDiffCanv = new TCanvas( "sseDiffGraph" , "" , 600 , 500 ) ;
      sseDiffCanv -> cd();
      sseDiffCanv -> SetLogx();
      sseDiffCanv -> SetLogy();
      
      sseDiffGraph -> SetMarkerStyle(20);
      sseDiffGraph -> Draw("PCOL");
      
      sseDiffGraph -> GetXaxis() -> SetTitle("Gamma");
      sseDiffGraph -> GetYaxis() -> SetTitle("C");
      sseDiffGraph -> GetZaxis() -> SetTitle("Data SSE - MC SSE");
   //     sseDiffGraph -> GetZaxis() -> SetRangeUser( 0.1 , 0.15 );
//       sseDiffGraph -> Draw("PCOL");
      sseDiffGraph -> Draw("same surf1");
//       sseDiffCanv -> Update();
      //       sseDiffGraph -> Draw("P");
//       sseDiffGraph -> Draw("same tri1");
      sseDiffCanv -> SaveAs("sseDiffGraph.png");
      sseDiffCanv -> SaveAs("sseDiffGraph.root");
   }
   
   cout << "End of main. " << endl;
   app->Run(kTRUE);
   return 0;
}