#include <TGraph2D.h>
#include <TTree.h>
#include <TFile.h>
#include <TCanvas.h>
#include <TAxis.h>

void testGraph2D(char* infilename){
	
	//load file and tree
	TFile* infile = new TFile(infilename, "read");
	TTree* conversationData = (TTree*)infile->Get("conversationData");

	//create variables and set branch address
	uint64_t root=0;
	int numNodes=0;
	int diameter=0;
	float density=0;
	float timeFromFirstToLast=0;
	float maxTimeBetweenTweets=0;
	float gap=0;
	
	conversationData->SetBranchAddress("root", &root);
	conversationData->SetBranchAddress("numNodes", &numNodes);
	conversationData->SetBranchAddress("diameter", &diameter);
	conversationData->SetBranchAddress("density", &density);
	conversationData->SetBranchAddress("timeFromFirstToLast", &timeFromFirstToLast);
	conversationData->SetBranchAddress("maxTimeBetweenTweets", &maxTimeBetweenTweets);
	conversationData->SetBranchAddress("gap", &gap);
	
	//create TGraph2D
	TGraph2D *density_gap_diameter = new TGraph2D();

	//loop over tree and fill graph
	int nPoints=0;
	for (int i=0; i<conversationData->GetEntries(); i++){
		conversationData->GetEntry(i);
		if(numNodes>1){			
			density_gap_diameter->SetPoint(nPoints, density, gap, diameter);
			nPoints++;						
		}
	}

	//create canvas, plot and set titles
	TCanvas *c1 = new TCanvas();
	density_gap_diameter->Draw("AP");
	density_gap_diameter->GetXaxis()->SetTitle("density");
	density_gap_diameter->GetYaxis()->SetTitle("gap");
	density_gap_diameter->GetZaxis()->SetTitle("diameter");
	
	c1->Modified();
	c1->Update();
	

	
	
	
	
}