#ifndef ALIANAPARTCORRBASECLASS_H
#define ALIANAPARTCORRBASECLASS_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice     */
/* $Id: $ */

//_________________________________________________________________________
// Base class for analysis algorithms
//-- Author: Gustavo Conesa (INFN-LNF)

#include <cstdlib>

//ROOT
class TClonesArray ;
class TObjArray ;
#include <TList.h> 
#include <TObject.h>

//Analysis
class AliESDCaloCluster;
class AliAODCaloCluster;
class AliAODCaloCells;
#include "AliCaloTrackReader.h"   
#include "AliCaloTrackESDReader.h" 
#include "AliCaloTrackAODReader.h" 
class AliCaloPID ;
class AliFiducialCut ;
class AliIsolationCut ;
class AliMCAnalysisUtils ;
class AliNeutralMesonSelection ;
class AliStack ; 
class AliHeader ; 
class AliGenEventHeader ; 
#include "AliAODPWG4ParticleCorrelation.h"

class AliAnaPartCorrBaseClass : public TObject {
	
public: 
	
  AliAnaPartCorrBaseClass() ; // default ctor
  AliAnaPartCorrBaseClass(const AliAnaPartCorrBaseClass & g) ; // cpy ctor
  AliAnaPartCorrBaseClass & operator = (const AliAnaPartCorrBaseClass & g) ;//cpy assignment
  virtual ~AliAnaPartCorrBaseClass() ; //virtual dtor
  
//	virtual void AddAODCaloCluster(AliAODCaloCluster calo) ;
  virtual void AddAODParticle(AliAODPWG4Particle part) ;
  
//	virtual void ConnectAODCaloClusters();
  virtual void ConnectAODPHOSCells();
  virtual void ConnectAODEMCALCells();
  virtual void ConnectInputOutputAODBranches();
  
  virtual TList * GetCreateOutputObjects()      { return (new TList) ;}
  //virtual TList * GetAnalysisOutputContainer()  { return fAnaOutContainer ;} 
	
  virtual void AddToHistogramsName(TString add) { fAddToHistogramsName = add; }  
  virtual TString GetAddedHistogramsStringToName() {return fAddToHistogramsName ;}
  
  virtual void Init() {;}
  virtual void InitParameters() ;
  
  virtual void Print(const Option_t * ) const ;
  
  virtual void MakeAnalysisFillAOD()  {;}
  
  virtual void MakeAnalysisFillHistograms() {;}
  
  virtual Int_t GetDebug() const  { return fDebug ; }
  virtual void SetDebug(Int_t d)   { fDebug = d ; }
  
  virtual Int_t GetEventNumber() const ;
		
  virtual AliCaloTrackReader * GetReader() const {return fReader ; }
  virtual void SetReader(AliCaloTrackReader * const reader) { fReader = reader ; }
  
  //Calorimeter indexes information
  Int_t GetModuleNumberCellIndexes(const Int_t absId, const TString calo, Int_t & icol, Int_t & irow, Int_t &iRCU) const {
	  return fReader->GetModuleNumberCellIndexes(absId, calo, icol, irow,iRCU);}
  Int_t GetModuleNumber(AliAODPWG4Particle * part) const {
	  return fReader->GetModuleNumber(part);}
  Int_t GetModuleNumber(AliESDCaloCluster * cluster) const {
	  return ((AliCaloTrackESDReader*)fReader)->GetModuleNumber(cluster);}
  Int_t GetModuleNumber(AliAODCaloCluster * cluster) const {
	  return ((AliCaloTrackAODReader*)fReader)->GetModuleNumber(cluster);}
	
  virtual void Terminate(TList * /*outputList*/) {;}
	
  //analysis AOD branch
  virtual TClonesArray * GetCreateOutputAODBranch() ;
  virtual TString GetInputAODName() const {return fInputAODName ; }
  virtual void SetInputAODName(TString name)   { fInputAODName = name; }	
  virtual TString GetOutputAODName()  const {return fOutputAODName ; }
  virtual void SetOutputAODName(TString name)   { fNewAOD = kTRUE ; fOutputAODName = name; }
  virtual Bool_t NewOutputAOD() const {return fNewAOD;}
  virtual TString GetOutputAODClassName() const {return fOutputAODClassName;}
  virtual void SetOutputAODClassName(TString name) {fOutputAODClassName = name; }
  
  virtual TString GetAODObjArrayName() const {return fAODObjArrayName;}
  virtual void SetAODObjArrayName(TString name) {fAODObjArrayName = name; }

  virtual TClonesArray* GetInputAODBranch() const {return fInputAODBranch ;}
  virtual TClonesArray* GetOutputAODBranch() const {if(fNewAOD) return fOutputAODBranch; else return fInputAODBranch ;}
  virtual TClonesArray* GetAODBranch(TString aodBranchName) const ;
	
//	virtual TClonesArray* GetAODCaloClusters() const {return fAODCaloClusters ;}
  virtual TClonesArray* GetAODCaloClusters() const ;
  virtual TClonesArray* GetAODTracks() const ;	
  virtual AliAODCaloCells* GetAODCaloCells() const {return fAODCaloCells ;}
  
  virtual TObjArray* GetAODCTS() const ;
  virtual TObjArray* GetAODEMCAL() const ;
  virtual TObjArray* GetAODPHOS() const ;
  
  virtual TString	GetBaseParametersList();
  
  virtual TNamed * GetEMCALCells() const ;
  virtual TNamed * GetPHOSCells() const ;
  
  virtual AliStack * GetMCStack() const ;
  virtual AliHeader* GetMCHeader() const ;
  virtual AliGenEventHeader* GetMCGenEventHeader() const ;
  
  //Analysis helpers classes pointers setters and getters
  virtual AliCaloPID * GetCaloPID() const {return  fCaloPID ;}
  virtual void SetCaloPID(AliCaloPID * const pid) { fCaloPID = pid ;}
  
  virtual AliFiducialCut * GetFiducialCut() const {return  fFidCut ;}
  virtual void SetFiducialCut(AliFiducialCut * const fc) { fFidCut = fc ;}
  
  virtual AliIsolationCut * GetIsolationCut() const {return  fIC ;}
  virtual void SetIsolationCut(AliIsolationCut * const ic) { fIC = ic ;}
  
  virtual AliMCAnalysisUtils * GetMCAnalysisUtils() const {return  fMCUtils ;}
  virtual void SetMCAnalysisUtils(AliMCAnalysisUtils * const mcutils) { fMCUtils = mcutils ;}	
  
  virtual AliNeutralMesonSelection * GetNeutralMesonSelection() const {return  fNMS ;}
  virtual void SetNeutralMesonSelection(AliNeutralMesonSelection * const nms) { fNMS = nms ;}
  
  virtual Bool_t     IsDataMC() const {return fDataMC ; }
  virtual void SwitchOnDataMC()    {fDataMC = kTRUE ; }
  virtual void SwitchOffDataMC()    {fDataMC = kFALSE ; }
  
  virtual Bool_t IsFiducialCutOn() const {return fCheckFidCut ; }
  virtual void SwitchOnFiducialCut() { fCheckFidCut = kTRUE;}
  virtual void SwitchOffFiducialCut() { fCheckFidCut = kFALSE;}
  
  virtual Bool_t IsCaloPIDOn() const {return fCheckCaloPID ; }
  virtual void SwitchOnCaloPID() { fCheckCaloPID = kTRUE;}
  virtual void SwitchOffCaloPID() { fCheckCaloPID = kFALSE;}
  
  virtual Bool_t IsCaloPIDRecalculationOn() const {return fRecalculateCaloPID ; }
  virtual void SwitchOnCaloPIDRecalculation() { fRecalculateCaloPID  = kTRUE;}
  virtual void SwitchOffCaloPIDRecalculation() { fRecalculateCaloPID  = kFALSE;}
  
  virtual Float_t    GetMaxPt()         const {return fMaxPt ; }
  virtual Float_t    GetMinPt()         const {return fMinPt ; }
  virtual void SetMaxPt(Float_t pt)              {fMaxPt = pt ; }
  virtual void SetMinPt(Float_t pt)              {fMinPt = pt ; }
  void SetPtCutRange(Double_t ptmin, Double_t ptmax)
  {  fMaxPt=ptmax;   fMinPt=ptmin;}
  
  //Histogrammes setters and getters
  //Pt, Energy 
  virtual void SetHistoPtRangeAndNBins(Float_t min, Float_t max, Int_t n) {
    fHistoPtBins = n ;
    fHistoPtMax = max ;
    fHistoPtMin = min ;
  }
  
  Int_t   GetHistoPtBins()  const { return fHistoPtBins; }
  Float_t GetHistoPtMin()   const { return fHistoPtMin ; }
  Float_t GetHistoPtMax()   const { return fHistoPtMax ; }
  
  //Azimuthal angle
  virtual void SetHistoPhiRangeAndNBins(Float_t min, Float_t max, Int_t n) {
    fHistoPhiBins  = n ;
    fHistoPhiMax   = max ;
    fHistoPhiMin   = min ;
  }
  
  Int_t   GetHistoPhiBins()  const { return fHistoPhiBins; }
  Float_t GetHistoPhiMin()   const { return fHistoPhiMin ; }
  Float_t GetHistoPhiMax()   const { return fHistoPhiMax ; }
  
  //Pseudorapidity-rapidity
  virtual void SetHistoEtaRangeAndNBins(Float_t min, Float_t max, Int_t n) {
    fHistoEtaBins = n ;
    fHistoEtaMax  = max ;
    fHistoEtaMin  = min ;
  }
  
  Int_t   GetHistoEtaBins()  const { return fHistoEtaBins; }
  Float_t GetHistoEtaMin()   const { return fHistoEtaMin ; }
  Float_t GetHistoEtaMax()   const { return fHistoEtaMax ; }
  
  //Mass
  virtual void SetHistoMassRangeAndNBins(Float_t min, Float_t max, Int_t n) {
	fHistoMassBins = n ;
	fHistoMassMax  = max ;
	fHistoMassMin  = min ;
  }
	
  Int_t   GetHistoMassBins()  const { return fHistoMassBins ; }
  Float_t GetHistoMassMin()   const { return fHistoMassMin ; }
  Float_t GetHistoMassMax()   const { return fHistoMassMax ; }
	
  //Asymetry
  virtual void SetHistoAsymmetryRangeAndNBins(Float_t min, Float_t max, Int_t n) {
	fHistoAsymBins = n ;
	fHistoAsymMax  = max ;
	fHistoAsymMin  = min ;
  }
	
  Int_t   GetHistoAsymmetryBins()  const { return fHistoAsymBins ; }
  Float_t GetHistoAsymmetryMin()   const { return fHistoAsymMin ; }
  Float_t GetHistoAsymmetryMax()   const { return fHistoAsymMax ; }	
	
 private:    
  
  Bool_t  fDataMC ;             // Flag to access MC data when using ESD or AOD     
  Int_t   fDebug ;              // Debug level
  Bool_t  fCheckFidCut ;        // Do analysis for clusters in defined region         
  Bool_t  fCheckCaloPID ;       // Do analysis for calorimeters
  Bool_t  fRecalculateCaloPID ; // Recalculate PID or use PID weights in calorimeters
  Float_t fMinPt ;              // Maximum pt of (trigger) particles in the analysis
  Float_t fMaxPt ;              // Minimum pt of (trigger) particles in the analysis
	
  AliCaloTrackReader * fReader; // Acces to ESD/AOD/MC data
  
  TClonesArray* fInputAODBranch ;    //! Selected input particles branch
  TString       fInputAODName ;      //  Name of input AOD branch;
  TClonesArray* fOutputAODBranch ;   //! Selected output particles branch
  Bool_t        fNewAOD ;            //  Flag, new aod branch added to the analysis or not.
  TString       fOutputAODName ;     //  Name of output AOD branch;
  TString       fOutputAODClassName; //  Type of aod objects to be stored in the TClonesArray (AliAODPWG4Particle, AliAODPWG4ParticleCorrelation ...)	
  TString       fAODObjArrayName ;   // Name of ref array kept in a TList in AliAODParticleCorrelation with clusters or track references.
  TString       fAddToHistogramsName;// Add this string to histograms name
  
  //TClonesArray* fAODCaloClusters ;  //! selected PHOS/EMCAL CaloClusters
  AliAODCaloCells * fAODCaloCells ; //! selected PHOS/EMCAL CaloCells
  
  //Analysis helper classes access pointers
  AliCaloPID               * fCaloPID; //! PID calculation
  AliFiducialCut           * fFidCut;  //! Acceptance cuts
  AliIsolationCut          * fIC;      //! Isolation cut 
  AliMCAnalysisUtils       * fMCUtils; //! MonteCarlo Analysis utils 
  AliNeutralMesonSelection * fNMS;     //! Neutral Meson Selection
  
  //TList * fAnaOutContainer;	// Temporal histogram output container, contents to be added to the main container passed to the main analysis frame

  //Histograms binning and range    
  Int_t   fHistoPtBins   ;  // Number of bins in pt axis
  Float_t fHistoPtMax    ;  // Maximum value of pt histogram range
  Float_t fHistoPtMin    ;  // Minimum value of pt histogram range
  Int_t   fHistoPhiBins  ;  // Number of bins in phi axis
  Float_t fHistoPhiMax   ;  // Maximum value of phi histogram range
  Float_t fHistoPhiMin   ;  // Minimum value of phi histogram range
  Int_t   fHistoEtaBins  ;  // Number of bins in eta axis
  Float_t fHistoEtaMax   ;  // Maximum value of eta histogram range
  Float_t fHistoEtaMin   ;  // Minimum value of eta histogram range
  Int_t   fHistoMassBins ;  // Number of bins in mass axis
  Float_t fHistoMassMax  ;  // Maximum value of mass histogram range
  Float_t fHistoMassMin  ;  // Minimum value of mass histogram range
  Int_t   fHistoAsymBins ;  // Number of bins in asymmetry axis
  Float_t fHistoAsymMax  ;  // Maximum value of asymmetry histogram range
  Float_t fHistoAsymMin  ;  // Minimum value of asymmetry histogram range
	
  ClassDef(AliAnaPartCorrBaseClass,6)
    } ;


#endif //ALIANAPARTCORRBASECLASS_H