//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
//%%%%%%% RIV analysis step - displays Chi2 results from multiple isotopes on one graph %%%%%%%%%%%%
//%%%%%%%%%%%%%%%%% - must be preceded by running HelioGkChi2 for each of the isotopes of interest %%%%%
//%%%%%%%Version last edited 20100101 - by Russell Leslie%%%%%%%%%%%%%%
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
//
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
//
#include "TF1.h"
#include "TCanvas.h"
#include "TGraph.h"
#include "TGraph2D.h"
#include "TFile.h"
#include "TTree.h"
#include "TLeaf.h"
#include "TH1F.h"
#include "TPaveText.h"
#include "TString.h"
#include "TText.h"
#include "TGraphErrors.h"
#include "TROOT.h"
#include "TStopwatch.h"
#include "Riostream.h"
#include <iostream>
#include <fstream>
#include <string>
//
using namespace std;


 void TestScript() {
//
TStopwatch stopwatch;
//
//
Int_t Choice=0;
Int_t i;
Int_t j;
Int_t k;
Int_t m=80;
Double_t md=80;
//
TCanvas *ChiCanvas;
//
	TGraph2D *gr0 = new TGraph2D();
//
	Double_t Chi2_value=0;
//
	for (i=0;i<m;i++){
		for (j=0;j<m;j++){
			gr0->SetPoint((i*md+j),(i/md),(j/md),(i/md)+(j/md));
			}
		}
	
//
   const UInt_t Number = 2;
   Double_t Red[Number]   = { 0.00, 1.00};
   Double_t Green[Number] = { 0.00, 1.00};
   Double_t Blue[Number]  = { 0.00, 1.00};
   Double_t Stops[Number] = { 0.00, 1.00};
//                                                                               
   Int_t nb=100;
   TColor::CreateGradientColorTable(Number,Stops,Red,Green,Blue,nb);
//
   TCanvas *ChiGroup = new TCanvas("ChiGroup","ChiGroup",10,10,1000,710);
    gr0->SetTitle("ChiGroup");
    gr0->GetZaxis()->SetTitle("-#chi^{2}");
    gr0->GetXaxis()->SetTitle("G_{2}");
    gr0->GetYaxis()->SetTitle("G_{4}");
    gr0->GetXaxis()->SetLabelSize(0.02);
    gr0->GetXaxis()->SetLabelOffset(0.02);
    gr0->GetYaxis()->SetLabelSize(0.02);
    gr0->GetYaxis()->SetLabelOffset(0);
    gr0->GetXaxis()->SetTickLength(0.01);
    gr0->GetYaxis()->SetTickLength(0.01);
    gr0->GetZaxis()->CenterTitle();
    gr0->GetXaxis()->CenterTitle();
    gr0->GetYaxis()->CenterTitle();
    gr0->GetZaxis()->SetTitleOffset(1.20);
    gr0->GetXaxis()->SetTitleOffset(1.45);

    ChiGroup->SetTheta(90);  //<======== 
    ChiGroup->SetPhi(0.00001);     //<======== 
//
    gr0->Draw("surf1");
//
//
stopwatch.Stop();
//
cout << "\nProgram run-time details \nRealtime = " << stopwatch.RealTime() << "\tCputime = " << stopwatch.CpuTime() << "\tLivetime = " << 100*stopwatch.CpuTime()/stopwatch.RealTime() <<"%"<< endl;
}