// -*- C++ -*-
#ifndef __BTJuly2006__
#define __BTJuly2006__

//////////////////////////////////////////////////////////////////////////
//                                                                      //
//                                                                      //
//                                                                      //
//////////////////////////////////////////////////////////////////////////
#include "/afs/f9.ijs.si//home/irena/Delo/phd_analiza/Beam-Tests/CERNbt_lib/CERNbt.h"
#include "TObject.h"
#include "TH1.h"
#include "TH2.h"
#include "TArray.h"
#include "TCanvas.h"
#include "TGraph.h"
#include "TGraphErrors.h"
#include "TMath.h"
#include "TNtuple.h"
#include "TF1.h"
#include "TFile.h"
#include <iostream>
#include <fstream>
#include <istream>
#include <stdio.h>
#include <stdlib.h>
#include "TROOT.h"
#include "TSystem.h"
#include <cstring>
#include <limits>

using namespace std;

struct EFFdat_t{
  Float_t RunIndx;
  Float_t XX;
  Float_t YY;
  Float_t All[2];
  Float_t Yes[2][6];
};



class BTJuly2006:public TObject
{
public:
  BTJuly2006();
  BTJuly2006(string config,string rootfile);
  ~BTJuly2006();
  TGraphErrors **Get_NINOefficiency(string cut_on_tracks,Int_t drw);
  TH2F ***Get_HitMap();
  TGraph **Draw_Box(Float_t x1, Float_t dx,Float_t y1,Float_t dy,Int_t col,Int_t draw);
  TH1F **Get_NINOSignalSpectra(string cut_on_tracks);
  TH1F **Get_FotecSignalSpectra(string cut_on_tracks);
  TH1F **Get_FotecNoiseSpectra(string cut_on_tracks);
  TH2F **Get_FotecHitMap(Float_t *amplcut);
  TH2F **Get_NINOHitMap(Float_t *amplcut);
  TH2F **ScanZ_HitMap();
  TH2F **Get_Fotec3Defficiency(Float_t *amplcut);
  TH2F ***Get_Fotec3Defficiency_v2(Float_t *amplcut,Float_t BoxWidth, Int_t &ngraphs);
  Double_t Project_Ampl(Int_t ich,Float_t amplcut,Int_t irun);
  TTree *Get_Fotec3Defficiency_v2_Tree(Float_t *MIP,Float_t BoxWidth);

  
  Int_t NCH;
  Int_t NF;
  Char_t path[1000]; 
  std::vector<Int_t> RUN;
  std::vector<Float_t> NINOthr;
  std::vector<Int_t> NINOch;
  std::vector<Int_t> FOTECch;
  CERNbt **datobj;      //!
  TFile **datfile;      //!
  string ROOTFILE;      // READ or WRITE
  string CONFIG;
  string ADC_ChInput[4];  //hwat's connectec to ADC channel, nino fotec moduel...
  Float_t DistZ[4];       //distance from module [ich] to first plane
  //for determinig the box inside the diamond, to do cuts on tracks 
  //when projectiing amplitudes:
  std::vector<std::vector<Float_t> > Pos_x1;
  std::vector<std::vector<Float_t> > Pos_y1;
  Float_t Pos_dx;
  Float_t Pos_dy;

private:
  std::vector<Int_t> StartEvent;
  Double_t Project_Ampl(Int_t ich,Float_t amplcut);
  Double_t Project_Ampl_v2(Int_t ich,Float_t amplcut,Int_t start_run,Int_t end_run);
  ClassDef(BTJuly2006,1) 
};

#endif