#ifndef  HMDCCPSELECTOR_H
#define  HMDCCPSELECTOR_H

#include "hreconstructor.h"
#include "TString.h"

class HGeantKine;
class HCategory;
class HMdcSegSim;
class HMdcHitSim;

class TArrayI;
class TNtuple;
class TH1F;

class HMdcCPSelector : public HReconstructor {
 protected: 
  Int_t theta_bin_size;
  Int_t phi_bin_size;
  
  TNtuple * mdc_singles;
  TNtuple * mdc_doubles;
  TH1F * hsystem;
 
  TArrayI* arr_doub_mdcseg;
  TArrayI* arr_sing_mdcseg;

  TArrayI* arr_doub_mdchit0;
  TArrayI* arr_doub_mdchit1;

  TArrayI* arr_sing_mdchit0;
  TArrayI* arr_sing_mdchit1;

  Int_t cout_singles;
  Int_t cout_doubles;

  Int_t cout_singles_m0;
  Int_t cout_singles_m1;

  Int_t cout_doubles_m0;
  Int_t cout_doubles_m1;

  HCategory *fMdcSegCat;  //! Pointer to the category
  HCategory *fMdcHitCat;  //! Pointer to the category
  HCategory *fGeantCat;   //! Pointer to the category
  HCategory *fGeantMdcCat;   //! Pointer to the category

  HCategory *fPidPart; 
  HIterator *fPidPartIter;
  
  HCategory *fClusInf; //! Pointer to the category
  HIterator *iterator_clusinf;  //! Iterator on HMdcClusInf
  HIterator *iterator_kine; //! Iterator on HGeantKine
  HIterator *iterator_mdcseg; //! Iterator on HMdcSegSim
  HIterator *iterator_mdchit; //! Iterator on HMdcHitSim
  
  Int_t  nThreshold,index,dummy,dummy1;
  Int_t pTrack_parent,nEvCounter;
  Int_t pID;
  Int_t special_mode;//=0 it select singles and doubles normaly ;=1 for doubles it consider also such singles that has coresponding pair lepton passing possition  of Hit in Geant   
  
  
  TString nameHistoFile;
  
  Int_t mod_mdcseg;

 public: 
  HMdcCPSelector(void);
  HMdcCPSelector(const Text_t *name,const Text_t *title,const Char_t *histoname="histograms_convanal.root",Int_t mode=1);
  ~HMdcCPSelector(void);
  void resetArray(TArrayI *ar); 
  Bool_t checkArray(Int_t member,TArrayI *ar);
  void put2Array(Int_t member,TArrayI *ar,Int_t index);
  
  Bool_t init(void);
  Bool_t finalize(void);
  Int_t execute(void);
  void findConvInMdcSeg();
  void findConvInMdcSeg(HMdcSegSim* mdcseg);
  void findConvPairsInPidParticle();
  void setHistoName(const Char_t* name){nameHistoFile=name;}
  void setSpecialMode(Int_t mode){special_mode=mode;}
  HCategory * getMdcHitCat(){return  fMdcHitCat;} 
  HCategory * getClusInfCat(){return  fClusInf;} 
  
  Bool_t fillHistoSingle(HMdcSegSim *mdchit); 
  Bool_t fillHistoDouble(HMdcSegSim *mdchit);

  Bool_t fillHistoSingle(HMdcHitSim *mdchit,HMdcSegSim *mdcseg); 
  Bool_t fillHistoDouble(HMdcHitSim *mdchit,HMdcSegSim *mdcseg);


  void writeHisto2File();
  ClassDef(HMdcCPSelector,0)
    };
    
#endif /* !HMDCCPSELECTOR_H */