StripPulseFitter.h

#pragma once
 
#include <vector>
 
#include "Tracking/Digitization/PulseShape.h"
 
namespace tracking::digitization {
 
class StripPulseFitter {
 public:
  struct FitResult {
    double amplitude{
        0};          
    double t0{0};    
    double chi2{0};  
    int ndf{0};      
    bool converged{false};  
  };
 
  StripPulseFitter(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;
 
 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 tracking::digitization