/* 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"

int main(int argc, char **argv) {  
  
  int Slot = 0;
  int DISPLAY = 1;

  int version = 1;

  char files[4][128];

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

      if (!strcmp(argv[i],"h")){
	cout<<"-P string  print plot to file \"string\"   "<<endl;
	cout<<"-V #     data sample number"<<endl;
	return 0;
      }
      
    }
  }


  sprintf(files[0],"tdcdata_ROCFDC1_%d_ALLslots.dat",version);
  sprintf(files[1],"tdcdata_ROCFDC14_%d_ALLslots.dat",version);
  sprintf(files[2],"tdcdata_ROCFDC4_%d_ALLslots.dat",version);
  sprintf(files[3],"tdcdata_ROCFDC13_%d_ALLslots.dat",version);


  int Plane[4][12] = { 1,2,2,3,3,4,1,2,2,3,3,4,
		       1,4,5,5,6,6,1,4,5,5,6,6,
		       1,2,2,3,3,4,1,2,2,3,3,4,
		       1,4,5,5,6,6,1,4,5,5,6,6};

  int Offsets[4][12] = { 0,0,2,0,2,2,0,0,2,0,2,2,
			 2,0,0,2,0,2,2,0,0,2,0,2,
			 0,0,2,0,2,2,0,0,2,0,2,2,
			 2,0,0,2,0,2,2,0,0,2,0,2};

  TApplication theApp("App", &argc, argv);
  
  char *val = getenv("USER");
  char username[128];
  strcpy(username,val);
  char DataDir[128];
  sprintf(DataDir,"/gluonraid1/Users/%s/fdcdata",username);
  char fnam[128];

  char str1[128],str2[128];
  TH1F *WirePlanes[4][6];

  for (int i=0;i<4;i++) {
    for (int k=0;k<6;k++) {
      sprintf(str1,"WirePlane%d%d",i+1,k+1);
      sprintf(str2,"FDC %d TDCDATA WirePlane%d",i+1,k+1);    
      WirePlanes[i][k] = new TH1F(str1, str2, 98,-0.5,97.5);
    }
  }

  for (int n=0;n<4;n++){ 

    sprintf(fnam,"%s/%s",DataDir,files[n]);  
    ifstream INF;
    INF.open(fnam, ios::in);
  
    if (!INF) {
      cout<<"file not found: "<<fnam<<endl;
      return 0;
    }


    int SLOT = -1;
    int Chan = -1;
    int EVENT = 0;
    int time;
    
    char lin[5000];
    while (!INF.eof()){
      
      INF.getline(lin,5000);
      string s = lin;
      istringstream iss (s);
      
      int NewSLOT;
      int NewChan;
      int NewEVENT;
      int time;
      iss>>NewEVENT>>NewSLOT>>NewChan>>time;
      
      if (NewChan != Chan){
	
	EVENT = NewEVENT;
	SLOT = NewSLOT;
	Chan = NewChan;
	int Package = 0;

	if ((n<2) && (SLOT>5)){
	  Package = 1;
	} else if ((n>1) && (SLOT<6)){
	  Package = 2;
	} else if ((n>1) && (SLOT>5)){
	  Package = 3;
	}

	int c = Chan + 1 + Offsets[n][SLOT]*24;

	if (c<25){
	  c += 24;
	} else if (c<49){
	  c -= 24;
	} else if (c<73){
	  c += 24;
	} else if (c<97){
	  c -= 24;
	}

	//cout<<Plane<<" "<<c<<"     "<<EVENT<<" "<<SLOT<<" "<<Chan<<endl;
	//cout<<SLOT<<"  "<<Chan<<"     "<<Plane<<"  "<<c<<endl;

	//if (EVENT>20000000)
	  WirePlanes[Package][Plane[n][SLOT]-1]->Fill((float)c);
	
      }
      
    }
    
    INF.close();
    
  }

  TFile *fout;
  sprintf(fnam,"%s/tdcs_ALLPackages_NEW%d.root",DataDir,version);
  fout = new TFile(fnam,"RECREATE");
  

  for (int n=0;n<4;n++) {
    for (int k=0;k<6;k++) {
      WirePlanes[n][k]->Write();
    }
  }
  fout->Close();

  return 0;
  
}