CompositePulse.h

#ifndef RECON_EVENT_COMPOSITEPULSE_H_
#define RECON_EVENT_COMPOSITEPULSE_H_
 
#include <Rtypes.h>  // For ClassDef
 
#include <vector>
 
#include "TF1.h"
 
namespace ldmx {
 
class CompositePulse {
 public:
  CompositePulse(TF1 func, const double& gain, const double& pedestal)
      : pulse_func_{func}, gain_{gain}, pedestal_{pedestal} {}
 
  CompositePulse() = default;
 
  virtual ~CompositePulse() = default;
 
  void clear() {};
 
  void addOrMerge(const std::pair<double, double>& hit, double hit_merge_ns);
 
  double findCrossing(double low, double high, double level,
                      double prec = 0.01);
 
  void setGainPedestal(double gain, double pedestal) {
    gain_ = gain;
    pedestal_ = pedestal;
  }
 
  double operator()(double time) const { return at(time); }
 
  double at(double time) const {
    double signal = gain_ * pedestal_;
    for (auto hit : hits_)
      signal += hit.first * pulse_func_.Eval(time - hit.second);
    return signal;
  };
 
  const std::vector<std::pair<double, double>>& hits() const { return hits_; }
 
 private:
  std::vector<std::pair<double, double>> hits_;
 
  TF1 pulse_func_;
 
  double gain_;
 
  double pedestal_;
 
  ClassDef(CompositePulse, 2);
 
};  // CompositePulse
 
}  // namespace ldmx
 
#endif