PndSttHitProducerRealFast.h

Go to the documentation of this file.

// Fast Digitization

#ifndef PNDSTTHITPRODUCERREALFAST_H
#define PNDSTTHITPRODUCERREALFAST_H 1


#include <PndPersistencyTask.h>
#include "PndGeoSttPar.h"

#include "TVector3.h"

class PndSttHit;
class PndSttHitInfo;
class TClonesArray;
class TObjectArray;

class PndSttHitProducerRealFast : public PndPersistencyTask
{

 public:

  PndSttHitProducerRealFast();


  virtual ~PndSttHitProducerRealFast();


  virtual InitStatus Init();


  virtual void Exec(Option_t* opt);

  PndSttHit* AddHit(Int_t detID, Int_t tubeID, Int_t iPoint, TVector3& pos, TVector3& dpos, Double_t p, Double_t rsim, Double_t closestDistanceError, Double_t depcharge);
  PndSttHit* AddHit(TClonesArray *hitarray, Int_t detID, Int_t tubeID, Int_t iPoint, TVector3& pos, TVector3& dpos, Double_t p, Double_t rsim, Double_t closestDistanceError, Double_t depcharge);

  PndSttHitInfo* AddHitInfo(Int_t fileNumber, Int_t eventNumber, Int_t trackID, Int_t pointID, Int_t nMerged, Bool_t isFake);

  void FoldZPosWithResolution(Double_t &zpos, Double_t &zposError, TVector3 localInPos, TVector3 localOutPos);
  
  Double_t GetError(Double_t, Int_t);

 
  void SetPersistence(Bool_t persistence) { SetPersistency(persistence); }
  void SwitchOnOverlap() { fOverlap = kTRUE; }

  void SeparateHits() { fSeparate = kTRUE; }

  private: 

  virtual void SetParContainers();

  TClonesArray* fPointArray;

  TClonesArray* fHitArray;  
  TClonesArray* fOverlapHitArray;  

  Bool_t fSeparate;


  TClonesArray* fSttParalHitArray;  
  TClonesArray* fSttParalOverlapHitArray;  
  TClonesArray* fSttSkewHitArray;  
  TClonesArray* fSttSkewOverlapHitArray;  



  TObjArray *fVolumeArray;

  TClonesArray* fHitInfoArray;
  Int_t fevtn;

  PndGeoSttPar*  fSttParameters;
  
  Bool_t  fPersistence; 
  Bool_t fOverlap; 
 //  CHECK added
  TClonesArray* fTubeArray;
  Int_t fGeoType;   //< Distinguishes between old geometry description (1) and new, root-based one (2)

 ClassDef(PndSttHitProducerRealFast,1);

};

#endif