////////////////////////////////////
// KRATTA SIM Prameters analysator - events counters
// for the Asy-Eos experiment
// TKratSimKapAna2 DECLARATION
// May 2013
// revison 05/2013
// E.d.F ver 1.0 
// sebastian.kupny@uj.edu.pl
// Changes:
////////////////////////////////////


#ifndef KRATTA_SIM_KAP_ANA2_H
#define KRATTA_SIM_KAP_ANA2_H

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

#include "TClonesArray.h"
#include "TF1.h"
#include "TGraph.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TLine.h"
#include "TMarker.h"

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

#include "R3BMCTrack.h"
#include "R3BMCStack.h"
#include "KrakowArrayPoint.h"


#include "TKratSimKapReaderDataContainer.h"
#include "TKratSimKapAna1DataContainer.h"


class TKratSimKapAna2 : public FairTask
{

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

    /** Destructor **/
    ~TKratSimKapAna2();

    /** Virtual method Init **/
    virtual InitStatus Init();
    
    /** Virtual method  ReInit **/
    virtual InitStatus ReInit();
    
    /** Virtual method Exec **/
    virtual void Exec(Option_t* opt);
    
    /** Virtual method Finish **/
    virtual void Finish();
    
    
    virtual void Reset();
    virtual void PrintResults();
    
    /** Input container #1 - TKratSimKapReaderDataContainer **/
    void SetTKratSimKapReaderContInName( TString TKratSimKapReaderDataContainer ){
      fTKratSimKapReaderDataContainer = TKratSimKapReaderDataContainer;
    }
    TString GetTKratSimKapReaderContInName( void ) const {
      return fTKratSimKapReaderDataContainer;
    }
    /** Input container #2 - TKratSimKapAna1DataContainer **/
    void SetTKratSimKapAna1ContInName( TString TKratSimKapAna1DataContainer ){; 
      fTKratSimKapAna1DataContainer = TKratSimKapAna1DataContainer;             
    }
    TString GetTKratSimKapAna1ContInName( void ) const {
      return fTKratSimKapAna1DataContainer;
    }
    
  protected:
 
    
    TKratSimKapReaderDataContainer *fKapReaderDataCont;
    TClonesArray                   *fKapReaderDataContClones;
    
    TKratSimKapAna1DataContainer   *fKapAna1DataCont;
    TClonesArray                   *fKapAna1DataContClones;
    
    TString fTKratSimKapReaderDataContainer;
    TString fTKratSimKapAna1DataContainer;  
        
    /// Variables used in analysis
    
  private:
    
    TKratSimKapAna2(const TKratSimKapAna2&);
    TKratSimKapAna2& operator=(const TKratSimKapAna2&){
      return *this;
    }
  
    Double_t  fKineticEnergy;
    Double_t  fMomentum;
    TString   fParticle;
    TString   fSimDescription;
    
    Long_t    fEventsInFile;
    Long_t    fCounterReact;
    Long_t    fCounterReactThin;
    Long_t    fCounterReactThick;

    /// Presentation
    TH1F *fHistEnergySumWhenReactionInThin;
    TH1F *fHistEnergySumWhenReactionInThick;
    TH1F *fHistEnergySumNoReaction;

    TH1F *fHistEnergyPrimWhenReactionInThin;
    TH1F *fHistEnergyPrimWhenReactionInThick;
    TH1F *fHistEnergyPrimNoReaction;
    
    TH1F *fHistReactVsEnInAny;
    TH1F *fHistReactVsEnInThin;
    TH1F *fHistReactVsEnInThick;
    TH1F *fHistReactVsEnNo;
    
    
        
    virtual void SetParContainers();
 
  ClassDef(TKratSimKapAna2,1);
  
};

#endif ///KRATTA_SIM_KAP_ANA2_H