PHG4InttDigitizer.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file PHG4InttDigitizer.h

#ifndef G4INTT_PHG4INTTDIGITIZER_H
#define G4INTT_PHG4INTTDIGITIZER_H

#include <phparameter/PHParameterInterface.h>

#include <fun4all/SubsysReco.h>

#include <gsl/gsl_rng.h>

#include <map>
#include <string>                              // for string
#include <utility>                             // for pair
#include <vector>

class PHCompositeNode;

class PHG4InttDigitizer : public SubsysReco, public PHParameterInterface
{
 public:
  PHG4InttDigitizer(const std::string &name = "PHG4InttDigitizer");
  ~PHG4InttDigitizer() override;

  int InitRun(PHCompositeNode *topNode) override;

  int process_event(PHCompositeNode *topNode) override;

  int End(PHCompositeNode *topNode) override;

  void SetDefaultParameters() override;

  void Detector(const std::string &d) { detector = d; }

  void set_adc_scale(const int &layer, const std::vector<double> &userrange);

 private:
  void CalculateLadderCellADCScale(PHCompositeNode *topNode);

  void DigitizeLadderCells(PHCompositeNode *topNode);

  std::string detector;
  // noise electrons
  float added_noise();

  float mNoiseMean;      // Mean of noise electron distribution
  float mNoiseSigma;     // Sigma of noise electron distribution
  float mEnergyPerPair;  // GeV/e-h pair

  // settings
  std::map<int, unsigned int> _max_adc;
  std::map<int, float> _energy_scale;

  // storage
  //SvtxHitMap *_hitmap;

  const unsigned int nadcbins = 8;
  std::map<int, std::vector<std::pair<double, double> > > _max_fphx_adc;

  unsigned int m_nCells;
  unsigned int m_nDeadCells;

  gsl_rng *RandomGenerator;
};

#endif