#ifndef PNDSTTRAWANA_H
#define PNDSTTRAWANA_H
	
#include "FairTask.h"
#include "PndSttHit.h"
#include "TMatrixT.h"

class TClonesArray;
class TH1F;
class TH2F;
class TCanvas;
	
	
class PndSttRawAna : public FairTask {
 public:
  PndSttRawAna();

  virtual ~PndSttRawAna();
	
  virtual InitStatus Init();
	
  virtual void Exec(Option_t *option);
	
  virtual void FinishTask();

  void SetPersistence(Bool_t persistence) { fPersistence = persistence; }

  void SetVerbose(int ver) { fVerbose = ver; }

  void Calibration();

  void ReverseTDCspectra();

  void FindDeltaT();

  void WriteOffset();

 private:
  TClonesArray *fRawData;
	
  TFile *fInFile;

  static const Int_t fNLayers = 8;
  static const Int_t fNTubes = 129;
  static const Double_t fRadiusS = 0.5; // cm
  static const Double_t fRadiusW = 0.001; // cm = 10 micro
  Int_t fevtn;
  Bool_t fPersistence;
  Int_t fVerbose;
  Int_t fTInfoBin[fNTubes][3];

  TMatrixT<Double_t> fTInfo;

  TH1F *hsumled;
  TH1F 
    *hrawled[fNTubes],
    *hrawzct[fNTubes],
    *hrawcfd[fNTubes],
    *hrawled0[fNTubes],
    *hrawzct0[fNTubes],
    *hrawcfd0[fNTubes],
    *hfpgaled[fNTubes],
    *hfpgazct[fNTubes],
    *hfpgacfd[fNTubes],
    *hfpgaled0[fNTubes],
    *hfpgazct0[fNTubes],
    *hfpgacfd0[fNTubes];

  TCanvas 
    *crled[fNLayers],
    *crzct[fNLayers],
    *crcfd[fNLayers],
    *cfled[fNLayers],
    *cfzct[fNLayers],
    *cfcfd[fNLayers];

  void CreateHistos();

  void WriteHistos();

  void WriteSumHisto();

  void ClearHistos();
	
 public:
  ClassDef(PndSttRawAna, 0)
    };
	
	
#endif