// File: hrichpedrawnew.h
//
// Author: Angelika Elhardt  <angelika.elhardt@physik.tu-muenchen.de>
// Last update by W.Schoen: 00/11/08 22:00:32
// Modified by martin.jurkovic@ph.tum.de 2005/09/11

#ifndef HRICHPEDRAWNEW_H
#define HRICHPEDRAWNEW_H

#include "TArrayS.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TH2.h"
#include "TNtuple.h"
#include "TStyle.h"
#include "TVector.h"
#include "TString.h"

#include "haddef.h"
#include "hades.h"
#include "hevent.h"
#include "hparset.h"
#include "hreconstructor.h"
#include "hrichcalpar.h"

class HCategory;
class HIterator;
class HRichPedRawNew : public HReconstructor {

 public:
  HRichPedRawNew(const Text_t *name,const Text_t *title,const Char_t *,const Char_t*);
  ~HRichPedRawNew();

  void createCalibration();
  void createPedestal();
  void setCalPar(HParSet* pCalPar) {if(fCalPar) delete fCalPar; fCalPar = pCalPar;}

  Bool_t finalize(void);
  Bool_t init(void);

  Int_t execute(void);
  Int_t getCol(int);
  Int_t getRow(int);

  HParSet *getCalPar() {return fCalPar;}

  TFile *dataFile;

  TH2F *meanHisto[6];
  TH1F *mean1dHisto[6];
  TH2F *sigmaHisto[6];
  TH1F *sigma1dHisto[6];
  TH2F *statHisto[6];
  TH2F *runAwayHisto[6];


 private:
  void bookCanvases();
  void bookHistos();
  void calculate();
  void createMeans();
  void createStatistics();
  void deleteMeans();
  void deleteStatistics();
  void initializeMeans();
  void initializeMeanSigma();
  void initializeStatistics();
  void initCalPar();
  void initMappingPar();
  void fillMask();
  void fill2DHistos();
  void fill1DHistos();
  void meanAndSigma(int,int,int,float);
  void outputFile();
  void printCalParCntr();
  void uipFile();

  Bool_t fillCalParCntr();
  Bool_t checkSectorList();

  Double_t (*pDiff)      [92][90];
  Double_t (*pMean)      [92][90];
  Double_t (*pSigma)     [92][90];

  Int_t    (*countArray) [92][90];
  Int_t    (*pRunAway)   [92][90];
  Int_t    (*pStat)      [92][90];

  Bool_t checkInput;
  Bool_t sigmaTransformation;

  TString pedPathName;
  Char_t sectorList[7];

  Int_t channels;
  Int_t colMax;
  Int_t hardwareMax;
  Int_t mask[6][2];
  Int_t nCounter;
  Int_t numSec;
  Int_t rowMax;
  Int_t sectorMax;

  Float_t epsilon;
  Float_t size;

  HCategory *pRichRaw;

  HIterator *pIter; 

  HParSet *fCalPar; //Pointer to the cal data category
  HParSet *fMappingPar;

  TVector *vec3;
  TVector *vec5;
  TVector *vec6;

  ClassDef(HRichPedRawNew,0)
};

#endif /* !HRICHPEDRAWNEW_H */