#include "Riostream.h"
#include "TString.h"
#include "TNtuple.h"
#include "TSystem.h"

void tidfit() {
	
	 ifstream in;
   in.open("Run_13-19.txt");
   
   
   Float_t Time,Vi_Brick1,Io_Brick1,Vo_Brick1,Vi_Brick2,Io_Brick2,Vo_Brick2,Vi_Brick3,Io_Brick3,Vo_Brick3,Vi_Brick4,Io_Brick4,Vo_Brick4,Vi_Brick1err,Vi_Brick2err,Vi_Brick3err,Vi_Brick4err,Io_Brick1err,Io_Brick2err,Io_Brick3err,Io_Brick4err,Vo_Brick1err,Vo_Brick2err,Vo_Brick3err,Vo_Brick4err;
   Int_t n1 = 12;
   auto f = TFile::Open("Run_13-19.root","RECREATE");
   
    TNtuple ntuple("ntuple","data from ascii file","Time:Vi_Brick1:Vi_Brick1err:Io_Brick1:Io_Brick1err:Vo_Brick1:Vo_Brick1err:Vi_Brick2:Vi_Brick2err:Io_Brick2:Io_Brick2err:Vo_Brick2:Vo_Brick2err:Vi_Brick3:Vi_Brick3err:Io_Brick3: Io_Brick3err:Vo_Brick3:Vo_Brick3err:Vi_Brick4:Vi_Brick4err:Io_Brick4:Io_Brick4err:Vo_Brick4:Vo_Brick4err");


    while (1) {
           
            in >>  Time >> Vi_Brick1 >> Vi_Brick1err >> Io_Brick1 >> Io_Brick1err >> Vo_Brick1 >> Vo_Brick1err >> Vi_Brick2 >> Vi_Brick2err >> Io_Brick2 >>Io_Brick2err >> Vo_Brick2 >>Vo_Brick2err >> Vi_Brick3 >>Vi_Brick3err >>  Io_Brick3 >>  Io_Brick3err >>  Vo_Brick3 >>  Vo_Brick3err >> Vi_Brick4 >> Vi_Brick4err >> Io_Brick4 >> Io_Brick4err >> Vo_Brick4 >> Vo_Brick4err;
          if(!in.good()) break;
	   
	}
   
   
  Double_t time1[] = {Time};
  Double_t Vi_brk[] = {Vi_Brick1,Vi_Brick2,Vi_Brick3,Vi_Brick4};
  Double_t Io_brk[] ={Io_Brick1,Io_Brick2,Io_Brick3,Io_Brick4};
  Double_t Vo_brk[] ={Vo_Brick1,Vo_Brick2,Vo_Brick3,Vo_Brick4};
  
  Double_t Vi_brkerr [] ={Vi_Brick1err,Vi_Brick2err,Vi_Brick3err,Vi_Brick4err};
  Double_t Io_brkerr [] = {Io_Brick1err,Io_Brick2err,Io_Brick3err,Io_Brick4err};
  Double_t Vo_brkerr [] = {Vo_Brick1err,Vo_Brick2err,Vo_Brick3err,Vo_Brick4err};
  
  
  TGraphErrors *g1 = new TGraphErrors(n1,time1,Vi_brk,0,Vi_brkerr);
   
   TGraphErrors *g2 = new TGraphErrors(n1,time1,Io_brk,0,Io_brkerr);
     
   TGraphErrors *g3 = new TGraphErrors(n1,time1,Vo_brk,0,Vo_brkerr);
 
   
   TCanvas *c1 = new TCanvas("c1","tidfit",700,500);
         g1->SetMarkerStyle(21);
         g1->SetMarkerColor(2); 
		 g1->Draw("ap");
		 
	 TCanvas *c2 = new TCanvas("c2","tidfit",700,500);
         g2->SetMarkerStyle(21);
         g2->SetMarkerColor(2); 
		 g2->Draw("ap");

 TCanvas *c3 = new TCanvas("c3","tidfit",700,500);
         g3->SetMarkerStyle(21);
         g3->SetMarkerColor(2); 
		 g3->Draw("ap");		 
    
   
       in.close();
}