{

#include <iostream.h>
#include <stdio.h>
#include <fcntl.h>

  FILE *fp;
  float ptmp,p[20];
  int i, iline=0;

  ntuple = new TNtuple("ntuple","NTUPLE","x:y:z");

  fp = fopen("xyz.dat","r");

  while ( fscanf(fp,"%f",&ptmp) != EOF ){
    p[i++]=ptmp;
    if (i==3){
      i=0; 
      iline++;
      ntuple->Fill(p[0],p[1],p[2]); 
    }
  }

  ntuple.Draw("z:y:x");

}