#ifndef MARS_MTFillMatrix
#define MARS_MTFillMatrix

#ifndef MARS_MH3
#include "MH3.h"
#endif

class MRead;
class MTask;
class MHMatrix;

class MTFillMatrix : public MParContainer
{
private:
    MH3      *fReference;
    MRead    *fReader;
    MHMatrix *fDestMatrix1;
    MHMatrix *fDestMatrix2;
    Int_t     fNumDestEvents1;
    Int_t     fNumDestEvents2;
    MTask    *fWriteFile1;
    MTask    *fWriteFile2;

    Bool_t CheckResult(MHMatrix *m, Int_t num) const;
    Bool_t WriteMatrix(MHMatrix *m, const TString &fname, Int_t i) const;

public:
    MTFillMatrix(const MH3 *ref=NULL);
    ~MTFillMatrix();

    void SetDestMatrix1(MHMatrix *matrix, UInt_t num=0)
    {
        fDestMatrix1 = matrix;
        if (num>0)
            fNumDestEvents1 = num;
    }
    void SetWriteFile1(MTask *write, UInt_t num=0)
    {
        fWriteFile1 = write;
        if (num>0)
            fNumDestEvents1 = num;
    }
    void SetNumDestEvents1(UInt_t num)
    {
        fNumDestEvents1 = num;
    }
    void SetDestMatrix2(MHMatrix *matrix, UInt_t num=0)
    {
        fDestMatrix2 = matrix;
        if (num>0)
            fNumDestEvents2 = num;
    }
    void SetWriteFile2(MTask *write, UInt_t num=0)
    {
        fWriteFile2 = write;
        if (num>0)
            fNumDestEvents2 = num;
    }
    void SetNumDestEvents2(UInt_t num)
    {
        fNumDestEvents2 = num;
    }

    void SetReader(MRead *task) { fReader = task; }

    Bool_t Process();
    Bool_t WriteMatrix1(const TString &fname) const { return WriteMatrix(fDestMatrix1, fname, 1); }
    Bool_t WriteMatrix2(const TString &fname) const { return WriteMatrix(fDestMatrix2, fname, 2); }
    Bool_t WriteMatrices(const TString &fname) const;

    ClassDef(MTFillMatrix, 0) // Tool to fill matrices (eg. trainings- and test-matrices)
};

#endif