////////////////////////////////////
// KRATTA CALIB analyser
// for the Asy-Eos experiment
// TKratCalibAna DECLARATION
// Jan 2013
// revison 01/2014
// E.d.F ver 1.0 
// sebastian.kupny@uj.edu.pl
// Changes:
////////////////////////////////////


#ifndef KRATTA_CALIB_ANALYSER_H
#define KRATTA_CALIB_ANALYSER_H 


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

#include "TCanvas.h"
#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 "TProfile.h"

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


#include "TRootKRATEvent.h"
#include "ArrayMeanValues.h"
 
class TKratCalibAna : public FairTask
{

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

    /** Destructor **/
    ~TKratCalibAna();

    /** 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 = "");
    
    void  SetModuleToAnalyse( Int_t moduleToAnalyse){ fModToAnalysisNumber = moduleToAnalyse; };
    Int_t GetModuleToAnalyse( void ){ return fModToAnalysisNumber; };
    
    /** KRATTA_IN **/
    void SetKratContInName( TString KratCalEvtClonArrName ){
      fKratCalEvtClonArrName = KratCalEvtClonArrName;
    }
    TString GetKratContInName( void ) const {
      return fKratCalEvtClonArrName;
    }
    
    //Int_t LoadCuts( TString A_fileName );
  protected:
  
    Int_t fModToAnalysisNumber;
    Int_t fStartRun;
    Int_t fStopRun;
         
    //ASYEvent      *fASYEvent;
    //TString        fASYEventName;
    
    //TClonesArray  *fASYEventCopy;
    //TClonesArray  *fKrattaAsyPeakClonesArray;
    
    
    TRootKRATEvent    *fKratCalEvt;
    //TString            fKratCalEvtName;
                      
    TClonesArray      *fKratCalEvtClonArr;
    TString            fKratCalEvtClonArrName;
    
    TRootKRATParticle *fKratCalPart;
    TClonesArray      *fKratCalPartClonArr;
    
    //TFile * fFileWithAnalysisResults;

    ///----------------------- Analysis tools ---------------------------
    //TH1F *fHistTimestamp;
    //TH1F *fHistEventNo;
    //TH1F *fHistPeakNo;

    //TH2F *fDeltaEE;
    //TH2F *fHistogram;
    //TH1F *fParam1;
    //TH1F *fParam2;
    //TH2F *fGammaLine;

    //TCanvas * fC1;

    ///----------------------- Analysis functions ----------------------

  private:
    TKratCalibAna(const TKratCalibAna&);
    TKratCalibAna& operator=(const TKratCalibAna&){
      return *this;
    }
    

    TProfile *fProf[fNMod][fNPar];
    TH2F     *fHsit[fNMod][fNPar];
    
    virtual void SetParContainers();

    Double_t  GetVal(Int_t paramNumber, TRootKRATEvent *ev, TRootKRATParticle * fk);
    
     Int_t fRunFi;
     Int_t fRunLa;
     
  ClassDef(TKratCalibAna,1);
  
};

#endif ///KRATTA_CALIB_ANALYSER_H