/* $Id: myplotting.cc 19 2009-07-16 22:28:31Z yangsuli $
   $Date: 2009-07-17 00:28:31 +0200 (Fri, 17 Jul 2009) $
 */

#include <TMath.h>
#include <TCanvas.h>
#include <TCanvas.h>
#include <TFile.h>
#include <TH1F.h>
#include <TTree.h>
#include <TChain.h>
#include <TF1.h>

#include <stdlib.h>
#include <stdio.h>
#include <iostream>
using namespace std;
#include <sys/stat.h>
#include <vector>
#include "../aux/OutputTree_.h"
#define EVE 10000000

OutputTree_ *p;
int  NEvent;

void get_event(int i) {
	if ( p->LoadTree(i) < 0) { 
		cout<<"\nProblem in LoadTree."
			<<"\nEntry: "<<i<<endl;
		exit(0);
	}
	p->fChain->GetEntry(i);
}

void make_histgram() {


	NEvent = p->fChain->GetEntries();
	cout<<"Found "<<NEvent<<" events.\n"; 

	int counter = 0;
	for (int i=0; i<NEvent; i++) {
		get_event(i);
		if(fabs(p->JetInvMasses_eta1stJ)>0.5&&fabs(p->JetInvMasses_eta1stJ)<5){
		cout<<( p->JetInvMasses_flagCenteralJetVeto)<<endl;
		cout<<( p->JetInvMasses_deltaEtaJ1J2)<<endl;
		cout<<( p->JetInvMasses_mJ1J2)<<endl;
		cout<<( p->JetInvMasses_eta1stJ)<<endl;
	//	getchar();
		counter ++;
		}
		}
	cout<<counter<<endl;

	}


	int main(int argc, char **argv) {
		if (argc < 2)
		{
			printf("\nUsage: %s *.root\n\n",argv[0]);
			exit(0);
		}

		int const index = argc - 1;

		TString path[index];

		for ( int j = 0; j < argc-1; j++)
		{
			path[j] = (argv[j+1]);
		}

		TChain *chain = new TChain("OutputTree_");

		for (int k = 0; k < argc-1 ; k++)
		{
			chain -> Add(path[k]); // Add the root file to the TChain chain
		}
	/*	TString inputfile(argv[1]);
		TFile f(inputfile,"read");
		TTree *t = (TTree *) f.Get("OutputTree_");

		if(!t){
			cout<<"No OutputTree_ tree found in "<<inputfile<<endl;
			exit(1);
		}
*/
		p = new OutputTree_(chain);

		make_histgram();

		delete chain;
		delete p;


	}