////////////////////////////////////
// KRATTA DST analyser
// for the Asy-Eos experiment
// TKratDstAna DECLARATION
// Feb 2013
// revison 02/2013
// E.d.F ver 1.0 
// sebastian.kupny@uj.edu.pl
// Changes:
////////////////////////////////////


#ifndef KRATTA_DST_ANALYSER_H
#define KRATTA_DST_ANALYSER_H 


#include <iostream>
#include <map>
#include <string>
#include <vector>

#include "TClonesArray.h"
#include "TF1.h"
#include "TFumili.h"
#include "TGraph.h"
#include "TH1D.h"
#include "TLine.h"
#include "TMarker.h"
#include "TMath.h"

#include "TGeoManager.h"
#include "FairRootManager.h"
#include "FairRunAna.h"
#include "FairRuntimeDb.h"
#include "FairTask.h"

#include "mktreeraw.h"
#include "mktreedstchain_sct_fitABC.h"

class TKratDstAna : public FairTask
{

  public:
    /** Default constructor **/  
    TKratDstAna();
    
    /** Constructor **/  
    TKratDstAna (const char* name, Int_t iVerbose = 1);
    
    TKratDstAna(Int_t verbose);

    /** Destructor **/
    ~TKratDstAna();

    /** Virtual method Init **/
    virtual InitStatus Init();

    /** Virtual method Exec **/
    virtual void Exec(Option_t* opt);
    virtual void Finish();
    virtual void Reset();
    TString BuildAnalysisName( Int_t startRun, Int_t lastRun, Int_t modNo, TString param1, TString param2 = "");
    /** KRATTA input data container (form the fairroot manager) **/
    void SetKratContInName( TString KratContInName ){
      fKratContInName = KratContInName;
    }
    TString GetKratContInName( void ) const {
      return fKratContInName;
    }    
    void  SetModuleToAnalyse( Int_t moduleToAnalyse){ fModToAnalysisNumber = moduleToAnalyse; };
    Int_t GetModuleToAnalyse( void ){ return fModToAnalysisNumber; };
    
    Int_t LoadCuts( TString A_fileName );
  protected:
  
    Int_t fModToAnalysisNumber;
    Int_t fStartRun;
    Int_t fStopRun;
         
    ASYEvent      *fASYEvent;
    TString        fKratContInName;
    
    TClonesArray  *fASYEventCopy;

    
    TClonesArray  *fKrattaAsyPeakClonesArray;
    
    TFile * fFileWithAnalysisResults;
    
    ///----------------------- Analysis tools ---------------------------
    //TH1F *fHistTimestamp;
    TH1F *fHistEventNo;
    TH1F *fHistPeakNo;
    
    TH2F *fDeltaEE;
    TH2F *fHistogram;
    TH1F *fParam1;
    TH1F *fParam2;
    TH2F *fGammaLine;
    
  

    TCanvas * fC1;
    
    struct {
      TCutG * p;
      TCutG * d;
      TCutG * t;		
      TCutG * He3;
      TCutG * He4;
      TCutG * GammaLine;
    } fCuts[35];

    ///----------------------- Analysis functions ----------------------
    
    bool CheckPeak( ASYFadcPeak *peak );
    
    
    
    
  private:
    TKratDstAna(const TKratDstAna&);
    TKratDstAna& operator=(const TKratDstAna&){
      return *this;
    }
      
    virtual void SetParContainers();

 
  ClassDef(TKratDstAna,1);
  
};

#endif ///KRATTA_DST_ANALYSER_H