#ifndef HEMCDIGIPAR_H
#define HEMCDIGIPAR_H

#include "hparcond.h"
#include "TArrayF.h"
#include "TArrayD.h"
#include "TF1.h"

class HEmcDigiPar : public HParCond {
protected:
  Float_t sigmaT;       // time resolution [ns]
  Float_t phot2E;       // mean energy deposition in 1.5" PMT per photon electron (1000./1306.) in MeV
  Float_t phot2E2;      // mean energy deposition in 3" PMT per photon electron (1000./1306.) in MeV
  Float_t sigmaEIntern; // internal energy resolution
  Float_t sigmaEReal;   // measured energy resolution in 1.5" PMT
  Float_t sigmaEReal2;  // measured energy resolution in 3" PMT
  TArrayF enCorrection; // cell's energy correction values: geant_en * (1.0 + value)
  TArrayD sigmaTFunc;   // TF1 pars for sigma Time in 1.5" PMT
  TArrayD sigmaTFunc2;  // TF1 pars for sigma Time in 3" PMT
  TF1*    fsigmaTFunc;  //! TF1 pars for sigma Time in 1.5" PMT : [p0]+[p1]/(x+[p2])
  TF1*    fsigmaTFunc2; //! TF1 pars for sigma Time in 3"   PMT : [p0]+[p1]/(x+[p2])

public:
  HEmcDigiPar(const char* name="EmcDigiPar",
              const char* title="Emc digitization parameters",
              const char* context="EmcDigiProduction");
  ~HEmcDigiPar() {;}
  Float_t getSigmaT()       const {return sigmaT;}
  Float_t getPhot2E()       const {return phot2E;}
  Float_t getPhot2E2()      const {return phot2E2;}
  Float_t getSigmaEIntern() const {return sigmaEIntern;}
  Float_t getSigmaEReal()   const {return sigmaEReal;}
  Float_t getSigmaEReal2()  const {return sigmaEReal2;}
  Float_t getEnCorr(Int_t cell) const {return enCorrection[cell];}
  Double_t getSigmaTPars (Int_t i) const {return (i >= 0 && i<sigmaTFunc .GetSize()-1) ? sigmaTFunc [i]: -1000.;}
  Double_t getSigmaTPars2(Int_t i) const {return (i >= 0 && i<sigmaTFunc2.GetSize()-1) ? sigmaTFunc2[i]: -1000.;}

  void setEnCorr(Int_t cell, Float_t corr) {enCorrection[cell] = corr;}
  void setSigmaTPars (Int_t n, Double_t* pars) { sigmaTFunc .Set(n,pars); }
  void setSigmaTPars2(Int_t n, Double_t* pars) { sigmaTFunc2.Set(n,pars); }
  Double_t calcSigmaT(Int_t type, Double_t energy);


  Bool_t init(HParIo*, Int_t*);
  void   clear(void);
//  void   fill(Float_t, Float_t, Float_t, Float_t, Float_t, Float_t);
  void   fill(Float_t sT, Float_t p2e, Float_t p2e2, Float_t seI, Float_t seR, Float_t seR2);
  void   putParams(HParamList*);
  Bool_t getParams(HParamList*);

  ClassDef(HEmcDigiPar,3) // Container for the EMC digitization parameters
};

#endif /* !HEMCDIGIPAR_H */