#ifndef ALICOMPARISONDCA_H
#define ALICOMPARISONDCA_H

//------------------------------------------------------------------------------
// Class to keep information from comparison of 
// reconstructed and MC particle tracks (DCA - Distance of Closest Approach 
// to the vertex).   
// 
// Author: J.Otwinowski 04/02/2008 
//------------------------------------------------------------------------------

class TFile;
class AliMCInfo;
class AliESDRecInfo;
class AliESDEvent; 
class AliESD;
class AliESDfriend;
class AliRecInfoCuts;
class AliMCInfoCuts;
class TH1I;
class TH3F;
class TH3;
class TProfile;
class TProfile2D;
class TString;
class AliESDVertex;

#include "TNamed.h"
#include "AliComparisonObject.h"

//class AliComparisonDCA : public TNamed {
class AliComparisonDCA : public AliComparisonObject {
public :
  AliComparisonDCA(); 
  ~AliComparisonDCA();

  // Init data members
  virtual void Init();

  // Execute analysis
  virtual void Exec(AliMCInfo* infoMC, AliESDRecInfo *infoRC);

  // Merge output objects (needed by PROOF) 
  virtual Long64_t Merge(TCollection* list);

  // Analyse output histograms
  virtual void Analyse();

  // Get analysis folder
  virtual TFolder* GetAnalysisFolder() {return fAnalysisFolder;}

  // Create folder for analysed histograms
  TFolder *CreateFolder(TString folder = "folderDCA",TString title = "Analysed DCA histograms");

  // Export objects to folder
  TFolder *ExportToFolder(TObjArray * array=0);

  // Process events
  void  Process(AliMCInfo* infoMC, AliESDRecInfo *infoRC);

  // Selection cuts
  void SetAliRecInfoCuts(AliRecInfoCuts* cuts=0) {fCutsRC = cuts;}
  void SetAliMCInfoCuts(AliMCInfoCuts* cuts=0) {fCutsMC = cuts;}  

  AliRecInfoCuts*  GetAliRecInfoCuts() const {return fCutsRC;}
  AliMCInfoCuts*   GetAliMCInfoCuts()  const {return fCutsMC;}

  // getters
  TH3F  *GetD0TanSPtTPCITS() {return fD0TanSPtTPCITS;}
  TH3F  *GetD1TanSPtTPCITS() {return fD1TanSPtTPCITS;}
  TH3F  *GetD0TanSPt() {return fD0TanSPt;}
  TH3F  *GetD1TanSPt() {return fD1TanSPt;}
  TH3F  *GetD0TanSPtTPC() {return fD0TanSPtTPC;}
  TH3F  *GetD1TanSPtTPC() {return fD1TanSPtTPC;}

private:
  // DCA resolution
  TH3F  *fD0TanSPtTPCITS; //-> distance to vertex y (TPC+ITS clusters) 
  TH3F  *fD1TanSPtTPCITS; //-> distance to vertex z (TPC+ITS clusters) 
  TH3F  *fD0TanSPt;     //-> distance to vertex y  
  TH3F  *fD1TanSPt;     //-> distance to vertex z 
  TH3F  *fD0TanSPtTPC;  //-> distance to vertex y (only TPC track parameters) 
  TH3F  *fD1TanSPtTPC;  //-> distance to vertex z (only TPC track parameters)

  // Global cuts objects
  AliRecInfoCuts*  fCutsRC; // selection cuts for reconstructed tracks
  AliMCInfoCuts*  fCutsMC;  // selection cuts for MC tracks

  // analysis folder 
  TFolder *fAnalysisFolder; // folder for analysed histograms

  AliComparisonDCA(const AliComparisonDCA&); // not implemented
  AliComparisonDCA& operator=(const AliComparisonDCA&); // not implemented

  ClassDef(AliComparisonDCA,1);
};

#endif