#ifndef HDEANATASK_H
#define HDEANATASK_H

#include "hlocation.h"
#include "hreconstructor.h"
#include "hcategory.h"
#include "hiterator.h"
#include "hgeantkine.h"

#include "HDEAnaHistManager.h"
#include "HDEAnaCuts.h"


#include "TCanvas.h"
#include "TFile.h"

#include <vector>
#include <string>

using namespace std;


class HDEAnaTask : public HReconstructor {

private:

    HCategory* fCatKine;
    HCategory* fCatParticleCand;
    HCategory* fCatRichHit;

    Int_t fEventNum;

    HDEAnaHistManager* fHM;

    string fOutputDir; // output directory for the images

    string fOutputFileName; // output root file with histogramms


    Double_t fSignalMBR; // signal Multiplicity * Branching Ratio

    HDEAnaCuts fCuts;


    void initHist();

    void fillKineHists();

    void fillOpenAngleToBackgroundHists();

    void fillParticleCandHists();

    void fillGeantPidHists();

    void fillQPHists();

    void fillRichHists();

    void fillDecayHists();

    void fillNofCandsHists();

    void fillNofRingsHists();

    void fillPairSources(double xBin, double yBin, bool isLeptonCut, bool isMomCut, bool isOACut);


    void drawHist();


public:
   HDEAnaTask(const string& outputFileName);
   ~HDEAnaTask();

   Bool_t init();
   Bool_t reinit();
   Int_t  execute();
   Bool_t finalize();

   void setOutputDir(const string& dir) {fOutputDir = dir;}


   ClassDef(HDEAnaTask, 0)

};

#endif // HDEANATASK_H