#include <stdlib.h>
#include <fstream>
#include <iostream>
#include <iomanip>
#include <math.h>
#include "TBenchmark.h"
#include "TFile.h"
#include "TTree.h"
#include "TChain.h"
#include "TH1F.h"

using namespace std;

void ntuple()
{
   gROOT->Reset();
   TFile f("/Volumes/MacHD/out_all_muon.root","recreate");
   TTree data("data"," muon data ");



   ifstream in;
   in.open("/Volumes/MacHD/jdsummary_2440000.dat");

   Double_t jd;
   Int_t runevent,ntrackf,cutf,ntracks,cuts,index;
   Float_t ra1,dec1,ra2,dec2,xlevent,xleventf,xlevents,hr,zen;
   Int_t nlines = 0;
   short levent,leventf,levents;
   data.Branch("runevent",&runevent,"runevent/I");
   data.Branch("ntrackf",&ntrackf,"ntrackf/b");
   data.Branch("ntracks",&ntracks,"ntracks/b");
   data.Branch("cutf",&cutf,"cutf/b");
   data.Branch("cuts",&cuts,"cut/b");
   data.Branch("jd",&jd,"jd/D");
   data.Branch("ra1",&ra1,"ra1/F");
   data.Branch("ra2",&ra2,"ra2/F");
   data.Branch("dec1",&dec1,"dec1/F");
   data.Branch("dec2",&dec2,"dec2/F");
   data.Branch("levent",&levent,"levent/S");
   data.Branch("leventf",&leventf,"leventf/S");   
   data.Branch("levents",&levents,"levents/S");




   for(index=0;index<1000000000;index++) {
      in>>runevent>>jd>>ra1>>dec1>>ra2>>dec2>>xlevent>>xleventf>>xlevents>>hr>>zen>>ntrackf>>cutf>>ntracks>>cuts;
      if (!in.good()) break;
      levent=xlevent+0.5;
      leventf=xleventf+0.5;
      levents=xlevents+0.5;
      if (nlines%500000==0){
      cout<<nlines<<endl;
      data.Print();
};
	  data.Fill();
      nlines++;
   }

   in.close();

   data.Write();
}