Loading and plotting arrays

ROOT
1

I have to load in 3 big arrays into ROOT these comprise double precision data and are maybe 1000x10000 in size. I then wish to produce a very simple contour plot using this data in Matlab I could do this using the following code:

load x.dat
load t.dat
load h.dat

contour(x,t,h)

Can anybody tell me how to do this in ROOT please?

Cheers

Dave

2

In order to read the data in I’ve written the following C++ type routine which seems to work (I can read t.dat in the same way):

[color=red]
void Contplotter()
{

//#include<fstream.h>

// Declaration ifstream
ifstream xdata_file;
ifstream hdata_file;

// Declaration local variable
int i(0),j(0);
const Int_t m = 200;
const Int_t n = 10000;
Double_t xd[m][n], hd[m][n];

// Opening file
xdata_file.open("/home/evxdmk/Desktop/Models/SM63/xdata.dat");
hdata_file.open("/home/evxdmk/Desktop/Models/SM63/hdata.dat");

// Read data
i = 0;
while (!xdata_file.eof()) {
j=0;
while (j<n) {

xdata_file >> xd[i][j];
hdata_file >> hd[i][j];
j = j+1;

} //while
i = i+1;
}

// Close file
xdata_file.close();
hdata_file.close();

// Exit and return
return;
}

[/color]

I still do not know how to produce a contour plot of this in ROOT though – any help on this would be much appreciated.

Regards,

Dave