#include <vector> 
#include "TTree.h"
#include "TFile.h"
#include <TROOT.h>
#include <iostream>
#ifdef __MAKECINT__
#pragma link C++ class std::vector < std::vector<int> >+;   
#endif 
using namespace std;
void add_var(){
   
  Int_t tau_n;
  
  vector<float> * tau_etEflow= 0;
  vector<float> *tau_jet_Et=0;
  vector<vector<int> > *tau_track_nTRTHits = 0;
  vector<vector<int> > *tau_track_nTRTHighTHits = 0;
  vector<vector<int> > *tau_track_nTRTOutliers = 0;
  vector<vector<int> > *tau_track_nTRTHighTOutliers = 0;
  TBranch        *b_tau_track_nTRTHighTHits;   //!
  TBranch        *b_tau_track_nTRTOutliers;   //!
  TBranch        *b_tau_track_nTRTHighTOutliers;   //!
  TBranch        *b_tau_n; 
  TBranch        *b_tau_track_nTRTHits;

  TBranch  *b_tau_etEflow;
  TBranch  *b_tau_jet_Et; 
  vector<Float_t> *TRThitsRatio = new vector<Float_t>;
  
  vector<Float_t> *etEflowOverjetEt  = new vector<Float_t>;

  TFile f("/data0/asircar/Slims/NewSlims/j4.root", "READ");
  TTree *t3 = (TTree*)f.Get("tauPerf");
  t3->SetBranchAddress("tau_track_nTRTHits", &tau_track_nTRTHits, &b_tau_track_nTRTHits);
  t3->SetBranchAddress("tau_track_nTRTHighTHits", &tau_track_nTRTHighTHits, &b_tau_track_nTRTHighTHits);
  t3->SetBranchAddress("tau_track_nTRTOutliers", &tau_track_nTRTOutliers, &b_tau_track_nTRTOutliers);
  t3->SetBranchAddress("tau_track_nTRTHighTOutliers", &tau_track_nTRTHighTOutliers, &b_tau_track_nTRTHighTOutliers);
  t3->SetBranchAddress("tau_n", &tau_n, &b_tau_n);
  
  t3->SetBranchAddress("tau_etEflow", &tau_etEflow, &b_tau_etEflow);
  t3->SetBranchAddress("tau_jet_Et",&tau_jet_Et, &b_tau_jet_Et); 
  
  TFile fnew("j4.new.root","RECREATE");
  TTree* t4 = t3->CloneTree(0);
  
  t4->Branch("tau_TRThitsRatio",&TRThitsRatio);      
  
  t4->Branch("tau_etEflowOverjetEt",&etEflowOverjetEt);

  Long64_t nentries = t3->GetEntries();
  //cout<<"entries"<< nentries<<endl;
    for (Long64_t i = 0; i < nentries; i++){
    //  for (Long64_t i = 0; i < 100; i++){
    TRThitsRatio->clear();
   
    etEflowOverjetEt->clear();
   
     t3->GetEntry(i);
      //cout<< "num tau"<<tau_n<<endl;

    for(int j = 0; j!= tau_n; j++){
      int    highthresh=0 ;
      int    lowthresh=0 ; 
      float  ratio =0;

      float  rationew=0;   
      
      vector<int> nTRTHits = tau_track_nTRTHits->at(j);
      vector<int> nTRTOutliers = tau_track_nTRTOutliers->at(j);
      vector<int> nTRTHighThreshHits = tau_track_nTRTHighTHits->at(j);
      vector<int> nTRTHighThreshOutliers = tau_track_nTRTHighTOutliers->at(j);

      //cout <<"len hits"<<nTRTHits.size()<<endl;
      
    float etflow = tau_etEflow->at(j);
    cout <<"et"<< etflow<<endl; 
   
   float jetet =  tau_jet_Et->at(j);      
   cout <<"jetEt"<< jetet<<endl; 
   if (etflow!=-1111 && jetet!=0){
   rationew = etflow/jetet;
   cout <<"ratio"<<rationew<<endl;
   // if(fabs(rationew)>1000000){
   //cout<<":::::::::::::::::::::::::::::::::::::::::::::::::: ratio fail"<<rationew<<endl;
   // }
    etEflowOverjetEt->push_back(rationew);  
   }

      int k = 0;
      if(nTRTHits.size()){


	lowthresh = nTRTHits.at(k) + nTRTOutliers.at(k);

	//cout<< "trthits"<< nTRTHits.at(k)<< endl;
	//cout<< "trtouthits"<<nTRTOutliers.at(k)<< endl;


	highthresh = nTRTHighThreshHits.at(k)+nTRTHighThreshOutliers.at(k);

	//cout<<"hihits"<<nTRTHighThreshHits.at(k)<<endl;
	//cout<<"hiout"<<nTRTHighThreshOutliers.at(k)<<endl;   

         

	//cout<<"lw"<<lowthresh <<endl;


	//cout<<"hi"<<highthresh <<endl;
	if( lowthresh > 0 ) ratio = (float(highthresh)/float(lowthresh));
	else ratio = -1111;
	//cout<<"ratio"<<ratio<<endl;
	
	TRThitsRatio->push_back(ratio);

      }
      else {TRThitsRatio->push_back(-1111);}
     

  
  }

    t4->Fill();
  }
 

  fnew.cd();
  fnew.Write();
  fnew.Close();
}