CbmPsdMCbmHitProducer.h

Go to the documentation of this file.

 
#ifndef CBM_PSD_MCBM_HIT_PRODUCER
#define CBM_PSD_MCBM_HIT_PRODUCER
 
#include "FairTask.h"
//#include "CbmRichRecGeoPar.h"
 
class TClonesArray;
class TVector3;
class CbmEvent;
class CbmDigiManager;
class CbmPsdDigi;
 
class CbmPsdMCbmHitProducer : public FairTask {
public:
  CbmPsdMCbmHitProducer();
 
  virtual ~CbmPsdMCbmHitProducer();
 
  virtual void SetParContainers();
 
  virtual InitStatus Init();
 
 
  virtual void Exec(Option_t* option);
 
  virtual void Finish();
 
 
  void ProcessData(CbmEvent* event);
 
 
  void ProcessDigi(CbmEvent* event, Int_t digiIndex);
 
 
  //void SetMappingFile(const string& mappingFile){fMappingFile = mappingFile;}
 
 
  void setEnRange(double low, double high) {
    fEnLimitLow  = low;
    fEnLimitHigh = high;
  }
 
  void applyEnCut() { fDoEnCut = true; }
 
private:
  CbmDigiManager* fDigiMan = nullptr;
  TClonesArray* fPsdHits;              // PSD hits
  TClonesArray* fCbmEvents = nullptr;  // CbmEvent for time-based simulations
  bool fDoEnCut            = false;
  double fEnLimitLow       = 0.;
  double fEnLimitHigh      = 100000.;
 
  //map<Int_t,CbmRichMCbmMappingData> fPsdMapping; //TODO
 
  Int_t fEventNum;  // event number
 
  Double_t fHitError;
 
  TString fMappingFile;
 
  void InitMapping();
 
  bool isInEnRange(const double energy);
 
  void AddHit(CbmEvent* event,
              Double_t time,
              Double_t energy,
              UInt_t moduleId,
              UInt_t sectionId,
              Int_t index);
 
  CbmPsdMCbmHitProducer(const CbmPsdMCbmHitProducer&);
 
  CbmPsdMCbmHitProducer& operator=(const CbmPsdMCbmHitProducer&);
 
  ClassDef(CbmPsdMCbmHitProducer, 1)
};
 
#endif