SimQIE.h

#ifndef TRIGSCINT_SIMQIE_H
#define TRIGSCINT_SIMQIE_H
 
#include <iostream>
 
#include "TMath.h"
#include "TRandom3.h"
#include "TrigScint/QIEInputPulse.h"
 
namespace trigscint {
 
class SimQIE {
 public:
  SimQIE();
 
  SimQIE(float pd, float sg, uint64_t seed);
 
  void setTDCThreshold(const float thr) { tdc_thr_ = thr; }
 
  void setGain(const float gg = 1e+6) { gain_ = gg * 16e-5; }
 
  void setFreq(const float sf = 40) { tau_ = 1000 / sf; }
 
  void setNTimeSamples(const int maxts = 5) { maxts_ = maxts; }
 
  float QErr(float Q);
 
  int Q2ADC(float Charge);
 
  float ADC2Q(int ADC);
 
  int TDC(QIEInputPulse* pp, float T0);
 
  std::vector<int> Out_ADC(QIEInputPulse* pp);
 
  std::vector<int> Out_TDC(QIEInputPulse* pp);
 
  std::vector<int> CapID(QIEInputPulse* pp);
 
  bool PulseCut(QIEInputPulse* pulse, float cut);
 
 private:
  int nbins_[5] = {0, 16, 36, 57, 64};
  float edges_[17] = {-16,   34,    158,    419,    517,   915,
                      1910,  3990,  4780,   7960,   15900, 32600,
                      38900, 64300, 128000, 261000, 350000};
  float sense_[16] = {3.1,   6.2,   12.4,  24.8, 24.8, 49.6, 99.2, 198.4,
                      198.4, 396.8, 793.6, 1587, 1587, 3174, 6349, 12700};
 
  float gain_{1};
  float tau_{25};
  int maxts_{0};
 
  float tdc_thr_{3.74};
 
  std::unique_ptr<TRandom3> rand_ptr{nullptr};
  TRandom3* trg_;
 
  float mu_{0};
  float sg_{0};
  bool isnoise_{false};
};
}  // namespace trigscint
#endif