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 "TRandom.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;
  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