{
		
	gROOT->Reset();

#include "Riostream.h"

	ifstream in;
	in.open("field.txt");

	Float_t x,y,z,a,b,c;
	Int_t nlines = 0;
	
	TFile *f = new TFile("MATRIX.root", "RECREATE");
	TH1F *h1= new TH1F("h1", "x distribution", 100, -4, 4);
	TNtuple *ntuple = new TNtuple("ntuple", "data from the MATRIX file", "x:y:z:a,b,c");
	
	while(1) {
		in >> x >> y >> z >> a >> b >> c;
		if (!in.good()) break;
		if (nlines < 0) printf("x=%8f, y=%8f, z=%8fn, a=%8f, b=%8f, c=%8fn", x,y,z,a,b,c);
		h1->Fill(x);
		ntuple->Fill(x,y,z,a,b,c);
		ntuple++;

}
	print("found %d pointsn", nlines);

	in.close();

	f->Write();


}