PulseFitter.h

#pragma once
 
#include <vector>
 
#include "Tools/PulseShape.h"
 
namespace ldmx {
 
class PulseFitter {
 public:
  // I like not having underscores for struct members,
  // but that's not settable until a later version of clang-tidy
  // NOLINTBEGIN(readability-identifier-naming)
  struct FitResult {
    double amplitude{
        0};          
    double t0{0};    
    double chi2{0};  
    int ndf{0};      
    bool converged{false};  
  };
  // NOLINTEND(readability-identifier-naming)
 
  PulseFitter(const PulseShape& shape, double t0_offset_ns,
              double sampling_interval_ns, double pedestal_adc,
              double noise_sigma_adc, double t_scan_min_ns = -50.0,
              double t_scan_max_ns = 150.0, double t_scan_step_ns = 1.0);
 
  FitResult fit(const std::vector<short>& samples) const;
 
  FitResult fitFixed(const std::vector<short>& samples, double T) const;
 
 private:
  std::pair<double, double> evalAtT(const std::vector<short>& samples,
                                    double T) const;
 
  const PulseShape* shape_;
  double t0_offset_ns_;
  double sampling_interval_ns_;
  double pedestal_adc_;
  double inv_sigma2_;  
  double t_scan_min_ns_;
  double t_scan_max_ns_;
  double t_scan_step_ns_;
};
 
}  // namespace ldmx