//--------------------------------------------------------------------------
// Description:
//	Class EmcERatioLocMaxMaxFinder./
//      Searches for local maxima in a cluster based on the ratio
//      between the energy of the maxima crystal and that of 
//      its neighbours 
//
// Environment:
//	Software developed for the BaBar Detector at the SLAC B-Factory.
//
// Adapted for the PANDA experiment at GSI
//
// Author List:
//      Phil Strother         
//      Helmut Schmuecker                   Ruhr-Uni Bochum      
//
// Copyright Information:
//	Copyright (C) 1997	            Imperial College
// Modified:
// M. Babai
//------------------------------------------------------------------------
#pragma once
#ifndef PNDEMC2DLOCALMAXFINDER_H
#define PNDEMC2DLOCALMAXFINDER_H

#include "FairTask.h"
#include "TObject.h"
#include "PndEmcDataTypes.h"

class PndEmcDigi;
class PndEmcCluster;
class PndEmcTwoCoordIndex;
class PndEmcGeoPar;
class PndEmcDigiPar;
class PndEmcRecoPar;

class PndEmc2DLocMaxFinder: public FairTask
{
  
  
 public:
  
  // Constructors
  PndEmc2DLocMaxFinder(Int_t verbose=0);
  
  // Destructor
  virtual ~PndEmc2DLocMaxFinder();
  
  /** Virtual method Init **/
  virtual InitStatus Init();

  /** Virtual method Exec **/
  virtual void Exec(Option_t* opt);
  
  void SetStorageOfData(Bool_t p = kTRUE) {fPersistance=p;};
  
 protected:
  
  // Helper functions
  virtual bool isALocalMax( const PndEmcDigi *const, const PndEmcCluster * const, 
			    const PndEmcCoordIndexSet &amongstTheseNeighbours ) const;
  
  // Test to see whether the given Digi is a local maxima, within the set of
  // crystals given in "amongstTheseNeighbours".  Note that this set need
  // all be crystals that actually were within the cluster.
  
 private:
  // Methods
  void getNeighbourDigis( PndEmcCoordIndexSet &, PndEmcCoordIndexSet &, int,
			  std::map<Int_t, Int_t>) const;
  
  /** Input array of PndEmcClusters **/
  TClonesArray* fClusterArray;
  /** Input array of PndEmcDigis **/
  TClonesArray* fDigiArray;
  
  PndEmcGeoPar*     fGeoPar;       /** Geometry parameter container **/
  PndEmcDigiPar*    fDigiPar;      /** Digitisation parameter container **/
  PndEmcRecoPar*    fRecoPar;      /** Reconstruction parameter container **/
  /** Get parameter containers **/
  virtual void SetParContainers();
  
  Bool_t fPersistance;
  
  // Data members
  Double_t fMaxECut;
  Double_t fNeighbourECut;
  Double_t fCutSlope;
  Double_t fCutOffset;
  Double_t fERatioCorr;
  Int_t fTheNeighbourLevel;
  /** Verbosity level **/
  Int_t fVerbose;
 
  PndEmc2DLocMaxFinder(const  PndEmc2DLocMaxFinder& L);
  PndEmc2DLocMaxFinder& operator= (const  PndEmc2DLocMaxFinder&) {return *this;};
 
  ClassDef(PndEmc2DLocMaxFinder,1);
};
#endif // PNDEMC2DLOCALMAXFINDER_HH