// CPP includes
#include <cstdlib>
#include <cstdio>
#include <iostream>
#include <cmath>
#include <vector>
#include <string>
#include <iostream>
#include <sstream>
#include <fstream>

// C includes
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>

// ROOT includes
#include <TSystem.h>
#include <TROOT.h>
#include <TApplication.h>
#include <TRint.h>
#include <TDirectory.h>
#include <TFile.h>
#include <TChain.h>
#include <TTree.h>
#include <TNtuple.h>
#include <TBranch.h>
#include <TH1.h>
#include <TAxis.h>
#include <TCanvas.h>
#include <TColor.h>
#include <TPad.h>
#include <TStyle.h>
#include <TKey.h>
#include <TLegend.h>
#include <TText.h>
#include <TPaveText.h>
#include <TLine.h>
#include <TMarker.h>
#include <TLatex.h>
#include <TMath.h>
#include <TH1.h>
#include <TH2.h>
#include <TH3.h>
#include <TF1.h>
#include <TF2.h>
#include <TString.h>
#include "TObjString.h"
#include "TRandom3.h"
#include "TObjArray.h"
#include "TFrame.h"


using namespace std;

#ifdef __MAKECINT__
#pragma link C++ class vector<float>+;
#endif
void run(){
    
   TFile *f = TFile::Open("Results.root","READ");
   f->ls();
   TTree *t= (TTree*)f->Get("tvec");

  
  vector<vector<double> > history_ps_x_all_part;
  vector<vector<double> > history_ps_y_all_part;
  vector<double> *x = 0;
  vector<double> *y = 0;
  
  TBranch *bx = 0;
  TBranch *by = 0;
  t->SetBranchAddress("x",&x,&bx);
  t->SetBranchAddress("y",&y,&by);
   


   for (Int_t i = 0; i < 1; i++) {
      
      Long64_t tentry = t->LoadTree(i);
      bx->GetEntry(tentry);
      history_ps_x_all_part.push_back(*x);
      
      by->GetEntry(tentry);
      history_ps_y_all_part.push_back(*y);

       
      
}

      for (Int_t i = 0; i < 1; i++) {
          for (UInt_t j = 0; j <history_ps_x_all_part[i].size() ; ++j) {
                 cout << history_ps_x_all_part[i][j]<< " " <<   history_ps_x_all_part[i][j]<<endl;
                 
          }}

   
   // Since we passed the address of a local variable we need
   // to remove it.
   t->ResetBranchAddresses();
    
   
}   
  
void run2(){
    
   TFile *f = TFile::Open("Results.root","READ");
   f->ls();
   TTree *t= (TTree*)f->Get("tvec");

   // Prepare matrixes to store the data read from the file containing the rsults.
   // The first matrix will contains, as first entry, all the x_positions of first particle followed,
   // as second entry, all the x_positions of the second particle followed and so on;
   // The second matrix will contains, as first entry,all the y_positions of first particle followed,
   // as second entry, all the y_positions of the second particle followed and so on;
   // The same for all othe matrixes down here. 

   /*vector<vector<double> > history_ps_x_all_part;
   vector<double> pos_x;
   vector<double> *x;
   double sigma_x=2e-3;
   double distance = 0.05;
   
   int nEntries=t->GetEntries();

   
   cout << t->GetEntries() << endl;
 

   for(int i=0;i<1;i++){
       
      // Long64_t tentry = t->LoadTree(i);
      // TBranch *bvpx = 0;
       t->SetBranchAddress("x",&x);
       //bvpx->GetEntry(tentry);
       //t->GetEntry(i);
       cout << x->size()<<endl;
       for (UInt_t j = 0; j < x->size(); ++j) 
         {
           cout <<j << " " <<x->at(j)<<endl;
           pos_x.push_back(x->at(j));
          }
       history_ps_x_all_part.push_back(pos_x);
       t->ResetBranchAddresses();         
                           } */
  vector<vector<double> > history_ps_x_all_part;
  vector<vector<double> > history_ps_y_all_part;
  vector<double> *x = 0;
  vector<double> *y = 0;
  
  TBranch *bx = 0;
  TBranch *by = 0;
  t->SetBranchAddress("x",&x,&bx);
  t->SetBranchAddress("y",&y,&by);
   


   for (Int_t i = 0; i < 2; i++) {
      
      Long64_t tentry = t->LoadTree(i);
      bx->GetEntry(tentry);
      history_ps_x_all_part.push_back(*x);
 
      //by->GetEntry(tentry);

       
      
}

      for (Int_t i = 0; i < 2; i++) {
         // for (UInt_t j = 0; j <history_ps_x_all_part[i].size() ; ++j) {
                 cout << history_ps_x_all_part[i][0]<<endl;
                 
          }//}

   
   // Since we passed the address of a local variable we need
   // to remove it.
   t->ResetBranchAddresses();
    
   
}