//
//  Tree_reading.cpp
//  
//
// 
//


#include <algorithm>
#include <TRandom3.h>
void Tree_reading(){
    TRandom3 ran3;
    TCanvas *c1= new TCanvas();
    TFile *Input =new TFile("xBQ2binning.root","read");
    
    double xBj,QSq;
    int ixB,jQ2;
    double zeta_min=0.0, zeta_max=1.0;
    double zeta;
    double W2min=4.0;
    double MpSq=0.938*0.938;
    
     TTree * tree= (TTree*) Input->Get("xBQ2binning_tree");
  
    tree->SetBranchAddress("xBj", &xBj);
    tree->SetBranchAddress("QSq", &QSq);
    tree->SetBranchAddress("ixB", &ixB);
    tree->SetBranchAddress("jQ2", &jQ2);
    
    
    
    int entries= tree->GetEntries();
    
//    cout << entries<< endl;
    double Q2Min[jQ2],Q2Max[jQ2],Q2;
    
    
    for(int i=0;i< entries;i++)
    {
        
        tree->GetEntry(i);
  
       
    
    for(int k=0; k< ixB; k++){
        for(int j=0; j < jQ2; j++){
           Q2Min[j]= tree->GetMinimum("QSq");
          //  Q2Max[j][k]= tree->GetMaximum("QSq");
            
            //Q2Max[j]= tree->GetMinimum("QSq");
         //   printf("Q2Min[j]=%f \n", Q2Min[j]);
            //Q2=(Q2Max[j]-Q2Min[j])*(ran3.Uniform())+ Q2Min[j];
            
            
        }
        
    }
    
    }
    

   
    
   
   
    
    Input->Close();
    
}