// read data from BLMLHC_M Class post mortem data file
// Authors: M. Stockner, D. Kramer
// modified:
//		16.11.2006

#include <stdio.h>
#include <stdlib.h>
#include <Riostream.h>
#include <TString.h>
#include <TROOT.h>
#include <TStyle.h>
#include <TSystem.h>
#include <TCanvas.h>
#include <TGraph.h>
#include <TMultiGraph.h>
#include <TLegend.h>
#include <TApplication.h>
#include <TObjArray.h>
#include <TObjString.h>

int main(int argc, char **argv)
{
  gROOT->SetStyle("Plain");
  gStyle->SetPalette(1);
	gROOT->ForceStyle();

  TString filename;// = "PMData-20061115080058.txt";
	
  if(argc>1) filename = argv[1];
  else {cout << "Filename ???" << endl; return -1;}

	// cout << filename <<endl;
  // data separator
  char splitter[] = ",";		//possilbe = " ,.";
  // read ascii file into 2 dim matrix
  int x,y,i;
  TString next;
  TString aLine;
  ifstream in;
  in.open(filename);
  if(!in.good()) {cout << "file error" <<endl; in.close(); return -1;}
  x=0; y=0;
  double myarray[20][48000];
	TObjArray* StringArray;

  TApplication myApp("App", &argc, argv);
	
	while(aLine.ReadLine(in))
	{
    if(aLine[0]=='P') continue;
    if(aLine[0]=='<') continue;		
		x=0;
		StringArray = aLine.Tokenize(splitter);
		for(i=0; i<StringArray->GetEntries(); i++){
			next = ((TObjString*)StringArray->At(i))->GetString();
			myarray[x][y] = next.Atof();
			x++;
		}
		y++;
//  		if(y>40000) break;    
	}	
	printf ("x: %d y: %d \n",x,y);
  // end of read ascii file into 2 dim matrix
  // results:
  // x ... ROWS, int
  // y ... LINES, int
  // myarray ... all data from file, double

//   for(int a=0;a<x;a++) {
//     for(int u=0;u<y;u++) {
//       cout << myarray[a][u] << " ";
//     }
// 		cout << endl;
// 	}

	TCanvas *myc = new TCanvas("myc","Simulation Response",1);
	myc->SetFillColor(0);
	myc->SetLeftMargin(0.11);  
//   myc->SetLogx();
//   myc->SetLogy();
//   myc->SetTickx();
//   myc->SetTicky();
  
  TMultiGraph *mg1 = new TMultiGraph();

  TLegend *leg1 = new TLegend(0.903736,0.790254,0.989943,0.902542);
  leg1->SetFillColor(0);
	TString grname[] = {"BLMS","BLMI"};
	
	int v = 0;
  double tmparray[y];
  for(int i=0;i<y;i++) tmparray[i] = (double)i*40./1e3;
	
  for(int a=0;a<x;a++) {
		if(a==1 || a==9){
			TGraph *gr = new TGraph(y, tmparray, myarray[a]);
			gr->SetLineColor(a);
			gr->SetMarkerColor(a); 
			// gr->SetMarkerStyle(markernb[a]);
			mg1->Add(gr);
			leg1->AddEntry(gr,grname[v++],"lp");
		}
  }

//   mg1->SetMinimum(1e-6);
//   mg1->SetMaximum(2e-3);
  mg1->Draw("ALP");
	mg1->SetTitle(filename);
	mg1->GetXaxis()->SetRangeUser(0.,tmparray[y-1]);	

  mg1->GetXaxis()->SetTitle("time [ms]");
  mg1->GetXaxis()->SetTitleOffset(1.1);
  mg1->GetYaxis()->SetTitle("counts + ADC");
  mg1->GetYaxis()->SetTitleOffset(1.5);
	leg1->Draw();

myApp.Run();

return 0;
}