/* 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) {  
  
  char CrateName[128] = "ROCFDC";
  int version = 1;
  int ROC = 2;
  int SLOT = 3;
  int CHANNEL = 1;
  int DISPLAY = 1;

  char InputFile[128];
  int ROCNum[10] = {2,3,5,6,7,8,9,10,11,12};

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

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

      if (!strcmp(argv[i],"-V")){
        version = atoi(argv[i+1]);
        i++;
      }
      if (!strcmp(argv[i],"-R")){
        ROC = atoi(argv[i+1]);
        i++;
      }
      if (!strcmp(argv[i],"-S")){
        SLOT = atoi(argv[i+1]);
        i++;
      }
      if (!strcmp(argv[i],"-C")){
        CHANNEL = atoi(argv[i+1]);
        i++;
      }
      if (!strcmp(argv[i],"-F")){
        sprintf(CrateName,"%s",argv[i+1]);
        i++;
      }

      if (!strcmp(argv[i],"h")){
	cout<<"-ND        do not display individual events"<<endl;
        cout<<"-R #     ROC number (default 2)"<<endl;
        cout<<"-S  #    SLOT number (default 3)"<<endl;
        cout<<"-C #     CHANNEL number (default 1)"<<endl;
        cout<<"-V #     data version number (default 1)"<<endl;
	cout<<"-F string  data label (default ROCFDC) \"string\"   "<<endl;
	return 0;
      }
      
    }
  }

  sprintf(InputFile,"adcdata_%s%d_%d_ALLslots_pedestals.dat",CrateName,ROC,version);

  char DataDir[128] = "/gluonraid1/Users/zihlmann/FDCDATA";
  char fnam[128];
  ifstream INF;
  
  char lin[1000];
  int NSlot = -1;
  int NCon = 0;
  int EventNew = 1;
  float x[100000];
  float y[100000];
  float dx[100000];
  float dy[100000];

  sprintf(fnam,"%s/%s",DataDir,InputFile);

  cout<<"open file "<<fnam<<endl;
  INF.open(fnam, ios::in);
  int Slot = SLOT;
  if (Slot<11){
    Slot -= 3;
  } else {
    Slot -= 5;
  }

  if (!INF) {
    cout<<"file not found: "<<fnam<<endl;
  } else {
    
    while (!INF.eof()){
      
      INF.getline(lin,5000);
      string s = lin;
      istringstream iss (s);
      
      int Sl;
      int ch;
      int EVENT;
      float p,dp;
      iss>>EVENT>>Sl>>ch;
      
      if ((Sl==Slot) && (ch+1 == CHANNEL)){
	iss>>p>>dp;
	x[NCon] = (float)EVENT;
	y[NCon] = p;
	dx[NCon] = 0.;
	dy[NCon++] = dp;
      }
      
    }
      
    INF.close();
  }
  
  TGraphErrors *graph;
  char t1[128];
  char t2[128];
  sprintf(t1,"graph%d_%d_%d",ROC,SLOT,CHANNEL);
  sprintf(t2,"Pedestals ROC %d SLOT %d CHANNEL %d",ROC,SLOT,CHANNEL);
  graph = new TGraphErrors(NCon,x,y,NULL,dy);
  graph->SetTitle(t2);
  graph->SetMarkerColor(4);
  graph->SetMarkerStyle(21);
  graph->GetYaxis()->SetTitle("Pedestal [ADC counts]"); 

  if (DISPLAY){
    TCanvas *myc1;
    char str1[128];
    sprintf(str1,"Pedestals ROC %d ADC125 Samples in SLOT %d CHANNEL %d",ROC,SLOT,CHANNEL);
    myc1 = new TCanvas("myc1", str1, 800, 400);
    myc1->SetFillColor(42);
    gPad->SetGrid();
    graph->Draw("AP");
    myc1->Update();

    cout<<"Continue: ";
    char inp[128];
    cin>>inp;
  }

  TFile *fout;
  sprintf(fnam,"%s/PedestalsROC%dSLOT%dCHANNEL%d.root",DataDir,ROC,SLOT,CHANNEL);
  fout = new TFile(fnam,"RECREATE");
  graph->Write();

  fout->Close();

  return 0;
  
}