CbmRichMCbmHitProducer.h

Go to the documentation of this file.

 
#ifndef CBM_RICH_MCBM_HIT_PRODUCER
#define CBM_RICH_MCBM_HIT_PRODUCER
 
#include "CbmRichRecGeoPar.h"
#include "FairTask.h"
 
#include "CbmDigiManager.h"  // for ROOTCLING
 
class TClonesArray;
class TVector3;
class CbmEvent;
class CbmRichDigi;
 
class CbmRichMCbmMappingData {
public:
  UInt_t fTrbId;
  UInt_t fChannel;
  Double_t fX;
  Double_t fY;
  Double_t fZ;
};
 
class CbmRichMCbmHitProducer : public FairTask {
public:
  CbmRichMCbmHitProducer();
 
  virtual ~CbmRichMCbmHitProducer();
 
  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 std::string& mappingFile) {
    fMappingFile = mappingFile;
  }
 
 
  void setToTLimits(double low, double high) {
    fToTLimitLow  = low;
    fToTLimitHigh = high;
  }
 
 
  void setToTLimitLow(double low) { fToTLimitLow = low; }
 
 
  void setToTLimitHigh(double high) { fToTLimitHigh = high; }
 
 
  void applyToTCut() { fDoToT = true; }
 
  void DoRestrictToAcc(bool val = true) { fRestrictToAcc = val; }
 
  void DoRestrictToFullAcc(bool val = true) { fRestrictToFullAcc = val; }
 
 
private:
  CbmDigiManager* fDigiMan = nullptr;
  TClonesArray* fRichHits;             // RICH hits
  TClonesArray* fCbmEvents = nullptr;  // CbmEvent for time-based simulations
  bool fDoToT              = false;
  bool fRestrictToAcc      = false;
  bool fRestrictToFullAcc  = false;
  bool fRestrictToAerogelAccDec2019 = false;
  double fToTLimitLow               = 0.;
  double fToTLimitHigh              = 1000.;
 
  std::map<Int_t, CbmRichMCbmMappingData> fRichMapping;
 
  Int_t fEventNum;  // event number
 
  Int_t fNofHits = 0;
 
  Double_t fHitError;
 
  std::string fMappingFile;
 
  void InitMapping();
 
  bool isInToT(const double ToT);
 
  bool RestrictToAcc(TVector3& pos);
  bool RestrictToAcc(Double_t x, Double_t y);
 
  bool RestrictToFullAcc(TVector3& pos);
  bool RestrictToFullAcc(Double_t x, Double_t y);
 
  bool RestrictToAerogelAccDec2019(TVector3& pos);
  bool RestrictToAerogelAccDec2019(Double_t x, Double_t y);
 
  void AddHit(CbmEvent* event,
              TVector3& posHit,
              const CbmRichDigi* digi,
              Int_t index,
              Int_t PmtId);
 
  CbmRichMCbmHitProducer(const CbmRichMCbmHitProducer&);
 
  CbmRichMCbmHitProducer& operator=(const CbmRichMCbmHitProducer&);
 
  ClassDef(CbmRichMCbmHitProducer, 1)
};
 
#endif