/* include files */
#include <iostream>
#include <sstream>   
#include <iomanip>
#include <stdio.h>
#include <stdlib.h>
#include <fstream>
#include <string>
using namespace std;

// root include files
#include "TApplication.h"  // needed to display canvas
#include "TSystem.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TH1.h"
#include "TH1F.h" 
#include "TH2F.h"
#include "TH3F.h"
#include "TF1.h"
#include "TMath.h"
#include "TGraph.h"
#include "TGraphErrors.h"
#include "FDCchannels.h"

int main(int argc, char **argv) {  
  
  int Slot = 0;
  int DISPLAY = 1;
  char CrateName[128] = "One";
  char filename[128] = "samples.pdf";
  int version = 1;
  int doprint=0;
  int TheRoc=0;
  char InputFiles[10][128];
  int ROCNum[10] = {2,3,5,6,7,8,9,10,11,12};
  int ADCSlotNum[15] = {3,4,5,6,7,8,9,10,13,14,15,16,17,18,19};
  int SCANNER=1;
  char detector[128] = "FDC";

  if (argc>1){
    
    for (int i=1;i<argc;i++){
      
      if (!strcmp(argv[i],"-ND")){
	DISPLAY = 0;
	SCANNER = 0;
      }

      if (!strcmp(argv[i],"-P")){
	sprintf(filename,"%s",argv[i+1]);
	i++;
	doprint = 1;
      }

      if (!strcmp(argv[i],"-V")){
        version = atoi(argv[i+1]);
        i++;
      }

      if (!strcmp(argv[i],"-R")){
        TheRoc = atoi(argv[i+1]);
        i++;
      }

      if (!strcmp(argv[i],"-S")){
        SCANNER = 0;
      }

      if (!strcmp(argv[i],"-D")){
	sprintf(detector,"%s",argv[i+1]);
	i++;
      }

      if (!strcmp(argv[i],"h")){
	cout<<"-ND        do not display individual events"<<endl;
	cout<<"-P string  print plot to file \"string\"   "<<endl;
	cout<<"-V #     data sample number"<<endl;
	cout<<"-R #     scan only ROC #"<<endl;
	return 0;
      }
      
    }
  }

  TApplication theApp("App", &argc, argv);

  ArrayInit();

  for (int n=0;n<10;n++){
    sprintf(InputFiles[n],"adcdata_ROC%s%d_%d_ALLslots.dat",detector,ROCNum[n],version);
  }
  
  char DataDir[128] = "/gluonraid1/Users/zihlmann/FDCDATA";
  char fnam[128];
  ifstream INF;
  int MAXAdcs = 15;
  TH1F *mult[10][15];
  TH1F *Strips[4][6][2];
  char t1[128];
  char t2[128];
  for (int m=0;m<10;m++) {
    for (int k=0;k<MAXAdcs;k++) {
      int s = ADCSlotNum[k];
      sprintf(t1,"mult%d_%d",m,s);
      sprintf(t2,"Hit multiplicity ROC %d SLOT %d",ROCNum[m],s);
      mult[m][k] = new TH1F(t1,t2,74,-1.5,72.5);  
    }
  }

  for (int m=0;m<4;m++) {
    for (int i=0;i<6;i++) {
      for (int j=0;j<2;j++) {
	sprintf(t1,"StripsFDC%dPlane%dU",m+1,i+1);
	sprintf(t2,"Hit multiplicity FDC %d Plane %d U",m+1,i+1);
	if (j){
	  sprintf(t1,"StripsFDC%dPlane%dD",m+1,i+1);
	  sprintf(t2,"Hit multiplicity FDC %d Plane %d D",m+1,i+1);
	}
	Strips[m][i][j] = new TH1F(t1,t2,24*9+2,-0.5,24.*9.+1.5);	
      }
    }
  }
  


  TCanvas *myc1;
  if (DISPLAY) {
    char str1[128];
    sprintf(str1,"ADC125 Samples");
    myc1 = new TCanvas("myc1", str1, 800, 400);
    myc1->SetFillColor(42);
    gPad->SetGrid();
  }
  
  double array[5000];
  double xaxis[5000];
  int loli = 0;
  int hili = 10;
  if (TheRoc){
    loli = TheRoc-2;
    if (TheRoc>3){
      loli = TheRoc-3;
    }
    hili = loli+1;
  }
  for (int n=loli;n<hili;n++) {
    
    char lin[5000];
    int NSlot = -1;
    int NCon = -1;
    int EventNew = -1;
    
    sprintf(fnam,"%s/%s",DataDir,InputFiles[n]);

    cout<<"open file "<<fnam<<endl;
    INF.open(fnam, ios::in);
    
    if (!INF) {
      cout<<"file not found: "<<fnam<<endl;
    } else {
    
      while (!INF.eof()){

	INF.getline(lin,5000);
	string s = lin;
	istringstream iss (s);
	
	int SLOT;
	int chan;
	int EVENT;
	iss>>EVENT>>SLOT>>chan;
	
	//if ((EVENT != EventNew) || (NSlot != SLOT)){
	if (1){
	  mult[n][SLOT]->Fill(chan);
	  EventNew = EVENT;
	  NSlot = SLOT;
	  int con = (chan)/24;
	  int FDC = Package[n][SLOT][con];
	  int P =  Plane[n][SLOT][con];
	  int C =  Connector[n][SLOT][con];
	  int L = Layer[n][SLOT][con];
	  float channel = (C-1)*24. + (23-(chan - con*24)) + 1; // Reverse order!=> 23-c
	  Strips[FDC-1][P-1][L]->Fill(channel);
	  
	}

	if ((DISPLAY) && (SLOT>-1)){
	  int Max = 0;
	  while (iss>>array[Max]){
	    xaxis[Max] = Max;
	    Max++;
	  }

	  //cout<<EVENT<<" "<<SLOT<<" "<<chan<<"  :"<<Max<<endl;
	  
	  char tit[128];
	  sprintf(tit,"Crate ROC %d : ADC SLOT %d Channel %d",
		  ROCNum[n], ADCSlotNum[SLOT],chan);
	  //cout<<Max<<"  "<<tit<<endl;
	  //if (gr){
	  //        delete gr;
	  //}
	  
	  if (Max>0){
	    TGraph *gr =  new TGraph(Max,xaxis,array);
	    gr->SetTitle(tit);
	    gr->SetMarkerColor(4);
	    gr->SetMarkerStyle(21);
	    gr->GetXaxis()->SetTitle("sample number");
	    gr->GetYaxis()->SetTitle("ADC counts");  
	    gr->Draw("ALP");
	    myc1->Update();
	    gPad->SetGrid();
	    gPad->Modified();
	    
	    if (doprint){
	      myc1->SaveAs(filename);
	    }
	    
	    if (SCANNER){
	      cout<<"continue?"<<EVENT<<" :";
	      char answ[128];
	      cin>>answ;
	    }
	  }
	}
	
	
      }
      
      INF.close();
    }
    
  }
  
  TFile *fout;
  sprintf(fnam,"%s/adcs_hitmultiplicities_%d_ALLROCS.root",DataDir,version);
  if (TheRoc){
    sprintf(fnam,"%s/adcs_hitmultiplicities_%d_ROC%d.root",DataDir,version,TheRoc);
  }
  fout = new TFile(fnam,"RECREATE");


  for (int m=0;m<10;m++){
    for (int k=0;k<MAXAdcs;k++){
      mult[m][k]->Write();
    }
  }
  for (int m=0;m<4;m++) {
    for (int i=0;i<6;i++) {
      for (int j=0;j<2;j++) {
	Strips[m][i][j]->Write();	
      }
    }
  }

  fout->Close();

  return 0;
  
}