CbmPsdMCbmHitProducer.h

Go to the documentation of this file.

/* Copyright (C) 2020 Facility for Antiproton and Ion Research in Europe, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Pierre-Alain Loizeau [committer] */
 
#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