/* CbmMuchFindHitsGem.h
 *
 *@author  V.Singhal <vikas@vecc.gov.in>
 *@version 2.0 (removeing fDaq Time stream mode as per CbmDaqNew Scheme)
 *@define  Clusterization and Hit creation is not dependend of EBE or Time Based mode.
 *         fClusterSeprationTime taken as 100 ns which will work with Time Based mode.
 *@since   22.08.17
 *
 *@author  E.Kryshen <e.kryshen@gsi.de>
 *@version 2.0
 *@since   31.07.12
 *
 *@author  M.Ryzhinskiy <m.ryzhinskiy@gsi.de>
 *@version 1.0
 *@since   11.02.08
 *
 *  Produces hits in GEM-like modules using different
 *  clustering algorithms. Clustering algorithms can be chosen
 *  by the switch SetAlgorithm(Int_t iAlgorithm)
 *  Several clustering algorithms are supported
 *  0 - one hit per fired pad
 *  1 - one hit per cluster
 *  2 - simple cluster deconvolution with charge threshold
 *  3 - local maxima finder
 */

#ifndef CBMMUCHFINDHITSGEM_H
#define CBMMUCHFINDHITSGEM_H 1
#include "FairTask.h"
#include "TString.h"
#include <vector>
#include "CbmMuchDigi.h"
#include "CbmEvent.h"
#include "CbmStsReco.h"

class CbmMuchPad;
class CbmMuchCluster;
class TClonesArray;
class CbmMuchGeoScheme;

// Root file for both the Event by Event Mode and Time Stream mode is same therefore 
// Time based data can be read directly from the Root Tree ( One Tree entry = One Time Slice)
// class CbmTimeSlice;

class CbmMuchFindHitsGem: public FairTask {
 public:
  CbmMuchFindHitsGem(const char* digiFileName, Int_t flag);
  virtual ~CbmMuchFindHitsGem() {}
  virtual void Exec(Option_t* opt);
  void ProcessData(CbmEvent*);
  void SetAlgorithm(Int_t iAlgorithm)             { fAlgorithm = iAlgorithm;}
  void SetThresholdRatio(Double_t thresholdRatio) {fThresholdRatio = thresholdRatio; }
  void SetClusterSeparationTime(Double_t time)    { fClusterSeparationTime = time; }
  // Removing SetDaq functionality as Cluster and Hit Finder algorithm is same for both the Time Based and Event Based mode.
  // void SetDaq(Bool_t daq)                         { fDaq = daq; }
  
 private:
  virtual InitStatus Init();
  void FindClusters(CbmEvent*);
  void CreateCluster(CbmMuchPad* pad);
  void ExecClusteringSimple(CbmMuchCluster* cluster,Int_t iCluster, CbmEvent* event);
  void ExecClusteringPeaks(CbmMuchCluster* cluster,Int_t iCluster, CbmEvent* event);
  void CreateHits(CbmMuchCluster* cluster,Int_t iCluster,CbmEvent* event);
  TString  fDigiFile;                     // Digitization file
  Int_t    fFlag;
  Int_t    fAlgorithm;                    // Defines which algorithm to use
  Double_t fClusterSeparationTime;        // Minimum required time between two clusters
  Double_t fThresholdRatio;               // Charge threshold ratio relative to max cluster charge
  Int_t fEvent;                           // Event counter
  Int_t fNofTimeslices;		    // TimeSlice Counter
  TClonesArray* fDigis;                   // Input array of CbmMuchDigi
  TClonesArray* fEvents;                   // Input array of CbmEvent after Event Building
  std::vector<Int_t> fClusterCharges;          //!
  std::vector<Bool_t> fLocalMax;               //!
  std::vector<CbmMuchPad*> fClusterPads;       //!
  std::vector<std::vector<Int_t> >fNeighbours;      //!
  
  TClonesArray* fClusters;                // Output array of CbmMuchCluster objects
  TClonesArray* fHits;                    // Output array of CbmMuchHit
  CbmMuchGeoScheme* fGeoScheme;           // Geometry scheme
  // auxiliary maps and vectors
  std::vector<Int_t> fDigiIndices;             //!
  std::vector<CbmMuchPad*> fFiredPads;         //!
  ECbmMode fMode = kCbmTimeslice;                 // To Select CbmEvent Branch in Time Based Event building mode
  // Removing SetDaq functionality as Cluster and Hit Finder algorithm is same for both the Time Based and Event Based mode.
  //Bool_t fDaq;                            // To daq or not to daq // TODO
  //CbmTimeSlice* fTimeSlice;               // Time slice object in the DAQ approach
  std::vector<CbmMuchDigi> fDigiData;          // Vector of digis
  UInt_t fuClusters = 0; 				//Number of Clusters.
  
  CbmMuchFindHitsGem(const CbmMuchFindHitsGem&);
  CbmMuchFindHitsGem operator=(const CbmMuchFindHitsGem&);
  
  ClassDef(CbmMuchFindHitsGem,1);
};

#endif