#ifndef MARS_MPad
#define MARS_MPad

#ifndef MARS_MTask
#include "MTask.h"
#endif

#ifndef MARS_MH
#include "MH.h"
#endif

class TH1D;
class TH2D;
class TH3D;

class MGeomCam;
class MCerPhotEvt;
class MPedPhotCam;
class MMcEvt;
class MSigmabar;
class MParList;
class MBlindPixels;
class MRead;
class MFilterList;


class MPad : public MTask
{
private:
    MGeomCam       *fCam;
    MCerPhotEvt    *fEvt; 
    MSigmabar      *fSigmabar;
    MMcEvt         *fMcEvt;
    MPedPhotCam    *fPed;
    MBlindPixels   *fBlinds;

    TString        fType;           // type of data to be padded
    TFile          *fInfile;        // input file containing padding histograms

    Int_t          fIter;

    Int_t          fErrors[9];
    Int_t          fWarnings[2];

    //----------------------------------
    // plots used for the padding
    // for all plots it is assumed that the pedestal RMS is given in units of "number of photons"
    TH2D  *fHBlindPixIdTheta;    // 2D-histogram (blind pixel Id vs. Theta)
    TH2D  *fHBlindPixIdThetaMC;  // 2D-histogram (blind pixel Id vs. Theta)
    TH2D  *fHBlindPixIdThetaON;  // 2D-histogram (blind pixel Id vs. Theta)
    TH2D  *fHBlindPixIdThetaOFF; // 2D-histogram (blind pixel Id vs. Theta)

    //---------------------
    TH2D  *fHBlindPixNTheta;    // 2D-histogram (no.of blind pixels vs. Theta)
    TH2D  *fHBlindPixNThetaMC;  // 2D-histogram (no.of blind pixels vs. Theta)
    TH2D  *fHBlindPixNThetaON;  // 2D-histogram (no.of blind pixels vs. Theta)
    TH2D  *fHBlindPixNThetaOFF; // 2D-histogram (no.of blind pixels vs. Theta)

    //---------------------
    TH2D  *fHSigmaTheta;       // 2D-histogram (sigmabar vs. Theta)
    TH2D  *fHSigmaThetaMC;     // 2D-histogram (sigmabar vs. Theta)
    TH2D  *fHSigmaThetaON;     // 2D-histogram (sigmabar vs. Theta)
    TH2D  *fHSigmaThetaOFF;    // 2D-histogram (sigmabar vs. Theta)

    //---------------------
    TH3D  *fHDiffPixTheta;     // 3D-histogram (Theta, pixel, sigma^2-sigmabar^2)
    TH3D  *fHDiffPixThetaMC;   // 3D-histogram (Theta, pixel, sigma^2-sigmabar^2)
    TH3D  *fHDiffPixThetaON;   // 3D-histogram (Theta, pixel, sigma^2-sigmabar^2)
    TH3D  *fHDiffPixThetaOFF;  // 3D-histogram (Theta, pixel, sigma^2-sigmabar^2)

    //---------------------
    TH3D  *fHSigmaPixTheta;     // 3D-histogram (Theta, pixel, sigma)
    TH3D  *fHSigmaPixThetaMC;   // 3D-histogram (Theta, pixel, sigma)
    TH3D  *fHSigmaPixThetaON;   // 3D-histogram (Theta, pixel, sigma)
    TH3D  *fHSigmaPixThetaOFF;  // 3D-histogram (Theta, pixel, sigma)

    //---------------------
    TH3D  *fHgMC;           // matrix (Theta, sigbarold, sigbarnew) for MC data
    TH3D  *fHgON;           // matrix (Theta, sigbarold, sigbarnew) for ON data
    TH3D  *fHgOFF;          // matrix (Theta, sigbarold, sigbarnew) for OFF data

    //-------------------------------
    // plots for checking the padding
    TH2D  *fHSigmaPedestal; // 2D-histogram : pedestal sigma after
                                     //                versus before padding
    TH2D  *fHPhotons;       // 2D-histogram : no.of photons after
                                     //                versus before padding
    TH1D  *fHNSB;           // 1D-histogram : additional NSB


    Bool_t Merge2Distributions(TH1D *hista, TH1D * histb, TH1D * histap,
                               TH2D *fHga,  TH2D * fHgb, Int_t nbinssig);


public:
    MPad(const char *name=NULL, const char *title=NULL);
    ~MPad();

    Bool_t MergeONOFFMC(
      TH2D *sigthmc,  TH3D *diffpixthmc, TH3D *sigmapixthmc,
      TH2D *blindidthmc,  TH2D *blindnthmc,
      TH2D *sigthon,  TH3D *diffpixthon, TH3D *sigmapixthon,
      TH2D *blindidthon,  TH2D *blindnthon,
      TH2D *sigthoff=NULL, TH3D *diffpixthoff=NULL,TH3D *sigmapixthoff=NULL,
      TH2D *blindidthoff=NULL, TH2D *blindnthoff=NULL);

    Bool_t MergeONMC(
      TH2D *sigthmc,  TH3D *diffpixthmc, TH3D *sigmapixthmc,
      TH2D *blindidthmc,  TH2D *blindnthmc,
      TH2D *sigthon,  TH3D *diffpixthon, TH3D *sigmapixthon,
      TH2D *blindidthon,  TH2D *blindnthon);

    Bool_t ReadPaddingDist(const char *filein);
    Bool_t WritePaddingDist(const char *fileout);

    void SetDataType(const char *type);   // type of data to be padded

    Int_t PreProcess(MParList *pList);
    Int_t Process();
    Int_t PostProcess();

    ClassDef(MPad, 0)    // task for the padding 
}; 

#endif