#ifndef __EEmcIsoSmdCal_h__
#define __EEmcIsoSmdCal_h__

#include <TDirectory.h>

class TH1F;
class TH2F;
class TProfile;

class EEezCluster;
class EEmcDb;

class EEmcIsoSmdCal : public TDirectory {
  
 public:

  EEmcIsoSmdCal( const Char_t *name, const Char_t *title = "SMD calibration using isolated tower response" );
  ~EEmcIsoSmdCal() { /* nada */ };

  Int_t Init();
  void Fill( EEezCluster *cluster );

  void RemakeGains ( Int_t sector, EEmcDb *db );
  void RandomizeGains ( Int_t sector, Float_t sigma, EEmcDb *db );

 private: 

  TH2F *mADC[2]; // ADC response of highest strip in plane
  TH2F *mMIP[6]; // Calibrated MIP response of highest strip in plane (w/ eta cuts)
  TH2F *mETow[2]; // Energy of tower used in single-strip response

  TH1F *mFreqLarge[2]; // Frequency plot of how often each strip is largest

  TProfile *mMIP_pfy[2]; // Profile plots of mMIP U and V

 protected:

  ClassDef(EEmcIsoSmdCal,1);

};

#endif