#ifndef AliAnalysisTaskEMCALPi0V2ShSh_h
#define AliAnalysisTaskEMCALPi0V2ShSh_h

// $Id: AliAnalysisTaskEMCALPi0V2ShSh.h$

class TH1F;
class TH1D;
class TH2F;
class THnSparse;
class TList;
class TObjArray;
class AliEMCALGeometry;
class AliESDEvent;
class AliESDtrack;
class AliVCluster;
class AliCentrality;

#include "AliAnalysisTaskSE.h"


class AliAnalysisTaskEMCALPi0V2ShSh : public AliAnalysisTaskSE {
 public:
  AliAnalysisTaskEMCALPi0V2ShSh();
  AliAnalysisTaskEMCALPi0V2ShSh(const char *name);
  virtual ~AliAnalysisTaskEMCALPi0V2ShSh() {}

  void                   UserCreateOutputObjects();
  void                   UserExec(Option_t *option);
  void                   VZEROEventPlane();
  void                   FillClusterHists();
  void                   FillTrackHists();
  void                   Terminate(Option_t *);

 protected:
  TRefArray             *fCaloClusters;          //!pointer to EMCal clusters
  AliEventplane         *fEventPlane;
  Double_t               fCentralityV0M;
  
  
 private:
  AliESDEvent *fESD;      //! ESD object
  TList       *fOutputList; //! General Output list

  Double_t    fEPTPC;
  Double_t    fEPTPCResolution;
  Double_t    fEPV0;
  Double_t    fEPV0A;
  Double_t    fEPV0C;
  Double_t    fEPV0Ar;
  Double_t    fEPV0Cr;
  Double_t    fEPV0r;
  Double_t    fEPV0A4r;
  Double_t    fEPV0A5r;
  Double_t    fEPV0A6r;
  Double_t    fEPV0A7r;
  Double_t    fEPV0C0r;
  Double_t    fEPV0C1r;
  Double_t    fEPV0C2r;
  Double_t    fEPV0C3r;

  //histograms
  TH1F        *fHistAllcentV0;
  TH1F        *fHistAllcentV0r;
  TH1F        *fHistAllcentV0A;
  TH1F        *fHistAllcentV0C;
  TH1F        *fHistAllcentTPC;

  TH2F        *fHistEPTPC;
  TH2F        *fHistEPTPCResolution;

  TH2F        *fHistEPV0;
  TH2F        *fHistEPV0A;
  TH2F        *fHistEPV0C;
  TH2F        *fHistEPV0Ar;
  TH2F        *fHistEPV0Cr;
  TH2F        *fHistEPV0r;
  TH2F        *fHistEPV0A4r;
  TH2F        *fHistEPV0A7r;
  TH2F        *fHistEPV0C0r;
  TH2F        *fHistEPV0C3r;

  TH2F        *fHistdifV0A_V0C0r;
  TH2F        *fHistdifV0A_V0C3r;
  TH2F        *fHistdifV0C0r_V0C3r;
  TH2F        *fHistdifV0C_V0A4r;
  TH2F        *fHistdifV0C_V0A7r;
  TH2F        *fHistdifV0A4r_V0A7r;
  TH2F        *fHistdifV0Ar_V0Cr;	

  TH1F        *fHistClusterEta;
  TH1F        *fHistClusterPhi;
  TH1F        *fHistClusterE;
  TH1F        *fHistClusterEt;
  TH1F        *fHistClusterN;
  TH1F        *fHistClusterM02;
  TH2F        *fHistClusterEN;
  TH2F        *fHistClusterEM02;
  TH2F        *fHistClusterPhiEta;
  TH2F	      *fHistClusterEtN;
  TH2F        *fHistClusterEtM02;
  TH1D        *fHistClusterdphiV0;

  TH1F        *fHistTrackPt;
  TH1F        *fHistTrackEta;
  TH1F        *fHistTrackPhi;
  TH2F        *fHistTrackPhiEta;

  THnSparse   *fClusterPbV0;
  THnSparse   *fClusterPbV0A;
  THnSparse   *fClusterPbV0C;
  THnSparse   *fClusterPbTPC;
 

  AliAnalysisTaskEMCALPi0V2ShSh(const AliAnalysisTaskEMCALPi0V2ShSh&); // not implemented
  AliAnalysisTaskEMCALPi0V2ShSh& operator=(const AliAnalysisTaskEMCALPi0V2ShSh&); // not implemented
  
  ClassDef(AliAnalysisTaskEMCALPi0V2ShSh, 1);
};
#endif