#ifndef HShowerHitFinder_H
#define HShowerHitFinder_H

#include "hreconstructor.h"
#include "hlocation.h"
#include "hparset.h"
#include "hratree.h"
#include "hspecgeompar.h"
#include "hgeomvector.h"

class HCategory;
class HIterator;
class HShowerCal;
class HShowerHit;
class HShowerPID;
class HShowerHitHeader;
class HShowerCriterium;
class HSpecGeomPar;

class HShowerHitFinder : public HReconstructor {
public:
        HShowerHitFinder();
        HShowerHitFinder(const Text_t *name,const Text_t *title);
       ~HShowerHitFinder();

        Int_t execute(void);
        Bool_t init(void);
        Bool_t finalize(void);
        Bool_t initParameters(void);

        Bool_t lookForHit(HShowerCal* cal, HLocation &fLoc);
        Float_t calculateSum(HLocation &fLoc, Int_t nRange, Int_t* pncs = NULL);
        Float_t calculateVar(HLocation &fLoc, Int_t nRange, Float_t avg);
        Bool_t isLocalMax(HLocation &fLoc);
        HShowerHitFinder &operator=(HShowerHitFinder &c);

        HCategory* getCalCat(){return m_pCalCat;}
        HCategory* getHitCat(){return m_pHitCat;}
        HCategory* getPIDCat(){return m_pPIDCat;}
        HParSet* getHitFPar(){return m_pHitFPar;}
        HParSet* getGeometry(){return m_pGeometry;}
	HSpecGeomPar* getHSpecGeomPar(){return m_pHSpecGeomPar;}
		

        void setCalCat(HCategory* pCalCat){m_pCalCat = pCalCat;}
        void setHitCat(HCategory* pHitCat){m_pHitCat = pHitCat;}
        void setPIDCat(HCategory* pPIDCat){m_pPIDCat = pPIDCat;}
        void setHitFPar(HParSet* pPar){m_pHitFPar= pPar;}
        void setGeometry(HParSet* pGeometry){m_pGeometry= pGeometry;}
        void setCriterium(HShowerCriterium* pCrit);
	void setHSpecGeomPar(HSpecGeomPar *pSpecGeometry){m_pHSpecGeomPar = pSpecGeometry;}

        Bool_t IsSortFlagSet(){return m_bIsSort;}
        Bool_t IsFillPID(){return m_bIsFillPID;}
        void setSortFlag(Bool_t bSort = kTRUE){m_bIsSort = bSort;}
        void setFillPID(Bool_t bIsFillPID = kTRUE){m_bIsFillPID = bIsFillPID;}

private:
        HLocation m_Loc;
        HRaTree* m_pCellArr;  //random access table for searching local maxima

        HCategory *m_pCalCat; //!Pointer to the cal data category
        HCategory *m_pHitCat; //!Pointer to the hit data category
        HCategory *m_pPIDCat; //!Pointer to the hit data category
        HCategory *m_pHitHdrCat; //!Pointer to the hit header data category

        HParSet *m_pHitFPar; //Pointer to the analysis parameters container
        HParSet *m_pGeometry; //!Pointer to geometry parameters container
	HSpecGeomPar* m_pHSpecGeomPar;

        HIterator *fIter; //!Iterator for calibrated fired pads
        Bool_t  m_bIsFillPID;//flag for filling PID level - default TRUE
        Bool_t  m_bIsSort;//flag for sorting data by m_nAddress - default FALSE

        HShowerCriterium* m_pCriterium;  //!definition of shower criterium

        void fillSums(HShowerHit* hit, HLocation &fLoc);
        virtual void calcCoord(HShowerHit* hit, HLocation &fLoc);
        virtual void calcCoordWithSigma(HShowerHit* hit, HLocation &fLoc, Int_t nRange);


        void fillPID(HShowerHit* hit, HShowerPID* pid);

        HShowerHitHeader* getHitHeader(HLocation &fLoc);
        void updateClusters(HLocation &fLoc);
        void updateLocalMax(HLocation &fLoc);
        void updateFiredCells(HLocation &fLoc);

        ClassDef(HShowerHitFinder,0) //ROOT extension
};

#endif