//
//  main.cpp
//  data analysis
//
//  Created by Fatimah awdhah on 7/30/15.
//  Copyright (c) 2015 Fatimah awdhah. All rights reserved.
//

#include "Riostream.h"
#include "TFile.h"
#include "TH1.h"
#include "TNtuple.h"


void basic(const char* fdata="basic.dat", const char* froot="basic.root")
{
    //example of macro to read data from an ascii file and create a root file with an histogram and an ntuple.
    //arguments:
    //fdata= nom du fichier ascii a lire (default:basic.dat)
    //froot= nom du fichier ROOT resultat (default:basic.root)
    
    ifstream in;
    //we assume a file basic.dat in the current directory
    //this file has 3 columns of float data
    in.open(fdata, ios::in);
    float x,y,z;
    int nlines =0.;
    TFile* f = new TFile (froot, "RRECREATE");
    TH1F *h1 = new TH1F("h1", "X Distribution",100,-4,4);
    TNtuple* ntuple= new TNtuple("ntuple","data from ascii file","x:y:z");
    
    while (in.good()) {
        in >>x >>y >>z;
        if (nlines <5) {
            cout<< "X= " <<x<< ", Y= "<< y<<", Z= "<<z<< endl;
        }
        h1->Fill(x);
        ntuple->Fill (x,y,z);
        nlines++;
    }
    cout<<"on a trouve " <<nlines << "lignes. " <<endl;
    in.close();
    f->Write ();
    delete f;
  
}