PndMdtHitProducer.h

Go to the documentation of this file.

#ifndef PNDMDTHITPRODUCER_H
#define PNDMDTHITPRODUCER_H 1

#include "FairTask.h"
#include "PndMdtHit.h"
#include "TVector3.h"
#include "FairTSBufferFunctional.h"

class TClonesArray;
class PndMdtDigi;

using std::map;
using std::vector;

class PndMdtHitProducer : public FairTask {

 public:
  PndMdtHitProducer();

  ~PndMdtHitProducer();

  virtual InitStatus Init();

  virtual void RunTimeBased() { fTimeOrderedDigi = kTRUE; }
  virtual void FinishTask();
  virtual void Exec(Option_t *opt);

  PndMdtHit *AddHit(Int_t detID, Int_t stripID, TVector3 &pos, TVector3 &dpos, Int_t bIndex, Int_t sIndex);

  //  Bool_t MdtMapping(); // Creates maps of MDT hits
  // void Reset();        // reset maps

 private:
  virtual void Exec_old(Option_t *opt);

  TClonesArray *fBoxArray;
  TClonesArray *fStripArray;
  TClonesArray *fDigiArray;

  TClonesArray *fBoxClusterArray;
  TClonesArray *fStripClusterArray;
  TClonesArray *fHitArray;

  Int_t fNumofBoxDigis;
  Int_t fNumofStripDigis;
  Int_t fNumofHits;
  Int_t fNumofGroupedBoxDigis;
  Int_t fNumofGroupedStripDigis;
  /*   map<Int_t, vector<Int_t> >mapBoxBarrel; */
  /*   map<Int_t, vector<Int_t> >mapBoxEndcap; */
  /*   map<Int_t, vector<Int_t> >mapBoxForward; */
  /*   map<Int_t, vector<Int_t> >mapStripBarrel; */
  /*   map<Int_t, vector<Int_t> >mapStripEndcap; */
  /*   map<Int_t, vector<Int_t> >mapStripForward; */
  std::map<Int_t, Int_t> fMatchMap;
  // statistics record
  std::map<Int_t, Int_t> fBoxHitMapofLayer;
  std::map<Int_t, Int_t> fStripHitMapofLayer;

  typedef std::vector<std::vector<PndMdtDigi *>> ClustersCollection;
  typedef ClustersCollection::iterator ClustersColIter;
  typedef std::vector<PndMdtDigi *> SingleCluster;
  typedef SingleCluster::iterator SingleDigiIter;

  typedef std::vector<std::vector<PndMdtDigi *>> TrackletCollection;
  typedef TrackletCollection::iterator TrackletColIter;
  typedef std::vector<PndMdtDigi *> Tracklet;
  typedef Tracklet::iterator TrackletDigiIter;

  typedef std::vector<PndMdtDigi *> VecDigi;
  typedef VecDigi::iterator VecDigiIter;

  Bool_t fTimeOrderedDigi;
  BinaryFunctor *fFunctor;
  Int_t fEventCounter;
  Double_t fTimeWindow;

  struct PndMdtDigiLess {
    bool operator()(const PndMdtDigi *lv, const PndMdtDigi *rv) const;
  };

  ClassDef(PndMdtHitProducer, 1);
};

#endif