//#include "/home/penning/diss/caf/myStyle.cpp"
#include <TH1F.h>
#include <TH2.h>
#include <TF1.h>
#include <TStyle.h>
#include <TCanvas.h>
#include <iostream>
#include <sstream>
#include <fstream>
#include <TLegend.h>
#include <TPaveText.h>
#include <TPostScript.h>
#include <TLatex.h>
#include <TFile.h>


#define theCut 3

#define binN 20000
#define binDown 0
#define binUp 20

#define  OutputDirectory "profiles/run"

void calibConst(int _selLayer=8, Int_t ieta=-13, TString _identifier="_aug_zb") {
  gROOT->SetStyle("Plain");
  //  setMyStyle();
  gStyle->SetOptStat(0000);
  // Prepare histograms
  
  TH1F* _hCalib=new TH1F("calibConst","calibConst",64,0,64);
  // Read information about stores
  TString inN="calib/calib";
  inN+=_identifier;
  inN+="/calibLayer";
  inN+=_selLayer;
  inN+=+"_";
  inN+=_identifier;
  inN+=+"_";
  inN+=ieta;
  inN+=".txt";
  float _calibs[65];
  float _calibs_error[65];
  _calibs[0]=-99;
  _calibs_error[0]=-99;
  int _phi, _Bin, _N;
  float _calib, _calib_error;
  cout<<inN<<endl;
  fstream in(inN,ios::in);
  //  assert(in);
  string line;
  while (in.good()) {
    in>>_phi>>_Bin>>_N>>_calib>>_calib_error;
    //    cout<<" "<<_phi<<" "<<_Bin<<" "<<_N<<" "<<_calib<<" "<<_calib_error<<endl;
    _calibs[_phi]=_calib;
    _calibs_error[_phi]=_calib_error;
    if (in.eof()) break;
  }
  in.close();
  // Loop over the 64 iphis
  for (Int_t i=1; i<=64; i++) {
    //    cout<<"phi: "<< i <<", _calib: "<<_calibs[i]<<" +/-" <<_calibs_error[i]<<endl;
    _hCalib->SetBinContent(i,_calibs[i]);
    _hCalib->SetBinError(i,_calibs_error[i]);
    // loop over data
  }
  
  TCanvas* c1=new TCanvas("c1");
  TString titleStr="Calibration constants for #eta="; titleStr+=ieta;
  _hCalib->SetTitle(titleStr);
  _hCalib->GetYaxis()->SetRangeUser(0.5,1.5);
  _hCalib->GetXaxis()->SetTitle("#phi");
  _hCalib->GetYaxis()->SetTitle("calib. const.");
  _hCalib->SetMarkerStyle(8);
  _hCalib->Draw();
  TString strPs="calib/calib"; strPs+=_identifier; strPs+="/calibConstLayer";strPs+=_selLayer; strPs+=+"_"; strPs+=_identifier;strPs+=+"_";strPs+=ieta;strPs+=".eps";
  //  TString strPs="calib/calib"; strPs+=_identifier; strPs+="/calibConstLayer"+_selLayer+"_"+_identifier+"_"+=ieta+=".txt";
  //  cout<<strPs<<endl;
  c1->Print(strPs);
}