#ifndef __EEsmdIsoClusters_h__
#define __EEsmdIsoClusters_h__

#include <TDirectory.h>

#include "EEezCluster.h"


class TH1F;
class TH2F;

class EEezClAnalysis;
class EEezAnalysis;
class EEmcSmdMap;

class EEsmdIsoClusters : public TDirectory {

 public:
  
  EEsmdIsoClusters( const Char_t *name, const Char_t *title = "ESMD calibration using isolated gammas" );
  EEsmdIsoClusters() { /* default */ Warning("EEsmdIsoClusters","Use named constructor!"); }
  ~EEsmdIsoClusters() { /* nada */ };

  Int_t Init();
  Int_t Make();

  void  Clear( Option_t *opts );

  void  setClusterFinder ( EEezClAnalysis *c ) { mEEmcClusters = c; }  
  void  setAnalysis ( EEezAnalysis *a ) { mEEmcAnalysis = a; }
  void  setMaxSum( Float_t max ) { mMaxSharedSum = max; }

  
  Int_t getNClusters() { return mClusters.size(); }
  EEezCluster *getClusterPtr( Int_t i ) { return mClusters[i]; }


 private:

  EEezAnalysis   *mEEmcAnalysis; // Reconstructed tower+pre+post+smd response
  EEezClAnalysis *mEEmcClusters; // Reconstructed clusters

  EEezClusterPtrVec_t mClusters;

  EEmcSmdMap *mEEsmdMap; // tower-to-smd (and vice versa) mapping

  Float_t mMaxSharedSum; // Maximum energy sum [GeV] of towers along shared smd strips

 protected:

  ClassDef(EEsmdIsoClusters,1);

};

#endif