PndSdsHybridHitProducer.h

Go to the documentation of this file.

#ifndef PNDSDSHYBRIDHITPRODUCER_H
#define PNDSDSHYBRIDHITPRODUCER_H

#include "PndSdsTask.h"
#include "PndSdsPixelDigiPar.h"
#include "PndSdsMCPoint.h"
#include "PndSdsPixel.h"
#include "PndSdsDigiPixel.h"
#include "FairGeoVector.h"
#include "FairGeoTransform.h"
#include "FairMCEventHeader.h"
#include "TVector3.h"
#include "TRandom.h"
#include "TGeoMatrix.h"
#include "TGeoBBox.h"
#include "PndGeoHandling.h"
#include "PndSdsChargeConversion.h"
#include "PndSdsTotDigiPar.h"
#include "PndSdsFE.h"

#include "PndSdsDigiPixelWriteoutBuffer.h"

#include <string>
#include <vector>

class TClonesArray;

class PndSdsHybridHitProducer : public PndSdsTask {
 public:
  PndSdsHybridHitProducer();

  PndSdsHybridHitProducer(PndSdsHybridHitProducer &other)
    : PndSdsTask(), fDigiPixelMCInfo(other.fDigiPixelMCInfo), fPointArray(other.fPointArray), fPixelArray(other.fPixelArray), fPixelMCArray(other.fPixelMCArray),
      fDigiPar(other.fDigiPar), fTotDigiPar(other.fTotDigiPar), fChargeConverter(other.fChargeConverter), fMCEventHeader(other.fMCEventHeader), fDataBuffer(other.fDataBuffer),
      flx(other.flx), fly(other.fly), fthreshold(other.fthreshold), fnoise(other.fnoise), fqsigma(other.fqsigma), fcols(other.fcols), frows(other.frows),
      fPixelHits(other.fPixelHits), fGeoH(other.fGeoH), fEventNr(other.fEventNr), fOverwriteParams(other.fOverwriteParams), fTimeOrderedDigi(other.fTimeOrderedDigi),
      fPixelList(other.fPixelList)
  {
  }
  PndSdsHybridHitProducer(const char *name);

  PndSdsHybridHitProducer(Double_t lx, Double_t ly, Double_t threshold, Double_t noise);

  virtual ~PndSdsHybridHitProducer();

  virtual void SetParContainers();
  virtual InitStatus Init();
  virtual InitStatus ReInit();

  virtual void Exec(Option_t *opt);

  virtual void FinishEvent();
  virtual void FinishTask();

  void RunTimeBased() { fTimeOrderedDigi = kTRUE; }

  PndSdsHybridHitProducer &operator=(PndSdsHybridHitProducer &other)
  {
    if (this != &other) // protect against invalid self-assignment
    {
      fDigiPixelMCInfo = other.fDigiPixelMCInfo;
      fPointArray = other.fPointArray;
      fPixelArray = other.fPixelArray;
      fPixelMCArray = other.fPixelMCArray;
      fDigiPar = other.fDigiPar;
      fTotDigiPar = other.fTotDigiPar;
      fChargeConverter = other.fChargeConverter;
      fMCEventHeader = other.fMCEventHeader;
      fDataBuffer = other.fDataBuffer;
      flx = other.flx;
      fly = other.fly;
      fthreshold = other.fthreshold;
      fnoise = other.fnoise;
      fqsigma = other.fqsigma;
      fcols = other.fcols;
      frows = other.frows;
      fPixelHits = other.fPixelHits;
      fGeoH = other.fGeoH;
      fEventNr = other.fEventNr;
      fOverwriteParams = other.fOverwriteParams;
      fTimeOrderedDigi = other.fTimeOrderedDigi;
      fPixelList = other.fPixelList;
    }
    return *this;
  };

 protected:
  Bool_t fPersistance;     // switch to turn on/off storing the arrays to a file
  Bool_t fDigiPixelMCInfo; // switch to turn on/off storing additional MC Info of Digis

  TClonesArray *fPointArray;

  TClonesArray *fPixelArray;

  TClonesArray *fPixelMCArray;

  PndSdsPixelDigiPar *fDigiPar;
  PndSdsTotDigiPar *fTotDigiPar;
  PndSdsChargeConversion *fChargeConverter;
  FairMCEventHeader *fMCEventHeader;

  PndSdsDigiPixelWriteoutBuffer *fDataBuffer;
  // fDetectorType fMCPointType;

  void Register();
  void Reset();
  void ProduceHits();

  TGeoHMatrix GetTransformation(Int_t sensorID);
  void GetLocalHitPoints(PndSdsMCPoint *myPoint, FairGeoVector &myHitIn, FairGeoVector &myHitOut);
  //  PndSdsHit CalcGlobalPoint(std::vector<PndSdsPixel> pixels);
  TVector3 GetSensorDimensions(Int_t sensorID);

  void AddHit(PndSdsPixel &hit, int mcIndex);
  void AddHits(std::vector<PndSdsPixel> *hitList, int mcIndex);

  // PndSdsFE* fFEModel;

  Double_t flx;        // pixel width in x;
  Double_t fly;        // pixel width in y;
  Double_t fthreshold; // pixel threshold in electrons
  Double_t fnoise;     // pixel noise in electrons
  Double_t fqsigma;    // gaussian charge cloud spread
  Int_t fcols;         // pixel columns in one FE
  Int_t frows;         // pixel rows in one FE
  Int_t fPixelHits;
  //  FairMCEventHeader* fMcEventHeader;
  PndGeoHandling *fGeoH; 
  Int_t fEventNr;
  Bool_t fOverwriteParams;
  Bool_t fTimeOrderedDigi; 

  std::vector<PndSdsPixel> fPixelList;

  ClassDef(PndSdsHybridHitProducer, 8);
};

#endif