PndGemDigitize.h

Go to the documentation of this file.

//* $Id: */

// -------------------------------------------------------------------------
// -----                    PndGemDigitize header file                 -----
// -----                  Created 12/02/2009 by R. Karabowicz          -----
// -------------------------------------------------------------------------


#ifndef PNDGEMDIGITIZE_H
#define PNDGEMDIGITIZE_H 1

#include <PndPersistencyTask.h>
#include "TStopwatch.h"
#include "TRandom2.h"
#include "PndGemMCPoint.h"
#include "PndGemSensor.h"

#include <list>
#include <map>

class TClonesArray;
class PndGemDigiPar;
class PndGemDigiWriteoutBuffer;


class PndGemDigitize : public PndPersistencyTask
{

 public:

  PndGemDigitize();


  PndGemDigitize(Int_t iVerbose);


  PndGemDigitize(const char* name, Int_t iVerbose);


  virtual ~PndGemDigitize();


  virtual void Exec(Option_t* opt);


  void SaveOutsideHits(Bool_t bt=kTRUE) {fSaveOutsideHits = bt;}


  void SetRealisticResponse(Bool_t bt=kTRUE) {fRealisticResponse = bt;}

  void RunTimeBased(Bool_t bt=kTRUE) {fTimeOrderedDigi   = bt;}

  Bool_t GetTimeBased()   const { return fTimeOrderedDigi; }

  void SetShowerSigma(Double_t dd) {fSSigma = dd;} //sigam for realistic response
  
 private:

  PndGemDigiPar* fDigiPar;
  TClonesArray*     fPoints;       
  TClonesArray*     fDigis;        
  TClonesArray*     fDigiMatches;  
  Double_t          fSSigma;
  Bool_t fSaveOutsideHits; 
  TClonesArray* fHitOutsideArray; 
  Bool_t fRealisticResponse;
  Int_t             fTNofEvents;
  Int_t             fTNofPoints;
  Int_t             fTNofDigis;
  TRandom2*         fRand;
  PndGemDigiWriteoutBuffer* fDataBuffer;
  Bool_t fTimeOrderedDigi;

  Int_t             fNPoints;
  Int_t             fNFailed;
  Int_t             fNOutside;
  Int_t             fNMulti;
  Int_t             fNDigis;

  TStopwatch        fTimer;

  std::map<std::pair<Int_t, Int_t>, Int_t> fChannelMap; 
 
 
  virtual void SetParContainers();


  virtual InitStatus Init();


  virtual InitStatus ReInit();


  void Reset();


  virtual void Finish();


  void DigitizeEvent();

  void DigitizeRealisticEvent();

  void SimulateRectangularResponse(Int_t sensorDetId, Int_t side, Double_t channelInd, Double_t stripWidth, Double_t showerSigma, Double_t showerStrength, Int_t iPoint);

  void SimulateGaussianResponse(PndGemSensor* sensor, Int_t side, PndGemMCPoint* gemPoint, Double_t showerSigma, Double_t showerStrength, Int_t iPoint);

  void ActivateChannel(Int_t sensorDetId, Int_t sensorSide, Int_t channelNumber, Double_t signalHeight, Double_t signalTime, Int_t iPoint);

  void PrintDigis();

  ClassDef(PndGemDigitize,1);

};

#endif