rootana/src/TAP_generators/template_fitting/TemplateFitter.h

Go to the documentation of this file.

#ifndef TemplateFitter_h_
#define TemplateFitter_h_

#include "TFitterMinuit.h"
#include "TPulseIsland.h"
#include "TTemplate.h"

#include "definitions.h"

class TemplateFitter {
 public:
  TemplateFitter(std::string detname, int refine_factor);
  ~TemplateFitter();

 private:
  TFitterMinuit* fMinuitFitter;
  void operator=(const TemplateFitter&);

 public:
  int FitPulseToTemplate(const TTemplate* hTemplate, const TH1D* hPulse, const std::string& bankname);
  double GetPedestalOffset() { return fPedestalOffset; }
  double GetAmplitudeScaleFactor() { return fAmplitudeScaleFactor; }
  double GetTimeOffset() { return fTimeOffset; }
  double GetChi2() { return fChi2; }
  double GetNDoF() { return fNDoF; }

 private:
  IDs::channel fChannel;
  double fPedestalOffset;
  double fAmplitudeScaleFactor;
  double fTimeOffset;
  double fChi2;
  int fNDoF; // number of degress of freedom in fit

  double fPedestalOffset_estimate;
  double fAmplitudeScaleFactor_estimate;
  double fTimeOffset_estimate;

  double fPedestalOffset_minimum;
  double fAmplitudeScaleFactor_minimum;
  double fTimeOffset_minimum;

  double fPedestalOffset_maximum;
  double fAmplitudeScaleFactor_maximum;
  double fTimeOffset_maximum;

 public:
  void SetInitialParameterEstimates(double pedestal, double amplitude, double time);

 private:
  int fRefineFactor;
};

#endif

---