#include "TComplex.h"
#include "TCanvas.h"
#include "TPad.h"
#include "TH1.h"
#include "TF1.h"
#include "TGraph.h"
#include "TStyle.h"
#include "TVirtualFitter.h"
#include "TDatime.h"
#include "TTree.h"
#include "TFile.h"
#include "Riostream.h"


using namespace std;

void Readdat(){  

//---------------------------------------------------------------------
//---------------------Readdata from .txt:   
//--------------------------------------------------------------------- 
  ifstream in_x,f;
 
  f.open("Envrconv_05112006_SAVE2.txt");
  string line;
  Int_t nlines=0;
  while(f.eof()==false && getline(f,line,'\n'))
  {
  //   getline(f,line,'\n');
  //   cout << line << endl;
    nlines++;
  }
  f.close();

  cout << "Number of lines: "<< nlines << endl;

  Int_t ELMB,channel,year,month,day,hours,min;
  Double_t sec,totalsec,val;
  char t;
  Int_t minsec=0;
  Int_t hourssec=0;
  Int_t minhoursec=0;
  Int_t lines=0;

  in_x.open("Envrconv_05112006_SAVE2.txt"); 
  TArrayD *ELMBarr = new TArrayD(nlines);
  TArrayD *Chan = new TArrayD(nlines);    
  TArrayD *time = new TArrayD(nlines);
  TArrayD *temp = new TArrayD(nlines);     

  while (in_x.eof()==false &&   in_x >> ELMB >> channel >> year >> t >> month >> t >> day >> hours >> t >> min >> t >>sec
    >>val) 
  {
     if (!in_x.good()) cout<< "Can't read file" << endl;

     hourssec = (hours* 60.0* 60.0); 
     minsec = (min* 60.0);
     minhoursec = (hourssec + minsec);
     totalsec = (minhoursec + sec);

     ELMBarr->AddAt(ELMB,lines);
     Chan->AddAt(channel,lines);
     time->AddAt(totalsec,lines);
     temp->AddAt(val,lines);
     lines++;      
  } 
 
 in_x.close();

//---------------------------------------------------------------------
//---------------------Read Mappingfile:   
//--------------------------------------------------------------------- 

  ifstream m,in_m;

  m.open("Mapping.txt");
  int mlines = 0;
  string ELMBmaps;
  string ELMBChs;
  string Disk;
  string Type;
  string Side;
  string Quadrant;
  string Num;
  string identifier;

  vector<string> Identivec;

  while(m.eof()==false &&  m >> ELMBmaps >> ELMBChs >> Disk >> Type >> Side >> Quadrant >> Num) 
  {
    identifier ="ELMB_"+ELMBmaps + "_Ch" + ELMBChs + "_D" + Disk +"_" + Type + "_" +Side
    +"_" + Quadrant + "_Num:" + Num; 

    Identivec.push_back(identifier);

    mlines++;
  }
  m.close();

  cout << "Number of mlines: "<< mlines << endl;
 /*
  for(int k=0;k<Identivec.size();k++) 
  { 
    cout<< Identivec[k] <<endl; 
  } 
 */
  in_m.open("Mapping.txt");
  int count = 0;
  Int_t ELMBmap;
  Int_t ELMBCh;

  TArrayD *ELMBmaparr = new TArrayD(mlines);
  TArrayD *ELMBCha = new TArrayD(mlines);     


  while(in_m.eof()== false &&  in_m >> ELMBmap >> ELMBCh >> Disk >> Type >> Side >> Quadrant >> Num) 
  {

    ELMBmaparr->AddAt(ELMBmap,count);
    ELMBCha->AddAt(ELMBCh,count);

    count++;
  }
  in_m.close();

//---------------------------------------------------------------------
//---------------------stopppoints to divide data per channel
//----------------------------------------------------------------------

  TArrayD *stoparr = new TArrayD(mlines);
  TArrayD *Chanstop = new TArrayD(mlines);
  TArrayD *ELMBstop = new TArrayD(mlines);
  Int_t p=1;
  stoparr->At(0)==0;

  for(int i=0; i<lines; i++){
    if(Chan->At(i)!= Chan->At(i+1)){
      stoparr -> AddAt(i,p);
      Chanstop -> AddAt(Chan->At(i),p);
      ELMBstop -> AddAt(ELMBarr->At(i),p);
      //   cout << " p: " << p << " stoparr: " << stoparr->At(p) <<endl;
      //   cout << " p: " << p << " Chanstop: " << Chanstop->At(p) << " ELMBstop: " << ELMBstop->At(p) <<endl;
      p++;  
    }

  }
//---------------------------------------------------------------------
//---------------------Mapping   
//----------------------------------------------------------------------


  TArrayD *idselectarr = new TArrayD(mlines);
  Int_t id=0;

  for(int i=0;i<p;i++){
    for(int j=0;j<mlines;j++){
      if(ELMBstop->At(i)== ELMBmaparr->At(j) && Chanstop->At(i)==ELMBCha->At(j)){
	idselectarr ->AddAt(j,id);
	id++;
      }   
    }
  }

   
  //for(int s=0; s<id; s++){
  // cout<< idselectarr->At(s) <<endl;    
  // cout<< Identivec[idselectarr->At(s)]<<endl;  
  //}


//---------------------------------------------------------------------
//---------------------create Tree and fill it:   
//--------------------------------------------------------------------- 

  Double_t tELMB, tchannel,ttsec,tval;
  Int_t tlines=0;

  TFile g("Envrconv_05112006.root","recreate");
  TTree *t1 = new TTree("Envronetree","HiTree");
  t1->Branch("tlines",&tlines, "tlines/I");
  t1->Branch("tELMB",&tELMB,"tELMB/D");
  t1->Branch("tchannel",&tchannel,"tchannel/D"); 
  t1->Branch("ttsec",&ttsec,"ttsec/D");
  t1->Branch("tval",&tval,"tval/D");

  for(Int_t i=0; i<nlines; i++){
    tELMB = ELMBarr->At(i);
    tchannel = Chan->At(i);
    ttsec = time->At(i);
    tval = temp->At(i); 
    tlines = i; 

    t1->Fill();
  }
 
// t1->Print();
  g.cd();
  t1->Write();
 
//-----------------------Tree for seperating channels--------------------------//   

  Int_t tstop,nstops;

  TFile fg("treefg.root","recreate");
  TTree *tfg = new TTree("tfg","Idtree");
  tfg->Branch("tstop",&tstop,"tstop/I");
  tfg->Branch("nstops",&nstops,"nstops/I");

  for(Int_t i=0; i<p; i++){
    tstop = stoparr->At(i);
    nstops = i;

    tfg->Fill();
  } 
  fg.cd();
  tfg->Write();

//-----------------------Tree for identifier--------------------------//   

  vector<string> *tidentivec = new vector<string>;

  TFile fi("identi.root","recreate");
  TTree *Treeidenti = new TTree("Treeidenti","Vectortree");
  Treeidenti->Branch("tidentivec","vector<string>",&tidentivec);
  
  for(Int_t i=0; i<id; i++){
    tidentivec->push_back(Identivec[idselectarr->At(i)]);
     
    Treeidenti->Fill();    
  } 
  fi.cd();
  Treeidenti->Write();
}