PndSdsNoiseProducer.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                PndSdsNoiseProducer header file                -----
// -----                  Created 01.07.08  by R.Kliemt                -----
// -------------------------------------------------------------------------

#ifndef PNDSDSNOISEPRODUCER_H
#define PNDSDSNOISEPRODUCER_H

#include "PndSdsTask.h"

#include <vector>

#include "TRandom.h"
//#include "PndSdsGeoPar.h"
#include "PndGeoHandling.h"
#include "PndSdsStripDigiPar.h"
#include "PndSdsPixelDigiPar.h"
#include "PndDetectorList.h"
#include "TString.h"

class TClonesArray;

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

  PndSdsNoiseProducer(PndSdsNoiseProducer &other)
    : PndSdsTask(), fBranchNameStrip(other.fBranchNameStrip), fBranchNamePixel(other.fBranchNamePixel), fDigiStripArray(other.fDigiStripArray),
      fDigiPixelArray(other.fDigiPixelArray), fDigiParRect(other.fDigiParRect), fDigiParTrap(other.fDigiParTrap), fDigiParPix(other.fDigiParPix), fMCPointType(other.fMCPointType),
      fGeoH(other.fGeoH), fPixelIds4(other.fPixelIds4), fPixelIds6(other.fPixelIds6), fPixelIds8(other.fPixelIds8), fPixelIds12(other.fPixelIds12),
      fStripRectIds(other.fStripRectIds), fStripTrapIds(other.fStripTrapIds), fNoiseSpread(other.fNoiseSpread), fThreshold(other.fThreshold),
      fIonizationEnergy(other.fIonizationEnergy){};

  virtual ~PndSdsNoiseProducer();

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

  virtual void SetBranchNames() = 0;
  virtual void SetMCPointType() = 0;

  virtual void Exec(Option_t *opt);

  Double_t CalcDistFraction(Double_t spread, Double_t threshold);
  //   Int_t CalcChanWhite(Int_t chanleft, Double_t frac);
  Int_t CalcChargeAboveThreshold(Double_t spread, Double_t threshold);
  void AddDigiStrip(Int_t &iStrip, Int_t iPoint, Int_t sensorID, Int_t fe, Int_t chan, Double_t charge);
  void AddDigiPixel(Int_t &noisies, Int_t iPoint, Int_t sensorID, Int_t fe, Int_t col, Int_t row, Double_t charge);

  void DiveDownNode(TGeoNode *fN);
  void FillSensorLists();
  //   void Finish();

  PndSdsNoiseProducer &operator=(PndSdsNoiseProducer &other)
  {
    if (this != &other) // protect against invalid self-assignment
    {
      fBranchNameStrip = other.fBranchNameStrip;
      fBranchNamePixel = other.fBranchNamePixel;
      fDigiStripArray = other.fDigiStripArray;
      fDigiPixelArray = other.fDigiPixelArray;
      fDigiParRect = other.fDigiParRect;
      fDigiParTrap = other.fDigiParTrap;
      fDigiParPix = other.fDigiParPix;
      fMCPointType = other.fMCPointType;
      fGeoH = other.fGeoH;
      fPixelIds4 = other.fPixelIds4;
      fPixelIds6 = other.fPixelIds6;
      fPixelIds8 = other.fPixelIds8;
      fPixelIds12 = other.fPixelIds12;
      fStripRectIds = other.fStripRectIds;
      fStripTrapIds = other.fStripTrapIds;
      fNoiseSpread = other.fNoiseSpread;
      fThreshold = other.fThreshold;
      fIonizationEnergy = other.fIonizationEnergy;
    }
    return *this;
  }

 protected:
  TString fBranchNameStrip;
  TString fBranchNamePixel;
  TClonesArray *fDigiStripArray;
  TClonesArray *fDigiPixelArray;

  PndSdsStripDigiPar *fDigiParRect;
  PndSdsStripDigiPar *fDigiParTrap;
  PndSdsPixelDigiPar *fDigiParPix;

  DetectorType fMCPointType;

  PndGeoHandling *fGeoH; 

  std::vector<Int_t> fPixelIds4;
  std::vector<Int_t> fPixelIds6;
  std::vector<Int_t> fPixelIds8;
  std::vector<Int_t> fPixelIds12;
  std::vector<Int_t> fStripRectIds;
  std::vector<Int_t> fStripTrapIds;

  Int_t fNoiseSpread;
  Int_t fThreshold;

  Double_t fIonizationEnergy; // electrons per GeV (Panda unified units)
                              //   void Register();
                              //   void Reset();
                              //   void ProduceHits();

  ClassDef(PndSdsNoiseProducer, 1);
};

#endif